In this study, mRNA-Seq was used to characterize and compare the leukocyte transcriptomes from two taurine breeds (Holstein and Jersey), and one indicine breed (Cholistani). At the genomic level, we identified breed-specific base changes in protein coding regions. Among 7,793,425 coding bases, only 165 differed between Holstein and Jersey, and 3,383 (0.04%) differed between Holstein and Cholistani, 817 (25%) of which resulted in amino acid changes in 627 genes. At the transcriptional level, we assembled transcripts and estimated their abundances including those from more than 3,000 unannotated intergeneic regions. Differential gene expression analysis showed a high similarity between Holstein and Jersey, and a much greater difference between the taurine breeds and the indicine breed. We identified gene ontology pathways that were systematically altered, including the electron transport chain and immune response pathways that may contribute to different levels of heat tolerance and disease resistance in taurine and indicine breeds. At the post-transcriptional level, sequencing mRNA allowed us to identify a number of genes undergoing differential alternative splicing among different breeds. This study provided a high-resolution survey of the variation between bovine transcriptomes at different levels and may provide important biological insights into the phenotypic differentiation among cattle breeds.
References
[1]
MacHugh DE, Shriver MD, Loftus RT, Cunningham P, Bradley DG (1997) Microsatellite DNA variation and the evolution, domestication and phylogeography of taurine and zebu cattle (Bos taurus and Bos indicus). Genetics 146: 1071–1086.
[2]
Fries R, Ruvinsky A (1999) The genetics of cattle. Fries R, Ruvinsky A, editors. New York, NY: CABI Publising.
[3]
Utech KB, Wharton RH, Kerr JD (1978) Resistance to Boophilus microplus (Canestrini) in different breeds of cattle. Australian Journal of Agricultural Research 29: 885–895.
[4]
Hansen PJ (2004) Physiological and cellular adaptations of zebu cattle to thermal stress. Anim Reprod Sci 82–83: 349–360.
[5]
Farooq U, Samad HA, Sher F, Asim M, Khan MA (2010) Cholistan and Cholistani breed of cattle. Pak Vet J 30: 126–130.
[6]
Elsik CG, Tellam RL, Worley KC, The Bovine Genome Sequencing and Analysis Consortium (2009) The genome sequence of taurine cattle: a window to ruminant biology and evolution. Science 324: 522–528.
[7]
The Bovine HapMap Consortium (2009) Genome-wide survey of SNP variation uncovers the genetic structure of cattle breeds. Science 324: 528–532.
[8]
Passador-Gurgel G, Hsieh WP, Hunt P, Deighton N, Gibson G (2007) Quantitative trait transcripts for nicotine resistance in Drosophila melanogaster. Nat Genet 39: 264–268.
[9]
Gilad Y, Oshlack A, Smyth GK, Speed TP, White KP (2006) Expression profiling in primates reveals a rapid evolution of human transcription factors. Nature 440: 242–245.
[10]
Oleksiak MF, Churchill GA, Crawford DL (2002) Variation in gene expression within and among natural populations. Nat Genet 32: 261–266.
[11]
Blekhman R, Marioni JC, Zumbo P, Stephens M, Gilad Y (2010) Sex-specific and lineage-specific alternative splicing in primates. Genome Res 20: 180–189.
[12]
Lin L, Shen S, Jiang P, Sato S, Davidson BL, et al. (2010) Evolution of alternative splicing in primate brain transcriptomes. Hum Mol Genet 19: 2958–2973.
[13]
Langmead B, Trapnell C, Pop M, Salzberg SL (2009) Ultrafast and memory-efficient alignment of short DNA sequences to the human genome. Genome Biol 10: R25.
Trapnell C, Williams BA, Pertea G, Mortazavi A, Kwan G, et al. (2010) Transcript assembly and quantification by RNA-Seq reveals unannotated transcripts and isoform switching during cell differentiation. Nat Biotechnol 28: 511–515.
[16]
Bullard JH, Purdom E, Hansen KD, Dudoit S (2010) Evaluation of statistical methods for normalization and differential expression in mRNA-Seq experiments. BMC Bioinformatics 11: 94.
[17]
Hansen KD, Brenner SE, Dudoit S (2010) Biases in Illumina transcriptome sequencing caused by random hexamer priming. Nucleic Acids Res 38: e131.
[18]
Roberts A, Trapnell C, Donaghey J, Rinn JL, Pachter L (2011) Improving RNA-Seq expression estimates by correcting for fragment bias. Genome Biol 12: R22.
[19]
Young MD, Wakefield MJ, Smyth GK, Oshlack A (2010) Gene ontology analysis for RNA-seq: accounting for selection bias. Genome Biol 11: R14.
[20]
Benjamini Y, Hochberg Y (1995) Controlling the False Discovery Rate: a Practical and Powerful Approach to Multiple Testing. J R Statist Soc B 57: 963–971.
[21]
Livak KJ, Schmittgen TD (2001) Analysis of relative gene expression data using real-time quantitative PCR and the 2(-Delta Delta C(T)) Method. Methods 25: 402–408.
[22]
Xia C, Ma W, Stafford LJ, Marcus S, Xiong WC, et al. (2001) Regulation of the p21-activated kinase (PAK) by a human Gbeta -like WD-repeat protein, hPIP1. Proc Natl Acad Sci U S A 98: 6174–6179.
[23]
Kumar R, Gururaj AE, Barnes CJ (2006) p21-activated kinases in cancer. Nat Rev Cancer 6: 459–471.
[24]
Piper EK, Jonsson NN, Gondro C, Lew-Tabor AE, Moolhuijzen P, et al. (2009) Immunological profiles of Bos taurus and Bos indicus cattle infested with the cattle tick, Rhipicephalus (Boophilus) microplus. Clin Vaccine Immunol 16: 1074–1086.
[25]
Hwang SO, Boswell SA, Seo JS, Lee SW (2008) Novel oxidative stress-responsive gene ERS25 functions as a regulator of the heat-shock and cell death response. J Biol Chem 283: 13063–13069.
[26]
Kamwanja LA, Chase CC Jr, Gutierrez JA, Guerriero V Jr, Olson TA, et al. (1994) Responses of bovine lymphocytes to heat shock as modified by breed and antioxidant status. J Anim Sci 72: 438–444.
[27]
Wang ET, Sandberg R, Luo S, Khrebtukova I, Zhang L, et al. (2008) Alternative isoform regulation in human tissue transcriptomes. Nature 456: 470–476.
