Ruminal Transcriptomic Analysis of Grass-Fed and Grain-Fed Angus Beef Cattle

Transcriptomic profiles of the ruminal wall in Italian Mediterranean dairy buffaloes fed green forage

BACKGROUND

Green feed diet in ruminants exerts a beneficial effect on rumen metabolism and enhances the content of milk nutraceutical quality. At present, a comprehensive analysis focused on the identification of genes, and therefore, biological processes modulated by the green feed in buffalo rumen has never been reported. We performed RNA-sequencing in the rumen of buffaloes fed a total mixed ration (TMR) + the inclusion of 30% of ryegrass green feed (treated) or TMR (control), and identified differentially expressed genes (DEGs) using EdgeR and NOISeq tools.

RESULTS

We found 155 DEGs using EdgeR (p-values < 0.05) and 61 DEGs using NOISeq (prob ≥0.8), 30 of which are shared. The rt-qPCR validation suggested a higher reliability of EdgeR results as compared with NOISeq data, in our biological context. Gene Ontology analysis of DEGs identified using EdgeR revealed that green feed modulates biological processes relevant for the rumen physiology and, then, health and well-being of buffaloes, such as lipid metabolism, response to the oxidative stress, immune response, and muscle structure and function. Accordingly, we found: (i) up-regulation of HSD17B13, LOC102410803 (or PSAT1) and HYKK, and down-regulation of CDO1, SELENBP1 and PEMT, encoding factors involved in energy, lipid and amino acid metabolism; (ii) enhanced expression of SIM2 and TRIM14, whose products are implicated in the immune response and defense against infections, and reduced expression of LOC112585166 (or SAAL1), ROR2, SMOC2, and S100A11, encoding pro-inflammatory factors; (iii) up-regulation of NUDT18, DNAJA4 and HSF4, whose products counteract stressful conditions, and down-regulation of LOC102396388 (or UGT1A9) and LOC102413340 (or MRP4/ABCC4), encoding detoxifying factors; (iv) increased expression of KCNK10, CACNG4, and ATP2B4, encoding proteins modulating Ca²⁺ homeostasis, and reduced expression of the cytoskeleton-related MYH11 and DES.

CONCLUSION

Although statistically unpowered, this study suggests that green feed modulates the expression of genes involved in biological processes relevant for rumen functionality and physiology, and thus, for welfare and quality production in Italian Mediterranean dairy buffaloes. These findings, that need to be further confirmed through the validation of additional DEGs, allow to speculate a role of green feed in the production of nutraceutical molecules, whose levels might be enhanced also in milk.

Impact of Cattle Feeding Strategy on the Beef Metabolome

The present study explored changes in the meat metabolome of animals subjected to different finishing systems and growth rates. Thirty-six Angus × Nellore crossbred steers were used in a completely randomized design with four treatments: (1) feedlot system with high average daily gain (ADG; FH); (2) feedlot system with low ADG (FL); (3) pasture system with high ADG (PH); and (4) pasture system with low ADG (PL). After harvest and chilling, Longissimus thoracis (LT) muscle samples were taken for metabolite profile analysis using nuclear magnetic resonance. Spectrum was analyzed using chenomx software, and multi- and mega-variate data analyses were performed. The PLS-DA showed clear separation between FH and PL groups and overlap among treatments with different finishing systems but similar for matching ADG (FL and PH) treatments. Using a VIP cut-off of around 1.0, ATP and fumarate were shown to be greater in meat from PL cattle, while succinate, leucine, AMP, glutamate, carnosine, inosine, methionine, G1P, and choline were greater in meat from FH. Comparing FL and PH treatments, glutamine, carnosine, urea, NAD+, malonate, lactate, isoleucine, and alanine were greater in the meat of PH cattle, while G6P and betaine were elevated in that of FL cattle. Relevant pathways were also identified by differences in growth rate (FH versus PL) and finishing system were also noted. Growth rate caused a clear difference in meat metabolism that was highlighted by energy metabolism and associated pathways, while the feeding system tended to alter protein and lipid metabolism.

