DoGSD: the dog and wolf genome SNP database (original) (raw)

Journal Article

,

1Laboratory for Conservation and Utilization of Bioresource & Key Laboratory for Microbial Resources of the Ministry of Education, Yunnan University, Kunming 650091, China

2State Key Laboratory of Genetic Resources and Evolution, and Yunnan Laboratory of Molecular Biology of Domestic Animals, Kunming Institute of Zoology, Chinese Academy of Sciences, Kunming 650223, China

4Kunming College of Life Science, University of Chinese Academy of Sciences,Kunming 650204, China

†

The authors wish it to be known that, in their opinion, the first three authors should be regarded as Joint First Authors.

Search for other works by this author on:

,

3Core Genomic Facility, Beijing Institute of Genomics, Chinese Academy of Sciences, Beijing 100101, China

†

The authors wish it to be known that, in their opinion, the first three authors should be regarded as Joint First Authors.

Search for other works by this author on:

,

3Core Genomic Facility, Beijing Institute of Genomics, Chinese Academy of Sciences, Beijing 100101, China

5CAS Key Laboratory of Genome Sciences and Information, Beijing Institute of Genomics, Chinese Academy of Sciences, Beijing 100101, China

†

The authors wish it to be known that, in their opinion, the first three authors should be regarded as Joint First Authors.

Search for other works by this author on:

,

3Core Genomic Facility, Beijing Institute of Genomics, Chinese Academy of Sciences, Beijing 100101, China

Search for other works by this author on:

,

1Laboratory for Conservation and Utilization of Bioresource & Key Laboratory for Microbial Resources of the Ministry of Education, Yunnan University, Kunming 650091, China

Search for other works by this author on:

,

5CAS Key Laboratory of Genome Sciences and Information, Beijing Institute of Genomics, Chinese Academy of Sciences, Beijing 100101, China

Search for other works by this author on:

,

2State Key Laboratory of Genetic Resources and Evolution, and Yunnan Laboratory of Molecular Biology of Domestic Animals, Kunming Institute of Zoology, Chinese Academy of Sciences, Kunming 650223, China

4Kunming College of Life Science, University of Chinese Academy of Sciences,Kunming 650204, China

6Department of Molecular and Cell Biology, School of Life Sciences, University of Science and Technology of China, Hefei 230026, China

Search for other works by this author on:

,

1Laboratory for Conservation and Utilization of Bioresource & Key Laboratory for Microbial Resources of the Ministry of Education, Yunnan University, Kunming 650091, China

Search for other works by this author on:

,

3Core Genomic Facility, Beijing Institute of Genomics, Chinese Academy of Sciences, Beijing 100101, China

Search for other works by this author on:

,

2State Key Laboratory of Genetic Resources and Evolution, and Yunnan Laboratory of Molecular Biology of Domestic Animals, Kunming Institute of Zoology, Chinese Academy of Sciences, Kunming 650223, China

7Department of Laboratory Medicine and Pathobiology, University of Toronto, Toronto, Canada

Search for other works by this author on:

... Show more

Accepted:

02 November 2014

Published:

17 November 2014

• PDF
• Split View
• • Article contents
  • Figures & tables
  • Video
  • Audio
  • Supplementary Data
• Cite

Cite

Bing Bai, Wen-Ming Zhao, Bi-Xia Tang, Yan-Qing Wang, Lu Wang, Zhang Zhang, He-Chuan Yang, Yan-Hu Liu, Jun-Wei Zhu, David M. Irwin, Guo-Dong Wang, Ya-Ping Zhang, DoGSD: the dog and wolf genome SNP database, Nucleic Acids Research, Volume 43, Issue D1, 28 January 2015, Pages D777–D783, https://doi.org/10.1093/nar/gku1174
Close

Navbar Search Filter Mobile Enter search term Search

Abstract

