A Population Genomics Analysis of the Native Irish Galway Sheep Breed.

The Galway sheep population is the only native Irish sheep breed and this livestock genetic resource is currently categorised as 'at-risk'. In the present study, comparative population genomics analyses of Galway sheep and other sheep populations of European origin were used to investigate the microevolution and recent genetic history of the breed. These analyses support the hypothesis that British Leicester sheep were used in the formation of the Galway. When compared to conventional and endangered breeds, the Galway breed was intermediate in effective population size, genomic inbreeding and runs of homozygosity. This indicates that, although the Galway breed is declining, it is still relatively genetically diverse and that conservation and management plans informed by genomic information may aid its recovery. The Galway breed also exhibited distinct genomic signatures of artificial or natural selection when compared to other breeds, which highlighted candidate genes that may be involved in production and health traits.

Domestication of cattle: Two or three events?

Cattle have been invaluable for the transition of human society from nomadic hunter-gatherers to sedentary farming communities throughout much of Europe, Asia and Africa since the earliest domestication of cattle more than 10,000 years ago. Although current understanding of relationships among ancestral populations remains limited, domestication of cattle is thought to have occurred on two or three occasions, giving rise to the taurine (Bos taurus) and indicine (Bos indicus) species that share the aurochs (Bos primigenius) as common ancestor ~250,000 years ago. Indicine and taurine cattle were domesticated in the Indus Valley and Fertile Crescent, respectively; however, an additional domestication event for taurine in the Western Desert of Egypt has also been proposed. We analysed medium density Illumina Bovine SNP array (~54,000 loci) data across 3,196 individuals, representing 180 taurine and indicine populations to investigate population structure within and between populations, and domestication and demographic dynamics using approximate Bayesian computation (ABC). Comparative analyses between scenarios modelling two and three domestication events consistently favour a model with only two episodes and suggest that the additional genetic variation component usually detected in African taurine cattle may be explained by hybridization with local aurochs in Africa after the domestication of taurine cattle in the Fertile Crescent. African indicine cattle exhibit high levels of shared genetic variation with Asian indicine cattle due to their recent divergence and with African taurine cattle through relatively recent gene flow. Scenarios with unidirectional or bidirectional migratory events between European taurine and Asian indicine cattle are also plausible, although further studies are needed to disentangle the complex human-mediated dispersion patterns of domestic cattle. This study therefore helps to clarify the effect of past demographic history on the genetic variation of modern cattle, providing a basis for further analyses exploring alternative migratory routes for early domestic populations.

Genomic Characterisation of the Indigenous Irish Kerry Cattle Breed.

Kerry cattle are an endangered landrace heritage breed of cultural importance to Ireland. In the present study we have used genome-wide SNP array data to evaluate genomic diversity within the Kerry population and between Kerry cattle and other European breeds. Patterns of genetic differentiation and gene flow among breeds using phylogenetic trees with ancestry graphs highlighted historical gene flow from the British Shorthorn breed into the ancestral population of modern Kerry cattle. Principal component analysis (PCA) and genetic clustering emphasised the genetic distinctiveness of Kerry cattle relative to comparator British and European cattle breeds. Modelling of genetic effective population size (Ne) revealed a demographic trend of diminishing Ne over time and that recent estimated Ne values for the Kerry breed may be less than the threshold for sustainable genetic conservation. In addition, analysis of genome-wide autozygosity (FROH) showed that genomic inbreeding has increased significantly during the 20 years between 1992 and 2012. Finally, signatures of selection revealed genomic regions subject to natural and artificial selection as Kerry cattle adapted to the climate, physical geography and agro-ecology of southwest Ireland.

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

The 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.

The application of transcriptomic data in the authentication of beef derived from contrasting production systems.