Plasma Amino Acid Appearance and Status of Appetite Following a Single Meal of Red Meat or a Plant-Based Meat Analog: A Randomized Crossover Clinical Trial

Background

Red meat is a nutrient-dense food and a dietary staple. A new generation of plant-based meat analogs (PBMAs) have been designed to mimic the experience of eating meat, but there is limited evidence about their digestive efficacy and nutritional quality.

Objectives

We compared the postprandial digestive response of a single meal containing meat commercially raised in New Zealand, including lamb, on-farm pasture-raised beef (Pasture), or grain-finished beef (Grain) with a PBMA (Beyond Burger; Beyond Meat) sold through consumer retail. The primary outcome was the appearance of amino acids in plasma. Secondary outcomes included glucose and insulin, appetite assessment, and anthropometry.

Methods

Thirty healthy men (20-34 y) participated in a double-blinded randomized crossover trial. Each consumed 1 of the 4 test meals on 4 occasions separated by a washout period of at least 1 wk, following an overnight fast. The meal was a burrito-style wrap containing meat or PBMAs, vegetables, salsa, and seasonings in a flour tortilla. The amount of Pasture, Grain, Lamb, or BB was 220 g raw (∼160 g cooked). Venous blood samples were collected over 4 h. Appetite and hunger status was scored with visual analog scales.

Results

Pre-meal amino acid concentrations in plasma did not differ by group (P > 0.9), although several nonessential amino acids differed strongly according to participant BMI. Postprandial amino acids peaked at 2-3 h in all groups. The BB meal produced significantly lower plasma concentrations of total, essential, branched-chain, and non-proteogenic amino acids than the Lamb, Pasture, or Grain meals, based on AUC. There were no significant differences between meal groups in scores for hunger, fullness, or cravings.

Conclusions

Red meat meals exhibited greater bioavailability of amino acids compared with the PBMA (BB). Pasture versus Grain origins of the beef had little influence on participants' responses. This trial was registered at ClinicalTrials.gov as NCT04545398.

Gossypol detoxification in the rumen and Helicoverpa armigera larvae: A review

Gossypol, a phenolic compound found in the cotton plant, is widely distributed in cottonseed by-products. Although ruminant animals are believed to be more tolerant of gossypol toxicity than monogastric animals due to rumen microbial fermentation, the actual mechanisms of detoxification remain unclear. In contrast, the metabolic detoxification of gossypol by Helicoverpa armigera (Lepidoptera: Noctuidae) larvae has achieved great advances. The present review discusses the clinical signs of gossypol in ruminant animals, as well as summarizing advances in the study of gossypol detoxification in the rumen. It also examines the regulatory roles of several key enzymes in gossypol detoxification and transformation known in H. armigera. With the rapid development of modern molecular biotechnology and -omics technology strategies, evidence increasingly indicates that research into the biological degradation of gossypol in H. armigera larvae and some microbes, in terms of these key enzymes, could provide scientific insights that would underpin future work on microbial gossypol detoxification in the rumen, with the ultimate aim of further alleviating gossypol toxicity in ruminant animals.

Metabolic Regulations by lncRNA, miRNA, and ceRNA Under Grass-Fed and Grain-Fed Regimens in Angus Beef Cattle