The rapid advancement of next-generation sequencing technology has generated a deluge of genomic data from domesticated dogs and their wild ancestor, grey wolves, which have simultaneously broadened our understanding of domestication and diseases that are shared by humans and dogs. To address the scarcity of single nucleotide polymorphism (SNP) data provided by authorized databases and to make SNP data more easily/friendly usable and available, we propose DoGSD (http://dogsd.big.ac.cn), the first canidae-specific database which focuses on whole genome SNP data from domesticated dogs and grey wolves. The DoGSD is a web-based, open-access resource comprising ∼19 million high-quality whole-genome SNPs. In addition to the dbSNP data set (build 139), DoGSD incorporates a comprehensive collection of SNPs from two newly sequenced samples (1 wolf and 1 dog) and collected SNPs from three latest dog/wolf genetic studies (7 wolves and 68 dogs), which were taken together for analysis with the population genetic statistics, Fst. In addition, DoGSD integrates some closely related information including SNP annotation, summary lists of SNPs located in genes, synonymous and non-synonymous SNPs, sampling location and breed information. All these features make DoGSD a useful resource for in-depth analysis in dog-/wolf-related studies.

INTRODUCTION

Dogs have been the dearest friends of humans as guardians, companions and working partners for thousands of years (1,2). Our natural curiosity, together with practical and aesthetical needs propelled the exploration of the genetic basis of the remarkable diversity of dog phenotypes (3–5). Furthermore, being one of the most thoroughly domesticated animals has made them a long-term focus and perfect model for domestication genetics studies (6–10). Researches into their parallel evolution with humans has facilitated our understanding of human evolution itself and, intriguingly, genes that cause diseases that are mutually shared, especially those related to cancer and neurological disorders (11,12). As a powerful genetic marker, single nucleotide polymorphism (SNP) plays a pivotal role in dog genetic researches referring to human selection (3,4,10), demographic studies (6), linkage and segregation analysis (5) and genome-wide association studies (13).

To date, the latest and most widely used dog SNP data set is from dbSNP (build139:ftp://ftp.ncbi.nlm.nih.gov/snp/), where ∼80% of the ∼2.7M SNPs were called from two dogs, one Boxer (14) and one Standard Poodle (15). It provides the SNP list from all samples without individual lists, which limits its usability in population genetic analysis. The rapid development of next-generation sequencing (NGS) technology has facilitated the generation of massive dog/wolf genome data sets (7,16). However, SNP calling from the massive data generated by NGS is very laborious and requires a large amount of computational resources.

Here, we develop DoGSD, the first web-based public database for accessing highly dense, broadly representative dog/wolf whole-genome SNP data. We collected SNPs from two unpublished dog/wolf genomes, three recently published works (8,9,11) and the latest dog SNP data set (dbSNP139). A total of 34 our samples are Chinese indigenous dogs which represent a key phase in dog domestication and are not included in any existing SNP databases. DoGSD provides a powerful SNP retrieving interface for each individual samples and a non-redundant data set. We made annotation to integrate information, such as SNP-related genes, transcripts, proteins and calculate allele frequencies. DoGSD has a functionality to search gene-related synonymous and non-synonymous SNPs. In addition, we incorporate the essential genetic statistics, viz. F- statistics (Fst) into DoGSD.

DATABASE CONSTRUCTION

Figure 1 shows the database construction pipeline. Details are described in Data sources, Data processing and Database implementation.

Figure 1.

The data sources and pipeline to construct DoGSD.

Data sources

Our data set integrates SNPs called from two newly sequenced wolf/dog samples and three published dog genetic works (8,9,11) with the newly released dog SNP data set (dbSNP139) (Table 1). The two newly sequenced samples include one grey wolf (solidGW1) sampled from Harbin, China and a Chinese indigenous dog (CID1) sampled from Ya'an, China. The 75 published canines include 7 grey wolves and 68 dogs. The 7 grey wolves are sampled across the Eurasian continent and America, with 33 indigenous dogs from China, another 34 from 5 breeds and 1 Dingo. Dogs from the 5 breeds include 11 German Shepherd Dogs (Germany), 10 Kunming Dogs (China), 11 Tibetan Mastiffs (China), 1 Belgian Malinois (Belgium) and 1 Basenji (Congo). We only used Solexa data of the 5 individuals from Freedman et al. In total, DoGSD includes 8 grey wolves, 34 Chinese indigenous dogs, 13 ancient breed dogs, 22 modern breed dogs and the dbSNP139.

Data sets in DoGSD

Data processing

Except for dbSNP139, we divided our samples into two groups. Group 1 is the two unpublished dog/wolf genomes (solidGW1 and CID1) generated by the SOLiD system and group 2 is from three recently published works (8,9,11), which were sequenced with the Illumina Solexa technology. All of the individual genomes were sequenced to an average of 15× coverage.

For the two samples sequenced by the SOLiD system, two DNA mate-paired libraries were prepared for solidGW1 with 5 and 1.5 kb insert size, which generated 72.93G bp of raw data. A single mate-paired library with 500 bp insert size was prepared for CID1, producing 45.96G bp raw data. Using the Bioscope pipeline, 42.01 Gb of solidGW1 was aligned to the Boxer genome assembly (canFam2) with 17.9× coverage, while 32.99 Gb of the CID1 was aligned to the canFam2 with 14× coverage.

We used two sets of methods and the camFam2/canFam3 reference genomes download from UCSC (http://hgdownload.soe.ucsc.edu/goldenPath/canFam2 or 3/bigZips/canFam2 or 3.fa.gz) to identify SNPs from these two groups, and integrated the SNPs from all 77 samples and the dbSNP139 data set to obtain a non-redundant SNP set. For group 1, SNPs were detected with the BioScope diBayes tool where valid adjacent two-base mismatches occurred between a read and the reference (canFam2). Depending on whether there was high or low read coverage on the reference, either a Frequentist or Bayesian algorithm was applied, respectively. The Frequentist algorithm is based on the null hypothesis that a given position is homozygous and any other valid adjacent mismatches are errors subject to a Poisson distribution. The Bayesian algorithm calculates the posterior probability of each site according to the expected polymorphism rate in the genome, GC content, coverage, position in a read and the quality value of the colour call and prior errors derived from the 6mer probe annealing error. SNP calling stringency was set as medium. Finally, we converted genome positions of all resulting SNPs from Canfam2.0 to Canfam3.0.

For group 2, we used Burrows–Wheeler Aligner (version: 0.6.2-r126) (17,18) to map Illumina short reads to the dog reference sequence (Canfam3) in the first step. Picard (version: 1.87) (downloaded from http://picard.sourceforge.net) was then used to eliminate duplicated reads generated in the library polymerase chain reaction construction. After that, we used the tools in GATK (version: 2.5–2-gf57256b) (19,20) to realign reads around known indels (downloaded from Ensembl ftp://ftp.ensembl.org/pub/release-73/variation/vcf/canis_familiaris/), and recalibrate base quality score to obtain more accurate quality score for each base. The refined data from all individuals were jointly used to call a raw SNPs set by GATK UnifiedGenotyper, and finally, we identified a high quality set of SNPs, using the variant quality score to recalibrate procedure in GATK.

For both groups and the dbSNP139 data set, we filtered SNPs with two-thirds derived alleles. SNPs were then annotated to genes, transcripts and proteins downloaded from the Ensembl FTP site (ftp://ftp.ensembl.org/pub/release-75/fasta/canis_familiaris/). Allele frequency was calculated for wolf/dog populations for each SNP from 77 sequenced samples. Also the Fst was calculated using vcftools (21) for the batch two data.

Database implementation

Several tables containing all resulting high-quality SNPs and their annotations were processed with Perl scripts and put into MySQL database. We use JSP/JAVA to implement data visualization, searching and downloading. GBrowse was also integrated for chromosome-based data visualization.

Data characteristics

In total, we obtained 19,333,098 non-redundant SNPs, approximately seven times greater than in dbSNP139. The average length of SNP intervals of dbSNP139 is 862.0 bp, with a median value of 241.0 bp, while the same statistics for DoGSD are 120.4 bp and 68.0 bp, respectively, showing that the SNP distribution of our database is denser and more uniform than that of dbSNP139. Ancestral LD/cross-breed LD blocks of the dog have an average size of 10 Kb (14). Thus, we compared the amount of intervals |$\geq$|10 Kb with no SNP of both data sets. Our results show that there were 333.4 Mb (covering 13.8% of the genome) without SNPs in dbSNP139 while only 7.4 Mb (0.3%) in DoGSD, clearly demonstrating that our database has fewer blocks without SNPs than dbSNP139 data set. This pattern is still true when calculating SNP intervals |$\geq$|10 Kb with only one (191.1 Mb/7.93% versus 1 Mb/0.04%) or two (148.6 Mb/6.17% versus 1.1 Mb/0.05%) SNPs. These superior characteristics are primarily due to the wider geographical sampling locations and the unbiased resequencing strategy. Together, higher density, more uniform distribution and fewer SNP rare blocks make our database more effective and accurate in specifying genome locations closely related to certain phenotypes and sequence sites under selection (22).

In the 77 samples collected in DoGSD, 34 are Chinese indigenous dogs, which have never been provided by any other databases. The Chinese indigenous dog has been demonstrated as the first stage in dog domestication (11) and its SNP data set should undoubtedly help researchers to profoundly understand the evolutionary gap between wolf and breed dog.

In addition to SNPs and their annotations, we also provide Fst for almost all SNP, which is an important and widely used population genetic statistic for selection detection. Users may browse and locate SNPs with potential selection signals. The Weir and Cocker-ham method (23) was used to calculate Fst between the wolf and dog populations using the inferred genotypes. Fst values provided in DoGSD are based on all 75 samples from the group 2 data set.

USAGE AND ACCESS

The DoGSD database can be accessed through a simple user interface. Online documentation is provided to help users to access the database. DoGSD has been designed with two main functionalities for data retrieval, Browse and Search. Users may browse the no-redundant and individual sample SNPs either with text format in tables (Figure 2A), or with a chromosome-based graphical GBrowse (24) interface. In the text format tables, SNP ID, sample and sequencing platform information, chromosome location, reference and derived alleles, three-fifths flank sequences are given for each SNP and annotated gene, transcript and protein, derived residual and counterparts from other samples are given if available (Figure 2B). In GBrowse interface, users can obtain Fst values, a pie chart showing the allele frequency of the dog and wolf populations, SNP density in 300 kb windows size, related gene and transcript information (Figure 2C).

Figure 2.

The screen dumps of the SNP lists of DoGSD and typical browse results. The curly arrow indicates two steps are not directly linked and the straight arrow means two steps are directly linked. (A) The non-redundant and individual sample SNP lists. (B) The text format browse result of an SNP. (C) GBrowse visualization of a SNP.

Users may type in chromosome numbers, start and end positions as input to retrieve data for either a single SNP or a group of SNPs. Comparative search of SNPs between two or more individuals were also implemented (Figure 3A). Figure 3B shows a typical comparative search result. Users can compare the SNPs of individuals interested them within a chromosome range.

Figure 3.

Example of a comparative search. (A) An example of search entry. (B) The comparative search result.

Considering that SNPs locate in genes and coding sequences (CDSs), and especially since non-synonymous SNPs are often of greater interests for further study, we provide several functions within the Browse pull-down menu: SNP Hit On Gene, SNP Coding, SNP Coding Synonymous and SNP Coding NonSynonymous. Users may obtain useful information by summaries of SNPs hit on genes or CDSs and examine whether they are synonymous or non-synonymous.

All the SNP data can be downloaded freely as tab-delimited files and bam/fastq/sra format files of the 75 cited samples are provided.

DISCUSSION

Nowadays, researches towards dog evolution, traits and human co-evolution are still hotspots in biological field. As the first database concentrated on dog/wolf whole genome SNPs, DoGSD stores huge amount of uniformly distributed high-quality SNPs, which compensates for the scarcity of wolf/dog SNPs provided by other databases. With high re-sequencing depth, sampling coverage and calling accuracy from 77 individual samples, our non-redundant SNP data set can be used as a dog SNP reference. Also, users can download individual sample SNP list from our database for population genetic analyses. Especially, SNPs indentified in the Chinese indigenous dog are valuable to researches about early phase of dog evolution. Fst value will also shed light on the selection signals in dog/wolf genomes, guiding researchers to in-depth investigation of these signals on a larger population scale and thoroughly explain the evolutionary driving forces between dog and wolf.

We will keep up with whole-genome SNP data releases and update DoGSD in a timely manner with new released data from population studies of dogs and wolves either by our own research groups or from publicly available resources. We would add uploading functionality to DoGSD by which users could submit whole-genome sequence data or SNP list directly. Once we get the whole-genome sequence, we would identify SNPs in it by group 2 pipeline described in Data processing. Interfaces for calculating Fst and θw will also be developed. We will add additional genetic statistical parameters, such as Tajima D, EHH and recombination maps to the database. Since detection of large chromosomal events promotes our understanding of complex traits and evolution from a perspective that complements SNP, we will further integrate structural variations of dog/wolf genomes into our database in the future.

FUNDING

973 program [2013CB835200 and 2013CB835202 to G.D.W.]; The Strategic Priority Research Program of the Chinese Academy of Sciences,Grant No.XDB13040500; Chinese Academy of Sciences [1731200000001 to W.Z.]; The Youth Innovation Promotion Association, Chinese Academy of Sciences [to G.D.W.]. Funding for open access charge: 973 program.

Conflict of interest statement. None declared.

†

The authors wish it to be known that, in their opinion, the first three authors should be regarded as Joint First Authors.

REFERENCES

, , , , .

Genetic evidence for an East Asian origin of domestic dogs

,

Science (New York, N.Y.)

,

2002

, vol.

298

(pg.

1610

-

1613

)

, , , , , , , , , , et al.

mtDNA data indicate a single origin for dogs south of Yangtze River, less than 16,300 years ago, from numerous wolves

,

Mol. Biol. Evol.

,

2009

, vol.

26

(pg.

2849

-

2864

)

, , , , , , , , , , et al.

Variation of BMP3 contributes to dog breed skull diversity

,

PLoS Genet.

,

2012

, vol.

8

pg.

e1002849

, , , , , , , .

Tracking footprints of artificial selection in the dog genome

,

Proc. Natl. Acad. Sci. U.S.A.

,

2010

, vol.

107

(pg.

1160

-

1165

)

, , , , , , , , , , et al.

Linkage and segregation analysis of black and brindle coat color in domestic dogs

,

Genetics

,

2007

, vol.

176

(pg.

1679

-

1689

)

, , , , , , , , , , et al.

Genome-wide SNP and haplotype analyses reveal a rich history underlying dog domestication

,

Nature

,

2010

, vol.

464

(pg.

898

-

902

)

, , , , , , , , , .

The genomic signature of dog domestication reveals adaptation to a starch-rich diet

,

Nature

,

2013

, vol.

495

(pg.

360

-

364

)

, , , , , , , , , .

Whole-genome sequencing of six dog breeds from continuous altitudes reveals adaptation to high-altitude hypoxia

Genome Res.

,

2014

, vol.

24

(pg.

1315

-

1308

)

, , , , , , , , , , et al.

Genome sequencing highlights the dynamic early history of dogs

,

PLoS Genet.

,

2014

, vol.

10

pg.

e1004016

, , , , , , , , , , et al.

Rethinking dog domestication by integrating genetics, archeology, and biogeography

,

Proc. Natl. Acad. Sci. U.S.A.

,

2012

, vol.

109

(pg.

8878

-

8883

)

, , , , , , , , , .

The genomics of selection in dogs and the parallel evolution between dogs and humans

,

Nat. Commun.

,

2013

, vol.

4

1860

, , , , , , , , , .

Genetic convergence in the adaptation of dogs and humans to the high altitude environment of the Tibetan plateau

Genome Biol. Evol.

,

2014

, vol.

6

(pg.

2122

-

2128

)

, , , , , , , , , , et al.

Coat variation in the domestic dog is governed by variants in three genes

,

Science (New York, N.Y.)

,

2009

, vol.

326

(pg.

150

-

153

)

, , , , , , , , , , et al.

Genome sequence, comparative analysis and haplotype structure of the domestic dog

,

Nature

,

2005

, vol.

438

(pg.

803

-

819

)

, , , , , , , , , , et al.

The dog genome: survey sequencing and comparative analysis

,

Science (New York, N.Y.)

,

2003

, vol.

301

(pg.

1898

-

1903

)

, , , , , , , , , , et al.

Genome analysis of the domestic dog (Korean Jindo) by massively parallel sequencing

,

DNA Res.

,

2012

, vol.

19

(pg.

275

-

287

)

, .

Fast and accurate long-read alignment with Burrows-Wheeler transform

,

Bioinformatics

,

2010

, vol.

26

(pg.

589

-

595

)

, .

Fast and accurate short read alignment with Burrows-Wheeler transform

,

Bioinformatics

,

2009

, vol.

25

(pg.

1754

-

1760

)

, , , , , , , , , , et al.

The Genome Analysis Toolkit: a MapReduce framework for analyzing next-generation DNA sequencing data

,

Genome Res.

,

2010

, vol.

20

(pg.

1297

-

1303

)

, , , , , , , , , , et al.

A framework for variation discovery and genotyping using next-generation DNA sequencing data

,

Nat. Genet.

,

2011

, vol.

43

(pg.

491

-

498

)

, , , , , , , , , , et al.

The variant call format and VCFtools

,

Bioinformatics

,

2011

, vol.

27

(pg.

2156

-

2158

)

, , , , .

Domestication genomics evidence with animals

,

Annu. Rev. Animal Biosci.

,

2014

, vol.

2

(pg.

1

-

24

)

.

Estimating F-statistics for the analysis of population structure

,

Evolution

,

1984

, vol.

38

(pg.

1358

-

1370

)

.

Using the Generic Genome Browser (GBrowse)

Curr. Protoc. Bioinformatics

,

2009

, vol.

28

(pg.

9.9.1

-

9.9.25

)

© The Author(s) 2014. Published by Oxford University Press on behalf of Nucleic Acids Research.

This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by-nc/4.0/), which permits non-commercial re-use, distribution, and reproduction in any medium, provided the original work is properly cited. For commercial re-use, please contact journals.permissions@oup.com

I agree to the terms and conditions. You must accept the terms and conditions.

Submit a comment

Name

Affiliations

Comment title

Comment

You have entered an invalid code

Thank you for submitting a comment on this article. Your comment will be reviewed and published at the journal's discretion. Please check for further notifications by email.

Citations

Views

Altmetric

Metrics

Total Views 6,143

5,049 Pageviews

1,094 PDF Downloads

Since 12/1/2016

Citations

64 Web of Science

×

Email alerts

Citing articles via

• Latest

• Most Read

• Most Cited

More from Oxford Academic