PMCCPMCCPMCC

Search tips
Search criteria 

Advanced


 
Formats
Article sections
Authors
Related links
Logo of bmcgenoBioMed Centralsearchsubmit a manuscriptregisterthis articleBMC Genomics
 
BMC Genomics. 2012; 13: 78.
Published online 2012 February 17. doi:  10.1186/1471-2164-13-78
PMCID: PMC3309927
Copyright ©2012 Doan et al; licensee BioMed Central Ltd.
Whole-Genome sequencing and genetic variant analysis of a quarter Horse mare
Ryan Doan,1 Noah D Cohen,2 Jason Sawyer,3,4 Noushin Ghaffari,5 Charlie D Johnson,5 and Scott V Dindotcorresponding author1,6
1Department of Veterinary Pathobiology, Texas A&M University College of Veterinary Medicine and Biomedical Sciences, College Station, TX, USA
2Department of Large Animal Clinical Sciences, Texas A&M University College of Veterinary Medicine, College Station, TX, USA
3Texas AgriLife Research Center, McGregor, TX, USA
4Department of Animal Science, Texas A&M University, College Station, TX, USA
5AgriLife Genomics and Bioinformatics Center, College Station, TX, USA
6Department of Molecular and Cellular Medicine, Texas A&M Health Science Center College of Medicine, College Station, TX, USA
corresponding authorCorresponding author.
Ryan Doan: rdoan/at/cvm.tamu.edu; Noah D Cohen: ncohen/at/cvm.tamu.edu; Jason Sawyer: j-sawyer/at/tamu.edu; Noushin Ghaffari: nghaffari/at/tamu.edu; Charlie D Johnson: charlie/at/ag.tamu.edu; Scott V Dindot: sdindot/at/cvm.tamu.edu
Received October 22, 2011; Accepted February 17, 2012.
This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.
Small right arrow pointing to: This article has been cited by other articles in PMC.
Abstract
Background
The catalog of genetic variants in the horse genome originates from a few select animals, the majority originating from the Thoroughbred mare used for the equine genome sequencing project. The purpose of this study was to identify genetic variants, including single nucleotide polymorphisms (SNPs), insertion/deletion polymorphisms (INDELs), and copy number variants (CNVs) in the genome of an individual Quarter Horse mare sequenced by next-generation sequencing.
Results
Using massively parallel paired-end sequencing, we generated 59.6 Gb of DNA sequence from a Quarter Horse mare resulting in an average of 24.7X sequence coverage. Reads were mapped to approximately 97% of the reference Thoroughbred genome. Unmapped reads were de novo assembled resulting in 19.1 Mb of new genomic sequence in the horse. Using a stringent filtering method, we identified 3.1 million SNPs, 193 thousand INDELs, and 282 CNVs. Genetic variants were annotated to determine their impact on gene structure and function. Additionally, we genotyped this Quarter Horse for mutations of known diseases and for variants associated with particular traits. Functional clustering analysis of genetic variants revealed that most of the genetic variation in the horse's genome was enriched in sensory perception, signal transduction, and immunity and defense pathways.
Conclusions
This is the first sequencing of a horse genome by next-generation sequencing and the first genomic sequence of an individual Quarter Horse mare. We have increased the catalog of genetic variants for use in equine genomics by the addition of novel SNPs, INDELs, and CNVs. The genetic variants described here will be a useful resource for future studies of genetic variation regulating performance traits and diseases in equids.
Background
The sequencing and assembly of a horse genome was a great achievement in equine genomics and veterinary medicine because of the broad range of potential applications of this information for improving health and performance and for understanding differences among species [1]. To date, however, only a single genome of a Thoroughbred mare has been sequenced and made publicly available [1]. The current catalog of genetic variants in the equine genome consists of 1,163,580 single nucleotide polymorphisms (SNPs; dbSNP: http://www.ncbi.nlm.nih.gov/projects/SNP/ [build 135]), with no insertion/deletion polymorphisms (INDELs) or copy number variants (CNVs) having been deposited into a publicly available database (dbSNP or dbVar, http://www.ncbi.nlm.nih.gov/dbvar). Of the known SNPs in horses, most (~64%) originate from the Thoroughbred mare used for the genome assembly [1].
Sequencing the genome of a Quarter Horse was considered important for several reasons. Almost 3 million Quarter Horses were registered in the United States in 2010 according to the American Quarter Horse Association, making it the single largest breed registry in the country [2,3]. When one includes non-registered Quarter Horses and other Quarter Horse-influenced breed registries (e.g., American Paint Horse Association), the Quarter Horse is by far the largest contributor to the population of horses in the United States. Although Quarter Horse breeding was, and continues to be, strongly influenced by Thoroughbred bloodlines, the 2 breeds were selectively bred to enhance different traits. Thoroughbred breeding has selected for speed over distances of 3/4 to 2 miles, whereas selection in the Quarter Horse has emphasized speed over shorter distances and compliant disposition suited to working cattle and other ranch-related duties [4]. Thus, there are important phenotypic differences between the breeds that have been achieved by selective breeding with a clear underlying genetic component [5,6]. Moreover, there are single-gene disorders in Quarter Horses, such as polysaccharide storage myopathy (PSSM) [7,8], hyperkalemic periodic paralysis (HYPP) [9], glycogen branching enzyme deficiency (GBED) [10], and hereditary equine regional dermal asthenia (HERDA) [11,12]. Most diseases and traits (such as predisposition to osteoarthritis or body-type) are complex, involving multiple genes, which may be modulated by environmental factors [5,13,14]. Thus, identifying genetic variants in the genome of a single Quarter Horse would provide a wealth of information for future genomic studies in equids.
Here we describe the whole-genome sequencing of an individual Quarter Horse mare using massively parallel paired-end sequencing. Sequence reads were mapped to the reference Thoroughbred nuclear and mitochondrial genomes. We developed a comprehensive list of genetic variants, including SNPs, INDELs, and CNVs. We annotated genetic variants and examined their impact on gene structure and function. The genomic sequence was also examined for mutations and polymorphisms associated with diseases and traits in horses. Furthermore, we examined biological processes enriched for genetic variants and compared these biological processes between this Quarter Horse mare and the reference Thoroughbred mare.
Results
Whole-genome sequencing, alignment, and identification of new genomic sequence
Genomic DNA was obtained from a single blood sample of a Quarter Horse mare owned by Texas A&M University. The mare's pedigree was heavily influenced by stock and racing Quarter Horses, with no introgression of Thoroughbred lines during the preceding 4 generations. A small insert library with an average size of approximately 270 base-pair (bp) was generated (Additional file 1, Figure S1), and 14 lanes of 75-bp paired-end sequencing were performed using the Illumina Genome Analyzer II (Additional file 2, Table S1) (Illumina, Inc., San Diego, CA). A schematic of the sequence mapping and variant identification, annotation, and analysis is presented in Figure ​Figure1.1. The sequencing reactions yielded 849,600,742 reads totaling 61,310,282,925 bases of DNA (Additional file 3, Table S2). Sequence reads were trimmed and mapped to the assembled autosomes, × chromosome, and mitochondrial genome of the reference Thoroughbred horse genome (EquCab2.0). A total of 817,470,439 reads (59,643,456,282 bases) were aligned to approximately 97% of the reference genome (53,783,377,648 bases were uniquely mapped), resulting in an average of 24.7X sequence coverage of the horse's genome (Additional file 3, Table S2). Sequence coverage of assembled chromosomes ranged from 13X to 55X (Additional file 4, Figure S2). Additionally, 81,993 reads mapped to 100% of the reference mitochondrial genome, resulting in an average of 355.6X sequence coverage. Reads not mapping to assembled chromosomes were subsequently mapped to the unassembled chromosomes (ChrUn), where 12,865,623 reads (935,090,584 bases) were mapped. Next, we de novo assembled 12,657,236 reads (707,253,380 bases) not mapping to the assembled or unassembled chromosomes. The de novo assembly yielded 35,540 contigs totaling 19.1 million bases (Mb) of sequence (Table ​(Table1).1). BLAT and BLAST analysis of the contigs to other genomes did not result in any significant alignment.
Figure 1
Figure 1
Overview of genome-wide variation analysis pipeline by next generation sequencing.
Table 1
Table 1
Sequence generation and mapping to equCab2 reference genome
Identification of genetic variants
We analyzed sequences mapped to the assembled chromosomes for putative SNPs, INDELs, and CNVs. SNPs were required to meet several stringent criteria (see Materials and Methods) and were further filtered to remove those near other variants (i.e., SNPs and INDELs) or within the pseudo-autosomal region (PAR) of the × chromosome. The resulting SNPs were then filtered using different minimum read-depths (5X and 10X [Additional file 5, Table S3 and Additional file 6, Table S4]). As expected, increasing the minimum read-depth resulted in fewer SNPs being detected (1,131,519 SNPs lost from 5X to 10X coverage; Additional file 7, Table S5 and Additional file 8, Figure S3A and B). Additionally, we found that a greater proportion of the total SNPs was lost than that of the total non-synonymous SNPs, with homozygous SNPs being lost at the highest rate (Additional file 8, Figure S3C and D). At 10X coverage, which was the minimum read-depth used for the subsequent annotation and downstream analysis, we identified 3,157,093 SNPs. Additionally, we characterized and annotated SNPs detected at lower minimum read-depths (Additional files 9, 10, 11, 12, 13, 14, 15, 16, 17, 18, 19, 20, 21, 22, and 23, Tables S6, S7, S8, S9, S10, S11, S12, S13, S14, S15, S16, S17, S18, S19, and S20). Intergenic SNPs, located within 1,000 bases (1 kb) upstream or downstream of a gene, were determined as well as genic SNPs located within introns, non-coding exons, 5' and 3' untranslated regions (UTR), intron and exon splice sites, and coding exons. Non-synonymous SNPs were further annotated as either radical or conservative (Table ​(Table2,2, and Additional files 9, 10, 11, 12, 13, 14, 15, 16, 17, 18, 19, 20, 21, 22, and 23, Tables S6, S7, S8, S9, S10, S11, S12, S13, S14, S15, S16, S17, S18, S19, and S20; see Materials and Methods for description of amino acid classification). Additionally, we identified SNPs located within the coding exons of 305 candidate imprinted genes for future studies in the horse ([15]; Additional file 24, Table S21). Analysis of the mitochondrial genome identified 71 SNPs, including 5 heteroplasmic SNPs (Additional file 5, Table S3 and Additional file 6, Table S4).
Table 2
Table 2
Annotation of SNPs in the Quarter Horse genome
Comparison of putative SNPs in the mapped sequences to the horse SNP database (dbSNP, http://www.ncbi.nlm.nih.gov/projects/SNP/) revealed that 342,726 were known and 2,814,367 were novel, including 18,140 non-synonymous SNPs (11,434 radical and 6,706 conservative) and 2,629 complex (i.e., tri-allelic) SNPs. Comparison of the putative SNPs to those identified in the sequencing and assembly of the reference Thoroughbred genome revealed that 244,669 SNPs were overlapping between the 2 horses, including 174,682 intergenic and 69,987 genic SNPs. By removing overlapping known SNPs, we found that the Quarter Horse and Thoroughbred genomes had 2,912,424 and 522,610 unique SNPs, respectively.
Using a minimum read-depth of 10X and similar filtering methods described for the SNP analysis, we identified 193,271 INDELs (1-8 bp insertion or deletions; Table ​Table3and3and Additional file 25, Table S22). Additionally, we identified and annotated INDELs using lower read-depth coverage (5X; Additional file 26, Table S23). Although our analysis was limited to INDELs no greater than 8 bp, the overwhelming majority were single bp gains and losses (Additional file 27, Figure S4). Of the INDELs identified, 2,574 caused a frameshift mutation (frameshift deletion + insertions) in a coding exon, whereas 35 did not (Table ​(Table3).3). In addition, 21 INDELs lead to the creation of a stop codon (stop - gain). Given the high number of INDELs affecting coding exons, particularly the large number of homozygous INDELs, we annotated them using the RefSeq annotation. The analysis revealed that 50 RefSeq genes contained INDELs, including 7 genes with loss of a stop codon and 47 with a frameshift.
Table 3
Table 3
Annotation of INDELs in the Quarter Horse genome
Copy number variants (CNVs) were identified using a read-depth algorithm that corrects for GC bias. CNVs were then filtered to remove those within telomeric regions, as investigation of these regions is associated with high false discovery rates [16]. The filtered analysis revealed 282 CNVs, including 274 gains, 6 losses, and 2 homozygous deletions (Additional file 28, Table S24, Additional file 29, Figure S5 and Additional file 30, Figure S6). The sizes ranged from 3.74 kb to 4.84 Mb, with an average length of 296.1 kb. Annotation of the CNVs indicated that 192 and 90 were genic and intergenic, respectively.
Pathway, trait, and disease analysis of identified genetic variants
We performed a functional annotation clustering analysis of genes containing SNPs, INDELs, and CNVs to identify biological processes (BP) enriched for these classes of genetic variants. Using the horse Ensembl gene identification (ID) for the functional annotation clustering, we found that only a small percentage of BP terms were returned (0.6%). Therefore, we decided to perform the analysis using the human orthologs of horse genes found to possess genetic variants. Of the 6,899 genes containing non-synonymous SNPs, 5,712 had human orthologs and 3,880 returned BP terms (67.9%). Clustering analysis indicated that non-synonymous SNPs were enriched (P < 0.05) in pathways for sensory perception (27%; P = 3.9 × 10-15), cellular processes (24%; P = 1.7 × 10-6), and signal transduction (16% P = 2.7 × 10-7; Figure ​Figure2A2A and Additional file 31, Table S25). Separate functional enrichment analysis of genes with either radical or conservative non-synonymous SNPs revealed that the biological processes of the enriched pathways remained unchanged, although there was an increase in enrichment of conservative amino acid changes in the immunity and defense pathway (3% radical and 14% conservative; Additional file 32, Figure S7A and B). Clustering analysis of INDELs indicated that genes containing variants were enriched in pathways for signal transduction (26%; P = 4.2 × 10-6), developmental processes (15%; P = 3.5 × 10-10), and cellular processes (14%; P = 1.4 × 10-5; Figure ​Figure2B).2B). Analysis of CNVs indicated that genes containing variants were enriched in pathways for sensory perception (60%; P = 2.1 × 10-61), immunity and defense (19%; P = 2.7 × 10-2), and signal transduction (13%; P = 1.3 × 10-3; Figure ​Figure2C2C).
Figure 2
Figure 2
Biological process enrichment analysis of identified genetic variants. (A) BP enrichment analysis of non-synonymous coding SNPs, (B) INDELs, and (C) CNVs.
Comparison of biological pathways containing heterozygous non-synonymous SNPs between the Quarter Horse and reference Thoroughbred (Additional file 33, Table S26) revealed that the Quarter Horse had SNPs enriched in pathways for sensory perception (36%; P = 1.3 × 10-36), cellular processes (18%; P = 4.8 × 10-5), and signal transduction (16%; P = 1.2 × 10-3; Additional file 34, Figure S8A), whereas the Thoroughbred had SNPs enriched in pathways for immunity and defense (36%; P = 5.2 × 10-8), sensory perception (23%; P = 6.4 × 10-14), and cellular processes (13%; P = 2.7 × 10-1; Additional file 34, Figure S8B). As the enrichment analysis reflected percentages of genes within pathways, and because the Quarter Horse contained nearly three times as many heterozygous SNPs as the Thoroughbred, we determined the total number of genes for each BP term that contained a non-synonymous heterozygous SNP. As expected, the Quarter Horse had more genes containing SNPs for each BP class, except for the immune and defense (P = 2.4 × 10-27) and nucleic acid metabolism (P = 1.1 × 10-11) pathways (Additional file 34, Figure S8C, Additional file 35, Tables S27 and Additional file 36, Table S28).
Next, we investigated the DNA sequence of the Quarter Horse mare for mutations of known diseases and for SNPs associated with various traits in horses (Table ​(Table4).4). We found the Quarter Horse to be heterozygous for a c.115G > A mutation in the cyclophilin B (PPIB) gene, which causes hereditary equine regional dermal asthenia (HERDA)[12], heterozygous for the g.901C > T SNP associated with chestnut coat color [17], and homozygous for the g.66493737C > T and g.22684390C > T SNPs associated with racing endurance in Thoroughbred horses [18,19]. No other mutations or polymorphisms underlying known diseases or traits were identified.
Table 4
Table 4
Analysis of genetic variants for known traits and diseases
Discussion
The current catalog of equine genetic variants is limited and primarily consists of those detected from the sequencing and assembly of a single Thoroughbred horse [1]. In the present study, we describe the whole-genome sequencing and identification of genetic variants in the genome of a Quarter Horse mare. To our knowledge, this is the first published report of a whole-genome sequence of a Quarter Horse and the only horse genome sequenced by next-generation sequencing. At 25X sequence coverage, we estimated that approximately 85-88% of the horse's genome could be genotyped [20]. Our analysis yielded 3.1 million SNPs, 193 thousand INDELs, and 282 CNVs. Despite the fact that approximately 10% of the novel homozygous SNPs likely reflect errors in the reference genome (given that the accuracy of the reference genome is 99.99%) and that our false discovery rate (FDR) for SNP detection is approximately 1.5% (Doan et al. unpublished studies), the genetic variants identified here represent a significant addition to what is currently available for studies in horses. It should be noted that the breeding structure of Quarter Horses and Thoroughbreds, as well as the differences in sequencing technology (i.e., next-generation sequencing of DNA fragments vs. Sanger sequencing of bacterial artificial clones) most likely contributed to the increased amount of genetic variation identified in this study. We suspect that many genetic variants were missed due to the parameters applied during our variant detection (e.g., requiring a minimum of 10X sequence coverage). To minimize this false negative rate (FNR) we provide a list of SNPs and INDELs using a less stringent calling criterion (5X minimum sequence coverage; Additional file 5, Table S3). However, caution should be used with these variants as the FDR is expected to increase with a reduced minimum coverage, although not directly proportional to the decrease of the FNR.
De novo assembly of reads not mapping to the assembled or un-assembled chromosomes led to the generation of 19.1 Mb of new horse genomic sequence. Our analysis of CNVs in horses and cattle using array comparative genomic hybridization (Doan et al., unpublished studies), as well as studies in the human and mouse show that large (Mb) deletions are common variants in the genome [21-24]. This presence/absence variation [25] is common, and thus we suspect that some percentage of the de novo assembled sequence represents sequence missing from the reference Thoroughbred genome due to homozygous deletions.
Functional annotation clustering analysis of genetic variants revealed that pathways for sensory perception, signal transduction, protein processing, cellular process, and immunity and defense were differentially affected by each type of genetic variant (i.e., SNPs, INDELs, and CNVs), suggesting varying degrees of tolerance and selection for genetic variants underlying these biological processes. Genes involved in sensory perception pathways contained most of the genetic variation, primarily SNPs (27%) and CNVs (60%). The observed enrichment in sensory perception genes may be related to selection of the Quarter Horse for a relatively calm disposition [26,27], although this is highly speculative at this point, and will only be determined by future population based-studies. Conversely, these variants could reflect misassemblies in the reference genome or misalignments in the Quarter Horse reads to the reference sequence, as these genes exist as large families with numerous pseudogenes.
Conclusions
This is the first sequencing of a horse genome by next-generation sequencing and the first genomic sequence of an individual Quarter Horse. The genetic variants identified in this study will be a useful resource for future studies to understand the genetic basis of phenotypic variation and disease in equids.
Methods
DNA sample
DNA was isolated from whole blood of a Quarter Horse mare using a standard phenol-chloroform extraction including 2 phenol-chloroform-isoamyl (PCI) steps, followed by rinses with chloroform, isopropanol, and 70% ethanol. The sample was suspended in Qiagen EB buffer (Qiagen Sciences, Germantown, MD). The Texas A&M University Institutional Animal Care and Use Committee approved this study.
Whole-genome sequencing
For the construction of sequencing libraries, we first sonicated high-quality genomic DNA by pulsing 3X for 15 sec/pulse at 14% using a Sonic Dismembrator 500 (Fisher Scientific, Pittsburg, PA) and purified with an Invitrogen PureLink PCR Kit (Invitrogen, Carlsbad, CA). The DNA was blunt end-repaired, adenylated, and ligated with paired-end adaptors, according to the manufacturer's protocol (Illumina, San Diego CA). The prepared library was resolved on a 2% low range agarose gel and a 2-mm section of DNA was isolated at 271 bp (Qiagen Sciences, Germantown, MD). The library was then enriched according to the manufacturer's protocol (Illumina, San Diego CA). The size and concentration of the sequencing library was determined by PCR, visualization of polyacrylamide gel electrophoresis (PAGE) gels, and through the use of the Agilent 2100 Bioanalyzer DNA kit (Agilent Technologies, San Diego CA). Cluster generation and paired-end sequencing was performed according to the manufacturer's protocols (Illumina) at the AgriLife Genomics and Bioinformatics Center (College Station, TX)
Sequence alignment
We used the trim function in the CLC Genomics Workbench (CLC Bio, Aarhus, Denmark) using the following parameters: ambiguous limit, 2; ambiguous trim, yes; quality limit, 0.1; quality trim, yes; and, remove 3' nucleotide, no; remove 5' nucleotide, no. The CLC Genomics Workbench Reference Mapping function was used to assemble the trimmed reads to the equCab2 reference assembly using the following parameters: similarity score = 0.8; and, length fraction = 0.5. Paired-end reads were mapped using an insert range of 180-bp to 500-bp and reads mapping to multiple places in the reference were mapped using the random setting. In order to determine the percent of the genome that was mapped by uniquely mapped read, we used the previously stated settings with nonspecific reads being ignored.
SNP detection
We used the SNP detection function in CLC Genomics Workbench, based on the neighborhood quality standard (NQS) algorithm, using the following parameters: minimum coverage = 5; minimum central base quality = 30; average base quality over a window length of 11 nucleotides = 15; and, minimum allele frequency = 35%. SNPs were filtered by removing those within 10 bp of another variant (both INDELs and SNPs). Also, SNPs located within the pseudo-autosomal region (PAR) were removed [28]. The remaining SNPs were filtered by read-depth to create SNP analyses with minimum depths of 5, 6, 7, 8, 9, and 10.
INDEL detection
We used the deletion and insertion polymorphism (DIP) function in CLC Genomics Workbench using the following parameters: minimum coverage = 5; minimum allele frequency = 35%; and, maximum expected variations = 2. INDELs were filtered in a similar manner as SNPs, with those near other variants and within the PAR being removed. Analyses were conducted at both 5X and 10X read-depths.
CNV detection
We used the Control-FREE copy number algorithm (FREEC) program to identify CNVs in the mapped sequence data [16]. We optimized the program for our data using a break-point threshold of -0.0013 and a coefficient of variation of 0.045. Additionally, we removed any CNVs located within 1 Mb of the beginning and end of all chromosomes in order to reduce the potential for erroneous calling.
Genetic variant annotation and analysis
We re-annotated the variant calls from CLC Genomics using both Galaxy http://galaxy.psu.edu/ and ANNOVAR software programs [29,30]. The join and merge functions in Galaxy were used to annotate the SNPs and INDELs. We used these functions to compare the SNPs to all known SNPs in dbSNP as well as those identified in the reference horse. The ANNOVAR program used the Ensembl annotation database to create an mRNA library, allowing for the determination of amino acid changes. The program was used to determine the locations of all variants within the genome. We also annotated the known SNPs from the reference Thoroughbred genome to determine amino acid changes as well as gene locations for a comparison to the Quarter Horse. The SNPs were divided into groups based on radical and conservative amino acid changes, where radical SNPs result in a difference in polarity or charge when the amino acid was changed while conservative SNPs cause no change in polarity or charge.
The gene lists for each group were converted to human Ensembl gene IDs using Ensembl Biomart. Biological function analysis was performed through the DAVID Functional Annotation Tool, using the default settings [31]. The resulting biological process terms were further grouped by similarities in function to determine enrichment for specific biological processes (Additional file 37, Table S29). Statistical significance (p-value) for each enriched group was determined using Fisher's combined probability test with the p-value created from the DAVID Functional Annotation Tool. The enriched genes for each biological process group were compared between the Quarter Horse and reference Thoroughbred using the Fisher's exact test for a 2 × 2 contingency table.
Analysis of de novo assembly of unmapped reads
We created a de novo assembly of all reads that did not align to the reference genome (including ChrMt and ChrUn) using de Bruijn Graphs (CLC Genomics). We were able to assemble 8,186,040 of the 12,657,236 unmapped reads into contigs with minimum and average contig lengths of 200 bp and 537 bp, respectively. We performed de novo assembly of reads not mapping to either assembled or unassembled chromosomes using the following parameters: similarity = 0.8; length fraction = 0.5; insertion cost = 3; deletion cost = 3; mismatch cost = 3; minimum paired distance = 180 bp; and, maximum paired distance = 500 bp. The resulting 35,540 contigs were further analyzed by BLAST, mapping all complete genomes and chromosomes from RefSeq.
Sequence data and genetic variants
The Illumina FASTQ data generated from this study has been submitted to the NCBI Sequence Read Archive http://www.ncbi.nlm.nih.gov/sra under accession number SRX110702. The mapped sequences and genetic variants have been submitted to Intrepid Bioinformatics and are available for public viewing http://server1.intrepidbio.com/FeatureBrowser/ngsdatasetrecord/record?ngsrecord=6197673305. A complete list of SNPs, INDELs, and CNVs are listed in Additional files 5, 6, and 24, Tables S3, S4 and S24.
Competing interests
The authors declare that they have no competing interests.
Authors' contributions
RD generated the sequencing libraries and analyzed the data. NG and CJ sequenced the libraries and analyzed the data. RD, NDC, JS and SVD conceived the study and wrote the manuscript. All authors read and approved the final manuscript.
Supplementary Material
Additional file 1
Distribution of paired-end sequencing library insert length. Figure indicating distribution of library insert length and number of sequence reads generated per length.
Click here for file(437K, TIFF)
Additional file 2
Sequencing reaction statistics. Table of raw summary statistics of each lane of Illumina sequencing.
Click here for file(11K, XLSX)
Additional file 3
Detailed overview of mapping results. Table of total sequence reads (paired and single end), the number of mapped reads and bp to assembled and unassembled chromosomes.
Click here for file(12K, XLSX)
Additional file 4
Average depth of coverage of assembled chromosomes. Figure of average mapping read coverage of each chromosome.
Click here for file(1.0M, TIFF)
Additional file 5
5X filtered SNPs and INDELS identified by whole genome sequencing. Table listing chromosomal coordinates of SNPs and INDELs identified using a minimum sequence read depth coverage of 5X.
Click here for file(19M, ZIP)
Additional file 6
10X filtered SNPs and INDELS identified by whole genome sequencing. Table listing chromosomal coordinates of SNPs and INDELs identified using a minimum sequence read depth coverage of 10X.
Click here for file(14M, ZIP)
Additional file 7
Loss of SNPs by post-filtering. Table describing the number of SNPs and INDELs removed by post-filtering.
Click here for file(8.4K, XLSX)
Additional file 8
Comparison of minimum depth of coverage SNP filters. Total number of (A) SNPs and (B) non-synonymous SNP remaining after increasing the minimum sequence read depth of coverage. Proportion of homozygous and heterozygous (C) SNPs and (D) non-synonymous SNPs lost by increasing minimum sequence read depth of coverage.
Click here for file(568K, TIFF)
Additional file 9
Annotation of SNPs with minimum read depth coverage of 5X. Table describing the genic annotation of SNPs identified at 5X coverage.
Click here for file(20K, XLSX)
Additional file 10
Annotation of SNPs with minimum read depth coverage of 5X relative to the Ensembl gene annotation. Table describing the intergenic and genic location of SNPs identified at 5X coverage.
Click here for file(21K, XLSX)
Additional file 11
Genes affected by SNPs with minimum read depth coverage of 5X. Table describing the number of genes identified with SNPs, including the annotation of the SNPs within the genetic elements.
Click here for file(13K, XLSX)
Additional file 12
Overview of SNP annotation with minimum read depth coverage of 5X. Table summarizing the number and annotation of SNPs with minimum read depth coverage of 5X.
Click here for file(9.8K, XLSX)
Additional file 13
Annotation of SNPs with minimum read depth coverage of 6X. Table describing the genic annotation of SNPs identified with minimum read depth coverage of 6X coverage.
Click here for file(20K, XLSX)
Additional file 14
Annotation of SNPs with minimum read depth coverage of 6X relative to the Ensembl gene annotation. Table describing the intergenic and genic location of SNPs identified with minimum read depth coverage of 6X coverage.
Click here for file(21K, XLSX)
Additional file 15
Annotation of SNPs with minimum read depth coverage of 7X. Table describing the genic annotation of SNPs identified with minimum read depth coverage of 7X coverage.
Click here for file(20K, XLSX)
Additional file 16
Annotation of SNPs with minimum read depth coverage of 7X relative to the Ensembl gene annotation. Table describing the intergenic and genic location of SNPs identified with minimum read depth coverage of 7X coverage.
Click here for file(21K, XLSX)
Additional file 17
Annotation of SNPs with minimum read depth coverage of 8X. Table describing the genic annotation of SNPs identified with minimum read depth coverage of 8X coverage.
Click here for file(20K, XLSX)
Additional file 18
Annotation of SNPs with minimum read depth coverage of 8X relative to the Ensembl gene annotation. Table describing the intergenic and genic location of SNPs identified with minimum read depth coverage of 8X coverage.
Click here for file(21K, XLSX)
Additional file 19
Annotation of SNPs with minimum read depth coverage of 9X. Table describing the genic annotation of SNPs identified with minimum read depth coverage of 9X coverage.
Click here for file(19K, XLSX)
Additional file 20
Annotation of SNPs with minimum read depth coverage of 9X relative to the Ensembl gene annotation. Table describing the intergenic and genic location of SNPs identified with minimum read depth coverage of 9X coverage.
Click here for file(20K, XLSX)
Additional file 21
Annotation of SNPs with minimum read depth coverage of 10X. Table describing the genic annotation of SNPs identified with minimum read depth coverage of 10X coverage.
Click here for file(19K, XLSX)
Additional file 22
Annotation of SNPs with minimum read depth coverage of 10X relative to the Ensembl gene annotation. Table describing the intergenic and genic location of SNPs identified with minimum read depth coverage of 10X coverage.
Click here for file(20K, XLSX)
Additional file 23
Genes affected by SNPs with minimum read depth coverage of 10X. Table describing the number of genes identified with SNPs, including the annotation of the SNPs within the genetic elements.
Click here for file(13K, XLSX)
Additional file 24
Imprinted genes containing SNPs in coding regions. List of known imprinted genes in humans containing coding SNPs in Quarter Horse.
Click here for file(22K, XLSX)
Additional file 25
Annotation of INDELs with minimum read depth coverage of 5X.
Click here for file(24K, XLSX)
Additional file 26
Annotation of INDELs with minimum read depth coverage of 10X. Table describing genic annotation and location of INDELs with minimum read depth coverage of 10X.
Click here for file(22K, XLSX)
Additional file 27
Distribution of total number of INDELs by length. Plot showing the number of INDELS by length identified with minimum read depth coverage of 10X.
Click here for file(784K, TIFF)
Additional file 28
Coordinates of CNVs. Table describing the location of CNVs identified by analysis of read depth.
Click here for file(51K, XLSX)
Additional file 29
Identification of a loss within an olfactory gene cluster by sequencing read-depth. Plot of read depth indicating a copy number loss within olfactory gene cluster.
Click here for file(302K, TIFF)
Additional file 30
Identification of a homozygous deletion by read-depth coverage. Plot of read depth indicating a homozygous deletion.
Click here for file(235K, TIFF)
Additional file 31
Biological process enrichments and associated p-values for SNPs, INDELS, and CNVs with minimum read depth coverage of 10X. Statistical analyses of biological process enrichments for SNP, INDELS, and CNVs with minimum read depth coverage of 10X.
Click here for file(11K, XLSX)
Additional file 32
Biological process enrichment analysis of conserved and radical non-synonymous SNPs. Statistical analyses of biological process enrichments for conserved and radical SNPs with minimum read depth coverage of 10X.
Click here for file(953K, TIFF)
Additional file 33
Annotation of SNPs in the reference Thoroughbred genome. Table describing genic annotation and location of known SNPs in the reference Thoroughbred genome.
Click here for file(12K, XLSX)
Additional file 34
Comparison of biological processes affected by SNPs in the Quarter Horse and Thoroughbred genomes. (A) BP enrichment of non-synonymous homozygous SNPs and (B) heterozygous SNPs in the Quarter Horse genome. (C) BP enrichment of non-synonymous heterozygous SNPs in the Thoroughbred genome. (D) Number of genes containing non-synonymous SNPs within each BP pathway.
Click here for file(1.4M, TIFF)
Additional file 35
Biological process enrichments and associated p-values for SNP comparison between Quarter Horse and reference Thoroughbred with minimum coverage of 10 reads. Comparison and statistical analysis of biological process enrichments between Quarter Horse and reference Thoroughbred.
Click here for file(10K, XLSX)
Additional file 36
Statistical analysis of enriched genes in the Quarter Horse and reference Thoroughbred. Statistical analysis of comparison between enriched genes in the Quarter Horse and reference Thoroughbred.
Click here for file(9.6K, XLSX)
Additional file 37
Biological process terms and classification. All terms used for biological processes enrichment and their associated category.
Click here for file(30K, XLSX)
Acknowledgements
We thank Sam Wigington at Texas A&M University, College of Veterinary Medicine and Biomedical Sciences for identification of the mare whose genome we sequenced. Funding for this study was provided by Texas AgriLife Research; the GW & G Pool Large Animal Hospital Endowment at the Department of Large Animal Clinical Sciences, College of Veterinary Medicine and Biomedical Sciences; the Department of Animal Science; and the Link Equine Research Endowment, Texas A&M University. This publication is based in part on work supported by Award No. KUS-C1-016-04, made by King Abdullah University of Science and Technology (KAUST).
References
  • Wade CM, Giulotto E, Sigurdsson S, Zoli M, Gnerre S, Imsland F, Lear TL, Adelson DL, Bailey E, Bellone RR, Blocker H, Distl O, Edgar RC, Garber M, Leeb T, Mauceli E, MacLeod JN, Penedo MC, Raison JM, Sharpe T, Vogel J, Andersson L, Antczak DF, Biagi T, Binns MM, Chowdhary BP, Coleman SJ, Della Valle G, Fryc S, Guerin G. et al. Genome sequence, comparative analysis, and population genetics of the domestic horse. Science. 2009;326(5954):865–867. doi: 10.1126/science.1178158. [PubMed] [Cross Ref]
  • http://www.aqha.com/About.aspx
  • 2010 Annual Report: American Quarter Horse Association. In
  • Ridgeway W. The origin and influence of the Thoroughbred horse. New York, NY: Benjamin Blom Inc; 1972.
  • Gu J, Orr N, Park SD, Katz LM, Sulimova G, MacHugh DE, Hill EW. A genome scan for positive selection in thoroughbred horses. PloS one. 2009;4(6):e5767. doi: 10.1371/journal.pone.0005767. [PMC free article] [PubMed] [Cross Ref]
  • Hill EW, Gu J, McGivney BA, MacHugh DE. Targets of selection in the Thoroughbred genome contain exercise-relevant gene SNPs associated with elite racecourse performance. Anim Genet. 2010;41(Suppl 2):56–63. [PubMed]
  • McCue ME, Valberg SJ, Miller MB, Wade C, DiMauro S, Akman HO, Mickelson JR. Glycogen synthase (GYS1) mutation causes a novel skeletal muscle glycogenosis. Genomics. 2008;91(5):458–466. doi: 10.1016/j.ygeno.2008.01.011. [PMC free article] [PubMed] [Cross Ref]
  • Tryon RC, Penedo MC, McCue ME, Valberg SJ, Mickelson JR, Famula TR, Wagner ML, Jackson M, Hamilton MJ, Nooteboom S, Bannasch DL. Evaluation of allele frequencies of inherited disease genes in subgroups of American Quarter Horses. J Am Vet Med Assoc. 2009;234(1):120–125. doi: 10.2460/javma.234.1.120. [PubMed] [Cross Ref]
  • Rudolph JA, Spier SJ, Byrns G, Rojas CV, Bernoco D, Hoffman EP. Periodic paralysis in quarter horses: a sodium channel mutation disseminated by selective breeding. Nat Genet. 1992;2(2):144–147. doi: 10.1038/ng1092-144. [PubMed] [Cross Ref]
  • Ward TL, Valberg SJ, Adelson DL, Abbey CA, Binns MM, Mickelson JR. Glycogen branching enzyme (GBE1) mutation causing equine glycogen storage disease IV. Mammalian genome: official journal of the International Mammalian Genome Society. 2004;15(7):570–577. [PubMed]
  • Finno CJ, Spier SJ, Valberg SJ. Equine diseases caused by known genetic mutations. Vet J. 2009;179(3):336–347. doi: 10.1016/j.tvjl.2008.03.016. [PubMed] [Cross Ref]
  • Tryon RC, White SD, Bannasch DL. Homozygosity mapping approach identifies a missense mutation in equine cyclophilin B (PPIB) associated with HERDA in the American Quarter Horse. Genomics. 2007;90(1):93–102. doi: 10.1016/j.ygeno.2007.03.009. [PubMed] [Cross Ref]
  • Panizzi L, Barber SM, Lang HM, Carmalt JL. Carpometacarpal osteoarthritis in thirty-three horses. Vet Surg. 2009;38(8):998–1005. doi: 10.1111/j.1532-950X.2009.00589.x. [PubMed] [Cross Ref]
  • Brooks SA, Makvandi-Nejad S, Chu E, Allen JJ, Streeter C, Gu E, McCleery B, Murphy BA, Bellone R, Sutter NB. Morphological variation in the horse: defining complex traits of body size and shape. Anim Genet. pp. 159–165. [PubMed]
  • Gregg C, Zhang J, Weissbourd B, Luo S, Schroth GP, Haig D, Dulac C. High-resolution analysis of parent-of-origin allelic expression in the mouse brain. Science. 2010;329(5992):643–648. doi: 10.1126/science.1190830. [PMC free article] [PubMed] [Cross Ref]
  • Boeva V, Zinovyev A, Bleakley K, Vert JP, Janoueix-Lerosey I, Delattre O, Barillot E. Control-free calling of copy number alterations in deep-sequencing data using GC-content normalization. Bioinformatics. pp. 268–269. [PMC free article] [PubMed]
  • Wagner HJ, Reissmann M. New polymorphism detected in the horse MC1R gene. Anim Genet. 2000;31(4):289–290. doi: 10.1046/j.1365-2052.2000.00655.x. [PubMed] [Cross Ref]
  • Hill EW, McGivney BA, Gu J, Whiston R, Machugh DE. A genome-wide SNP-association study confirms a sequence variant (g.66493737C > T) in the equine myostatin (MSTN) gene as the most powerful predictor of optimum racing distance for Thoroughbred racehorses. BMC Genomics. p. 552. [PMC free article] [PubMed]
  • Gu J, MacHugh DE, McGivney BA, Park SD, Katz LM, Hill EW. Association of sequence variants in CKM (creatine kinase, muscle) and COX4I2 (cytochrome c oxidase, subunit 4, isoform 2) genes with racing performance in Thoroughbred horses. Equine Vet J. pp. 569–575. [PubMed]
  • Ajay SS, Parker SC, Abaan HO, Fajardo KV, Margulies EH. Accurate and comprehensive sequencing of personal genomes. Genome Res. 2011;21(9):1498–1505. doi: 10.1101/gr.123638.111. [PubMed] [Cross Ref]
  • McCarroll SA, Hadnott TN, Perry GH, Sabeti PC, Zody MC, Barrett JC, Dallaire S, Gabriel SB, Lee C, Daly MJ, Altshuler DM. Common deletion polymorphisms in the human genome. Nat Genet. 2006;38(1):86–92. doi: 10.1038/ng1696. [PubMed] [Cross Ref]
  • Egan CM, Sridhar S, Wigler M, Hall IM. Recurrent DNA copy number variation in the laboratory mouse. Nat Genet. 2007;39(11):1384–1389. doi: 10.1038/ng.2007.19. [PubMed] [Cross Ref]
  • Graubert TA, Cahan P, Edwin D, Selzer RR, Richmond TA, Eis PS, Shannon WD, Li X, McLeod HL, Cheverud JM, Ley TJ. A high-resolution map of segmental DNA copy number variation in the mouse genome. PLoS Genet. 2007;3(1):e3. doi: 10.1371/journal.pgen.0030003. [PubMed] [Cross Ref]
  • Pinto D, Marshall C, Feuk L, Scherer SW. Copy-number variation in control population cohorts. Hum Mol Genet. 2007;16(Spec No. 2):R168–173. [PubMed]
  • Springer NM, Ying K, Fu Y, Ji T, Yeh CT, Jia Y, Wu W, Richmond T, Kitzman J, Rosenbaum H, Iniguez AL, Barbazuk WB, Jeddeloh JA, Nettleton D, Schnable PS. Maize inbreds exhibit high levels of copy number variation (CNV) and presence/absence variation (PAV) in genome content. PLoS Genet. 2009;5(11):e1000734. doi: 10.1371/journal.pgen.1000734. [PMC free article] [PubMed] [Cross Ref]
  • Oki H, Kusunose R, Nakaoka H, Nishiura A, Miyake T, Sasaki Y. Estimation of heritability and genetic correlation for behavioural responses by Gibbs sampling in the Thoroughbred racehorse. J Anim Breed Genet. 2007;124(4):185–191. doi: 10.1111/j.1439-0388.2007.00659.x. [PubMed] [Cross Ref]
  • Momozawa Y, Takeuchi Y, Kusunose R, Kikusui T, Mori Y. Association between equine temperament and polymorphisms in dopamine D4 receptor gene. Mamm Genome. 2005;16(7):538–544. doi: 10.1007/s00335-005-0021-3. [PubMed] [Cross Ref]
  • Raudsepp T, Chowdhary BP. The horse pseudoautosomal region (PAR): characterization and comparison with the human, chimp and mouse PARs. Cytogenetic and genome research. 2008;121(2):102–109. doi: 10.1159/000125835. [PubMed] [Cross Ref]
  • Goecks J, Nekrutenko A, Taylor J. Galaxy: a comprehensive approach for supporting accessible, reproducible, and transparent computational research in the life sciences. Genome Biol. p. R86. [PMC free article] [PubMed]
  • Wang K, Li M, Hakonarson H. ANNOVAR: functional annotation of genetic variants from high-throughput sequencing data. Nucleic Acids Res. p. e164. [PMC free article] [PubMed]
  • Huang da W, Sherman BT, Zheng X, Yang J, Imamichi T, Stephens R, Lempicki RA. Extracting biological meaning from large gene lists with DAVID. Curr Protoc Bioinformatics. 2009;Chapter 13:Unit 13 11. [PubMed]
  • Bellone RR, Forsyth G, Leeb T, Archer S, Sigurdsson S, Imsland F, Mauceli E, Engensteiner M, Bailey E, Sandmeyer L, Grahn B, Lindblad-Toh K, Wade CM. Fine-mapping and mutation analysis of TRPM1: a candidate gene for leopard complex (LP) spotting and congenital stationary night blindness in horses. Briefings in functional genomics. 2010;9(3):193–207. doi: 10.1093/bfgp/elq002. [PubMed] [Cross Ref]
  • Brooks SA, Gabreski N, Miller D, Brisbin A, Brown HE, Streeter C, Mezey J, Cook D, Antczak DF. Whole-genome SNP association in the horse: identification of a deletion in myosin Va responsible for Lavender Foal Syndrome. PLoS Genet. 2010;6(4):e1000909. doi: 10.1371/journal.pgen.1000909. [PMC free article] [PubMed] [Cross Ref]
  • Marklund L, Moller MJ, Sandberg K, Andersson L. A missense mutation in the gene for melanocyte-stimulating hormone receptor (MC1R) is associated with the chestnut coat color in horses. Mammalian genome: official journal of the International Mammalian Genome Society. 1996;7(12):895–899. doi: 10.1007/s003359900264. [PubMed] [Cross Ref]
  • Brooks SA, Bailey E. Exon skipping in the KIT gene causes a Sabino spotting pattern in horses. Mammalian genome: official journal of the International Mammalian Genome Society. 2005;16(11):893–902. doi: 10.1007/s00335-005-2472-y. [PubMed] [Cross Ref]
  • Brooks SA, Lear TL, Adelson DL, Bailey E. A chromosome inversion near the KIT gene and the Tobiano spotting pattern in horses. Cytogenetic and genome research. 2007;119(3-4):225–230. doi: 10.1159/000112065. [PubMed] [Cross Ref]
  • Spirito F, Charlesworth A, Linder K, Ortonne JP, Baird J, Meneguzzi G. Animal models for skin blistering conditions: absence of laminin 5 causes hereditary junctional mechanobullous disease in the Belgian horse. The Journal of investigative dermatology. 2002;119(3):684–691. doi: 10.1046/j.1523-1747.2002.01852.x. [PubMed] [Cross Ref]
  • Brunberg E, Andersson L, Cothran G, Sandberg K, Mikko S, Lindgren G. A missense mutation in PMEL17 is associated with the Silver coat color in the horse. BMC Genet. 2006;7:46. [PMC free article] [PubMed]
  • Graves KT, Henney PJ, Ennis RB. Partial deletion of the LAMA3 gene is responsible for hereditary junctional epidermolysis bullosa in the American Saddlebred Horse. Anim Genet. 2009;40(1):35–41. doi: 10.1111/j.1365-2052.2008.01795.x. [PubMed] [Cross Ref]
  • Shin EK, Perryman LE, Meek K. A kinase-negative mutation of DNA-PK(CS) in equine SCID results in defective coding and signal joint formation. J Immunol. 1997;158(8):3565–3569. [PubMed]
  • Aleman M, Riehl J, Aldridge BM, Lecouteur RA, Stott JL, Pessah IN. Association of a mutation in the ryanodine receptor 1 gene with equine malignant hyperthermia. Muscle & nerve. 2004;30(3):356–365. doi: 10.1002/mus.20084. [PubMed] [Cross Ref]
  • Cannon SC, Hayward LJ, Beech J, Brown RH Jr. Sodium channel inactivation is impaired in equine hyperkalemic periodic paralysis. Journal of neurophysiology. 1995;73(5):1892–1899. [PubMed]
  • Cook D, Brooks S, Bellone R, Bailey E. Missense mutation in exon 2 of SLC36A1 responsible for champagne dilution in horses. PLoS Genet. 2008;4(9):e1000195. doi: 10.1371/journal.pgen.1000195. [PMC free article] [PubMed] [Cross Ref]
  • Yang GC, Croaker D, Zhang AL, Manglick P, Cartmill T, Cass D. A dinucleotide mutation in the endothelin-B receptor gene is associated with lethal white foal syndrome (LWFS); a horse variant of Hirschsprung disease. Hum Mol Genet. 1998;7(6):1047–1052. doi: 10.1093/hmg/7.6.1047. [PubMed] [Cross Ref]
  • Mariat D, Taourit S, Guerin G. A mutation in the MATP gene causes the cream coat colour in the horse. Genetics, selection, evolution: GSE. 2003;35(1):119–133. doi: 10.1186/1297-9686-35-1-119. [PMC free article] [PubMed] [Cross Ref]
  • Rieder S, Taourit S, Mariat D, Langlois B, Guerin G. Mutations in the agouti (ASIP), the extension (MC1R), and the brown (TYRP1) loci and their association to coat color phenotypes in horses (Equus caballus) Mammalian genome: official journal of the International Mammalian Genome Society. 2001;12(6):450–455. doi: 10.1007/s003350020017. [PubMed] [Cross Ref]
  • Rosengren Pielberg G, Golovko A, Sundstrom E, Curik I, Lennartsson J, Seltenhammer MH, Druml T, Binns M, Fitzsimmons C, Lindgren G, Sandberg K, Baumung R, Vetterlein M, Stromberg S, Grabherr M, Wade C, Lindblad-Toh K, Ponten F, Heldin CH, Solkner J, Andersson L. A cis-acting regulatory mutation causes premature hair graying and susceptibility to melanoma in the horse. Nat Genet. 2008;40(8):1004–1009. doi: 10.1038/ng.185. [PubMed] [Cross Ref]
Articles from BMC Genomics are provided here courtesy of
BioMed Central