In utero hyperthermia in late gestation derails dairy calf early-life mammary development

Prenatal hyperthermia has immediate and long-term consequences on dairy cattle growth, immunity, and productivity. While changes in the molecular architecture are reported in the mature mammary gland (MG), any influence on early-life mammary development is unknown. Herein, we characterize the impact of late-gestation in utero heat stress on heifer mammary gross and cellular morphology at early-life developmental stages (i.e., birth and weaning). During summer, pregnant dams were exposed to environmental heat stress (shade of a free-stall barn) or offered active cooling (shade, fans, and water soakers) for 54 ± 5 d before parturition (avg. temperature-humidity index = 79). Heifer calves born to these dams were either in utero heat-stressed (IU-HT; n = 36) or in utero cooled (IU-CL; n = 37) and were managed as a single cohort thereafter. A subset of heifers was euthanized at birth (d0; n = 8/treatment; 4.6 ± 2.3 h after birth) and after weaning (d63; n = 8/treatment; 63.0 ± 1.5 d) to harvest the whole MG. An ultrasound of rear mammary parenchyma (MPAR) was taken prior to d63 and correlated to harvested MPAR cross-sectional area and weight. Portions of mammary fat pad (MFP) and MPAR were preserved for compositional and histological analysis, including ductal structure number and cross-sectional area, connective tissue area, and adipocyte number and cross-sectional area. Cellular proliferation in MPAR was assessed via Ki-67 immunohistochemistry. Relative to IU-CL heifers, the MGs of IU-HT heifers were shorter in length at d0 and d63 (P ≤ 0.02). There were moderate correlations between d63 ultrasound and harvest measures. The IU-HT heifers had reduced MG and MFP mass at d0 and d63 (P ≤ 0.05), whereas MPAR mass was reduced only at d0 (P = 0.01). IU-HT heifers had greater MPAR protein and DNA content at d63 (P ≤ 0.04), but there were no MFP compositional differences (P ≥ 0.12). At d0, IU-HT heifers had fewer MPAR ductal structures (P ≤ 0.06), but there were no differences at d63. Yet, MPAR luminal and total ductal structure cross-sectional areas of IU-HT heifers were reduced at both d0 and d63 (P ≤ 0.01). The MFP adipocytes of IU-HT heifers were smaller at d0 (P ≤ 0.01), but differences were not detected at d63. The IU-HT heifers had diminished MPAR total, stromal, and epithelial cellular proliferation at both d0 and d63 (P < 0.01). Prenatal hyperthermia derails dairy calf early-life mammary development with potential carry-over consequences on future synthetic capacity.

Late gestation hyperthermia: epigenetic programming of daughter's mammary development and function

Exposure to stressors during early developmental windows, such as prenatally (i.e., in utero), can have life-long implications for an animal's health and productivity. The mammary gland starts developing in utero and, like other developing tissues and organs, may undergo fetal programming. Previous research has implicated factors, such as prenatal exposure to endocrine disruptors or alterations in maternal diet (e.g., maternal over or undernutrition), that can influence the developmental trajectory of the offspring mammary gland in postnatal life. However, the direct links between prenatal insults and future productive outcomes are less documented in livestock species. Research on in utero hyperthermia effects on early-life mammary development is scarce. This review will provide an overview of key developmental milestones taking place in the bovine mammary gland during the pre- and postnatal stages. We will showcase how intrauterine hyperthermia, experienced by the developing fetus during the last trimester of gestation, derails postnatal mammary gland development and impairs its synthetic capacity later in life. We will provide insights into the underlying histological, cellular, and molecular mechanisms taking place at key postnatal developmental life stages, including birth, weaning and the first lactation, that might explain permanent detriments in productivity long after the initial exposure to hyperthermia. Collectively, our studies indicate that prenatal hyperthermia jeopardizes the normal developmental trajectory of the mammary gland from fetal development to lactation. Further, in utero hyperthermia epigenetically programs the udder, and possibly other organs critical to lactation, yielding a less resilient and less productive cow for multiple lactations.

Late-Gestation in utero Heat Stress Limits Dairy Heifer Early-Life Growth and Organ Development

Dairy calves exposed to late-gestation heat stress weigh less, have impaired immunity, produce less milk across multiple lactations, and have reduced productive life. However, less is known about the relationship between in utero heat stress and organ morphology and development. Herein, we characterized the consequences of late-gestation in utero heat stress on body and organ growth trajectories during early-life development. Holstein heifers were either in utero heat-stressed (IU-HT, n = 36, dams exposed to THI &gt; 68) or cooled (IU-CL, n = 37, dams exposed to THI &gt; 68 with access to active cooling) during late gestation (54 ± 5 d prepartum). All heifers were reared identically from birth to weaning. Upon birth, calves were weighed and fed 3.78 L of colostrum followed by 0.87 kg DM/d milk replacer (MR) over two feedings and ad libitum starter concentrate daily. Weaning began at 49 d and ended at 56 d of age. Feed intake was recorded daily, and body weight (BW) and other growth measures were recorded at 0, 28, 56, and 63 d. Blood was collected at d 1 then weekly. Subsets of heifers were selected for euthanasia at birth and 7 d after complete weaning (n = 8 per group each) to harvest and weigh major organs. Reduced BW and stature measures persisted in IU-HT heifers from 0 to 63 d of age with a 7% lower average daily gain and reduced starter consumption relative to IU-CL heifers. IU-HT heifers had lower hematocrit percentages and reduced apparent efficiency of absorption of IgG relative to IU-CL heifers. Additionally, IU-HT heifers had reduced gross thymus, spleen, thyroid gland, and heart weight at birth and larger adrenal glands and kidneys but smaller ovaries relative to BW at 63 d. The mammary gland of IU-HT heifers was smaller relative to IU-CL heifers at birth and 63 d adjusted for BW, suggesting mechanisms leading to impaired milk yield in mature IU-HT cows are initiated early in development. In summary, in utero heat stress reduces whole-body size and limits development of key organs with potential repercussions on dairy calf metabolic adaptation, immune function, and future productivity.

Growth performance and mammary gland development of Holstein calves fed milk with 15% solids and treated with growth promoters

