1. Introduction
Swine production is facing new challenges to improve its efficiency. One of the current objectives is to improve prolificacy without penalizing the offspring survival and development. The possible damaging effect of the intrauterine growth restriction process (IUGR) on the offspring, mainly linked to large litters, is well-known, as well as the negative consequences on their postnatal development [
1,
2,
3,
4,
5]. Many studies on selected swine breeds have shown the IUGR impact on pig production by inducing lower homogeneity and low birth-weight (BIW) pigs [
6,
7,
8]. These effects have also been described in traditional swine breeds, even with a greater impact than in modern breeds [
9,
10]. Although traditional breeds show lower reproductive parameters [
11,
12,
13,
14], the increase in litter size results in lower BIW and higher variability than in selected swine breeds [
10]. Even the delay of days to market is longer in traditional breeds [
10]. This consequence affects farm profitability and makes clear the wide range of improvement to work on useful strategies, particularly in traditional and fatty breeds, which is a growing industry due to its high-quality cured pork products.
At the sow level, the nutritional strategies are essential to diminish negative consequences in piglets throughout their productive lifetime. Among the prenatal nutritional strategies, amino acids (AA) clearly stand out. Particularly some functional AA, such as arginine [Arg], glutamine [Gln, AA abbreviation], and proline [Pro], have been the target of the largest number of studies on pregnancy in selected swine breeds [
15,
16,
17,
18,
19,
20]. They have shown benefits such as increased fetal-placental unit growth and reduced variation in BIW, with Arg being the most commonly AA used in previous studies [
15,
19,
21,
22,
23,
24,
25,
26,
27]. Glutamine is particularly interested because it is one of the most abundant AA in fetal tissue protein and is a primary energy source for the fetal small intestine [
18,
19]. However, there remain significant factors to be assessed before translating this knowledge into farm strategy to guarantee the best equilibrium between advantages and investments. These factors include the parity number (Pa), the supplementation period, and the current farm conditions, particularly for Gln. Moreover, previous studies were only developed in selected swine breeds. There are no studies on the effect of the prenatal supplement of AA on traditional breeds, although their metabolic differences could produce differences in the response [
28,
29,
30,
31,
32,
33,
34].
Hence, the present study aimed at determining the effects of prenatal Gln supplementation on the prenatal and postnatal developments of the offspring and their carcass and meat quality under farm conditions in a traditional swine breed. The second objective was to examine the interaction of prenatal Gln supplementation and the parity of sow.
2. Materials and Methods
2.1. Animals and Design
The experiment was carried out in a commercial farm, Ibéricos de Arauzo 2004 S.L. (Zorita de la Frontera, Salamanca, Spain), according to the European Union Directive and the Spanish Policy for Animal Protection RD53/2013. As a traditional breed, the Iberian pig was used in this trial. The management of sows and their offspring followed standard commercial farm practices, housing them indoors under controlled temperature, and electronic chip identification was used. Moreover, sows and offspring were fed with standard grain-based diets for Iberian pigs (diets shown in
Table S1; [
35]).
A total of 78 Iberian sows (Retinto variety), from parity (Pa) 1 to 3, were inseminated with cooled semen from Duroc PIC boars (Genus plc, Worcester, UK). Day 36 of pregnancy, 46 pregnant sows (Pa1: 21 gilts, Pa2: 16 sows, Pa3: 9 sows; treated group), randomly chosen, were supplemented with 1% L-glutamine (GLN [treated group abbreviation]; S.A. Ajinomoto OmniChem N.V., Wetteren, Belgium) on the gestation diet up to delivery. The 32 sows of the control group (C; Pa1: 10 gilts, Pa2: 13 sows, Pa3: 9 sows) were fed with the same diet without supplementation. Average daily feed intakes obtained during pregnancy, after Day 35, were between 1 and 1.5 kg with no difference by treatments. Sows were housed in groups from Day 35 to Day 101 of pregnancy and, afterward, individually allocated in pens until weaning. After birth, piglets were sexed and weighed, and a total of 486 alive piglets (C: 198, 53% males; GLN: 288, 49% males) were measured and allocated to mothers until weaning. Males were surgically castrated within the two first days of birth.
Piglets were classified into two birth-weight (BIW) categories, Low and Normal BIW (LBIW and NBIW), for experimental purposes. The classification cut-off value relied on previous studies on the same farm (BIW ≤ 0.99 kg; [
9,
10]). At weaning (average age, 24 days), 136 control (C; 51% males) and 151 treated (GLN; 50% males) piglets randomly selected were weighed and measured. Later, piglets were housed, distributed by treatment and sex in groups of 12 piglets/pen maximum. Pigs were monitored until slaughter. At 215 days-old, 96 control (37% male) and 103 treated pigs (47% males) were sampled. Finally, 54 control (61% males) and 79 treated pigs (62% males) were sampled at the slaughterhouse.
2.2. Birth Data and Offspring Development
Birth data were recorded per sow. The BIW mean and its SD and coefficient of variation (CoV) were calculated per sow using the total piglets born (without stillbirths) and statistically analyzed by parity and treatment. The remaining birth data were assessed based on piglets born alive.
Pigs were individually weighed at birth, weaning, and slaughterhouse. Pigs were shipped to the slaughterhouse in three batches. The first and second batches were slaughtered when pigs reached the minimum market carcass weight established for the Iberian breed (115 kg; Day 244 to 267 of average age). The remaining pigs were sent to market regardless of their body weight on Day 270 of average age. Average daily weight gain (ADWG) was individually calculated for the lactation phase, for the following period (growing and fattening) until the slaughterhouse and for the whole productive life.
At birth and weaning, morphological measurements (occipito-nasal length, biparietal diameter, trunk length, maximum thoracic diameter, and abdominal and thoracic circumferences) were recorded with a measuring tape. At weaning and 215 days-old, backfat thickness, distinguishing the inner and outer layers, and longissimus dorsi (LD, loin) muscle diameter were determined with ultrasound (SonoSite Inc., Bothell, WA, USA) at the P2 point (level of the head of the last rib). At the slaughterhouse, the length of carcasses (from the posterior edge of the symphysis pubica to the anterior edge of the first rib) and the backfat thickness (at the last rib) were measured with a tape. Carcass yield was calculated individually (carcass weight/body weight).
2.3. Tissue Sampling at Slaughter
At birth, 24 control and 24 treated piglets (6 LBIW and 6 NBIW piglets per sex) were euthanized and sampled. To avoid reducing the sample size at weaning, only eight control newborns (4 LBIW and 4 NBIW) were selected from control sows to slaughter at birth. The remaining euthanized control newborns came from sows with similar characteristics to our control group. Selected neonates were euthanized by stunning and exsanguination according to RD53/2013 standard procedures. Immediately head, carcass, and total and individual viscerae (brain, heart, intestine, kidneys, liver, lungs, pancreas, and spleen) were weighed. Later, weight-ratios of carcass and individual viscerae to BIW were calculated. Duodenal samples for gene expression analysis were immediately submerged in RNAlater (Invitrogen, Carlsbad, CA, USA) and stored at −20 °C. Samples of brain, the right lateral lobe of the liver, and LD muscle were also biobanked at −20 °C until fatty acids (FA) composition analysis.
At the slaughterhouse, samples of the right lateral lobe of the liver, LD muscle, and subcutaneous backfat fat at the measurement level (P2 point) were biobanked until FA composition analysis. On the same sampling day, a second sample of LD muscle was also used for the pH and drip-loss analyses [
36]. Duodenal samples were also collected and stored at −80 °C.
2.4. Amino Acids, Glucose and Lipid Metabolism and Oxidative Status
Blood samples were drawn with vacuum tubes (Vacutainer Systems Europe, Meylan, France) from euthanized neonates at birth and after fasting at the slaughterhouse. Samples were centrifuged at 1500×
g for 15 min, and plasma was separated and stored at −80 °C until analyses of AA and glucose (fructosamine) and lipid metabolism. The AA analysis of plasma samples was carried out by the Instrumental Techniques laboratory of the Universidad de Valladolid (UVA, Valladolid, Spain). It was performed using the ZORBAX Eclipse Plus method with an Agilent 1100 HPLC system [
37]. Parameters of glucose and lipid (total cholesterol, high-density lipoprotein cholesterol [HDL-c], low-density lipoprotein cholesterol [LDL-c], and triglycerides) profiles were assessed by a clinical analyzer (Saturno 300 plus, Crony Instruments Srl, Rome, Italy), according to the manufacturer’s instructions. Values for total antioxidant capacity were determined by FRAP (ferric reducing antioxidant power assay; [
38]), while lipids oxidative damage was assessed by MDA (malondialdehyde; [
39]).
2.5. Gene Expression by Quantitative PCR
The following steps described, from the RNA extraction to the raw expression data, were carried out by the Unit of Genomics service of the Universidad Complutense de Madrid (UCM, Spain). First, available duodenum tissue samples (~10 to 20 mg) obtained from the birth sampling (48 samples) and slaughterhouse (52 samples; C: 12 males and 9 females, GLN: 15 males and 16 females) were used for total RNA extraction using RNeasy PowerLyzer Tissue and Cells isolation kit (Qiagen, Hilden, Germany) and the TyssueLyser LT equipment (Qiagen), following the manufacturer’s recommendations. The obtained RNA was quantified using NanoDrop ONE (NanoDrop Technologies, Wilmington, DE, USA), and the RNA quality was assessed with an Agilent bioanalyzer device (Agilent Technologies, Palo Alto, CA, USA). Moreover, 2.2 µg of obtained RNA were treated with RNAse-Free DNAse Set (Qiagen), following the manufacturer’s instructions. Second, the synthesis of cDNA was carried out with High Capacity RNA-to-cDNA kit (Applied Biosystems, Foster City, CA, USA) using 9 µl of RNA treated with con DNAse I (~1 µg of total RNA) in a total volume of 20 µL containing, following the supplier’s instructions. Nine target genes related to mechanistic target of rapamycin (mTOR) pathway (linked to Gln;
EIF4EBP1,
HIFIA,
MTOR,
PPARG,
RPS6,
RPS6KB1,
RPTOR, and
SREBF1 [
40,
41]) and intestinal tight junctions (
OCLN; Occludin) and two reference genes (
ACTB and
GAPDH) were analyzed. Primer pairs were designed using Universal ProbeLibrary Assay Design Center (Roche Life Science, Basel, Switzerland) from the available ENSEMBL sequences and covered different exons to assure the cDNA amplification (
Table S2). Third, the transcript quantification was performed using Power SYBRGreen PCR Master Mix (Applied Biosystems, Warrington, UK) in a QuantStudio 12 K Flex (Applied Biosystems). The qPCR reactions were prepared in a total volume of 10 µL containing 2.5 µL of cDNA (1/10 dilution) and forward and reverse primers (concentration of 300 nM), following the manufacturer’s recommendations. All points and samples were duplicated as technical replicates, and mixes without cDNA were used as negative controls. Cycling conditions were 95 °C for 10 min, followed by 40 cycles of 95 °C (15 s) and 60 °C (1 min). Data were extracted with the QuantStudio 12 K Flex software v1.2.2 (Applied Biosystems).
For statistical analysis of gene expression, the influence of maternal treatment and BIW classification were analyzed with a linear model fitting factors and their interactions as fixed effects with sow (at birth) as a random effect. The method employed for the statistical analysis of the gene expression data [
42] simultaneously analyzed the Cp values for the target and endogenous genes using a linear mixed model. The model used and posterior calculations were performed as previously described by our group [
43]. Treatment and sexes were assessed at both ages, but BIW was only assessed at birth. The adjusted
p-values are indicated as
q-values.
2.6. Fat Content and Fatty Acid Composition of Diets and Tissue Samples
The one-step procedure [
44] was used for the extraction and methylation of the diet fatty acids (FA). Gas chromatography (Hewlett Packard HP-6890, Palo Alto, CA, USA) was used to identify fatty acid methyl esters by a flame ionization detector and a capillary column (HP-Innowax, 30 m × 0.32 mm i.d., and 0.25 µm polyethylene glycol-film thickness; [
45]).
The lipids from the brain, intramuscular fat (IMF) at the LD muscle, and liver fat were extracted [
46]. Fat content was expressed as a percentage (%) of dry matter (DM). Afterward, total lipids at IMF and liver fat were separated into the neutral lipid (NL; in fat storage such as triglycerides) and polar lipid fractions (PL; in cell membranes such as phospholipids; [
47]). Subcutaneous fat was individually analyzed in outer and inner layers. Extracts were methylated [
48] and analyzed using protocols developed at our laboratory [
45]. The individual FA percentages for saturated, monounsaturated, and polyunsaturated FA (SFA, MUFA, and PUFA) were calculated. Total n-3, total n-6 FA, the Σn-6/Σn-3 ratio, and the unsaturated index (UI) were also calculated [
49]. The activity of stearoyl-CoA desaturase enzyme 1 (SCD1) was estimated as C18:1/C18:0 and MUFA/SFA ratios (desaturation indexes; [
50]).
2.7. Statistical Analysis
Data were analyzed by the SAS version 9.4 (Statistical Analysis System Institute Inc., Cary, NC, USA) to assess the maternal treatment effect (C vs. GLN) and its interactions with sex (female vs. male) and BIW (LBIW vs. NBIW). Dependent variables were assessed using two-way ANOVA in a general linear model, including maternal treatment and sex. Except at birth and weaning, when the BIW classification was also added to the model (three-way ANOVA) and all maternal treatment interactions. No significant triple and double interactions were removed from models, and significant interactions were studied individually. Changes over time in body weight, backfat depth, and meat color traits were assessed using a repeated measures test with the Greenhouse–Geisser correction. Reproductive data and birth data were analyzed by parity and maternal treatment, and chi-square was used to assess the percentage of LBIW piglets and percentages of parity. Sow was used as a random effect in the birth and weaning analysis to consider the common maternal environment. Litter size was categorized and used as a random effect for birth data [
10]. For performance parameters, the respective age was used as a covariate. Finally, pig was the experimental unit for all variables studied except for the reproductive data, where sow was the unit. Results were expressed as mean ± SE. Statistical significance was accepted from
p < 0.05 and statistical trend was defined as 0.05 <
p < 0.1.
4. Discussion
The industry of traditional swine breeds is currently increasing its census, production, and economic influence. Batch homogeneity is one of the main production challenges, so the search for useful solutions for its improvement from birth is critical. Thus, testing feed supplements showing positive results in other swine breeds, such as Gln, is required, particularly considering the metabolic differences between traditional and selected swine breeds [
28,
29,
30,
31,
32,
33,
34]. Glutamine is one of the most interesting AA in fetuses and adults because it is one of the most used AA by pig enterocytes as an energy source and in all tissues for nucleic acid synthesis [
18,
19]. The current study is not only important because it is the first one focused on the effect of prenatal Gln supplementation on traditional breeds, but because it is also the first trial of a prenatal supplementation of any AA in this kind of swine breeds.
Among the valuable data of this study is the first assessment of plasma AA levels in newborns of the Iberian breed. Some AA levels (i.e., Gln, Ala, and Asn) were greater in treated newborns. Moreover, the BIW was also studied, and the treated LBIW group showed higher Arg than the corresponding control group. The transfer of Gln to fetal blood is critical for maximizing fetal growth by synthesizing other AA, such as Glu, Asp, and Asn [
51]. In a previous study, the maternal concentrations of Gln in Meishan gilts (Pa1, nulliparous), a highly prolific breed with a remarkable homogeneity in the offspring, was higher than in Large White × Landrace gilts [
52,
53,
54,
55]. The Meishan newborns also had greater Gln, Ala, and Asp concentrations, similarly to Gln supplemented newborns in a previous study, which showed higher Gln values [
27,
55]. However, only the proportional weight of viscerae was higher in our treated newborns related to prenatal body growth. Glutamine is recognized as a nutritionally essential AA for gestating, especially after the early gestation period and in gilts, so more improvements in litters from gilts were expected [
15,
56]. On the contrary, the lack of effects on the development of litters from sows supplemented with Gln has been previously described in selected swine breeds [
57,
58]. However, the promising potential of the Gln maternal supplementation to reduce negative effects of the IUGR process, according to scientific literature [
15,
19,
27,
59,
60], led to expect more effects than those found under our conditions. Nevertheless, future studies with larger sample size, particularly gilts, could increase the detection power of productive effects, based on the result of previous studies carried out in other swine breeds.
Regarding the FA profile of different tissues, there were significant effects of the maternal treatment at birth. First, the control group showed differences between LBIW and NBIW newborns, but not in the treated group. Thus, the FA composition at birth would be more homogeneous in newborns from pregnancies supplemented with Gln. Second, treated newborns showed greater amounts of Σn-3 FA in brain and muscle cellular membranes, related to a protective effect because of the improvement of pro-/anti-inflammatory status and the reduction of pathological risks [
61,
62]. Finally, the treated group had a higher amount of C18:1 FA in the liver, which was previously found in Meishan newborns [
54]. This result could be associated with a possible adaptation to increase the survival capacity [
54], although it was not possible to corroborate it with our results.
On the other hand, the supplementation of Gln also activates, directly and indirectly, the protein synthesis in fetal skeletal muscle and the mTOR signaling pathway, linked to cell growth and proliferation [
19,
23,
40,
63]. Our study has shown higher expression in components of the mTOR protein complex 1 (mTORC1;
MTOR and
RPTOR [raptor]) in treated NBIW newborns than in their counterparts of the control group. The mTORC1 is related to the regulation of protein and lipid synthesis, autophagy, and energy metabolism, but none of its downstream genes studied showed differences between treatment groups [
40,
41,
63]. This finding would consist of the lack of evidence of greater prenatal growth in supplemented pregnancies. Furthermore, differences in the
MTOR expression between LBIW and NBIW of the control group, but not in the treated group, were also found. Thus, treated newborns showed, as in the FA profile, more homogeneity than the control group. Finally, the expression related to intestinal tight junctions was also not different by treatment, although a previous study found benefits from the Gln supplementation in LBIW and NBIW [
27]. So, physiological mechanisms related to Gln could improve prenatal development that would improve productive parameters at birth. However, 1% Gln supplementation in Iberian gilts and sows under our conditions seems not to be enough to trigger it.
The current study is the first to show the offspring development after weaning in pigs from pregnancies supplemented with Gln. Although it is always a challenge to find medium or long-term effects of maternal supplementations on the offspring, this valuable data allows us to assess the possibility of physiological or cellular changes with postnatal effects due to the prenatal Gln supplementation. However, no beneficial effects of the maternal treatment were found at weaning, at 215 days-old, or at the slaughterhouse under our conditions. Furthermore, some differences in backfat depth, carcass traits, and tissue FA profiles were affected by sex, mainly in females, although no sex-related effect was found at birth. Nevertheless, sex is a well-known important factor in swine production and, especially, in the Iberian breed [
10,
64,
65,
66].