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

peer-reviewedBackground: Recovering high quality intact RNA from post-mortem tissue is of major concern for gene expression studies in animals and humans. Since the availability of post-mortem tissue is often associated with substantial delay, it is important that we understand the temporal variation in the stability of total RNA and of individual gene transcripts so as to be able to appropriately interpret the data generated from such studies. Hence, the objective of this experiment was to qualitatively and quantitatively assess the integrity of total and messenger RNA extracted from bovine skeletal muscle, subcutaneous adipose tissue and liver stored at 4°C at a range of time points up to 22 days post-mortem. These conditions were designed to mimic the environment prevailing during the transport of beef from the abattoir to retail outlets. Results: The 28S and 18S rRNA molecules of total RNA were intact for up to 24 h post-mortem in liver and adipose tissues and up to 8 days post-mortem in skeletal muscle. The mRNA of housekeeping genes (GAPDH and ACTB) and two diet-related genes (RBP5 and SCD) were detectable up to 22 days post-mortem in skeletal muscle. While the mRNA stability of the two housekeeping genes was different in skeletal muscle and liver, they were similar to each other in adipose tissue. After 22 days post-mortem, the relative abundance of RBP5 gene was increased in skeletal muscle and in adipose tissue and decreased in liver. During this period, the relative abundance of SCD gene also increased in skeletal muscle whereas it decreased in both adipose tissue and liver. Conclusion: Stability of RNA in three tissues (skeletal muscle, subcutaneous adipose tissue and liver) subjected to long-term post-mortem storage at refrigeration temperature indicated that skeletal muscle can be a suitable tissue for recovering biologically useful RNA for gene expression studies even if the tissue is subjected to post-mortem storage for weeks, whereas adipose tissue and liver should be processed within 24 hours post-mortem.This research was part funded by the Irish Government under the National Development Plan 2000-2006. Support to BB under the Walsh Fellowship Programme of Teagasc is acknowledged

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

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

Polymorphism discovery and allele frequency estimation using high-throughput DNA sequencing of target-enriched pooled DNA samples

BACKGROUND: The central role of the somatotrophic axis in animal post-natal growth, development and fertility is well established. Therefore, the identification of genetic variants affecting quantitative traits within this axis is an attractive goal. However, large sample numbers are a pre-requisite for the identification of genetic variants underlying complex traits and although technologies are improving rapidly, high-throughput sequencing of large numbers of complete individual genomes remains prohibitively expensive. Therefore using a pooled DNA approach coupled with target enrichment and high-throughput sequencing, the aim of this study was to identify polymorphisms and estimate allele frequency differences across 83 candidate genes of the somatotrophic axis, in 150 Holstein-Friesian dairy bulls divided into two groups divergent for genetic merit for fertility. RESULTS: In total, 4,135 SNPs and 893 indels were identified during the resequencing of the 83 candidate genes. Nineteen percent (n = 952) of variants were located within 5' and 3' UTRs. Seventy-two percent (n = 3,612) were intronic and 9% (n = 464) were exonic, including 65 indels and 236 SNPs resulting in non-synonymous substitutions (NSS). Significant (P < 0.01) mean allele frequency differentials between the low and high fertility groups were observed for 720 SNPs (58 NSS). Allele frequencies for 43 of the SNPs were also determined by genotyping the 150 individual animals (Sequenom(® )MassARRAY). No significant differences (P > 0.1) were observed between the two methods for any of the 43 SNPs across both pools (i.e., 86 tests in total). CONCLUSIONS: The results of the current study support previous findings of the use of DNA sample pooling and high-throughput sequencing as a viable strategy for polymorphism discovery and allele frequency estimation. Using this approach we have characterised the genetic variation within genes of the somatotrophic axis and related pathways, central to mammalian post-natal growth and development and subsequent lactogenesis and fertility. We have identified a large number of variants segregating at significantly different frequencies between cattle groups divergent for calving interval plausibly harbouring causative variants contributing to heritable variation. To our knowledge, this is the first report describing sequencing of targeted genomic regions in any livestock species using groups with divergent phenotypes for an economically important trait

Tracing Cattle Breeds with Principal Components Analysis Ancestry Informative SNPs

The recent release of the Bovine HapMap dataset represents the most detailed survey of bovine genetic diversity to date, providing an important resource for the design and development of livestock production. We studied this dataset, comprising more than 30,000 Single Nucleotide Polymorphisms (SNPs) for 19 breeds (13 taurine, three zebu, and three hybrid breeds), seeking to identify small panels of genetic markers that can be used to trace the breed of unknown cattle samples. Taking advantage of the power of Principal Components Analysis and algorithms that we have recently described for the selection of Ancestry Informative Markers from genomewide datasets, we present a decision-tree which can be used to accurately infer the origin of individual cattle. In doing so, we present a thorough examination of population genetic structure in modern bovine breeds. Performing extensive cross-validation experiments, we demonstrate that 250-500 carefully selected SNPs suffice in order to achieve close to 100% prediction accuracy of individual ancestry, when this particular set of 19 breeds is considered. Our methods, coupled with the dense genotypic data that is becoming increasingly available, have the potential to become a valuable tool and have considerable impact in worldwide livestock production. They can be used to inform the design of studies of the genetic basis of economically important traits in cattle, as well as breeding programs and efforts to conserve biodiversity. Furthermore, the SNPs that we have identified can provide a reliable solution for the traceability of breed-specific branded products

