Transcriptomic analysis of skeletal muscle from beef cattle exposed to illicit schedules containing dexamethasone: identification of new candidate biomarkers and their validation using samples from a field monitoring trial

The Near-Infrared Spectroscopy of Ethanol-Fixed Tissues to Detect Illicit Treatments with Glucocorticoids in Bulls

This study aimed to set up indirect, rapid methods involving near infrared (NIR) spectroscopy analysis, to detect illicit treatments with glucocorticoids in bull. The ethanol fixation method (EtOH) was applied to 7 different tissues obtained from 20 Friesian bulls, 12 of which were experimentally administered with dexamethasone as part of a growth-promoting protocol for 60 days and slaughtered 26 days after the end of the treatment. A perfect discrimination was obtained for the 7 sampled tissues, considering a full UV-Vis-NIR range (350 ÷ 2500 nm), for both false positive and negative animals. The validated true positive and negative errors were zero for the longissimus thoracis muscle, 10% for the skin-dermis, 15% for the fat, 25% for the thymus gland and the semitendinosus muscle, 30% for the sternomandibularis muscle and 35% for the skin-hair. A multiple test on the most accessible tissues, that is, the thymus gland, the sternomandibularis muscle and fat, can be used as an alternative to provide indications about animals that have been subjected to illicit treatments. In the short space of three days from the slaughter, NIR spectroscopy of ETOH fixed tissues, would allow at least cost the detection of a probable illicit which could eventually be reported to health authorities for specific investigation in the frame of official controls.

Transcriptional Biomarkers and Immunohistochemistry for Detection of Illicit Dexamethasone Administration in Veal Calves

Corticosteroids such as Dexamethasone (DEX) are commonly licensed for therapy in meat animals due to their known pharmacological properties. However, their misuse aimed to achieve anabolic effects is often found by National Residues Control Plans. The setup of a complementary “biomarker based” methods to unveil such illicit practices is encouraged by current European legislation. In this study, the combined use of molecular and histological quantitative techniques was applied on formalin fixed paraffin embedded (FFPE) muscle samples to assess the effects of illicit DEX treatment on veal calves. A PCR array, including 28 transcriptional biomarkers related to DEX exposure, was combined with a histochemical analysis of muscle fiber. An analysis based on unsupervised (PCA) and supervised (PLS-DA and Kohonen’s SOM) methods, was applied in order to define multivariate models able to classify animals suspected of illicit treatment by DEX. According to the conventional univariate approach, a not-significant reduction in type I fibres was recorded in the DEX-treated group, and only 12 out of 28 targeted genes maintained their expected differential expression, confirming the technical limitations of a quantitative analysis on FFPE samples. However, the multivariate models developed highlighted the possibility to establish complementary screening strategies, particularly when based on transcriptional biomarkers characterised by low expression profiles.

FKBP5 gene expression in skeletal muscle as a potential biomarker for illegal glucocorticoid treatment in veal calves

For the current European legislation, the chemical analysis of drug residues is the exclusive accepted method to identify animals illicitly treated with growth promoters. Glucocorticoids and their metabolites are no detectable by LC/MS-MS methods in biological fluids when the growth promoter administration is discontinued several days prior to the slaughtering. The aim of this study was to elucidate the effect on the expression of genes belonging to the glucocorticoid pathway in three types of skeletal muscle of calves treated with prednisolone or dexamethasone in combination with estradiol. A gene expression change of glucocorticoid receptors (NR3C1 and NR3C2), their chaperones molecules (FKBP prolyl isomerase 4 and 5, FKBP4 and 5) and pre-receptor system (hydroxysteroid 11-beta dehydrogenases 1 and 2, HSD11B1 and 2) may indicate potential biomarkers of glucocorticoid treatment. In the biceps brachii muscle, the administration of dexamethasone with estradiol increased HSD11B2 (P < 0.01) and NR3C2 (P < 0.01) gene expression, whereas prednisolone administration increased HSD11B1 transcript levels (P < 0.05). In the longissimus lumborum muscle, NR3C2 gene expression decreased following prednisolone administration (P < 0.05). FKBP5 gene expression decreased in all considered muscles of calves administered with dexamethasone and estradiol (P < 0.01), whereas increased in the longissimus lumborum (P < 0.01) and vastus lateralis (P < 0.05) muscle of prednisolone-treated group (P < 0.05). The opposite effect of dexamethasone and prednisolone appears very promising to develop a low-cost screening test, because the expression analysis of a unique gene in a given tissue may distinguish the dispensed molecules.

