Search

Infection with the zoonotic trematode Fasciola hepatica, common in many regions with a temperate climate, leads to delayed growth and loss of productivity in cattle, while infection in sheep can have more severe effects, potentially leading to death. Previous transcriptomic analyses revealed upregulation of TGFB1, cell death and Toll-like receptor signalling, T-cell activation, and inhibition of nitric oxide production in macrophages in response to infection. However, the differences between ovine and bovine responses have not yet been explored. The objective of this study was to further investigate the transcriptomic response of ovine peripheral blood mononuclear cells (PBMC) to F. hepatica infection, and to elucidate the differences between ovine and bovine PBMC responses.Sixteen male Merino sheep were randomly assigned to infected or control groups (n = 8 per group) and orally infected with 120 F. hepatica metacercariae. Transcriptomic data was generated from PBMC at 0, 2 and 16 weeks post-infection (wpi), and analysed for differentially expressed (DE) genes between infected and control animals at each time point (analysis 1), and for each group relative to time 0 (analysis 2). Analysis 2 was then compared to a similar study performed previously on bovine PBMC.A total of 453 DE genes were found at 2 wpi, and 2 DE genes at 16 wpi (FDR < 0.1, analysis 1). Significantly overrepresented biological pathways at 2 wpi included role of PKR in interferon induction and anti-viral response, death receptor signalling and RIG-I-like receptor signalling, which suggested that an activation of innate response to intracellular nucleic acids and inhibition of cellular apoptosis were taking place. Comparison of analysis 2 with the previous bovine transcriptomic study revealed that anti-inflammatory response pathways which were significantly overrepresented in the acute phase in cattle, including IL-10 signalling, Th2 pathway, and Th1 and Th2 activation were upregulated only in the chronic phase in sheep. We propose that the earlier activation of anti-inflammatory responses in cattle, as compared with sheep, may be related to the general absence of acute clinical signs in cattle. These findings offer scope for “smart vaccination” strategies for this important livestock parasite.

SUMMARYThe phenotypic diversity of African cattle reflects adaptation to a wide range of agroecological conditions, human-mediated selection preferences, and complex patterns of admixture between the humpless Bos taurus (taurine) and humped Bos indicus (zebu) subspecies, which diverged 150-500 thousand years ago. Despite extensive admixture, all African cattle possess taurine mitochondrial haplotypes, even populations with significant zebu biparental and male uniparental nuclear ancestry. This has been interpreted as the result of a human-mediated dispersal ultimately stemming from zebu bulls imported from South Asia during the last three millennia. Here we assess whether ancestry at mitochondrially-targeted nuclear genes in African admixed cattle is impacted by mitonuclear functional interactions. Using high-density SNP data, we find evidence for mitonuclear coevolution across hybrid African cattle populations with a significant increase of taurine ancestry at mitochondrially-targeted nuclear genes. Our results, therefore, support the hypothesis of incompatibility between the taurine mitochondrial genome and the zebu nuclear genome.GRAPHICAL SUMMARYHighlights•Using high-density genome-wide SNP data, we present evidence for mitonuclear coevolution in hybrid African cattle.•We observe a significant increase of taurine ancestry across multiple hybrid populations at mitochondrially-targeted nuclear genes.•Our results provide support for the hypothesis of mitonuclear incompatibility between the zebu nuclear genome and the taurine mitochondrial genome.

The Mongolian horse represents one of the most ancient extant horse populations. In this study we determined the male-specific region of the Y chromosome (MSY) haplotype distribution in 60 Chinese Mongolian horses representing five distinct populations. Cosmopolitan male lineages were predominant in horses from one improved (Sanhe), one Chinese Mongolian subtype (Baicha Iron Hoof) and one indigenous (Abaga Black) population. In contrast, autochthonous Y chromosome diversity was evident among the two landrace populations (Wushen and Wuzhumuqin), as the majority of their MSY haplotypes were situated at root nodes in a network. Our results also suggest gene flow between Chinese Mongolian and Arabian horses, as an appreciable number of Wuzhumuqin horses carried haplotypes that are typically observed in Arabian horses. Although most horses carried modern haplotypes as a direct result of recent breed improvement, authentic Chinese Mongolian horses retain an ancient signature of paternal lineages that has not previously been described in extant horse populations. Therefore, further characterization of MSY variation in these populations will be important for the discovery of lost diversity in modern domestic horses and also for understanding the evolutionary history of equine paternal lineages.

Bovine tuberculosis, caused by infection with members of the Mycobacterium tuberculosis complex, particularly 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 RNA sequencing, 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 analysed the transcriptome of bovine whole peripheral blood samples collected at −1 week pre-infection and +1, +2, +6, +10, and +12 weeks post-infection time points. Differentially expressed genes were catalogued and evaluated at each post-infection time point relative to the −1 week pre-infection time point and used for the identification of putative candidate host transcriptional biomarkers for M. bovis infection. Differentially expressed gene sets were also used for examination of cellular pathways associated with the host response to M. bovis infection, construction of de novo gene interaction networks enriched for host differentially expressed genes, and time-series analyses to identify functionally important groups of genes displaying similar patterns of expression across the infection time course. A notable outcome of these analyses was identification of a 19-gene transcriptional biosignature of infection consisting of genes increased in expression across the time course from +1 week to +12 weeks post-infection.

Fasciola hepatica is a trematode parasite responsible for major economic losses in livestock production, and is also a food-borne zoonotic agent in developing rural regions. For years, the immunoregulatory mechanisms employed by the parasite have hampered efforts to develop a successful vaccine candidate. Given that a comprehensive understanding of the immune response to infection is needed, we investigated the gene expression changes in ovine hepatic lymph nodes after experimental infection with F. hepatica. Lymph nodes from uninfected and infected animals were processed for RNA sequencing (RNA-seq) at 16 weeks post-infection. Comparison of groups revealed 5,132 differentially-expressed genes (DEGs). An inhibition of pro-inflammatory pathways, which has previously been described during fasciolosis, was evident in our data. However, other signals previously identified in ruminant peripheral blood mononuclear cells (PBMC) or liver tissue, such as activation of TGF-β or apoptosis-related pathways were not detected. We found inhibition of some key immunological pathways, including natural killer (NK) cell activity and IgE-mediated signaling. These may point to additional some as yet unrecognized mechanisms employed by the parasite to evade the host immune response. Understanding these, and leveraging information from this and other omics studies, will be important for the development of future vaccine prototypes against this parasite.