The objective of this study was to evaluate growth performance, and mammary gland parenchyma (PAR) development in Holstein female calves fed whole milk combined with milk replacer for 15% total solids (15TS) or only milk (11% TS; 11TS). The effects of zilpaterol hydrochloride (ZH) and an estrogenic implant (EI) on growth performance, serum metabolites, and PAR development in post-weaning calves (63-90 days of age) fed 15TS or 11TS were also investigated. In the first phase, 78 calves were randomly divided into two groups (n = 39/group), 15TS, or 11TS (5 to 8 L/day; starter ad libitum). After weaning, calves were randomly assigned to four treatments, no EI, EI, ZH, and ZH-EI. In phase 1, 15TS calves had greater (P < 0.01) average body weight at weaning than 11TS calves (82.0 ± 5.8 vs. 74.1 ± 5.7 kg). Pre-weaning average daily gain (718 ± 67 vs. 576 ± 64 g/day) was higher (P < 0.01) for 15TS calves than 11TS calves. Feed/gain (1.5 ± 0.3 vs. 1.8 ± 0.3 kg of DMI/kg of gain) was lower (P < 0.01) for 15TS calves than 11TS calves. PAR area (8.8 ± 3.3 vs. 5.0 ± 1.6 cm²) was greater (P < 0.01) for 15TS calves than 11TS calves. No differences (P > 0.10) in average daily gain (986 ± 79 g/day across treatments) were observed in calves receiving ZH, EI, or ZH-EI fed 15TS or 11TS diets previously. However, the PAR area of 15TS calves administered ZH-EI was 2.2 times greater (P < 0.01) than 11TS calves with no ZH and CI. This study demonstrated the benefits of supplying milk with 15% TS for enhancing overall animal growth and PAR development. Further use of ZH and EI during the post-weaning period did improve mammary PAR growth.

Effects of sodium salicylate and time postpartum on mammary tissue proliferation, gene transcript profile, and DNA methylation

Previous studies have demonstrated nonsteroidal antiinflammatory drug treatment in early lactation had a positive impact on whole-lactation milk production in older cows. The objective of this study was to evaluate proliferative, transcriptional, and epigenetic changes in the mammary gland that could explain increased production responses due to nonsteroidal antiinflammatory drug treatment. Sodium salicylate (SAL; 125 g/d) or water (CON) were administered via oral drench to multiparous Holstein cows (n = 8/treatment) once daily for 3 d beginning approximately 24 h after parturition, and mammary tissue was collected on d 1, 4, and 45 postpartum. Day 1 tissue was collected immediately preceding the initial drench, and d 4 tissue was collected 24 h following the final drench. Blood was collected twice weekly and analyzed for plasma glucose, insulin, β-hydroxybutyrate, free fatty acids, and prolactin. Cows were milked twice daily until d 7 of lactation, and thrice daily for the remainder of the study. Total RNA extracted from tissue was deep-sequenced and analyzed for differential gene expression using DESeq2. We detected no treatment effect on milk yield or plasma metabolites through 45 d of lactation; additionally, no change in mammary epithelial cell proliferation was detected when assessed by Ki67 labeling. Comparison of SAL versus CON revealed that only 16 of 18,286 genes were differentially expressed (false discovery rate <0.1) in mammary tissue collected on d 45, whereas no differentially expressed genes due to treatment were detected on d 1 or 4. Analysis of transcriptional differences over time showed downregulation of pathways related to immune cell recruitment and differentiation, and extensive overlap with pathways related to cholesterol synthesis and liver X receptor signaling. Global DNA methylation of mammary tissue was decreased for CON compared with SAL. Transcriptome analysis emphasized extensive involvement of immune-related signaling pathways in the switch from lactogenesis to galactopoiesis, and changes in methylation with SAL treatment merit future investigation into epigenetic effects on milk production.

Organization of mammary blood vessels as affected by mammary parenchymal region and estradiol administration in Holstein heifer calves

Mammary blood flow is central to mammary growth, development, and productivity, but the development of the vasculature network is poorly understood. The objective of this study was to determine how the vascular system adapts to mammary growth by inducing different levels of mammary growth and examining 2 regions of mammary parenchymal tissue. Holstein heifer calves (n = 12) received daily injections on the days immediately preceding euthanasia at 82 d of age. Treatments were control (CON), short-term estradiol (STE), and long-term estradiol (LTE). The CON calves received corn oil injections, the STE calves received 9 injections of corn oil followed by 3 injections of estradiol, and the LTE calves received 12 estradiol injections. Mammary tissues were collected from the center and edge parenchymal regions of all right rear mammary glands to quantify the tissue area of various tissue structures, the percentage of proliferating epithelial cells, and the number and form of blood vessels. Results showed that LTE calves had a greater tissue area occupied by epithelium than CON and STE calves, and the epithelial area in CON and STE calves was similar. Edge parenchyma had a greater percentage of proliferating epithelial cells than center parenchyma across all treatment groups. In the edge region, LTE calves had the greatest percentage of proliferating epithelial cells, coinciding with greater epithelial area. The number of blood vessels per unit of tissue area was greater in center than in edge parenchyma; the corresponding vessel surface area per unit of tissue area followed the same pattern. Mammary blood vessel measures were not markedly influenced by estradiol treatment. These results highlight the marked difference in the number and organization of blood vessels in different mammary parenchyma regions but indicate that the effects of estradiol on stimulating mammary epithelial proliferation does not directly translate to increasing numbers of blood vessels.

Review: Importance of colostrum supply and milk feeding intensity on gastrointestinal and systemic development in calves

