Genome-Wide microRNA Binding Site Variation between Extinct Wild Aurochs and Modern Cattle Identifies Candidate microRNA-Regulated Domestication Genes

peer-reviewedThe domestication of cattle from the now-extinct wild aurochs (Bos primigenius) involved selection for physiological and behavioral traits, with underlying genetic factors that remain largely unknown. Non-coding microRNAs have emerged as key regulators of the spatio-temporal expression of target genes controlling mammalian growth and development, including in livestock species. During the domestication process, selection of mutational changes in miRNAs and/or miRNA binding sites could have provided a mechanism to generate some of the traits that differentiate domesticated cattle from wild aurochs. To investigate this, we analyzed the open reading frame DNA sequence of 19,994 orthologous protein-coding gene pairs from extant Bos taurus genomes and a single extinct B. primigenius genome. We identified miRNA binding site polymorphisms in the 3′ UTRs of 1,620 of these orthologous genes. These 1,620 genes with altered miRNA binding sites between the B. taurus and B. primigenius lineages represent candidate domestication genes. Using a novel Score Site ratio metric we have ranked these miRNA-regulated genes according to the extent of divergence between miRNA binding site presence, frequency and copy number between the orthologous genes from B. taurus and B. primigenius. This provides an unbiased approach to identify cattle genes that have undergone the most changes in miRNA binding (i.e., regulation) between the wild aurochs and modern-day cattle breeds. In addition, we demonstrate that these 1,620 candidate domestication genes are enriched for roles in pigmentation, fertility, neurobiology, metabolism, immunity and production traits (including milk quality and feed efficiency). Our findings suggest that directional selection of miRNA regulatory variants was important in the domestication and subsequent artificial selection that gave rise to modern taurine cattle

Long-term stability of RNA in post-mortem bovine skeletal muscle, liver and subcutaneous adipose tissues

peer-reviewedBackground: 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

RNA Sequencing (RNA-Seq) Reveals Extremely Low Levels of Reticulocyte-Derived Globin Gene Transcripts in Peripheral Blood From Horses (Equus caballus) and Cattle (Bos taurus)

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

Antigen stimulation of peripheral blood mononuclear cells from Mycobacterium bovis infected cattle yields evidence for a novel gene expression program

<p>Abstract</p> <p>Background</p> <p>Bovine tuberculosis (BTB) caused by <it>Mycobacterium bovis </it>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 <it>in vivo </it>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 <it>M. bovis</it>. A functional genomics approach was used to examine the immune response of BTB-infected (<it>n </it>= 6) and healthy control (<it>n </it>= 6) cattle to stimulation with bovine tuberculin (purified protein derivative – PPD-b) <it>in vitro</it>. PBMC were harvested before, and at 3 h and 12 h post <it>in vitro </it>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.</p> <p>Results</p> <p>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 (<it>P </it>≤ 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 (<it>n </it>= 8) versus control animals (<it>n </it>= 8) after stimulation with bovine tuberculin.</p> <p>Conclusion</p> <p>The PBMC from BTB-infected animals exhibit different transcriptional profiles compared with PBMC from healthy control animals in response to <it>M. bovis </it>antigen stimulation, providing evidence of a novel gene expression program due to <it>M. bovis </it>exposure.</p

Key Hub and Bottleneck Genes Differentiate the Macrophage Response to Virulent and Attenuated Mycobacterium bovis

Mycobacterium bovis is an intracellular pathogen that causes tuberculosis in cattle. Following infection, the pathogen resides and persists inside host macrophages by subverting host immune responses via a diverse range of mechanisms. Here, a high-density bovine microarray platform was used to examine the bovine monocyte-derived macrophage transcriptome response to M. bovis infection relative to infection with the attenuated vaccine strain, M. bovis Bacille Calmette–Guérin. Differentially expressed genes were identified (adjusted P-value ≤0.01) and interaction networks generated across an infection time course of 2, 6, and 24 h. The largest number of biological interactions was observed in the 24-h network, which exhibited scale-free network properties. The 24-h network featured a small number of key hub and bottleneck gene nodes, including IKBKE, MYC, NFKB1, and EGR1 that differentiated the macrophage response to virulent and attenuated M. bovis strains, possibly via the modulation of host cell death mechanisms. These hub and bottleneck genes represent possible targets for immuno-modulation of host macrophages by virulent mycobacterial species that enable their survival within a hostile environment

Global Gene Expression and Systems Biology Analysis of Bovine Monocyte-Derived Macrophages in Response to In Vitro Challenge with Mycobacterium bovis

peer-reviewedBackground 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)

Innate gene repression associated with Mycobacterium bovis infection in cattle: toward a gene signature of disease