Background Bovine TB (bTB), caused by infection with Mycobacterium bovis, is a major endemic disease affecting global cattle production. The key innate immune cell that first encounters the pathogen is the alveolar macrophage, previously shown to be substantially reprogrammed during intracellular infection by the pathogen. Here we use differential expression, and correlation- and interaction-based network approaches to analyse the host response to infection with M. bovis at the transcriptome level to identify core infection response pathways and gene modules. These outputs were then integrated with genome-wide association study (GWAS) data sets to enhance detection of genomic variants for susceptibility/resistance to M. bovis infection. Results The host gene expression data consisted of RNA-seq data from bovine alveolar macrophages (bAM) infected with M. bovis at 24 and 48 h post-infection (hpi) compared to non-infected control bAM. These RNA-seq data were analysed using three distinct computational pipelines to produce six separate gene sets: 1) DE genes filtered using stringent fold-change and P-value thresholds (DEG-24: 378 genes, DEG-48: 390 genes); 2) genes obtained from expression correlation networks (CON-24: 460 genes, CON-48: 416 genes); and 3) genes obtained from differential expression networks (DEN-24: 339 genes, DEN-48: 495 genes). These six gene sets were integrated with three bTB breed GWAS data sets by employing a new genomics data integration tool—gwinteR. Using GWAS summary statistics, this methodology enabled detection of 36, 102 and 921 prioritised SNPs for Charolais, Limousin and Holstein-Friesian, respectively. Conclusions The results from the three parallel analyses showed that the three computational approaches could identify genes significantly enriched for SNPs associated with susceptibility/resistance to M. bovis infection. Results indicate distinct and significant overlap in SNP discovery, demonstrating that network-based integration of biologically relevant transcriptomics data can leverage substantial additional information from GWAS data sets. These analyses also demonstrated significant differences among breeds, with the Holstein-Friesian breed GWAS proving most useful for prioritising SNPS through data integration. Because the functional genomics data were generated using bAM from this population, this suggests that the genomic architecture of bTB resilience traits may be more breed-specific than previously assumed.

BackgroundMastitis in dairy cows caused byStaphylococcus aureusis a major problem hindering economic growth in dairy farms worldwide. It is difficult to prevent or eliminate due to its asymptomatic nature and long persistence of infection. Although transcriptomic responses of bovine mammary gland cells to pathogens that cause mastitis have been studied, the common responses of peripheral blood leukocytes toS. aureusinfection across two consecutive generations of dairy cattle have not been investigated.MethodsIn the current study, RNA-Seq was used to profile the transcriptomes of peripheral blood leukocytes sampled fromS. aureus-infected mothers and theirS. aureus-infected daughters, and also healthy non-infected mothers and their healthy daughters. Differential gene expression was evaluated as follows: 1)S. aureus-infected cows versus healthy non-infected cows (S vs. H, which include all the mothers and daughters), 2)S. aureus-infected mothers versus healthy non-infected mothers (SM vs. HM), and 3)S. aureus-infected daughters versus healthy non-infected daughters (SMD vs. HMD).ResultsAnalysis of all identified expressed genes in the four groups (SM, SMD, HM, and HMD) showed thatEPOR,IL9,IFNL3,CCL26,IL26were exclusively expressed in both the HM and HMD groups, and that they were significantly (P Padj. 1.2) were detected in the three comparisons, respectively. DEGs withP 2 were used for functional enrichment analyses. For the S vs. H comparison, DEGs detected includedCCL20,IL13andMMP3, which are associated with the IL-17 signaling pathway. In the SM vs. HM and SMD vs. HMD comparisons, five (BLA-DQB,C1R,C2,FCGR1A, andKRT10) and six (BLA-DQB,C3AR1,CFI,FCAR,FCGR3A, andLOC10498484) genes, respectively, were involved in theS. aureusinfection pathway.ConclusionsOur study provides insights into the transcriptomic responses of bovine peripheral blood leukocytes across two generations of cattle naturally infected withS. aureus. The genes highlighted in this study could serve as expression biomarkers for mastitis and may also contain sequence variation that can be used for genetic improvement of dairy cattle for resilience to mastitis.

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

Adaptation to early training and racing (i.e. precocity), which is highly variable in racing Thoroughbreds, has implications for the selection and training of horses. We hypothesised that precocity in Thoroughbred racehorses is heritable. Age at first sprint training session (work day), age at first race and age at best race were used as phenotypes to quantify precocity. Using high-density SNP array data, additive SNP heritability (hSNP2) was estimated to be 0.17, 0.14 and 0.17 for the three traits respectively. In genome-wide association studies (GWAS) for age at first race and age at best race, a 1.98-Mb region on equine chromosome 18 (ECA18) was identified. The most significant association was with the myostatin (MSTN) g.66493737C>T SNP (P = 5.46 × 10-12 and P = 1.89 × 10-14 respectively). In addition, two SNPs on ECA1 (g.37770220G>A and g.37770305T>C) within the first intron of the serotonin receptor gene HTR7 were significantly associated with age at first race and age at best race. Although no significant associations were identified for age at first work day, the MSTN:g.66493737C>T SNP was among the top 20 SNPs in the GWAS (P = 3.98 × 10-5 ). Here we have identified variants with potential roles in early adaptation to training. Although there was an overlap in genes associated with precocity and distance aptitude (i.e. MSTN), the HTR7 variants were more strongly associated with precocity than with distance. Because HTR7 is closely related to the HTR1A gene, previously implicated in tractability in young Thoroughbreds, this suggests that behavioural traits may influence precocity.