Expression of corticosteroid hormone receptors, prereceptors, and molecular chaperones in hypothalamic-pituitary-adrenal axis and adipose tissue after the administration of growth promoters in veal calves

The action of glucocorticoids on target tissues is regulated by the glucocorticoid and mineralocorticoid receptors (codified by the NR3C1 and NR3C2 gene, respectively). Moreover, the prereceptor system, represented by the hydroxysteroid 11-beta dehydrogenases (HSD11Bs), catalyzes the interconversion from active glucocorticoids into inactive compounds. This study aimed to determine whether the expression of the prereceptor system, the corticosteroid receptors, and the molecules regulating their intracellular trafficking (FKBP prolyl isomerase 4 and FKBP prolyl isomerase 5) could be regulated in the hypothalamic-pituitary-adrenal axis and in different type of adipose tissue of calves by the administration of dexamethasone in combination with estradiol or prednisolone. Research about the glucocorticoid effects on bovine target tissues may allow development of new diagnostic methods that use potential molecular biomarkers of glucocorticoid treatment. The administration of dexamethasone in combination with estradiol increased the gene expression of HSD11B1 (P < 0.01), HSD11B2 (P < 0.05), NR3C1 (P < 0.01), and NR3C2 (P < 0.01) in the adrenal glands; NR3C2 in the intramuscular adipose tissue (P < 0.01), and HSD11B1 in the subcutaneous adipose tissue (P < 0.01). Prednisolone administration increased the gene expression of HSD11B1 (P < 0.01), NR3C1 (P < 0.05), and NR3C2 (P < 0.05) in the adrenal glands and HSD11B1 (P < 0.01) in the subcutaneous adipose tissue. Interestingly, most of the examined tissues/organs showed a significant variation of FKBP5 gene expression after the administration of dexamethasone in combination with estradiol. So, these changes suggest that the FKBP5 gene expression could be a possible biomarker of the illegal dexamethasone administration in calves.

Role of FKBP51 in the modulation of the expression of the corticosteroid receptors in bovine thymus following glucocorticoid administration

The aim of this work was to study the transcriptional effects of glucocorticoids on corticosteroid hormone receptors, prereceptors (11β-hydroxysteroid dehydrogenase 1 and 2, 11β-HSD1 and 2), and chaperones molecules regulating intracellular trafficking of the receptors (FKBP51 and FKBP52) in thymus of veal calves. Moreover, the expression of FKBP51 and FKBP52 gene were investigated in beef cattle thymus. In the cervical thymus of veal calves, dexamethasone administration in combination with estradiol decreased FKBP51 expression (P < 0.01). The same treatment increased mineralocorticoid receptor (MR) (P < 0.01) and 11β-HSD1 expression (P < 0.05) compared to control group in the cervical thymus of veal calves. The thoracic thymus of veal calves treated with dexamethasone and estradiol showed a decrease of FKBP51 (P < 0.05), FKBP52 (P < 0.05), glucocorticoid receptor (P < 0.05), and MR expression (P < 0.05) compared to control group in the thoracic thymus of veal calves. The gene expression of FKBP51 decreased both in cervical (P < 0.01) and thoracic thymus (P < 0.01) of beef cattle treated with dexamethasone and estradiol. In addition, also prednisolone administration reduced FKBP51 expression in the cervical thymus (P < 0.01) and in the thoracic thymus of beef cattle (P < 0.01). The gene expression of FKBP52 increased only in the cervical thymus following dexamethasone administration (P < 0.01). The decrease of FKBP51 gene expression in thymus could be a possible biomarker of illicit dexamethasone administration in bovine husbandry. Moreover, so far, an effective biomarker of prednisolone administration is not identified. In this context, the decrease of FKBP51 gene expression in thymus of beef cattle following prednisolone administration could play an important role in the indirect identification of animals illegally treated with prednisolone.