Feeding management of the postnatal and preweaning calf has an important impact on calf growth and development during this critical period and affects the health and well-being of the calves. After birth, an immediate and sufficient colostrum supply is a prerequisite for successful calf rearing. Colostrum provides high amounts of nutrient as well as non-nutrient factors that promote the immune system and intestinal maturation of the calf. The maturation and function of the neonatal intestine enable the calf to digest and absorb the nutrients provided by colostrum and milk. Therefore, colostrum intake supports the start of anabolic processes in several tissues, stimulating postnatal body growth and organ development. After the colostrum feeding period, an intensive milk feeding protocol, that is, at least 20% of BW milk intake/day, is required to realise the calf potential for growth and organ development during the preweaning period. Insufficient milk intake delays postnatal growth and may have detrimental effects on organ development, for example, the intestine and the mammary gland. The somatotropic axis as the main postnatal endocrine regulatory system for body growth is stimulated by the intake of high amounts of colostrum and milk and indicates the promotion of anabolic metabolism in calves. The development of the forestomach is an important issue during the preweaning period in calves, and forestomach maturation is best achieved by solid feed intake. Unfortunately, intensive milk-feeding programmes compromise solid feed intake during the first weeks of life. In the more natural situation for beef calves, when milk and solid feed intake occurs at the same time, calves benefit from the high milk intake as evidenced by enhanced body growth and organ maturation without impaired forestomach development during weaning. To realise an intensive milk-feeding programme, it is recommended that the weaning process should not start too early and that solid feed intake should be at a high extent despite intensive milk feeding. A feeding concept based on intensive milk feeding prevents hunger and abnormal behaviour of the calves and fits the principles of animal welfare during preweaning calf rearing. Studies on milk performance in dairy cows indicate that feeding management during early calf rearing influences lifetime performance. Therefore, an intensive milk-feeding programme affects immediate as well as long-term performance, probably by programming metabolic pathways during the preweaning period.

Effect of whole-milk allowance on liveweight gain and growth of parenchyma and fat pads in the mammary glands of dairy heifers at weaning

The rates of development of 2 tissues in mammary glands, parenchyma (PAR) and the mammary fat pad (MFP), in response to nutrition in early life might have a major bearing on lifetime milk production. Historical studies reported that feeding greater amounts of dietary nutrients from postweaning to puberty increased growth rates of heifers and stimulated the growth of MFP at the expense of PAR, which might suggest compromised mammary development and future milk production. The current study sought to determine if a higher volume of whole milk (8 vs. 4 L/d) offered to calves would increase rates of growth and development of PAR in mammary glands at weaning (1 to 12 wk). To measure these tissues, we developed 2 simple methods to assess the size of PAR and MFP at the time of screening using ultrasound. We report that calves offered 8 L/d of whole milk had greater rates of growth until weaning (0.86 ± 0.06 vs. 0.81 ± 0.09 kg/d), compared with calves offered 4 L/d. Ultrasonography showed that despite the faster rates of growth in calves offered 8 L/d of milk/d, the ratio of PAR:MFP depth was 40% less at weaning in the front glands (34%) compared with calves offered 4 L of milk/d. Rear glands were less impaired. The ultrasound methods developed here might be useful to monitor the development of mammary glands in response to different nutritional regimens during the preweaning period.

Form of calf diet and the rumen. I: Impact on growth and development