Background A single bout of exercise induces changes in gene expression in skeletal muscle. Regular exercise results in an adaptive response involving changes in muscle architecture and biochemistry, and is an effective way to manage and prevent common human diseases such as obesity, cardiovascular disorders and type II diabetes. However, the biomolecular mechanisms underlying such responses still need to be fully elucidated. Here we performed a transcriptome-wide analysis of skeletal muscle tissue in a large cohort of untrained Thoroughbred horses (n = 51) before and after a bout of high-intensity exercise and again after an extended period of training. We hypothesized that regular high-intensity exercise training primes the transcriptome for the demands of high-intensity exercise. Results An extensive set of genes was observed to be significantly differentially regulated in response to a single bout of high-intensity exercise in the untrained cohort (3241 genes) and following multiple bouts of high-intensity exercise training over a six-month period (3405 genes). Approximately one-third of these genes (1025) and several biological processes related to energy metabolism were common to both the exercise and training responses. We then developed a novel network-based computational analysis pipeline to test the hypothesis that these transcriptional changes also influence the contextual molecular interactome and its dynamics in response to exercise and training. The contextual network analysis identified several important hub genes, including the autophagosomal-related gene GABARAPL1, and dynamic functional modules, including those enriched for mitochondrial respiratory chain complexes I and V, that were differentially regulated and had their putative interactions ‘re-wired’ in the exercise and/or training responses. Conclusion Here we have generated for the first time, a comprehensive set of genes that are differentially expressed in Thoroughbred skeletal muscle in response to both exercise and training. These data indicate that consecutive bouts of high-intensity exercise result in a priming of the skeletal muscle transcriptome for the demands of the next exercise bout. Furthermore, this may also lead to an extensive ‘re-wiring’ of the molecular interactome in both exercise and training and include key genes and functional modules related to autophagy and the mitochondrion. Electronic supplementary material The online version of this article (doi:10.1186/s12864-017-4007-9) contains supplementary material, which is available to authorized users.

During the last decade, ancient DNA research has been revolutionized by the availability of increasingly powerful DNA sequencing and ancillary genomics technologies, giving rise to the new field of paleogenomics. In this review, we show how our understanding of the genetic basis of animal domestication and the origins and dispersal of livestock and companion animals during the Upper Paleolithic and Neolithic periods is being rapidly transformed through new scientific knowledge generated with paleogenomic methods. These techniques have been particularly informative in revealing high-resolution patterns of artificial and natural selection and evidence for significant admixture between early domestic animal populations and their wild congeners.

Animal domestication has fascinated biologists since Charles Darwin first drew the parallel between evolution via natural selection and human-mediated breeding of livestock and companion animals. In this review we show how studies of ancient DNA from domestic animals and their wild progenitors and congeners have shed new light on the genetic origins of domesticates, and on the process of domestication itself. High-resolution paleogenomic data sets now provide unprecedented opportunities to explore the development of animal agriculture across the world. In addition, functional population genomics studies of domestic and wild animals can deliver comparative information useful for understanding recent human evolution.

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

The Mongolian horse is one of the oldest extant horse populations and although domesticated, most animals are free-ranging and experience minimal human intervention. As an ancient population originating in one of the key domestication centers, the Mongolian horse may play a key role in understanding the origins and recent evolutionary history of horses. Here we describe an analysis of high-density genome-wide single-nucleotide polymorphism (SNP) data in 40 globally dispersed horse populations (n = 895). In particular, we have focused on new results from Chinese Mongolian horses (n = 100) that represent 5 distinct populations. These animals were genotyped for 670K SNPs and the data were analyzed in conjunction with 35K SNP data for 35 distinct breeds. Analyses of these integrated SNP data sets demonstrated that the Chinese Mongolian populations were genetically distinct from other modern horse populations. In addition, compared to other domestic horse breeds, the Chinese Mongolian horse populations exhibited relatively high genomic diversity. These results suggest that, in genetic terms, extant Chinese Mongolian horses may be the most similar modern populations to the animals originally domesticated in this region of Asia. Chinese Mongolian horse populations may therefore retain ancestral genetic variants from the earliest domesticates. Further genomic characterization of these populations in conjunction with archaeogenetic sequence data should be prioritized for understanding recent horse evolution and the domestication process that has led to the wealth of diversity observed in modern global horse breeds.

The Thoroughbred horse is a highly valued domestic animal population under strong selection for athletic phenotypes. Here we present a high resolution genomics-based analysis of inbreeding in the population that may form the basis for evidence-based discussion amid concerns in the breeding industry over the increasing use of small numbers of popular sire lines, which may accelerate a loss of genetic diversity. In the most comprehensive globally representative sample of Thoroughbreds to-date (n = 10,118), including prominent stallions (n = 305) from the major bloodstock regions of the world, we show using pan-genomic SNP genotypes that there has been a highly significant decline in global genetic diversity during the last five decades (FISR2 = 0.942, P = 2.19 × 10−13; FROHR2 = 0.88, P = 1.81 × 10−10) that has likely been influenced by the use of popular sire lines. Estimates of effective population size in the global and regional populations indicate that there is some level of regional variation that may be exploited to improve global genetic diversity. Inbreeding is often a consequence of selection, which in managed animal populations tends to be driven by preferences for cultural, aesthetic or economically advantageous phenotypes. Using a composite selection signals approach, we show that centuries of selection for favourable athletic traits among Thoroughbreds acts on genes with functions in behaviour, musculoskeletal conformation and metabolism. As well as classical selective sweeps at core loci, polygenic adaptation for functional modalities in cardiovascular signalling, organismal growth and development, cellular stress and injury, metabolic pathways and neurotransmitters and other nervous system signalling has shaped the Thoroughbred athletic phenotype. Our results demonstrate that genomics-based approaches to identify genetic outcrosses will add valuable objectivity to augment traditional methods of stallion selection and that genomics-based methods will be beneficial to actively monitor the population to address the marked inbreeding trend.

