Search

SummaryWhile over ten thousand genetic loci have been associated with phenotypic traits and inherited diseases in genome-wide association studies, in most cases only a relatively small proportion of the trait heritability is explained and biological mechanisms underpinning these traits have not been clearly identified. Expression quantitative trait loci (eQTL) are subsets of genomic loci shown experimentally to influence gene expression. Since gene expression is one of the primary determinants of phenotype, the identification of eQTL may reveal biologically relevant loci and provide functional links between genomic variants, gene expression and ultimately phenotype. Skeletal muscle (gluteus medius) gene expression was quantified by RNA-seq for 111 Thoroughbreds (47 male, 64 female) in race training at a single training establishment sampled at two time-points: at rest (n = 92) and four hours after high-intensity exercise (n = 77); n = 60 were sampled at both time points. Genotypes were generated from the Illumina Equine SNP70 BeadChip. Applying a Benjamini-Hochberg corrected P-value threshold (PFDR < 0.05), association tests identified 3,583 cis-eQTL associated with expression of 1,456 genes at rest; 4,992 cis-eQTL associated with the expression of 1,922 genes post-exercise; 1,703 trans-eQTL associated with 563 genes at rest; and 1,219 trans-eQTL associated with 425 genes post-exercise. The gene with the highest cis-eQTL association at both time-points was the endosome-associated-trafficking regulator 1 gene (ENTR1; Rest: PFDR = 3.81 × 10−27, Post-exercise: PFDR = 1.66 × 10−24), which has a potential role in the transcriptional regulation of the solute carrier family 2 member 1 glucose transporter protein (SLC2A1). Functional analysis of genes with significant eQTL revealed significant enrichment for cofactor metabolic processes. These results suggest heritable variation in genomic elements such as regulatory sequences (e.g. gene promoters, enhancers, silencers), microRNA and transcription factor genes, which are associated with metabolic function and may have roles in determining end-point muscle and athletic performance phenotypes in Thoroughbred horses.

BackgroundBovine tuberculosis is caused by infection withMycobacterium bovis, which can also cause disease in a range of other mammals, including humans. Alveolar macrophages are the key immune effector cells that first encounterM. bovisand how the macrophage epigenome responds to mycobacterial pathogens is currently not well understood.ResultsHere, we have used chromatin immunoprecipitation sequencing (ChIP-seq), RNA-seq and miRNA-seq to examine the effect ofM. bovisinfection on the bovine alveolar macrophage (bAM) epigenome. We show that H3K4me3 is more prevalent, at a genome-wide level, in chromatin fromM. bovis-infected bAM compared to control non-infected bAM; this was particularly evident at the transcriptional start sites of genes that determine programmed macrophage responses to mycobacterial infection (e.g. M1/M2 macrophage polarisation). This pattern was also supported by the distribution of RNA Polymerase II (PolII) ChIP-seq results, which highlighted significantly increased transcriptional activity at genes demarcated by permissive chromatin. Identification of these genes enabled integration of high-density GWAS data, which revealed genomic regions associated with resilience to infection withM. bovisin cattle.ConclusionsThrough integration of these data, we show that bAM transcriptional reprogramming occurs through differential distribution of H3K4me3 and PolII at key immune genes. Furthermore, this subset of genes can be used to prioritise genomic variants from a relevant GWAS data set.

RNA-seq has emerged as an important technology for measuring gene expression in peripheral blood samples collected from humans and other vertebrate species. In particular, transcriptomics analyses of whole blood can be used to study immunobiology and develop novel biomarkers of infectious disease. However, an obstacle to these methods in many mammalian species is the presence of reticulocyte-derived globin mRNAs in large quantities, which can complicate RNA-seq library sequencing and impede detection of other mRNA transcripts. A range of supplementary procedures for targeted depletion of globin transcripts have, therefore, been developed to alleviate this problem. Here, we use comparative analyses of RNA-seq data sets generated from human, porcine, equine and bovine peripheral blood to systematically assess the impact of globin mRNA on routine transcriptome profiling of whole blood in cattle and horses. The results of these analyses demonstrate that total RNA isolated from equine and bovine peripheral blood contains very low levels of globin mRNA transcripts, thereby negating the need for globin depletion and greatly simplifying blood-based transcriptomic studies in these two domestic species.

microRNAs (miRNAs) are a class of small non-coding endogenous RNA molecules that regulate a wide range of biological processes by post-transcriptionally regulating gene expression. Thousands of these molecules have been discovered to date, and multiple miRNAs have been shown to coordinately fine-tune cellular processes key to organismal development, homeostasis, neurobiology, immunobiology, and control of infection. The fundamental regulatory role of miRNAs in a variety of biological processes suggests that differential expression of these transcripts may be exploited as a novel source of molecular biomarkers for many different disease pathologies or abnormalities. This has been emphasized by the recent discovery of remarkably stable miRNAs in mammalian biofluids, which may originate from intracellular processes elsewhere in the body. The potential of circulating miRNAs as biomarkers of disease has mainly been demonstrated for various types of cancer. More recently, however, attention has focused on the use of circulating miRNAs as diagnostic/prognostic biomarkers of infectious disease; for example, human tuberculosis caused by infection with Mycobacterium tuberculosis, sepsis caused by multiple infectious agents, and viral hepatitis. Here, we review these developments and discuss prospects and challenges for translating circulating miRNA into novel diagnostics for infectious disease.

MicroRNAs (miRNAs) are important regulators of gene expression and are known to play a key role in regulating both adaptive and innate immunity. Bovine alveolar macrophages (BAMs) help maintain lung homeostasis and constitute the front line of host defense against several infectious respiratory diseases, such as bovine tuberculosis. Little is known, however, about the role miRNAs play in these cells. In this study, we used a high-throughput sequencing approach, RNA-seq, to determine the expression levels of known and novel miRNAs in unchallenged BAMs isolated from lung lavages of eight different healthy Holstein-Friesian male calves. Approximately 80 million sequence reads were generated from eight BAM miRNA Illumina sequencing libraries, and 80 miRNAs were identified as being expressed in BAMs at a threshold of at least 100 reads per million (RPM). The expression levels of miRNAs varied over a large dynamic range, with a few miRNAs expressed at very high levels (up to 800,000RPM), and the majority lowly expressed. Notably, many of the most highly expressed miRNAs in BAMs have known roles in regulating immunity in other species (e.g. bta-let-7i, bta-miR-21, bta-miR-27, bta-miR-99b, bta-miR-146, bta-miR-147, bta-miR-155 and bta-miR-223). The most highly expressed miRNA in BAMs was miR-21, which has been shown to regulate the expression of antimicrobial peptides in Mycobacterium leprae-infected human monocytes. Furthermore, the predicted target genes of BAM-expressed miRNAs were found to be statistically enriched for roles in innate immunity. In addition to profiling the expression of known miRNAs, the RNA-seq data was also analysed to identify potentially novel bovine miRNAs. One putatively novel bovine miRNA was identified. To the best of our knowledge, this is the first RNA-seq study to profile miRNA expression in BAMs and provides an important reference dataset for investigating the regulatory roles miRNAs play in this important immune cell type.

