ORIGINAL RESEARCH article

Front. Vet. Sci., 11 September 2020

Sec. Livestock Genomics

Volume 7 - 2020 | https://doi.org/10.3389/fvets.2020.00615

Synthetic Alkaloid Treatment Influences the Intestinal Epithelium and Mesenteric Adipose Transcriptome in Holstein Steers

  • 1. Ruminant Nutrition Laboratory, Department of Animal and Food Sciences, University of Kentucky, Lexington, KY, United States

  • 2. Department of Animal Science, University of Tennessee Institute of Agriculture, Knoxville, TN, United States

  • 3. Beltsville Agricultural Research Center, Agricultural Research Service, United States Department of Agriculture, Beltsville, MD, United States

  • 4. Forage-Animal Production Research Unit, Agricultural Research Service, United States Department of Agriculture, Lexington, KY, United States

Abstract

Holstein steers (n = 16) were used to determine if a synthetic alkaloid, bromocriptine, would alter the transcriptome of the small intestine and adjacent mesenteric adipose. On d 0, steers were assigned to one of two treatments: control (CON; saline only) or bromocriptine (BROMO; 0.1 mg/kg BW bromocriptine mesylate injected intramuscularly every 3 d for 30 d). Steers were slaughtered and midpoint sections of jejunal epithelium and associated mesenteric fat were collected for RNA isolation. Transcriptome analysis was completed via RNA-Seq to determine if BROMO differed compared with CON within intestinal epithelium or mesenteric adipose mRNA isolates. Differential expression thresholds were set at a significant P-value (P < 0.05) and a fold change ≥ 1.5. Only two genes were differentially expressed within the intestinal epithelium but there were 20 differentially expressed genes in the mesenteric adipose tissue (six up regulated and 14 down regulated). Functions related to cell movement, cell development, cell growth and proliferation, cell death, and overall cellular function and maintenance were the top five functional molecular categories influenced by BROMO treatment within the intestinal epithelium. The top molecular categories within mesenteric adipose were antigen presentation, protein synthesis, cell death, cell movement, and cell to cell signaling and interaction. In conclusion, BROMO treatment influenced the intestinal epithelium and mesenteric adipose transcriptome and identified genes and pathways influential to the effects associated with alkaloid exposure which are important to beef production.

Introduction

Tall fescue (Lolium arundinaceum) makes up a sizable portion of the grazing land in the United States and can reduce costs due to increased persistence and drought tolerance. The hardiness of tall fescue is attributed to the presence of an endophyte infection (Epichloë coenophiala). This symbiotic relationship between the plant and endophyte increases the sustainability of pastures and increases tall fescues usefulness as a grazing feed source (). However, tall fescue infected with endophyte decreases animal productivity due to ergot alkaloid ingestion from the endophyte infected plants (, ). Ingestion of ergot alkaloids from the plant itself (, ) and alkaloids produced by the fungi that infect tall fescue () have been implicated in fescue toxicosis. Animals experiencing fescue toxicosis exhibit a variety of symptoms such as increased respiration, impaired immune function, reduced feed intake, and decreased weight gain (). Cattle grazing tall fescue exhibit an unhealthy appearance even after access to fescue has ceased but can experience compensatory gains with increased diet quality (). Stuedemann and Hoveland () proposed that fescue toxicosis can alter lipid metabolism, adipose composition, and increase the occurrence of necrotic fat depots. However, instances of necrotic fat have only been reported in the abdominal fat depots () and mainly in older cattle which have grazed endophyte infected fescue over long periods of time.

Ergot alkaloids have been reported to inhibit prolactin secretion in cattle (, ). However, ergot alkaloid levels vary greatly in feedstuffs and make consistent exposure difficult in experimental procedures. Bromocriptine is a dopamine receptor agonist that effectively inhibits prolactin secretion () similar to ergot alkaloids. Lactating cows fed endophyte-infected fescue or treated with bromocriptine had over 90% of differentially expressed genes (>850 genes) influenced in a similar manner when compared with control cows (). Doses of bromocriptine have ranged from 15 to 80 mg/animal, and treatment frequencies have ranged from two times a day to three times per week () but provide a model to decrease prolactin and initiate fescue toxicosis symptoms in cattle with much greater control by not relying on animal intake. Thus, the hypothesis of this study was that ergot alkaloid exposure via bromocriptine treatment would impact the transcriptome of the intestinal epithelium and mesenteric adipose tissue relevant to beef production of growing Holstein steers without impacting feed intake.

Materials and Methods