The Thoroughbred horse is a highly valued domestic animal population under strong selection for athletic phenotypes. Here we present a high resolution genomics-based analysis of inbreeding in the population that may form the basis for evidence-based discussion amid concerns in the breeding industry over the increasing use of small numbers of popular sire lines, which may accelerate a loss of genetic diversity. In the most comprehensive globally representative sample of Thoroughbreds to-date (n= 10,118), including prominent stallions (n= 305) from the major bloodstock regions of the world, we show using pan-genomic SNP genotypes that there has been a highly significant decline in global genetic diversity during the last five decades (FISR2= 0.942,P= 2.19 × 10−13;FROHR2= 0.88,P= 1.81 × 10−10) that has likely been influenced by the use of popular sire lines. Estimates of effective population size in the global and regional populations indicate that there is some level of regional variation that may be exploited to improve global genetic diversity. Inbreeding is often a consequence of selection, which in managed animal populations tends to be driven by preferences for cultural, aesthetic or economically advantageous phenotypes. Using a composite selection signals approach, we show that centuries of selection for favourable athletic traits among Thoroughbreds acts on genes with functions in behaviour, musculoskeletal conformation and metabolism. As well as classical selective sweeps at core loci, polygenic adaptation for functional modalities in cardiovascular signalling, organismal growth and development, cellular stress and injury, metabolic pathways and neurotransmitters and other nervous system signalling has shaped the Thoroughbred athletic phenotype. Our results demonstrate that genomics-based approaches to identify genetic outcrosses will add valuable objectivity to augment traditional methods of stallion selection and that genomics-based methods will be beneficial to actively monitor the population to address the marked inbreeding trend.Author SummaryIn the highly valuable global Thoroughbred horse industry, there is no systematic industry-mediated genetic population management. Purposeful inbreeding is common practice and there is an increasing use of popular sires. Inbreeding can lead to population health and fertility decline, but there is little objective genomics-based data for the Thoroughbred to catalyse action and support changes in breeding practices. Here, we describe the most comprehensive genetic analysis in the population among 10,000 Thoroughbreds from the major bloodstock regions of the world and reveal a highly significant increase in inbreeding during the last five decades. The main drivers of genetic diversity are the most influential ‘breed-shaping’ sire lines,Sadler’s Wells, DanehillandA.P. Indy. We identified genomic regions subject to positive selection containing genes for athletic traits. Our results highlight the need for population-wide efforts to proactively avert the potential increase of deleterious alleles that may impact on animal health in order to safeguard the future of a breed that is admired for its athleticism and enjoyed for sport worldwide.

The parasitic helminth Fasciola hepatica (liver fluke) causes economic loss to the livestock industry globally and also causes zoonotic disease. New control strategies such as vaccines are urgently needed, due to the rise of drug resistance in parasite populations. Vaccine development requires a comprehensive understanding of the immunological events during infection. Previous in vivo studies by our group have investigated global differentially expressed genes (DEGs) in ovine peripheral blood mononuclear cells (PBMC) in response to both acute and chronic F. hepatica infection. This work demonstrated that pathways involved in the pathogenesis of ovine fasciolosis included fibrosis, inhibition of macrophage nitric oxide production, and antibody isotype switching, among others. Transcriptomic changes in PBMC populations following F. hepatica infection in cattle, in which the disease phenotype is quite different, have not yet been examined. Using RNA sequencing we investigated gene expression changes in PBMC isolated from 9 non-infected and 11 F. hepatica-experimentally-infected calves immediately before infection, at 1 and at 14 weeks post-infection. Longitudinal time-course comparisons between groups revealed 21 and 1,624 DEGs driven exclusively by F. hepatica infection in cattle at acute and chronic stages, respectively. These results show that fewer DEGs at the acute stage of infection can be identified in cattle, as compared with sheep. In addition, the log2 fold-changes of these DEGs were relatively low (−1 to 3) reflecting the different clinical presentation of F. hepatica infection in cattle. Gene pathways for hepatic fibrosis and hepatic cholestasis along with apoptosis of antigen-presenting cells were enriched at chronic stages. Our results reflect the major differences in the disease phenotype between cattle and sheep and may indicate pathways to target in vaccine development.

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

Despite strong selection for athletic traits in Thoroughbred horses, there is marked variation in speed and aptitude for racing performance within the breed. Using global positioning system monitoring during exercise training, we measured speed variables and temporal changes in speed with age to derive phenotypes for GWAS. The aim of the study was to test the hypothesis that genetic variation contributes to variation in end‐point physiological traits, in this case galloping speed measured during field exercise tests. Standardisation of field‐measured phenotypes was attempted by assessing horses exercised on the same gallop track and managed under similar conditions by a single trainer. PCA of six key speed indices captured 73.9% of the variation with principal component 1 (PC1). Verifying the utility of the phenotype, we observed that PC1 (median) in 2‐year‐old horses was significantly different among elite, non‐elite and unraced horses (P < 0.001) and the temporal change with age in PC1 varied among horses with different myostatin (MSTN) g.66493737C>T SNP genotypes. A GWAS for PC1 in 2‐year‐old horses (n = 122) identified four SNPs reaching the suggestive threshold for association (P < 4.80 × 10−5), defining a 1.09 Mb candidate region on ECA8 containing the myosin XVIIIB (MYO18B) gene. In a GWAS for temporal change in PC1 with age (n = 168), five SNPs reached the suggestive threshold for association and defined candidate regions on ECA2 and ECA11. Both regions contained genes that are significantly differentially expressed in equine skeletal muscle in response to acute exercise and training stimuli, including MYO18A. As MYO18A plays a regulatory role in the skeletal muscle response to exercise, the identified genomic variation proximal to the myosin family genes may be important for the regulation of the response to exercise and training.