Beef cattle raised under grass-fed and grain-fed have many differences, including metabolic efficiency and meat quality. To investigate these two regimens' intrinsic influence on beef cattle, we used high-throughput sequencing and metabolomics analyses to explore differentially expressed genes (DEGs) and metabolimic networks in the liver. A total of 200 DEGs, 76 differentially expressed miRNAs (DEmiRNAs), and two differentially expressed lncRNAs (DElncRNAs) were detected between regimen groups. Metabolic processes and pathways enriched functional genes including target genes of miRNAs and lncRNAs. We found that many genes were involved in energy, retinol and cholesterol metabolism, and bile acid synthesis. Combined with metabolites such as low glucose concentration, high cholesterol concentration, and increased primary bile acid concentration, these genes were mainly responsible for lowering intramuscular fat, low cholesterol, and yellow meat in grass-fed cattle. Additionally, we identified two lncRNAs and eight DEGs as potential competing endogenous RNAs (ceRNAs) to bind miRNAs by the interaction network analysis. These results revealed that the effects of two feeding regimens on beef cattle were mainly induced by gene expression changes in metabolic pathways mediated via lncRNAs, miRNAs, and ceRNAs, and contents of metabolites in the liver. It may provide a clue on feeding regimens inducing the metabolic regulations.

Diet induced the change of mtDNA copy number and metabolism in Angus cattle

Background

Grass-fed and grain-fed Angus cattle differ in the diet regimes. However, the intricate mechanisms of different beef quality and other phenotypes induced by diet differences are still unclear. Diet affects mitochondrial function and dynamic behavior in response to changes in energy demand and supply. In this study, we examined the mtDNA copy number, mitochondria-related genes expression, and metabolic biomarkers in grass-fed and grain-fed Angus cattle.

Results

We found that the grass-fed group had a higher mtDNA copy number than the grain-fed group. Among different tissues, the mtDNA copy number was the highest in the liver than muscle, rumen, and spleen. Based on the transcriptome of the four tissues, a lower expression of mtDNA-encoded genes in the grass-fed group compared to the grain-fed group was discovered. For the mitochondria-related nuclear genes, however, most of them were significantly down-regulated in the muscle of the grass-fed group and up-regulated in the other three tissues. In which, COX6A2, POLG2, PPIF, DCN, and NDUFA12, involving in ATP synthesis, mitochondrial replication, transcription, and maintenance, might contribute to the alterations of mtDNA copy number and gene expression. Meanwhile, 40 and 23 metabolic biomarkers were identified in the blood and muscle of the grain-fed group compared to a grass-fed group, respectively. Integrated analysis of the altered metabolites and gene expression revealed the high expression level of MDH1 in the grain-fed group might contribute to the mitochondrial NADH oxidation and spermidine metabolism for adapting the deletion mtDNA copy number.

Conclusions

Overall, the study may provide further deep insight into the adaptive and regulatory modulations of the mitochondrial function in response to different feeding systems in Angus cattle.

DNA methylation, microRNA expression profiles and their relationships with transcriptome in grass-fed and grain-fed Angus cattle rumen tissue

Rumen is an organ for supplying nutrients for the growth and production of bovine, which might function differently under grass-fed and grain-fed regimens considering the association of gene expression, DNA methylation, and microRNA expression. The objective of this study was to explore the potential mechanism influencing rumen function of grass-fed and grain-fed animals. Methylated DNA binding domain sequencing (MBD-Seq) and microRNA-Seq were respectively utilized to detect the DNA methylation and microRNA expression in rumen tissue of grass-fed and grain-fed Angus cattle. Combined analysis revealed that the expression of the differentially expressed genes ADAMTS3 and ENPP3 was correlated with the methylation abundance of the corresponding differentially methylated regions (DMRs) inside these two genes, and these two genes were reported to be respectively involved in biosynthesis and regulation of glycosyltransferase activity; the differentially expressed microRNA bta-mir-122 was predicted to possibly target the differentially expressed genes OCLN and RBM47, potentially affecting the rumen function; the microRNA bta-mir-655 was exclusively detected in grain-fed group; its targets were significantly enriched in insulin and TGF-beta signaling pathways, which might worked together to regulate the function of rumen, resulting in different characteristics between grass-fed and grain-fed cattle. Collectively, our results provided insights into understanding the mechanisms determining rumen function and unraveled the biological basis underlying the economic traits to improve the productivity of animals.