All animal procedures were conducted with approval from the Institutional Animal Care and Use Committee at the University of Kentucky (Protocol # 01065A2006).

Animals

Holstein steers [n = 16; initial body weight (BW) 313 ± 16 kg] were used to determine if the synthetic alkaloid, bromocriptine, at levels known to induce fescue toxicosis symptoms () would alter the transcriptome of the jejunal epithelium and the adjacent mesenteric adipose. Steers were housed in individual pens (3 × 3 m) within a temperature-controlled research facility and allowed 3 d to adjust to being housed indoors prior to the initiation of the study (d 0). Steers were maintained on a total mixed ration of 61% corn silage, 24% cracked corn, 5% dry distiller's grains with solubles, and 10% soybean meal supplement with 150 mg/head/d rumensin. The overall diet was 49% dry matter, had a net energy for maintenance of 0.367 Mcal/kg, contained 13% crude protein, and was fed to provide 1.5 × net energy for maintenance daily. Maintenance was determined by multiplying empty BW raised to the 0.75 power by 0.77 and increased by 20% to address the increased maintenance energy requirement for dairy animals ().

On d 0 steers were randomly assigned to one of two treatments: control (CON) or bromocriptine (BROMO). Control animals were intramuscularly injected in the injection triangle of the neck with saline and ethanol (60 and 40%, respectively) at a similar volume to BROMO (~5 mL) every 3 d to coincide with bromocriptine injections. Animals who were assigned to BROMO received bromocriptine mesylate (Santa Cruz Biotechnology; Dallas, TX) at 0.1 mg/kg BW every 3 d beginning on d 0 and ending on d 27 resulting in 10 intramuscular injections. Bromocriptine was reconstituted as working stock in 95% ethanol and diluted with saline so no more than 40% ethanol was given in a single injection. The concentration of BROMO was determined in a previous study where BROMO give at this dosage elicited the same gene expression response in 90% of differentially expressed genes in mammary tissue (n = 866) identified comparing the BROMO treatment response with the response observed in cows fed endophyte-infected fescue seed (). Body weights were taken weekly to ensure feed was provided at 1.5 × maintenance and injection volumes calculated properly. Feed and water intake was individually provided and recorded daily. Orts were weighed daily prior to feeding to calculate daily intake for the previous 24 h. Steers were kept on automatic waterers that were equipped with a water meter to capture water intake. Daily water intake was recorded prior to feeding. Circulating levels of prolactin were used to confirm treatment effect on animal did in fact induce similar effects to fescue toxicosis. Blood samples were taken prior to feeding via jugular venipuncture with 10 mL vacutainer tubes (Becton Dickinson Co.; Franklin Lakes, NJ) upon entry into the housing facility (d −3) and d 0, 7, 14, 21, and 28 to ensure BROMO treated steers had decreased systemic prolactin concentrations. Serum prolactin concentrations were quantified via radioimmunoassay according to procedures previously described (). The intra and inter assay CVs were 5.35 and 4.37%, respectively, with a sensitivity of 0.05 ng/mL. Prolactin concentrations were analyzed within the mixed procedure of SAS (SAS Institute, Cary, NC) with repeated measures and a compound symmetry covariance structure. The model included treatment, day and the interaction of treatment×day and block as a random covariate. The repeated statement had day as the variable and animal within block as the subject. Least square means were used to separate main effect means with significance set as P < 0.05.

On d 30, at 07:00 h the steers were transported 26.7 km from the animal research center beef unit to the meats laboratory at the University of Kentucky. After exsanguination the gastrointestinal tract was removed from the body cavity as quickly as possible (<10 min). The gastrointestinal tract was separated into the forestomach complex (rumen, reticulum, omasum, and abomasum), small intestine, and large intestine. The small intestine was looped over a peg board and a 10 cm section was removed 7 m from the cranial end. Mesenteric adipose adjacent to the 10 cm intestinal section was also removed for further processing. Following separation from the entire tract, 50 mg of adipose tissue submerged in RNAlater (Thermo Fisher Scientific; Waltham, MA) and flash frozen in liquid nitrogen. The small intestinal lumen was exposed, the epithelium was scraped away from the other intestinal tissue layers, and 50 mg of epithelium was submerged in RNAlater (Thermo Fisher Scientific) and flash frozen in liquid nitrogen.

RNA Isolation and Library Construction

Both jejunum and mesenteric adipose samples were isolated using RNeasy Plus Mini Kit (Qiagen, Valenica, CA). Quality check was performed using Tapestation RNA HS Assay (Agilent Technologies, Santa Clara, CA) and quantified by Qubit RNA HS assay (ThermoFisher Scientific, Waltham, MA). The RNA integrity number for samples was lower and more variable than expected (2.6 ± 0.6 for intestine and 5.0 ± 1.3 for adipose); thus, ribosomal RNA depletion was performed with Ribo-zero Magnetic Gold Kit (Illumina Inc., San Diego, CA) to ensure RNA isolated was of the highest possible quality. Samples were randomly primed and fragmented based on manufacturer's recommendation (NEBNext® Ultra™ RNA Library Prep Kit for Illumina® Illumina Inc., San Diego, CA). The first strand was synthesized with the Protoscript II Reverse Transcriptase for 40 min at 42°C. Synthesis of second strand cDNA was then carried out and products were amplified by PCR to create the cDNA library for each sample. All remaining steps for library construction were done according to the NEBNext® Ultra™ RNA Library Prep Kit for Illumina® (Illumina Inc.). Illumina 8-nt dual-indices were used. Samples were pooled by animal and sequencing was done on a HiSeq with a read length configuration of 150 paired end base pairs with up to 50 million reads per sample.

Bioinformatic Analysis

Intestinal and adipose library preparation and sequencing, was performed using FASTQC (version v0.11.3). A splice aware aligner, STAR (version 2.5.2b), was used to perform RNA-Seq alignment. The UMD3.1 and UMD3.1.90 from Ensembl were used as genome reference and annotation reference, respectively, during the alignment. Then dupRader and Picard CollectRnaSeqMetrics (version 2.10.5) were used to evaluate duplicates level and overall alignment performance. HT-Seq (version 0.0.6) was used to calculate the per gene expression count. DE-Seq2 was used to estimate differentially expressed genes, all genes are listed in Supplementary Table 1. This list of genes was uploaded into Ingenuity Pathway Analysis (IPA, Qiagen) for further analyses (Supplementary Table 2).

Ingenuity Pathway Analysis

Differentially expressed genes from RNA-Seq were analyzed via IPA according to previously described methods (). The RNA-Seq data was used by IPA to identify new gene targets, biological mechanisms, pathways, and functions, and candidate biomarkers. Datasets uploaded into the IPA software consisted of all genes identified as significant with a P ≤ 0.10 (Supplementary Table 2). From these datasets the most relevant genes were identified as differentially expressed genes (DEG). To qualify as a DEG, genes had to have a significant P-value and a positive or negative fold change of >1.5. Canonical pathway analysis identified the most significantly influenced canonical pathways from the IPA library within each tissue. Genes related to each canonical pathway were considered for further analysis. Two measurements are evaluated for canonical pathway analysis (1) the ratio of the number of genes from the data set that map to the pathway divided by the total number of genes that map to the canonical pathway and (2) the Fisher's exact test was used to calculate a P-value for the genes in the dataset and the identified canonical pathway. Networks of DEG were then algorithmically generated based on their connectivity. These genes were overlaid onto a global molecular network from the IPA library. Functional analysis identified biological functions associated with DEG that were most significantly influenced by treatment within a tissue. A right-tailed Fisher's exact test was used to calculate a P-value for each biological function and/or disease.

Results

Circulating prolactin concentrations were similar (P > 0.11) between BROMO and Con steers prior to treatment initiation, d −3 and 0 (Figure 1). However, as expected, plasma concentrations were dramatically reduced (P = 0.003) in BROMO steers compared with CON steers at all-time points following bromocriptine administration, d 7, 14, 21, and 28.

Figure 1

Threshold for determination of DEG was set at a significant P-value (P < 0.05) and a fold change of >1.5. Based on selection criteria, RNA-Seq analysis revealed 15 DEG that were down regulated in BROMO compared with CON; one DEG within the intestinal epithelium, fatty acid binding protein 4, and 16 DEG in the mesenteric adipose tissue (Table 1). Analysis identified 7 DEG were identified to be up regulated in BROMO compared with CON; 6 in the mesenteric and UDP-GlcNAc:betaGal Beta-1,3-N-Acetylglucosaminyltransferase in the intestinal epithelium (Table 2).

Table 1

SymbolGene nameFold changeTissue typeP-valueLocationType
FABP4Fatty Acid Binding Protein 41.842Intestinal0.020CytoplasmTransporter
CYP4F2Cytochrome P450 Family 4 Subfamily F Member 25.448Adipose0.048CytoplasmEnzyme
C8AComplement C8 Alpha Chain5.14Adipose0.019ExtracellularOther
ASIC5Acid Sensing Ion Channel Subunit Family Member 55.021Adipose0.041PMIon Channel
GLT1D1Glycosyltransferase 1 Domain Containing 14.965Adipose0.0278ExtracellularEnzyme
SERPINB2Serpin Family B Member 24.924Adipose0.053ExtracellularOther
AICDAActivation Induced Cytidine Deaminase4.902Adipose0.040CytoplasmEnzyme
CDH16Cadherin 164.83Adipose0.040PMEnzyme
CXCL2C-X-C Motif Chemokine Ligand 24.436Adipose0.039ExtracellularCytokine
CHGBChromogranin B4.356Adipose0.023ExtracellularOther
SLC26A5Solute Carrier Family 26 Member 54.037Adipose0.029PMTransporter
HEPACAM2HEPACAM Family Member 24.014Adipose0.050CytoplasmOther
NR1H4Nuclear Receptor Subfamily 1 Group H Member 43.862Adipose0.019NucleusLDNR
SLC9A2Solute Carrier Family 9 Member A23.628Adipose0.027PMTransporter
LEAP2Liver Enriched Antimicrobial Peptide 23.473Adipose0.051ExtracellularOther
DTHD1Death Domain Containing 12.326Adipose0.054OtherOther
CD8ACD8a Molecule1.535Adipose0.042PMOther

Differentially expressed genes* that were down regulated after exposure to Bromocriptine in the intestinal epithelium and mesenteric adipose tissue of Holstein steers.

*

Differentially expressed genes listed above were limited to those with a P-value <0.05 and >1.5 fold change.

PM, Plasma Membrane.

LDNR, ligand-dependent nuclear receptor.

Table 2

SymbolGene nameFold changeTissue typeP-valueLocationType
B3GNT6UDP-GlcNAc:betaGal Beta-1,3-N-Acetylglucosaminyltransferase 62.853Intestinal0.006CytoplasmEnzyme
MOXD1Monooxygenase DBH Like 11.772Adipose0.040CytoplasmEnzyme
SORCS2Sortilin Related VPS10 Domain Containing Receptor 21.833Adipose0.047PMTransporter
SLC11A1Solute Carrier Family 11 Member 11.988Adipose0.043PMTransporter
LY6G5BLymphocyte Antigen 6 Family Member G5B2.042Adipose0.050OtherOther
DNERDelta/Notch Like EGF Repeat Containing3.558Adipose0.027PMTransmembrane Receptor
PF4Platelet Factor 45.393Adipose0.007Extracellular SpaceCytokine

Differentially expressed genes* that were up regulated after exposure to Bromocriptine in the intestinal epithelium and mesenteric adipose tissue of Holstein steers.

*

Differentially expressed genes listed above were limited to those with a P-value <0.05 and >1.0 fold change.

PM, Plasma Membrane.

The canonical pathways influenced by BROMO treatment are listed in Table 3. These pathways were determined to be influenced from the all significant (P < 0.05) genes with and without a 1.5 fold change within the datasets (BROMO and CON) via IPA library. In the intestinal epithelium and mesenteric adipose, LXR/RXR activation, granulocyte adhesion, and diapedesis were depressed (P < 0.001) in BROMO steers. Atherosclerosis and high mobility group box 1 (HMGB1) signaling and the inflammasome pathway were also depressed (P < 0.001) in the intestinal epithelium of BROMO steers compare with CON steers. In the mesenteric adipose, LPS/IL-1 mediated inhibition of RXR function, FXR/RXR activation, antigen presentation, IL-6 signaling, dendritic cell maturation, and serotonin degradation were all depressed (P < 0.001) by BROMO treatment (Table 3).

Table 3

PathwayP-valueRatio*Genes
Intestinal epithelium
LXR/RXR Activation<0.00010.033↓CCL2, ↓IL18, ↓LPL, ↓MMP9
Atherosclerosis Signaling<0.00010.031↓CCL2, ↓IL18, ↓LPL, ↓MMP9
Inflammasome Pathway0.00010.100↓IL18 and ↑NLRP1
HMGB1 Signaling0.00010.023↓CCL2, ↓IL18, ↑PLAT
Granulocyte Adhesion and Diapedesis0.00030.017↓CCL2, ↓IL18, ↓MMP9
Mesenteric adipose
LPS/IL-1 Mediated Inhibition of RXR Function<0.00010.054↓ABCB11; ↓ALAS1; ↓ALDH1L1, 8A1; ↓FABP7; ↓IL1A; ↓NR1H4; ↓PAPSS2; ↓PLTP; ↓SULT2A1;
↑IL1RL1; ↑SULT1A1
FXR/RXR Activation<0.00010.071↓ABCB11; ↓APOA4, B; ↓IL1A; ↓NR1H4; ↓PLTP; ↓SDC1; ↓SULT2A1; ↑SERPINF1
Antigen Presentation Pathway<0.00010.132↓B2M; ↓PSMB8, 9; ↓TAP1, BP
Granulocyte Adhesion and Diapedesis<0.00010.050↓CCL8; ↓CXCL2; ↓IL1A; ↓MMP1; ↓SDC1; ↓XCL2; ↑CCL11; ↑IL1RL1; ↑PF4
LXR/RXR Activation0.00010.058↓APOA4, B; ↓IL1A; ↓NR1H4; ↓PLTP; ↑IL1RL1; ↑SERPINF1
Serotonin Degradation0.00510.052↓ADH4; ↓ALDH1L1; ↓SULT2A1; ↑SULT1A1
Dendritic Cell Maturation0.03100.026↓B2M; ↓FCGR3A, B; ↓IL1A; ↑COL1A1; ↑NFKBIA
IL-6 Signaling0.02840.031↓IL1A; ↑COL1A1; ↑IL1RL1; ↑NFKBIA

Canonical pathways influenced by treatment with bromocriptine in the intestinal epithelium and mesenteric adipose in Holstein steers.

*

Ratio: The number of genes differentially expressed within the pathway divided by the total number of genes within the pathway.

Genes that exhibited directional influence by treatment within a pathway.

HMGB1: High Mobility Group Box 1. The ↑ arrow is for genes that were up-regulated and the ↓ arrow is for a gene that was down-regulated.

The top 5 upstream regulators were different between the intestinal epithelium and mesenteric adipose (Tables 4, 5). However, upstream regulation by IL 1A, IL 2, IL 4, IL 6, IL 18, IGF 1, growth hormone (GH), extracellular signal–regulated kinases 1/2 (ERK1/2), and toll-like receptor 4 (TLR 4) were found in the top 30 upstream regulators for both tissue types. In the intestinal epithelium, regulator of G protein signaling (RGS2), MAPK 9, prostaglandin-endoperoxide synthase 2 (PTGS2), proteinase 3 (PRTN3), and histidine rich glycoprotein (HRG) were found to be the top 5 upstream regulators (Table 4). In the mesenteric adipose, TNF, interferon α receptor (IFNαR), IL 1B, NF-κβ, and IFNγ were found to be the top five upstream regulators (Table 5). Also, prolactin, which is depressed by BROMO treatment (Figure 1), was identified as an upstream regulation of the transcriptome in the mesenteric adipose of BROMO steers (Table 5).

Table 4

GenesUpstream Regulators*
Top 5IL1AIL2IL4IL6IL18ERK 1/2IGF 1GHIFNαRTLR4IFNγ
RGS2MAPK 9PTGS 2PRTN 3HRG
CCL2
FABP4
IL 18
LPL
MMP9
PLAT

Upstream regulators in the intestinal epithelium and how they influence gene expression.

*

Abbreviations for upstream regulators: RGS2, regulator of G protein signaling; PTGS2, Prostaglandin-endoperoxide synthase 2; PRTN3, Proteinase 3; HRG, Histidine Rich Glycoprotein; ERK1/2, extracellular signal–regulated kinases 1/2; GH, Growth Hormone; IFNαR, Interferon α Receptor; TLR4, Toll-like receptor 4.

Abbreviations for genes: CCL2, C-C motif ligand 2; FABP4, fatty acid binding protein 4; LPL, lipoprotein lipase; MMP9, Matrix metallopeptidase 9; PLAT, Plasminogen Activator. The ↑ arrow is for genes that were up-regulated and the ↓ arrow is for a gene that was down-regulated.

Table 5

GenesUpstream Regulators*
Top 5IL1AIL2IL4IL6IL18IGF-1GHTLR 4PRLERK 1/2
TNFIFNARIL-1BNF-κβIFNγ
B2M
CCL11
CCND1
COL1α1
CXCL2
GADD45
GBP2
IDO1
IL1A
IRF1
JUNB
MMP1
MYC
NFκβIα
PSMB8
PSMB9
RSAD2
TAP1

Upstream regulators in mesenteric adipose tissue and how they influence gene expression.

*

Abbreviations for upstream regulators: IFNαR, interferon α receptor; GH, growth hormone; TLR 4, toll-like receptor 4; PRL, prolactin; ERK1/2, extracellular signal–regulated kinases 1/2.

Abbreviations for genes: B2M, beta-2-microglobulin; CCL11, C-C motif ligand 11; CCND1, cyclin D1; COL1α1, collagen type 1 α chain 1; C-X-C motif, chemokine; CXCL2, ligand 2; GADD45, growth arrest and DNA damage 45; GBP2, guanine binding protein 2; IDO1, indoleamine 2,3-dioxygenase 1; IRF1, interferon regulatory factor 1; MMP1, Matrix metallopeptidase 1; NFκβIα, NF-κβ inhibitor α; PSMB8, proteasome subunit beta 8; PSMB9, proteasome subunit beta 9; RSAD2, radical S-adenosyl methionine domain containing 2; TAP1, transporter 1. The ↑ arrow is for genes that were up-regulated and the ↓ arrow is for a gene that was down-regulated.

Functional analysis of all significant (P < 0.05) genes, with and without 1.5 fold change within the datasets (BROMO and CON) was done via IPA library. Functions related to cell movement, cell development, cell growth and proliferation, cell death, and overall cellular function and maintenance were the top 5 functional molecular categories (P < 0.001) affected by BROMO treatment within the intestinal epithelium (Table 6). Whereas, hematological system development and function, immune cell trafficking, tissue morphology, skeletal and muscular system development and overall tissue development were the top physiological function categories within the epithelium of the intestine depressed (P < 0.001) in BROMO steers (Table 6). In the intestinal epithelium all individual functions were depressed (P < 0.001) in BROMO steers compared with CON steers (Table 7). Most of these functions were related to immune cell activation and immune response such as: cell movement of neutrophils, infiltration, activation, and accumulation of phagocytes, activation and accumulation of myeloid cells, quantity of macrophages, and overall inflammation. Other functions were related to metabolism, quantity of carbohydrates and proteins, synthesis of lipids, organismal survival, development of epithelial tissue, and vasculogenesis (Table 7).

Table 6

Top 5 molecular and cellular functionsTop 5 physiological system development and functions
Cellular MovementHematological System Development and Function
Cellular DevelopmentImmune Cell Trafficking
Cellular Growth and ProliferationTissue Morphology
Cell Death and SurvivalSkeletal and Muscular System Development and Function
Cellular Function and MaintenanceTissue Development

IPA* bio-functions enriched by significantly expressed genes in the intestinal epithelium of bromocriptine treated Holstein steers.

*

IPA, Ingenuity Pathway Analysis.

Table 7

IPA Bio-functionFunctional effectP-valuez-Score of activationGenes, n
Carbohydrate MetabolismQuantity of carbohydrate0.0017−1.1764
Cardiovascular System DevelopmentVasculogenesis0.0009−1.0735
Cell-To-Cell SignalingActivation of cells0.0015−1.5565
Activation of myeloid cells<0.0001−1.2824
Activation of phagocytes0.0001−1.2824
Cellular MovementInfiltration by myeloid cells<0.0001−1.4455
Chemotaxis0.0018−1.1334
Cell movement of neutrophils<0.0001−1.6955
Cellular infiltration by phagocytes<0.0001−1.0675
Migration of smooth muscle cells<0.0001−1.2824
Immune Cell Trafficking and Inflammatory ResponseAccumulation of phagocytes<0.0001−1.0674
Quantity of macrophages<0.0001−1.6124
Accumulation of myeloid cells<0.0001−1.1924
Inflammatory response0.0008−1.675
Inflammation of absolute anatomical region<0.0001−1.5337
Inflammation of organ<0.0001−1.41410
Synthesis of lipid0.0043−1.9634
Organismal SurvivalSurvival of organism0.00091.1385
Protein SynthesisMetabolism of protein0.0026−1.9275
Quantity of protein in blood0.0018−1.064
Tissue DevelopmentDevelopment of epithelial tissue0.0001−1.1035

Individual functions* activated or depressed within the intestinal epithelium of Holstein steers treated with bromocriptine compared to saline treated steers.

*

Functions were selected upon the criteria of having a P value < 0.05 and an absolute value for the z score of activation of > 1.

PA, Ingenuity Pathway Analysis.

A negative z score of activation indicates the function was depressed and a positive z score of activation indicates the functions was activated.

Ingenuity pathway analysis of mesenteric adipose determined that the top molecular categories were antigen presentation, protein synthesis, cell death, cell movement, and cell to cell signaling and interaction (Table 8). Whereas, the top physiological function categories were organismal survival, hematological system development and function, tissue morphology, organismal development, and immune cell trafficking (Table 8). The complete list of individual functions influenced by BROMO treatment in the mesenteric adipose is provided in Supplementary Table 1. Physiologically relevant functions to the hypothesis of the study (P < 0.001) related to metabolism, tissue and cell health and growth, cell movement, and some immune and inflammatory responses were down regulated (P < 0.001; Table 9). The few immune response functions that showed activation (P < 0.001) with BROMO treatment were recruitment of macrophages, activation of myeloid cells and leukocytes, and proliferation of mononuclear leukocytes and lymphatic system cells. There were multiple gastrointestinal disease functions, focused largely around inflammation, activated due to exposure to BROMO including enteritis and colitis (Table 9). Functions that were related to metabolism consisted of mainly lipid metabolism functions, transport of lipid, cholesterol, and steroids and the quantity of HDL in the blood. Cellular movement and migration of both epithelial and smooth muscles cells was depressed in BROMO steers compared with CON steers (Table 9).

Table 8

Top 5 molecular and cellular functionsTop 5 physiological system development and functions
Antigen PresentationOrganismal Survival
Protein SynthesisHematological System Development and Function
Cell Death and SurvivalTissue Morphology
Cellular MovementOrganismal Development
Cell to Cell Signaling and InteractionImmune Cell Trafficking

IPA* bio-functions enriched by significantly expressed genes in mesenteric adipose of bromocriptine treated Holstein steers.

*

IPA, Ingenuity Pathway Analysis.

Table 9

IPAFunctional effectP-valuez-Score of activationGenes, n
Bio-function
Antigen PresentationActivation of leukocytes<0.00011.02819
Proliferation of mononuclear leukocytes0.00412.01218
Proliferation of lymphatic system cells0.00401.819
Gastrointestinal DiseaseGastroenteritis0.00141.58813
Enteritis0.00251.35212
Colitis0.00311.74111
Inflammation of gastrointestinal tract0.00161.58813
Abnormality of large intestine0.00241.96912
Hematological Disease, Development, and FunctionBlood protein disorder0.00081.95416
Lymphoproliferative disorder<0.00011.47547
Quantity of hematopoietic progenitor cells0.0012−1.22614
Quantity of phagocytes<0.0001−1.3516
Quantity of macrophages<0.0001−1.11711
Hematological System Development and FunctionQuantity of T lymphocytes0.0035−1.47615
Quantity of lymphocytes0.0011−1.08420
Quantity of blood cells<0.0001−2.14834
Quantity of granulocytes0.0004−1.55112
Quantity of leukocytes<0.0001−1.50929
Quantity of mononuclear leukocytes0.0007−1.39921
Quantity of myeloid cells<0.0001−1.14418
Humoral Immune ResponseQuantity of immunoglobulin0.0022−1.26110
Inflammatory Disease and ResponseInflammation of body cavity<0.00011.25826
Inflammation of absolute anatomical region<0.00011.9733
Organismal DevelopmentGrowth of organism0.00061.1622
Organismal Injury and AbnormalitiesCytosis0.0011−1.36212
Tissue MorphologyQuantity of cells<0.0001−1.51949

Individual functions* activated or depressed within mesenteric adipose of Holstein steers treated with bromocriptine compared to saline treated Holstein steers.

*

Functions were selected upon the criteria of having a P value < 0.05, more than 10 genes identified within, and an absolute value for the z score of activation of > 1.

IPA, Ingenuity Pathway Analysis.

A negative z score of activation indicates the function was depressed and a positive z score of activation indicates the functions was activated.

Discussion

The widespread usage of tall fescue for grazing animals makes the decrease in productivity due to fescue toxicosis economically important to beef production. Tall fescue (Lolium arundinaceum) has a symbiotic relationship to a fungal endophyte (Epichloë coenophiala) which infects mainly the seed head of the plant but also confers growth and heartiness advantages which increase the usefulness of tall fescue as a grazing feed source (). The endophytic infection creates a source of ergot alkaloids of which ergovaline is found in the greatest concentrations (). Consumption of these associated ergot alkaloids in sufficient amounts to cause fescue toxicosis in cattle results in a dramatic reduction in circulating prolactin concentration; a diagnostic tool used to confirm exposure to detrimental ergot alkaloids (, , ). Bromocriptine, a synthetic ergopeptine and dopamine receptor agonist, shares a high degree of structural homology with ergovaline and has also been shown to reduce circulating prolactin concentration in cattle (, ). More recently we have shown that 90% of DEG (n = 866) in mammary tissue respond similarly for both fed fescue-derived alkaloids and BROMO treated cows (); indicating that BROMO provides an adequate model for studying physiological changes in cattle grazing endophyte-infected tall fescue. In the current study we extend our finding concerning the effects of BROMO on the transcriptome to the intestinal epithelium and mesenteric adipose in growing ruminant animals.

Intestinal Epithelium

Based on the current dataset and previous work (, ) in dairy cows BROMO impacted the intestinal transcriptome less than that observed in the transcriptome of mesenteric adipose and the mammary gland. We cannot eliminate the possibility that the intestinal transcriptome could be differentially influenced by luminal vs. systemic administration of ergot alkaloids, however, previous studies have shown that subcutaneous injections of bromocriptine alter intestinal transit in mice () and both intravenous and intraluminal bromocriptine stimulated small intestinal absorption of nutrients in multiple mammalian species (). Based on selection criteria, RNA-Seq analysis revealed two DEG within the intestinal epithelium. These genes were fatty acid binding protein 4 (FABP4) which was down regulated and UDP-GlcNAc:betaGal Beta-1,3-N-Acetylglucosaminyltransferase 6 (B3GNT6) which was up regulated in steers treated with BROMO. Fatty acid binding protein is known to be involved in transport and metabolism of long chain fatty acids and overall lipid metabolism (, , ). Up regulation of B3GNT6 is indicative of increased lymphocyte trafficking and homing which may suggested increased immune activity as reported in other published data (). Interestingly, although a limited number of DEG were detected both fit with previously reported effects of ergot alkaloid exposure.

There were two canonical pathways (LXR/RXR activation, granulocyte adhesion and diapedesis pathways) depressed in both tissue types. These pathways provide further evidence that alkaloid exposure negatively impacts immune function. Granulocyte (neutrophils, basophils, and eosinophils) adhesion and diapedesis are more involved in the efficiency and activity of the immune response within these tissues. This is supported by depressed movement of neutrophil identified by individual function analysis within IPA. Most of the other intestine specific pathways influenced by alkaloid exposure have a role in the inflammatory response: the inflammasome pathway, HMGB1, and atherosclerosis signaling. Interestingly, HMGB1 signaling and monocyte infiltration can also influence expression of adhesion molecules (). While this was not evident in DEGs, upon further analysis we found both IL18 and C-C Motif Chemokine Ligand 2 (CCL2) tended to be down regulated and plasminogen activator (PLAT) tended to be up regulated. These genes were sufficient to identify the HMGB1 pathway as significantly affected and may suggest that BROMO treatment may impact barrier function of the intestine similar to results in mice ().

Cellular and physiological functions identified from these datasets were indicative of effects on cell movement, cell and tissue growth and development, and cell survival. Functional analysis using IPA indicated migration of smooth muscle cells was influenced in BROMO steers. This is meaningful in that intravenous administration of ergot alkaloids has also resulted in an immediate inhibition of muscle contractions in the forestomach of ruminants (). Thus, these effects could be due to direct actions of alkaloids on myenteric neurons and on smooth muscle cells (). The most significant physiological function identified was hematological system development and function. Initially this was thought to be changes in vascularity which would agree with previously published research (, , ). However, a closer look at individual functions revealed that circulating immune surveillance rather than angiogenesis was impacted within the intestinal epithelium. Because vasoactivity was not measured in the current study it is equivocal whether fescue alkaloid-induced vasoconstriction, i.e., predominantly ergovaline, fails to result in DEG or whether the dose of BROMO used in the current study was below the threshold to induce vasoconstriction in visceral blood vessels as previously reported. The latter can be supported by the fact that the dose of BROMO used in the current study had no impact on the ability of treated steers to consume 1.5 × net energy for maintenance daily. Whereas, previous studies using fescue seed have shown much more severe effects on feed intake in cattle (, , ).

Mesenteric Adipose

There was a greater number of DEG due to BROMO in the mesenteric adipose with 20 DEG (14 down regulated and 6 up regulated). Down regulated genes functions included enzymes, transporters, ion channel, cytokines, and other immune response molecules. The 4 DEG identified as an enzymes were Cytochrome P450 Family 4 Subfamily F Member 2 (CYP4F2), involved in fatty acid and cholesterol synthesis; Activation Induced Cytidine Deaminase (AICDA), regulation of immunoglobulin somatic hypermutation and class switch recombination (CSR), processes required for B-cell development of high affinity antibodies; cadherin 16 (CDH16) an adhesion molecule, and Glycosyltransferase 1 Domain Containing 1 (GLT1D1) for which an exact function is unknown. The transporters down regulated in adipose were both members of the solute carrier family, SLC9A2 and SLC26A5, and are involved in cation and anion transport, respectively. In addition to SLC9A2, which is a sodium/hydrogen exchanger, Acid Sensing Ion Channel Subunit Family Member 5 (ASIC5) a Na+ ion channel was also down regulated. The down regulation of SLC9A2 which functions to counter adverse environmental conditions, specifically pH regulation is intriguing but the exact role cannot be determined in the current dataset.

Maruo et al. () suggested activation of a negative feedback mechanism in the intestine and liver of pigs fed ergot alkaloids which also had a down regulation of toll-like receptors and a variety for cytokines including IL-6 and Tumor Necrosis Factor α (TNFα). In pigs fed ergot alkaloids tissues may be attempting to restore barrier function by up regulating mRNA expression of adhesion proteins such as occludins, claudins 3 and 4, and others (). Accordingly, the down regulation of most genes in the current study may exemplify this purposed feedback mechanism within the adipose and reflect the drive to return to normal physiological conditions. In contrast, alkaloid ingestion has been reported to cause increased inflammation () and decrease overall immune competency (). These conflicting results were also apparent in our data. However, an explanation is not immediately apparent but there are still clear alterations in immunological response and function.

The six up regulated DEG had identified functions of enzyme activity [Monooxygenase DBH Like 1 (MOXD1)], molecular transport [Sortilin Related VPS10 Domain Containing Receptor 2 (SORCS2); Solute Carrier Family 11 Member 1 (SLC11A1)], transmembrane receptor [Delta/Notch Like EGF Repeat (DNER)], cytokine [Platelet Factor 4 (PF4)] and other unknown functions [Lymphocyte Antigen 6 Family Member G5B (LY6G5B)]. In contrast, mRNA and protein expression in the intestine and liver of swine fed ergot alkaloids exhibited up regulation of adhesion molecules (); however, CSH16 was the only DEG identified that has adhesion properties and it was down regulated. This may be a differential effect of alkaloid exposure between species but may simply confirm the wide variety of physiological functions influenced by ergot alkaloids.

The canonical pathways influenced by BROMO treatment were largely focused on inflammation, immune response, and lipid metabolism. Pathway identification from IPA utilized all significantly expressed genes (P < 0.05) with and without a fold change >1.5. The LXR/RXR activation and granulocyte adhesion and diapedesis pathways were depressed in the mesenteric adipose as well as the intestine. Inhibition of the LXR/RXR heterodimer would directly modulate the initiation of immune and inflammatory responses in macrophages (). Individual functions support the influence on macrophage function from recruitment and movement to proliferation and overall quantity. Whereas, granulocyte (neutrophils, basophils, and eosinophils) adhesion and diapedesis are more involved in the efficiency and activity of the immune response within these tissues. Unlike in the intestine, overall granulocyte quantity and recruitment were depressed in adipose tissue of these steers rather than neutrophils specifically. Genes involved with the compliment system (SERPINB2 and C8A), chemokine secretion (CXCL2 and PF4), and immunoglobulin function (AICDA and HEPACAM2) were all down regulated DEG except for PF4 which was one of a fed up regulated DEG in adipose tissue. Further work needs to be completed to fully understand the role and mechanisms of DEG within these pathways and how ergot alkaloid exposure changes the function of these pathways and overall physiology of the animal.

Accumulation of ergot alkaloids in adipose tissue has been reported to disrupt or alter lipid metabolism (). Both CYP4F2 and NR1H4 were down regulated in the current dataset. Lipid and bile acid synthesis were depressed with down regulation of CYP4F2 and NR1H4 and further supported by identification of HDL and LDL inhibition within LXR/RXR pathways. These data begin to identify specific genes and a putative mechanism to explain how alkaloids, specifically BROMO, influence adipocytes and lipid metabolism.

LPS/IL-1 mediated inhibition of RXR and IL-6 signaling were activated and FXR/RXR activation, antigen presentation, and dendritic cell maturation were depressed by alkaloid exposure. Canonical pathways influenced by BROMO treatment, like DEG, indicate opposing effects on immune response. Specifically, the influence of ergot alkaloids on immune cell trafficking, cellular movement, and immune response within abdominal tissues. Previously reported data demonstrated that ergot alkaloids increased inflammation but decreased overall immune function in steers that grazed endophyte-infected tall fescue as evidenced by lower Cu status, depressed immune-competency, decreased major histocompatibility complex class II expression, and macrophage phagocytosis activity (). The current dataset agrees with previous literature () that exposed cattle could have a decreased immune response through inhibition of associated genes and other pathways associated with a variety of mechanisms. However, our data doesn't supply any further evidence for increased inflammation (, ). This dichotomy in the literature is similar to the conflicting DEG, pathways, and functions observed in the current dataset and thus requires more work to completely elucidate the mechanisms behind this contradiction.

In the current dataset, the serotonin degradation pathway (Table 2) was depressed by BROMO treatment in the adipose which is consistent with the known interaction between bromocriptine () or alkaloids () and serotonin metabolism. Specifically, ergot alkaloids have been reported to interact with serotonin receptors in vascular smooth muscle (, , ). Expression of the serotonin receptor 5HTR2A and 5HTR4 in gastrointestinal smooth muscle was reduced in steers treated with ergot alkaloids compared with control steers (). Serotonin receptor down regulation is likely due to alkaloids binding to the receptor. Bromocriptine has a strong affinity for the 5HTR2A receptor (). Therefore, BROMO is likely inhibiting the serotonin degradation pathway through antagonistic activity within the mesenteric adipose after exposure to BROMO. These data and that from Klotz et al. () suggest an interactive role of the adipose, intestine, and ergot alkaloids throughout the serotonin pathway from ligand to receptor.

Surprisingly, there were limited pathways identified by the current dataset specifically involved with angiogenesis or vascularity. As discussed earlier, alkaloid exposure has been reported to constrict blood vessels () and flow (, ) in other datasets. None of the DEGs identified had a direct effect on vascularity but CD8A, CXCL2, C8A, SERPINB2, AICDA, HEPCAM2, and PF4 are directly involved in immune system response and function. This may suggest that abdominal effects, such as necrotic fat, are a result of persistent activation of the immune system vs. vasoconstriction as seen in peripheral tissues (, , , ).

In the mesenteric adipose compared with the intestinal epithelium there was a greater number and variety of functions influenced by BROMO treatment; however most functions identified represented cell movement and survival, immune cell trafficking, tissue/organismal development, morphology, and survival. Mulac and Humpf () demonstrated apoptotic effects of ergot alkaloids on human primary cells in culture which may help to explain the inhibition of cell survival. Functions related to overall quantity, recruitment, and movement of macrophages along with other leukocytes were identified by functional analysis in the current dataset. Macrophages within adipose tissue are generally found scattered or in crown-like structures around dead adipocytes (). These crown-like structures were differentially expressed in abdominal fat depots compared with subcutaneous and pelvic adipose of obese () or leptin and leptin receptor deficient mice (). These structures within adipose may be a mechanism responsible for necrotic fat depot development but more likely are functioning to clear necrotic adipocytes. In the current dataset, more work needs to be completed to know if the macrophages were in crown-like structures associated with deceased adipocytes and if this a possible mechanism in the formation of necrotic fat depots rather than vasoconstriction.

The influence of BROMO on lipid metabolism is likely due, indirectly, to the decrease in prolactin (). Connor et al. () found exposure to shortened day lengths decreased circulating prolactin which led to decreased acyl-CoA dehydrogenase. As the first step in β-oxidation acyl-CoA dehydrogenase would reduce fatty acids () but have also been associated with increased hepatic triacylglycerol () and non-esterified fatty acids (, ). Alterations in lipid metabolism due to bromoergocriptine did not impact milk fat in lactating dairy cows (). Thus, indicating prolactin alters fatty acid () and cholesterol (, ) metabolism. However, in addition to those mediated through prolactin, it is possible that BROMO treatment and ergot alkaloids have a more direct effect on lipid metabolism. These data also indicate BROMO can alter adipose metabolism and function which would change the morphology of adipose tissue and influence immune surveillance by leukocytes located within this tissue.

In summary, exposure of Holstein steers to systemic bromocriptine resulted in decreased plasma prolactin concentration which is what occurs when steers graze endophyte-infected tall fescue or are administered ergot alkaloids. RNA-Seq analysis revealed that exposure to a synthetic alkaloid was influential on pathways and functions involved with immune cell trafficking, inflammatory responses, and lipid metabolism compared with other physiological pathways and functions. These data provide insight into the specific tissue responses and potential pathways and genes affected and thereby begin the process of elucidating the mechanisms behind compromised health and productivity of animals grazing ergot alkaloids. However, further and more specific analyses need to be conducted to completely understand the mechanisms behind the negative effects of animals exposed to ergot alkaloids.

Statements

Data availability statement

The raw data reads were lost and are not available for upload, however the extended and comprehensive findings have been provided in the supplemental tables. Any further information that we can provide will be available upon discussion with the authors.

Ethics statement

The animal study was reviewed and approved by University of Kentucky Animal Care and Use Committee.

Author contributions

KMcLea and KMcLeo performed experiment and conception and design of research. KMcLea, RB, CL, JE, and KMcLeo analyzed data. KMcLea, RB, CL, and KMcLeo interpreted results of experiments. KMcLea prepared figures and drafted manuscript. KMcLea, RB, CL, JK, and KMcLeo edited and revised manuscript. KMcLea, RB, CL, JK, JE, and KMcLeo approved final version of manuscript. All authors contributed to the article and approved the submitted version.

Funding

The information reported in this paper (No. 18-01-077) is part of a project of the Kentucky Agricultural Experiment Station and is published with the approval of the Director. Mention of trade name, proprietary product, or specified equipment does not constitute a guarantee or warranty by the University of Kentucky and does not imply approval to the exclusion of other products that may be available.

Acknowledgments

Researchers would like to thank Kirk Vanzant, Lauren Clark, and the rest of the beef unit staff, Brock Billingsley, and Adam Bohannon for their support during sample collection and completion of this project. We would also like to extend our appreciation to Rebecca Payton for running plasma samples for prolactin quantification.

Conflict of interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Supplementary material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fvets.2020.00615/full#supplementary-material

References

Summary

Keywords

adipose, alkaloids, intestine, ruminants, transcriptome

Citation

McLean KJ, Baldwin VI RL, Li C, Klotz JL, Edwards JL and McLeod KR (2020) Synthetic Alkaloid Treatment Influences the Intestinal Epithelium and Mesenteric Adipose Transcriptome in Holstein Steers. Front. Vet. Sci. 7:615. doi: 10.3389/fvets.2020.00615

Received

02 October 2019

Accepted

29 July 2020

Published

11 September 2020

Volume

7 - 2020

Edited by

Nguyen Hong Nguyen, University of the Sunshine Coast, Australia

Reviewed by

Yuanning Li, Yale University, United States; Stephen Brent Smith, Texas A&M University, United States

Updates

Copyright

*Correspondence: Kyle R. McLeod

This article was submitted to Livestock Genomics, a section of the journal Frontiers in Veterinary Science

Disclaimer

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.

Outline

Figures

Cite article

Copy to clipboard


Export citation file


Share article

Article metrics