Durability traits in Thoroughbred horses are heritable, economically valuable and may affect horse welfare. The aims of this study were to test the hypotheses that (i) durability traits are heritable and (ii) genetic data may be used to predict a horse's potential to have a racecourse start. Heritability for the phenotype ‘number of 2‐ and 3‐year‐old starts’ was estimated to be urn:x-wiley:02689146:media:age12798:age12798-math-0001 = 0.11 ± 0.02 (n = 4499). A genome‐wide association study identified SNP contributions to the trait. The neurotrimin (NTM), opioid‐binding protein/cell adhesion molecule like (OPCML) and prolylcarboxypeptidase (PRCP) genes were identified as candidate genes associated with the trait. NTM functions in brain development and has been shown to have been selected during the domestication of the horse. PRCP is an established expression quantitative trait locus involved in the interaction between voluntary exercise and body composition in mice. We hypothesise that variation at these loci contributes to the motivation of the horse to exercise, which may influence its response to the demands of the training and racing environment. A random forest with mixed effects (RFME) model identified a set of SNPs that contributed to 24.7% of the heritable variation in the trait. In an independent validation set (n = 528 horses), the cohort with high genetic potential for a racecourse start had significantly fewer unraced horses (16% unraced) than did low (27% unraced) potential horses and had more favourable race outcomes among those that raced. Therefore, the information from SNPs included in the model may be used to predict horses with a greater chance of a racecourse start.

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

Cattle were domesticated ∼10,000 years ago, but analysis of modern breeds has not elucidated their origins. Verdugo et al. performed genome-wide analysis of 67 ancient Near Eastern Bos taurus DNA samples. Several populations of ancient aurochs were progenitors of domestic cows. These genetic lineages mixed ∼4000 years ago in a region around the Indus Valley. Interestingly, mitochondrial analysis indicated that genetic material likely derived from arid-adapted Bos indicus (zebu) bulls was introduced by introgression.Science, this issue p. 173Genome-wide analysis of 67 ancient Near Eastern cattle, Bos taurus, remains reveals regional variation that has since been obscured by admixture in modern populations. Comparisons of genomes of early domestic cattle to their aurochs progenitors identify diverse origins with separate introgressions of wild stock. A later region-wide Bronze Age shift indicates rapid and widespread introgression of zebu, Bos indicus, from the Indus Valley. This process was likely stimulated at the onset of the current geological age, ~4.2 thousand years ago, by a widespread multicentury drought. In contrast to genome-wide admixture, mitochondrial DNA stasis supports that this introgression was male-driven, suggesting that selection of arid-adapted zebu bulls enhanced herd survival. This human-mediated migration of zebu-derived genetics has continued through millennia, altering tropical herding on each continent.

SUMMARYThe Galway sheep population is the only native Irish sheep breed and represents an important livestock genetic resource, which 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 breed and suggest more recent gene flow from the Suffolk sheep breed. 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 meat and wool production.

Mycobacterial infections still represent a major burden on human and animal health, particularly in the context of human and bovine tuberculosis caused by infection with Mycobacterium tuberculosis and M. bovis, respectively. It is now well understood that differentiation of immune cells into their various lineages is influenced by microRNA transcriptional regulatory networks, and that dynamic and orchestrated control of miRNA and mRNA expression profiles is critical in mounting an effective host immune response to infection. In this chapter, we review studies that demonstrate how mycobacterial infections alter host miRNA expression profiles to evade immune responses and promote pathogen survival. We also discuss the potential diagnostic and therapeutic applications of miRNAs for mycobacterial infections.

Background: Race distance aptitude in Thoroughbred horses is highly heritable and is influenced largely by variation at the myostatin gene (MSTN). Objectives: In addition to MSTN, we hypothesised that other modifying loci contribute to best race distance. Study design: Using 3006 Thoroughbreds, including 835 ‘elite’ horses, which were >3 years old, had race records and were sampled from Europe/ Middle-East, Australia/New Zealand, North America and South Africa, we performed genome-wide association (GWA) tests and separately developed a genomic prediction algorithm to comprehensively catalogue additive genetic variation contributing to best race distance. Methods: 48,896 single‐nucleotide polymorphism (SNP) genotypes were generated from high‐density SNP genotyping arrays. Heritability estimates, tests of GWA and genomic prediction models were derived for the phenotypes: average race distance, best race distance for elite, nonelite and all winning horses. Results: Heritability estimates were high (urn:x-wiley:04251644:media:evj13058:evj13058-math-0001 = 0.51, best race distance – elite; urn:x-wiley:04251644:media:evj13058:evj13058-math-0002 = 0.42, best race distance – nonelite; urn:x-wiley:04251644:media:evj13058:evj13058-math-0003 = 0.40, best race distance – all) and most of the variation was attributed to the MSTN gene. MSTN locus SNPs were the most strongly associated with the trait and included BIEC2‐438999 (ECA18:66913090; P = 4.51 × 10−110, average race distance; P = 2.33 × 10−42, best race distance – elite). The genomic prediction algorithm enabled the inclusion of variation from all SNPs in a model that partitioned horses into short and long cohorts following assignment of MSTN genotype. Additional genes with minor contributions to best race distance were identified. Main limitations: The nongenetic influence of owner/trainer decisions on placement of horses in suitable races could not be controlled. Conclusions: MSTN is the single most important genetic contributor to best race distance in the Thoroughbred. Employment of genetic prediction models will lead to more accurate placing of horses in races that are best suited to their inherited genetic potential for distance aptitude.

Epigenetic modifications, such as DNA methylation, are pivotal in orchestrating gene expression patterns in various mammalian biological processes. They represent compelling mechanisms for conveying cellular response to external environmental stimuli more durably than transient transcriptional bursts/regulation events. Perturbation of the bovine alveolar macrophage (bAM) transcriptome due to Mycobacterium bovis infection on transcriptomic patterns has been well documented; however, the impact of this intracellular pathogen on the bAM epigenome has not been determined. For the current investigation, whole genome bisulfite sequencing (WGBS) was used to assess the effect of M. bovis infection on the bAM DNA methylome (5- methylcytosine). The methylomes of bAM infected with M. bovis (n = 8) were compared to those of non-infected control bAM (n = 8) at 24 hours post-infection (hpi). No differences in DNA methylation (CpG or non-CpG) were observed between control and infected bAM. Analysis of DNA methylation at proximal promoter regions uncovered >250 genes harbouring intermediately methylated (IM) promoters (average methylation = 33-66%). Gene ontology (GO) analysis, focusing on genes with low, intermediate or highly methylated promoters, revealed that genes with IM promoters were enriched for immune-related GO categories; this enrichment was not observed for genes in the high or low methylation groups. Targeted analysis of two non-imprinted genes in the IM category, C1QB and IL2RA, confirmed the WGBS observation. This study is the first in cattle to examine genome-wide DNA methylation at single nucleotide resolution in an important bovine cellular host-pathogen interaction model and provides evidence for intermediate promoter methylation in bAM.