A Nutrigenomics Approach Using RNA Sequencing Technology to Study Nutrient-Gene Interactions in Agricultural Animals

Thorough understanding of animal gene expression driven by dietary nutrients can be regarded as a bottom line of advanced animal nutrition research. Nutrigenomics (including transcriptomics) studies the effects of dietary nutrients on cellular gene expression and, ultimately, phenotypic changes in living organisms. Transcriptomics can be applied to investigate animal tissue transcriptomes at a defined nutritional state, which can provide a holistic view of intracellular RNA expression. As a novel transcriptomics approach, RNA sequencing (RNA-Seq) technology can monitor all gene expressions simultaneously in response to dietary intervention. The principle and history of RNA-Seq are briefly reviewed, and its 3 principal steps are described in this article. Application of RNA-Seq in different areas of animal nutrition research is summarized. Lastly, the application of RNA-Seq in swine science and nutrition is also reviewed. In short, RNA-Seq holds significant potential to be employed for better understanding the nutrient-gene interactions in agricultural animals.

MicroRNA Sequencing Reveals the Effect of Different Levels of Non-Fibrous Carbohydrate/Neutral Detergent Fiber on Rumen Development in Calves

Rumen development in calves is affected by many factors, including dietary composition. MicroRNAs (miRNAs) are known to function in the development of the rumen in cattle, what is not known is how these miRNAs function in rumen development of calves fed with high and low ratios of non-fibrous carbohydrate (NFC)/neutral detergent fiber (NDF). A total of six healthy Charolais hybrids bull calves of similar weight were divided into two groups; three calves were fed a mixed diet with NFC/NDF = 1.35 (H group), and three were fed a mixed diet with NFC/NDF = 0.80 (L group). After 105 days on the diet, calves were sacrificed and rumen tissues were collected. Tissues were subjected to histological observation and miRNA expression analysis. Functional enrichment analysis was conducted on the target genes of the miRNAs. Targeting and regulatory relationships were verified by luciferase reporter assay and quantitative PCR (qPCR). We found that the length of rumen papilla in the L group was significantly greater than that in the H group, while the width of rumen papilla in H group was significantly greater than that that in L group. We identified 896 miRNAs; 540 known miRNAs, and 356 novel predicted miRNAs. After statistical testing, we identified 24 differentially expressed miRNAs (DEmiRNAs). miRNA-mRNA-cluster network analysis and literature reviews revealed that cell proliferation, differentiation, physical and nutrient stimuli processes participate in rumen development under different NFC/NDF levels. The regulatory relationships between three DEmiRNAs and five target genes were verified by examining the levels of expression. The binding sites on bta-miR-128 for the peroxisome proliferator activated receptor gamma (PPARG) and solute carrier family 16 member 1 (SLC16A1) genes were investigated using a dual luciferase assay. The results of this study provide insight into the role of miRNAs in rumen development in calves under different NFC/NDF levels.

Comprehensive Transcriptional Profiling of the Gastrointestinal Tract of Ruminants from Birth to Adulthood Reveals Strong Developmental Stage Specific Gene Expression

One of the most significant physiological challenges to neonatal and juvenile ruminants is the development and establishment of the rumen. Using a subset of RNA-Seq data from our high-resolution atlas of gene expression in sheep (Ovis aries) we have provided the first comprehensive characterization of transcription of the entire gastrointestinal (GI) tract during the transition from pre-ruminant to ruminant. The dataset comprises 164 tissue samples from sheep at four different time points (birth, one week, 8 weeks and adult). Using network cluster analysis we illustrate how the complexity of the GI tract is reflected in tissue- and developmental stage-specific differences in gene expression. The most significant transcriptional differences between neonatal and adult sheep were observed in the rumen complex. Comparative analysis of gene expression in three GI tract tissues from age-matched sheep and goats revealed species-specific differences in genes involved in immunity and metabolism. This study improves our understanding of the transcriptomic mechanisms involved in the transition from pre-ruminant to ruminant by identifying key genes involved in immunity, microbe recognition and metabolism. The results form a basis for future studies linking gene expression with microbial colonization of the developing GI tract and provide a foundation to improve ruminant efficiency and productivity through identifying potential targets for novel therapeutics and gene editing.