<p>Abstract</p> <p>Background</p> <p>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 <it>Mycobacterium bovis </it>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 <it>M. bovis </it>infected animals <it>in vivo</it>.</p> <p>Results</p> <p>In total, 378 gene features were differentially expressed at the <it>P </it>≤ 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 (<it>TLR2</it>) and <it>TLR4 </it>genes, lack of differential expression of indicator adaptive immune gene transcripts (<it>IFNG, IL2, IL4</it>), and lower <it>BOLA </it>major histocompatibility complex – class I (<it>BOLA</it>) and class II (<it>BOLA-DRA</it>) 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 (<it>n </it>= 8 per group) by real time quantitative reverse transcription PCR.</p> <p>Conclusion</p> <p>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 <it>in vivo </it>infection, in the absence of exogenous antigenic stimulation.</p

DNA sequence polymorphisms in a panel of eight candidate bovine imprinted genes and their association with performance traits in Irish Holstein-Friesian cattle

peer-reviewedBackground: Studies in mice and humans have shown that imprinted genes, whereby expression from one of the two parentally inherited alleles is attenuated or completely silenced, have a major effect on mammalian growth, metabolism and physiology. More recently, investigations in livestock species indicate that genes subject to this type of epigenetic regulation contribute to, or are associated with, several performance traits, most notably muscle mass and fat deposition. In the present study, a candidate gene approach was adopted to assess 17 validated single nucleotide polymorphisms (SNPs) and their association with a range of performance traits in 848 progeny-tested Irish Holstein-Friesian artificial insemination sires. These SNPs are located proximal to, or within, the bovine orthologs of eight genes (CALCR, GRB10, PEG3, PHLDA2, RASGRF1, TSPAN32, ZIM2 and ZNF215) that have been shown to be imprinted in cattle or in at least one other mammalian species (i.e. human/mouse/pig/sheep). Results: Heterozygosities for all SNPs analysed ranged from 0.09 to 0.46 and significant deviations from Hardy-Weinberg proportions (P ≤ 0.01) were observed at four loci. Phenotypic associations (P ≤ 0.05) were observed between nine SNPs proximal to, or within, six of the eight analysed genes and a number of performance traits evaluated, including milk protein percentage, somatic cell count, culled cow and progeny carcass weight, angularity, body conditioning score, progeny carcass conformation, body depth, rump angle, rump width, animal stature, calving difficulty, gestation length and calf perinatal mortality. Notably, SNPs within the imprinted paternally expressed gene 3 (PEG3) gene cluster were associated (P ≤ 0.05) with calving, calf performance and fertility traits, while a single SNP in the zinc finger protein 215 gene (ZNF215) was associated with milk protein percentage (P ≤ 0.05), progeny carcass weight (P ≤ 0.05), culled cow carcass weight (P ≤ 0.01), angularity (P ≤ 0.01), body depth (P ≤ 0.01), rump width (P ≤ 0.01) and animal stature (P ≤ 0.01). Conclusions: Of the eight candidate bovine imprinted genes assessed, DNA sequence polymorphisms in six of these genes (CALCR, GRB10, PEG3, RASGRF1, ZIM2 and ZNF215) displayed associations with several of the phenotypes included for analyses. The genotype-phenotype associations detected here are further supported by the biological function of these six genes, each of which plays important roles in mammalian growth, development and physiology. The associations between SNPs within the imprinted PEG3 gene cluster and traits related to calving, calf performance and gestation length suggest that this domain on chromosome 18 may play a role regulating pre-natal growth and development and fertility. SNPs within the bovine ZNF215 gene were associated with bovine growth and body conformation traits and studies in humans have revealed that the human ZNF215 ortholog belongs to the imprinted gene cluster associated with Beckwith-Wiedemann syndrome--a genetic disorder characterised by growth abnormalities. Similarly, the data presented here suggest that the ZNF215 gene may have an important role in regulating bovine growth. Collectively, our results support previous work showing that (candidate) imprinted genes/loci contribute to heritable variation in bovine performance traits and suggest that DNA sequence polymorphisms within these genes/loci represents an important reservoir of genomic markers for future genetic improvement of dairy and beef cattle populations.Department of Agriculture, Fisheries and Food Ireland - Research Stimulus Fund (project numbers: RSF-06-406, RSF-06-0353 and RSF-06-0409); Science Foundation Ireland - Investigator Programme Grants(SFI/01/F.1/B028; SFI/08/IN.1/B1931, 07/SRC/B1156 (MPM)

MicroRNA profiling of the bovine alveolar macrophage response to Mycobacterium bovis infection suggests pathogen survival is enhanced by microRNA regulation of endocytosis and lysosome trafficking

peer-reviewedMycobacterium 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