Transcriptomic characterization of bovine primary cultured hepatocytes; a cross-comparison with a bovine liver and the Madin-Darby bovine kidney cells

Bovine primary cultured hepatocytes (CHs) are widely used in vitro models for liver toxicity testing. However, little is known about their whole-transcriptome profile and its resemblance to the normal liver tissue. In the present study, we profiled - by microarray - the whole-transcriptome of bovine CHs (n=4) and compared it with the transcriptomic landscape of control liver samples (n=8), as well the Madin-Darby bovine kidney (MDBK) cells (n=4). Compared with liver tissue, the bovine CHs relatively expressed (fold change >2, P<0.05) about 2155 and 2073 transcripts at a lower and higher abundance, respectively. Of those expressed at a lower abundance, many were drug biotransformation enzyme-coding genes, such as the cytochrome P450 family (CYPs), sulfotransferases, methyltransferases, and glutathione S-transferases. Also, several drug transporters and solute carriers were expressed at a lower abundance in bovine CHs. 'Drug metabolism', 'PPAR signaling', and 'metabolism of xenobiotics by CYPs' were among the most negatively-enriched pathways in bovine CHs compared with liver. A qPCR cross-validation using 8 selected genes evidenced a high correlation (r=0.95, P=0.001) with the corresponding microarray results. Although from a kidney origin, and albeit to a lower extent compared to bovine CHs, the MDBK cells showed a basal expression of many CYP-coding genes. Our study provides a whole-transcriptome-based evidence for the bovine CHs and hepatic tissue resemblance. Overall, the bovine CHs' transcriptomic profile might render it unreliable as an in vitro model to study drug metabolism.

Transcriptomic signature of high dietary organic selenium supplementation in sheep: A nutrigenomic insight using a custom microarray platform and gene set enrichment analysis

The objective of this study was to investigate the effect of a high dietary Se supplementation on the whole transcriptome of sheep. A custom sheep whole-transcriptome microarray, with more than 23,000 unique transcripts, was designed and then used to profile the global gene expression of sheep after feeding a high dietary supplementation of organic Se. Lactating crossbred ewes ( = 10; 3 to 4 yr of age and 55 to 65 kg BW) at late lactation (100 ± 8 d in milk) were acclimated to indoor individual pen feeding of a basal control diet (0.40 mg Se/d, sodium selenite) for 4 wk. Sheep were then kept on a diet with an extra (high) supplementation of organic Se (1.45 mg Se/d as Sel-Plex; Alltech Biotechnology Pty Ltd, Dandenong, Victoria, Australia) for 40 d. Whole blood was collected at 2 time points (last day of the acclimatization period [T0] and after 40 d of the organic Se supplementation [T40]), and then total RNA was isolated and labeled for the subsequent microarray analysis. Significance Analysis of Microarrays, using the -statistic, of the microarray data (T40 versus T0) evidenced the up- and downregulation of 942 and 244 transcripts (false discovery rate < 0.05), respectively. Seven genes showed the same trend of expression (up- or downregulation) when tested by quantitative real-time PCR (qPCR) in a cross-validation step. The microarray showed significant upregulation of the following selenoproteins at T40: selenium binding protein 1 (SELENBP1), selenoprotein W1 (SEPW1), glutathione peroxidase 3 (GPX3), and septin 8 (SEPT8). And the expression trends for SEPW1 and SEPT8 were validated using qPCR. Functional annotation of the differentially expressed genes showed the enrichment of several immune system-related biological processes (lymphocyte activation, cytokine binding, leukocyte activation, T cell differentiation, and B cell activation) and pathways (cytokine and interleukin signaling). Moreover, Gene Set Enrichment Analysis evidenced the enrichment of B and T cell receptors signaling pathways, with an enrichment score of 0.63 and 0.59, respectively. Overall, from a global gene expression (whole-transcriptome) point of view, short-term supplementation of a high dietary organic Se to Se-nondeficient sheep results in a transcriptomic signature that mainly reflects an induced immune system and a modulation of transcription effect. Also, the present study provides a custom whole-transcriptome microarray platform that can be used in further global gene expression studies in the ovine species.