Members of the Mycobacterium tuberculosis complex (MTBC) are the causative agents of tuberculosis in a range of mammals, including humans. A key feature of MTBC pathogens is their high degree of genetic identity, yet distinct host tropism. Notably, while Mycobacterium bovis is highly virulent and pathogenic for cattle, the human pathogen M. tuberculosis is attenuated in cattle. Previous research also suggests that host preference amongst MTBC members has a basis in host innate immune responses. To explore MTBC host tropism, we present in-depth profiling of the MTBC reference strains M. bovis AF2122/97 and M. tuberculosis H37Rv at both the global transcriptional and translational level via RNA-sequencing and SWATH mass spectrometry. Furthermore, a bovine alveolar macrophage infection time course model was employed to investigate the shared and divergent host transcriptomic response to infection with M. tuberculosis or M. bovis. Significant differential expression of virulence-associated pathways between the two bacilli was revealed, including the ESX-1 secretion system. A divergent transcriptional response was observed between M. tuberculosis and M. bovis infection of bovine alveolar macrophages, in particular cytosolic DNA-sensing pathways at 48 hours post-infection, and highlights a distinct engagement of M. bovis with the bovine innate immune system. The work presented here therefore provides a basis for the identification of host innate immune mechanisms subverted by virulent host-adapted mycobacteria to promote their survival during the early stages of infection.ImportanceThe Mycobacterium tuberculosis complex (MTBC) includes the most important global pathogens for humans and animals, namely Mycobacterium tuberculosis and Mycobacterium bovis, respectively. These two exemplar mycobacterial pathogens share a high degree of genetic identity, but the molecular basis for their distinct host preference is unknown. In this work we integrated transcriptomic and proteomic analyses of the pathogens to elucidate global quantitative differences between them at the mRNA and protein level. We then integrated this data with transcriptome analysis of the bovine macrophage response to infection with either pathogen. Increased expression of the ESX-1 virulence system in M. bovis appeared a key driver of an increased cytosolic nucleic acid sensing and interferon response in bovine macrophages infected with M. bovis compared to M. tuberculosis. Our work demonstrates the specificity of host-pathogen interaction and how the subtle interplay between mycobacterial phenotype and host response may underpin host specificity amongst MTBC members.

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.

Kerry cattle are an endangered landrace heritage breed of cultural importance to Ireland. In the present study we have used genome-wide SNP data (Illumina®BovineSNP50 array) to evaluate genomic diversity within the Kerry cattle population and between Kerry cattle and other European cattle breeds. Visualisation of patterns of genetic differentiation and gene flow among cattle breeds using phylogenetic trees with ancestry graphs highlighted, in particular, 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 diminishingNeover time and that recent estimatedNevalues 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.Note 1:This is an Associate Editor (D.E.M) Inaugural Article submission toFrontiers in Genetics: Livestock GenomicsNote 2:British English language style preferred for publication of this article.

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.

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

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 (ORF) DNA sequence of 19,994 orthologous protein-coding gene pairs from extant Bos taurus genomes and a single extinct B. primigenius genome. Using this data from 19,994 gene pairs, we identified miRNA binding site polymorphisms in the 3' UTRs of 1,620 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 (SSr) 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.

Under selection pressure from pathogens, variable NK cell receptors that recognize polymorphic MHC class I evolved convergently in different species of placental mammal. Unexpectedly, diversified killer cell Ig–like receptors (KIRs) are shared by simian primates, including humans, and cattle, but not by other species. Whereas much is known of human KIR genetics and genomics, knowledge of cattle KIR is limited to nine cDNA sequences. To facilitate comparison of the cattle and human KIR gene families, we determined the genomic location, structure, and sequence of two cattle KIR haplotypes and defined KIR sequences of aurochs, the extinct wild ancestor of domestic cattle. Larger than its human counterpart, the cattle KIR locus evolved through successive duplications of a block containing ancestral KIR3DL and KIR3DX genes that existed before placental mammals. Comparison of two cattle KIR haplotypes and aurochs KIR show the KIR are polymorphic and the gene organization and content appear conserved. Of 18 genes, 8 are functional and 10 were inactivated by point mutation. Selective inactivation of KIR3DL and activating receptor genes leaves a functional cohort of one inhibitory KIR3DL, one activating KIR3DX, and six inhibitory KIR3DX. Functional KIR diversity evolved from KIR3DX in cattle and from KIR3DL in simian primates. Although independently evolved, cattle and human KIR gene families share important function-related properties, indicating that cattle KIR are NK cell receptors for cattle MHC class I. Combinations of KIR and MHC class I are the major genetic factors associated with human disease and merit investigation in cattle.

The analysis of mitochondrial and nuclear DNA sequence polymorphisms from modern cattle populations has had a profound impact on our understanding of the events surrounding the domestication of cattle. From these studies, it has been possible to distinguish between pre- and post-domestic genetic differentiation, supporting previous assertions from archaeological studies and, in some cases, revealing novel aspects of the demographic history of cattle.\ud Analyses of genetic material retrieved from the remains of extinct ancestral wild cattle have also added valuable layers of information pertaining to cattle domestic origins; however, information from these investigations have, in general, been limited to small, variable portions of the mitochondrial genome owing to technical challenges associated with the retrieval and amplification of ancient DNA.\ud In recent years, however, new high-throughput, massively parallel genomics technology platforms, such as single-nucleotide polymorphism (SNP) genotyping arrays and next-generation sequencing (NGS), have provided a new impetus to the studies of genetic variation in extant and ancient cattle.\ud Arrays of SNP have facilitated high-resolution genetic surveys of global cattle populations and detection of ancient and recent genomic selective sweeps. Next-generation sequencing analyses of modern and ancient cattle hold great promise for identifying and cataloging of pre- and post-domestication patterns of genomic variation and correlating this with natural and artificial selection processes.

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