Background Identification of gene expression profiles that differentiate experimental groups is critical for discovery and analysis of key molecular pathways and also for selection of robust diagnostic or prognostic biomarkers. While integration of differential expression statistics has been used to refine gene set enrichment analyses, such approaches are typically limited to single gene lists resulting from simple two-group comparisons or time-series analyses. In contrast, functional class scoring and machine learning approaches provide powerful alternative methods to leverage molecular measurements for pathway analyses, and to compare continuous and multi-level categorical factors. Results We introduce GOexpress, a software package for scoring and summarising the capacity of gene ontology features to simultaneously classify samples from multiple experimental groups. GOexpress integrates normalised gene expression data (e.g., from microarray and RNA-seq experiments) and phenotypic information of individual samples with gene ontology annotations to derive a ranking of genes and gene ontology terms using a supervised learning approach. The default random forest algorithm allows interactions between all experimental factors, and competitive scoring of expressed genes to evaluate their relative importance in classifying predefined groups of samples. Conclusions GOexpress enables rapid identification and visualisation of ontology-related gene panels that robustly classify groups of samples and supports both categorical (e.g., infection status, treatment) and continuous (e.g., time-series, drug concentrations) experimental factors. The use of standard Bioconductor extension packages and publicly available gene ontology annotations facilitates straightforward integration of GOexpress within existing computational biology pipelines. Electronic supplementary material The online version of this article (doi:10.1186/s12859-016-0971-3) contains supplementary material, which is available to authorized users.

Interleukin 8 (IL-8) is a major mediator of the innate immune response and polymorphisms in this gene are associated with susceptibility to inflammatory disease in humans. The aim of this study was to characterise the promoter region of the bovine IL8 gene towards understanding its regulation and the effect of promoter polymorphisms on gene expression levels. Twenty-nine polymorphic sites were identified across a 2.1kb upstream promoter region of the IL8 gene including two insertion/deletion polymorphisms. Sequence analysis and SNP genotyping identified two distinct promoter haplotypes (IL8-h1 and IL8-h2), which were present at significantly different frequencies in two divergently selected cattle breeds - Holstein-Friesian and Norwegian Red (IL8-h1 at 48% and 80% respectively). IL8-h1 was functionally less responsive in unstimulated mammary epithelial cells and in response to stimulation with LPS or bovine TNF. Serial deletion analysis and in silico transcription-factor binding site analysis indicated that allele specific binding of the transcriptional repressor Oct-1 may account for the reduced sensitivity of IL8-h1. Our finding of genetic variation in the bovine IL8 promoter that differentially regulates its expression has significant functional implications for IL8 expression in vitro and which may impact on susceptibility to bovine infectious disease and inflammation.

The growth hormone gene (GH1) and its polypeptide product (GH) have a crucial role in reproduction, embryogenesis and general development. A polymorphism present in the fifth exon of the bovine GH1 gene (GH1 p.Leu127Val) has been associated with GH release and milk production in cattle. The objective of the present study was to examine the genotype frequencies of the GH1 p.Leu127Val polymorphism in bovine blastocysts produced in vitro and in vivo to determine if allelic variation of the GH1 gene affects embryo development and survival. A heterozygous (p.Leu127/Val127) sire was used for in vitro fertilization of oocytes of unknown maternal genotype (n = 104) and known maternal genotype (n = 115). PCR amplification and genotyping of the GH1 gene from Day 8 blastocysts derived from these fertilized oocytes demonstrated that there was significant over-representation from the expected Mendelian ratio of GH1 p.Leu127/Leu127 homozygotes from oocytes of known maternal genotype (P = 0.006). Contrary to this, analysis of in vivo-produced bovine blastocysts of known parental GH1 genotype (n = 69) did not reveal an overrepresentation of GH1 p.Leu127/Leu127 homozygotes. These results suggest that developing in vitro-produced embryos are exposed to a selection process, probably due to a less favorable culture environment, that acts to increase the number of GH1 p.Leu127/Leu127 homozygotes, thereby giving rise to the observed transmission ratio distortion (TRD) of GH1 genotypes when compared to in vivo produced embryos.