Performance Measurement and Comparative Transcriptome Analysis Revealed the Efforts on Hybrid Improvement of Qinchuan Cattle

Crossbreeding can provide productive gains through heterosis, however, surveys about the effects of crossbreeding through global transcriptomic sequencing are few. This study revealed that Angus × Qinchuan cattle (AQF) have improved performance characteristics compared to Qinchuan cattle (QCF). We performed RNA-seq on the subcutaneous fat tissue of QCF and AQF. More than 42.2 million clean reads were obtained in each sample. We detected 40 and 21 breed-specific highly expressed genes (FPKM > 500) in QCF and AQF, respectively. Furthermore, a total of 353 differentially expressed genes (DEGs, |log2 ratio| ≥ 1 and Probability ≥ 0.8) were found between these two groups, of which 227 genes were upregulated in AQF and 126 genes were upregulated in QCF. Functional enrichment analyses showed that breed-specific highly expressed genes and DEGs were closely related to terms such as development in AQF, and adaption or immune in QCF. In addition, we also identified the novel transcript units, alternative splicing events, single-nucleotide polymorphisms and Indels. Our results revealed differences in inherent characteristics and genetic differences when comparing QCF with AQF.

Infrequently expressed miRNAs influence survival after diagnosis with colorectal cancer

Half of miRNAs expressed in colorectal tissue are expressed < 50% of the population. Many infrequently expressed miRNAs have low levels of expression. We hypothesize that less frequently expressed miRNAs, when expressed at higher levels, influence both disease stage and survival after diagnosis with colorectal cancer (CRC); low levels of expression may be background noise.

We examine 304 infrequently expressed miRNAs in 1893 population-based cases of CRC with paired carcinoma and normal mucosa miRNA profiles. We evaluate miRNAs with disease stage and survival after adjusting for age, study center, sex, MSI status, and AJCC stage. These miRNAs were further evaluated with RNA-Seq data to identify miRNA::mRNA associations that may provide insight into the functionality of miRNAs.

Eleven miRNAs were associated with advanced disease stage among colon cancer patients (Q value = 0.10). Eight infrequently expressed miRNAs influenced survival if highly expressed in overall CRC. Of these, five increased likelihood of dying if they were highly expressed, i.e. miR-124-3p, miR-143-5p, miR-145-3p, miR31-5p, and miR-99b-5p, while three were associated with better survival if highly expressed, i.e. miR-362-5p, miR-374a-5p, and miR-590-5p. Thirteen miRNAs infrequently expressed in colon-specific carcinoma tissue were associated with CRC survival if highly expressed. Evaluation of miRNAs::mRNA associations showed that mRNA expression influenced by infrequently expressed miRNA contributed to networks and pathways shown to influence disease progression and prognosis.

Our large study enabled us to examine the implications of infrequently expressed miRNAs after removal of background noise. These results require replication in other studies. Confirmation of our findings in other studies could lead to important markers for prognosis.

Infrequently expressed miRNAs in colorectal cancer tissue and tumor molecular phenotype