MHC Class I Bound to an Immunodominant Theileria parva Epitope Demonstrates Unconventional Presentation to T Cell Receptors

T cell receptor (TCR) recognition of peptide-MHC class I (pMHC) complexes is a crucial event in the adaptive immune response to pathogens. Peptide epitopes often display a strong dominance hierarchy, resulting in focusing of the response on a limited number of the most dominant epitopes. Such T cell responses may be additionally restricted by particular MHC alleles in preference to others. We have studied this poorly understood phenomenon using Theileria parva, a protozoan parasite that causes an often fatal lymphoproliferative disease in cattle. Despite its antigenic complexity, CD8+ T cell responses induced by infection with the parasite show profound immunodominance, as exemplified by the Tp1214–224 epitope presented by the common and functionally important MHC class I allele N*01301. We present a high-resolution crystal structure of this pMHC complex, demonstrating that the peptide is presented in a distinctive raised conformation. Functional studies using CD8+ T cell clones show that this impacts significantly on TCR recognition. The unconventional structure is generated by a hydrophobic ridge within the MHC peptide binding groove, found in a set of cattle MHC alleles. Extremely rare in all other species, this feature is seen in a small group of mouse MHC class I molecules. The data generated in this analysis contribute to our understanding of the structural basis for T cell-dependent immune responses, providing insight into what determines a highly immunogenic p-MHC complex, and hence can be of value in prediction of antigenic epitopes and vaccine design

Trans-Species Polymorphism and Selection in the MHC Class II DRA Genes of Domestic Sheep

Highly polymorphic genes with central roles in lymphocyte mediated immune surveillance are grouped together in the major histocompatibility complex (MHC) in higher vertebrates. Generally, across vertebrate species the class II MHC DRA gene is highly conserved with only limited allelic variation. Here however, we provide evidence of trans-species polymorphism at the DRA locus in domestic sheep (Ovis aries). We describe variation at the Ovar-DRA locus that is far in excess of anything described in other vertebrate species. The divergent DRA allele (Ovar-DRA*0201) differs from the sheep reference sequences by 20 nucleotides, 12 of which appear non-synonymous. Furthermore, DRA*0201 is paired with an equally divergent DRB1 allele (Ovar-DRB1*0901), which is consistent with an independent evolutionary history for the DR sub-region within this MHC haplotype. No recombination was observed between the divergent DRA and B genes in a range of breeds and typical levels of MHC class II DR protein expression were detected at the surface of leukocyte populations obtained from animals homozygous for the DRA*0201, DRB1*0901 haplotype. Bayesian phylogenetic analysis groups Ovar-DRA*0201 with DRA sequences derived from species within the Oryx and Alcelaphus genera rather than clustering with other ovine and caprine DRA alleles. Tests for Darwinian selection identified 10 positively selected sites on the branch leading to Ovar-DRA*0201, three of which are predicted to be associated with the binding of peptide antigen. As the Ovis, Oryx and Alcelaphus genera have not shared a common ancestor for over 30 million years, the DRA*0201 and DRB1*0901 allelic pair is likely to be of ancient origin and present in the founding population from which all contemporary domestic sheep breeds are derived. The conservation of the integrity of this unusual DR allelic pair suggests some selective advantage which is likely to be associated with the presentation of pathogen antigen to T-cells and the induction of protective immunity

A genome-wide SNP-association study confirms a sequence variant (g.66493737C>T) in the equine myostatin (MSTN) gene as the most powerful predictor of optimum racing distance for Thoroughbred racehorses