Preweaning diet is known to affect rumen tissue appearance at the gross level. The objectives of this experiment were to investigate effects of different preweaning diets on the growth and development of the rumen epithelium and on putative rumen epithelial stem and progenitor cell measurements at the gene and cell levels. Neonatal Holstein bull calves (n = 11) were individually housed and randomly assigned to 1 of 2 diets. The diets were milk replacer only (MRO; n = 5) or milk replacer with starter (MRS; n = 6). Diets were isoenergetic (3.87 ± 0.06 Mcal of metabolizable energy per day) and isonitrogenous (0.17 ± 0.003 kg/d of apparent digestible protein). Milk replacer was 22% crude protein, 21.5% fat (dry matter basis). The textured calf starter was 21.5% crude protein (dry matter basis). Water was available ad libitum and feed and water intake were recorded daily. Putative stem and progenitor cells were labeled by administering a thymidine analog (5-bromo-2'-deoxyuridine, BrdU; 5 mg/kg of body weight in sterile saline) for 5 consecutive days and allowed a 25-d washout period. Calves were killed at 43 ± 1 d after a 6 h exposure to a defined concentration of volatile fatty acids. We obtained rumen tissue from the ventral sac and used it for immunohistochemical analyses of BrdU (putative stem and progenitor cells) and Ki67 (cell proliferation), gene expression analysis, and morphological measurements via hematoxylin and eosin staining. Epithelial stem and progenitor cell gene markers of interest, analyzed by real-time quantitative PCR, were β1-integrin, keratin-14, notch-1, tumor protein p63, and leucine-rich repeat-containing G protein-coupled receptor 5. Body growth did not differ by diet, but empty reticulorumens were heavier in MRS calves (MRS: 0.67 ± 0.04 kg; MRO: 0.39 ± 0.04 kg). The percentage of label-retaining BrdU basale cells was higher in MRO calves than in MRS calves (2.0 ± 0.3% vs. 0.3 ± 0.2%, respectively). We observed a higher percentage of basale cells undergoing proliferation in MRS calves than in MRO calves (18.4 ± 2.6% vs. 10.8 ± 2.8%, respectively). Rumen epithelial gene expression was not affected by diet, but the submucosa was thicker in MRO calves and the epithelium and corneum/keratin layers were thicker in MRS calves. Presumptive stem and progenitor cells in the rumen epithelium were identifiable by their ability to retain labeled DNA in the long term, changed proliferative status in response to diet, and likely contributed to observed treatment differences in rumen tissue thickness.

Form of calf diet and the rumen. II: Impact on volatile fatty acid absorption

Diet is known to affect rumen growth and development. Calves fed an all-liquid diet have smaller and less developed rumens and a decreased ability to absorb volatile fatty acids (VFA) compared to calves fed both liquid and dry feed. However, it is unknown how rumens respond when challenged with a defined concentration of VFA. The objective of this study was to assess the effects of 2 different feeding programs on VFA absorption in preweaned calves. Neonatal Holstein bull calves were individually housed and randomly assigned to 1 of 2 diets. The diets were milk replacer only (MRO; n = 5) or milk replacer with starter (MRS; n = 6). Diets were isoenergetic (3.87 ± 0.06 Mcal of metabolizable energy per day) and isonitrogenous (0.17 ± 0.003 kg/d of apparent digestible protein). Milk replacer was 22% crude protein, 21.5% fat (dry matter basis). The textured calf starter was 21.5% crude protein (dry matter basis). Feed and ad libitum water intakes were recorded daily. Calves were exposed to a defined concentration of VFA buffer (acetate 143 mM, propionate 100 mM, butyrate 40.5 mM) 6 h before euthanasia on d 43 ± 1. Rumen fluid samples were obtained every 15 to 30 min for 6 h to measure the rate of VFA absorption. Rumen tissues were obtained from the ventral sac region and processed for morphological and immunohistochemical analyses of the VFA transporters monocarboxylate transporter 1 (MCT1) and 4 (MCT4). Body growth did not differ between diets, but empty reticulorumens were heavier in MRS than MRO calves (0.67 vs. 0.39 ± 0.04 kg) and MRS calves had larger papillae areas (0.76 vs. 15 ± 0.08 mm²). We observed no differences between diets in terms of the abundance of MCT1 and MCT4 per unit area. These results indicate that the extrapolated increase in total abundance of MCT1 or MCT4 in MRS calves was not due to increased transporter density per unit area. Modeled VFA absorption metrics (flux, mmol/h, or 6 h absorbed VFA in mmol) were not different across diets. These results demonstrate that the form of calfhood diet, whether solely MR or MR and starter, does not alter VFA absorption capacity when the rumen is exposed to a defined concentration of VFA at 6 wk of age.

Effect of weaning age on growth, mammary gland development, and immune function in Holstein Friesian calves fed conserved alfalfa (FiberStart)

This study aimed to evaluate intake, body growth, and the development of the rumen, mammary gland, and immune system in Holstein Friesian calves reared for 100 d on the commercially available feed FiberStart (conserved alfalfa, Medicago sativa; Fiber Fresh Feeds Ltd., Reporoa, New Zealand) and fed calf milk replacer (CMR) for either 56 or 91 d. Eighty calves (40 bulls and 40 heifer calves) were reared indoors in groups (n = 5 of the same sex/pen). All calves were fed 4 L of CMR/d (175 g/L of CMR) in 2 feeds/d for the first 10 d and then 1 feed/d until d 49 or 84. The calves were gradually weaned by d 56 (earlier weaned; n = 8 pens) and d 91 (later weaned; n = 8 pens). All calves were fed FiberStart ad libitum as the only solid feed source from d 1 to 100 of the study. Irrespective of treatment, all calves had similar body weights at d 0 (40.9 ± 3.0 kg) and d 49 (74.2 ± 5.1 kg) of the study. Calf sex had no effect on intake, growth, blood, and immune system parameters. Earlier-weaned calves consumed 18% more solid feed dry matter but had 16% lower body weight gain (28.9 vs. 38.5 kg, respectively) from d 56 to 100 relative to later-weaned calves, resulting in different body weight at 100 d (104 vs. 121 ± 1.3 kg). Although earlier-weaned calves could compensate for the loss of CMR dry matter and crude protein intake from d 56 to 100 by increasing forage intake, they were unable to compensate for the loss of energy from the CMR by increasing solid feed consumption. Plasma β-hydroxybutyrate concentrations were 52% greater in earlier-weaned calves than in later-weaned calves at d 84, indicating greater metabolic activity of the rumen wall. The duration of CMR feeding had no influence on humoral or cell-mediated immune functions of the calves, as evidenced by a lack of effect on antivaccine antibody responses as well as on immune gene expression. Earlier- versus later-weaned heifer calves had 5% lower mammary gland mass, indicating that greater energy supply increased mammary mass. The results of this experiment demonstrate the ability to artificially rear dairy calves on a conserved alfalfa as the only solid feed. Furthermore, earlier weaning off CMR promotes solid feed intake and an associated increase in blood β-hydroxybutyrate, an indicator of rumen development, but increasing the duration of CMR feeding improves growth and mammary gland mass by d 100. The implications of these findings on lifetime growth, health, and milk production in dairy heifers warrant further investigation.

Preweaning nutrient supply alters mammary gland transcriptome expression relating to morphology, lipid accumulation, DNA synthesis, and RNA expression in Holstein heifer calves

The objective of this study was to analyze the mammary gland transcriptome to determine how preweaning nutrient supply alters the molecular mechanisms that regulate preweaning mammary development. Holstein heifers were fed via milk replacer (MR) either an elevated level of nutrient intake (ELE; on average, 5.9 ± 0.2 Mcal of ME in 8.4 L of MR/d, n = 6) or a restricted amount of nutrients (RES; 2.8 ± 0.2 Mcal of ME in 4 L of MR/d, n = 5) for 54 d after birth, at which point they were slaughtered and samples of mammary parenchyma tissue were obtained. Parenchymal mRNA was analyzed, and the fold change (FC) of 18,111 genes (ELE relative to RES) was uploaded to Ingenuity Pathway Analysis (IPA) software (Qiagen Bioinformatics, Redwood City, CA) for transcriptomic analysis. Using a threshold of P < 0.05, IPA identified that the FC of 1,931 of 18,811 differentially expressed genes (DEG) could be used for the analysis. A total of 18 molecular and cellular functions were relevant to DEG arising from the treatments; the 5 functions most associated with DEG were cell death and survival, cellular movement, cellular development, cellular growth and proliferation, and lipid metabolism. Based on the directional FC of DEG, the mammary gland of ELE heifers was predicted to have increased epithelial-mesenchymal transition (Z = 2.685) and accumulation of lipid (Z = 2.322), whereas the synthesis of DNA (Z = -2.137), transactivation of RNA (Z = -2.254), expression of RNA (Z = -2.405), transcription (Z = -2.482), and transactivation (Z = -2.611) were all predicted to be decreased. Additionally, IPA predicted the activation status of 13 upstream regulators with direct influence on DEG as affected by ELE feeding that were ligand-dependent nuclear receptors (n = 2), enzymes (n = 1), or transcription regulators (n = 10). Of these, 6 were activated (Z > 2) and 7 were inhibited (Z < -2). In summary, feeding ELE preweaning altered the mammary transcriptome of Holstein heifers, affecting cell functions involved in the morphological and physiological development of the mammary gland.

Higher plane of nutrition pre-weaning enhances Holstein calf mammary gland development through alterations in the parenchyma and fat pad transcriptome

Background

To reduce costs of rearing replacement heifers, researchers have focused on decreasing age at breeding and first calving. To increase returns upon initiation of lactation the focus has been on increasing mammary development prior to onset of first lactation. Enhanced plane of nutrition pre-weaning may benefit the entire replacement heifer operation by promoting mammary gland development and greater future production.

Methods

Twelve Holstein heifer calves (< 1 week old) were reared on 1 of 2 dietary treatments (n = 6/group) for 8 weeks: a control group fed a restricted milk replacer at 0.45 kg/d (R, 20% crude protein, 20% fat), or an accelerated group fed an enhanced milk replacer at 1.13 kg/d (EH, 28% crude protein, 25% fat). At weaning (8 weeks), calves were euthanized and sub-samples of mammary parenchyma (PAR) and mammary fat pad (MFP) were harvested upon removal from the body. Total RNA from both tissues was extracted and sequenced using the Illumina HiSeq2500 platform. The Dynamic Impact Approach (DIA) and Ingenuity Pathway Analysis (IPA) were used for pathway analysis and functions, gene networks, and cross-talk analyses of the two tissues.

Results

When comparing EH vs R 1561 genes (895 upregulated, 666 downregulated) and 970 genes (506 upregulated, 464 downregulated) were differentially expressed in PAR and MFP, respectively. DIA and IPA results highlight a greater proliferation and differentiation activity in both PAR and MFP, supported by an increased metabolic activity. When calves were fed EH, the PAR displayed transcriptional signs of greater overall organ development, with higher ductal growth and branching, together with a supportive blood vessel and nerve network. These activities were mediated by intracellular cascades, such as AKT, SHH, MAPK, and Wnt, probably activated by hormones, growth factors, and endogenous molecules. The analysis also revealed strong communication between MFP and PAR.

Conclusion

The transcriptomics and bioinformatics approach highlighted key mechanisms that mediate the mammary gland response to a higher plane of nutrition in the pre-weaning period.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-5303-8) contains supplementary material, which is available to authorized users.