BACKGROUND: Differences between cattle production systems can influence the nutritional and sensory characteristics of beef, in particular its fatty acid (FA) composition. As beef products derived from pasture-based systems can demand a higher premium from consumers, there is a need to understand the biological characteristics of pasture produced meat and subsequently to develop methods of authentication for these products. Here, we describe an approach to authentication that focuses on differences in the transcriptomic profile of muscle from animals finished in different systems of production of practical relevance to the Irish beef industry. The objectives of this study were to identify a panel of differentially expressed (DE) genes/networks in the muscle of cattle raised outdoors on pasture compared to animals raised indoors on a concentrate based diet and to subsequently identify an optimum panel which can classify the meat based on a production system. RESULTS: A comparison of the muscle transcriptome of outdoor/pasture-fed and Indoor/concentrate-fed cattle resulted in the identification of 26 DE genes. Functional analysis of these genes identified two significant networks (1: Energy Production, Lipid Metabolism, Small Molecule Biochemistry; and 2: Lipid Metabolism, Molecular Transport, Small Molecule Biochemistry), both of which are involved in FA metabolism. The expression of selected up-regulated genes in the outdoor/pasture-fed animals correlated positively with the total n-3 FA content of the muscle. The pathway and network analysis of the DE genes indicate that peroxisome proliferator-activated receptor (PPAR) and FYN/AMPK could be implicit in the regulation of these alterations to the lipid profile. In terms of authentication, the expression profile of three DE genes (ALAD, EIF4EBP1 and NPNT) could almost completely separate the samples based on production system (95 % authentication for animals on pasture-based and 100 % for animals on concentrate- based diet) in this context. CONCLUSIONS: The majority of DE genes between muscle of the outdoor/pasture-fed and concentrate-fed cattle were related to lipid metabolism and in particular ß-oxidation. In this experiment the combined expression profiles of ALAD, EIF4EBP1 and NPNT were optimal in classifying the muscle transcriptome based on production system. Given the overall lack of comparable studies and variable concordance with those that do exist, the use of transcriptomic data in authenticating production systems requires more exploration across a range of contexts and breeds.

Genome sequencing of the extinct Eurasian wild aurochs, Bos primigenius, illuminates the phylogeography and evolution of cattle.

BACKGROUND: Domestication of the now-extinct wild aurochs, Bos primigenius, gave rise to the two major domestic extant cattle taxa, B. taurus and B. indicus. While previous genetic studies have shed some light on the evolutionary relationships between European aurochs and modern cattle, important questions remain unanswered, including the phylogenetic status of aurochs, whether gene flow from aurochs into early domestic populations occurred, and which genomic regions were subject to selection processes during and after domestication. Here, we address these questions using whole-genome sequencing data generated from an approximately 6,750-year-old British aurochs bone and genome sequence data from 81 additional cattle plus genome-wide single nucleotide polymorphism data from a diverse panel of 1,225 modern animals. RESULTS: Phylogenomic analyses place the aurochs as a distinct outgroup to the domestic B. taurus lineage, supporting the predominant Near Eastern origin of European cattle. Conversely, traditional British and Irish breeds share more genetic variants with this aurochs specimen than other European populations, supporting localized gene flow from aurochs into the ancestors of modern British and Irish cattle, perhaps through purposeful restocking by early herders in Britain. Finally, the functions of genes showing evidence for positive selection in B. taurus are enriched for neurobiology, growth, metabolism and immunobiology, suggesting that these biological processes have been important in the domestication of cattle. CONCLUSIONS: This work provides important new information regarding the origins and functional evolution of modern cattle, revealing that the interface between early European domestic populations and wild aurochs was significantly more complex than previously thought.

Analysis of the Bovine Monocyte-Derived Macrophage Response to Mycobacterium avium Subspecies Paratuberculosis Infection Using RNA-seq.

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 6-h infection time course with non-infected controls. We observed 245 and 574 differentially expressed (DE) genes in MAP-infected versus non-infected control samples (adjusted P value ≤0.05) at 2 and 6 h post-infection, respectively. Functional analyses of these DE 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.

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.

RNA-seq Transcriptional Profiling of Peripheral Blood Leukocytes from Cattle Infected with Mycobacterium bovis.

Bovine tuberculosis, 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.

Whole-transcriptome, high-throughput RNA sequence analysis of the bovine macrophage response to Mycobacterium bovis infection in vitro.

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.

Global gene expression and systems biology analysis of bovine monocyte-derived macrophages in response to in vitro challenge with Mycobacterium bovis.

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.

Genome-wide transcriptional profiling of peripheral blood leukocytes from cattle infected with Mycobacterium bovis reveals suppression of host immune genes.

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.

A complete mitochondrial genome sequence from a mesolithic wild aurochs (Bos primigenius).

