A large and diverse collection of bovine genome sequences from the Canadian Cattle Genome Project
GigaScience volume 4, Article number: 49 (2015)
The Canadian Cattle Genome Project is a large-scale international project that aims to develop genomics-based tools to enhance the efficiency and sustainability of beef and dairy production. Obtaining DNA sequence information is an important part of achieving this goal as it facilitates efforts to associate specific DNA differences with phenotypic variation. These associations can be used to guide breeding decisions and provide valuable insight into the molecular basis of traits.
We describe a dataset of 379 whole-genome sequences, taken primarily from key historic Bos taurus animals, along with the analyses that were performed to assess data quality. The sequenced animals represent ten populations relevant to beef or dairy production. Animal information (name, breed, population), sequence data metrics (mapping rate, depth, concordance), and sequence repository identifiers (NCBI BioProject and BioSample IDs) are provided to enable others to access and exploit this sequence information.
The large number of whole-genome sequences generated as a result of this project will contribute to ongoing work aiming to catalogue the variation that exists in cattle as well as efforts to improve traits through genotype-guided selection. Studies of gene function, population structure, and sequence evolution are also likely to benefit from the availability of this resource.
The sequencing dataset presented here primarily consists of key influential sires of ten cattle populations. Seven of these are the Simmental, Limousin, Angus, Charolais, Hereford, Gelbvieh, and Holstein purebreds. Holstein is an important dairy breed, whereas the others are used in beef production. The remaining three populations, referred to as Alberta, Guelph, and Beefbooster, are crossbred and composite animals with purebred, crossbred or composite sires. The Alberta population consists of research animals located at the Agriculture and Agri-Food Canada research station at Lacombe, and the University of Alberta research station at Kinsella, Alberta, Canada. Sires of these animals are mostly purebred (all breeds mentioned previously plus Red Angus) and Beefbooster animals (see below). Three of the sequenced sires are referred to as Kinsella composites (KC) of many different breeds, created at the station in the 1960s. The Guelph population consists of animals in the University of Guelph beef research herd, located at New Liskeard and Elora, Ontario, Canada. Although the herds have origins going back to the 1970s, since 1995 Angus (black and red) and Simmental sires, or sires that are a combination of these breeds (composite or hybrid sires) have predominantly been used. These research populations usually have recorded phenotypes such as feed efficiency, and have been largely selected in parallel with commercial selection trends. Beefbooster is a population (including the sequenced sires) of hybrid animals created in the 1970s. They are a composite of many other breeds including Angus, Simmental, Gelbvieh, Limousin, and Shorthorn. The purpose of creating this population was to capitalize on heterosis by crossbreeding. After generations of crossbreeding and selection, Beefbooster cattle are now, for all intents and purposes, their own breed (although not registered). Beefbooster Inc. is located in Calgary, Alberta. Breed and population information is provided for each sequenced animal in Additional file 1.
For the seven purebred populations, animals were chosen for sequencing following an in-depth pedigree analysis performed using pedigree files obtained from the respective breed associations. Analysis of the pedigrees revealed good pedigree completeness and depth . The programs CFC  and Pedig  were used to identify the most influential animals, and these were then ranked based on their total and marginal genetic contributions. Animals’ relationships with each other were also taken into account to avoid sequencing closely related animals. The top 30 ancestors were chosen for sequencing, with the goal of representing approximately 50 % of the effective genome. Younger bulls with a high number of progeny were also considered for sequencing. For the Alberta, Guelph, and Beefbooster populations, animals were selected based on number of progeny, and on relationships to animals with feed efficiency and meat quality measurements. To avoid duplication of sequencing efforts, animal selections were registered with the ongoing 1000 Bull Genomes Project .
The complete set of 379 genome sequences is described in Additional file 1, which includes the name, breed, and relevant population of each individual. An International Bull Identification (Interbull ID) number is also given to each animal. This ID consists of a three-letter breed code, followed by a three-letter country code, followed by a single letter to indicate the sex of the animal (M or F), and lastly a 12-character animal identifier. Thus the ID itself directly conveys country of registration, breed, and sex. The animal identifier can be used to retrieve additional animal information from various breed-specific databases. Identifiers for sequence retrieval from the NCBI Sequence Read Archive (SRA), and quality control measures and comments are also included in Additional file 1.
Sequencing and sequence data quality assessment
The project used the SOLiD 5500xl system to sequence 85 animals, and the HiSeq 2000 platform for the remaining 294. Standard filtering criteria were used to remove low quality reads prior to alignment. Further quality checking was performed using FastQC version 0.10.1 . Reads were mapped against the bovine genome assembly UMD 3.1 , including unassembled contigs using BWA version 0.5.9 . Following read alignment, local realignment was performed using GATK version 2.4 , and duplicates were then marked using Picard version 1.54 . Read mapping rate, duplication rate, and genome coverage without duplicate reads were determined and recorded for each sample (Additional file 1). The majority of samples yielded a high mapping rate (more than 95 % of reads mapped) and low duplicate reads rate (less than 10 %). IGV version 2.3 was used to visually inspect the read alignments for selected samples . This visualization allowed us to identify and correct a software bug in the mapping pipeline for SOLiD paired-end data.
As a further data quality check, SNP calling was performed using Samtools-0.1.18 mpileup . Sequencing genotypes were then compared to those obtained using the Illumina BovineHD BeadChip array (770 K) (Fig. 1) to establish the concordance rate. The median and mean SNP concordance rates for sequencing and genotyping are, respectively: 91.6 and 90.8 % for the SOLiD platform; and 99.4 and 98.2 % for the Illumina platform. This comparison allowed us to further assess data quality and to identify sample-handling errors. Finally, principal component analysis of genotypes from purebred animals was used to visualize the clustering of animals and to identify possible outliers (Fig. 2). Principal components 1 and 2 explain 75.9 % of the total variance. Quality control metrics and comments are included for each sequenced animal in Additional file 1.
Availability of supporting data and materials
Sequence data is available from the NCBI SRA repository under BioProjects PRJNA176557 and PRJNA256210. Supporting data are also available from the GigaScience GigaDB database .
Coancestry, inbreeding (F) and contribution
Genome analysis toolkit
Integrative genomics viewer
Single nucleotide polymorphism
University of Maryland
Stachowicz K, Larmer S, Jamrozik J, Moore SS, Miller SP. Sequencing and genotyping for the whole genome selection in Canadian beef populations. In: Proceedings of the Twentieth Conference of the Association for the Advancement of Animal Breeding and Genetics. Napier, New Zealand; 2013. p. 344–7.
Sargolzaei M, Iwaisaki H, Colleau JJ. CFC: a tool for monitoring genetic diversity. In: Proceedings of the 8th World Congress on Genetics Applied to Livestock Production. Belo Horizonte, Brazil: 2006; paper 27–28.
Boichard D. Pedig: a Fortran package for pedigree analysis suited for large populations. In: Proceedings of the 7th World Congress on Genetics Applied to Livestock Production. Montpellier, France: 2002; paper 28-13.
Daetwyler HD, Capitan A, Pausch H, Stothard P, van Binsbergen R, Brøndum RF, et al. Whole-genome sequencing of 234 bulls facilitates mapping of monogenic and complex traits in cattle. Nat Genet. 2014;46:858–65.
FastQC. http://www.bioinformatics.babraham.ac.uk/projects/fastqc/Accessed 30 Aug 2012.
Zimin AV, Delcher AL, Florea L, Kelley DR, Schatz MC, Puiu D, et al. A whole-genome assembly of the domestic cow, Bos taurus. Genome Biol. 2009;10:R42.
Li H, Durbin R. Fast and accurate short read alignment with Burrows-Wheeler transform. Bioinformatics. 2009;25:1754–60.
McKenna A, Hanna M, Banks E, Sivachenko A, Cibulskis K, Kernytsky A, et al. The Genome Analysis Toolkit: a MapReduce framework for analyzing next-generation DNA sequencing data. Genome Res. 2010;20:1297–303.
Picard. http://broadinstitute.github.io/picard/Accessed 30 Aug 2012.
Thorvaldsdóttir H, Robinson JT, Mesirov JP. Integrative Genomics Viewer (IGV): high-performance genomics data visualization and exploration. Brief Bioinform. 2013;14:178–92.
Li H, Handsaker B, Wysoker A, Fennell T, Ruan J, Homer N, et al. The Sequence Alignment/Map format and SAMtools. Bioinformatics. 2009;25:2078–9.
Stothard P, Liao X, Arantes AS, Pauw MD, Coros C, Plastow GS, et al. Bovine whole-genome sequence alignments from the Canadian Cattle Genome Project. GigaScience Database. 2015. http://dx.doi.org/10.5524/100157
Funding from the following sources supported this work: Genome Canada, Genome Alberta, Alberta Livestock and Meat Agency Ltd., Western Economic Diversification, and Alberta Innovates Bio Solutions. We would also like to acknowledge our industry partners, including the Canadian Angus, Simmental, Hereford, Charolais, Limousin and Gelbvieh associations, as well as Beefbooster Inc. for contributing samples and our international collaborators: AgResearch, BeefCRC, BIO, Semex, SRUC, Teagasc and the USDA–ARS. This research was enabled in part by support provided by WestGrid (http://www.westgrid.ca) and Compute Canada Calcul Canada (http://www.computecanada.ca).
The authors declare that they have no competing interests.
SM, SPM, PS, FS, JAB, and GSP designed the study. SPM, MS, and JAB identified animals for sequencing. XL, PS, and ASA developed the sequence analysis and management pipelines and analyzed the data. MDP, JJC, and CC assisted with sample acquisition, sample processing, and quality control. PS, ASA, and MDP drafted the manuscript. XL prepared the dataset for submission to NCBI. All authors read, revised and approved the manuscript.
Bovine genome sequences from the Canadian Cattle Genome Project. Description: Additional file 1 provides detailed information about each sequenced animal, including NCBI identifiers, Interbull ID, animal name, breed, population, sequencing platform, sequencing depth, and concordance rate. (XLSX 99 kb)
About this article
Cite this article
Stothard, P., Liao, X., Arantes, A.S. et al. A large and diverse collection of bovine genome sequences from the Canadian Cattle Genome Project. GigaSci 4, 49 (2015). https://doi.org/10.1186/s13742-015-0090-5
- Whole-genome sequencing
- Bos taurus