TRIENNIAL LACTATION SYMPOSIUM/BOLFA: Plasticity of mammary development in the prepubertal bovine mammary gland

Although peripubertal mammary development represents only a small fraction of the total mass of mammary parenchyma present in the udder at the end of gestation and into lactation, there is increasing evidence that the tissue foundations created in early life can affect future mammary development and function. Studies on expression of estrogen and progesterone receptors seem to confirm the relevance of these steroids in prepubertal mammary development, but connections with other growth factors, hormones, and local tissue factors remain elusive. Enhanced preweaning feeding in the bovine appears to enhance the capacity of mammary tissue to response to mammogenic stimulation. This suggests the possibility that improved early nutrition might allow for creation of stem or progenitor cell populations to better support the massive ductal growth and lobulo-alveolar development during gestation. Increasing evidence that immune cells are involved in mammary development suggests there are unexpected and poorly understood connections between the immune system and mammary development. This is nearly unexplored in ruminants. Development of new tools to identify, isolate, and characterize cell populations within the developing bovine mammary gland offer the possibility of identifying and perhaps altering populations of mammary stem cells or selected progenitor cells to modulate mammary development and, possibly, mammary function.

Effects of rumen-undegradable protein on intake, performance, and mammary gland development in prepubertal and pubertal dairy heifers

The objective of this study was to evaluate the influence of different amounts of rumen-undegradable protein (RUP) on intake, N balance, performance, mammary gland development, carcass traits, and hormonal status of Holstein heifers at different physiological stages (PS). Sixteen prepubertal (PRE) heifers (initial BW = 106 ± 7.6 kg; age = 4.3 ± 0.46 mo) and 16 pubertal (PUB) heifers (initial BW = 224 ± 7.9 kg; age = 12.6 ± 0.45 mo) were used in an experiment over a period of 84 d. Four diets with increasing RUP contents (38, 44, 51, and 57% of dietary crude protein) and heifers at 2 PS (PRE or PUB) were used in a 4 × 2 factorial arrangement of treatments in a completely randomized design. Throughout the experiment, 2 digestibility trials were performed over 5 consecutive days (starting at d 36 and 78) involving feed and ort sampling and spot collections of feces and urine. At d 0 and 83, body ultrasound images were obtained for real-time carcass trait evaluation. The mammary gland was ultrasonically scanned at d 0 and every 3 wk during the experiment. Blood samples were taken at d 0 and 84 to determine serum concentrations of progesterone, estrogen, insulin-like growth factor I (IGF-I), and insulin. No interaction between PS and the level of RUP was found for any trait. Apparent digestibility of dry matter, organic matter, and neutral detergent fiber corrected for ash and protein was not affected by RUP level but was lower for PRE compared with PUB heifers. Sorting against neutral detergent fiber corrected for ash and protein (tendency only) and for crude protein was greater for PUB than PRE heifers. Pubertal heifers had greater average daily gain (905 vs. 505 g/d) and N retention (25.9 vs. 12.5 g/d) than PRE heifers. In addition, average daily gain and N retention were greatest at 51% RUP of dietary protein. Mammary ultrasonography indicated no effects of RUP amounts on mammary gland composition, whereas PRE heifers had greater pixel values than PUB, indicating higher contents of fat rather than protein in the mammary glands of PRE heifers. Serum progesterone and IGF-I concentration was affected only by PS, and PRE heifers had greater values of progesterone and IGF-I concentrations than PUB heifers. Serum insulin concentration was unaffected by PS but tended to be higher at 51% of RUP. In conclusion, an RUP level of 51% increases body weight, average daily gain, feed efficiency, and N retention in heifers regardless of the PS. In addition, PRE heifers have a lower sorting ability and reduced intake, total-tract digestibility, and N retention. They also have higher amounts of fat in their mammary glands, even at moderate growth rates.

Different milk diets have substantial effects on the jejunal mucosal immune system of pre-weaning calves, as demonstrated by whole transcriptome sequencing

There is increasing evidence that nutrition during early mammalian life has a strong influence on health and performance in later life. However, there are conflicting data concerning the appropriate milk diet. This discrepancy particularly applies to ruminants, a group of mammals that switch from monogastric status to rumination during weaning. Little is known regarding how the whole genome expression pattern in the juvenile ruminant gut is affected by alternative milk diets. Thus, we performed a next-generation-sequencing-based holistic whole transcriptome analysis of the jejunum in male pre-weaned German Holstein calves fed diets with restricted or unlimited access to milk during the first 8 weeks of life. Both groups were provided hay and concentrate ad libitum. The analysis of jejunal mucosa samples collected 80 days after birth and four weeks after the end of the feeding regimes revealed 275 differentially expressed loci. While the differentially expressed loci comprised 67 genes encoding proteins relevant to metabolism or metabolic adaptation, the most distinct difference between the two groups was the consistently lower activation of the immune system in calves that experienced restricted milk access compared to calves fed milk ad libitum. In conclusion, different early life milk diets had significant prolonged effects on the intestinal immune system.

