1. Introduction
The first days postnatal are a critical period of metabolic-nutritional programming in pigs. Of particular interest to this study is the relationship between early nutritional environment and future lactation performance. In swine, greater colostrum intake by gilts was related to earlier puberty and better lactation performance (as sows) than the low colostrum intake counterparts [
1]. Studies in sheep showed the nutrition of ewes during pregnancy affected the yield and composition of milk produced by the offspring [
2,
3]. The preweaning growth rate of heifers was positively related to their milk production as cows [
4]. Heifer dairy calves fed two liters of colostrum produced less milk in their first and second lactations than calves who were fed four liters [
5]. The lower milk production was related to greater rates of morbidity and lower body weights of calves. Heifer calves fed restricted versus ad libitum intake of milk replacer had less mammary gland mass, mammary parenchyma, fat pad mass, and lower expansion of epithelium into the adjacent stromal tissue [
6].
The number of mammary epithelial cells is highly correlated to milk production [
7]. Mammary cell number is established during the development of the gland, which begins in utero [
8]. The development of mammary glands, from birth to the peripubertal period, in pigs is characterized by ductal elongation and formation of a lumen [
9]. The pattern of ductal branching in swine is similar to the human breast, which is characterized by terminal ductal lobular units (TDLU). Prior to the onset of puberty, mammary branching and organization increases in complexity from primarily TDLU-1 to TDLU-2 through the proliferation of lobular buds and ductal elongation.
Knowing that nutrition in early life impacts future milk production and that milk production is determined by the number of mammary epithelial cells led us to hypothesize that the level of colostrum intake during the first 24 h postnatal affects the proliferation rate of mammary epithelial cells in the gland, which in the long term would affect the number of milk-producing epithelial cells. Stem and progenitor cells populate tissues by asymmetric cell division [
10]. The proliferation and self-renewal of stem-progenitor cells are balanced in the tissue by daughter cells exiting the cell cycle and beginning the differentiation process. As cells differentiate, they become more specialized, and this specialization is marked by protein synthesis [
11]. In the gland, at any point in time, there is an array of cellular states as cells progress towards differentiated states [
12]. Accounting for these parameters, we posited that the ratio of newly synthesized protein to newly synthesized DNA could be used as an indicator of the relative amount of cellular differentiation to cell division.
The overall goal of this study was to test this hypothesis and determine the relationship between the level of colostrum intake and 24 h level of circulating amino acid, glucose and insulin with mammary parenchyma histological features, as well as cell division and protein synthesis over the first postnatal week. For the described studies, one of two standardized doses of a homogenate colostrum sample, 10% and 20% of birth bodyweight, was fed to study animals [
13]. Gilts were given a bolus of heavy water (deuterium oxide) immediately after birth and daily for seven days to label newly synthesized DNA and protein over this time period. DNA and protein were isolated from mammary parenchymal tissue. DNA and protein fractional synthesis (f) and fractional synthetic rate (fraction per day, FSR) were calculated using mass isotopomer distribution analysis (MIDA) [
14,
15,
16,
17]. The ratio of protein f and FSR to DNA f and FSR were calculated and used to indicate the relative amount of differentiation to cell division. Mammary morphological development was also analyzed by measuring the epithelial area of parenchymal tissue and the stromal and epithelial proliferation index on postnatal day seven.
2. Materials and Methods
2.1. Animals and Study Design
Animals and Study Design (
Figure 1). Prior to the start of this study, all animal procedures were reviewed and approved by the Institutional Animal Care and Use Committee (Protocol # 1907001920). Animals used for the study were born at the Purdue University Animal Sciences Research and Education Center Swine Farm to eight York (×) Landrace multiparous (3.25 ± 1.67 parities) sows bred to terminal sire Duroc boars. Sixteen animals (
n = 16) used for the investigations were part of a larger study [
13]. Animals were selected at birth from across eight different litters (
n = 2 gilts/litter; birth litter size was 12.3 ± 1.3 live piglets/sow). At birth, gilts were towel-dried, weighed and assigned to one of two study treatments (1 gilt/litter/treatment); colostrum fed at a rate of 24 h intake of 10% of birth body weight (BW; COL10) or 20% of BW (COL20). Gilts that were selected weighed between 1.2 and 1.8 kg and treatments were matched across litters by weight. All gilts were bottle-fed with a pooled colostrum sample every 2 h from birth until 24 h of age. Between feedings, they were returned to a nursery area, which was maintained at 40 °C. Nursery temperature was chosen to prevent any risk of hyperthermia at this age, as this temperature is the same as the offsprings’ in utero environment.
A blood sample was collected from each gilt at 24 h postnatal, via jugular venipuncture, using a 22 gauge × 2.5 cm needle into a 2 mL potassium-EDTA coated vacutainer tube (BD367841, BD, Franklin Lakes, NJ, USA). At 24 h, body weight, body temperature and blood samples for plasma preparation were collected from each gilt to measure blood glucose, insulin, total protein and amino acid. Plasma was prepared by immediately centrifuging blood at 2000× g for 15 min (E8 Centrifuge, LW Scientific Inc., Lawrenceville, GA, USA).
Gilts were returned to their birth dam to be nursed for the remaining days of the study. Nursing litter size was standardized to 12–14 piglets per sow. Body weights were recorded daily for each gilt. Gilts in the COL10 and COL20 groups were euthanized on day seven postnatal. Animals were euthanized using CO2 inhalation, and then mammary tissue was collected by doing an incision longitudinally along both sides of the left mammary chain. Additionally, six gilts weighing between 1.2 and 1.8 kg were identified at birth, immediately euthanized and used to evaluate mammary morphology at birth.
The whole mammary chain was removed by dissecting through the subcutaneous tissue. For histology, a square of 1.5 × 1.5 cm of skin around the nipple was sectioned and placed in 10% buffered formalin. After 24 h, histology samples were transferred to 1X PBS. Mammary tissue used for MIDA analysis of DNA and protein synthesis was collected from the thoracic mammary glands, and mammary parenchyma was grossly dissected to remove muscle, skin and other tissue not associated with gland parenchyma, snap-frozen in liquid nitrogen and stored at −80 °C.
A sample size power analysis was performed prior to the start of the study. The power of the study with six animals per treatment with an alpha-error of 0.05, 1.5-fold difference between treatments and 0.25 standard deviation was 0.95. If the difference dropped to 1.4-fold, the power of the study was 0.8 with six animals. Since we anticipated the potential of loss of piglets, the study was begun with eight animals per treatment. Following tissue and plasma collection, all researchers were blinded to treatment during the experimental analysis portion of the study. The treatment groups were revealed for data compilation and statistical analysis on the effect of treatment.
2.2. Colostrum Sample and Analysis
Approximately 50 mL of colostrum was collected from multiple sows (~250) over the course of 7 mo. Colostrum collection was done manually during active farrowing when oxytocin levels are naturally high. Following collection, colostrum was frozen and stored at −80 °C until the day prior to the start of the study. A homogenate-pooled sample was prepared following overnight thawing of colostrum at 4 °C. Piglets were fed this homogenate sample, and several aliquots were collected and stored at −80 °C for subsequent composition analysis.
Colostrum composition was analyzed for percent fat, protein, and insulin concentration. Fat percentage was determined using the creamatocrit approach by centrifuging homogenate samples at 12,000× g for 10 min in a non-heparinized hematocrit tube (3 tubes per sample). Fat percentage was calculated as the ratio of the length of fat to total sample length measured with a caliper and then multiplied by 100.
The protein content of colostrum samples was measured using the Bradford Assay Kit (Pierce Coomassie Plus Assay Kit, Thermo Fisher Scientific, Waltham, MA, USA). Samples were diluted at 1:100 in phosphate buffer manufacturer’s instructions were followed. A plate spectrophotometer (Sparks 10M multimode microplate reader, Tecan) was used to analyze absorbance at 495 nm wavelength.
Colostrum composition was analyzed for percent fat, protein and insulin concentration. Fat percentage was determined using the creamatocrit approach by centrifuging homogenate samples at 12,000× g for 10 min in a non-heparinized hematocrit tube (3 tubes per sample). Fat percentage was calculated as the ratio of the length of fat to total sample length measured with a caliper then multiplied by 100. The protein content of colostrum samples was measured using a Bradford Assay Kit (Pierce Coomassie Plus Assay Kit, Thermo Fisher Scientific; Waltham, MA, USA). Samples were diluted at 1:100 in phosphate buffer, and the manufacturer’s instructions were followed. A plate spectrophotometer (Sparks 10M multimode microplate reader, Tecan Trading AG, Mannedorf, Switzerland) was used to analyze absorbance at 495 nm wavelength.
Colostrum insulin was analyzed in duplicate samples using a porcine insulin ELISA kit (cat no. 10-1200-01; Mercodia AB; Winston Salem, NC, USA). Insulin was measured in both homogenate and skimmed colostrum samples. Intraplate variation was 4.75%.
2.3. Neonate Plasma
2.3.1. Protein
Plasma protein was measured in duplicate using the Bradford Assay (Bio-Rad Laboratories, Inc., Hercules, CA, USA) following manufacturer instructions. Prior to analysis, plasma was diluted 1:100 with phosphate-buffered saline. Intraplate CV was 3.65%.
2.3.2. Insulin
Plasma insulin was analyzed in duplicate samples using a porcine insulin ELISA kit (cat no. 10-1200-01; Mercodia AB; Winston Salem, NC, USA), following manufacturer instructions. Intraplate variation was 4.75%.
2.3.3. Glucose
Plasma glucose was determined using Autokit Glucose (Fujifilm Wako Diagnostics USA Corporation, Mountain View, CA, USA) following manufacturer instructions. Intraplate CV was 4.84%.
2.3.4. Free Amino Acids
Free amino acid content of neonate plasma was analyzed using liquid chromatography-tandem mass spectrometry (LC/MS-MS) in Purdue University’s Bindley Biosciences Metabolite Profiling Facility. Briefly, 10 μL of amino-butyric acid at a concentration of 1 μg/uL and 25 μL of 100% trichloroacetic acid (TCA) solution were added to 100 µL of plasma. Samples were incubated for 10 min at 4 °C followed by centrifugation at 14,000× g for 10 min. The supernatant was collected and stored at −20 °C until analysis. Just prior to liquid chromatography, 100 µL of acetonitrile (ACN) was mixed with 100 µL of supernatant. Liquid chromatography was performed using Intrada Amino Acid 3 μm, 2 × 150 mm column (Imtrakt USA, Portland, OR, USA) connected to an Agilent 6470 QQQ LC-MS/MS system (Agilent, Santa Clara, CA, USA). Acetonitrile with 0.3% of formic acid and acetonitrile with 100 mM ammonium formate solution (20:80 v/v) were used as mobile phases.
2.4. Histological Analysis of Mammary Gland Development
All tissue preparations for histological analysis were done by the Purdue University Histology Research Laboratory. Mammary tissues were fixed in 10% neutral buffered formalin for 24 h and transferred to PBS until processing for paraffin embedding. Paraffin processing was done in a Sakura Tissue-Tek VIP6 tissue processor for dehydration through graded ethanols, clearing in xylene and infiltration with Leica Paraplast Plus paraffin. After processing, tissues were embedded in Leica Paraplast Plus paraffin. Tissue sections were taken at a thickness of 4 µm using a Thermo HM355S microtome. Sections were mounted on charged slides and dried for 30–60 min in a 60 °C oven. After drying, all slides were deparaffinized through 3 changes of xylene and rehydrated through graded ethanols to water in a Leica Autostainer XL. For hematoxylin and eosin (H&E) staining of tissues, the Leica Autostainer XL was used. Tissue sections were stained in Gill’s II hematoxylin, blued and counterstained in an eosin/phloxine B mixture. Finally, tissues were dehydrated, cleared in xylene and cover-slipped in a toluene-based mounting media (Leica MM24).
H&E-stained tissues were used to measure the proportion of epithelial tissue within the parenchymal compartment. First, ImagePro Plus 5.1 (Media Cybernetics) was used to capture histological images in conjunction with a Nikon Eclipse 50i microscope (Nikon Inc., New York, NY, USA; Evolution MP, Media Cybernetics Inc., Rockville, MD, USA). Multiple images of H&E stained tissue were captured at 10× magnification to encompass the entire parenchymal area of the gland for each animal. The parenchymal area was defined for this study as the epithelial cells of the terminal ductal lobular units (TDLU) and associated ducts along with intralobular and interlobular stroma. To create a panorama of the entire parenchymal area of the cross-section, images were merged into a single image using Adobe Photoshop (V 22.1.0, Adobe). ImageJ was used to measure the area in the tissue section (
Figure 2). The “Draw/Merge: Trace” tool was used to first select parenchymal tissue and calculate the area, then to trace and calculate the entire epithelial area of TDLU (epithelium plus lumen) and finally to trace around the lumen and calculate that area. The ratio of epithelium within parenchyma was calculated by subtracting the lumen from the epithelial area of the TDLU and then dividing this by parenchyma area, and this was defined as parenchymal epithelial area (PEA).
Tissue sections were also immunostained with KI67 to mark proliferating populations of cells. After deparaffinization, antigen retrieval was done with a TRIS/EDTA pH 9.0 solution in a BioCare decloaking chamber (Pacheco, CA, USA) at a temperature of 95 °C for 20 min. Slides were cooled for 20 min at room temperature and transferred to TRIS buffer with Tween 20 detergent (TBST). The rest of the staining was carried out at room temperature using a BioCare Intellipath stainer. Slides were incubated with 3% hydrogen peroxide in water for 5 min. Slides were rinsed with TBST and incubated in 2.5% normal goat serum for 20 min. Excess reagent was blown off, and Ki67 primary antibody (Cell Marque, 275R-16, Rocklin, CA, USA) was applied at a dilution of 1:100 (0.364ug/mL) for 30 min. The negative control slide was stained with Rabbit IgG (Vector Labs, I-1000, Burlingame, CA, USA) at a concentration of 1:5000 (1 µg/mL) for 30 min. Slides were rinsed twice in TBST, and a goat anti-rabbit secondary antibody (Vector Labs, MP-7451) was applied for 30 min. Slides were rinsed twice in TBST, and Vector ImmPACT DAB (Vector Labs, SK-4105) was applied for 5 min. Slides were rinsed in water and transferred to a Leica Autostainer XL(Wetzlar, Germany)or hematoxylin counterstain, dehydration and cover-slipping.
Five images per gilt were taken at 200× magnification. Sections of jejunum tissue were used as a positive control for the specificity of KI67 staining for proliferating populations of cells. To determine the proliferation index of mammary epithelial cells and proliferating intralobular stroma cells in parenchymal tissue, an ImageJ plugin called Cell Count by GNU General Public License was utilized. As above, the parenchymal area was defined for this study as the epithelial cells of the TDLU with ducts and associated intralobular and interlobular stroma. The proliferation index of epithelial and stromal cells within parenchymal tissue was determined. All epithelial cells were positively stained for KI67, the five sections were counted, and the epithelial cells without staining were counted. Similarly, intralobular and interlobular stromal cells that were immunostained for KI67 were counted, and all cells not stained were counted. The total number of each cell type was determined, and then the number of proliferating epithelial or stromal cells was divided by the total of each type to determine the percent of proliferating cells.
All research assistants that analyzed histomorphology were blinded to treatment and day and trained by one individual on the approach to conducting analyses. For each animal, three research assistants analyzed histomorphic features, and data across the three researchers were averaged for final counts. Tissue was available for all animals that survived to postnatal day 7 of COL10 (n = 7). However, the quality of tissue collected for one COL20 piglet was not representative of parenchyma, so only six animals in this treatment were used for histological analysis. To determine relative changes in the proliferating index of epithelial and intralobular stromal cells in the parenchyma and parenchymal epithelial area between birth and postnatal day 7, tissue from the baseline group (n = 6) of gilts was also analyzed.
2.5. Mass Isotopomer Distribution Analysis (MIDA) of DNA and Protein Synthesis (f) and Fractional Synthetic Rate
2.5.1. Metabolic Labelling with Deuterium Oxide (D2O)
In order to obtain good quality data for mass isotopomer distribution analysis for proteome studies, greater than 2.5% steady-state enrichment of heavy water-deuterium oxide (D
2O) is recommended [
14]. To achieve this, the doses of D
2O used were adapted from Lam et al. [
17]. Metabolic labeling of newly synthesized deoxyribose and amino acid molecules was begun immediately after birth by administering a deuterium oxide bolus to animals via intraperitoneal injection (IP) of 0.9% NaCl in D
2O (Millipore Sigma, Burlington, MA, USA), at 20 mL/kg of BW. A second bolus was given to animals 4 h later, using the same dose and route of administration. At 24 h postnatal, and then daily until postnatal day 7, piglets were orally gavaged with 10 mL/kg of BW deuterium oxide at 0600 h to maintain enrichment. Six unlabeled animals were administered equivalent doses of saline. Tissue and blood of the unlabeled group were used to determine the baseline mass isotopomer distribution. At birth and postnatal days 1, 3, 5 and 7, blood samples were collected by jugular puncture and plasma was isolated for analysis of D
2O enrichment.
2.5.2. Determination of D2O Enrichment in Plasma Using Gas Chromatography-Tandem Mass Spectrometry (GC-MS/MS)
To determine the percent of body water that was D
2O, the approach of base-catalyzed exchange of hydrogen (deuterium) between water and acetone was used [
18], with modifications by Purdue University’s Bindley Biosciences Metabolite Profiling Facility. Briefly, piglet plasma (20 μL) was mixed with 2 μL of 10 N NaOH and 4 μL of 5% (
v/
v) acetone in acetonitrile. The sample and standard curve (ranged from 0.5 to 16% of D
2O in H
2O) mixtures were incubated at room temperature overnight, after which the acetone portion was extracted by the introduction of 500 μL of chloroform and 0.5 g of anhydrous sodium sulfate (Na
2SO4). Samples were centrifuged at 14,000×
g for 1 min to precipitate the Na
2SO4. The chloroform solution was transferred to a vial, and acetone was measured using GC/MS in Purdue University’s Metabolomics Core. Gas chromatography was carried out using an Agilent Select FAME column (CP7419, Agilent, Santa Clara, CA, USA) attached to a TSQ 8000 triple quadrupole GC-MS/MS (Thermo Fisher Scientific, Waltham, MA). The intensity of acetone was measured at 58 and 59
m/
z and was used to calculate D
2O percentage in plasma.
2.5.3. Isolation of DNA from Tissue and DNA Hydrolysis
The procedures for DNA isolation and hydrolysis described were modified from the approach described by others [
19] in the following way. Approximately 25 mg of mammary tissue from the parenchymal area was used for DNA isolation. Tissue was ground using a tissue homogenizer in DNA extraction buffer from the gMAX Mini Genomic DNA Kit (IBI Scientific, Dubuque, IA, USA). DNA was eluted from the column using 100 µL of TE buffer (tris-HCl + EDTA), and the concentration was measured with a nanodrop system. DNA was diluted with hydrolysis buffer (20 mM Tris-HCl, 100 mM NaCl, 20 mM MgCl
2, pH 7.9) to 1 µg in a total volume of 50 µL, and 50 µL of hydrolysis enzyme cocktail [benzonase (E-1014); phosphodiesterase I (P-3243) and alkaline phosphatase (P-7923); MilliporeSigma, Burlington, MA, USA) was added. The hydrolysis reaction was carried at 37 °C for 6 h. Then, samples were dried overnight and stored at −20 °C until analysis.
2.5.4. Liquid Chromatography-Tandem Mass Spectrometry (LC-MS/MS) Analysis of Adenosine Isotopomer Distribution
D2O labels the deoxyribose moiety of dNTPs in replicating DNA through the de novo nucleotide synthesis pathway. The isotopic enrichment of the purine deoxyribonucleoside adenosine is then determined by LC-MS/MS. Briefly, samples were reconstituted in 100 µL of 5% MeOH/95% 5 mM ammonium formate. Molecule separation was carried out with 5 mM ammonium fumarate and 100% methanol as mobile phases in a Waters Atlantis T3, 3 μm, 2.1 × 50 mm column (186003717, Waters Corp., Milford, MA, USA) connected to an Agilent 6470 QQQ LC-MS/MS system (Agilent, Santa Clara, CA, USA). Multiple reaction monitoring (MRM) of the ribose portion of adenosine (dA) was measured based on the parental and product ions 251 → 117 m/z (M0). Ion combinations for M+1 and M+2 were identified and measured based on the identifications of 252 → 118 m/z and 253 → 119 m/z, respectively.
2.5.5. Protein Hydrolysis
Preparation of protein hydrolysate for measuring global protein synthesis was done as described [
15] with some modifications. Briefly, approximately 25 mg of parenchymal mammary tissue were placed in a 5 mL amber glass vial (Fisherbrand, Thermo Fisher Scientific, Waltham, MA, USA), and 1 mL of 6 M HCl was added under the fume hood. Samples were homogenized using the Fisherbrand 150 handheld tissue homogenizer (Thermo Fisher Scientific, Waltham, MA). The probe of the homogenizer was washed with sterile water between samples. Caps were placed in vials and incubated at 120 °C in a forced air oven (Model 414004-576, VWR International, West Chester PA, USA) for 24 h. Following incubation, samples were transferred to a 1.5 mL tube and centrifuged at 14,000×
g for 10 min. The supernatant was transferred to a 1.5 mL tube and dried in a savant SPD 2010 speedvac concentrator (Waltham, MA, USA) overnight. The dried samples were stored at −20 °C until amino acid extraction.
2.5.6. Amino Acid Extraction LC/MS Analysis of Isotopomer Distribution of Alanine
Dried protein hydrolysates were reconstituted by adding 300 µL of PBS and vortexing the samples, and 100 µL was transferred to a new 1.5 mL tube. Twenty-five µL of TCA (trichloroacetic acid, saturated solution, 1000 mg of TCA + 700 µL H2O) was added and samples vortexed to mix. Samples were then centrifuged at 14,000× g for 10 min, and 50 µL were transferred to a new tube, being careful to avoid black precipitate. Then 50 µL of acetonitrile was added, and samples were mixed well by vortexing. One hundred µL of this extract was used for LC/MS analysis of alanine.
The method used to determine the isotopomers of alanine was developed by Purdue University’s Metabolite Profiling Facility, Bindley Bioscience Center, through modification of the methods used to measure amino acids. In this method, an Intrada Amino Acid column was used for the liquid chromatography (LC), followed by a quadrupole mass spectrometer (MS). Alanine is retained to ~11.5 min of the run, and the mass spectrometry returns a precursor ion of 90
m/
z and a product ion of 44
m/
z. The fragment of 44
m/
z (with chemical formula C
2H
6N) contains four hydrogens that can potentially be replaced by deuterium during the synthesis process. The precursor (alanine, C
3H
7NO
2) and product (C
2H
6N) will increase mass equally as deuterium is added to the molecule. For this method, the LC/MS machine and software is programmed to measure the intensity/area of the peaks of molecules with precursor → product ion pair of 90
m/
z → 44
m/
z; 91
m/
z → 45
m/
z; 92
m/
z → 46
m/
z; 93
m/
z → 47
m/
z; and 94
m/
z → 48
m/
z; in order to measure the intensity/area of isotopomer (M) with no heavy isotopes (M0), one (M+1), two (M+2), three (M+3) and four (M+4), respectively.
Supplementary Materials Information S1 shows the distribution of alanine M0, M+1, etc., in a sample from an unlabeled animal (blank) and eight samples from D
2O labeled animals, with corresponding LC/MS spectra of samples.
2.5.7. Mass Isotopomer Distribution Analysis (MIDA) of Adenosine and Alanine for Calculation of the Fraction (f) of DNA and Protein Newly Synthesize and the Fractional Synthetic Rate (FSR) of DNA and Protein
To determine the percent of newly synthesized DNA and proteins, the percent of enrichment (
p) was calculated as the mean percentage of D
2O enrichment from postnatal day one to day seven in each piglet. LC-MS/MS analysis of adenosine and alanine provided the intensities for M0 to M+2 and M0 to M+4, respectively. The percentage of M0 (%M0) at birth and day seven (%M0
7) was calculated as the intensity M0 isotopomer divided by the total of the intensities for all isotopomers measured for adenosine and alanine. The difference between %M0 and %M0t is defined as EM0t [
14].
EM0t can be defined as the modifications on %M0 after a specific time (t) of exposure to D
2O and is, in part, determined by the synthetic ratio (balance between synthesis and degradation) of DNA, as tissue cells are in a constant turnover of division and death. Each cell has a different turnover based on its biological functions and properties. Moreover, the speed of synthesis of dA, or FSR, can be modified by physiological events or experimental treatments. The %M
0* (max %M
0 when 100% of the DNA are de novo synthesized) can be calculated based on %M
0,
p and
n. Where
p is the probability of being labeled, which is the percent enrichment of D2O, and
n is the number of hydrogens in a molecule that can be labeled (‘replace’) hydrogen. The following formulas for DNA synthetic rate were adapted from [
16,
17]:
After we know the %M
0*, we are able to elucidate EM
0*, similar to EM
0t.
EM0* is the %M0 when 100% of the DNA is turned over.
Once EM
0t and EM
0* are determined, fractional synthesis (
f) is calculated.
2.6. Statistical Analysis
All statistical analyses were performed in SAS (version 9.4; Cary, NC, USA). Amino acids that were below the detectable limit were assigned a value that was one-fifth the highest value, which was below the lowest value. Procs Mixed model was used with treatment run as the class. Normality was checked by running the residuals in the model. All statistical models included treatment and birthweight. Birthweight was excluded from the model if not significant (p < 0.05). The least-square means were performed using the Tukey-Kramer method, with a pairwise comparison. p < 0.05 was considered significant, with p-value > 0.05, but ≤ 0.1 discussed as a tendency. Bodyweight and average daily gain were run with Procs Mixed with treatment and day as fixed effects with the day as a repeated measurement. Treatment by day interactions and birthweights were included in the model and removed if not significant (p < 0.05). Proc corr was used for correlations.
4. Discussion
The data collected supports the relationship between factors indicative of perinatal nutritional environment and mammary growth and development over the first week postnatal. In particular, plasma lysine level at 24 h postnatal was positively related to average daily gain, the fraction of newly synthesized DNA (f) in mammary parenchymal tissue over the first seven days postnatal, and the fractional synthetic rate of DNA in mammary parenchyma. Plasma lysine was also inversely related to the ratio of protein to DNA f and FSR. This relationship, as posited in the introduction, may reflect that higher lysine levels favored a greater degree of cell division versus cells leaving the cell cycle and differentiating. The relationships between nutritional environment and mammary development were found despite the fact that colostrum dose was not related to any of the variables used to evaluate mammary development. There may not be an effect of colostrum dose on variables measured. In light of this possibility, it is interesting to note that the amount of DNA isolated per unit of mammary parenchymal tissue was numerically higher in COL20 versus COL10 animals. This finding suggests that the level of colostrum intake may affect the number of cells in parenchyma. Analysis of DNA content at an earlier time point is needed to determine this. Moreover, future studies using tools like single-cell RNA-seq would help in understanding whether the amount of colostrum consumed affects the developmental program of subpopulations of cells within the gland. The lack of an effect may also have been related to the study design. Different doses of colostrum resulted in COL20 animals weighing significantly more after the 24 h of colostrum feeding, and these differences were maintained to postnatal day seven [
13]. However, returning piglets to birth litters likely had unmeasured impacts on perinatal nutrition. One piglet in each group died by crushing, and the growth rates were highly variable after return to litters. Bottle feeding and returning neonates to litters where they competed for access to milk, likely differentially stressed animals and contributed to piglets’ nutritional environment. Moreover, sow milk quality likely varied across litters. With this in mind, circulating levels of nutrients and gross measures of growth were used as corollaries of the postnatal nutritional environment.
Variables used to assess mammary development over the first week postnatal correlated with each other, supporting the potential for underlying relationships. PEA reflects the relative amount of epithelial tissue in the gland and was positively correlated with the percent of epithelial and stromal cells proliferating on postnatal day seven, as well as the average daily gain. Percent of proliferating epithelial cells was also related to protein synthesis over the first week postnatal. The different degrees of mammary development between animals in our study were most likely attributed to the postnatal environment and not prenatal factors, such as mammary histomorphology and DNA and protein f and FSR variables were not correlated with gilt birthweight. Whereas PEA and f and FSR of mammary DNA and protein were correlated with average daily gain. Moreover, crown-rump length on postnatal day seven was correlated with the fraction of newly synthesized DNA and DNA FSR. Although it is important to note that the nutritional environment in late gestation carries over to postnatal growth performance, as discussed below, in regard to lysine. Since gilts were matched by birthweight when assigned to treatments, the relationship to these variables would be minimally expected.
The relationship between average daily gain and mammary growth metrics supports mammary growth as isometric to body growth in the first postnatal week. The relationship between average daily gain to the mammary parenchymal epithelial area is particularly intriguing with regard to our hypothesis. This relationship likely indicates that the more adequately nourished the gilt is, the more resources are partitioned to growing the secretory component of the gland. A higher proliferative index of epithelial cells was found in mammary tissue isolated at birth compared to tissue from COL10 and COL20 animals on postnatal day seven. The proliferative population at birth may be particularly sensitive to nutritional growth cues in the gilt’s postnatal environment. Studies of mature pigs support that the nutritional environment affects mammary development [
20,
21].
As an essential amino acid, lysine availability is limiting to porcine growth. Studies of the effect of lysine availability on suckling neonates through maternal milk found diets with a 20% deficiency in lysine content had reduced litter growth by 8–10% [
20]. At the same time, litter growth increased by 2.35 times with a 2.90-fold increase in amino acids in the sow’s diet. Improved growth was directly related to the increased intake of lysine and other amino acids by neonates [
20]. Thus, our finding that plasma lysine levels at 24 h were related to the average daily gain of neonates is consistent with its availability to neonates limiting growth. Currently, it is not known what led to the varying concentrations of circulating lysine across gilts, as levels of lysine were not related to the 24 h dose of colostrum. Increasing lysine content of sow diet in late gestation increased the total number born alive and birth weight of piglets [
22]. Similarly, increasing the lysine and fat content of sow diets in late gestation diets improved the overall performance of litters [
23]. It also increased sow colostrum production. Although all sows in the present study were on the same diet, the efficiency of absorption of lysine by sows and placental lysine transfer to gilts during their late fetal growth may be a factor. Further research in this area is needed, as lysine levels at 24 h postnatal were also positively related to the fraction of newly synthesized DNA and the fractional synthetic rate of DNA in mammary parenchymal tissue over the first seven days postnatal. At the level of the cell, nutrients, including amino acids, regulate gene expression [
24]. The mTORC1 signaling pathway in cells functions to integrate nutrient availability, growth factor signaling and developmental cues to regulate growth [
25]. The production of proteins, lipids and nucleotides need to increase for cells to grow and divide, while catabolic pathways such as autophagy need to be suppressed. mTORC1 regulates all of these processes [
25]. Amino acids and positive cellular energy status activate mTORC1, whereas lack of amino acids or energy inhibits its activity [
26]. Thus, the association of lysine with the total amount of DNA synthesized and DNA synthetic fractional synthetic rate likely reflects that it is indicative of nutrient-amino acid availability that can be used for cellular growth, and in particular, nucleotide synthesis.
The body does not store amino acids, so muscle protein is mobilized to produce free amino acids [
24]. Protein undernutrition decreases the plasma level of most essential amino acids and causes adjustments in physiological functions, with a primary consequence of feeding a low protein diet resulting in the inhibition of growth [
24]. The lower levels of essential amino acids and lower growth of COL10 gilts relative to COL20 animals likely reflect adaptations of COL10 animals to the undernutrition they experienced over the first 24 h of the experiment. Furthermore, lower in COL10 animals, there were several non-essential amino acids, including glutamine. For maximal growth performance, pigs require dietary glutamine [
27]. Studies of neonatal pigs found that nearly all glutamate and glutamine feed was metabolized by the gut, so that glutamate and glutamine in the body must derive almost entirely from synthesis de novo [
28]. Glutamine and glutamate are precursors and products of each other, with reactions catalyzed by glutamine synthetase and glutaminase. Circulating glutamate levels appear to be resistant to large variations, with levels relatively constant in experimental manipulations of diet and energy [
29]. Part of this persistence may be due to the sensitivity of hepatic glutaminase to metabolic state and its transcriptional stimulation induced by starvation. Meanwhile, there is little effect of metabolic state on glutamine synthetase [
29]. This may explain the significantly lower levels of glutamine in COL10 versus COL20 animals. Moreover, although seemingly counterintuitive, the negative relationship between glutamate and PEA may reflect that lower energy in the postnatal environment negatively impacts mammary epithelial expansion.
In the present study, the effect of returning piglets to the birth litter after 24 h of bottle feeding was not controlled, other than by matching treatments by litter. Differences in milk composition between sows and competition between piglets could influence developmental trajectory in the mammary tissue. In addition, this study did not consider other bioactive factors in colostrum, like fatty acids or hormones. Future studies aiming to understand the programming effects of colostrum need to control for maternal effects on piglet development after colostrum feeding by returning piglets used in the study to a common sow. This approach would limit competition of study piglets with established piglets that have experience nursing from the dam. Moreover, to control for nourishment versus bioactive factors in milk, future studies should use a nutrient-matched formula that represents the energy provided by colostrum but does not contain bioactive components. Finally, future studies should include more time points for sample and tissue collection to increase the understanding of the mechanisms involved in neonatal programing.