Mycobacteria are a genus of bacteria that range from the non-pathogenic Mycobacterium smegmatis to Mycobacterium tuberculosis, the causative agent of tuberculosis in humans. Mycobacteria primarily infect host tissues through inhalation or ingestion. They are phagocytosed by host macrophages and dendritic cells. Here, conserved pathogen-associated molecular patterns (PAMPs) on the surface of mycobacteria are recognized by phagocytic pattern recognition receptors (PRRs). Several families of PRRs have been shown to non-opsonically recognize mycobacterial PAMPs, including membrane-bound C-type lectin receptors, membrane-bound and cytosolic Toll-like receptors and cytosolic NOD-like receptors. Recently, a possible role for intracellular cytosolic PRRs in the recognition of mycobacterial pathogens has been proposed. Here, we discuss currentideas on receptor-mediated recognition of mycobacterial pathogens by macrophages and dendritic cells.

Variation at the myostatin (MSTN) gene locus has been shown to influence racing phenotypes in Thoroughbred horses, and in particular, early skeletal muscle development and the aptitude for racing at short distances. Specifically, a single nucleotide polymorphism (SNP) in the first intron of MSTN (g.66493737C/T) is highly predictive of best race distance among Flat racing Thoroughbreds: homozygous C/C horses are best suited to short distance races, heterozygous C/T horses are best suited to middle distance races, and homozygous T/T horses are best suited to longer distance races. Patent applications for this gene marker association, and other linked markers, have been filed. The information contained within the patent applications is exclusively licensed to the commercial biotechnology company Equinome Ltd, which provides a DNA-based test to the international Thoroughbred horse racing and breeding industry. The application of this information in the industry enables informed decision making in breeding and racing and can be used to assist selection to accelerate the rate of change of genetic types among distinct populations (Case Study 1) and within individual breeding operations (Case Study 2).

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

Sequence variation at the equine myostatin gene ( MSTN) locus has previously been shown to have a singular genomic influence on optimum race distance in Thoroughbred racehorses. Myostatin, encoded by the MSTN gene, is a member of the TGF-β superfamily that regulates skeletal muscle development in a range of mammalian species including the horse. In the Thoroughbred, the C-allele at the g.66493737C/T SNP has been found at significantly higher frequency in subgroups of the population that are suited to fast, short distance, sprint races and also influences body composition phenotypes. We investigated the influence of the g.66493737C/T SNP on speed indexes measured in a cohort of n = 85 Thoroughbred horses-in-training. We found significant associations between genotypes at the g.66493737C/T SNP and all measured speed variables: Dist6 [distance travelled during 6 s before and after maximal velocity (Vmax); P = 0.0040], Vmaxt (duration at Vmax; P = 0.0249), Vmax ( P = 0.0265), Dist6b (distance travelled during 6 s before Vmax; P = 0.0032), and Dist6a(distance travelled during 6 s after Vmax; P = 0.0317). For each measure, horses with the C/C and C/T genotypes outperformed T/T horses, indicating the requirement for at least one C-allele to improve speed. For the most significantly associated variables (Dist6 and Dist6b) the C/C cohort performed better than the T/T cohort with the heterozygotes intermediate, indicating a dose-dependent manifestation. These findings clearly indicate that variation at the MSTN gene influences speed in Thoroughbred horses.

Genomic imprinting is an epigenetic non-Mendelian phenomenon found predominantly in placental mammals. Imprinted genes display differential expression in the offspring depending on whether the gene is maternally or paternally inherited. Currently, some 100 imprinted genes have been reported in mammals, and while some of these genes are imprinted across most mammalian species, others have been shown to be imprinted in only a few species. The PHLDA2 gene that codes for a pleckstrin homology-like domain, family A (member 2), protein has to date been shown to be a maternally expressed imprinted gene in humans, mice and pigs. Genes subject to imprinting can have major effects on mammalian growth, development and disease. For instance, disruption of imprinted genes can lead to aberrant growth syndromes in cloned domestic mammals, and it has been demonstrated that PHLDA2 mRNA expression levels are aberrant in the placenta of somatic clones of cattle. In this study, we demonstrate that PHLDA2 is expressed across a range of cattle foetal tissues and stages and provide the first evidence that PHLDA2 is a monoallelically expressed imprinted gene in cattle foetal tissues, and also in the bovine placenta.

Summary The regulation of the bioavailability of insulin-like growth factors (IGFs) is critical for normal mammalian growth and development. The imprinted insulin-like growth factor 2 receptor gene (IGF2R) encodes a transmembrane protein receptor that acts to sequester and degrade excess circulating insulin-like growth factor 2 (IGF-II) – a potent foetal mitogen – and is considered an important inhibitor of growth. Consequently, IGF2R may serve as a candidate gene underlying important growth- and body-related quantitative traits in domestic mammalian livestock. In this study, we have quantified genotype–phenotype associations between three previously validated intronic bovine IGF2R single nucleotide polymorphisms (SNPs) (IGF2R:g.64614T>C, IGF2R:g.65037T>C and IGF2R:g. 86262C>T) and a range of performance traits in 848 progeny-tested Irish Holstein-Friesian artificial insemination sires. Notably, all three polymorphisms analysed were associated (P £ 0.05) with at least one of a number of performance traits related to animal body size: angularity, body depth, chest width, rump width, and animal stature. In addition, the C-toT transition at the IGF2R:g.65037T>C polymorphism was positively associated with cow carcass weight and angularity. Correction for multiple testing resulted in the retention of two genotype–phenotype associations (animal stature and rump width). None of the SNPs analysed were associated with any of the milk traits examined. Analysis of pairwise r 2 measures of linkage disequilibrium between all three assayed SNPs ranged between 0.41 and 0.79, suggesting that some of the observed SNP associations with performance may be independent. To our knowledge, this is one of the first studies demonstrating associations between IGF2R polymorphisms and growth- and body-related traits in cattle. These results also support the increasing body of evidence that imprinted genes harbour polymorphisms that contribute to heritable variation in phenotypic traits in domestic livestock species.