<p>Abstract</p> <p>Background</p> <p>Thoroughbred horses have been selected for traits contributing to speed and stamina for centuries. It is widely recognized that inherited variation in physical and physiological characteristics is responsible for variation in individual aptitude for race distance, and that muscle phenotypes in particular are important.</p> <p>Results</p> <p>A genome-wide SNP-association study for optimum racing distance was performed using the EquineSNP50 Bead Chip genotyping array in a cohort of <it>n </it>= 118 elite Thoroughbred racehorses divergent for race distance aptitude. In a cohort-based association test we evaluated genotypic variation at 40,977 SNPs between horses suited to short distance (≤ 8 f) and middle-long distance (> 8 f) races. The most significant SNP was located on chromosome 18: BIEC2-417495 ~690 kb from the gene encoding myostatin (<it>MSTN</it>) [<it>P</it>_unadj.= 6.96 × 10^-6]. Considering best race distance as a quantitative phenotype, a peak of association on chromosome 18 (chr18:65809482-67545806) comprising eight SNPs encompassing a 1.7 Mb region was observed. Again, similar to the cohort-based analysis, the most significant SNP was BIEC2-417495 (<it>P</it>_unadj.= 1.61 × 10^-9; <it>P</it>_Bonf.= 6.58 × 10^-5). In a candidate gene study we have previously reported a SNP (g.66493737C>T) in <it>MSTN </it>associated with best race distance in Thoroughbreds; however, its functional and genome-wide relevance were uncertain. Additional re-sequencing in the flanking regions of the <it>MSTN </it>gene revealed four novel 3' UTR SNPs and a 227 bp SINE insertion polymorphism in the 5' UTR promoter sequence. Linkage disequilibrium was highest between g.66493737C>T and BIEC2-417495 (<it>r</it>²= 0.86).</p> <p>Conclusions</p> <p>Comparative association tests consistently demonstrated the g.66493737C>T SNP as the superior variant in the prediction of distance aptitude in racehorses (g.66493737C>T, <it>P </it>= 1.02 × 10^-10; BIEC2-417495, <it>P</it>_unadj.= 1.61 × 10^-9). Functional investigations will be required to determine whether this polymorphism affects putative transcription-factor binding and gives rise to variation in gene and protein expression. Nonetheless, this study demonstrates that the g.66493737C>T SNP provides the most powerful genetic marker for prediction of race distance aptitude in Thoroughbreds.</p

Whole mitochondrial genomes unveil the impact of domestication on goat matrilineal variability

Background: The current extensive use of the domestic goat (Capra hircus) is the result of its medium size and high adaptability as multiple breeds. The extent to which its genetic variability was influenced by early domestication practices is largely unknown. A common standard by which to analyze maternally-inherited variability of livestock species is through complete sequencing of the entire mitogenome (mitochondrial DNA, mtDNA). Results: We present the first extensive survey of goat mitogenomic variability based on 84 complete sequences selected from an initial collection of 758 samples that represent 60 different breeds of C. hircus, as well as its wild sister species, bezoar (Capra aegagrus) from Iran. Our phylogenetic analyses dated the most recent common ancestor of C. hircus to ~460,000 years (ka) ago and identified five distinctive domestic haplogroups (A, B1, C1a, D1 and G). More than 90 % of goats examined were in haplogroup A. These domestic lineages are predominantly nested within C. aegagrus branches, diverged concomitantly at the interface between the Epipaleolithic and early Neolithic periods, and underwent a dramatic expansion starting from ~12–10 ka ago. Conclusions: Domestic goat mitogenomes descended from a small number of founding haplotypes that underwent domestication after surviving the last glacial maximum in the Near Eastern refuges. All modern haplotypes A probably descended from a single (or at most a few closely related) female C. aegagrus. Zooarchaelogical data indicate that domestication first occurred in Southeastern Anatolia. Goats accompanying the first Neolithic migration waves into the Mediterranean were already characterized by two ancestral A and C variants. The ancient separation of the C branch (~130 ka ago) suggests a genetically distinct population that could have been involved in a second event of domestication. The novel diagnostic mutational motifs defined here, which distinguish wild and domestic haplogroups, could be used to understand phylogenetic relationships among modern breeds and ancient remains and to evaluate whether selection differentially affected mitochondrial genome variants during the development of economically important breeds

Large-Scale Mitochondrial DNA Analysis of the Domestic Goat Reveals Six Haplogroups with High Diversity

Background. From the beginning of domestication, the transportation of domestic animals resulted in genetic and demographic processes that explain their present distribution and genetic structure. Thus studying the present genetic diversity helps to better understand the history of domestic species. Methodology/Principal Findings. The genetic diversity of domestic goats has been characterized with 2430 individuals from all over the old world, including 946 new individuals from regions poorly studied until now (mainly the Fertile Crescent). These individuals represented 1540 haplotypes for the HVI segment of the mitochondrial DNA (mtDNA) control region. This large-scale study allowed the establishment of a clear nomenclature of the goat maternal haplogroups. Only five of the six previously defined groups of haplotypes were divergent enough to be considered as different haplogroups. Moreover a new mitochondrial group has been localized around the Fertile Crescent. All groups showed very high haplotype diversity. Most of this diversity was distributed among groups and within geographic regions. The weak geographic structure may result from the worldwide distribution of the dominant A haplogroup (more than 90% of the individuals). The large-scale distribution of other haplogroups (except one), may be related to human migration. The recent fragmentation of local goat populations into discrete breeds is not detectable with mitochondrial markers. The estimation of demographic parameters from mismatch analyses showed that all groups had a recent demographic expansion corresponding roughly to the period when domestication took place. But even with a large data set it remains difficult to give relative dates of expansion for different haplogroups because of large confidence intervals. Conclusions/Significance. We propose standard criteria for the definition of the different haplogroups based on the result of mismatch analysis and on the use of sequences of reference. Such a method could be also applied for clarifying the nomenclature of mitochondrial haplogroups in other domestic species