Effect of fescue toxicosis on hepatic gene expression in mice1

Gene expression profiling indicates an increased capacity for proline, serine, and ATP synthesis and mitochondrial mass by the liver of steers grazing high vs. low endophyte-infected tall fescue

Grazing -infected forages results in a variety of reduced animal performance parameters, collectively known as "fescue toxicosis." The initial, limited evaluations of hepatic mechanisms affected by fescue toxicosis have used transcriptomic expression profiling of experimental phenotypes developed by short-term feeding of concentrated ergot alkaloids or fescue seeds to rodents and steers. To assess the effects of fescue toxicosis in growing cattle using a commercially relevant phenotype, we induced fescue toxicosis in beef steers by summer-long grazing (89 to 105 d) of a single high toxic endophyte-infected tall fescue pasture (HE; 0.746 μg/g ergot alkaloids; 5.7 ha; = 10; BW = 267 ± 14.5 kg) vs. a low toxic endophyte tall fescue-mixed pasture (LE; 0.023 μg/g ergot alkaloids; 5.7 ha; = 9; BW = 266 ± 10.9 kg). High toxic endophyte tall fescue-mixed pasture steers had decreased BW (313 vs. 338 kg) and an increased potential for hepatic gluconeogenesis from AA-derived carbons. To gain a greater perspective into fescue toxicosis-induced hepatic metabolism and identify candidate regulatory mechanisms, the goal of the current research was to examine liver samples for changes in gene (mRNA) expression profiles using a Bovine Affymetrix microarray and selected reverse-transcription PCR and immunoblot analyses. The expression (false discovery rate < 10%; < 0.01) of 147 genes was increased (7 to 268%) and that of 227 was decreased (4 to 87%) in livers of HE vs. LE steers. The top (1) functional gene category was cell-mediated immune response (33 genes; ≤ 0.012), (2) canonical cell signaling pathway was primary immunodeficiency signaling (8 genes; ≤ 0.0003), and (3) canonical metabolic pathways were oxidative phosphorylation (5 genes; ≤ 0.016) and purine metabolism (8 genes; ≤ 0.029). High toxic endophyte tall fescue-mixed pasture steers had increased ( ≤ 0.022) expression of genes critical for increased (1) Pro () and Ser () synthesis, (2) shunting of AA carbons into pyruvate () and ATP synthesis (, , , COX4, , and ), and (3) mitochondrial mass (COX4). Targeted reverse-transcribed PCR or immunoblot assays corroborated ( ≤ 0.035) these latter microarray findings for , , , , and COX4. Moreover, network analysis identified glucocorticoid receptor-mediated signaling as the most probable mechanism to coordinate the above findings. These results greatly extend our knowledge of the consequences of summer-long grazing of endophyte-infected tall fescue to the hepatic metabolism of growing steers.

NutriPhysioGenomics applications to identify adaptations of cattle to consumption of ergot alkaloids and inorganic versus organic forms of selenium: altered nutritional, physiological and health states?

NutriPhysioGenomics (NPG) seeks to measure nutrition-responsive genome expression during specific physiological states, thus defining how a given challenge alters the ‘basal’ transcriptome. Application of NGS regimens (transcriptome and bioinformatics analyses) in combination with targeted-gene approaches has revealed cellular mechanisms putatively responsible for complex, whole-animal, metabolic syndromes such as heat stress and subacute ruminal acidosis. Using similar approaches, our laboratory sought to understand how the basal physiology of developing cattle adapted to two prevalent health challenges of forage-based beef cattle production in south-eastern USA: fescue toxicosis and selenium (Se) deficiency. In Model 1, pituitary and hepatic genomic expression profiles of growing beef steers grazing high (E+, n = 9) versus low (E–, n = 10) endophyte-infected tall fescue pastures for 85 days, and consuming sodium selenite (ISe) as a free-choice Se supplement, were compared by a combination of microarray, bioinformatic, and targeted-gene/protein (real-time reverse transcription–PCR, Nanostring, immunoblot) analyses. In Model 2, hepatic genomic expression profiles of growing beef heifers (0.5 kg gain/day) fed a cotton seed hull-based diet and different sources (n = 9) of dietary Se (3 mg/day) supplements (no supplement, Control; inorganic Se, sodium selenite, ISe; organic Se, Sel-Plex; OSe), or a 1.5 mg : 1.5 mg blend of ISe and OSe, MIX) were compared after 168 days of supplementation, as described for Model 1. The results for Model 1 showed, that in the pituitary of E+ steers, expression of genes for prolactin signalling; redox capacity; regulation of lactotroph, gonadotroph, and thyrotroph proliferation; gonadotropin-releasing hormone-mediated signalling; and Se-based metabolism was impaired. Concomitantly, the livers of E+ steers had an increased level of expression of genes encoding proteins responsible for shunting of amino acid carbons into pyruvate and ATP synthesis capacity (oxidative phosphorylation pathway, mitochondrial mass), increased serine and proline biosynthesis, and reduced selenoprotein-mediated metabolism. Result for Model 2 showed that, overall, there were clear differences in the profiles of differentially expressed genes (DEG) among the four Se treatment groups, with the form of Se administered being more reflective of DEG profiles than the total amount of Se assimilated. Moreover, hepatic transcriptomes profiles of MIX heifers revealed an increased potential for selenoprotein synthesis and selenoprotein-mediated metabolism. In addition, several genes involved with increased redox capacity were upregulated in MIX versus ISe heifers. Taken together, our NGS approach characterised adaptation to physiological challenges and, serendipitously, identified suppression of several metabolic pathways by consumption of ergot alkaloid consumption that have the potential to be increased with supplementation of the MIX form of Se.

