David A. Magee

Genetic evidence for Near-Eastern origins of European cattle

The limited ranges of the wild progenitors of many of the primary European domestic species point to their origins further east in Anatolia or the fertile crescent. The wild ox (Bos primigenius), however, ranged widely and it is unknown whether it was domesticated within Europe as one feature of a local contribution to the farming economy. Here we examine mitochondrial DNA control-region sequence variation from 392 extant animals sampled from Europe, Africa and the Near East, and compare this with data from four extinct British wild oxen. The ancient sequences cluster tightly in a phylogenetic analysis and are clearly distinct from modern cattle. Network analysis of modern Bos taurus identifies four star-like clusters of haplotypes, with intra-cluster diversities that approximate to that expected from the time depth of domestic history. Notably, one of these clusters predominates in Europe and is one of three encountered at substantial frequency in the Near East. In contrast, African diversity is almost exclusively composed of a separate haplogroup, which is encountered only rarely elsewhere. These data provide strong support for a derived Near-Eastern origin for European cattle.

Zebu Cattle Are an Exclusive Legacy of the South Asia Neolithic

Animal domestication was a major step forward in human prehistory, contributing to the emergence of more complex societies. At the time of the Neolithic transition, zebu cattle (Bos indicus) were probably the most abundant and important domestic livestock species in Southern Asia. Although archaeological evidence points toward the domestication of zebu cattle within the Indian subcontinent, the exact geographic origins and phylogenetic history of zebu cattle remains uncertain. Here, we report evidence from 844 zebu mitochondrial DNA (mtDNA) sequences surveyed from 19 Asiatic countries comprising 8 regional groups, which identify 2 distinct mitochondrial haplogroups, termed I1 and I2. The marked increase in nucleotide diversity (P < 0.001) for both the I1 and I2 haplogroups within the northern part of the Indian subcontinent is consistent with an origin for all domestic zebu in this area. For haplogroup I1, genetic diversity was highest within the Indus Valley among the three hypothesized domestication centers (Indus Valley, Ganges, and South India). These data support the Indus Valley as the most likely center of origin for the I1 haplogroup and a primary center of zebu domestication. However, for the I2 haplogroup, a complex pattern of diversity is detected, preventing the unambiguous pinpointing of the exact place of origin for this zebu maternal lineage. Our findings are discussed with respect to the archaeological record for zebu domestication within the Indian subcontinent.

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.

Polymorphisms in bovine immune genes and their associations with somatic cell count and milk production in dairy cattle

BackgroundMastitis, an inflammation of the mammary gland, is a major source of economic loss on dairy farms. The aim of this study was to quantify the associations between two previously identified polymorphisms in the bovine toll-like receptor 2 (TLR2) and chemokine receptor 1 (CXCR1) genes and mammary health indictor traits in (a) 246 lactating dairy cow contemporaries representing five breeds from one research farm and (b) 848 Holstein-Friesian bulls that represent a large proportion of the Irish dairy germplasm. To expand the study, a further 14 polymorphisms in immune genes were included for association studies in the bull population.ResultsTLR4-2021 associated (P < 0.05) with both milk protein and fat percentage in late lactation (P < 0.01) within the cow cohort. No association was observed between this polymorphism and either yield or composition of milk within the bull population. CXCR1-777 significantly associated (P < 0.05) with fat yield in the bull population and tended to associate (P < 0.1) with somatic cell score (SCS) in the cows genotyped. CD14-1908 A allele was found to associate with increased (P < 0.05) milk fat and protein yield and also tended to associate with increased (P < 0.1) milk yield. A SERPINA1 haplotype with superior genetic merit for milk protein yield and milk fat percentage (P < 0.05) was also identified.ConclusionOf the sixteen polymorphisms in seven immune genes genotyped, just CXCR1-777 tended to associate with SCS, albeit only in the on-farm study. The lack of an association between the polymorphisms with SCS in the Holstein-Friesian data set would question the potential importance of these variants in selection for improved mastitis resistance in the Holstein-Friesian cow.

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.

Single nucleotide polymorphisms at the imprinted bovine insulin-like growth factor 2 (IGF2) locus are associated with dairy performance in Irish Holstein-Friesian cattle

The imprinted insulin-like growth factor 2 gene (IGF2) encodes a fetal mitogenic hormone protein (IGF-II) and has previously been shown to be associated with performance in dairy cattle. In this study we assessed genotype-phenotype associations between four single nucleotide polymorphisms (SNPs) located within the bovine IGF2 locus on chromosome 29 and a range of performance traits related to milk production, animal growth and body size, fertility and progeny survival in 848 progeny-tested Irish Holstein-Friesian sires. Two of the four SNPs (rs42196909 and IGF2.g-3815A>G), which were in strong linkage disequilibrium (r2 = 0·995), were associated with milk yield (P ≤ 0·01) and milk protein yield (P ≤ 0·05); the rs42196901 SNP was also associated (P ≤ 0·05) with milk fat yield. Associations (P ≤ 0·05) with milk fat percentage and milk protein percentage were observed at the rs42196901 and IGF2.g-3815A>G SNPs, respectively. The rs42196909 and IGF2.g-3815A>G SNPs were also associated with progeny carcass conformation (P ≤ 0·05), while an association (P ≤ 0·01) with progeny carcass weight was observed at the rs42194733 SNP locus. None of the four SNPs were associated with body size, fertility and progeny survival. These findings support previous work which suggests that the IGF2 locus is an important biological regulator of milk production in dairy cattle and add to an accumulating body of research showing that imprinted genes influence many complex performance traits in cattle.

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

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

Transcriptional profiling of immune genes in bovine monocyte-derived macrophages exposed to bacterial antigens.

The involvement of Toll-like receptors (TLRs) and other immune signalling genes during challenge of bovine macrophages with bacterial products derived from disease-causing bacteria in cattle was investigated. An in vitro cell culture model of bovine monocyte derived macrophages (MDM) was established and these cells were exposed to purified protein derivative (PPD-b) derived from Mycobacterium bovis and to lipopolysaccharide (LPS) derived from Escherichia coli. Following 24h incubation, total RNA was extracted and expression of immune related genes was determined by real time quantitative reverse transcription PCR (qRT-PCR). Expression of a selection of genes spanning the TLR-2 and TLR-4 pathways, from the initial activation of the receptors to the production of pro-inflammatory cytokines and chemokines was determined. Results from repeat experiments using MDM from seven different age-matched dairy cattle showed that PPD-b treatment caused significant up-regulation of the TLR2 and TLR4 genes and the expression profile of TLR adaptor molecules suggested that this signalling is MYD88-dependent. Conversely, LPS caused significant up-regulation of TLR4 via a MYD88-independent signalling pathway. Significant up-regulation of genes involved with NF-κB signalling was also detected in PPD-b- and LPS-treated samples accompanied by the expression of pro-inflammatory cytokine (TNF, IL1B, IL6) and chemokine genes (IL8, CCL5, CCL3). Overall, LPS challenge resulted in a more marked up-regulation of immune-related genes. Furthermore, the magnitude fold-change difference in gene expression suggests, at least in part, that bovine macrophages produce IFN-γ as a result of LPS challenge.

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.