Performance strategies affect mammary gland development in prepubertal heifers

In Brazil, the majority of dairy cattle are Holstein × Gyr (H×G). It is unknown whether excessive energy intake negatively affects their mammary development to the same extent as in purebred Holsteins. We hypothesized that mammary development of H×G heifers can be affected by dietary energy supply. We evaluated the effect of different average daily gains (ADG) achieved by feeding different amounts of a standard diet during the growing period on biometric measurements, development of mammary parenchyma (PAR) and mammary fat pad (MFP), and blood hormones. At the outset of this 84-d experiment, H×G heifers (n = 18) weighed 102.2 ± 3.4 kg and were 3 to 4 mo of age. Heifers were randomly assigned to 1 of 3 ADG programs using a completely randomized design. Treatments were high gain (HG; n = 6), where heifers were fed to gain 1 kg/d; low gain (LG; n = 6), where heifers were fed to gain 0.5 kg/d; and maintenance (MA; n = 6), where heifers were fed to gain a minimal amount of weight per day. Heifers were fed varying amounts of a single TMR to support desired BW gains. Over the 84 d, periodic biometric and blood hormone measurements were obtained. On d 84, all heifers were slaughtered and carcass and mammary samples were collected. At the end, HG heifers weighed the most (181 ± 7.5 kg), followed by LG (146 ± 7.5 kg) and MA (107 ± 7.5 kg) heifers. The ADG were near expected values and averaged 0.907, 0.500, and 0.105 ± 0.03 kg/d for HG, LG, and MA, respectively. In addition, body lengths, heart girths, and withers heights were affected by dietary treatment, with MA heifers generally being the smallest and HG heifers generally being the largest. Body condition scores differed by treatment and were highest in HG and lowest in MA heifers; in vivo subcutaneous fat thickness measurement and direct analysis of carcass composition supported this. The HG heifers had the heaviest MFP, followed by LG and then MA heifers. Amount of PAR was highest in LG heifers and was the same for HG and MA heifers. The percentage of udder mass occupied by PAR was lowest in HG heifers, differing from LG and MA heifers. Composition of MFP was not evaluated. Regarding PAR composition, no differences in ash or DM were found. On the other hand, CP concentration of PAR for HG heifers was lower than that for LG heifers, which was lower than that for MA heifers. Regarding the fat content, HG treatment was higher than LG and MA treatment, which did not differ from each other. In PAR, differences in relative abundance of genes related to both stimulation and inhibition of mammary growth were observed to depend on dietary treatment, sampling day, or both. The same can be said for most of the blood hormones that were measured in this experiment. In this experiment, high ADG achieved by feeding different amounts of a standard diet during the growing period negatively affected mammary development.

Ad libitum milk replacer feeding, but not butyrate supplementation, affects growth performance as well as metabolic and endocrine traits in Holstein calves

The enhanced growth performance of calves fed a higher plane of nutrition pre-weaning is well documented, and the effect of butyrate on the development of the gastrointestinal tract in calves has been evaluated. The aim of this study was to examine the synergistic effects of ad libitum milk replacer (MR) feeding and butyrate supplementation on growth performance and energy metabolism in calves. Sixty-four (32 male, 32 female) Holstein calves were examined from birth until wk 11 of life. Calves received MR either ad libitum (Adl) or restrictively (Res) with (AdlB+, ResB+) or without (AdlB-, ResB-) 0.24% butyrate supplementation. Colostrum and transition milk were fed in predefined amounts (Res or Adl) for the first 3 d postpartum. Ad libitum and restrictive MR feeding with or without butyrate was performed from d 4 until wk 8 of age. From wk 9 to 10, all calves were gradually weaned and were fed 2 L/d until the end of the trial. Concentrate (CON), hay, and water were freely available. Intakes of MR and CON were measured daily. Calves were weighed at birth and weekly thereafter. Blood was drawn on d 1 before the first colostrum intake; on d 2, 4, and 7; and weekly thereafter until the end of the study to measure plasma concentrations of metabolites and hormones. Liver samples were taken at d 50 and at the end of the study to determine gene expression related to glucose metabolism. Milk, MR, and total nutrient intake were greater, but CON intake was lower in Adl than in Res calves, resulting in a greater body weight, but partially lower gain to feed ratio in Adl than in Res. Plasma concentrations of glucose and insulin were higher during the ad libitum milk-feeding period, whereas plasma β-hydroxybutyrate was lower in Adl than in Res. Plasma concentrations of nonesterified fatty acids, lactate, total bilirubin, and cortisol were lower, but triglyceride and cholesterol concentrations were higher in Adl than in Res at specific time points. Feed intake, growth performance, and metabolic and endocrine changes were insignificantly affected by butyrate, and hepatic gene expression of enzymes related to endogenous glucose production was barely influenced by ad libitum MR feeding and butyrate supplementation. Intensive MR feeding indicated greater stimulation of growth and anabolic metabolism, but butyrate supplementation did not further improve postnatal growth or anabolic processes either in intensive or restrictive MR-fed calves.

Growth, intake, and health of Holstein heifer calves fed an enhanced preweaning diet with or without postweaning exogenous estrogen

