Whole-genome resequencing of Hanwoo (Korean cattle) and insight into regions of homozygosity
- Kyung-Tai Lee†1,
- Won-Hyong Chung†2,
- Sung-Yeoun Lee†3,
- Jung-Woo Choi†4,
- Jiwoong Kim2,
- Dajeong Lim1,
- Seunghwan Lee5,
- Gul-Won Jang1,
- Bumsoo Kim1,
- Yun Ho Choy6,
- Xiaoping Liao4,
- Paul Stothard4,
- Stephen S Moore4, 7,
- Sang-Heon Lee2, 8,
- Sungmin Ahn3Email author,
- Namshin Kim2, 8Email author and
- Tae-Hun Kim1Email author
© Lee et al.; licensee BioMed Central Ltd. 2013
Received: 14 September 2012
Accepted: 18 June 2013
Published: 30 July 2013
Hanwoo (Korean cattle), which originated from natural crossbreeding between taurine and zebu cattle, migrated to the Korean peninsula through North China. Hanwoo were raised as draft animals until the 1970s without the introduction of foreign germplasm. Since 1979, Hanwoo has been bred as beef cattle. Genetic variation was analyzed by whole-genome deep resequencing of a Hanwoo bull. The Hanwoo genome was compared to that of two other breeds, Black Angus and Holstein, and genes within regions of homozygosity were investigated to elucidate the genetic and genomic characteristics of Hanwoo.
The Hanwoo bull genome was sequenced to 45.6-fold coverage using the ABI SOLiD system. In total, 4.7 million single-nucleotide polymorphisms and 0.4 million small indels were identified by comparison with the Btau4.0 reference assembly. Of the total number of SNPs and indels, 58% and 87%, respectively, were novel. The overall genotype concordance between the SNPs and BovineSNP50 BeadChip data was 96.4%. Of 1.6 million genetic differences in Hanwoo, approximately 25,000 non-synonymous SNPs, splice-site variants, and coding indels (NS/SS/Is) were detected in 8,360 genes. Among 1,045 genes containing reliable specific NS/SS/Is in Hanwoo, 109 genes contained more than one novel damaging NS/SS/I. Of the genes containing NS/SS/Is, 610 genes were assigned as trait-associated genes. Moreover, 16, 78, and 51 regions of homozygosity (ROHs) were detected in Hanwoo, Black Angus, and Holstein, respectively. ‘Regulation of actin filament length’ was revealed as a significant gene ontology term and 25 trait-associated genes for meat quality and disease resistance were found in 753 genes that resided in the ROHs of Hanwoo. In Hanwoo, 43 genes were located in common ROHs between whole-genome resequencing and SNP chips in BTA2, 10, and 13 coincided with quantitative trait loci for meat fat traits. In addition, the common ROHs in BTA2 and 16 were in agreement between Hanwoo and Black Angus.
We identified 4.7 million SNPs and 0.4 million small indels by whole-genome resequencing of a Hanwoo bull. Approximately 25,000 non-synonymous SNPs, splice-site variants, and coding indels (NS/SS/Is) were detected in 8,360 genes. Additionally, we found 25 trait-associated genes for meat quality and disease resistance among 753 genes that resided in the ROHs of Hanwoo. These findings will provide useful genomic information for identifying genes or casual mutations associated with economically important traits in cattle.
KeywordsHanwoo Resequencing NS/SS/I ROH
The bovine genome was one of the first mammalian genomes sequenced, likely because cattle are important farm animals serving as major nutritional sources for humans and because of their evolutionary position as a representative of the Ruminantia, a phylogenetically distant clade to humans and rodents . The bovine genome sequencing consortium sequenced a single inbred female Hereford cow and her sire using a combination of hierarchical sequencing and whole-genome shotgun sequencing ; the data have been assembled into two reference genomes, Btau and UMD [3, 4].
After the bovine reference genome was assembled, several bovine genomes were resequenced, providing more insight into the genetic diversity of cattle that may be associated with phenotypic differences between breeds. In 2008, Van Tassell et al. reported more than 60,000 putative single-nucleotide polymorphisms (SNPs) identified from a reduced representation DNA library of 66 cattle representing three populations . In 2009, Eck et al. performed the first single cattle whole-genome resequencing and reported more than 2 million novel SNPs in a Fleckvieh bull . In 2011, Kawahara-Miki et al. resequenced the genome of a single Kuchinoshima-Ushi bull, a Japanese native cattle breed whose lineage has been strictly maintained in a small island secluded from mainland Japan . In that study, more than 5.5 million novel SNPs were reported, and the Kuchinoshima-Ushi bull was determined to be genetically distinct from European domestic cattle breeds. Most recently, Stothard and colleagues reported whole-genome resequencing of Black Angus and Holstein, representative beef and dairy breeds, respectively, in North America, leading to the identification of substantial numbers of SNPs and copy number variants (CNVs) that could potentially be used as genetic markers across the genome .
When high-density genome-wide SNP data are available, analyses can identify genetic differences between similar populations. Understanding the genetic mechanisms leading to phenotypic differentiation requires identification of the genomic regions that have been under artificial selection in cattle breeds. For example, strong artificial selection will increase the frequency of favorable alleles at loci affecting meat quality traits in meat-producing breeds such as Hanwoo or Black Angus. In this process, a small region of the genome surrounding the mutations is also selected, resulting in a small genome region that shows reduced variation. Many methods have been developed for the detection of selection signatures from genome analyses, such as the use of regions of homozygosity (ROHs) , the integrated haplotype score (iHS) , FST, and the extended haplotype homozygosity (EHH) statistic , according to the detection of the timescale for selection signatures. ROHs are without heterozygosity in the diploid states and provide association evidence at the genome-wide scale for complex traits.
Hanwoo, a Korean cattle breed, is reported to have originated from crossbreeding between taurine and zebu cattle and migrated to the Korean peninsula through North China; their history as a draft animal dates back at least 5,000 years [13, 14]. Afterward, Hanwoo was maintained without the introduction of additional germplasm. Hanwoo was raised as a draft animal until the 1970s. In the late 1970s, the Korean government initiated a Hanwoo genetic breeding program to improve meat quantity and quality.
In this study, we sequenced the genome of a Hanwoo breeding bull and identified single nucleotide polymorphisms (SNPs) based on the Bos taurus reference genome assembly (Btau4.0). SNPs of Hanwoo were compared with those of Black Angus and Holstein. Moreover, functional annotation was carried out for SNPs. We also investigated genomic regions of homozygosity in Hanwoo, Black Angus, and Holstein.
Results and discussion
Genome sequencing, SNP/indel detection, and genotype concordance
Summary of the sequenced reads for Hanwoo, Black Angus, and Holstein
Number of total reads
Number of mapped reads
Total bases of mapped reads
Genome coverage (%)
In total, 4,781,758 SNPs were identified in the Hanwoo genome using the Genome Analysis Tool Kit (GATK) 1.0.5974 [15, 16]. Among them, 2,327,616 SNPs (48.8%) were found in the single-nucleotide polymorphism database (dbSNP, build 133) while the remaining 2,454,142 SNPs (51.2%) were novel; 3,104,888 (64.9%) were heterozygous and 1,676,870 (35.1%) were homozygous, with a ratio of 1:1.85 (homozygous:heterozygous) (Additional file 3). Using UnifiedGenotyper in GATK, we identified 391,512 small indels (−14 to +22bp); 228,121 (58.3%) were heterozygous and 163,391 (41.7%) were homozygous (160,316 insertions and 231,196 deletions). Of these indels, 49,225 were found in dbSNP (build 133) while the remaining 342,287 indels (87.4%) were novel. All SNPs and indels identified in Hanwoo were submitted to the dbSNP at NCBI under the handle NIAS_AGBSGL.
To evaluate the SNP calling from our high-throughput genome sequencing data, concordance analysis was performed between Hanwoo genome resequencing and the SNP chip data. The same genomic DNA from Hanwoo used for deep resequencing was genotyped for 54,001 SNPs using BovineSNP50 BeadChip (Illumina). All probe sequences were mapped against the Btau4.0 reference genome assembly, and 50,411 positions were identified as unique genomic loci. In total, 1,061 (2.8%) of 38,049 homozygous calls by the SNP chip have been identified as heterozygous by NGS. In total, 526 (4.3%) of 12,362 heterozygous calls by the SNP chip were identified as homozygous by NGS (Additional file 4). The overall genotype concordance was 96.2%. The non-reference sensitivity and non-reference discrepancy rates were 97.1% and 7.0%, respectively. Non-reference sensitivity is the fraction of sites called variants (A/B or B/B) in comparison to those that are also called variants in evaluation data. The non-reference discrepancy rate, which is a good measure for testing the accuracy of genotype calls, can show the accuracy of genotype calling at sites called by both sites by excluding concordant genotypes (http://gatkforums.broadinstitute.org/discussion/48/using-varianteval).
Functional annotation of genomic variation
Next, we investigated whether NS/SS/I-containing genes were associated with economic traits and then categorized them into meat, disease resistance, growth, milk, and fecundity. We used previously reported information on trait-associated genes [7, 20, 21]. In total, 619 genes were assigned as trait-associated genes: 464 genes for meat quality, 144 genes for disease resistance, 25 genes for milk production, 8 genes for fecundity, and 6 genes for growth rate (Additional file 5). Of the 464 genes for meat quality, 228 contained more than one NS/SS/I. The titin (TTN) gene has 62 NS/SS/Is, the largest number among the genes related to meat quality. The bovine major histocompatibility complex (MHC) class I heavy chain isoform 1 precursor (BOLA, ENSBTAG00000002069) gene, which contain 32 NS/SS/Is, has the largest number of NS/SS/Is among the genes related to disease resistance (Additional file 5). The higher number of NS/SS/Is in TTN than BOLA may be due to the difference of gene size; 274.866 kb and 3.788 kb of TTN and BOLA genes, respectively. The TTN gene encodes the titin protein, the largest protein, which consists of 317 exons in 274.866 kb of genomic DNA in BTA2 (Ensembl database UMD3.1). The TTN gene plays a role in myofibrillogenesis and is associated with marbling . Moreover, the 231054C>T variant within the promoter region of TTN is associated with a marbling trait and is differentially expressed between high- and low-marbling muscle samples . However, many NS/SS/Is likely affect the function of titin, which acts as a molecular spring for the passive elasticity of muscles . These NS/SS/Is within the TTN gene may be informative variants for understanding the effects of steric changes in the TTN protein. Of the 144 genes for disease resistance, 74 also contained NS/SS/Is, and many novel damaging NS/SS/Is were detected in several genes including BCL2-like 1 (BCL2L1), nitric oxide synthase 1 (NOS1), nucleotide-binding oligomerization domain-containing protein 2 (NOD2), granzyme A (GZMA), and semaphorin-5A (SEMA5A), as well as the BOLA gene (Additional file 5). Among 109 genes containing more than one novel damaging specific NS/SS/Is in Hanwoo, the BCL2L1, GZMA, and CD5 genes are known as candidate genes for the disease resistance trait (Additional file 6). We suggest that the exonic variation identified in this study will provide valuable information for functional studies as well as marker development associated with economic traits in cattle.
Regions of homozygosity within the three breeds
A ROH is a continuous or uninterrupted stretch of DNA without heterozygosity in the diploid state. A discrepancy has existed in the minimum standard of definition of ROH among the groups that have been studied for ROHs to date . Most previous ROH studies have been performed with SNP chip results, an average of 50 SNPs of 5 Mb in size with an average distance of 100 kb between them, and an allowance of up to 2% heterozygous SNPs within a ROH . However, at present, no standardized criteria have been established for defining ROHs . In this study, using mass genotype data derived from whole-genome resequencing, we shortened the detection window of ROHs and loosened the permissible ratio of heterozygous SNPs (Additional file 7).
Summary of specific regions of homozygosity (ROHs) in Hanwoo, Black Angus, and Holstein
No. of ROHs
ROH size (Mb)
Average ROH size (Mb)
A total of 753 genes resided in the ROHs of Hanwoo, whereas 1,320 and 2,482 genes existed in the ROHs of Black Angus and Holstein, respectively. Among them, 77 and 30 common genes were located within the overlapping ROHs between Hanwoo and Angus and between Hanwoo and Holstein, respectively (Additional file 9). Among the 753 genes in the ROHs of Hanwoo, 505 genes contained no NS/SS/Is. A total of 2,158 genes (10.2%) contained no NS/SS/Is in the ROHs in any of the three breeds (Additional file 5). Moreover, we performed functional enrichment analysis using gene ontology (GO) for genes in the ROHs of the three breeds (Additional file 10). In Hanwoo, one significant GO term was ‘regulation of actin filament length related to muscle metabolism’ (GO:0030832, p-value = 0.044), including actin-related protein 3 homolog (ACTR3), actin-related protein 2/3 complex, subunit 2 (ARPC2), villin 1 (VIL1), and destrin (DSTN) genes. Meat tenderness is generated by the disruption of actin filaments and by breaking down the interaction between the actin and myosin filaments . Notably, a significant GO term of ‘striated muscle cell differentiation’ (GO:0051146, p-value = 0.034) was found in Holstein, including the retinoid X receptor, alpha (RXRA) gene, which inhibits adipogenesis  and plays a negative role in marbling in Hanwoo . Because Hanwoo and Black Angus were bred as beef cattle, 77 genes in overlapping ROHs between the breeds were used to analyze GO and the KEGG pathway. Although eight significant GO terms were detected, most were related to the immune system, such as T cell activation and lymphocyte activation, rather than meat traits. The presence of many immune system-related genes in the identified ROHs could reflect selection (natural or artificial) for disease resistance. According to functional enrichment analysis using KEGG pathway terms, vitamin B6 metabolism (bta00750, p-value = 0.025) was significantly enriched, including the aldehyde oxidase 1 (AOX1) gene in Hanwoo (Additional file 10). Vitamin B6 induces the differentiation of adipocytes from pre-adipocytes and facilitates fat accumulation . In particular, the AOX1 gene is a target of peroxisome proliferator-activated receptors alpha and gamma (PPARα and PPAR γ) as a key gene in adipogenesis . Melanogenesis (bta04916, p-value = 0.009) was also detected in Holstein. Sixteen genes, including tyrosinase (TYR) and the melanocortin 1 receptor (MC1R), exist in this term. TYR is the rate-limiting enzyme in the melanogenesis pathway. Tyrosinase activity is regulated by the MC1R. Recently, Kuhn and Weikard reported the dilution of black pigment (eumelanin) in an F2 Holstein × Charolais population . These genes may also be partially responsible for the coat color pattern of Holsteins. These observations suggest that genes within the ROH accumulate the biological functions for characteristics of each breed during a process of artificial selection.
Trait-associated genes in ROHs of Hanwoo, Black Angus, and Holstein
PDK1, HAT1, TTN, CHRNA1
DBI, CTLA4, CD28, IDH1, ALS2, MYL1, ACADL
PRKAG3, DES, SLC11A1
SP1, RARG, SOCS2
COL9A3, SIRPA, ATRN
CTLA4, CD28, ALS2
SREBF2, FGF6, MCHR1
MT, DR, MK3
CCL3, ACACA, CCL5
CYP1A1, CTSG, GZMH, GZMB, LOXL1
IGFBP2, PRKAG3, SLC11A1, IGFBP5
PRKAG1, WNT1, SLC16A7, IFNG, PFKM, GPD1
CSF2RB, IGF1, MB, PMCH, PVALB, CACNG2
ANXA1, VLDLR, CD72, TPM2, SMARCA2, ALDH1A1
PSMB7, CRAT, NOTCH1, ENG
ITGB1, NEBL, PHYH
MT, DR, MK
MC1R, CYBA, CDH15, FOXC2
SERPINA1, PSMA6, CFL2, SERPINA14
ELN, HSPB1, PLOD3, GPC2, MDH2, SERPINE1
ADRB3, CASP3, FGFR1
In addition, to detect common ROHs between the genome sequence and SNP chip, we calculated the ROHs from 40 Hanwoo bulls as well as 20 Angus and 19 Holstein individuals using the Bovine 50K SNP chip (Figure 2). We identified four, eight, and eight common ROHs between both data sets in Hanwoo, Black Angus, and Holstein, respectively (Additional file 9). In Hanwoo, 43 genes located in common ROHs were shared between genome sequencing and SNP chips in BTA2, 10, and 13 (Figure 2). Of 43 genes, 22 genes contained no NS/SS/I (Additional file 5). Moreover, four common ROHs in Hanwoo coincided with quantitative trait loci (QTLs) for meat fat traits (Figure 2). Specifically, two regions in BTA2 (95.3-96.4 Mb and 100.9-101.4 Mb in Btau4.0) were common ROHs between genome sequencing and SNP chips in Hanwoo and Black Angus. Of the 18 genes that resided in these regions, WD repeat domain 12 (WDR12), amyotrophic lateral sclerosis 2 (juvenile) chromosome region, candidate 8 (ALS2CR8), cytochrome P450, family 20, subfamily A, polypeptide 1 (CYP20A1), and cAMP responsive element binding protein 1 (CREB1) genes belonged to a significant GO term of metabolic processes in Hanwoo. Among them, the CREB1 gene has been shown to be related to fat metabolism. In 2012, Lee et al. reported that the expression of the cAMP responsive element binding protein (CREB1) gene is higher in muscle with high IMF content in Hanwoo . CREB1 is a transcription factor containing a basic leucine zipper. The CREB protein is phosphorylated in response to increased cAMP, allowing it to efficiently interact with the transcriptional co-activator protein, CREB binding protein, to stimulate the transcription of cAMP target genes . Moreover, Casimir and Ntambi reported that intracellular cAMP activates the expression of the stearoyl-CoA desaturase gene, a key enzyme involved in monounsaturated fatty acid synthesis through activation of the CREB protein . In 2009, Wang et al. observed that messenger RNA expression of a lipogenesis-related gene, stearoyl-coA desaturase (SCD), peaked at 20 to 25 months in crosses between Wagyu and Hereford, which was highly correlated with intramuscular fat content in these animals . These findings suggest that elevated CREB expression may stimulate genes involved in the lipid biosynthesis pathway such as SCD  and HMG-Co synthase , resulting in an increase in IMF content within muscles. Also, the ALS2 gene, which is related to meat traits, as well as cytotoxic T-lymphocyte-associated protein 4 (CTLA4) and CD28 molecule (CD28) genes for disease resistance, resided in a common ROH (BTA2: 94.8-96.9 Mb in Btau4.0) in Hanwoo and Black Angus according to genome sequencing (Figure 2). In livestock animals bred by an improvement scheme for economic traits, the use of ROHs will be a good genomic strategy for tracking and planning improvements in breeding.
In this study, we sequenced the whole genome of a Hanwoo bull and newly identified 2,454,142 SNPs and 342,287 small indels by comparison with the Hereford reference genome sequence. We also found 1,663,599 SNPs and 142,297 indels that were located in genic regions of 20,955 genes in the NCBI Reference Sequence Database (RefSeq), of which 22,752 SNPs and 2,163 indels were non-synonymous, frameshift, nonsense, or splice-site SNPs potentially capable of affecting protein functions. This suggests that genes containing several NS/SS/Is may have evolved into multi-copy genes for environmental adaptation, or that NS/SS/Is may be distorted due to an incorrect reference genome sequence. A ROH is a continuous or uninterrupted stretch of DNA without heterozygosity in the diploid state. In this study, we defined 16 ROHs in Hanwoo using a detection window of 400 kb and 20% of heterozygous SNPs using genotype data derived from whole-genome resequencing. The cumulative lengths of ROHs per genome, as well as the number of ROHs in Hanwoo, were smaller than those in Black Angus and Holstein. This suggests that the total lengths of ROHs in Holstein and Black Angus are longer than those of Hanwoo due to a longer period of time for artificial selection in those breeds. In addition, the dispersing pattern of ROHs in chromosomes was different between breeds. We suggest that these patterns would result from the different origins and breeding strategies among these three breeds. Moreover, 753 genes were observed in the ROHs of Hanwoo, of which 25 genes were associated with meat quality and disease resistance traits. In addition, we observed common ROHs between the genome sequence and high-density SNP chip data. This combinatorial ROH survey approach may be another effective method for identifying domestication genes. The findings of this study will provide valuable information for functional studies, as well as for marker development associated with economically important traits in cattle.
We sequenced the genome of a proven Hanwoo bull (27223) obtained from the Hanwoo Experiment Station, National Institute of Animal Science, Rural Development Administration, Korea. Bull 27223 was selected for mapping for its representativeness of the population at the Hanwoo Experiment station. Bull 27223 is a descendent of KPN369, which was one of the most frequently used Hanwoo bulls for artificial insemination in Korea during the early 2010s. Also, bull 27223 was selected for its superiority in growth performance with superior genetic potential in carcass quality. Therefore, many calves born since 2010 have been sired by this bull. The study protocol and standard operating procedures were reviewed and approved by the Institutional Animal Care and Use Committee of the National Institute of Animal Science (Suwon, Republic of Korea).
Whole-genome sequencing library preparation
Genomic DNA (gDNA) was extracted from whole blood with a QIAamp DNA Blood Maxi Kit according to the manufacturer’s instructions (Qiagen). Libraries were prepared according to the SOLiD System Mate-paired Library Preparation protocol of the Applied Biosystems SOLiD System: Library Preparation Guide (02/2009 & 10/2009 editions).
Briefly, gDNA was fragmented using Covaris S2 (Covaris) and HydroShear (Genomic Solutions) at the proper settings for targeted sizes. A QIAquick Gel Extraction Kit (Qiagen) was used for subsequent purification of sheared DNA, enzymatic reactions, and size-selected DNA in agarose gels according to the manufacturer’s instructions. To repair damaged DNA ends and obtain 5′-phosphorylated blunt-ends (5′P), the fragments were end-repaired using the End-It DNA End-Repair Kit (Epicentre Biotechnologies) according to the manufacturer’s instructions. Ligations for the adaptor attachment and circularization were accomplished using the Quick Ligation Kit (New England BioLabs). DNA quantitations were performed using a NanoDrop ND 1000 Spectrophotometer (Thermo Fisher Scientific), except for those followed by library amplification for emulsion PCR (ePCR).
In chronological order, the sheared gDNA fragments were end-repaired and the LMP CAP Adaptors (missing the 5′ phosphate from one oligonucleotide resulting in a nick on each strand when the DNA is circularized at a later step) were ligated to the end-repaired DNA fragments. The adaptor-ligated products were separated on a 1% agarose gel and excised from the gel at the appropriate positions for span size ranges (600–700 bp, 1–2 kb, and 0.6-2.2 kb). Size-selected DNA fragments were circularized with a biotinylated internal adaptor. Uncircularized DNA fragments were eliminated using Plasmid-Safe ATP-Dependent DNase (Epicentre Biotechnologies). Nick translation was performed for 14 min at 0°C in an ice-water bath using Escherichia coli DNA polymerase I with the circularized DNA fragments. The nick-translated products were cleaved at the nicks using T7 exonuclease and S1 nuclease, and end-repaired as described above. P1 and P2 adaptors (used for library amplification, ePCR, and ligation sequencing) were ligated to ends of the end-repaired DNA. Then the ligated DNA underwent nick translation with DNA polymerase I. The completed library was amplified using Library PCR primers 1 and 2 with Cloned Pfu polymerase (Stratagene) or Platinum® PCR Amplification Mix (SOLiD Long Mate-Paired Library Construction Kit, ABI). The amplified library was ran on a 4% agarose gel and the correct-sized band (275–300 bp) was excised and eluted, and quantitated by Qubit IT (Invitrogen). ePCR was carried out according to the Applied Biosystems SOLiD System: Template Bead Preparation Guide. The concentration of each library for ePCR was designed to range from 1.0 to 1.5 pM.
Library sequencing of template beads
Sequencing was performed according to the Applied Biosystems SOLiD System: Instrument Operation Guide. Templated beads were deposited onto two slides and sequencing was carried out to 50 bases using SOLiD v3.0 chemistry, with the exception that the library prepared from 0.6-2.2 kb-sheared DNA fragments was used for four slides and sequencing was carried out to 50 bases using SOLiD v3 plus chemistry.
Short-read alignment, variant calling, and annotation
Paired-end 50 bp reads from Hanwoo, Black Angus, and Holstein were mapped to the Btau4.0 reference genome assembly using BFAST 0.7.0a , with options bfast match “-A 1 -z -K 100 -M 500,” bfast localalign “-A 1 -o 10,” and bfast postprocess “-A 1 -a 3 -Y 2 -z -O 1.” Aligned reads considered to be PCR duplicates were removed using the MarkDuplicates algorithm in Picard tools 1.57. This algorithm identifies the 5′ coordinates and mapping orientations of each read pair by considering gaps and jumps. The reads that mapped to the same position and orientation are marked as duplicates except the best scored read pair. The score of a read pair is defined as the sum of base qualities >15. Next, the IndelRealigner module in the Genome Analysis Toolkit (GATK) 1.0.5974  was used to perform local realignment around indels to produce an accurate alignment and CountCovariates and TableRecalculation modules to recalibrate the base quality score. An in-house script was applied to modify the read quality, which was generated by BFAST before the GATK recalibration step. The quality scale generated by BFAST presented up to ~63 and was skewed to the maximum value. Such an overestimated quality scale prevented the filtration of false-positive variations while GATK runs genotyping. The in-house script scaled down the overestimated quality values to ~40. SNP and small indel calling were performed using GATK UnifiedGenotyper  with a minimum base quality of Q17 (phred score 17) with “--stand_call_conf 0 --stand_emit_conf 0 --max_deletion_fraction 1.00” and a minimum mapping quality of Q30 (phred score 30) with “--stand_call_conf 0 --stand_emit_conf 0 --genotype_likelihoods_model INDEL --minIndelCnt 3”. Hanwoo, Black Angus, and Holstein were genotyped separately using GATK UnifiedGenotyper. Then, the variants identified in three breeds were merged by genomic position for downstream analysis. A novel variant was defined as one that was not present in the cattle dbSNP 133.
Annotations of variants were based on the 34,577 Cow RefSeq in NCBI (downloaded April 2, 2012). The cattle RefSeqs were aligned against Btau4.0 using BLAT with the ‘fine’ option to obtain the genomic positions of genes, exons, and coding regions. In total, 33,080 RefSeqs were aligned against the reference genome. Among the aligned RefSeqs, the sequences with >90% coverage and a <1% error rate were selected. Then one representative RefSeq was selected from the RefSeqs derived from the same gene. As the result, we selected 29,197 RefSeqs for variant annotation. We identified 2-base canonical splice sites (GU/AG) at the end of an intron as a splice site. The genomic locations of some trait-associated genes that were not obtained from NCBI RefSeqs were defined from previously reported gene information . The selected genes were used to predefine the annotation data of all possible variants and pre-calculate the SIFT  predictions and scores. We selected the coding indels, splice-site variants, and non-synonymous SNPs (NS/SS/Is) that showed SIFT scores of <0.05 as the potentially damaging variants.
Specific NS/SS/I variants were detected by the following criteria: We first selected the NS/SS/Is for which at least 10 reads were aligned and an allele was 50% more abundant than the other alleles for all three breeds at the position.
To measure the genome-wide pattern of selection of a breed, we defined a ROH as follows. The minimum ROH size was set to 400 kb; each chromosome was divided into 400 kb bins, and the ratio of homozygous SNPs per bin was employed as the degree of homozygosity of the bin. To look for a series of high-degree bins rather than separated, one-point peak bins, a degree was smoothed by an average of the two neighbor bins on each side. A continuous extension of bins with a high degree of homozygosity was defined as a ROH. In this study, a 0.8 degree was imposed to determine the ROHs. One breed may contain a ROH that shows a high degree of homozygosity while the others do not. This helps to explain the breed-specific selection pressure. We defined a subset of ROHs that was not duplicated in the other breed’s ROHs as specific ROHs (sROHs). ROHs were identified from SNP chip data using HomozygosityMapper .
To evaluate the accuracy of SNP calling from resequencing of the Hanwoo genome, the same genomic DNA sample was applied to SNP chip analysis. We used BovineSNP50 BeadChip (Illumina)  to genotype the Hanwoo genome. In total, 40 proven bulls in the 45th Hanwoo Performance and Progeny Test Program in Korea, as well as 20 Angus and 19 Holstein individuals, were used for SNP genotyping with the same platform to investigate the ROHs. A consensus SNP genotype was obtained by selecting a maximally expressed genotype from the same location in a breed. Over 90% of the consensus genotypes appeared in more than half of the individuals for all three breeds (data not shown). A ROH was computed from the consensus SNP genotypes with the same criteria that were applied to calculate the ROH of whole-genome SNPs.
Trait-associated genes and QTL regions
We obtained information on trait-associated genes from previous reports to analyze the Kuchinoshima-Ushi breed genome . The genes were categorized into five economic traits: meat, disease resistance, growth, milk, and fecundity. Some genes that did not appear in NCBI RefSeq were added to the gene set for further analysis. QTL regions were identified from information on Cattle QTLs in the Animal QTLdb (Release 17; http://www.animalgenome.org/cgi-bin/QTLdb/BT/index) [20, 21]. QTL locations by bp (UMD 3.1) were downloaded and three types of QTLs were selected: meat fat, meat tenderness, and milk traits. The associated names of these three QTL types described in QTLdb are as follows: intramuscular fat, marbling score, and marbling score (EBV) for meat fat; shear force and tenderness score for meat tenderness; and milk yield, milk yield (daughter deviation), milk yield (EBV), milk yield (PTA), dairy capacity composite index, and dairy form for milk traits.
Because the QTL were based on the UMD 3.1 genome, we converted the locations to coordinates from the Btau4.0 genome. Sequences of the selected QTL were extracted from UMD 3.1 genome sequences and aligned to the Btau4.0 reference genome sequences using LASTZ with the following options: seed = 14 of 22; chain = gapped, step = 5. The alignments were filtered with a minimum of 1,000 bases, 99% average identity, and 5% coverage. The syntenic locations were merged into a large location allowing gaps of 10% at the syntenic locations at most.
Functional enrichment analysis
We determined genes whose genomic positions overlapped partially or completely with the ROH for each breed. We performed functional enrichment analysis against the candidate genes that were within a ROH region within the Gene Ontology and KEGG pathway terms using the Database for Annotation Visualization and Integrated Discovery (DAVID) tool (http://david.abcc.ncifcrf.gov/). Only the enriched GO terms with raw p-values <0.05 were used for further interpretation in this study. The functional relationships of the genes of interest were used in the Pathway studio program (Stratagene) .
Bos taurus reference genome assembly build 4.0
Copy number variant
Database for annotation visualization and integrated discovery
Estimated breeding value
Extended haplotype homozygosity
Genome analysis tool kit
Integrated haplotype score
Kyoto encyclopedia of genes and genomes
Predicted transmitting ability
Quantitative trait loci
Reference sequence database
Region of homozygosity
This work was supported by a grant from the BioGreen21 program, Rural Development Administration (RDA), Korea (grant no. PJ007197), a grant from the Korea Research Institute of Bioscience and Biotechnology Research Initiative Program, and through funds from the University of Alberta program provided to Professor Moore. Xiaoping Liao is supported by the Genome Canada project titled “Whole Genome Selection through Genome Wide Imputation in Beef Cattle.”
- Tellam RL, Lemay DG, Van Tassell CP, Lewin HA, Worley KC, Elsik CG: Unlocking the bovine genome. BMC Genomics. 2009, 10: 193-10.1186/1471-2164-10-193.PubMed CentralView ArticlePubMedGoogle Scholar
- Elsik CG, Tellam RL, Worley KC, Gibbs RA, Muzny DM, Weinstock GM, Adelson DL, Eichler EE, Elnitski L, Guigo R: The genome sequence of taurine cattle: a window to ruminant biology and evolution. Science. 2009, 324 (5926): 522-528.PubMed CentralView ArticlePubMedGoogle Scholar
- Liu Y, Qin X, Song XZ, Jiang H, Shen Y, Durbin KJ, Lien S, Kent MP, Sodeland M, Ren Y: Bos taurus genome assembly. BMC Genomics. 2009, 10: 180-10.1186/1471-2164-10-180.PubMed CentralView ArticlePubMedGoogle Scholar
- Zimin AV, Delcher AL, Florea L, Kelley DR, Schatz MC, Puiu D, Hanrahan F, Pertea G, Van Tassell CP, Sonstegard TS: A whole-genome assembly of the domestic cow, Bos taurus. Genome Biol. 2009, 10 (4): R42-10.1186/gb-2009-10-4-r42.PubMed CentralView ArticlePubMedGoogle Scholar
- Van Tassell CP, Smith TP, Matukumalli LK, Taylor JF, Schnabel RD, Lawley CT, Haudenschild CD, Moore SS, Warren WC, Sonstegard TS: SNP discovery and allele frequency estimation by deep sequencing of reduced representation libraries. Nat Methods. 2008, 5 (3): 247-252. 10.1038/nmeth.1185.View ArticlePubMedGoogle Scholar
- Eck SH, Benet-Pages A, Flisikowski K, Meitinger T, Fries R, Strom TM: Whole genome sequencing of a single Bos taurus animal for single nucleotide polymorphism discovery. Genome Biol. 2009, 10 (8): R82-10.1186/gb-2009-10-8-r82.PubMed CentralView ArticlePubMedGoogle Scholar
- Kawahara-Miki R, Tsuda K, Shiwa Y, Arai-Kichise Y, Matsumoto T, Kanesaki Y, Oda S, Ebihara S, Yajima S, Yoshikawa H: Whole-genome resequencing shows numerous genes with nonsynonymous SNPs in the Japanese native cattle Kuchinoshima-Ushi. BMC Genomics. 2011, 12: 103-10.1186/1471-2164-12-103.PubMed CentralView ArticlePubMedGoogle Scholar
- Stothard P, Choi JW, Basu U, Sumner-Thomson JM, Meng Y, Liao X, Moore SS: Whole genome resequencing of black Angus and Holstein cattle for SNP and CNV discovery. BMC Genomics. 2011, 12: 559-10.1186/1471-2164-12-559.PubMed CentralView ArticlePubMedGoogle Scholar
- Gibson J, Morton NE, Collins A: Extended tracts of homozygosity in outbred human populations. Hum Mol Genet. 2006, 15 (5): 789-795. 10.1093/hmg/ddi493.View ArticlePubMedGoogle Scholar
- Voight BF, Kudaravalli S, Wen X, Pritchard JK: A map of recent positive selection in the human genome. PLoS Biol. 2006, 4 (3): e72-10.1371/journal.pbio.0040072.PubMed CentralView ArticlePubMedGoogle Scholar
- Weir BS, Cardon LR, Anderson AD, Nielsen DM, Hill WG: Measures of human population structure show heterogeneity among genomic regions. Genome Res. 2005, 15 (11): 1468-1476. 10.1101/gr.4398405.PubMed CentralView ArticlePubMedGoogle Scholar
- Sabeti PC, Reich DE, Higgins JM, Levine HZ, Richter DJ, Schaffner SF, Gabriel SB, Platko JV, Patterson NJ, McDonald GJ: Detecting recent positive selection in the human genome from haplotype structure. Nature. 2002, 419 (6909): 832-837. 10.1038/nature01140.View ArticlePubMedGoogle Scholar
- Lee C, Pollak EJ: Genetic antagonism between body weight and milk production in beef cattle. J Anim Sci. 2002, 80 (2): 316-321.PubMedGoogle Scholar
- Han SW: The breed of cattle. Breeds of Livestock. 1996, Seoul: Sun-Jin publishing, 148-160. 1Google Scholar
- McKenna A, Hanna M, Banks E, Sivachenko A, Cibulskis K, Kernytsky A, Garimella K, Altshuler D, Gabriel S, Daly M: The genome analysis toolkit: a MapReduce framework for analyzing next-generation DNA sequencing data. Genome Res. 2010, 20 (9): 1297-1303. 10.1101/gr.107524.110.PubMed CentralView ArticlePubMedGoogle Scholar
- DePristo MA, Banks E, Poplin R, Garimella KV, Maguire JR, Hartl C, Philippakis AA, del Angel G, Rivas MA, Hanna M: A framework for variation discovery and genotyping using next-generation DNA sequencing data. Nat Genet. 2011, 43 (5): 491-498. 10.1038/ng.806.PubMed CentralView ArticlePubMedGoogle Scholar
- Decker JE, Pires JC, Conant GC, McKay SD, Heaton MP, Chen K, Cooper A, Vilkki J, Seabury CM, Caetano AR: Resolving the evolution of extant and extinct ruminants with high-throughput phylogenomics. Proc Natl Acad Sci USA. 2009, 106 (44): 18644-18649. 10.1073/pnas.0904691106.PubMed CentralView ArticlePubMedGoogle Scholar
- Marquez B, Ameye G, Vallet CM, Tulkens PM, Poirel HA, Van Bambeke F: Characterization of Abcc4 gene amplification in stepwise-selected mouse J774 macrophages resistant to the topoisomerase II inhibitor ciprofloxacin. PLoS One. 2011, 6 (12): e28368-10.1371/journal.pone.0028368.PubMed CentralView ArticlePubMedGoogle Scholar
- Schaeren-Wiemers N, Bonnet A, Erb M, Erne B, Bartsch U, Kern F, Mantei N, Sherman D, Suter U: The raft-associated protein MAL is required for maintenance of proper axon–glia interactions in the central nervous system. J Cell Biol. 2004, 166 (5): 731-742. 10.1083/jcb.200406092.PubMed CentralView ArticlePubMedGoogle Scholar
- Hu ZL, Fritz ER, Reecy JM: AnimalQTLdb: a livestock QTL database tool set for positional QTL information mining and beyond. Nucleic Acids Res. 2007, 35 (Database issue): D604-D609.PubMed CentralView ArticlePubMedGoogle Scholar
- Hu ZL, Reecy JM: Animal QTLdb: beyond a repository. A public platform for QTL comparisons and integration with diverse types of structural genomic information. Mamm Genome. 2007, 18 (1): 1-4. 10.1007/s00335-006-0105-8.View ArticlePubMedGoogle Scholar
- Sasaki Y, Nagai K, Nagata Y, Doronbekov K, Nishimura S, Yoshioka S, Fujita T, Shiga K, Miyake T, Taniguchi Y: Exploration of genes showing intramuscular fat deposition-associated expression changes in musculus longissimus muscle. Anim Genet. 2006, 37 (1): 40-46. 10.1111/j.1365-2052.2005.01380.x.View ArticlePubMedGoogle Scholar
- Yamada T, Sasaki S, Sukegawa S, Yoshioka S, Takahagi Y, Morita M, Murakami H, Morimatsu F, Fujita T, Miyake T: Association of a single nucleotide polymorphism in titin gene with marbling in Japanese Black beef cattle. BMC Res Notes. 2009, 2: 78-10.1186/1756-0500-2-78.PubMed CentralView ArticlePubMedGoogle Scholar
- Labeit S, Kolmerer B: Titins: giant proteins in charge of muscle ultrastructure and elasticity. Science. 1995, 270 (5234): 293-296. 10.1126/science.270.5234.293.View ArticlePubMedGoogle Scholar
- Ku CS, Naidoo N, Teo SM, Pawitan Y: Regions of homozygosity and their impact on complex diseases and traits. Hum Genet. 2011, 129 (1): 1-15. 10.1007/s00439-010-0920-6.View ArticlePubMedGoogle Scholar
- Weidemann JF, Kaess G, Carrljthers LD: The histology of pre-rigor and post-rigor ox muscle before and after cooking and its relation to tenderness. J Food Sci. 1967, 32 (1): 7-13. 10.1111/j.1365-2621.1967.tb01946.x.View ArticleGoogle Scholar
- Solomon C, White JH, Kremer R: Mitogen-activated protein kinase inhibits 1,25-dihydroxyvitamin D3-dependent signal transduction by phosphorylating human retinoid X receptor alpha. J Clin Invest. 1999, 103 (12): 1729-1735. 10.1172/JCI6871.PubMed CentralView ArticlePubMedGoogle Scholar
- Lim D, Kim NK, Park HS, Lee SH, Cho YM, Oh SJ, Kim TH, Kim H: Identification of candidate genes related to bovine marbling using protein-protein interaction networks. Int J Biol Sci. 2011, 7 (7): 992-1002.PubMed CentralView ArticlePubMedGoogle Scholar
- Huq MD, Tsai NP, Lin YP, Higgins L, Wei LN: Vitamin B6 conjugation to nuclear corepressor RIP140 and its role in gene regulation. Nat Chem Biol. 2007, 3 (3): 161-165. 10.1038/nchembio861.View ArticlePubMedGoogle Scholar
- Brandes R, Arad R, Bar-Tana J: Inducers of adipose conversion activate transcription promoted by a peroxisome proliferators response element in 3T3-L1 cells. Biochem Pharmacol. 1995, 50 (11): 1949-1951. 10.1016/0006-2952(95)02082-9.View ArticlePubMedGoogle Scholar
- Kuhn C, Weikard R: An investigation into the genetic background of coat colour dilution in a Charolais x German Holstein F2 resource population. Anim Genet. 2007, 38 (2): 109-113. 10.1111/j.1365-2052.2007.01569.x.View ArticlePubMedGoogle Scholar
- Abu Safieh L, Aldahmesh MA, Shamseldin H, Hashem M, Shaheen R, Alkuraya H, Al Hazzaa SA, Al-Rajhi A, Alkuraya FS: Clinical and molecular characterisation of Bardet-Biedl syndrome in consanguineous populations: the power of homozygosity mapping. J Med Genet. 2010, 47 (4): 236-241. 10.1136/jmg.2009.070755.View ArticlePubMedGoogle Scholar
- Collin RW, Safieh C, Littink KW, Shalev SA, Garzozi HJ, Rizel L, Abbasi AH, Cremers FP, den Hollander AI, Klevering BJ: Mutations in C2ORF71 cause autosomal-recessive retinitis pigmentosa. Am J Hum Genet. 2010, 86 (5): 783-788. 10.1016/j.ajhg.2010.03.016.PubMed CentralView ArticlePubMedGoogle Scholar
- Harville HM, Held S, Diaz-Font A, Davis EE, Diplas BH, Lewis RA, Borochowitz ZU, Zhou W, Chaki M, MacDonald J: Identification of 11 novel mutations in eight BBS genes by high-resolution homozygosity mapping. J Med Genet. 2010, 47 (4): 262-267. 10.1136/jmg.2009.071365.PubMed CentralView ArticlePubMedGoogle Scholar
- Iseri SU, Wyatt AW, Nurnberg G, Kluck C, Nurnberg P, Holder GE, Blair E, Salt A, Ragge NK: Use of genome-wide SNP homozygosity mapping in small pedigrees to identify new mutations in VSX2 causing recessive microphthalmia and a semidominant inner retinal dystrophy. Hum Genet. 2010, 128 (1): 51-60. 10.1007/s00439-010-0823-6.View ArticlePubMedGoogle Scholar
- Lapunzina P, Aglan M, Temtamy S, Caparros-Martin JA, Valencia M, Leton R, Martinez-Glez V, Elhossini R, Amr K, Vilaboa N: Identification of a frameshift mutation in Osterix in a patient with recessive osteogenesis imperfecta. Am J Hum Genet. 2010, 87 (1): 110-114. 10.1016/j.ajhg.2010.05.016.PubMed CentralView ArticlePubMedGoogle Scholar
- Nicolas E, Poitelon Y, Chouery E, Salem N, Levy N, Megarbane A, Delague V: CAMOS, a nonprogressive, autosomal recessive, congenital cerebellar ataxia, is caused by a mutant zinc-finger protein, ZNF592. Eur J Hum Genet. 2010, 18 (10): 1107-1113. 10.1038/ejhg.2010.82.PubMed CentralView ArticlePubMedGoogle Scholar
- Pang J, Zhang S, Yang P, Hawkins-Lee B, Zhong J, Zhang Y, Ochoa B, Agundez JA, Voelckel MA, Fisher RB: Loss-of-function mutations in HPSE2 cause the autosomal recessive urofacial syndrome. Am J Hum Genet. 2010, 86 (6): 957-962. 10.1016/j.ajhg.2010.04.016.PubMed CentralView ArticlePubMedGoogle Scholar
- Uz E, Alanay Y, Aktas D, Vargel I, Gucer S, Tuncbilek G, von Eggeling F, Yilmaz E, Deren O, Posorski N: Disruption of ALX1 causes extreme microphthalmia and severe facial clefting: expanding the spectrum of autosomal-recessive ALX-related frontonasal dysplasia. Am J Hum Genet. 2010, 86 (5): 789-796. 10.1016/j.ajhg.2010.04.002.PubMed CentralView ArticlePubMedGoogle Scholar
- Walsh T, Shahin H, Elkan-Miller T, Lee MK, Thornton AM, Roeb W, Abu Rayyan A, Loulus S, Avraham KB, King MC: Whole exome sequencing and homozygosity mapping identify mutation in the cell polarity protein GPSM2 as the cause of nonsyndromic hearing loss DFNB82. Am J Hum Genet. 2010, 87 (1): 90-94. 10.1016/j.ajhg.2010.05.010.PubMed CentralView ArticlePubMedGoogle Scholar
- Lencz T, Lambert C, DeRosse P, Burdick KE, Morgan TV, Kane JM, Kucherlapati R, Malhotra AK: Runs of homozygosity reveal highly penetrant recessive loci in schizophrenia. Proc Natl Acad Sci USA. 2007, 104 (50): 19942-19947. 10.1073/pnas.0710021104.PubMed CentralView ArticlePubMedGoogle Scholar
- Nalls MA, Guerreiro RJ, Simon-Sanchez J, Bras JT, Traynor BJ, Gibbs JR, Launer L, Hardy J, Singleton AB: Extended tracts of homozygosity identify novel candidate genes associated with late-onset Alzheimer’s disease. Neurogenetics. 2009, 10 (3): 183-190. 10.1007/s10048-009-0182-4.PubMed CentralView ArticlePubMedGoogle Scholar
- Yang TL, Guo Y, Zhang LS, Tian Q, Yan H, Papasian CJ, Recker RR, Deng HW: Runs of homozygosity identify a recessive locus 12q21.31 for human adult height. J Clin Endocrinol Metab. 2010, 95 (8): 3777-3782. 10.1210/jc.2009-1715.PubMed CentralView ArticlePubMedGoogle Scholar
- Barendse W, Bunch RJ, Thomas MB, Harrison BE: A splice site single nucleotide polymorphism of the fatty acid binding protein 4 gene appears to be associated with intramuscular fat deposition in longissimus muscle in Australian cattle. Anim Genet. 2009, 40 (5): 770-773. 10.1111/j.1365-2052.2009.01913.x.View ArticlePubMedGoogle Scholar
- Hoashi S, Hinenoya T, Tanaka A, Ohsaki H, Sasazaki S, Taniguchi M, Oyama K, Mukai F, Mannen H: Association between fatty acid compositions and genotypes of FABP4 and LXR-alpha in Japanese black cattle. BMC Genet. 2008, 9: 84-PubMed CentralView ArticlePubMedGoogle Scholar
- Cho S, Park TS, Yoon DH, Cheong HS, Namgoong S, Park BL, Lee HW, Han CS, Kim EM, Cheong IC: Identification of genetic polymorphisms in FABP3 and FABP4 and putative association with back fat thickness in Korean native cattle. BMB Rep. 2008, 41 (1): 29-34. 10.5483/BMBRep.2008.41.1.029.View ArticlePubMedGoogle Scholar
- Lee SH, van der Werf JH, Park EW, Oh SJ, Gibson JP, Thompson JM: Genetic polymorphisms of the bovine fatty acid binding protein 4 gene are significantly associated with marbling and carcass weight in Hanwoo (Korean Cattle). Anim Genet. 2010, 41 (4): 442-444.PubMedGoogle Scholar
- Laliotis GP, Bizelis I, Rogdakis E: Comparative approach of the de novo fatty acid synthesis (Lipogenesis) between ruminant and non ruminant mammalian species: from biochemical level to the main regulatory lipogenic genes. Curr Genomics. 2010, 11 (3): 168-183. 10.2174/138920210791110960.PubMed CentralView ArticlePubMedGoogle Scholar
- Lee SH, Kim SC, Choi BH, Lim D, Kim NK, Lee JH, Kim OH, Lee CS, Kim HC, Yang BS: mt-COX1, mt-ND1 and CREBP are indicators of intramuscular fat content in Hanwoo (Korean cattle). Livest Sci. 2012, 146: 160-167. 10.1016/j.livsci.2012.03.003.View ArticleGoogle Scholar
- Yamamoto KK, Gonzalez GA, Biggs WH, Montminy MR: Phosphorylation-induced binding and transcriptional efficacy of nuclear factor CREB. Nature. 1988, 334 (6182): 494-498. 10.1038/334494a0.View ArticlePubMedGoogle Scholar
- Casimir DA, Ntambi JM: cAMP activates the expression of stearoyl-CoA desaturase gene 1 during early preadipocyte differentiation. J Biol Chem. 1996, 271 (47): 29847-29853. 10.1074/jbc.271.47.29847.View ArticlePubMedGoogle Scholar
- Wang YH, Bower NI, Reverter A, Tan SH, De Jager N, Wang R, McWilliam SM, Cafe LM, Greenwood PL, Lehnert SA: Gene expression patterns during intramuscular fat development in cattle. J Anim Sci. 2009, 87 (1): 119-130.View ArticlePubMedGoogle Scholar
- Dooley KA, Bennett MK, Osborne TF: A critical role for CREB as a co-activator in sterol regulated transcription of HMG CoA synthase promoter. J Biol Chem. 1999, 274: 5285-5291. 10.1074/jbc.274.9.5285.View ArticlePubMedGoogle Scholar
- Homer N, Merriman B, Nelson SF: BFAST: an alignment tool for large scale genome resequencing. PLoS One. 2009, 4 (11): e7767-10.1371/journal.pone.0007767.PubMed CentralView ArticlePubMedGoogle Scholar
- Kumar P, Henikoff S, Ng PC: Predicting the effects of coding non-synonymous variants on protein function using the SIFT algorithm. Nat Protoc. 2009, 4 (7): 1073-1081.View ArticlePubMedGoogle Scholar
- Seelow D, Schuelke M, Hildebrandt F, Nurnberg P: HomozygosityMapper--an interactive approach to homozygosity mapping. Nucleic Acids Res. 2009, 37 (Web Server issue): W593-W599.PubMed CentralView ArticlePubMedGoogle Scholar
- Matukumalli LK, Lawley CT, Schnabel RD, Taylor JF, Allan MF, Heaton MP, O’Connell J, Moore SS, Smith TP, Sonstegard TS: Development and characterization of a high density SNP genotyping assay for cattle. PLoS One. 2009, 4 (4): e5350-10.1371/journal.pone.0005350.PubMed CentralView ArticlePubMedGoogle Scholar
- Nikitin A, Egorov S, Daraselia N, Mazo I: Pathway studio–the analysis and navigation of molecular networks. Bioinformatics. 2003, 19 (16): 2155-10.1093/bioinformatics/btg290.View ArticlePubMedGoogle Scholar
This article is published under license to BioMed Central Ltd. 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.