Summary Variations in the growth hormone receptor (GHR) gene sequence are associated with performance traits in cattle. For example, the single nucleotide polymorphism (SNP) F279Y in transmembrane exon 8 has a strong association with milk yield. In this study, 32 previously unreported, putative novel SNPs (31 in the 5′ non-coding region) were identified by resequencing ∼19 kb of the GHR gene in genomic DNA from 22 cattle of multiple breeds. A population of 848 Holstein–Friesian AI sires was subsequently genotyped for the 32 putative novel SNPs and seven published SNPs (including F279Y, one in exon 1A promoter and five in exon 10). Associations between each segregating SNP and genetic merit for performance were quantified in the 848 Holstein–Friesians using weighted animal linear mixed models. Six of the published SNPs and seven of the novel SNPs were associated with at least one of the traits – milk yield, fat yield, protein yield, fat percentage, protein percentage, somatic cell score, calving interval, survival and growth and size traits. Even when the allelic substitution effect (P

Athletic performance is influenced by a complex interplay among the environment and a suite of genes, which contributes to system-wide structure and function. In a panel of elite and non-elite Thoroughbred horses (n=148), we genotyped 68 SNPs in 17 putative exercise-relevant genes chosen from a genome scan for selection. We performed a series of case-control and quantitative association tests for relationships with racecourse performance. Thirteen SNPs in nine genes were significantly (P

Summary The recent completion of the horse genome and commercial availability of an equine SNP genotyping array has facilitated the mapping of disease genes. We report putative localization of the gene responsible for dwarfism, a trait in Friesian horses that is thought to have a recessive mode of inheritance, to a 2- MB region of chromosome 14 using just 10 affected animals and 10 controls. We successfully genotyped 34 429 SNPs that were tested for association with dwarfism using chi-square tests. The most significant SNP in our study, BIEC2-239376 (P2df = 4.54 · 10 )5 , Prec = 7.74 · 10 )6 ), is located close to a gene implicated in human dwarfism. Fine-mapping and resequencing analyses did not aid in further localization of the causative variant, and replication of our findings in independent sample sets will be necessary to confirm these results.

Summary Reasons for performing study: The role of molecular signalling pathways in the phenotypic adaptation of skeletal muscle to different exercise stimuli in the Thoroughbred horse has not been reported previously. Objective: To examine CKM, COX4I1, COX4I2 and PDK4 gene expression following high intensity sprint and moderate intensity treadmill exercise stimuli in skeletal muscle of Thoroughbred horses. Materials and methods: Two groups of trained 3-year-old Thoroughbred horses participated. Group A (n = 6 females, n = 3 males) participated in an incremental step test (moderate intensity) to fatigue or HRmax on a Sato high speed treadmill (distance = 5418.67m 343.21). Group B (n = 8 females) participated in routine 'work' (sprint) on an all- weather gallop (distance = 812.83 m 12.53). Biopsy samples were obtained from the gluteus medius pre-exercise (T0), immediately post exercise (T1) and 4 h post exercise (T2). For physiological relevance venous blood samples were collected to measure plasma lactate and creatine kinase concentrations. Changes in mRNA expression were determined by quantitative real-time RT-PCR for creatine kinase muscle (CKM), cytochrome c oxidase subunit IV isoform 1 (COX4I1), cytochrome c oxidase subunit IV isoform 2 (COX4I2) and pyruvate dehydrogenase kinase, isozyme 4 (PDK4) genes. Statistical significance (a

Summary 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 T was confirmed (P

Single nucleotide polymorphisms (SNPs) represent the most common form of DNA sequence variation in mammalian livestock genomes. While the past decade has witnessed major advances in SNP genotyping technologies, genotyping errors caused, in part, by the biochemistry underlying the genotyping platform used, can occur. These errors can distort project results and conclusions and can result in incorrect decisions in animal management and breeding programs; hence, SNP genotype calls must be accurate and reliable. In this study, 263 Bos spp. samples were genotyped commercially for a total of 16 SNPs. Of the total possible 4,208 SNP genotypes, 4,179 SNP genotypes were generated, yielding a genotype call rate of 99.31% (standard deviation ± 0.93%). Between 110 and 263 samples were subsequently re-genotyped by us for all 16 markers using a custom-designed SNP genotyping platform, and of the possible 3,819 genotypes a total of 3,768 genotypes were generated (98.70% genotype call rate, SD ± 1.89%). A total of 3,744 duplicate genotypes were generated for both genotyping platforms, and comparison of the genotype calls for both methods revealed 3,741 concordant SNP genotype call rates (99.92% SNP genotype concordance rate). These data indicate that both genotyping methods used can provide livestock geneticists with reliable, reproducible SNP genotypic data for in-depth statistical analysis.

Advances in high-throughput geno- typing technologies have afforded researchers the op- portunity to study ever-increasing numbers of SNP in animal genomes. However, many studies encounter difficulties in obtaining sufficient quantities of high- quality DNA for such analyses, particularly when the source biological material is limited or degraded. The recent development of in vitro whole-genome amplifica- tion approaches has permitted researchers to circum- vent these challenges by increasing the amount of us- able DNA in normally small-quantity samples. Here, we assess the performance of whole-genome amplification products generated from ovine genomic DNA using a high-throughput SNP genotyping platform, the newly developed Illumina ovineSNP50 BeadChip. Our results demonstrate a high genotype call rate for conventional genomic DNA and whole-genome amplified genomic DNA. The data also reveal an exceptionally high con- cordance rate (≥99%) between the genotypes generated from whole-genome amplified products and their con- ventional genomic DNA counterparts. This study sup- ports the use of whole-genome amplification as a viable solution for the analysis of high-density SNP genotypic data using compromised or limited starting material.