Research has shown that changes in nutrition both before and after weaning can affect mammary development. Additionally, estrogen is known to be a potent mammogenic stimulant. Our objectives were to determine effects of altered preweaning feeding and exogenous estradiol postweaning on growth, intake, and health. Thirty-six Holstein heifer calves were reared on (1) a restricted milk replacer (MR) diet fed at 0.44kg powder dry matter (DM)/day [R; 20.9% crude protein (CP), 19.8% fat, DM basis], or (2) an enhanced MR fed at 1.08kg powder DM/d (EH; 28.9% CP, 26.2% fat, DM basis). The MR feeding was reduced 50% during wk 8 to prepare for weaning. Starter was offered after wk 4 but balanced between treatments. Body weight and frame were measured weekly with intakes and health monitored daily. At weaning, a subset of calves were slaughtered (n=6/diet). Enhanced-fed calves had greater carcass, thymus, liver, spleen, and mammary gland (parenchyma and mammary fat pad) weights. The EH calves also had greater average daily gain (ADG) starting during wk 1 (0.36 vs. -0.06kg/d) and lasting through wk 7 (1.00 vs. 0.41kg/d). Remaining calves received estrogen implants or placebo and were slaughtered at the end of wk 10, creating 4 treatments: (1) R, (2) R + estrogen (R-E2), (3) EH, and (4) EH + estrogen (EH-E2). Postweaning ADG was similar between R, EH, and EH-E2 calves, but greater in R-E2 calves than E calves. The EH-E2 calves had the heaviest mammary glands, and R-E2 calves had heavier mammary glands than R calves. The EH calves consumed more MR DM, CP, and fat preweaning. The R-fed calves consumed more starter DM preweaning. Fecal score was greater for EH calves (1.74 vs. 1.50) preweaning, but days medicated did not differ. Fecal scores were lower for R-E2 calves postweaning. Improved preweaning feeding of calves increased body weights and frame measures. Differences in body weights remained postweaning. Enhanced-fed calves showed greater ADG during the preweaning period but not postweaning. Exogenous estrogen may elicit diet-dependent growth responses. Analysis of collected samples will allow determination of cellular and molecular processes responsible for the marked differences in mammary development observed.

Tamoxifen impairs prepubertal mammary development and alters expression of estrogen receptor α (ESR1) and progesterone receptors (PGR)

Research has shown that prepubertal heifers experience allometric mammary growth that is influenced by the ovaries. Our purpose was to determine the role of estrogen in prepubertal mammary gland development. Sixteen Holstein calves were randomly assigned to 1 of 2 treatment groups: tamoxifen-injected (TAM) or control (CON). Calves were administered the antiestrogen tamoxifen (0.3 mg kg(1) d(1)) or placebo from 28 to 120 d of age. At 120 d, calves were euthanized and udders removed. Weight and DNA content of trimmed parenchymal tissue were halved (P ≤ 0.0001) in TAM compared with CON calves. Parenchymal samples from 3 zones of the left rear mammary gland (lower, middle, and outer regions) were processed for immunohistochemical staining for estrogen receptor α (ESR1) and progesterone receptor (PGR), Ki67-positive cells, and 5-bromo-2'-deoxyuridine label retaining cells (LRCs). Overall, neither the percentage nor location within the epithelial tissue layer of either ESR1- or PGR-positive cells was impacted by TAM treatment. However, image analysis indicated a 6.2-fold lower (P = 0.0001) level of ESR1 protein expression in TAM calves. Similarly, messenger RNA expression of ESR1 was also reduced (P = 0.0001) in TAM heifers. In contrast, expression of PGR protein was greater by 43% (P = 0.03) in TAM calves, but messenger RNA expression did not differ between treatments. Overall, TAM calves had a higher (P ≤ 0.03) percentage and density (cells per tissue area) of Ki67-positive cells. Irrespective of treatment, there were also more Ki67-labeled cells in the outer zones of the mammary gland (P ≤ 0.001). We were able to effectively use multispectral imaging to identify positive cells and quantify the expression of ESR1 and PGR protein. We also identified and counted the proportion of label retaining cells (LCR) (putative epithelial stem cells). We noted an overall 2.9-fold greater number of LRCs in TAM heifers and more LRCs in the outer sampling zones. This suggests that a cohort of LCR cells in TAM remained inactivated in comparison with CON heifers, which exhibited markedly increased growth of the mammary parenchymal tissue over the treatment period. These results suggest that the impacts of ovariectomy are partially explained by loss of ESR1 expression and/or estrogen receptor signaling in the prepubertal bovine mammary gland. The significance of mammary expression of PGR in control of prepubertal bovine mammary development remains unresolved.

Level of nutrient intake affects mammary gland gene expression profiles in preweaned Holstein heifers

Bovine mammary parenchyma (PAR) and fat pad (MFP) development are responsive to preweaning level of nutrient intake. We studied transcriptome alterations in PAR and MFP from Holstein heifer calves (n=6/treatment) fed different nutrient intakes from birth to ca. 65 d age. Conventional nutrient intake received 441 g of dry matter (DM)/d of a control milk replacer (MR) [CON; 20% crude protein (CP), 20% fat, DM basis]. Calves in the accelerated nutrition groups received 951 g/d of high-protein/low-fat MR (HPLF; 28% CP, 20% fat, DM basis), 951 g/d of high-protein/high-fat MR (HPHF; 28% CP, 28% fat, DM basis), or 1,431 g/d of HPHF (HPHF+) MR. Out of 13,000 genes evaluated, over 1,500 differentially expressed genes (DEG) were affected (false discovery rate <0.10) by level of nutrient intake in PAR or MFP. Feeding HPLF versus CON resulted in the most dramatic changes in gene expression, with 278 and 588 DEG having ≥1.5-fold change in PAR and MFP. In PAR, the most-altered molecular functions were associated with metabolism of the cell (molecular transport and lipid metabolism) with most of the genes downregulated in HPLF versus CON. In MFP, DEG also were primarily associated with metabolism but changes also occurred in genes linked to cell morphology, cell-to-cell signaling, and immune response. Compared with CON, feeding HPHF or HPHF+ did not result in substantial additional effects on DEG beyond those observed with HPLF. The pentose phosphate, mitochondrial dysfunction, and ubiquinone biosynthesis pathways were among the most enriched due to HPLF versus CON in PAR and were inhibited, whereas glycosphingolipid biosynthesis, arachidonic acid metabolism, and eicosanoid synthesis pathways were among the most enriched due to HPLF versus CON in MFP and were inhibited. These responses suggest that, in PAR, doubling nutrient intake from standard feeding rates inhibited energy metabolism and activity of oxidative pathways that partly serve to protect cells against oxidative stress. The MFP in those heifers appeared to decrease production of lipid-derived metabolites that may play roles in signaling pathways within the adipocyte. Overall, results indicated that prepubertal/preweaned mammary transcriptome is responsive to long-term enhanced nutrient supply to achieve greater growth rates before weaning. The biological significance of these results to future milk production remains to be elucidated.

Effect of staged ovariectomy on measures of mammary growth and development in prepubertal dairy heifers