Excess iodine intake by the pregnant dam reduces lamb serum antibody concentration, specifically immunoglobulin G (IgG). An experiment was conducted to investigate the mechanisms under pinning the reduced serum IgG concentration at 24 h postpartum in the progeny of iodine supplemented dams. Forty-five mature twin bearing ewes (n=15/treatment) were allocated to one of three dietary treatments as follows: basal diet (Control); basal diet plus 26.6 mg of iodine per ewe per day as calcium iodate (CaIO3); or potassium iodide (KI). Ewes were individually housed and fed from d 119 of gestation until parturition. All lambs received colostrum at 1, 10 and 18 h postpartum via stomach tube. At 1 h postpartum lambs from the control and an iodine supplemented treatment (n=10 per treatment from control and CaIO3) were euthanised before colostrum consumption and ileal segments isolated to determine the gene expression profile of a panel of genes identified as having a role in antibody transfer. Preceding euthanasia, lambs were blood sampled for determination of serum IgG, total thyroxine and free tri-iodothyronine concentrations. Progeny of CaIO3 supplemented dams had lower tri-iodothyronine concentrations (P

Mycobacterium bovis, the agent of bovine tuberculosis, causes an estimated $3 billion annual losses to global agriculture due, in part, to the limitations of current diagnostics. Development of next-generation diagnostics requires a greater understanding of the interaction between the pathogen and the bovine host. Therefore, to explore the early response of the alveolar macrophage to infection, we report the first application of RNA-sequencing to define, in exquisite detail, the transcriptomes of M. bovis-infected and non-infected alveolar macrophages from ten calves at 2, 6, 24 and 48 hours post-infection. Differentially expressed sense genes were detected at these time points that revealed enrichment of innate immune signalling functions and transcriptional suppression of host defence mechanisms (e.g., lysosome maturation). We also detected differentially expressed natural antisense transcripts, which may play a role in subverting innate immune mechanisms following infection. Furthermore, we report differential expression of novel bovine genes, some of which have immune-related functions based on orthology with human proteins. This is the first in-depth transcriptomics investigation of the alveolar macrophage response to the early stages of M. bovis infection and reveals complex patterns of gene expression and regulation that underlie the immunomodulatory mechanisms used by M. bovis to evade host defence mechanisms.

Johne’s disease, caused by infection with Mycobacterium avium subsp. paratuberculosis, (MAP), is a chronic intestinal disease of ruminants with serious economic consequences for cattle production in the United States and elsewhere. During infection, MAP bacilli are phagocytosed and subvert host macrophage processes, resulting in subclinical infections that can lead to immunopathology and dissemination of disease. Analysis of the host macrophage transcriptome during infection can therefore shed light on the molecular mechanisms and host-pathogen interplay associated with Johne’s disease. Here we describe results of an in vitro study of the bovine monocyte-derived macrophage (MDM) transcriptome response during MAP infection using RNA-seq. MDM were obtained from seven age- and sex-matched Holstein-Friesian cattle and were infected with MAP across a six-hour infection time course with non-infected controls. We observed 245 and 574 differentially expressed genes in MAP-infected versus non-infected control samples (adjusted P value ≤ 0.05) at 2 and 6 hours post-infection, respectively. Functional analyses of these differentially expressed genes, including biological pathway enrichment, highlighted potential functional roles for genes that have not been previously described in the host response to infection with MAP bacilli. In addition, differential expression of pro- and anti-inflammatory cytokine genes, such as those associated with the IL-10 signaling pathway, and other immune-related genes that encode proteins involved in the bovine macrophage response to MAP infection emphasize the balance between protective host immunity and bacilli survival and proliferation. Systematic comparisons of RNA-seq gene expression results with Affymetrix® microarray data generated from the same experimental samples also demonstrated that RNA-seq represents a superior technology for studying host transcriptional responses to intracellular infection.

To examine differences in cytokine profiles that may confer tolerance/susceptibility to bovine African trypanosomiasis, N'Dama (trypanotolerant, n = 8) and Boran (trypanosusceptible, n = 8) cattle were experimentally challenged with Trypanosoma congolense. Blood samples were collected over a 34-day period, and RNA was extracted from peripheral blood mononuclear cells. The expression levels of a panel of 14 cytokines were profiled over the time course of infection and between breeds. Messenger RNA (mRNA) transcript levels for the IL2, IL8, and IL1RN genes were significantly downregulated across the time course of infection in both breeds. There was an early increase in transcripts for genes encoding proinflammatory mediators ( IFNG, IL1A, TNF, and IL12) in N'Dama by 14 days postinfection (dpi) compared with preinfection levels that was not detected in the susceptible Boran breed. By the time of peak parasitemia, a type 2 helper T cells (TH2)-like cytokine environment was prevalent that was particularly evident in the Boran. Increases in transcripts for the IL6 (29 and 34 dpi) and IL10 (21, 25, and 29 dpi) genes were detected that were higher in the Boran compared with N'Dama. These findings highlight the implications for using murine models to study the bovine immune response to trypanosomiasis, where in some cases cytokine expression patterns differ. Overall, these data suggest that the trypanotolerant N'Dama are more capable of responding very early in infection with proinflammatory and TH1 type cytokines than the trypanosusceptible Boran and may explain why N'Dama control parasitemia more efficiently than Boran during the early stages of infection.

Microarray analysis of messenger RNA (mRNA) abundance was used to investigate the gene expression program of peripheral blood mononuclear cells (PBMC) from cattle infected with Mycobacterium bovis, the causative agent of bovine tuberculosis. An immunospecific bovine microarray platform (BOTL-4) with spot features representing 1336 genes was used for transcriptional profiling of PBMC from six M. bovis-infected cattle stimulated in vitro with bovine purified protein derivative of tuberculin (PPD-bovine). Cells were harvested at four time points (3 h, 6 h, 12 h and 24 h post-stimulation) and a split-plot design with pooled samples was used for the microarray experiment to compare gene expression between PPD-bovine stimulated PBMC and unstimulated controls for each time point. Statistical analyses of these data revealed 224 genes (approximately 17% of transcripts on the array) differentially expressed between stimulated and unstimulated PBMC across the 24 h time course (P0.05). Of the 224 genes, 87 genes were significantly upregulated and 137 genes were significantly downregulated in M. bovis-infected PBMC stimulated with PPD-bovine across the 24 h time course. However, perturbation of the PBMC transcriptome was most apparent at time points 3 h and 12 h post-stimulation, with 81 and 84 genes differentially expressed, respectively. In addition, a more stringent statistical threshold (P0.01) revealed 35 genes (approximately 3%) that were differentially expressed across the time course. Real-time quantitative reverse transcription PCR (qRT-PCR) of selected genes validated the microarray results and demonstrated a wide range of differentially expressed genes in PPD-bovine-, PPD-avian- and Concanavalin A (ConA) stimulated PBMC, including the interferon-gamma gene (IFNG), which was upregulated in PBMC stimulated with PPD-bovine (40-fold), PPD-avian (10-fold) and ConA (8-fold) after in vitro culture for 12 h. The pattern of expression of these genes in PPD-bovine stimulated PBMC provides the first description of an M. bovis-specific signature of infection that may provide insights into the molecular basis of the host response to infection. Although the present study was carried out with mixed PBMC cell populations, it will guide future studies to dissect immune cell-specific gene expression patterns in response to M. bovis infection.

Myostatin, encoded by the MSTN gene, is a member of the TGF-β superfamily that regulates skeletal muscle development. A MSTN SNP significantly associated with Thoroughbred horse racing phenotypes has recently been identified as well as significant reductions in Thoroughbred skeletal muscle gene expression for three transcripts 400-1500 base pairs downstream of the MSTN gene following a period of training. Together, these findings indicate that MSTN genotypes may influence MSTN gene expression. To investigate this, MSTN mRNA expression was measured in biopsies from the middle gluteal muscle from 60 untrained yearling Thoroughbreds (C/C, n = 15; C/T, n = 28; T/T, n = 17) using two independent real-time qRT-PCR assays. MSTN gene expression was also evaluated in a subset (N = 33) of these animals using samples collected after a ten-month period of training. A significant association was observed between genotype and mRNA abundance for the untrained horses (assay I, P = 0.0237; assay II, P = 0.003559), with the C/C cohort having the highest MSTN mRNA levels, the T/T group the lowest levels and the C/T group intermediate levels. Following training, there was a significant decrease in MSTN mRNA (-3.35-fold; P = 6.9 × 10(-7) ), which was most apparent for the C/C cohort (-5.88-fold, P = 0.001). These data demonstrate the tight relationship between phenotype, genotype and gene expression at the MSTN gene in Thoroughbred racehorses.

Profitable milk production in dairy cows requires good reproductive performance. Calving interval is a trait used to measure reproductive efficiency. Herein we used a novel lactating Holstein cow model of fertility that displayed genetic and phenotypic divergence in calving interval, a trait used to define reproductive performance using a national breeding index in Ireland. Cows had similar genetic merit for milk production traits, but either very good genetic merit for fertility (Fert+; n = 7) or very poor genetic merit for fertility (Fert–; n = 6). We tested the hypothesis that Fert+ cows would have a corresponding detectable difference in endometrial gene expression compared with the Fert– cows. To do this, we sequenced the transcriptome of endometrial biopsies collected on Day 7 of the oestrous cycle (non-pregnant). This is an important stage for uterine remodelling and initiation of histotroph secretion. Significant differential expression (false discovery rate-adjusted P

peer-reviewed Mycobacterium bovis, the causative agent of bovine tuberculosis, a major problem for global agriculture, spreads via an airborne route and is taken up by alveolar macrophages (AM) in the lung. Here, we describe the first next-generation sequencing (RNA-seq) approach to temporally profile miRNA expression in primary bovine AMs post-infection with M. bovis. One, six, and forty miRNAs were identified as significantly differentially expressed at 2, 24 and 48 h post-infection, respectively. The differential expression of three miRNAs (bta-miR-142-5p, bta-miR-146a, and bta-miR-423-3p) was confirmed by RT-qPCR. Pathway analysis of the predicted mRNA targets of differentially expressed miRNAs suggests that these miRNAs preferentially target several pathways that are functionally relevant for mycobacterial pathogenesis, including endocytosis and lysosome trafficking, IL-1 signalling and the TGF-β pathway. Over-expression studies using a bovine macrophage cell-line (Bomac) reveal the targeting of two key genes in the innate immune response to M. bovis, IL-1 receptor-associated kinase 1 (IRAK1) and TGF-β receptor 2 (TGFBR2), by miR-146. Taken together, our study suggests that miRNAs play a key role in tuning the complex interplay between M. bovis survival strategies and the host immune response.

Bovine tuberculosis (BTB), caused by infection with Mycobacterium bovis, is a major endemic disease affecting cattle populations worldwide, despite the implementation of stringent surveillance and control programs in many countries. The development of high-throughput functional genomics technologies, including gene expression microarrays and RNA sequencing (RNA-seq), has enabled detailed analysis of the host transcriptome to M. bovis infection, particularly at the macrophage and peripheral blood level. In the present study, we have analyzed the peripheral blood leukocyte (PBL) transcriptome of eight natural M. bovis-infected and eight age- and sex-matched non-infected control Holstein-Friesian animals using RNA-seq. In addition, we compared gene expression profiles generated using RNA-seq with those previously generated using the high-density Affymetrix® GeneChip® Bovine Genome Array platform from the same PBL-extracted RNA. A total of 3,250 differentially expressed (DE) annotated genes were detected in the M. bovis-infected samples relative to the controls (adjusted P-value ≤ 0.05), with the number of genes displaying decreased relative expression (1,671) exceeding those with increased relative expression (1,579). Ingenuity® Systems Pathway Analysis (IPA) of all DE genes revealed enrichment for genes with immune function. Notably, transcriptional suppression was observed among several of the top ranking canonical pathways including Leukocyte Extravasation Signaling. Comparative platform analysis demonstrated that RNA-seq detected a larger number of annotated DE genes (3,250) relative to the microarray (1,398), of which 917 genes were common to both technologies and displayed the same direction of expression. Finally, we show that RNA-seq had an increased dynamic range compared to the microarray for estimating differential gene expression.

Background Mycobacterium bovis, the causative agent of bovine tuberculosis, is an intracellular pathogen that can persist inside host macrophages during infection via a diverse range of mechanisms that subvert the host immune response. In the current study, we have analysed and compared the transcriptomes of M. bovis-infected monocyte-derived macrophages (MDM) purified from six Holstein-Friesian females with the transcriptomes of non-infected control MDM from the same animals over a 24 h period using strand-specific RNA sequencing (RNA-seq). In addition, we compare gene expression profiles generated using RNA-seq with those previously generated by us using the high-density Affymetrix® GeneChip® Bovine Genome Array platform from the same MDM-extracted RNA. Results A mean of 7.2 million reads from each MDM sample mapped uniquely and unambiguously to single Bos taurus reference genome locations. Analysis of these mapped reads showed 2,584 genes (1,392 upregulated; 1,192 downregulated) and 757 putative natural antisense transcripts (558 upregulated; 119 downregulated) that were differentially expressed based on sense and antisense strand data, respectively (adjusted P-value ≤ 0.05). Of the differentially expressed genes, 694 were common to both the sense and antisense data sets, with the direction of expression (i.e. up- or downregulation) positively correlated for 693 genes and negatively correlated for the remaining gene. Gene ontology analysis of the differentially expressed genes revealed an enrichment of immune, apoptotic and cell signalling genes. Notably, the number of differentially expressed genes identified from RNA-seq sense strand analysis was greater than the number of differentially expressed genes detected from microarray analysis (2,584 genes versus 2,015 genes). Furthermore, our data reveal a greater dynamic range in the detection and quantification of gene transcripts for RNA-seq compared to microarray technology. Conclusions This study highlights the value of RNA-seq in identifying novel immunomodulatory mechanisms that underlie host-mycobacterial pathogen interactions during infection, including possible complex post-transcriptional regulation of host gene expression involving antisense RNA.

The potential for dietary supplementation with n-3 polyunsaturated fatty acids (n-3 PUFA) to improve reproductive efficiency in cattle has received much interest. The mechanisms by which n-3 PUFA may affect physiological and biochemical processes in key reproductive tissues are likely to be mediated by significant alterations in gene expression. The objective of this study was to examine the effects of dietary n-3 PUFA supplementation on global uterine endometrial gene expression in cattle. Beef heifers were supplemented with a rumen protected source of either a saturated fatty acid (CON; palmitic acid) or high n-3 PUFA (n-3 PUFA; 275 g) diet per animal per day for 45 days and global gene expression was determined in uterine endometrial tissue using an Affymetrix oligonucleotide bovine array. A total of 1,807 (946 up- and 861 downregulated) genes were differentially expressed following n-3 PUFA supplementation. Dietary n-3 PUFA altered numerous cellular processes potentially important in the control of reproduction in cattle. These included prostaglandin biosynthesis, steroidogenesis and transcriptional regulation, while effects on genes involved in maternal immune response and tissue remodeling were also observed. This study provides new insights into the effects of n-3 PUFA supplementation on the regulation of gene expression in the bovine uterus.

Mycobacterium avium subspecies paratuberculosis is the causative agent of Johne’s disease, an intestinal disease of ruminants with major economic consequences. Infectious bacilli are phagocytosed by host macrophages upon exposure where they persist, resulting in lengthy subclinical phases of infection that can lead to immunopathology and disease dissemination. Consequently, analysis of the macrophage transcriptome in response to M. avium subsp. paratuberculosis infection can provide valuable insights into the molecular mechanisms that underlie Johne’s disease. Here, we investigate pan-genomic gene expression in bovine monocyte-derived macrophages (MDM) purified from seven age-matched females, in response to in vitro infection with M. avium subsp. paratuberculosis (multiplicity of infection 2:1) at intervals of 2 hours, 6 hours and 24 hours post-infection (hpi). Differentially expressed genes were identified by comparing the transcriptomes of the infected MDM to the non-infected control MDM at each time point (adjusted P-value threshold ≤ 0.10). 1050 differentially expressed unique genes were identified 2 hpi, with 974 and 78 differentially expressed unique genes detected 6 and 24 hpi, respectively. Furthermore, in the infected MDM the number of upregulated genes exceeded the number of downregulated genes at each time point, with the fold-change in expression for the upregulated genes markedly higher than that for the downregulated genes. Inspection and systems biology analysis of the differentially expressed genes revealed an enrichment of genes involved in the inflammatory response, cell signalling pathways and apoptosis. The transcriptional changes associated with cellular signalling and the inflammatory response may reflect different immuno-modulatory mechanisms that underlie host-pathogen interactions during infection.

peer-reviewed Background Mycobacterium bovis, the causative agent of bovine tuberculosis, is a major cause of mortality in global cattle populations. Macrophages are among the first cell types to encounter M. bovis following exposure and the response elicited by these cells is pivotal in determining the outcome of infection. Here, a functional genomics approach was undertaken to investigate global gene expression profiles in bovine monocyte-derived macrophages (MDM) purified from seven age-matched non-related females, in response to in vitro challenge with M. bovis (multiplicity of infection 2:1). Total cellular RNA was extracted from non-challenged control and M. bovis-challenged MDM for all animals at intervals of 2 hours, 6 hours and 24 hours post-challenge and prepared for global gene expression analysis using the Affymetrix® GeneChip® Bovine Genome Array. Results Comparison of M. bovis-challenged MDM gene expression profiles with those from the non-challenged MDM controls at each time point identified 3,064 differentially expressed genes 2 hours post-challenge, with 4,451 and 5,267 differentially expressed genes detected at the 6 hour and 24 hour time points, respectively (adjusted P-value threshold ≤0.05). Notably, the number of downregulated genes exceeded the number of upregulated genes in the M. bovis-challenged MDM across all time points; however, the fold-change in expression for the upregulated genes was markedly higher than that for the downregulated genes. Systems analysis revealed enrichment for genes involved in: (1) the inflammatory response; (2) cell signalling pathways, including Toll-like receptors and intracellular pathogen recognition receptors; and (3) apoptosis. Conclusions The increased number of downregulated genes is consistent with previous studies showing that M. bovis infection is associated with the repression of host gene expression. The results also support roles for MyD88-independent signalling and intracellular PRRs in mediating the host response to M. bovis. Science Foundation Ireland (www.sfi.ie) Investigator grants (Nos: SFI/01/F.1/B028 and SFI/08/IN.1/B2038); Department of Agriculture, Fisheries and Food (www.agriculture.ie) Research Stimulus Grant (No: RSF 06 405); European Union Framework 7 (http://cordis.europa.eu/fp7) Project Grant (No: KBBE-211602-MACROSYS); Irish Research Council for Science, Engineering and Technology (IRCSET) funded Bioinformatics and Systems Biology PhD Programme (http://bioinfo-casl.ucd.ie/PhD).

RT-qPCR can be used to accurately determine expression levels of genes following RNA extraction from tissue samples. If blood is the source of total RNA, it is often desirable to process the samples immediately following collection because delays in processing for RNA extraction may influence mRNA expression estimates obtained from RT-qPCR analyses. However, this may not be feasible if the site of blood collection is distant from the processing laboratory. In the present study, the effects of delays in the processing of blood samples on mRNA expression data was investigated using a panel of 23 functionally diverse genes from five different gene ontology (GO) categories in peripheral blood sampled from ten age-matched healthy cattle. Venous blood was collected in Tempus™ Blood RNA tubes, which contain reagents that lyse blood cells immediately and stabilise the RNA signature (T0). Blood was also collected in conventional lithium heparin collection tubes, and stored at ambient temperature for T4, T6 and T8 h, prior to total RNA extraction. The mRNA expression profiles of these 23 genes were determined by RT-qPCR and compared across the time course. Thirteen genes showed significant up- or down-fold changes in mRNA expression over the 8 h time course. Among the GO categories, genes in the Immune response category showed the most differential expression. These results also demonstrated that the changes in mRNA expression for the IFNG gene, which encodes the cytokine IFN-γ, did not correspond to IFN-γ protein levels estimated using ELISA.

The involvement of Toll-like receptors (TLRs) and other immune signalling genes during challenge of bovine macrophages with bacterial products derived from disease-causing bacteria in cattle was investigated. An in vitro cell culture model of bovine monocyte derived macrophages (MDM) was established and these cells were exposed to purified protein derivative (PPD-b) derived from Mycobacterium bovis and to lipopolysachharide (LPS) derived from Escherichia coli. Following 24 h incubation, total RNA was extracted and expression of immune related genes was determined by real time quantitative reverse transcription PCR (qRT-PCR). Expression of a selection of genes spanning the TLR-2 and TLR-4 pathways, from the initial activation of the receptors to the production of pro-inflammatory cytokines and chemokines was determined. Results from repeat experiments using MDM from seven different age-matched dairy cattle showed that PPD-b treatment caused significant up-regulation of the TLR2 and TLR4 genes and the expression profile of TLR adaptor molecules suggested that this signalling is MYD88-dependent. Conversely, LPS caused significant up-regulation of TLR4 via a MYD88-independent signalling pathway. Significant up-regulation of genes involved with NF-κB signalling was also detected in PPD-b- and LPS-treated samples accompanied by the expression of pro-inflammatory cytokine (TNF, IL1B, IL6) and chemokine genes (IL8, CCL5, CCL3). Overall, LPS challenge resulted in a more marked up-regulation of immune-related genes. Furthermore, the magnitude fold-change difference in gene expression suggests, at least in part, that bovine macrophages produce IFN-γ as a result of LPS challenge.

Previous studies show that DNA sequence variation within the mammalian DLK1-DIO3 imprinted domain influences production traits in domestic livestock, most notably the ovine callipyge phenotype. We assessed genotype-phenotype associations between 7 single nucleotide polymorphisms (SNPs) within the orthologous bovine DLK1-DIO3 domain and performance traits in 848 progeny-tested Holstein-Friesian dairy sires. One SNP (MEG3_01) located proximal to the maternally expressed 3 (MEG3/Gtl2) gene was associated with milk yield, subcutaneous fat levels, and progeny carcass conformation (P < 0.01) and also tended to be associated with milk fat and protein yield (P < 0.10). A single SNP (CLPG_01) within the putative CLPG1 locus was associated with progeny carcass fat (P < 0.05), whereas a single SNP (PEG11_01) located proximal to the paternally expressed 11 (PEG11/Rtl) gene was associated with progeny carcass weight (P < 0.05). The MEG3_01 SNP together with an additional 2 SNPs (MEG8_01 and MEG8_02) located proximal to the putative maternally expressed 8 (MEG8/Rian) ortholog were associated (P < 0.05) with perinatal mortality. Finally, one SNP (MEG3_03) was associated (P < 0.05) with gestation length, whereas both the CLPG_01 and MEG8_01 SNPs also tended to be associated with calving interval (P < 0.10). Linkage disequilibrium analysis suggests that some phenotypic associations observed at these loci are independent. To our knowledge, this is one of the first studies demonstrating associations between the bovine DLK1-DIO3 domain and milk, carcass, fertility and, health traits in cattle. This imprinted domain may serve as a potential target for future genetic selection strategies.

Background Selection for exercise-adapted phenotypes in the Thoroughbred racehorse has provided a valuable model system to understand molecular responses to exercise in skeletal muscle. Exercise stimulates immediate early molecular responses as well as delayed responses during recovery, resulting in a return to homeostasis and enabling long term adaptation. Global mRNA expression during the immediate-response period has not previously been reported in skeletal muscle following exercise in any species. Also, global gene expression changes in equine skeletal muscle following exercise have not been reported. Therefore, to identify novel genes and key regulatory pathways responsible for exercise adaptation we have used equine-specific cDNA microarrays to examine global mRNA expression in skeletal muscle from a cohort of Thoroughbred horses (n = 8) at three time points (before exercise, immediately post-exercise, and four hours post-exercise) following a single bout of treadmill exercise. Results Skeletal muscle biopsies were taken from the gluteus medius before (T0), immediately after (T1) and four hours after (T2) exercise. Statistically significant differences in mRNA abundance between time points (T0 vs T1 and T0 vs T2) were determined using the empirical Bayes moderated t-test in the Bioconductor package Linear Models for Microarray Data (LIMMA) and the expression of a select panel of genes was validated using real time quantitative reverse transcription PCR (qRT-PCR). While only two genes had increased expression at T1 (P < 0.05), by T2 932 genes had increased (P < 0.05) and 562 genes had decreased expression (P < 0.05). Functional analysis of genes differentially expressed during the recovery phase (T2) revealed an over-representation of genes localized to the actin cytoskeleton and with functions in the MAPK signalling, focal adhesion, insulin signalling, mTOR signaling, p53 signaling and Type II diabetes mellitus pathways. At T1, using a less stringent statistical approach, we observed an over-representation of genes involved in the stress response, metabolism and intracellular signaling. These findings suggest that protein synthesis, mechanosensation and muscle remodeling contribute to skeletal muscle adaptation towards improved integrity and hypertrophy. Conclusions This is the first study to characterize global mRNA expression profiles in equine skeletal muscle using an equine-specific microarray platform. Here we reveal novel genes and mechanisms that are temporally expressed following exercise providing new knowledge about the early and late molecular responses to exercise in the equine skeletal muscle transcriptome.

Intense selection for elite racing performance in the Thoroughbred horse ( Equus caballus) has resulted in a number of adaptive physiological phenotypes relevant to exercise; however, the underlying molecular mechanisms responsible for these characteristics are not well understood. Adaptive changes in mRNA expression in equine skeletal muscle were investigated by real-time qRT-PCR for a panel of candidate exercise-response genes following a standardized incremental-step treadmill exercise test in eight untrained Thoroughbred horses. Biopsy samples were obtained from the gluteus medius before, immediately after, and 4 h after exercise. Significant ( P < 0.05) differences in gene expression were detected for six genes ( CKM, COX4I1, COX4I2, PDK4, PPARGC1A, and SLC2A4) 4 h after exercise. Investigation of relationships between mRNA and velocity at maximum heart rate (VHRmax) and peak postexercise plasma lactate concentration ([La]T1) revealed significant ( P < 0.05) associations with postexercise COX4I1 and PPARCG1A expression and between [La]T1and basal COX4I1 expression. Gene expression changes were investigated in a second cohort of horses after a 10 mo period of training. In resting samples, COX4I1 gene expression had significantly increased following training, and, after exercise, significant differences were identified for COX4I2, PDK4, and PPARGC1A. Significant relationships with VHRmaxand [La]T1were detected for PPARGC1A and COX4I1. These data highlight the roles of genes responsible for the regulation of oxygen-dependent metabolism, glucose metabolism, and fatty acid utilization in equine skeletal muscle adaptation to exercise.

Summary Bovine tuberculosis (BTB), caused by Mycobacterium bovis, continues to pose a threat to livestock worldwide and, as a zoonotic infection, also has serious implications for human health. The implementation of comprehensive surveillance programmes to detect BTB has been successful in reducing the incidence of infection in many countries, yet BTB has remained recalcitrant to eradication in several EU states, particularly in Ireland and the UK. There are well-recognized limitations in the use of the current diagnostics to detect all infected animals and this has led to renewed efforts to uncover novel diagnostic biomarkers that may serve to enhance the performance of the tests. Studies of single immunological parameters have so far been unable to unlock the complexities of the immune response to mycobacterial infection. However, the development of high-throughput methods including pan-genomic gene expression technologies such as DNA microarrays has facilitated the simultaneous identification and analysis of thousands of genes and their interactions during the immune response. In addition, the application of these new genomic technologies to BTB has identified pathogen-associated immune response signatures of host infection. The objective of these investigations is to understand the changing profile of immune responses throughout the course of infection and to identify biomarkers for sensitive diagnosis, particularly during the early stages of infection. Transcriptional profiling via microarray and more recently via next-generation sequencing technologies may lead to the development of specific and sensitive diagnostics for M. bovis infection and will enhance the prospect of eradication of tuberculosis from cattle populations.

Background Bovine tuberculosis (BTB) caused by Mycobacterium bovis continues to cause substantial losses to global agriculture and has significant repercussions for human health. The advent of high throughput genomics has facilitated large scale gene expression analyses that present a novel opportunity for revealing the molecular mechanisms underlying mycobacterial infection. Using this approach, we have previously shown that innate immune genes in peripheral blood mononuclear cells (PBMC) from BTB-infected animals are repressed in vivo in the absence of exogenous antigen stimulation. In the present study, we hypothesized that the PBMC from BTB-infected cattle would display a distinct gene expression program resulting from exposure to M. bovis. A functional genomics approach was used to examine the immune response of BTB-infected (n = 6) and healthy control (n = 6) cattle to stimulation with bovine tuberculin (purified protein derivative – PPD-b) in vitro. PBMC were harvested before, and at 3 h and 12 h post in vitro stimulation with bovine tuberculin. Gene expression changes were catalogued within each group using a reference hybridization design and a targeted immunospecific cDNA microarray platform (BOTL-5) with 4,800 spot features representing 1,391 genes. Results 250 gene spot features were significantly differentially expressed in BTB-infected animals at 3 h post-stimulation contrasting with only 88 gene spot features in the non-infected control animals (P ≤ 0.05). At 12 h post-stimulation, 56 and 80 gene spot features were differentially expressed in both groups respectively. The results provided evidence of a proinflammatory gene expression profile in PBMC from BTB-infected animals in response to antigen stimulation. Furthermore, a common panel of eighteen genes, including transcription factors were significantly expressed in opposite directions in both groups. Real-time quantitative reverse transcription PCR (qRT-PCR) demonstrated that many innate immune genes, including components of the TLR pathway and cytokines were differentially expressed in BTB-infected (n = 8) versus control animals (n = 8) after stimulation with bovine tuberculin. Conclusion The PBMC from BTB-infected animals exhibit different transcriptional profiles compared with PBMC from healthy control animals in response to M. bovis antigen stimulation, providing evidence of a novel gene expression program due to M. bovis exposure.

Background African animal trypanosomiasis (AAT) caused by tsetse fly-transmitted protozoa of the genus Trypanosoma is a major constraint on livestock and agricultural production in Africa and is among the top ten global cattle diseases impacting on the poor. Here we show that a functional genomics approach can be used to identify temporal changes in host peripheral blood mononuclear cell (PBMC) gene expression due to disease progression. We also show that major gene expression differences exist between cattle from trypanotolerant and trypanosusceptible breeds. Using bovine long oligonucleotide microarrays and real time quantitative reverse transcription PCR (qRT-PCR) validation we analysed PBMC gene expression in naïve trypanotolerant and trypanosusceptible cattle experimentally challenged with Trypanosoma congolense across a 34-day infection time course. Results Trypanotolerant N'Dama cattle displayed a rapid and distinct transcriptional response to infection, with a ten-fold higher number of genes differentially expressed at day 14 post-infection compared to trypanosusceptible Boran cattle. These analyses identified coordinated temporal gene expression changes for both breeds in response to trypanosome infection. In addition, a panel of genes were identified that showed pronounced differences in gene expression between the two breeds, which may underlie the phenomena of trypanotolerance and trypanosusceptibility. Gene ontology (GO) analysis demonstrate that the products of these genes may contribute to increased mitochondrial mRNA translational efficiency, a more pronounced B cell response, an elevated activation status and a heightened response to stress in trypanotolerant cattle. Conclusion This study has revealed an extensive and diverse range of cellular processes that are altered temporally in response to trypanosome infection in African cattle. Results indicate that the trypanotolerant N'Dama cattle respond more rapidly and with a greater magnitude to infection compared to the trypanosusceptible Boran cattle. Specifically, a subset of the genes analyzed by real time qRT-PCR, which display significant breed differences, could collectively contribute to the trypanotolerance trait in N'Dama.

peer-reviewed Background: Recovering high quality intact RNA from post-mortem tissue is of major concern for gene expression studies in animals and humans. Since the availability of post-mortem tissue is often associated with substantial delay, it is important that we understand the temporal variation in the stability of total RNA and of individual gene transcripts so as to be able to appropriately interpret the data generated from such studies. Hence, the objective of this experiment was to qualitatively and quantitatively assess the integrity of total and messenger RNA extracted from bovine skeletal muscle, subcutaneous adipose tissue and liver stored at 4°C at a range of time points up to 22 days post-mortem. These conditions were designed to mimic the environment prevailing during the transport of beef from the abattoir to retail outlets. Results: The 28S and 18S rRNA molecules of total RNA were intact for up to 24 h post-mortem in liver and adipose tissues and up to 8 days post-mortem in skeletal muscle. The mRNA of housekeeping genes (GAPDH and ACTB) and two diet-related genes (RBP5 and SCD) were detectable up to 22 days post-mortem in skeletal muscle. While the mRNA stability of the two housekeeping genes was different in skeletal muscle and liver, they were similar to each other in adipose tissue. After 22 days post-mortem, the relative abundance of RBP5 gene was increased in skeletal muscle and in adipose tissue and decreased in liver. During this period, the relative abundance of SCD gene also increased in skeletal muscle whereas it decreased in both adipose tissue and liver. Conclusion: Stability of RNA in three tissues (skeletal muscle, subcutaneous adipose tissue and liver) subjected to long-term post-mortem storage at refrigeration temperature indicated that skeletal muscle can be a suitable tissue for recovering biologically useful RNA for gene expression studies even if the tissue is subjected to post-mortem storage for weeks, whereas adipose tissue and liver should be processed within 24 hours post-mortem. This research was part funded by the Irish Government under the National Development Plan 2000-2006. Support to BB under the Walsh Fellowship Programme of Teagasc is acknowledged.

PUBLISHED This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited., Bovine tuberculosis is an enduring disease of cattle that has significant repercussions for human health. The advent of high-throughput functional genomics technologies has facilitated large-scale analyses of the immune response to this disease that may ultimately lead to novel diagnostics and therapeutic targets. Analysis of mRNA abundance in peripheral blood mononuclear cells (PBMC) from six Mycobacterium bovis infected cattle and six non-infected controls was performed. A targeted immunospecific bovine cDNA microarray with duplicated spot features representing 1,391 genes was used to test the hypothesis that a distinct gene expression profile may exist in M. bovis infected animals in vivo. Results: In total, 378 gene features were differentially expressed at the P ? 0.05 level in bovine tuberculosis (BTB)-infected and control animals, of which 244 were expressed at lower levels (65%) in the infected group. Lower relative expression of key innate immune genes, including the Toll-like receptor 2 (TLR2) and TLR4 genes, lack of differential expression of indicator adaptive immune gene transcripts (IFNG, IL2, IL4), and lower BOLA major histocompatibility complex ? class I (BOLA) and class II (BOLA-DRA) gene expression was consistent with innate immune gene repression in the BTB-infected animals. Supervised hierarchical cluster analysis and class prediction validation identified a panel of 15 genes predictive of disease status and selected gene transcripts were validated (n = 8 per group) by real time quantitative reverse transcription PCR. Conclusion: These results suggest that large-scale expression profiling can identify gene signatures of disease in peripheral blood that can be used to classify animals on the basis of in vivo infection, in the absence of exogenous antigenic stimulation., We would like to thank the staff at the UCD Lyons Research Farm for assistance with cattle blood sampling and Grace O'Gorman, Stephen Park and Emmeline Hill for technical assistance and advice. We also thank Frances Quigley and staff at the Mycobacteriology Laboratory, CVRL, Backweston for their assistance in the bacterial culture of M. bovis from infected cattle and Prof. Paul Coussens of the MSU Centre for Animal Functional Genomics for advice and provision of BOTL-5 microarrays. This work was supported by an Investigator Grant from Science Foundation Ireland (Grant No: SFI/01/F.1/B028)

Background Mycobacterium bovis is the causative agent of bovine tuberculosis (BTB), a pathological infection with significant economic impact. Recent studies have highlighted the role of functional genomics to better understand the molecular mechanisms governing the host immune response to M. bovis infection. Furthermore, these studies may enable the identification of novel transcriptional markers of BTB that can augment current diagnostic tests and surveillance programmes. In the present study, we have analysed the transcriptome of peripheral blood leukocytes (PBL) from eight M. bovis-infected and eight control non-infected age-matched and sex-matched Holstein-Friesian cattle using the Affymetrix® GeneChip® Bovine Genome Array with 24,072 gene probe sets representing more than 23,000 gene transcripts. Results Control and infected animals had similar mean white blood cell counts. However, the mean number of lymphocytes was significantly increased in the infected group relative to the control group (P = 0.001), while the mean number of monocytes was significantly decreased in the BTB group (P = 0.002). Hierarchical clustering analysis using gene expression data from all 5,388 detectable mRNA transcripts unambiguously partitioned the animals according to their disease status. In total, 2,960 gene transcripts were differentially expressed (DE) between the infected and control animal groups (adjusted P-value threshold ≤ 0.05); with the number of gene transcripts showing decreased relative expression (1,563) exceeding those displaying increased relative expression (1,397). Systems analysis using the Ingenuity® Systems Pathway Analysis (IPA) Knowledge Base revealed an over-representation of DE genes involved in the immune response functional category. More specifically, 64.5% of genes in the affects immune response subcategory displayed decreased relative expression levels in the infected animals compared to the control group. Conclusions This study demonstrates that genome-wide transcriptional profiling of PBL can distinguish active M. bovis-infected animals from control non-infected animals. Furthermore, the results obtained support previous investigations demonstrating that mycobacterial infection is associated with host transcriptional suppression. These data support the use of transcriptomic technologies to enable the identification of robust, reliable transcriptional markers of active M. bovis infection.

Background Digital gene expression profiling was used to characterize the assembly of genes expressed in equine skeletal muscle and to identify the subset of genes that were differentially expressed following a ten-month period of exercise training. The study cohort comprised seven Thoroughbred racehorses from a single training yard. Skeletal muscle biopsies were collected at rest from the gluteus medius at two time points: T1 - untrained, (9 ± 0.5 months old) and T2 - trained (20 ± 0.7 months old). Results The most abundant mRNA transcripts in the muscle transcriptome were those involved in muscle contraction, aerobic respiration and mitochondrial function. A previously unreported over-representation of genes related to RNA processing, the stress response and proteolysis was observed. Following training 92 tags were differentially expressed of which 74 were annotated. Sixteen genes showed increased expression, including the mitochondrial genes ACADVL, MRPS21 and SLC25A29 encoded by the nuclear genome. Among the 58 genes with decreased expression, MSTN, a negative regulator of muscle growth, had the greatest decrease. Functional analysis of all expressed genes using FatiScan revealed an asymmetric distribution of 482 Gene Ontology (GO) groups and 18 KEGG pathways. Functional groups displaying highly significant (P < 0.0001) increased expression included mitochondrion, oxidative phosphorylation and fatty acid metabolism while functional groups with decreased expression were mainly associated with structural genes and included the sarcoplasm, laminin complex and cytoskeleton. Conclusion Exercise training in Thoroughbred racehorses results in coordinate changes in the gene expression of functional groups of genes related to metabolism, oxidative phosphorylation and muscle structure.