Stephen D. E. Park

Integration of Proteomics and Genomics in Platelets A PROFILE OF PLATELET PROTEINS AND PLATELET-SPECIFIC GENES

Platelets, while anucleate, contain RNA, some of which is translated into protein upon activation. Hypothesising that the platelet proteome is reflected in the transcriptome, we identified 82 proteins secreted from activated platelets and compared these, as well as published proteomic data, to the transcriptional profile. We also compared the transcriptome of platelets to other tissues to identify platelet-specific genes and used ontology to determine gene categories over-represented in platelets. RNA was isolated from highly pure platelet preparations for hybridization to Affymetrix oligonucleotide arrays. We identified 2,928 distinct messages as being present in platelets. The platelet transcriptome was compared with the proteome by relating both to UniGene clusters. Platelet proteomic data correlated well with the transcriptome, with 69% of secreted proteins detectable at the mRNA level, and similar concordance was obtained using two published datasets. While many of the most abundant mRNAs are for known platelet proteins, messages were detected for proteins not previously reported in platelets. Some of these may represent residual megakaryocyte messages; however, proteomic analysis confirmed the expression of many previously unreported genes in platelets. Transcripts for well-described platelet proteins are among the most platelet-specific messages. Ontological categories related to signal transduction, receptors, ion channels, and membranes are over-represented in platelets, while categories involved in protein synthesis are depleted. Despite the absence of gene transcription, the platelet proteome is mirrored in the transcriptome. Conversely, transcriptional analysis predicts the presence of novel proteins in the platelet. Transcriptional analysis is relevant to platelet biology, providing insights into platelet function and the mechanisms of platelet disorders.

A Genome Scan for Positive Selection in Thoroughbred Horses

Thoroughbred horses have been selected for exceptional racing performance resulting in system-wide structural and functional adaptations contributing to elite athletic phenotypes. Because selection has been recent and intense in a closed population that stems from a small number of founder animals Thoroughbreds represent a unique population within which to identify genomic contributions to exercise-related traits. Employing a population genetics-based hitchhiking mapping approach we performed a genome scan using 394 autosomal and X chromosome microsatellite loci and identified positively selected loci in the extreme tail-ends of the empirical distributions for (1) deviations from expected heterozygosity (Ewens-Watterson test) in Thoroughbred (n = 112) and (2) global differentiation among four geographically diverse horse populations (FST). We found positively selected genomic regions in Thoroughbred enriched for phosphoinositide-mediated signalling (3.2-fold enrichment; P<0.01), insulin receptor signalling (5.0-fold enrichment; P<0.01) and lipid transport (2.2-fold enrichment; P<0.05) genes. We found a significant overrepresentation of sarcoglycan complex (11.1-fold enrichment; P<0.05) and focal adhesion pathway (1.9-fold enrichment; P<0.01) genes highlighting the role for muscle strength and integrity in the Thoroughbred athletic phenotype. We report for the first time candidate athletic-performance genes within regions targeted by selection in Thoroughbred horses that are principally responsible for fatty acid oxidation, increased insulin sensitivity and muscle strength: ACSS1 (acyl-CoA synthetase short-chain family member 1), ACTA1 (actin, alpha 1, skeletal muscle), ACTN2 (actinin, alpha 2), ADHFE1 (alcohol dehydrogenase, iron containing, 1), MTFR1 (mitochondrial fission regulator 1), PDK4 (pyruvate dehydrogenase kinase, isozyme 4) and TNC (tenascin C). Understanding the genetic basis for exercise adaptation will be crucial for the identification of genes within the complex molecular networks underlying obesity and its consequential pathologies, such as type 2 diabetes. Therefore, we propose Thoroughbred as a novel in vivo large animal model for understanding molecular protection against metabolic disease.

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

BackgroundDomestication 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.ResultsPhylogenomic 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.ConclusionsThis 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.

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.

Association of sequence variants in CKM (creatine kinase, muscle) and COX4I2 (cytochrome c oxidase, subunit 4, isoform 2) genes with racing performance in Thoroughbred horses

REASONS FOR PERFORMING STUDY The wild progenitors of the domestic horse were subject to natural selection for speed and stamina for millennia. Uniquely, this process has been augmented in Thoroughbreds, which have undergone at least 3 centuries of intense artificial selection for athletic phenotypes. While the phenotypic adaptations to exercise are well described, only a small number of the underlying genetic variants contributing to these phenotypes have been reported. OBJECTIVES A panel of candidate performance-related genes was examined for DNA sequence variation in Thoroughbreds and the association with racecourse performance investigated. MATERIALS AND METHODS Eighteen candidate genes were chosen for their putative roles in exercise. Re-sequencing in Thoroughbred samples was successful for primer sets in 13 of these genes. SNPs identified in this study and from the EquCab2.0 SNP database were genotyped in 2 sets of Thoroughbred samples (n = 150 and 148) and a series of population-based case-control investigations were performed by separating the samples into discrete cohorts on the basis of retrospective racecourse performance. RESULTS Twenty novel SNPs were detected in 3 genes: ACTN3, CKM and COX4I2. Genotype frequency distributions for 3 SNPs in CKM and COX4I2 were significantly (P < 0.05) different between elite Thoroughbreds and racehorses that had never won a race. These associations were not validated when an additional (n = 130) independent set of samples was genotyped, but when analyses included all samples (n = 278) the significance of association at COX4I2 g.22684390C > T was confirmed (P < 0.02). CONCLUSIONS While molecular genetic information has the potential to become a powerful tool to make improved decisions in horse industries, it is vital that rigour is applied to studies generating these data and that adequate and appropriate sample sets, particularly for independent replication, are used.

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

BackgroundSelection 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.ResultsSkeletal 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 (T0vs T1 and T0vs 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.ConclusionsThis 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.

Alterations in oxidative gene expression in equine skeletal muscle following exercise and training

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 (VHR(max)) and peak postexercise plasma lactate concentration ([La]T(1)) revealed significant (P < 0.05) associations with postexercise COX4I1 and PPARCG1A expression and between [La]T(1) and 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 VHR(max) and [La]T(1) were 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.

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.

Pan-genomic analysis of bovine monocyte-derived macrophage gene expression in response to in vitro infection with Mycobacterium avium subspecies paratuberculosis

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.

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.