We have previously shown that commonly expressed miRNAs influenced tumor molecular phenotype in colorectal cancer. We hypothesize that infrequently expressed miRNAs, when showing higher levels of expression, help to define tumor molecular phenotype. In this study, we examine 304 miRNAs expressed in at least 30 individuals, but in <50% of the population and with a mean level of expression above 1.0 relative florescent unit. We examine associations in 1893 individuals who have the tumor molecular phenotype data as well as miRNA expression levels for both carcinoma and normal colorectal tissue. We compare miRNAs uniquely associated with tumor molecular phenotype to the RNAseq data to identify genes associated with these miRNAs. This information is used to further identify unique pathways associated with tumor molecular phenotypes of TP53-mutated, KRAS-mutated, CpG island methylator phenotype and microsatellite instability tumors. Thirty-seven miRNAs were uniquely associated with TP53-mutated tumors; 30 of these miRNAs had higher level of expression in TP53-mutated tumors, while seven had lower levels of expression. Of the 34 miRNAs associated with CpG island methylator phenotype-high tumors, 16 were more likely to have a CpG island methylator phenotype-high tumor and 19 were less likely to be CpG island methylator phenotype-high. For microsatellite instability, 13 of the 22 infrequently expressed miRNAs were significantly less likely to be expressed in microsatellite unstable tumors. KRAS-mutated tumors were not associated with any miRNAs after adjustment for multiple comparisons. Of the dysregulated miRNAs, 17 were more likely to be TP53-mutated tumors while simultaneously being less likely to be CpG island methylator phenotype-high and/or microsatellite instability tumors. Genes regulated by these miRNAs were involved in numerous functions and pathways that influence cancer risk and progression. In summary, some infrequently expressed miRNAs, when expressed at higher levels, appear to have significant biological meaning in terms of tumor molecular phenotype and gene expression profiles.

Dietary intake alters gene expression in colon tissue: possible underlying mechanism for the influence of diet on disease

Supplemental Digital Content is available in the text.

Background

Although the association between diet and disease is well documented, the biologic mechanisms involved have not been entirely elucidated. In this study, we evaluate how dietary intake influences gene expression to better understand the underlying mechanisms through which diet operates.

Methods

We used data from 144 individuals who had comprehensive dietary intake and gene expression data from RNAseq using normal colonic mucosa. Using the DESeq2 statistical package, we identified genes that showed statistically significant differences in expression between individuals in high-intake and low-intake categories for several dietary variables of interest adjusting for age and sex. We examined total calories, total fats, vegetable protein, animal protein, carbohydrates, trans-fatty acids, mutagen index, red meat, processed meat, whole grains, vegetables, fruits, fiber, folate, dairy products, calcium, and prudent and western dietary patterns.

Results

Using a false discovery rate of less than 0.1, meat-related foods were statistically associated with 68 dysregulated genes, calcium with three dysregulated genes, folate with four dysregulated genes, and nonmeat-related foods with 65 dysregulated genes. With a more stringent false discovery rate of less than 0.05, there were nine meat-related dysregulated genes and 23 nonmeat-related genes. Ingenuity pathway analysis identified three major networks among genes identified as dysregulated with respect to meat-related dietary variables and three networks among genes identified as dysregulated with respect to nonmeat-related variables. The top networks (Ingenuity Pathway Analysis network score >30) associated with meat-related genes were (i) cancer, organismal injury, and abnormalities, tumor morphology, and (ii) cellular function and maintenance, cellular movement, cell death, and survival. Among genes related to nonmeat consumption variables, the top networks were (i) hematological system development and function, nervous system development and function, tissue morphology and (ii) connective tissue disorders, organismal injury, and abnormalities.

Conclusion

Several dietary factors were associated with gene expression in our data. These findings provide insight into the possible mechanisms by which diet may influence disease processes.

Large-scale gene co-expression network as a source of functional annotation for cattle genes

Background

Genome sequencing and subsequent gene annotation of genomes has led to the elucidation of many genes, but in vertebrates the actual number of protein coding genes are very consistent across species (~20,000). Seven years after sequencing the cattle genome, there are still genes that have limited annotation and the function of many genes are still not understood, or partly understood at best. Based on the assumption that genes with similar patterns of expression across a vast array of tissues and experimental conditions are likely to encode proteins with related functions or participate within a given pathway, we constructed a genome-wide Cattle Gene Co-expression Network (CGCN) using 72 microarray datasets that contained a total of 1470 Affymetrix Genechip Bovine Genome Arrays that were retrieved from either NCBI GEO or EBI ArrayExpress.

