1
|
Chen N, Zhang Z, Hou J, Chen J, Gao X, Tang L, Wangdue S, Zhang X, Sinding MHS, Liu X, Han J, Lü H, Lei C, Marshall F, Liu X. Evidence for early domestic yak, taurine cattle, and their hybrids on the Tibetan Plateau. SCIENCE ADVANCES 2023; 9:eadi6857. [PMID: 38091398 PMCID: PMC10848728 DOI: 10.1126/sciadv.adi6857] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 05/16/2023] [Accepted: 11/07/2023] [Indexed: 12/18/2023]
Abstract
Domestic yak, cattle, and their hybrids are fundamental to herder survival at high altitudes on the Tibetan Plateau. However, little is known about their history. Bos remains are uncommon in this region, and ancient domestic yak have not been securely identified. To identify Bos taxa and investigate their initial management, we conducted zooarchaeological analyses of 193 Bos specimens and sequenced five nuclear genomes from recently excavated assemblages at Bangga. Morphological data indicated that more cattle than yak were present. Ancient mitochondrial DNA and nuclear genome sequences identified taurine cattle and provided evidence for domestic yak and yak-cattle hybridization ~2500 years ago. Reliance on diverse Bos species and their hybrid has increased cattle adaptation and herder resilience to plateau conditions. Ancient cattle and yak at Bangga were closely related to contemporary livestock, indicating early herder legacies and the continuity of cattle and yak husbandry on the Tibetan Plateau.
Collapse
Affiliation(s)
- Ningbo Chen
- Key Laboratory of Animal Genetics Breeding and Reproduction of Shaanxi Province, College of Animal Science and Technology, Northwest A&F University, Yangling 712100, P. R. China
| | - Zhengwei Zhang
- Center for Archaeological Science, Sichuan University, Chengdu 610065, P. R. China
| | - Jiawen Hou
- Key Laboratory of Animal Genetics Breeding and Reproduction of Shaanxi Province, College of Animal Science and Technology, Northwest A&F University, Yangling 712100, P. R. China
| | - Jialei Chen
- Key Laboratory of Animal Genetics Breeding and Reproduction of Shaanxi Province, College of Animal Science and Technology, Northwest A&F University, Yangling 712100, P. R. China
| | - Xuan Gao
- Center for Archaeological Science, Sichuan University, Chengdu 610065, P. R. China
| | - Li Tang
- Department of Archaeology, Max Planck Institute of Geoanthropology, Jena 07745, Germany
- Domestication and Anthropogenic Evolution Research Group, Max Planck Institute of Geoanthropology, Jena 07745, Germany
| | - Shargan Wangdue
- Institute for Conservation and Research of Cultural Relics of Tibet Autonomous Region, Lhasa 850000, China
| | - Xiaoming Zhang
- State Key Laboratory of Genetic Resources and Evolution, Kunming Institute of Zoology, Chinese Academy of Sciences (CAS), Kunming 650201, P. R. China
| | - Mikkel-Holger S. Sinding
- Section for Computational and RNA Biology, Department of Biology, University of Copenhagen, Copenhagen DK-1350, Denmark
| | - Xuexue Liu
- National Germplasm Centre of Domestic Animal Resources, Institute of Animal Science, Chinese Academy of Agricultural Sciences (CAAS), Beijing 100193, P. R. China
- Centre for Anthropobiology and Genomics of Toulouse, CNRS UMR 5288, Université de Toulouse, Université Paul Sabatier, Toulouse 31000, France
| | - Jianlin Han
- Yazhouwan National Laboratory, Sanya 572024, P. R. China
- CAAS-ILRI Joint Laboratory on Livestock and Forage Genetic Resources, Institute of Animal Science, Chinese Academy of Agricultural Sciences (CAAS), Beijing 100193, P. R. China
| | - Hongliang Lü
- Center for Archaeological Science, Sichuan University, Chengdu 610065, P. R. China
| | - Chuzhao Lei
- Key Laboratory of Animal Genetics Breeding and Reproduction of Shaanxi Province, College of Animal Science and Technology, Northwest A&F University, Yangling 712100, P. R. China
| | - Fiona Marshall
- Department of Anthropology, Washington University in St. Louis, St. Louis, MO 63130, USA
| | - Xinyi Liu
- Department of Anthropology, Washington University in St. Louis, St. Louis, MO 63130, USA
| |
Collapse
|
2
|
Liu X, Liu W, Lenstra JA, Zheng Z, Wu X, Yang J, Li B, Yang Y, Qiu Q, Liu H, Li K, Liang C, Guo X, Ma X, Abbott RJ, Kang M, Yan P, Liu J. Evolutionary origin of genomic structural variations in domestic yaks. Nat Commun 2023; 14:5617. [PMID: 37726270 PMCID: PMC10509194 DOI: 10.1038/s41467-023-41220-x] [Citation(s) in RCA: 13] [Impact Index Per Article: 13.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/12/2023] [Accepted: 08/23/2023] [Indexed: 09/21/2023] Open
Abstract
Yak has been subject to natural selection, human domestication and interspecific introgression during its evolution. However, genetic variants favored by each of these processes have not been distinguished previously. We constructed a graph-genome for 47 genomes of 7 cross-fertile bovine species. This allowed detection of 57,432 high-resolution structural variants (SVs) within and across the species, which were genotyped in 386 individuals. We distinguished the evolutionary origins of diverse SVs in domestic yaks by phylogenetic analyses. We further identified 334 genes overlapping with SVs in domestic yaks that bore potential signals of selection from wild yaks, plus an additional 686 genes introgressed from cattle. Nearly 90% of the domestic yaks were introgressed by cattle. Introgression of an SV spanning the KIT gene triggered the breeding of white domestic yaks. We validated a significant association of the selected stratified SVs with gene expression, which contributes to phenotypic variations. Our results highlight that SVs of different origins contribute to the phenotypic diversity of domestic yaks.
Collapse
Affiliation(s)
- Xinfeng Liu
- State Key Laboratory of Herbage Improvement and Grassland Agro-ecosystem, College of Ecology, Lanzhou University, Lanzhou, 730000, China
- Key Laboratory of Animal Genetics and Breeding on Tibetan Plateau, Ministry of Agriculture and Rural Affairs, Lanzhou Institute of Husbandry and Pharmaceutical Sciences, Chinese Academy of Agricultural Sciences, Lanzhou, 730050, China
- Academy of Plateau Science and Sustainability, Qinghai Normal University, Xining, 810016, China
| | - Wenyu Liu
- State Key Laboratory of Herbage Improvement and Grassland Agro-ecosystem, College of Ecology, Lanzhou University, Lanzhou, 730000, China
| | - Johannes A Lenstra
- Faculty of Veterinary Medicine, Utrecht University, Utrecht, 3508 TD, The Netherlands
| | - Zeyu Zheng
- State Key Laboratory of Herbage Improvement and Grassland Agro-ecosystem, College of Ecology, Lanzhou University, Lanzhou, 730000, China
| | - Xiaoyun Wu
- Key Laboratory of Animal Genetics and Breeding on Tibetan Plateau, Ministry of Agriculture and Rural Affairs, Lanzhou Institute of Husbandry and Pharmaceutical Sciences, Chinese Academy of Agricultural Sciences, Lanzhou, 730050, China
| | - Jiao Yang
- State Key Laboratory of Herbage Improvement and Grassland Agro-ecosystem, College of Ecology, Lanzhou University, Lanzhou, 730000, China
| | - Bowen Li
- State Key Laboratory of Herbage Improvement and Grassland Agro-ecosystem, College of Ecology, Lanzhou University, Lanzhou, 730000, China
| | - Yongzhi Yang
- State Key Laboratory of Herbage Improvement and Grassland Agro-ecosystem, College of Ecology, Lanzhou University, Lanzhou, 730000, China
| | - Qiang Qiu
- State Key Laboratory of Herbage Improvement and Grassland Agro-ecosystem, College of Ecology, Lanzhou University, Lanzhou, 730000, China
| | - Hongyu Liu
- Anhui Provincial Laboratory of Local Livestock and Poultry Genetical Resource Conservation and Breeding, College of Animal Science and Technology, Anhui Agricultural University, Hefei, 230036, China
| | - Kexin Li
- State Key Laboratory of Herbage Improvement and Grassland Agro-ecosystem, College of Ecology, Lanzhou University, Lanzhou, 730000, China
| | - Chunnian Liang
- Key Laboratory of Animal Genetics and Breeding on Tibetan Plateau, Ministry of Agriculture and Rural Affairs, Lanzhou Institute of Husbandry and Pharmaceutical Sciences, Chinese Academy of Agricultural Sciences, Lanzhou, 730050, China
| | - Xian Guo
- Key Laboratory of Animal Genetics and Breeding on Tibetan Plateau, Ministry of Agriculture and Rural Affairs, Lanzhou Institute of Husbandry and Pharmaceutical Sciences, Chinese Academy of Agricultural Sciences, Lanzhou, 730050, China
| | - Xiaoming Ma
- Key Laboratory of Animal Genetics and Breeding on Tibetan Plateau, Ministry of Agriculture and Rural Affairs, Lanzhou Institute of Husbandry and Pharmaceutical Sciences, Chinese Academy of Agricultural Sciences, Lanzhou, 730050, China
| | - Richard J Abbott
- School of Biology, University of St Andrews, St Andrews, KY16 9AJ, UK
| | - Minghui Kang
- State Key Laboratory of Herbage Improvement and Grassland Agro-ecosystem, College of Ecology, Lanzhou University, Lanzhou, 730000, China.
| | - Ping Yan
- Key Laboratory of Animal Genetics and Breeding on Tibetan Plateau, Ministry of Agriculture and Rural Affairs, Lanzhou Institute of Husbandry and Pharmaceutical Sciences, Chinese Academy of Agricultural Sciences, Lanzhou, 730050, China.
| | - Jianquan Liu
- State Key Laboratory of Herbage Improvement and Grassland Agro-ecosystem, College of Ecology, Lanzhou University, Lanzhou, 730000, China.
- Academy of Plateau Science and Sustainability, Qinghai Normal University, Xining, 810016, China.
| |
Collapse
|
3
|
Brunson K, Witt KE, Monge S, Williams S, Peede D, Odsuren D, Bukhchuluun D, Cameron A, Szpak P, Amartuvshin C, Honeychurch W, Wright J, Pleuger S, Erdene M, Tumen D, Rogers L, Khatanbaatar D, Batdalai B, Galdan G, Janz L. Ancient Mongolian aurochs genomes reveal sustained introgression and management in East Asia. BIORXIV : THE PREPRINT SERVER FOR BIOLOGY 2023:2023.08.10.552443. [PMID: 37609302 PMCID: PMC10441390 DOI: 10.1101/2023.08.10.552443] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 08/24/2023]
Abstract
Societies in East Asia have utilized domesticated cattle for over 5000 years, but the genetic history of cattle in East Asia remains understudied. Genome-wide analyses of 23 ancient Mongolian cattle reveal that East Asian aurochs and ancient East Asian taurine cattle are closely related, but neither are closely related to any modern East Asian breeds. We observe binary variation in aurochs diet throughout the early Neolithic, and genomic evidence shows millennia of sustained male-dominated introgression. We identify a unique connection between ancient Mongolian aurochs and the European Hereford breed. These results point to the likelihood of human management of aurochs in Northeast Asia prior to and during the initial adoption of taurine cattle pastoralism.
Collapse
Affiliation(s)
| | - Kelsey E. Witt
- Department of Genetics and Biochemistry and Center for Human Genetics, Clemson University; Clemson, South Carolina 29634, USA
- Center for Computational Molecular Biology, Brown University; Providence 02912, USA
- Department of Ecology, Evolution, and Organismal Biology, Brown University; Providence 02912, USA
| | - Susan Monge
- Department of Anthropology, University of Illinois Chicago, Chicago, IL 60607, USA
| | - Sloan Williams
- Department of Anthropology, University of Illinois Chicago, Chicago, IL 60607, USA
| | - David Peede
- Center for Computational Molecular Biology, Brown University; Providence 02912, USA
- Department of Ecology, Evolution, and Organismal Biology, Brown University; Providence 02912, USA
- Institute at Brown for Environment and Society, Brown University; Providence 02912, USA
| | - Davaakhuu Odsuren
- Department of History, Mongolian National University of Education; Ulaanbaatar, Sukhbaatar district, 210648, Mongolia
- Institute of Archaeology, Mongolian Academy of Sciences, Ulaanbaatar-51, Mongolia
| | - Dashzeveg Bukhchuluun
- Department of Anthropology, Yale University, 10 Sachem St., New Haven, CT 06511, USA
| | - Asa Cameron
- Department of Anthropology, Yale University, 10 Sachem St., New Haven, CT 06511, USA
| | - Paul Szpak
- Department of Anthropology, Trent University; Peterborough K9J 6Y1, Canada
| | - Chunag Amartuvshin
- Department of Anthropology and Archaeology, National University of Mongolia; Ulaanbaatar-51, Mongolia
| | - William Honeychurch
- Department of Anthropology, Yale University, 10 Sachem St., New Haven, CT 06511, USA
| | - Joshua Wright
- Department of Archaeology, University of Aberdeen, King’s College; Aberdeen, AB24 3FX, UK
| | - Sarah Pleuger
- School of History, Classics and Archaeology, University of Edinburgh; Edinburgh EH8 9AG, UK
| | - Myagmar Erdene
- Department of Anthropology and Archaeology, National University of Mongolia; Ulaanbaatar-51, Mongolia
| | - Dashtseveg Tumen
- Department of Anthropology and Archaeology, National University of Mongolia; Ulaanbaatar-51, Mongolia
| | - Leland Rogers
- Department of Anthropology, University of North Carolina Wilmington; Wilmington, NC 28403, USA
| | - Dorjpurev Khatanbaatar
- School of Business Administration and Humanities, The Mongolian University of Science and Technology; Mongolia
| | - Byambatseren Batdalai
- Archaeological Research Center, National University of Mongolia; Ulaanbaatar-51, Mongolia
| | - Ganbaatar Galdan
- Institute of Archaeology, Mongolian Academy of Sciences, Ulaanbaatar-51, Mongolia
| | - Lisa Janz
- Department of Anthropology, University of Toronto Scarborough; Scarborough, ON M1C 1A4, Canada
| |
Collapse
|
4
|
Kelly CJ, Chitko-McKown CG, Chuong EB. Ruminant-specific retrotransposons shape regulatory evolution of bovine immunity. Genome Res 2022; 32:1474-1486. [PMID: 35948370 PMCID: PMC9435751 DOI: 10.1101/gr.276241.121] [Citation(s) in RCA: 7] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/01/2021] [Accepted: 05/05/2022] [Indexed: 02/03/2023]
Abstract
Cattle are an important livestock species, and mapping the genomic architecture of agriculturally relevant traits such as disease susceptibility is a major challenge in the bovine research community. Lineage-specific transposable elements (TEs) are increasingly recognized to contribute to gene regulatory evolution and variation, but this possibility has been largely unexplored in ruminant genomes. We conducted epigenomic profiling of the type II interferon (IFN) response in bovine cells and found thousands of ruminant-specific TEs including MER41_BT and Bov-A2 elements predicted to act as IFN-inducible enhancer elements. CRISPR knockout experiments in bovine cells established that critical immune factors including IFNAR2 and IL2RB are transcriptionally regulated by TE-derived enhancers. Finally, population genomic analysis of 38 individuals revealed that a subset of polymorphic TE insertions may function as enhancers in modern cattle. Our study reveals that lineage-specific TEs have shaped the evolution of ruminant IFN responses and potentially continue to contribute to immune gene regulatory differences across modern breeds and individuals. Together with previous work in human cells, our findings demonstrate that lineage-specific TEs have been independently co-opted to regulate IFN-inducible gene expression in multiple species, supporting TE co-option as a recurrent mechanism driving the evolution of IFN-inducible transcriptional networks.
Collapse
Affiliation(s)
- Conor J Kelly
- Department of Molecular, Cellular, and Developmental Biology and BioFrontiers Institute, University of Colorado Boulder, Boulder, Colorado 80309, USA
| | - Carol G Chitko-McKown
- USDA, ARS, Roman L. Hruska US Meat Animal Research Center (MARC), Clay Center, Nebraska 68933, USA
| | - Edward B Chuong
- Department of Molecular, Cellular, and Developmental Biology and BioFrontiers Institute, University of Colorado Boulder, Boulder, Colorado 80309, USA
| |
Collapse
|
5
|
He Y, Huang Y, Wang S, Zhang L, Gao H, Zhao Y, E G. Hereditary Basis of Coat Color and Excellent Feed Conversion Rate of Red Angus Cattle by Next-Generation Sequencing Data. Animals (Basel) 2022; 12:1509. [PMID: 35739846 PMCID: PMC9219544 DOI: 10.3390/ani12121509] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/01/2022] [Revised: 05/24/2022] [Accepted: 06/07/2022] [Indexed: 12/03/2022] Open
Abstract
Angus cattle have made remarkable contributions to the livestock industry worldwide as a commercial meat-type breed. Some evidence supported that Angus cattle with different coat colors have different feed-to-meat ratios, and the genetic basis of their coat color is inconclusive. Here, genome-wide association study was performed to investigate the genetic divergence of black and red Angus cattle with 63 public genome sequencing data. General linear model analysis was used to identify genomic regions with potential candidate variant/genes that contribute to coat color and feed conversion rate. Results showed that six single nucleotide polymorphisms (SNPs) and two insertion−deletions, which were annotated in five genes (ZCCHC14, ANKRD11, FANCA, MC1R, and LOC532875 [AFG3-like protein 1]), considerably diverged between black and red Angus cattle. The strongest associated loci, namely, missense mutation CHIR18_14705671 (c.296T > C) and frameshift mutation CHIR18_12999497 (c.310G>-), were located in MC1R. Three consecutive strongly associated SNPs were also identified and located in FANCA, which is widely involved in the Fanconi anemia pathway. Several SNPs of highly associated SNPs was notably enriched in ZCCHC14 and ANKRD11, which are related to myofiber growth and muscle development. This study provides a basis for the use of potential genetic markers to be used in future breeding programs to improve cattle selection in terms of coat color and meat phenotype. This study is also helpful to understand the hereditary basis of different coat colors and meat phenotypes. However, the putative candidate genes or markers identified in this study require further investigation to confirm their phenotypic causality and potential effective genetic relationships.
Collapse
Affiliation(s)
- Yongmeng He
- College of Animal Science and Technology, Southwest University, Chongqing 400716, China; (Y.H.); (Y.H.); (S.W.); (Y.Z.)
| | - Yongfu Huang
- College of Animal Science and Technology, Southwest University, Chongqing 400716, China; (Y.H.); (Y.H.); (S.W.); (Y.Z.)
| | - Shizhi Wang
- College of Animal Science and Technology, Southwest University, Chongqing 400716, China; (Y.H.); (Y.H.); (S.W.); (Y.Z.)
| | - Lupei Zhang
- Institute of Animal Science, Chinese Academy of Agricultural Sciences (CAAS), Beijing 100193, China; (L.Z.); (H.G.)
| | - Huijiang Gao
- Institute of Animal Science, Chinese Academy of Agricultural Sciences (CAAS), Beijing 100193, China; (L.Z.); (H.G.)
| | - Yongju Zhao
- College of Animal Science and Technology, Southwest University, Chongqing 400716, China; (Y.H.); (Y.H.); (S.W.); (Y.Z.)
| | - Guangxin E
- College of Animal Science and Technology, Southwest University, Chongqing 400716, China; (Y.H.); (Y.H.); (S.W.); (Y.Z.)
| |
Collapse
|
6
|
Chen SY, Schenkel FS, Melo ALP, Oliveira HR, Pedrosa VB, Araujo AC, Melka MG, Brito LF. Identifying pleiotropic variants and candidate genes for fertility and reproduction traits in Holstein cattle via association studies based on imputed whole-genome sequence genotypes. BMC Genomics 2022; 23:331. [PMID: 35484513 PMCID: PMC9052698 DOI: 10.1186/s12864-022-08555-z] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/05/2021] [Accepted: 04/12/2022] [Indexed: 02/06/2023] Open
Abstract
Background Genetic progress for fertility and reproduction traits in dairy cattle has been limited due to the low heritability of most indicator traits. Moreover, most of the quantitative trait loci (QTL) and candidate genes associated with these traits remain unknown. In this study, we used 5.6 million imputed DNA sequence variants (single nucleotide polymorphisms, SNPs) for genome-wide association studies (GWAS) of 18 fertility and reproduction traits in Holstein cattle. Aiming to identify pleiotropic variants and increase detection power, multiple-trait analyses were performed using a method to efficiently combine the estimated SNP effects of single-trait GWAS based on a chi-square statistic. Results There were 87, 72, and 84 significant SNPs identified for heifer, cow, and sire traits, respectively, which showed a wide and distinct distribution across the genome, suggesting that they have relatively distinct polygenic nature. The biological functions of immune response and fatty acid metabolism were significantly enriched for the 184 and 124 positional candidate genes identified for heifer and cow traits, respectively. No known biological function was significantly enriched for the 147 positional candidate genes found for sire traits. The most important chromosomes that had three or more significant QTL identified are BTA22 and BTA23 for heifer traits, BTA8 and BTA17 for cow traits, and BTA4, BTA7, BTA17, BTA22, BTA25, and BTA28 for sire traits. Several novel and biologically important positional candidate genes were strongly suggested for heifer (SOD2, WTAP, DLEC1, PFKFB4, TRIM27, HECW1, DNAH17, and ADAM3A), cow (ANXA1, PCSK5, SPESP1, and JMJD1C), and sire (ELMO1, CFAP70, SOX30, DGCR8, SEPTIN14, PAPOLB, JMJD1C, and NELL2) traits. Conclusions These findings contribute to better understand the underlying biological mechanisms of fertility and reproduction traits measured in heifers, cows, and sires, which may contribute to improve genomic evaluation for these traits in dairy cattle. Supplementary Information The online version contains supplementary material available at 10.1186/s12864-022-08555-z.
Collapse
Affiliation(s)
- Shi-Yi Chen
- Department of Animal Sciences, Purdue University, 270 S. Russell Street, West Lafayette, IN, 47907-2041, USA.,Farm Animal Genetic Resources Exploration and Innovation Key Laboratory of Sichuan Province, Sichuan Agricultural University, Chengdu, 611130, Sichuan, China
| | - Flavio S Schenkel
- Centre for Genetic Improvement of Livestock, Department of Animal Biosciences, University of Guelph, Guelph, ON, N1G 2W1, Canada
| | - Ana L P Melo
- Department of Reproduction and Animal Evaluation, Rural Federal University of Rio de Janeiro, Seropédica, RJ, 23897-000, Brazil
| | - Hinayah R Oliveira
- Department of Animal Sciences, Purdue University, 270 S. Russell Street, West Lafayette, IN, 47907-2041, USA.,Centre for Genetic Improvement of Livestock, Department of Animal Biosciences, University of Guelph, Guelph, ON, N1G 2W1, Canada
| | - Victor B Pedrosa
- Department of Animal Sciences, Purdue University, 270 S. Russell Street, West Lafayette, IN, 47907-2041, USA.,Department of Animal Sciences, State University of Ponta Grossa, Ponta Grossa, PR, 84030-900, Brazil
| | - Andre C Araujo
- Department of Animal Sciences, Purdue University, 270 S. Russell Street, West Lafayette, IN, 47907-2041, USA
| | - Melkaye G Melka
- Department of Animal and Food Science, University of Wisconsin River Falls, River Falls, WI, 54022, USA
| | - Luiz F Brito
- Department of Animal Sciences, Purdue University, 270 S. Russell Street, West Lafayette, IN, 47907-2041, USA. .,Centre for Genetic Improvement of Livestock, Department of Animal Biosciences, University of Guelph, Guelph, ON, N1G 2W1, Canada.
| |
Collapse
|
7
|
Sinding MHS, Ciucani MM, Ramos-Madrigal J, Carmagnini A, Rasmussen JA, Feng S, Chen G, Vieira FG, Mattiangeli V, Ganjoo RK, Larson G, Sicheritz-Pontén T, Petersen B, Frantz L, Gilbert MTP, Bradley DG. Kouprey ( Bos sauveli) genomes unveil polytomic origin of wild Asian Bos. iScience 2021; 24:103226. [PMID: 34712923 PMCID: PMC8531564 DOI: 10.1016/j.isci.2021.103226] [Citation(s) in RCA: 8] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/13/2021] [Revised: 08/11/2021] [Accepted: 10/01/2021] [Indexed: 12/30/2022] Open
Abstract
The evolution of the genera Bos and Bison, and the nature of gene flow between wild and domestic species, is poorly understood, with genomic data of wild species being limited. We generated two genomes from the likely extinct kouprey (Bos sauveli) and analyzed them alongside other Bos and Bison genomes. We found that B. sauveli possessed genomic signatures characteristic of an independent species closely related to Bos javanicus and Bos gaurus. We found evidence for extensive incomplete lineage sorting across the three species, consistent with a polytomic diversification of the major ancestry in the group, potentially followed by secondary gene flow. Finally, we detected significant gene flow from an unsampled Asian Bos-like source into East Asian zebu cattle, demonstrating both that the full genomic diversity and evolutionary history of the Bos complex has yet to be elucidated and that museum specimens and ancient DNA are valuable resources to do so. We generated two genomes from the likely extinct kouprey (Bos sauveli) Extensive mt and nuclear-genome-wide incomplete lineage sorting across wild Asian Bos Initial polytomic diversification of the wild Asian Bos—kouprey, banteng, and gaur
Collapse
Affiliation(s)
| | | | | | - Alberto Carmagnini
- School of Biological and Chemical Sciences, Queen Mary University of London, London, UK
| | - Jacob Agerbo Rasmussen
- Globe Institute, University of Copenhagen, Copenhagen, Denmark
- Laboratory of Genomics and Molecular Medicine, Department of Biology, University of Copenhagen, Copenhagen, Denmark
| | - Shaohong Feng
- BGI-Shenzhen, Beishan Industrial Zone, Shenzhen, China
| | - Guangji Chen
- BGI-Shenzhen, Beishan Industrial Zone, Shenzhen, China
- University of Chinese Academy of Sciences, Beijing, China
| | | | | | | | - Greger Larson
- The Palaeogenomics and Bio-Archaeology Research Network, Research Laboratory for Archaeology and History of Art, University of Oxford, Oxford, UK
| | - Thomas Sicheritz-Pontén
- Globe Institute, University of Copenhagen, Copenhagen, Denmark
- Centre of Excellence for Omics-Driven Computational Biodiscovery (COMBio), Faculty of Applied Sciences, AIMST University, Kedah, Malaysia
| | - Bent Petersen
- Globe Institute, University of Copenhagen, Copenhagen, Denmark
- Centre of Excellence for Omics-Driven Computational Biodiscovery (COMBio), Faculty of Applied Sciences, AIMST University, Kedah, Malaysia
| | - Laurent Frantz
- School of Biological and Chemical Sciences, Queen Mary University of London, London, UK
- Palaeogenomics Group, Department of Veterinary Sciences, Ludwig Maximilian University, Munich, Germany
| | - M. Thomas P. Gilbert
- Globe Institute, University of Copenhagen, Copenhagen, Denmark
- Center for Evolutionary Hologenomics, University of Copenhagen, Copenhagen, Denmark
- Norwegian University of Science and Technology, University Museum, Trondheim, Norway
| | - Daniel G. Bradley
- Smurfit Institute of Genetics, Trinity College Dublin, Dublin, Ireland
| |
Collapse
|
8
|
Dixit SP, Bhatia AK, Ganguly I, Singh S, Dash S, Sharma A, Anandkumar N, Dang AK, Jayakumar S. Genome analyses revealed genetic admixture and selection signatures in Bos indicus. Sci Rep 2021; 11:21924. [PMID: 34753978 PMCID: PMC8578574 DOI: 10.1038/s41598-021-01144-2] [Citation(s) in RCA: 10] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/01/2021] [Accepted: 10/21/2021] [Indexed: 11/20/2022] Open
Abstract
The genomic diversity and relationship among seven diverse cattle breeds viz. Sahiwal, Tharparkar, Gir, Vechur, Ongole, Kangayam and Hariana were investigated in 132 random samples based on high density SNP array comprising > 777 K SNPs. A total of 1993 SNPs (0.25% of the total) having greater power (FST ≥ 0.20) to differentiate these cattle populations were identified, and utilized to partition genome of each animal into a predefined number of clusters. The structure of these cattle indicated shared ancestry of dairy breeds viz. Gir, Tharparkar and Sahiwal. Most of the animals (> 76%) of different populations under study except Vechur clustered into their own group of animals called breed. Vechur population retained highest rate of admixture, consistent with its crossing with other breeds. Ongole, Kangayam and Hariana shared comparatively less of their genome (≤ 15%) with other breeds. The study indicated that all seven breeds evolved from their independent ancestry but there was intermixing of these breeds in the recent past. The selection signatures identified between draft (Kangayam) and dairy breeds included several genes like FAM19A2, RAB31P, BEST3, DGKA, AHCY, PIGU and PFKP which are involved in immune response, metabolic pathway, transportation of glucose and sugars, signaling pathways, cellular processes, cell division and glycolysis regulation, respectively. Moreover, these genomic regions also harbour QTLs affecting milk performance traits. The signatures were also identified even between the dairy breeds. In comparison to large-sized cattle, there were significant differences in the number of QTLs affecting production (body weight, growth rate etc.) and morphological traits (height) in short-statured Vechur breed. The presence of HMGA2 gene in the selection signature on chromosome 5 may explain the variations in stature between these cattle.
Collapse
Affiliation(s)
- S P Dixit
- ICAR - National Bureau of Animal Genetic Resources, Karnal, Haryana, 132001, India.
| | - A K Bhatia
- ICAR - National Bureau of Animal Genetic Resources, Karnal, Haryana, 132001, India
| | - Indrajit Ganguly
- ICAR - National Bureau of Animal Genetic Resources, Karnal, Haryana, 132001, India
| | - Sanjeev Singh
- ICAR - National Bureau of Animal Genetic Resources, Karnal, Haryana, 132001, India
| | - Soumya Dash
- ICAR - National Bureau of Animal Genetic Resources, Karnal, Haryana, 132001, India
| | - Anurodh Sharma
- ICAR - National Bureau of Animal Genetic Resources, Karnal, Haryana, 132001, India
| | - N Anandkumar
- ICAR - National Bureau of Animal Genetic Resources, Karnal, Haryana, 132001, India
| | - A K Dang
- ICAR - National Dairy Research Institute, Karnal, Haryana, 132001, India
| | - S Jayakumar
- ICAR - National Bureau of Animal Genetic Resources, Karnal, Haryana, 132001, India
| |
Collapse
|
9
|
Holečková B, Schwarzbacherová V, Galdíková M, Koleničová S, Halušková J, Staničová J, Verebová V, Jutková A. Chromosomal Aberrations in Cattle. Genes (Basel) 2021; 12:1330. [PMID: 34573313 PMCID: PMC8468509 DOI: 10.3390/genes12091330] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/07/2021] [Revised: 08/17/2021] [Accepted: 08/20/2021] [Indexed: 02/04/2023] Open
Abstract
Chromosomal aberrations and their mechanisms have been studied for many years in livestock. In cattle, chromosomal abnormalities are often associated with serious reproduction-related problems, such as infertility of carriers and early mortality of embryos. In the present work, we review the mechanisms and consequences of the most important bovine chromosomal aberrations: Robertsonian translocations and reciprocal translocations. We also discuss the application of bovine cell cultures in genotoxicity studies.
Collapse
Affiliation(s)
- Beáta Holečková
- Department of Biology and Physiology, University of Veterinary Medicine and Pharmacy, Komenského 73, 041 81 Košice, Slovakia; (V.S.); (M.G.); (S.K.); (J.H.); (A.J.)
| | - Viera Schwarzbacherová
- Department of Biology and Physiology, University of Veterinary Medicine and Pharmacy, Komenského 73, 041 81 Košice, Slovakia; (V.S.); (M.G.); (S.K.); (J.H.); (A.J.)
| | - Martina Galdíková
- Department of Biology and Physiology, University of Veterinary Medicine and Pharmacy, Komenského 73, 041 81 Košice, Slovakia; (V.S.); (M.G.); (S.K.); (J.H.); (A.J.)
| | - Simona Koleničová
- Department of Biology and Physiology, University of Veterinary Medicine and Pharmacy, Komenského 73, 041 81 Košice, Slovakia; (V.S.); (M.G.); (S.K.); (J.H.); (A.J.)
| | - Jana Halušková
- Department of Biology and Physiology, University of Veterinary Medicine and Pharmacy, Komenského 73, 041 81 Košice, Slovakia; (V.S.); (M.G.); (S.K.); (J.H.); (A.J.)
| | - Jana Staničová
- First Faculty of Medicine, Charles University in Prague, Salmovská 1, 121 08 Prague, Czech Republic;
- Department of Chemistry, Biochemistry and Biophysics, University of Veterinary Medicine and Pharmacy, Komenského 73, 041 81 Košice, Slovakia;
| | - Valéria Verebová
- Department of Chemistry, Biochemistry and Biophysics, University of Veterinary Medicine and Pharmacy, Komenského 73, 041 81 Košice, Slovakia;
| | - Annamária Jutková
- Department of Biology and Physiology, University of Veterinary Medicine and Pharmacy, Komenského 73, 041 81 Košice, Slovakia; (V.S.); (M.G.); (S.K.); (J.H.); (A.J.)
| |
Collapse
|
10
|
Terry SA, Basarab JA, Guan LL, McAllister TA. Strategies to improve the efficiency of beef cattle production. CANADIAN JOURNAL OF ANIMAL SCIENCE 2021. [DOI: 10.1139/cjas-2020-0022] [Citation(s) in RCA: 14] [Impact Index Per Article: 4.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 12/16/2022]
Abstract
Globally, there are approximately one billion beef cattle, and compared with poultry and swine, beef cattle have the poorest conversion efficiency of feed to meat. However, these metrics fail to consider that beef cattle produce high-quality protein from feeds that are unsuitable for other livestock species. Strategies to improve the efficiency of beef cattle are focusing on operational and breeding management, host genetics, functional efficiency of rumen and respiratory microbiomes, and the structure and composition of feed. These strategies must also consider the health and immunity of the herd as well as the need for beef cattle to thrive in a changing environment. Genotyping can identify hybrid vigor with positive consequences for animal health, productivity, and environmental adaptability. The role of microbiome–host interactions is key in efficient nutrient digestion and host health. Microbial markers and gene expression patterns within the rumen microbiome are being used to identify hosts that are efficient at fibre digestion. Plant breeding and processing are optimizing the feed value of both forages and concentrates. Strategies to improve the efficiency of cattle production are a prerequisite for the sustainable intensification needed to satisfy the future demand for beef.
Collapse
Affiliation(s)
- Stephanie A. Terry
- Agriculture and Agri-Food Canada, Lethbridge Research and Development Centre, 5403 1st Avenue South, Lethbridge, AB T1J 4B1, Canada
| | - John A. Basarab
- Alberta Agriculture and Forestry, Lacombe Research and Development Centre, 6000 C&E Trail, Lacombe, AB T4L 1W1, Canada
| | - Le Luo Guan
- Department of Agricultural, Food and Nutritional Science, University of Alberta, Edmonton, AB T6G 2P5, Canada
| | - Tim A. McAllister
- Agriculture and Agri-Food Canada, Lethbridge Research and Development Centre, 5403 1st Avenue South, Lethbridge, AB T1J 4B1, Canada
| |
Collapse
|
11
|
Zhang K, Lenstra JA, Zhang S, Liu W, Liu J. Evolution and domestication of the Bovini species. Anim Genet 2020; 51:637-657. [PMID: 32716565 DOI: 10.1111/age.12974] [Citation(s) in RCA: 28] [Impact Index Per Article: 7.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Accepted: 06/03/2020] [Indexed: 12/17/2022]
Abstract
Domestication of the Bovini species (taurine cattle, zebu, yak, river buffalo and swamp buffalo) since the early Holocene (ca. 10 000 BCE) has contributed significantly to the development of human civilization. In this study, we review recent literature on the origin and phylogeny, domestication and dispersal of the three major Bos species - taurine cattle, zebu and yak - and their genetic interactions. The global dispersion of taurine and zebu cattle was accompanied by population bottlenecks, which resulted in a marked phylogeographic differentiation of the mitochondrial and Y-chromosomal DNA. The high diversity of European breeds has been shaped through isolation-by-distance, different production objectives, breed formation and the expansion of popular breeds. The overlapping and broad ranges of taurine and zebu cattle led to hybridization with each other and with other bovine species. For instance, Chinese gayal carries zebu mitochondrial DNA; several Indonesian zebu descend from zebu bull × banteng cow crossings; Tibetan cattle and yak have exchanged gene variants; and about 5% of the American bison contain taurine mtDNA. Analysis at the genomic level indicates that introgression may have played a role in environmental adaptation.
Collapse
Affiliation(s)
- K Zhang
- State Key Laboratory of Grassland Agro-ecosystem, Institute of Innovation Ecology and College of Life Sciences, Lanzhou University, Lanzhou, 730000, China
| | - J A Lenstra
- Faculty of Veterinary Medicine, Utrecht University, Utrecht Yalelaan 104, Utrecht, 3584 CM, The Netherlands
| | - S Zhang
- State Key Laboratory of Grassland Agro-ecosystem, Institute of Innovation Ecology and College of Life Sciences, Lanzhou University, Lanzhou, 730000, China
| | - W Liu
- State Key Laboratory of Grassland Agro-ecosystem, Institute of Innovation Ecology and College of Life Sciences, Lanzhou University, Lanzhou, 730000, China
| | - J Liu
- State Key Laboratory of Grassland Agro-ecosystem, Institute of Innovation Ecology and College of Life Sciences, Lanzhou University, Lanzhou, 730000, China
- Key Laboratory for Bio-resource and Eco-environment of Ministry of Education, College of Life Sciences, Sichuan University, Chengdu, 610065, China
| |
Collapse
|
12
|
Sun HZ, Zhao K, Zhou M, Chen Y, Guan LL. Landscape of multi-tissue global gene expression reveals the regulatory signatures of feed efficiency in beef cattle. Bioinformatics 2020; 35:1712-1719. [PMID: 30329014 DOI: 10.1093/bioinformatics/bty883] [Citation(s) in RCA: 12] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/26/2018] [Revised: 10/01/2018] [Accepted: 10/16/2018] [Indexed: 01/08/2023] Open
Abstract
MOTIVATION Feed efficiency is an important trait for sustainable beef production that is regulated by the complex biological process, but the mode of action behinds it has not been clearly defined. Here, we aimed to elucidate the regulatory mechanisms of this trait through studying the landscape of the genome-wide gene expression of rumen, liver, muscle and backfat tissues, the key ones involved in the energy metabolism. RESULTS The transcriptome of 189 samples across four tissues from 48 beef steers with varied feed efficiency were generated using Illumina HiSeq4000. The analysis of global gene expression profiles of four tissues, functional analysis of tissue-shared and -unique genes, co-expressed network construction of tissue-shared genes, weighted correlations analysis between gene modules and feed efficiency-related traits in each tissue were performed. Among four tissues, the transcriptome of muscle tissue was distinctive from others, while those of rumen and backfat tissues were similar. The associations between co-expressed genes and feed efficiency related traits at single or all tissues level exhibited that the gene expression in the rumen, liver, muscle and backfat were the most correlated with feed conversion ratio, dry matter intake, average daily gain and residual feed intake, respectively. The 19 overlapped genes identified from the strongest module-trait relationships in four tissues are potential generic gene markers for feed efficiency. AVAILABILITY AND IMPLEMENTATION The distribution of gene expression data can be accessed at https://www.cattleomics.com/transcriptome. SUPPLEMENTARY INFORMATION Supplementary data are available at Bioinformatics online.
Collapse
Affiliation(s)
- Hui-Zeng Sun
- Department of Agricultural, Food and Nutritional Science, University of Alberta, Edmonton, AB, Canada
| | - Ke Zhao
- Department of Agricultural, Food and Nutritional Science, University of Alberta, Edmonton, AB, Canada.,Department of Food Quality and Safety, College of Food Engineering and Nutritional Science, Shaanxi Normal University, Xian, China
| | - Mi Zhou
- Department of Agricultural, Food and Nutritional Science, University of Alberta, Edmonton, AB, Canada
| | - Yanhong Chen
- Department of Agricultural, Food and Nutritional Science, University of Alberta, Edmonton, AB, Canada
| | - Le Luo Guan
- Department of Agricultural, Food and Nutritional Science, University of Alberta, Edmonton, AB, Canada
| |
Collapse
|
13
|
Butty AM, Chud TCS, Miglior F, Schenkel FS, Kommadath A, Krivushin K, Grant JR, Häfliger IM, Drögemüller C, Cánovas A, Stothard P, Baes CF. High confidence copy number variants identified in Holstein dairy cattle from whole genome sequence and genotype array data. Sci Rep 2020; 10:8044. [PMID: 32415111 PMCID: PMC7229195 DOI: 10.1038/s41598-020-64680-3] [Citation(s) in RCA: 13] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/19/2019] [Accepted: 04/15/2020] [Indexed: 12/15/2022] Open
Abstract
Multiple methods to detect copy number variants (CNV) relying on different types of data have been developed and CNV have been shown to have an impact on phenotypes of numerous traits of economic importance in cattle, such as reproduction and immunity. Further improvements in CNV detection are still needed in regard to the trade-off between high-true and low-false positive variant identification rates. Instead of improving single CNV detection methods, variants can be identified in silico with high confidence when multiple methods and datasets are combined. Here, CNV were identified from whole-genome sequences (WGS) and genotype array (GEN) data on 96 Holstein animals. After CNV detection, two sets of high confidence CNV regions (CNVR) were created that contained variants found in both WGS and GEN data following an animal-based (n = 52) and a population-based (n = 36) pipeline. Furthermore, the change in false positive CNV identification rates using different GEN marker densities was evaluated. The population-based approach characterized CNVR, which were more often shared among animals (average 40% more samples per CNVR) and were more often linked to putative functions (48 vs 56% of CNVR) than CNV identified with the animal-based approach. Moreover, false positive identification rates up to 22% were estimated on GEN information. Further research using larger datasets should use a population-wide approach to identify high confidence CNVR.
Collapse
Affiliation(s)
- Adrien M Butty
- Centre for Genetic Improvement of Livestock, University of Guelph, Guelph, ON, Canada
| | - Tatiane C S Chud
- Centre for Genetic Improvement of Livestock, University of Guelph, Guelph, ON, Canada
| | - Filippo Miglior
- Centre for Genetic Improvement of Livestock, University of Guelph, Guelph, ON, Canada
| | - Flavio S Schenkel
- Centre for Genetic Improvement of Livestock, University of Guelph, Guelph, ON, Canada
| | - Arun Kommadath
- Department of Agricultural, Food and Nutritional Science, University of Alberta, Edmonton, AB, Canada.,Lacombe Research and Development Centre, Agriculture and Agri-Food Canada, Lacombe, AB, Canada
| | - Kirill Krivushin
- Department of Agricultural, Food and Nutritional Science, University of Alberta, Edmonton, AB, Canada
| | - Jason R Grant
- Department of Agricultural, Food and Nutritional Science, University of Alberta, Edmonton, AB, Canada
| | - Irene M Häfliger
- Institute of Genetics, Vetsuisse Faculty, University of Bern, Bern, BE, Switzerland
| | - Cord Drögemüller
- Institute of Genetics, Vetsuisse Faculty, University of Bern, Bern, BE, Switzerland
| | - Angela Cánovas
- Centre for Genetic Improvement of Livestock, University of Guelph, Guelph, ON, Canada
| | - Paul Stothard
- Department of Agricultural, Food and Nutritional Science, University of Alberta, Edmonton, AB, Canada
| | - Christine F Baes
- Centre for Genetic Improvement of Livestock, University of Guelph, Guelph, ON, Canada. .,Institute of Genetics, Vetsuisse Faculty, University of Bern, Bern, BE, Switzerland.
| |
Collapse
|
14
|
Szmatoła T, Gurgul A, Jasielczuk I, Fu W, Ropka-Molik K. A detailed characteristics of bias associated with long runs of homozygosity identification based on medium density SNP microarrays. J Genomics 2020; 8:43-48. [PMID: 32328205 PMCID: PMC7171384 DOI: 10.7150/jgen.39147] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/08/2019] [Accepted: 02/18/2020] [Indexed: 01/18/2023] Open
Abstract
In the present study, runs of homozygosity (ROH) detected with the use of a standard bovine 54k single nucleotide polymorphism (SNP) genotyping assay and two different ROH detection approaches, based on 50 (M1) or 15 (M2) consecutive SNPs, were compared with results of whole genome sequencing. Both microarray-based methods accurately recognised medium-sized ROH, however, it was found that M2 method seemed to better than M1 identify short ROH, but highly overestimated their number, leading to numerous false positive calls. Moreover, long ROH identified with microarray data tended to break into shorter segments in sequencing data because of the presence of regions with high heterozygosity within the ROH sequences. This may indicate, that these long ROH are formed by closely positioned shorter homozygous segments that may be of older origin or may be created by two similar but not identical haplotypes, showing minor internal recombination signs. Such finding also suggests that at least some of the results of previous studies in regard to long ROH may be biased leading to inaccurate estimations of genomes autozygosity via ROH classification into length categories.
Collapse
Affiliation(s)
- Tomasz Szmatoła
- University Centre of Veterinary Medicine, University of Agriculture in Kraków, Al. Mickiewicza 24/28, 30-059 Kraków, Poland.,National Research Institute of Animal Production, Department of Animal Molecular Biology, Krakowska 1, 32-083 Balice, Poland
| | - Artur Gurgul
- University Centre of Veterinary Medicine, University of Agriculture in Kraków, Al. Mickiewicza 24/28, 30-059 Kraków, Poland.,National Research Institute of Animal Production, Department of Animal Molecular Biology, Krakowska 1, 32-083 Balice, Poland
| | - Igor Jasielczuk
- University Centre of Veterinary Medicine, University of Agriculture in Kraków, Al. Mickiewicza 24/28, 30-059 Kraków, Poland.,National Research Institute of Animal Production, Department of Animal Molecular Biology, Krakowska 1, 32-083 Balice, Poland
| | - Weiwei Fu
- College of Animal Science and Technology, Northwest A&F University, Yangling, Shaanxi 712100, China
| | - Katarzyna Ropka-Molik
- National Research Institute of Animal Production, Department of Animal Molecular Biology, Krakowska 1, 32-083 Balice, Poland
| |
Collapse
|
15
|
Chen N, Fu W, Zhao J, Shen J, Chen Q, Zheng Z, Chen H, Sonstegard TS, Lei C, Jiang Y. BGVD: An Integrated Database for Bovine Sequencing Variations and Selective Signatures. GENOMICS, PROTEOMICS & BIOINFORMATICS 2020; 18:186-193. [PMID: 32540200 PMCID: PMC7646086 DOI: 10.1016/j.gpb.2019.03.007] [Citation(s) in RCA: 42] [Impact Index Per Article: 10.5] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 06/30/2018] [Revised: 11/16/2018] [Accepted: 03/29/2019] [Indexed: 11/30/2022]
Abstract
Next-generation sequencing has yielded a vast amount of cattle genomic data for global characterization of population genetic diversity and identification of genomic regions under natural and artificial selection. However, efficient storage, querying, and visualization of such large datasets remain challenging. Here, we developed a comprehensive database, the Bovine Genome Variation Database (BGVD). It provides six main functionalities: gene search, variation search, genomic signature search, Genome Browser, alignment search tools, and the genome coordinate conversion tool. BGVD contains information on genomic variations comprising ~60.44 M SNPs, ~6.86 M indels, 76,634 CNV regions, and signatures of selective sweeps in 432 samples from modern cattle worldwide. Users can quickly retrieve distribution patterns of these variations for 54 cattle breeds through an interactive source of breed origin map, using a given gene symbol or genomic region for any of the three versions of the bovine reference genomes (ARS-UCD1.2, UMD3.1.1, and Btau 5.0.1). Signals of selection sweep are displayed as Manhattan plots and Genome Browser tracks. To further investigate and visualize the relationships between variants and signatures of selection, the Genome Browser integrates all variations, selection data, and resources, from NCBI, the UCSC Genome Browser, and Animal QTLdb. Collectively, all these features make the BGVD a useful archive for in-depth data mining and analyses of cattle biology and cattle breeding on a global scale. BGVD is publicly available at http://animal.nwsuaf.edu.cn/BosVar.
Collapse
Affiliation(s)
- Ningbo Chen
- Key Laboratory of Animal Genetics, Breeding and Reproduction Shaanxi Province, College of Animal Science and Technology, Northwest A&F University, Yangling 712100, China
| | - Weiwei Fu
- Key Laboratory of Animal Genetics, Breeding and Reproduction Shaanxi Province, College of Animal Science and Technology, Northwest A&F University, Yangling 712100, China
| | - Jianbang Zhao
- College of Information Engineering, Northwest A&F University, Yangling 712100, China
| | - Jiafei Shen
- Key Laboratory of Animal Genetics, Breeding and Reproduction Shaanxi Province, College of Animal Science and Technology, Northwest A&F University, Yangling 712100, China
| | - Qiuming Chen
- Key Laboratory of Animal Genetics, Breeding and Reproduction Shaanxi Province, College of Animal Science and Technology, Northwest A&F University, Yangling 712100, China
| | - Zhuqing Zheng
- Key Laboratory of Animal Genetics, Breeding and Reproduction Shaanxi Province, College of Animal Science and Technology, Northwest A&F University, Yangling 712100, China
| | - Hong Chen
- Key Laboratory of Animal Genetics, Breeding and Reproduction Shaanxi Province, College of Animal Science and Technology, Northwest A&F University, Yangling 712100, China
| | | | - Chuzhao Lei
- Key Laboratory of Animal Genetics, Breeding and Reproduction Shaanxi Province, College of Animal Science and Technology, Northwest A&F University, Yangling 712100, China
| | - Yu Jiang
- Key Laboratory of Animal Genetics, Breeding and Reproduction Shaanxi Province, College of Animal Science and Technology, Northwest A&F University, Yangling 712100, China.
| |
Collapse
|
16
|
Shen J, Hanif Q, Cao Y, Yu Y, Lei C, Zhang G, Zhao Y. Whole Genome Scan and Selection Signatures for Climate Adaption in Yanbian Cattle. Front Genet 2020; 11:94. [PMID: 32180793 PMCID: PMC7059643 DOI: 10.3389/fgene.2020.00094] [Citation(s) in RCA: 21] [Impact Index Per Article: 5.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/26/2019] [Accepted: 01/28/2020] [Indexed: 12/15/2022] Open
Abstract
Yanbian cattle is inhabitant of North of China, exhibiting many phenotypic features, such as long, dense body hair, and abundant intramuscular fat, designed to combat the extreme cold climate adaption. In the current study, we studied the cold tolerance of nine Yanbian cattle by whole genome resequencing and compared with African tropical cattle, N’Dama, as a control group. Yanbian cattle was aligned to the Bos taurus reference genome (ARS-UCD1.2) yielding an average of 10.8 fold coverage. The positive selective sweep analysis for the cold adaption in Yanbian cattle were analyzed using composite likelihood ratio (CLR) and nucleotide diversity (θπ), resulting in 292 overlapped genes. The strongest selective signal was found on BTA16 with potential mutation in CORT gene, a regulatory gene of primary hormone in the hypothalamic-pituitary-adrenal (HPA) axis, is reported to be associated with the cold stress, representedfour missense mutations (c.269C > T, p.Lys90Ile; c.251A > G, p.Glu84Gly; c.112C > T, p.Pro38Ser; c.86G > A, p.Pro29His). Meanwhile another gene on BTA6, showed significantly higher selective sweep signals for a cold adapted trait for hair follicle and length development, FGF5 (fibroblast growth factor 5) with a missense mutation (c.191C > T, p.Ser64Phe). Moreover, cold adapted Yanbian cattle was statistically compared with the hot adapted N’Dama cattle, a taurine cattle reported to show superior heat tolerance than zebu cattle, making them better adapted to the hot regions of Africa. XP-CLR, Fst, and θπ ratio were used to compare both breeds, yielding 487, 924, and 346 genes respectively. Among the 12 overlapped genes, (CD36) (c.638A > G, p.Lys 213Arg) involved in fat digestion and absorption plays an important role in membrane transport of long-chain fatty acid and its expression could increase in cold exposure. Henceforth, our study provides a novel genetic insights into the cold climate adaptation of Yanbian cattle and identified three candidate genes (CORT, FGF5, and CD36), which can add to an understanding of the cold climate adaptation of Yanbian cattle.
Collapse
Affiliation(s)
- Jiafei Shen
- Key Laboratory of Beef Cattle Genetics and Breeding in Ministry of Agriculture and Rural Agriculture, Branch of Animal Husbandry, Jilin Academy of Agricultural Sciences, Changchun, China.,College of Animal Science and Technology, Northwest A&F University, Yangling, China
| | - Quratulain Hanif
- Computational Biology Laboratory, Department of Agricutural Biotechnology, National Institute for Biotechnology and Genetic Engineering, Faisalabad, Pakistan.,Pakistan Institute of Engineering and Applied Sciences, Nilore, Pakistan
| | - Yang Cao
- Key Laboratory of Beef Cattle Genetics and Breeding in Ministry of Agriculture and Rural Agriculture, Branch of Animal Husbandry, Jilin Academy of Agricultural Sciences, Changchun, China
| | - Yongsheng Yu
- Key Laboratory of Beef Cattle Genetics and Breeding in Ministry of Agriculture and Rural Agriculture, Branch of Animal Husbandry, Jilin Academy of Agricultural Sciences, Changchun, China
| | - Chuzhao Lei
- College of Animal Science and Technology, Northwest A&F University, Yangling, China
| | - Guoliang Zhang
- Key Laboratory of Beef Cattle Genetics and Breeding in Ministry of Agriculture and Rural Agriculture, Branch of Animal Husbandry, Jilin Academy of Agricultural Sciences, Changchun, China
| | - Yumin Zhao
- Key Laboratory of Beef Cattle Genetics and Breeding in Ministry of Agriculture and Rural Agriculture, Branch of Animal Husbandry, Jilin Academy of Agricultural Sciences, Changchun, China
| |
Collapse
|
17
|
Kommadath A, Grant JR, Krivushin K, Butty AM, Baes CF, Carthy TR, Berry DP, Stothard P. A large interactive visual database of copy number variants discovered in taurine cattle. Gigascience 2020; 8:5523204. [PMID: 31241156 PMCID: PMC6593363 DOI: 10.1093/gigascience/giz073] [Citation(s) in RCA: 14] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/10/2018] [Revised: 02/27/2019] [Accepted: 05/28/2019] [Indexed: 12/19/2022] Open
Abstract
BACKGROUND Copy number variants (CNVs) contribute to genetic diversity and phenotypic variation. We aimed to discover CNVs in taurine cattle using a large collection of whole-genome sequences and to provide an interactive database of the identified CNV regions (CNVRs) that includes visualizations of sequence read alignments, CNV boundaries, and genome annotations. RESULTS CNVs were identified in each of 4 whole-genome sequencing datasets, which together represent >500 bulls from 17 breeds, using a popular multi-sample read-depth-based algorithm, cn.MOPS. Quality control and CNVR construction, performed dataset-wise to avoid batch effects, resulted in 26,223 CNVRs covering 107.75 unique Mb (4.05%) of the bovine genome. Hierarchical clustering of samples by CNVR genotypes indicated clear separation by breeds. An interactive HTML database was created that allows data filtering options, provides graphical and tabular data summaries including Hardy-Weinberg equilibrium tests on genotype proportions, and displays genes and quantitative trait loci at each CNVR. Notably, the database provides sequence read alignments at each CNVR genotype and the boundaries of constituent CNVs in individual samples. Besides numerous novel discoveries, we corroborated the genotypes reported for a CNVR at the KIT locus known to be associated with the piebald coat colour phenotype in Hereford and some Simmental cattle. CONCLUSIONS We present a large comprehensive collection of taurine cattle CNVs in a novel interactive visual database that displays CNV boundaries, read depths, and genome features for individual CNVRs, thus providing users with a powerful means to explore and scrutinize CNVRs of interest more thoroughly.
Collapse
Affiliation(s)
- Arun Kommadath
- Department of Agricultural, Food and Nutritional Science (AFNS), University of Alberta, Edmonton, AB, Canada.,Lacombe Research and Development Centre, Agriculture and Agri-Food Canada, Lacombe, Alberta, Canada
| | - Jason R Grant
- Department of Agricultural, Food and Nutritional Science (AFNS), University of Alberta, Edmonton, AB, Canada
| | - Kirill Krivushin
- Department of Agricultural, Food and Nutritional Science (AFNS), University of Alberta, Edmonton, AB, Canada
| | - Adrien M Butty
- Centre for Genetic Improvement of Livestock, Department of Animal Biosciences, University of Guelph, Guelph, ON, Canada
| | - Christine F Baes
- Centre for Genetic Improvement of Livestock, Department of Animal Biosciences, University of Guelph, Guelph, ON, Canada.,Institute of Genetics, Vetsuisse Faculty, University of Bern, Bern, Switzerland
| | - Tara R Carthy
- Teagasc, Animal & Grassland Research and Innovation Centre, Moorepark, Fermoy, Ireland
| | - Donagh P Berry
- Teagasc, Animal & Grassland Research and Innovation Centre, Moorepark, Fermoy, Ireland
| | - Paul Stothard
- Department of Agricultural, Food and Nutritional Science (AFNS), University of Alberta, Edmonton, AB, Canada
| |
Collapse
|
18
|
Zhang F, Wang Y, Mukiibi R, Chen L, Vinsky M, Plastow G, Basarab J, Stothard P, Li C. Genetic architecture of quantitative traits in beef cattle revealed by genome wide association studies of imputed whole genome sequence variants: I: feed efficiency and component traits. BMC Genomics 2020; 21:36. [PMID: 31931702 PMCID: PMC6956504 DOI: 10.1186/s12864-019-6362-1] [Citation(s) in RCA: 25] [Impact Index Per Article: 6.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/21/2019] [Accepted: 12/02/2019] [Indexed: 01/27/2023] Open
Abstract
BACKGROUND Genome wide association studies (GWAS) on residual feed intake (RFI) and its component traits including daily dry matter intake (DMI), average daily gain (ADG), and metabolic body weight (MWT) were conducted in a population of 7573 animals from multiple beef cattle breeds based on 7,853,211 imputed whole genome sequence variants. The GWAS results were used to elucidate genetic architectures of the feed efficiency related traits in beef cattle. RESULTS The DNA variant allele substitution effects approximated a bell-shaped distribution for all the traits while the distribution of additive genetic variances explained by single DNA variants followed a scaled inverse chi-squared distribution to a greater extent. With a threshold of P-value < 1.00E-05, 16, 72, 88, and 116 lead DNA variants on multiple chromosomes were significantly associated with RFI, DMI, ADG, and MWT, respectively. In addition, lead DNA variants with potentially large pleiotropic effects on DMI, ADG, and MWT were found on chromosomes 6, 14 and 20. On average, missense, 3'UTR, 5'UTR, and other regulatory region variants exhibited larger allele substitution effects in comparison to other functional classes. Intergenic and intron variants captured smaller proportions of additive genetic variance per DNA variant. Instead 3'UTR and synonymous variants explained a greater amount of genetic variance per DNA variant for all the traits examined while missense, 5'UTR and other regulatory region variants accounted for relatively more additive genetic variance per sequence variant for RFI and ADG, respectively. In total, 25 to 27 enriched cellular and molecular functions were identified with lipid metabolism and carbohydrate metabolism being the most significant for the feed efficiency traits. CONCLUSIONS RFI is controlled by many DNA variants with relatively small effects whereas DMI, ADG, and MWT are influenced by a few DNA variants with large effects and many DNA variants with small effects. Nucleotide polymorphisms in regulatory region and synonymous functional classes play a more important role per sequence variant in determining variation of the feed efficiency traits. The genetic architecture as revealed by the GWAS of the imputed 7,853,211 DNA variants will improve our understanding on the genetic control of feed efficiency traits in beef cattle.
Collapse
Affiliation(s)
- Feng Zhang
- Lacombe Research and Development Centre, Agriculture and Agri-Food Canada, Lacombe, AB, Canada.,Department of Agricultural, Food and Nutritional Science, University of Alberta, Edmonton, AB, Canada.,State Key Laboratory for Swine Genetics, Breeding and Production Technology, Jiangxi Agricultural University, Nanchang, Jiangxi, China.,Present Address: Institute of Translational Medicine, Nanchang University, Nanchang, Jiangxi, China
| | - Yining Wang
- Lacombe Research and Development Centre, Agriculture and Agri-Food Canada, Lacombe, AB, Canada.,Department of Agricultural, Food and Nutritional Science, University of Alberta, Edmonton, AB, Canada
| | - Robert Mukiibi
- Department of Agricultural, Food and Nutritional Science, University of Alberta, Edmonton, AB, Canada
| | - Liuhong Chen
- Lacombe Research and Development Centre, Agriculture and Agri-Food Canada, Lacombe, AB, Canada.,Department of Agricultural, Food and Nutritional Science, University of Alberta, Edmonton, AB, Canada
| | - Michael Vinsky
- Lacombe Research and Development Centre, Agriculture and Agri-Food Canada, Lacombe, AB, Canada
| | - Graham Plastow
- Department of Agricultural, Food and Nutritional Science, University of Alberta, Edmonton, AB, Canada
| | - John Basarab
- Alberta Agriculture and Forestry, Lacombe Research and Development Centre, 6000 C&E Trail, Lacombe, AB, Canada
| | - Paul Stothard
- Department of Agricultural, Food and Nutritional Science, University of Alberta, Edmonton, AB, Canada
| | - Changxi Li
- Lacombe Research and Development Centre, Agriculture and Agri-Food Canada, Lacombe, AB, Canada. .,Department of Agricultural, Food and Nutritional Science, University of Alberta, Edmonton, AB, Canada.
| |
Collapse
|
19
|
Wang Y, Zhang F, Mukiibi R, Chen L, Vinsky M, Plastow G, Basarab J, Stothard P, Li C. Genetic architecture of quantitative traits in beef cattle revealed by genome wide association studies of imputed whole genome sequence variants: II: carcass merit traits. BMC Genomics 2020; 21:38. [PMID: 31931697 PMCID: PMC6958779 DOI: 10.1186/s12864-019-6273-1] [Citation(s) in RCA: 21] [Impact Index Per Article: 5.3] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/26/2019] [Accepted: 11/12/2019] [Indexed: 12/12/2022] Open
Abstract
BACKGROUND Genome wide association studies (GWAS) were conducted on 7,853,211 imputed whole genome sequence variants in a population of 3354 to 3984 animals from multiple beef cattle breeds for five carcass merit traits including hot carcass weight (HCW), average backfat thickness (AFAT), rib eye area (REA), lean meat yield (LMY) and carcass marbling score (CMAR). Based on the GWAS results, genetic architectures of the carcass merit traits in beef cattle were elucidated. RESULTS The distributions of DNA variant allele substitution effects approximated a bell-shaped distribution for all the traits while the distribution of additive genetic variances explained by single DNA variants conformed to a scaled inverse chi-squared distribution to a greater extent. At a threshold of P-value < 10-5, 51, 33, 46, 40, and 38 lead DNA variants on multiple chromosomes were significantly associated with HCW, AFAT, REA, LMY, and CMAR, respectively. In addition, lead DNA variants with potentially large pleiotropic effects on HCW, AFAT, REA, and LMY were found on chromosome 6. On average, missense variants, 3'UTR variants, 5'UTR variants, and other regulatory region variants exhibited larger allele substitution effects on the traits in comparison to other functional classes. The amounts of additive genetic variance explained per DNA variant were smaller for intergenic and intron variants on all the traits whereas synonymous variants, missense variants, 3'UTR variants, 5'UTR variants, downstream and upstream gene variants, and other regulatory region variants captured a greater amount of additive genetic variance per sequence variant for one or more carcass merit traits investigated. In total, 26 enriched cellular and molecular functions were identified with lipid metabolisms, small molecular biochemistry, and carbohydrate metabolism being the most significant for the carcass merit traits. CONCLUSIONS The GWAS results have shown that the carcass merit traits are controlled by a few DNA variants with large effects and many DNA variants with small effects. Nucleotide polymorphisms in regulatory, synonymous, and missense functional classes have relatively larger impacts per sequence variant on the variation of carcass merit traits. The genetic architecture as revealed by the GWAS will improve our understanding on genetic controls of carcass merit traits in beef cattle.
Collapse
Affiliation(s)
- Yining Wang
- Lacombe Research and Development Centre, Agriculture and Agri-Food Canada, Lacombe, AB Canada
- Department of Agricultural, Food and Nutritional Science, University of Alberta, Edmonton, AB Canada
| | - Feng Zhang
- Lacombe Research and Development Centre, Agriculture and Agri-Food Canada, Lacombe, AB Canada
- Department of Agricultural, Food and Nutritional Science, University of Alberta, Edmonton, AB Canada
- State Key Laboratory for Swine Genetics, Breeding and Production Technology, Jiangxi Agricultural University, Nanchang, Jiangxi China
- Present Address: Institute of Translational Medicine, Nanchang University, Nanchang, Jiangxi China
| | - Robert Mukiibi
- Department of Agricultural, Food and Nutritional Science, University of Alberta, Edmonton, AB Canada
| | - Liuhong Chen
- Lacombe Research and Development Centre, Agriculture and Agri-Food Canada, Lacombe, AB Canada
- Department of Agricultural, Food and Nutritional Science, University of Alberta, Edmonton, AB Canada
| | - Michael Vinsky
- Lacombe Research and Development Centre, Agriculture and Agri-Food Canada, Lacombe, AB Canada
| | - Graham Plastow
- Department of Agricultural, Food and Nutritional Science, University of Alberta, Edmonton, AB Canada
| | - John Basarab
- Alberta Agriculture and Forestry, Lacombe Research and Development Centre, 6000 C&E Trail, Lacombe, AB Canada
| | - Paul Stothard
- Department of Agricultural, Food and Nutritional Science, University of Alberta, Edmonton, AB Canada
| | - Changxi Li
- Lacombe Research and Development Centre, Agriculture and Agri-Food Canada, Lacombe, AB Canada
- Department of Agricultural, Food and Nutritional Science, University of Alberta, Edmonton, AB Canada
| |
Collapse
|
20
|
Tijjani A, Utsunomiya YT, Ezekwe AG, Nashiru O, Hanotte O. Genome Sequence Analysis Reveals Selection Signatures in Endangered Trypanotolerant West African Muturu Cattle. Front Genet 2019; 10:442. [PMID: 31231417 PMCID: PMC6558954 DOI: 10.3389/fgene.2019.00442] [Citation(s) in RCA: 31] [Impact Index Per Article: 6.2] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/01/2018] [Accepted: 04/29/2019] [Indexed: 01/01/2023] Open
Abstract
Like most West African Bos taurus, the shorthorn Muturu is under threat of replacement or crossbreeding with zebu. Their populations are now reduced to a few hundred breeding individuals and they are considered endangered. So far, the genetic variation and genetic basis of the trypanotolerant Muturu environmental adaptation have not been assessed. Here, we present genome-wide candidate positive selection signatures in Muturu following within-population iHS and between population Rsb signatures of selection analysis. We compared the results in Muturu with the ones obtained in N’Dama, a West African longhorn trypanotolerant taurine, and in two European taurine (Holstein and Jersey). The results reveal candidate signatures of selection regions in Muturu including genes linked to the innate (e.g., TRIM10, TRIM15, TRIM40, and TRIM26) and the adaptive (e.g., JSP.1, BOLA-DQA2, BOLA-DQA5, BOLA-DRB3, and BLA-DQB) immune responses. The most significant regions are identified on BTA 23 at the bovine major histocompatibility complex (MHC) (iHS analysis) and on BTA 12 at genes including a heat tolerance gene (INTS6) (Rsb analysis). Other candidate selected regions include genes related to growth traits/stature (e.g., GHR and GHRHR), coat color (e.g., MITF and KIT), feed efficiency (e.g., ZRANB3 and MAP3K5) and reproduction (e.g., RFX2, SRY, LAP3, and GPX5). Genes under common signatures of selection regions with N’Dama, including for adaptive immunity and heat tolerance, suggest shared mechanisms of adaptation to environmental challenges for these two West African taurine cattle. Interestingly, out of the 242,910 SNPs identified within the candidate selected regions in Muturu, 917 are missense SNPs (0.4%), with an unequal distribution across 273 genes. Fifteen genes including RBBP8, NID1, TEX15, LAMA3, TRIM40, and OR12D3 comprise 220 missense variants, each between 11 and 32. Our results, while providing insights into the candidate genes under selection in Muturu, are paving the way to the identification of genes and their polymorphisms linked to the unique tropical adaptive traits of the West Africa taurine, including trypanotolerance.
Collapse
Affiliation(s)
- Abdulfatai Tijjani
- Cells, Organisms and Molecular Genetics, School of Life Sciences, University Park Campus, University of Nottingham, Nottingham, United Kingdom.,Center for Genomics Research and Innovation, National Biotechnology Development Agency, Abuja, Nigeria.,International Livestock Research Institute, Addis Ababa, Ethiopia
| | - Yuri Tani Utsunomiya
- Department of Support, Production and Animal Health, School of Veterinary Medicine, São Paulo State University, São Paulo, Brazil
| | - Arinze G Ezekwe
- Department of Animal Science, Faculty of Agriculture, University of Nigeria, Nsukka, Nigeria
| | - Oyekanmi Nashiru
- Center for Genomics Research and Innovation, National Biotechnology Development Agency, Abuja, Nigeria
| | - Olivier Hanotte
- Cells, Organisms and Molecular Genetics, School of Life Sciences, University Park Campus, University of Nottingham, Nottingham, United Kingdom.,International Livestock Research Institute, Addis Ababa, Ethiopia
| |
Collapse
|
21
|
Crysnanto D, Wurmser C, Pausch H. Accurate sequence variant genotyping in cattle using variation-aware genome graphs. Genet Sel Evol 2019; 51:21. [PMID: 31092189 PMCID: PMC6521551 DOI: 10.1186/s12711-019-0462-x] [Citation(s) in RCA: 23] [Impact Index Per Article: 4.6] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/02/2018] [Accepted: 05/03/2019] [Indexed: 12/22/2022] Open
Abstract
BACKGROUND Genotyping of sequence variants typically involves, as a first step, the alignment of sequencing reads to a linear reference genome. Because a linear reference genome represents only a small fraction of all the DNA sequence variation within a species, reference allele bias may occur at highly polymorphic or divergent regions of the genome. Graph-based methods facilitate the comparison of sequencing reads to a variation-aware genome graph, which incorporates a collection of non-redundant DNA sequences that segregate within a species. We compared the accuracy and sensitivity of graph-based sequence variant genotyping using the Graphtyper software to two widely-used methods, i.e., GATK and SAMtools, which rely on linear reference genomes using whole-genome sequencing data from 49 Original Braunvieh cattle. RESULTS We discovered 21,140,196, 20,262,913, and 20,668,459 polymorphic sites using GATK, Graphtyper, and SAMtools, respectively. Comparisons between sequence variant genotypes and microarray-derived genotypes showed that Graphtyper outperformed both GATK and SAMtools in terms of genotype concordance, non-reference sensitivity, and non-reference discrepancy. The sequence variant genotypes that were obtained using Graphtyper had the smallest number of Mendelian inconsistencies between sequence-derived single nucleotide polymorphisms and indels in nine sire-son pairs. Genotype phasing and imputation using the Beagle software improved the quality of the sequence variant genotypes for all the tools evaluated, particularly for animals that were sequenced at low coverage. Following imputation, the concordance between sequence- and microarray-derived genotypes was almost identical for the three methods evaluated, i.e., 99.32, 99.46, and 99.24% for GATK, Graphtyper, and SAMtools, respectively. Variant filtration based on commonly used criteria improved genotype concordance slightly but it also decreased sensitivity. Graphtyper required considerably more computing resources than SAMtools but less than GATK. CONCLUSIONS Sequence variant genotyping using Graphtyper is accurate, sensitive and computationally feasible in cattle. Graph-based methods enable sequence variant genotyping from variation-aware reference genomes that may incorporate cohort-specific sequence variants, which is not possible with the current implementation of state-of-the-art methods that rely on linear reference genomes.
Collapse
|
22
|
Abo-Ismail MK, Lansink N, Akanno E, Karisa BK, Crowley JJ, Moore SS, Bork E, Stothard P, Basarab JA, Plastow GS. Development and validation of a small SNP panel for feed efficiency in beef cattle. J Anim Sci 2018; 96:375-397. [PMID: 29390120 DOI: 10.1093/jas/sky020] [Citation(s) in RCA: 32] [Impact Index Per Article: 5.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/26/2017] [Accepted: 01/17/2018] [Indexed: 12/11/2022] Open
Abstract
The objective of this study was to develop and validate a customized cost-effective single nucleotide polymorphism (SNP) panel for genetic improvement of feed efficiency in beef cattle. The SNPs identified in previous association studies and through extensive analysis of candidate genomic regions and genes, were screened for their functional impact and allele frequency in Angus and Hereford breeds used as validation candidates for the panel. Association analyses were performed on genotypes of 159 SNPs from new samples of Angus (n = 160), Hereford (n = 329), and Angus-Hereford crossbred (n = 382) cattle using allele substitution and genotypic models in ASReml. Genomic heritabilities were estimated for feed efficiency traits using the full set of SNPs, SNPs associated with at least one of the traits (at P ≤ 0.05 and P < 0.10), as well as the Illumina bovine 50K representing a widely used commercial genotyping panel. A total of 63 SNPs within 43 genes showed association (P ≤ 0.05) with at least one trait. The minor alleles of SNPs located in the GHR and CAST genes were associated with decreasing effects on residual feed intake (RFI) and/or RFI adjusted for backfat (RFIf), whereas minor alleles of SNPs within MKI67 gene were associated with increasing effects on RFI and RFIf. Additionally, the minor allele of rs137400016 SNP within CNTFR was associated with increasing average daily gain (ADG). The SNPs genotypes within UMPS, SMARCAL, CCSER1, and LMCD1 genes showed significant over-dominance effects whereas other SNPs located in SMARCAL1, ANXA2, CACNA1G, and PHYHIPL genes showed additive effects on RFI and RFIf. Gene enrichment analysis indicated that gland development, as well as ion and cation transport are important physiological mechanisms contributing to variation in feed efficiency traits. The study revealed the effect of the Jak-STAT signaling pathway on feed efficiency through the CNTFR, OSMR, and GHR genes. Genomic heritability using the 63 significant (P ≤ 0.05) SNPs was 0.09, 0.09, 0.13, 0.05, 0.05, and 0.07 for ADG, dry matter intake, midpoint metabolic weight, RFI, RFIf, and backfat, respectively. These SNPs contributed to genetic variation in the studied traits and thus can potentially be used or tested to generate cost-effective molecular breeding values for feed efficiency in beef cattle.
Collapse
Affiliation(s)
- M K Abo-Ismail
- Livestock Gentec at University of Alberta, Edmonton, AB, Canada
- Animal and Poultry Production Department, Damanhour University, Damanhour, Egypt
| | - N Lansink
- Livestock Gentec at University of Alberta, Edmonton, AB, Canada
| | - E Akanno
- Livestock Gentec at University of Alberta, Edmonton, AB, Canada
| | - B K Karisa
- Livestock Gentec at University of Alberta, Edmonton, AB, Canada
| | - J J Crowley
- Livestock Gentec at University of Alberta, Edmonton, AB, Canada
- Canadian Beef Breeds Council, Calgary, AB, Canada
| | - S S Moore
- Centre for Animal Science, University of Queensland, St Lucia, Australia
| | - E Bork
- Rangeland Research Institute, Agriculture/Forestry Center, University of Alberta, Edmonton, AB, Canada
| | - P Stothard
- Livestock Gentec at University of Alberta, Edmonton, AB, Canada
| | - J A Basarab
- Alberta Agriculture and Forestry, Lacombe Research Centre, Lacombe, AB, Canada
| | - G S Plastow
- Livestock Gentec at University of Alberta, Edmonton, AB, Canada
| |
Collapse
|
23
|
Invited review: Improving feed efficiency of beef cattle – the current state of the art and future challenges. Animal 2018; 12:1815-1826. [DOI: 10.1017/s1751731118000976] [Citation(s) in RCA: 81] [Impact Index Per Article: 13.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/07/2023] Open
|
24
|
Gonen S, Ros-Freixedes R, Battagin M, Gorjanc G, Hickey JM. A method for the allocation of sequencing resources in genotyped livestock populations. Genet Sel Evol 2017; 49:47. [PMID: 28521728 PMCID: PMC5437657 DOI: 10.1186/s12711-017-0322-5] [Citation(s) in RCA: 22] [Impact Index Per Article: 3.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/10/2016] [Accepted: 05/12/2017] [Indexed: 11/18/2022] Open
Abstract
Background This paper describes a method, called AlphaSeqOpt, for the allocation of sequencing resources in livestock populations with existing phased genomic data to maximise the ability to phase and impute sequenced haplotypes into the whole population. Methods We present two algorithms. The first selects focal individuals that collectively represent the maximum possible portion of the haplotype diversity in the population. The second allocates a fixed sequencing budget among the families of focal individuals to enable phasing of their haplotypes at the sequence level. We tested the performance of the two algorithms in simulated pedigrees. For each pedigree, we evaluated the proportion of population haplotypes that are carried by the focal individuals and compared our results to a variant of the widely-used key ancestors approach and to two haplotype-based approaches. We calculated the expected phasing accuracy of the haplotypes of a focal individual at the sequence level given the proportion of the fixed sequencing budget allocated to its family. Results AlphaSeqOpt maximises the ability to capture and phase the most frequent haplotypes in a population in three ways. First, it selects focal individuals that collectively represent a larger portion of the population haplotype diversity than existing methods. Second, it selects focal individuals from across the pedigree whose haplotypes can be easily phased using family-based phasing and imputation algorithms, thus maximises the ability to impute sequence into the rest of the population. Third, it allocates more of the fixed sequencing budget to focal individuals whose haplotypes are more frequent in the population than to focal individuals whose haplotypes are less frequent. Unlike existing methods, we additionally present an algorithm to allocate part of the sequencing budget to the families (i.e. immediate ancestors) of focal individuals to ensure that their haplotypes can be phased at the sequence level, which is essential for enabling and maximising subsequent sequence imputation. Conclusions We present a new method for the allocation of a fixed sequencing budget to focal individuals and their families such that the final sequenced haplotypes, when phased at the sequence level, represent the maximum possible portion of the haplotype diversity in the population that can be sequenced and phased at that budget. Electronic supplementary material The online version of this article (doi:10.1186/s12711-017-0322-5) contains supplementary material, which is available to authorized users.
Collapse
Affiliation(s)
- Serap Gonen
- The Roslin Institute and Royal (Dick) School of Veterinary Studies, The University of Edinburgh, Easter Bush, Midlothian, Scotland, UK
| | - Roger Ros-Freixedes
- The Roslin Institute and Royal (Dick) School of Veterinary Studies, The University of Edinburgh, Easter Bush, Midlothian, Scotland, UK
| | - Mara Battagin
- The Roslin Institute and Royal (Dick) School of Veterinary Studies, The University of Edinburgh, Easter Bush, Midlothian, Scotland, UK
| | - Gregor Gorjanc
- The Roslin Institute and Royal (Dick) School of Veterinary Studies, The University of Edinburgh, Easter Bush, Midlothian, Scotland, UK
| | - John M Hickey
- The Roslin Institute and Royal (Dick) School of Veterinary Studies, The University of Edinburgh, Easter Bush, Midlothian, Scotland, UK.
| |
Collapse
|
25
|
Population Structure Analysis of Bull Genomes of European and Western Ancestry. Sci Rep 2017; 7:40688. [PMID: 28084449 PMCID: PMC5234001 DOI: 10.1038/srep40688] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/17/2016] [Accepted: 12/08/2016] [Indexed: 11/09/2022] Open
Abstract
Since domestication, population bottlenecks, breed formation, and selective breeding have radically shaped the genealogy and genetics of Bos taurus. In turn, characterization of population structure among diverse bull (males of Bos taurus) genomes enables detailed assessment of genetic resources and origins. By analyzing 432 unrelated bull genomes from 13 breeds and 16 countries, we demonstrate genetic diversity and structural complexity among the European/Western cattle population. Importantly, we relaxed a strong assumption of discrete or admixed population, by adapting latent variable models for individual-specific allele frequencies that directly capture a wide range of complex structure from genome-wide genotypes. As measured by magnitude of differentiation, selection pressure on SNPs within genes is substantially greater than that on intergenic regions. Additionally, broad regions of chromosome 6 harboring largest genetic differentiation suggest positive selection underlying population structure. We carried out gene set analysis using SNP annotations to identify enriched functional categories such as energy-related processes and multiple development stages. Our population structure analysis of bull genomes can support genetic management strategies that capture structural complexity and promote sustainable genetic breadth.
Collapse
|
26
|
Boussaha M, Michot P, Letaief R, Hozé C, Fritz S, Grohs C, Esquerré D, Duchesne A, Philippe R, Blanquet V, Phocas F, Floriot S, Rocha D, Klopp C, Capitan A, Boichard D. Construction of a large collection of small genome variations in French dairy and beef breeds using whole-genome sequences. Genet Sel Evol 2016; 48:87. [PMID: 27846802 PMCID: PMC5111192 DOI: 10.1186/s12711-016-0268-z] [Citation(s) in RCA: 11] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/03/2016] [Accepted: 11/04/2016] [Indexed: 12/22/2022] Open
Abstract
Background In recent years, several bovine genome sequencing projects were carried out with the aim of developing genomic tools to improve dairy and beef production efficiency and sustainability. Results In this study, we describe the first French cattle genome variation dataset obtained by sequencing 274 whole genomes representing several major dairy and beef breeds. This dataset contains over 28 million single nucleotide polymorphisms (SNPs) and small insertions and deletions. Comparisons between sequencing results and SNP array genotypes revealed a very high genotype concordance rate, which indicates the good quality of our data. Conclusions To our knowledge, this is the first large-scale catalog of small genomic variations in French dairy and beef cattle. This resource will contribute to the study of gene functions and population structure and also help to improve traits through genotype-guided selection. Electronic supplementary material The online version of this article (doi:10.1186/s12711-016-0268-z) contains supplementary material, which is available to authorized users.
Collapse
Affiliation(s)
- Mekki Boussaha
- GABI, INRA, AgroParisTech, Université Paris-Saclay, 78350, Jouy-en-Josas, France.
| | - Pauline Michot
- GABI, INRA, AgroParisTech, Université Paris-Saclay, 78350, Jouy-en-Josas, France.,Allice, Maison Nationale des Eleveurs, 75012, Paris, France
| | - Rabia Letaief
- GABI, INRA, AgroParisTech, Université Paris-Saclay, 78350, Jouy-en-Josas, France
| | - Chris Hozé
- GABI, INRA, AgroParisTech, Université Paris-Saclay, 78350, Jouy-en-Josas, France.,Allice, Maison Nationale des Eleveurs, 75012, Paris, France
| | - Sébastien Fritz
- GABI, INRA, AgroParisTech, Université Paris-Saclay, 78350, Jouy-en-Josas, France.,Allice, Maison Nationale des Eleveurs, 75012, Paris, France
| | - Cécile Grohs
- GABI, INRA, AgroParisTech, Université Paris-Saclay, 78350, Jouy-en-Josas, France
| | - Diane Esquerré
- GenPhySE, INRA, INPT, ENVT, Université de Toulouse, Castanet Tolosan, France
| | - Amandine Duchesne
- GABI, INRA, AgroParisTech, Université Paris-Saclay, 78350, Jouy-en-Josas, France
| | - Romain Philippe
- GMA, INRA, Université de Limoges, 87060, Limoges Cedex, France
| | | | - Florence Phocas
- GABI, INRA, AgroParisTech, Université Paris-Saclay, 78350, Jouy-en-Josas, France
| | - Sandrine Floriot
- GABI, INRA, AgroParisTech, Université Paris-Saclay, 78350, Jouy-en-Josas, France
| | - Dominique Rocha
- GABI, INRA, AgroParisTech, Université Paris-Saclay, 78350, Jouy-en-Josas, France
| | | | - Aurélien Capitan
- GABI, INRA, AgroParisTech, Université Paris-Saclay, 78350, Jouy-en-Josas, France.,Allice, Maison Nationale des Eleveurs, 75012, Paris, France
| | - Didier Boichard
- GABI, INRA, AgroParisTech, Université Paris-Saclay, 78350, Jouy-en-Josas, France
| |
Collapse
|
27
|
Heaton MP, Smith TPL, Carnahan JK, Basnayake V, Qiu J, Simpson B, Kalbfleisch TS. Using diverse U.S. beef cattle genomes to identify missense mutations in EPAS1, a gene associated with pulmonary hypertension. F1000Res 2016; 5:2003. [PMID: 27746904 PMCID: PMC5040160 DOI: 10.12688/f1000research.9254.2] [Citation(s) in RCA: 20] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Submit a Manuscript] [Subscribe] [Scholar Register] [Accepted: 10/04/2016] [Indexed: 01/08/2023] Open
Abstract
The availability of whole genome sequence (WGS) data has made it possible to discover protein variants
in silico. However, existing bovine WGS databases do not show data in a form conducive to protein variant analysis, and tend to under represent the breadth of genetic diversity in global beef cattle. Thus, our first aim was to use 96 beef sires, sharing minimal pedigree relationships, to create a searchable and publicly viewable set of mapped genomes relevant for 19 popular breeds of U.S. cattle. Our second aim was to identify protein variants encoded by the bovine endothelial PAS domain-containing protein 1 gene (
EPAS1), a gene associated with pulmonary hypertension in Angus cattle. The identity and quality of genomic sequences were verified by comparing WGS genotypes to those derived from other methods. The average read depth, genotype scoring rate, and genotype accuracy exceeded 14, 99%, and 99%, respectively. The 96 genomes were used to discover four amino acid variants encoded by
EPAS1 (E270Q, P362L, A671G, and L701F) and confirm two variants previously associated with disease (A606T and G610S). The six
EPAS1 missense mutations were verified with matrix-assisted laser desorption/ionization time-of-flight mass spectrometry assays, and their frequencies were estimated in a separate collection of 1154 U.S. cattle representing 46 breeds. A rooted phylogenetic tree of eight polypeptide sequences provided a framework for evaluating the likely order of mutations and potential impact of
EPAS1 alleles on the adaptive response to chronic hypoxia in U.S. cattle. This public, whole genome resource facilitates
in silico identification of protein variants in diverse types of U.S. beef cattle, and provides a means of translating WGS data into a practical biological and evolutionary context for generating and testing hypotheses.
Collapse
Affiliation(s)
| | | | | | | | | | | | - Theodore S Kalbfleisch
- Department of Biochemistry and Molecular Genetics, School of Medicine, University of Louisville, Louisville, USA
| |
Collapse
|
28
|
Heaton MP, Smith TP, Carnahan JK, Basnayake V, Qiu J, Simpson B, Kalbfleisch TS. Using diverse U.S. beef cattle genomes to identify missense mutations in EPAS1, a gene associated with pulmonary hypertension. F1000Res 2016; 5:2003. [PMID: 27746904 PMCID: PMC5040160 DOI: 10.12688/f1000research.9254.1] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Submit a Manuscript] [Subscribe] [Scholar Register] [Accepted: 10/04/2016] [Indexed: 08/24/2023] Open
Abstract
The availability of whole genome sequence (WGS) data has made it possible to discover protein variants in silico. However, existing bovine WGS databases do not show data in a form conducive to protein variant analysis, and tend to under represent the breadth of genetic diversity in global beef cattle. Thus, our first aim was to use 96 beef sires, sharing minimal pedigree relationships, to create a searchable and publicly viewable set of mapped genomes relevant for 19 popular breeds of U.S. cattle. Our second aim was to identify protein variants encoded by the bovine endothelial PAS domain-containing protein 1 gene ( EPAS1), a gene associated with pulmonary hypertension in Angus cattle. The identity and quality of genomic sequences were verified by comparing WGS genotypes to those derived from other methods. The average read depth, genotype scoring rate, and genotype accuracy exceeded 14, 99%, and 99%, respectively. The 96 genomes were used to discover four amino acid variants encoded by EPAS1 (E270Q, P362L, A671G, and L701F) and confirm two variants previously associated with disease (A606T and G610S). The six EPAS1 missense mutations were verified with matrix-assisted laser desorption/ionization time-of-flight mass spectrometry assays, and their frequencies were estimated in a separate collection of 1154 U.S. cattle representing 46 breeds. A rooted phylogenetic tree of eight polypeptide sequences provided a framework for evaluating the likely order of mutations and potential impact of EPAS1 alleles on the adaptive response to chronic hypoxia in U.S. cattle. This public, whole genome resource facilitates in silico identification of protein variants in diverse types of U.S. beef cattle, and provides a means of translating WGS data into a practical biological and evolutionary context for generating and testing hypotheses.
Collapse
Affiliation(s)
| | | | | | | | | | | | - Theodore S. Kalbfleisch
- Department of Biochemistry and Molecular Genetics, School of Medicine, University of Louisville, Louisville, USA
| |
Collapse
|
29
|
Raszek MM, Guan LL, Plastow GS. Use of Genomic Tools to Improve Cattle Health in the Context of Infectious Diseases. Front Genet 2016; 7:30. [PMID: 27014337 PMCID: PMC4780072 DOI: 10.3389/fgene.2016.00030] [Citation(s) in RCA: 20] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/09/2015] [Accepted: 02/18/2016] [Indexed: 12/15/2022] Open
Abstract
Although infectious diseases impose a heavy economic burden on the cattle industry, the etiology of many disorders that affect livestock is not fully elucidated, and effective countermeasures are often lacking. The main tools available until now have been vaccines, antibiotics and antiparasitic drugs. Although these have been very successful in some cases, the appearance of parasite and microbial resistance to these treatments is a cause of concern. Next-generation sequencing provides important opportunities to tackle problems associated with pathogenic illnesses. This review describes the rapid gains achieved to track disease progression, identify the pathogens involved, and map pathogen interactions with the host. Use of novel genomic tools subsequently aids in treatment development, as well as successful creation of breeding programs aimed toward less susceptible livestock. These may be important tools for mitigating the long term effects of combating infection and helping reduce the reliance on antibiotic treatment.
Collapse
Affiliation(s)
- Mikolaj M Raszek
- Livestock Gentec, Department of Agricultural, Food and Nutritional Science, University of Alberta Edmonton, AB, Canada
| | - Le L Guan
- Livestock Gentec, Department of Agricultural, Food and Nutritional Science, University of Alberta Edmonton, AB, Canada
| | - Graham S Plastow
- Livestock Gentec, Department of Agricultural, Food and Nutritional Science, University of Alberta Edmonton, AB, Canada
| |
Collapse
|