BACKGROUND: The derivation of domestic cattle from the extinct wild aurochs (Bos primigenius) has been well-documented by archaeological and genetic studies. Genetic studies point towards the Neolithic Near East as the centre of origin for Bos taurus, with some lines of evidence suggesting possible, albeit rare, genetic contributions from locally domesticated wild aurochsen across Eurasia. Inferences from these investigations have been based largely on the analysis of partial mitochondrial DNA sequences generated from modern animals, with limited sequence data from ancient aurochsen samples. Recent developments in DNA sequencing technologies, however, are affording new opportunities for the examination of genetic material retrieved from extinct species, providing new insight into their evolutionary history. Here we present DNA sequence analysis of the first complete mitochondrial genome (16,338 base pairs) from an archaeologically-verified and exceptionally-well preserved aurochs bone sample. METHODOLOGY: DNA extracts were generated from an aurochs humerus bone sample recovered from a cave site located in Derbyshire, England and radiocarbon-dated to 6,738+/-68 calibrated years before present. These extracts were prepared for both Sanger and next generation DNA sequencing technologies (Illumina Genome Analyzer). In total, 289.9 megabases (22.48%) of the post-filtered DNA sequences generated using the Illumina Genome Analyzer from this sample mapped with confidence to the bovine genome. A consensus B. primigenius mitochondrial genome sequence was constructed and was analysed alongside all available complete bovine mitochondrial genome sequences. CONCLUSIONS: For all nucleotide positions where both Sanger and Illumina Genome Analyzer sequencing methods gave high-confidence calls, no discrepancies were observed. Sequence analysis reveals evidence of heteroplasmy in this sample and places this mitochondrial genome sequence securely within a previously identified aurochsen haplogroup (haplogroup P), thus providing novel insights into pre-domestic patterns of variation. The high proportion of authentic, endogenous aurochs DNA preserved in this sample bodes well for future efforts to determine the complete genome sequence of a wild ancestor of domestic cattle.

Transcriptional profiling of cattle infected with Trypanosoma congolense highlights gene expression signatures underlying trypanotolerance and trypanosusceptibility.

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.

Transcriptional adaptations following exercise in thoroughbred horse skeletal muscle highlights molecular mechanisms that lead to muscle hypertrophy.

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 (T(0)), immediately after (T(1)) and four hours after (T(2)) exercise. Statistically significant differences in mRNA abundance between time points (T(0) vs T(1) and T(0) vs T(2)) 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 T(1) (P < 0.05), by T(2) 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 (T(2)) 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 T(1), 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.

Transmission ratio distortion at the growth hormone gene (GH1) in bovine preimplantation embryos: An in vitro culture-induced phenomenon?

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.

Linkage mapping of the locus for inherited ovine arthrogryposis (IOA) to sheep chromosome 5.

Arthrogryposis is a congenital malformation affecting the limbs of newborn animals and infants. Previous work has demonstrated that inherited ovine arthrogryposis (IOA) has an autosomal recessive mode of inheritance. Two affected homozygous recessive (art/art) Suffolk rams were used as founders for a backcross pedigree of half-sib families segregating the IOA trait. A genome scan was performed using 187 microsatellite genetic markers and all backcross animals were phenotyped at birth for the presence and severity of arthrogryposis. Pairwise LOD scores of 1.86, 1.35, and 1.32 were detected for three microsatellites, BM741, JAZ, and RM006, that are located on sheep Chr 5 (OAR5). Additional markers in the region were identified from the genetic linkage map of BTA7 and by in silico analyses of the draft bovine genome sequence, three of which were informative. Interval mapping of all autosomes produced an F value of 21.97 (p < 0.01) for a causative locus in the region of OAR5 previously flagged by pairwise linkage analysis. Inspection of the orthologous region of HSA5 highlighted a previously fine-mapped locus for human arthrogryposis multiplex congenita neurogenic type (AMCN). A survey of the HSA5 genome sequence identified plausible candidate genes for both IOA and human AMCN.

Bioinformatic discovery of novel bioactive peptides.

Short synthetic oligopeptides based on regions of human proteins that encompass functional motifs are versatile reagents for understanding protein signaling and interactions. They can either mimic or inhibit the parent protein's activity and have been used in drug development. Peptide studies typically either derive peptides from a single identified protein or (at the other extreme) screen random combinatorial peptides, often without knowledge of the signaling pathways targeted. Our objective was to determine whether rational bioinformatic design of oligopeptides specifically targeted to potentially signaling-rich juxtamembrane regions could identify modulators of human platelet function. High-throughput in vitro platelet function assays of palmitylated cell-permeable oligopeptides corresponding to these regions identified many agonists and antagonists of platelet function. Many bioactive peptides were from adhesion molecules, including a specific CD226-derived inhibitor of inside-out platelet signaling. Systematic screens of this nature are highly efficient tools for discovering short signaling motifs in molecular signaling pathways.

Cytokine mRNA profiling of peripheral blood mononuclear cells from trypanotolerant and trypanosusceptible cattle infected with Trypanosoma congolense.

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 (T(H)2)-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 T(H)1 type cytokines than the trypanosusceptible Boran and may explain why N'Dama control parasitemia more efficiently than Boran during the early stages of infection.

Gene expression profiling of peripheral blood mononuclear cells (PBMC) from Mycobacterium bovis infected cattle after in vitro antigenic stimulation with purified protein derivative of tuberculin (PPD).

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 (P<0.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 (P<0.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.