Results

The total of 16,607 probe sets, which represented 11,397 genes, with unique Entrez ID were consolidated into 32 co-expression modules that contained between 29 and 2569 probe sets. All of the identified modules showed strong functional enrichment for gene ontology (GO) terms and Reactome pathways. For example, modules with important biological functions such as response to virus, response to bacteria, energy metabolism, cell signaling and cell cycle have been identified. Moreover, gene co-expression networks using “guilt-by-association” principle have been used to predict the potential function of 132 genes with no functional annotation. Four unknown Hub genes were identified in modules highly enriched for GO terms related to leukocyte activation (LOC509513), RNA processing (LOC100848208), nucleic acid metabolic process (LOC100850151) and organic-acid metabolic process (MGC137211). Such highly connected genes should be investigated more closely as they likely to have key regulatory roles.

Conclusions

We have demonstrated that the CGCN and its corresponding regulons provides rich information for experimental biologists to design experiments, interpret experimental results, and develop novel hypothesis on gene function in this poorly annotated genome. The network is publicly accessible at http://www.animalgenome.org/cgi-bin/host/reecylab/d.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-016-3176-2) contains supplementary material, which is available to authorized users.

Epithelial, metabolic and innate immunity transcriptomic signatures differentiating the rumen from other sheep and mammalian gastrointestinal tract tissues

Background. Ruminants are successful herbivorous mammals, in part due to their specialized forestomachs, the rumen complex, which facilitates the conversion of feed to soluble nutrients by micro-organisms. Is the rumen complex a modified stomach expressing new epithelial (cornification) and metabolic programs, or a specialised stratified epithelium that has acquired new metabolic activities, potentially similar to those of the colon? How has the presence of the rumen affected other sections of the gastrointestinal tract (GIT) of ruminants compared to non-ruminants? Methods. Transcriptome data from 11 tissues covering the sheep GIT, two stratified epithelial and two control tissues, was analysed using principal components to cluster tissues based on gene expression profile similarity. Expression profiles of genes along the sheep GIT were used to generate a network to identify genes enriched for expression in different compartments of the GIT. The data from sheep was compared to similar data sets from two non-ruminants, pigs (closely related) and humans (more distantly related). Results. The rumen transcriptome clustered with the skin and tonsil, but not the GIT transcriptomes, driven by genes from the epidermal differentiation complex, and genes encoding stratified epithelium keratins and innate immunity proteins. By analysing all of the gene expression profiles across tissues together 16 major clusters were identified. The strongest of these, and consistent with the high turnover rate of the GIT, showed a marked enrichment of cell cycle process genes (P = 1.4 E-46), across the whole GIT, relative to liver and muscle, with highest expression in the caecum followed by colon and rumen. The expression patterns of several membrane transporters (chloride, zinc, nucleosides, amino acids, fatty acids, cholesterol and bile acids) along the GIT was very similar in sheep, pig and humans. In contrast, short chain fatty acid uptake and metabolism appeared to be different between the species and different between the rumen and colon in sheep. The importance of nitrogen and iodine recycling in sheep was highlighted by the highly preferential expression of SLC14A1-urea (rumen), RHBG-ammonia (intestines) and SLC5A5-iodine (abomasum). The gene encoding a poorly characterized member of the maltase-glucoamylase family (MGAM2), predicted to play a role in the degradation of starch or glycogen, was highly expressed in the small and large intestines. Discussion. The rumen appears to be a specialised stratified cornified epithelium, probably derived from the oesophagus, which has gained some liver-like and other specialized metabolic functions, but probably not by expression of pre-existing colon metabolic programs. Changes in gene transcription downstream of the rumen also appear have occurred as a consequence of the evolution of the rumen and its effect on nutrient composition flowing down the GIT.