Microarray evaluation of bovine hepatic gene response to fescue toxicosis

‘Fescue toxicosis’ is a disease in livestock caused by ingestion of ergot alkaloids produced by the fungal endophyte Neotyphodium coenophialum in tall fescue; it is estimated to cost 1 billion USD in damages per year to the beef industry alone. Clinical signs include decreased reproductive fitness, necrosis of extremities, and reduced average daily gain and milk production. Little is known about the cellular mechanisms that mediate these toxic sequelae. We evaluated the effects of ergovaline-based fescue toxicosis on gene expression via oligonucleotide microarray. Liver biopsies were obtained from steers (n=4) pre- and post-exposure (0 and 29 days) to feed containing 579 ng/g ergovaline. Analyses were performed using both ANOVA with false discovery rate correction and Storey's optimal discovery procedure. Overall, down-regulation of gene expression was observed; heart contraction and cardiac development, apoptosis, cell cycle control, and RNA processing genes represented the bulk of differentially expressed transcripts. 2 CYPs (CYP2E1 and CYP4F6) were amongst the significantly upregulated results. Thus, exposure of cattle to toxic levels of ergovaline caused widespread changes in hepatic gene expression, which can both help explain macroscopic clinical signs observed in ruminant animals, and reinforce previous findings in monogastric models.

Interpretation of ANOVA models for microarray data using PCA

MOTIVATION

ANOVA is a technique, which is frequently used in the analysis of microarray data, e.g. to assess the significance of treatment effects, and to select interesting genes based on P-values. However, it does not give information about what exactly is causing the effect. Our purpose is to improve the interpretation of the results from ANOVA on large microarray datasets, by applying PCA on the individual variance components. Interaction effects can be visualized by biplots, showing genes and variables in one plot, providing insight in the effect of e.g. treatment or time on gene expression. Because ANOVA has removed uninteresting sources of variance, the results are much more interpretable than without ANOVA. Moreover, the combination of ANOVA and PCA provides a simple way to select genes, based on the interactions of interest.

RESULTS

It is shown that the components from an ANOVA model can be summarized and visualized with PCA, which improves the interpretability of the models. The method is applied to a real time-course gene expression dataset of mesenchymal stem cells. The dataset was designed to investigate the effect of different treatments on osteogenesis. The biplots generated with the algorithm give specific information about the effects of specific treatments on genes over time. These results are in agreement with the literature. The biological validation with GO annotation from the genes present in the selections shows that biologically relevant groups of genes are selected.

AVAILABILITY

R code with the implementation of the method for this dataset is available from http://www.cac.science.ru.nl under the heading "Software".

Expression profiling of heat stress effects on mice fed ergot alkaloids

Fescue toxicosis affects wild and domestic animals consuming ergot alkaloids contained in tall fescue forage infected with the endophytic fungus, Neotyphodium coenophialum. When animals are consuming infected fescue (E+) forage during periods of elevated ambient temperatures (summer), a range of phenotypic disorders collectively called summer slump is observed. It is characterized by hyperthermia, with an accompanying decrease in feed intake, growth, milk yield, and reproductive fitness. Laboratory mice also exhibit symptoms of fescue toxicosis at thermoneutral (TN) temperature, as indicated by reduced growth rate and reproductive fitness. Our goal was to characterize the differences in gene expression in liver of mice exposed to summer-type heat stress (HS) and E+ when compared to mice fed E+ at TN temperature. Mice were fed E+ diet under HS (34 +/- 1 degrees C; n = 13; E+HS) or TN conditions (24 +/- 1 degrees C; n = 14; E+TN) for a period of 2 weeks between 47 and 60 days of age. Genes differentially expressed between E+HS versus E+TN were identified using DNA microarrays. Forty-one genes were differentially expressed between treatment groups. Expressions of eight genes were measured using quantitative real-time PCR. Genes coding for phase I detoxification enzymes were upregulated in E+HS mouse liver. This detoxification pathway is known to produce reactive oxidative species. We observed an upregulation of genes involved in the protection against reactive oxidative species. Key genes involved in de novo lipogenesis and lipid transport were also upregulated. Finally, genes involved in DNA damage control and unfolded protein responses were downregulated.