Previous studies in prepubertal heifers suggest that the magnitude of reduction in mammary parenchymal growth in response to ovariectomy varies with the age at which surgery is performed. We hypothesized that ovarian secretions are essential for initiating mammary development but not required to maintain allometric mammary growth in prepubertal dairy heifers. The objectives of this study were to determine the effect of staged ovariectomy during the prepubertal period on mammary growth and tissue composition and the expression of selected genes. Prepubertal Holstein heifers at 2, 3 or 4 months of age were randomly assigned to one of two treatments, ovariectomized (OVX; n = 12) or sham operated (INT; n = 12). Mammary parenchyma (PAR) and fat pad (MFP) were harvested 30 days after surgery. Proximate composition of PAR and MFP (DNA, protein and lipid) as well as expression of the selected estrogen-responsive genes stanniocalcin1 (STC1), tissue factor pathway inhibitor precursor (TFPI) and proliferating cell nuclear antigen (PCNA) were determined in PAR and MFP by quantitative real-time PCR. The relative amount of epithelium and proportion of epithelia cell nuclei expressing the proliferation marker Ki67 were determined by histological and immunohistochemical analyses, respectively. MFP mass was not impacted by treatment but was decreased with age as was lipid content and concentration (P ⩽ 0.01). The mass of mammary PAR was reduced in OVX and increased with age (P ⩽ 0.01). Parenchymal tissue tended to have less total DNA, protein and lipid in OVX heifers. Parenchymal tissue concentrations of protein and DNA were increased with age and there was an age × treatment interaction. Treatment had no effect on either the Ki67 labeling index or percent epithelial area. The relative abundances of STC1, TFPI and PCNA mRNA in PAR were reduced in OVX. We did not find a significant impact of ovariectomy on mRNA expression when surgery was performed at 2 months compared with surgery at 3 or 4 months of age. However, having nearly undetectable PAR in two heifers ovariectomized at the earliest period (2 months of age) suggests that early ovariectomy is especially detrimental to subsequent parenchymal development.

Bovine mammary stem cells: cell biology meets production agriculture

Mammary stem cells (MaSC) provide for net growth, renewal and turnover of mammary epithelial cells, and are therefore potential targets for strategies to increase production efficiency. Appropriate regulation of MaSC can potentially benefit milk yield, persistency, dry period management and tissue repair. Accordingly, we and others have attempted to characterize and alter the function of bovine MaSC. In this review, we provide an overview of current knowledge of MaSC gained from studies using mouse and human model systems and present research on bovine MaSC within that context. Recent data indicate that MaSC retain labeled DNA for extended periods because of their selective segregation of template DNA strands during mitosis. Relying on this long-term retention of bromodeoxyuridine-labeled DNA, we identified putative bovine MaSC. These label-retaining epithelial cells (LREC) are in low abundance within mammary epithelium (<1%). They are predominantly estrogen receptor (ER)-negative and localized in a basal or suprabasal layer of the epithelium throughout the gland. Thus, the response of MaSC to estrogen, the major mitogen in mammary gland, is likely mediated by paracrine factors released by cells that are ER-positive. This is consistent with considerable evidence for cross-talk within and between epithelial cells and surrounding stromal cells. Excision of classes of cells by laser microdissection and subsequent microarray analysis will hopefully provide markers for MaSC and insights into their regulation. Preliminary analyses of gene expression in laser-microdissected LREC and non-LREC are consistent with the concept that LREC represent populations of stem cells and progenitor cells that differ with regard to their properties and location within the epithelial layer. We have attempted to modulate the MaSC number by infusing a solution of xanthosine through the teat canal and into the ductal network of the mammary glands of prepubertal heifers. This treatment increased the number of putative stem cells, as evidenced by an increase in the percentage of LREC and increased telomerase activity within the tissue. The exciting possibility that stem cell expansion can influence milk production is currently under investigation.

Functional and gene network analyses of transcriptional signatures characterizing pre-weaned bovine mammary parenchyma or fat pad uncovered novel inter-tissue signaling networks during development

BACKGROUND

The neonatal bovine mammary fat pad (MFP) surrounding the mammary parenchyma (PAR) is thought to exert proliferative effects on the PAR through secretion of local modulators of growth induced by systemic hormones. We used bioinformatics to characterize transcriptomics differences between PAR and MFP from approximately 65 d old Holstein heifers. Data were mined to uncover potential crosstalk through the analyses of signaling molecules preferentially expressed in one tissue relative to the other.

RESULTS

Over 9,000 differentially expressed genes (DEG; False discovery rate <or= 0.05) were found of which 1,478 had a >or=1.5-fold difference between PAR and MFP. Within the DEG highly-expressed in PAR vs. MFP (n = 736) we noted significant enrichment of functions related to cell cycle, structural organization, signaling, and DNA/RNA metabolism. Only actin cytoskeletal signaling was significant among canonical pathways. DEG more highly-expressed in MFP vs. PAR (n = 742) belong to lipid metabolism, signaling, cell movement, and immune-related functions. Canonical pathways associated with metabolism and signaling, particularly immune- and metabolism-related were significantly-enriched. Network analysis uncovered a central role of MYC, TP53, and CTNNB1 in controlling expression of DEG highly-expressed in PAR vs. MFP. Similar analysis suggested a central role for PPARG, KLF2, EGR2, and EPAS1 in regulating expression of more highly-expressed DEG in MFP vs. PAR. Gene network analyses revealed putative inter-tissue crosstalk between cytokines and growth factors preferentially expressed in one tissue (e.g., ANGPTL1, SPP1, IL1B in PAR vs. MFP; ADIPOQ, IL13, FGF2, LEP in MFP vs. PAR) with DEG preferentially expressed in the other tissue, particularly transcription factors or pathways (e.g., MYC, TP53, and actin cytoskeletal signaling in PAR vs. MFP; PPARG and LXR/RXR Signaling in MFP vs. PAR).

CONCLUSIONS

Functional analyses underscored a reciprocal influence in determining the biological features of MFP and PAR during neonatal development. This was exemplified by the potential effect that the signaling molecules (cytokines, growth factors) released preferentially (i.e., more highly-expressed) by PAR or MFP could have on molecular functions or signaling pathways enriched in the MFP or PAR. These bidirectional interactions might be required to coordinate mammary tissue development under normal circumstances or in response to nutrition.