Using online tools at the Bovine Genome Database to manually annotate genes in the new reference genome.

With the availability of a new highly contiguous Bos taurus reference genome assembly (ARS-UCD1.2), it is the opportune time to upgrade the bovine gene set by seeking input from researchers. Furthermore, advances in graphical genome annotation tools now make it possible for researchers to leverage sequence data generated with the latest technologies to collaboratively curate genes. For many years the Bovine Genome Database (BGD) has provided tools such as the Apollo genome annotation editor to support manual bovine gene curation. The goal of this paper is to explain the reasoning behind the decisions made in the manual gene curation process while providing examples using the existing BGD tools. We will describe the sources of gene annotation evidence provided at the BGD, including RNA-seq and Iso-Seq data. We will also explain how to interpret various data visualizations when curating gene models, and will demonstrate the value of manual gene annotation. The process described here can be applied to manual gene curation for other species with similar tools. With a better understanding of manual gene annotation, researchers will be encouraged to edit gene models and contribute to the enhancement of livestock gene sets.

Symposium review: Advances in sequencing technology herald a new frontier in cattle genomics and genome-enabled selection.

The cattle reference genome assembly has underpinned major innovations in beef and dairy genetics through genome-enabled selection, including removal of deleterious recessive variants and selection for favorable alleles affecting quantitative production traits. The initial reference assemblies, up to and including UMD3.1 and Btau4.1, were based on a combination of clone-by-clone sequencing of bacterial artificial chromosome clones generated from blood DNA of a Hereford bull and whole-genome shotgun sequencing of blood DNA from his inbred daughter/granddaughter named L1 Dominette 01449 (Dominette). The approach introduced assembly gaps, misassemblies, and errors, and it limited the ability to assemble regions that undergo rearrangement in blood cells, such as immune gene clusters. Nonetheless, the reference supported the creation of genotyping tools and provided a basis for many studies of gene expression. Recently, long-read sequencing technologies have emerged that facilitated a re-assembly of the reference genome, using lung tissue from Dominette to resolve many of the problems and providing a bridge to place historical studies in common context. The new reference, ARS-UCD1.2, successfully assembled germline immune gene clusters and improved overall continuity (i.e., reduction of gaps and inversions) by over 250-fold. This reference properly places nearly all of the legacy genetic markers used for over a decade in the industry. In this review, we discuss the improvements made to the cattle reference; remaining issues present in the assembly; tools developed to support genome-based studies in beef and dairy cattle; and the emergence of newer genome assembly methods that are producing even higher-quality assemblies for other breeds of cattle at a fraction of the cost. The new frontier for cattle genomics research will likely include a transition from the individual Hereford reference genome, to a "pan-genome" reference, representing all the DNA segments existing in commonly used cattle breeds, bringing the cattle reference into line with the current direction of human genome research.

A homozygous ADAMTS2 nonsense mutation in a Doberman Pinscher dog with Ehlers Danlos syndrome and extreme skin fragility.

An eight-week old Doberman Pinscher was diagnosed with Ehlers Danlos syndrome based on the dog's hyper-mobile carpal, tarsal and stifle joints and abnormal skin. The skin was loose and hyper-elastic with several wounds and large atrophic scars. The dog was euthanized after a severe degloving injury from minimal trauma. A whole-genome sequence, generated with DNA from the dog's blood, contained a rare, homozygous C-to-T transition at position 2408978 on chromosome 11. This transition is predicted to alter the ADAMTS2 transcript (ADAMTS2:c.769C>T) and encode a nonsense mutation (p.Arg257Ter). Biallelic ADAMTS2 mutations have caused a type of Ehlers Danlos syndrome known as dermatosparaxis in other species.

GM2 Gangliosidosis in Shiba Inu Dogs with an In-Frame Deletion in HEXB.

Consistent with a tentative diagnosis of neuronal ceroid lipofuscinosis (NCL), autofluorescent cytoplasmic storage bodies were found in neurons from the brains of 2 related Shiba Inu dogs with a young-adult onset, progressive neurodegenerative disease. Unexpectedly, no potentially causal NCL-related variants were identified in a whole-genome sequence generated with DNA from 1 of the affected dogs. Instead, the whole-genome sequence contained a homozygous 3 base pair (bp) deletion in a coding region of HEXB. The other affected dog also was homozygous for this 3-bp deletion. Mutations in the human HEXB ortholog cause Sandhoff disease, a type of GM2 gangliosidosis. Thin-layer chromatography confirmed that GM2 ganglioside had accumulated in an affected Shiba Inu brain. Enzymatic analysis confirmed that the GM2 gangliosidosis resulted from a deficiency in the HEXB encoded protein and not from a deficiency in products from HEXA or GM2A, which are known alternative causes of GM2 gangliosidosis. We conclude that the homozygous 3-bp deletion in HEXB is the likely cause of the Shiba Inu neurodegenerative disease and that whole-genome sequencing can lead to the early identification of potentially disease-causing DNA variants thereby refocusing subsequent diagnostic analyses toward confirming or refuting candidate variant causality.

Homozygous PPT1 Splice Donor Mutation in a Cane Corso Dog With Neuronal Ceroid Lipofuscinosis.

A 10-month-old spayed female Cane Corso dog was evaluated after a 2-month history of progressive blindness, ataxia, and lethargy. Neurologic examination abnormalities indicated a multifocal lesion with primarily cerebral and cerebellar signs. Clinical worsening resulted in humane euthanasia. On necropsy, there was marked astrogliosis throughout white matter tracts of the cerebrum, most prominently in the corpus callosum. In the cerebral cortex and midbrain, most neurons contained large amounts of autofluorescent storage material in the perinuclear area of the cells. Cerebellar storage material was present in the Purkinje cells, granular cell layer, and perinuclear regions of neurons in the deep nuclei. Neuronal ceroid lipofuscinosis (NCL) was diagnosed. Whole genome sequencing identified a PPT1c.124 + 1G>A splice donor mutation. This nonreference assembly allele was homozygous in the affected dog, has not previously been reported in dbSNP, and was absent from the whole genome sequences of 45 control dogs and 31 unaffected Cane Corsos. Our findings indicate a novel mutation causing the CLN1 form of NCL in a previously unreported dog breed. A canine model for CLN1 disease could provide an opportunity for therapeutic advancement, benefiting both humans and dogs with this disorder.

Australian Cattle Dogs with Neuronal Ceroid Lipofuscinosis are Homozygous for a CLN5 Nonsense Mutation Previously Identified in Border Collies.

BACKGROUND: Neuronal ceroid lipofuscinosis (NCL), a fatal neurodegenerative disease, has been diagnosed in young adult Australian Cattle Dogs. OBJECTIVE: Characterize the Australian Cattle Dog form of NCL and determine its molecular genetic cause. ANIMALS: Tissues from 4 Australian Cattle Dogs with NCL-like signs and buccal swabs from both parents of a fifth affected breed member. Archived DNA samples from 712 individual dogs were genotyped. METHODS: Tissues were examined by fluorescence, electron, and immunohistochemical microscopy. A whole-genome sequence was generated for 1 affected dog. A TaqMan allelic discrimination assay was used for genotyping. RESULTS: The accumulation of autofluorescent cytoplasmic storage material with characteristic ultrastructure in tissues from the 4 affected dogs supported a diagnosis of NCL. The whole-genome sequence contained a homozygous nonsense mutation: CLN5:c.619C>T. All 4 DNA samples from clinically affected dogs tested homozygous for the variant allele. Both parents of the fifth affected dog were heterozygotes. Archived DNA samples from 346 Australian Cattle Dogs, 188 Border Collies, and 177 dogs of other breeds were homozygous for the reference allele. One archived Australian Cattle Dog sample was from a heterozygote. CONCLUSIONS AND CLINICAL IMPORTANCE: The homozygous CLN5 nonsense is almost certainly causal because the same mutation previously had been reported to cause a similar form of NCL in Border Collies. Identification of the molecular genetic cause of Australian Cattle Dog NCL will allow the use of DNA tests to confirm the diagnosis of NCL in this breed.

Origins of cattle on Chirikof Island, Alaska, elucidated from genome-wide SNP genotypes.

Feral livestock may harbor genetic variation of commercial, scientific, historical or esthetic value. The origins and uniqueness of feral cattle on Chirikof Island, Alaska, are uncertain. The island is now part of the Alaska Maritime Wildlife Refuge and Federal wildlife managers want grazing to cease, presumably leading to demise of the cattle. Here we characterize the cattle of Chirikof Island relative to extant breeds and discern their origins. Our analyses support the inference that Yakut cattle from Russia arrived first on Chirikof Island, then ~120 years ago the first European taurine cattle were introduced to the island, and finally a large wave of Hereford cattle were introduced on average 40 years ago. In addition, this mixture of European and East-Asian cattle is unique compared with other North American breeds and we find evidence that natural selection in the relatively harsh environment of Chirikof Island has further impacted their genetic architecture. These results provide an objective basis for decisions regarding conservation of the Chirikof Island cattle.

Golden Retriever dogs with neuronal ceroid lipofuscinosis have a two-base-pair deletion and frameshift in CLN5.

We studied a recessive, progressive neurodegenerative disease occurring in Golden Retriever siblings with an onset of signs at 15 months of age. As the disease progressed these signs included ataxia, anxiety, pacing and circling, tremors, aggression, visual impairment and localized and generalized seizures. A whole genome sequence, generated with DNA from one affected dog, contained a plausibly causal homozygous mutation: CLN5:c.934_935delAG. This mutation was predicted to produce a frameshift and premature termination codon and encode a protein variant, CLN5:p.E312Vfs*6, which would lack 39 C-terminal amino acids. Eighteen DNA samples from the Golden Retriever family members were genotyped at CLN5:c.934_935delAG. Three clinically affected dogs were homozygous for the deletion allele; whereas, the clinically normal family members were either heterozygotes (n = 11) or homozygous for the reference allele (n = 4). Among archived Golden Retrievers DNA samples with incomplete clinical records that were also genotyped at the CLN5:c.934_935delAG variant, 1053 of 1062 were homozygous for the reference allele, 8 were heterozygotes and one was a deletion-allele homozygote. When contacted, the owner of this homozygote indicated that their dog had been euthanized because of a neurologic disease that progressed similarly to that of the affected Golden Retriever siblings. We have collected and stored semen from a heterozygous Golden Retriever, thereby preserving an opportunity for us or others to establish a colony of CLN5-deficient dogs.

A homozygous KCNJ10 mutation in Jack Russell Terriers and related breeds with spinocerebellar ataxia with myokymia, seizures, or both.

BACKGROUND: Juvenile-onset spinocerebellar ataxia has been recognized in Jack Russell Terriers and related Russell group terriers (RGTs) for over 40 years. Ataxia occurs with varying combinations of myokymia, seizures, and other signs of neurologic disease. More than 1 form of the disease has been suspected. HYPOTHESIS/OBJECTIVES: The objective was to identify the mutation causing the spinocerebellar ataxia associated with myokymia, seizures, or both and distinguish the phenotype from other ataxias in the RGTs. ANIMALS: DNA samples from 16 RGTs with spinocerebellar ataxia beginning from 2 to 12 months of age, 640 control RGTs, and 383 dogs from 144 other breeds along with the medical records of affected dogs were studied. METHODS: This case-control study compared the frequencies of a KCNJ10 allele in RGTs with spinocerebellar ataxia versus control RGTs. This allele was identified in a whole-genome sequence of a single RGT with spinocerebellar ataxia and myokymia by comparison to whole-genome sequences from 81 other canids that were normal or had other diseases. RESULTS: A missense mutation in the gene coding for the inwardly rectifying potassium channel Kir4.1 (KCNJ10:c.627C>G) was significantly (P < .001) associated with the disease. Dogs homozygous for the mutant allele all had spinocerebellar ataxia with varying combinations of myokymia and seizures. CONCLUSIONS AND CLINICAL IMPORTANCE: Identification of the KCNJ10 mutation in dogs with spinocerebellar ataxia with myokymia, seizures, or both clarifies the multiple forms of ataxia seen in these breeds and provides a DNA test to identify carriers.

Genome scan for meat quality traits in Nelore beef cattle.

Meat quality traits are economically important because they affect consumers' acceptance, which, in turn, influences the demand for beef. However, selection to improve meat quality is limited by the small numbers of animals on which meat tenderness can be evaluated due to the cost of performing shear force analysis and the resultant damage to the carcass. Genome wide-association studies for Warner-Bratzler shear force measured at different times of meat aging, backfat thickness, ribeye muscle area, scanning parameters [lightness, redness (a*), and yellowness] to ascertain color characteristics of meat and fat, water-holding capacity, cooking loss (CL), and muscle pH were conducted using genotype data from the Illumina BovineHD BeadChip array to identify quantitative trait loci (QTL) in all phenotyped Nelore cattle. Phenotype count for these animals ranged from 430 to 536 across traits. Meat quality traits in Nelore are controlled by numerous QTL of small effect, except for a small number of large-effect QTL identified for a*fat, CL, and pH. Genomic regions harboring these QTL and the pathways in which the genes from these regions act appear to differ from those identified in taurine cattle for meat quality traits. These results will guide future QTL mapping studies and the development of models for the prediction of genetic merit to implement genomic selection for meat quality in Nelore cattle.

The accuracies of DNA-based estimates of genetic merit derived from Angus or multibreed beef cattle training populations.

Several organizations have developed prediction models for molecular breeding values (MBV) for quantitative growth and carcass traits in beef cattle using Bovine SNP50 genotypes and phenotypic or EBV data. Molecular breeding values for Angus cattle have been developed by IGENITY, Pfizer Animal Genetics, and a collaboration between researchers from Iowa State University and the University of Missouri-Columbia (ISU/UMC). The U.S. Meat Animal Research Center (USMARC; Clay Center, NE) has also developed MBV for 16 cattle breeds using 2 multibreed populations, the Germplasm Evaluation (GPE) Program and the 2,000 Bull Project (2K(ALL)), and 2 single breed subpopulations of the 2,000 Bull Project, Angus (2K(AN)) and Hereford (2K(HH)). In this study, these MBV were assessed relative to commercial ranch EBV estimated from the progeny phenotypes of Angus bulls naturally mated in multisire breeding pastures to commercial cows: 121 for USMARC MBV, 99 for ISU/UMC MBV, and 29 for IGENITY and Pfizer MBV (selected based on number of progeny carcass records). Five traits were analyzed: weaning weight (WW), HCW, marbling score (MS), rib-eye muscle area (RE), and, for IGENITY and Pfizer only, feedlot ADG. The average accuracies of MBV across traits were 0.38 ± 0.05 for IGENITY, 0.61 ± 0.12 for Pfizer, 0.46 ± 0.12 for ISU/UMC, 0.16 ± 0.04 for GPE, 0.26 ± 0.05 for 2K(ALL), 0.24 ± 0.04 for 2K(AN), and 0.02 ± 0.12 for 2K(HH). Angus-based MBV (IGENITY, Pfizer, ISU/UMC, and 2K(AN)) explained larger proportions of genetic variance in this population than GPE, 2K(ALL), or 2K(HH) MBV for the same traits. In this data set, IGENITY, Pfizer, and ISU/UMC MBV were predictive of realized performance of progeny, and incorporation of that information into national genetic evaluations would be expected to improve EPD accuracy, particularly for young animals.

Genome-wide association analysis for quantitative trait loci influencing Warner-Bratzler shear force in five taurine cattle breeds.

We performed a genome-wide association study for Warner-Bratzler shear force (WBSF), a measure of meat tenderness, by genotyping 3360 animals from five breeds with 54 790 BovineSNP50 and 96 putative single-nucleotide polymorphisms (SNPs) within µ-calpain [HUGO nomenclature calpain 1, (mu/I) large subunit; CAPN1] and calpastatin (CAST). Within- and across-breed analyses estimated SNP allele substitution effects (ASEs) by genomic best linear unbiased prediction (GBLUP) and variance components by restricted maximum likelihood under an animal model incorporating a genomic relationship matrix. GBLUP estimates of ASEs from the across-breed analysis were moderately correlated (0.31-0.66) with those from the individual within-breed analyses, indicating that prediction equations for molecular estimates of breeding value developed from across-breed analyses should be effective for genomic selection within breeds. We identified 79 genomic regions associated with WBSF in at least three breeds, but only eight were detected in all five breeds, suggesting that the within-breed analyses were underpowered, that different quantitative trait loci (QTL) underlie variation between breeds or that the BovineSNP50 SNP density is insufficient to detect common QTL among breeds. In the across-breed analysis, CAPN1 was followed by CAST as the most strongly associated WBSF QTL genome-wide, and associations with both were detected in all five breeds. We show that none of the four commercialized CAST and CAPN1 SNP diagnostics are causal for associations with WBSF, and we putatively fine-map the CAPN1 causal mutation to a 4581-bp region. We estimate that variation in CAST and CAPN1 explains 1.02 and 1.85% of the phenotypic variation in WBSF respectively.

Genome-wide association analysis for feed efficiency in Angus cattle.

Estimated breeding values for average daily feed intake (AFI; kg/day), residual feed intake (RFI; kg/day) and average daily gain (ADG; kg/day) were generated using a mixed linear model incorporating genomic relationships for 698 Angus steers genotyped with the Illumina BovineSNP50 assay. Association analyses of estimated breeding values (EBVs) were performed for 41,028 single nucleotide polymorphisms (SNPs), and permutation analysis was used to empirically establish the genome-wide significance threshold (P < 0.05) for each trait. SNPs significantly associated with each trait were used in a forward selection algorithm to identify genomic regions putatively harbouring genes with effects on each trait. A total of 53, 66 and 68 SNPs explained 54.12% (24.10%), 62.69% (29.85%) and 55.13% (26.54%) of the additive genetic variation (when accounting for the genomic relationships) in steer breeding values for AFI, RFI and ADG, respectively, within this population. Evaluation by pathway analysis revealed that many of these SNPs are in genomic regions that harbour genes with metabolic functions. The presence of genetic correlations between traits resulted in 13.2% of SNPs selected for AFI and 4.5% of SNPs selected for RFI also being selected for ADG in the analysis of breeding values. While our study identifies panels of SNPs significant for efficiency traits in our population, validation of all SNPs in independent populations will be necessary before commercialization.

Mapping of fertility traits in Finnish Ayrshire by genome-wide association analysis.

A whole-genome scan using single marker association was used to detect chromosome regions associated with seven female fertility traits in Finnish Ayrshire dairy cattle. The phenotypic data consisted of de-regressed estimated breeding values for 340 bulls which were estimated using a single trait model. Genotypes were obtained with the Illumina BovineSNP50 panel and a total of 35 630 informative, high-quality single nucleotide polymorphism (SNP) markers were used. The association analysis was performed using a mixed-model approach which fitted a fixed effect for each SNP and a random polygenic effect. We detected eleven genome-wide significant associations on eight different chromosomes. With at least chromosome-wise significance after Bonferroni correction, sixteen SNPs on nine chromosomes showed significant associations with one or more fertility traits. The results confirmed quantitative trait loci on three chromosomes (1, 2 and 20) for fertility traits previously reported for the same breed and one on chromosome four previously detected in Holstein cattle.

A truncated retrotransposon disrupts the GRM1 coding sequence in Coton de Tulear dogs with Bandera's neonatal ataxia.

BACKGROUND: Bandera's neonatal ataxia (BNAt) is an autosomal recessive cerebellar ataxia that affects members of the Coton de Tulear dog breed. OBJECTIVE: To identify the mutation that causes BNAt. ANIMALS: The study involved DNA from 112 Cotons de Tulear (including 15 puppies with signs of BNAt) and 87 DNA samples from dogs of 12 other breeds. METHODS: The BNAt locus was mapped with a genome-wide association study (GWAS). The coding exons of positional candidate gene GRM1, which encodes metabotropic glutamate receptor 1, were polymerase chain reaction (PCR)-amplified and resequenced. A 3-primer PCR assay was used to genotype individual dogs for a truncated retrotransposon inserted into exon 8 of GRM1. RESULTS: The GWAS indicated that the BNAt locus was in a canine chromosome 1 region that contained candidate gene GRM1. Resequencing this gene from BNAt-affected puppies indicated that exon 8 was interrupted by the insertion of a 5'-truncated retrotransposon. All 15 BNAt-affected puppies were homozygous for the insert, whereas all other Cotons de Tulear were heterozygotes (n = 43) or homozygous (n = 54) for the ancestral allele. None of the 87 dogs from 12 other breeds had the insertion allele. CONCLUSIONS AND CLINICAL IMPORTANCE: BNAt is caused by a retrotransposon inserted into exon 8 of GRM1. A DNA test for the GRM1 retrotransposon insert can be used for genetic counseling and to confirm the diagnosis of BNAt.

Characterization and mode of inheritance of a paroxysmal dyskinesia in Chinook dogs.

BACKGROUND: Paroxysmal dyskinesias are episodes of abnormal, involuntary movement or muscle tone, distinguished from seizures by the character of the episode and lack of seizure activity on ictal EEG. HYPOTHESIS: Paroxysmal dyskinesia is an inherited, autosomal recessive disorder in Chinook dogs. ANIMALS: Families of Chinook dogs with paroxysmal dyskinesia. METHODS: Pedigrees and medical histories were reviewed for 299 Chinook dogs. A family of 51 dogs was used for analysis. Episodes were classified as seizures, paroxysmal dyskinesia, or unknown, and segregation analysis was performed. RESULTS: Paroxysmal dyskinesia was identified in 16 of 51 dogs and characterized by an inability to stand or ambulate, head tremors, and involuntary flexion of 1 or multiple limbs, without autonomic signs or loss of consciousness. Episode duration varied from minutes to an hour. Inter-ictal EEGs recorded on 2 dogs with dyskinesia were normal. Three dogs with dyskinesia also had generalized tonic-clonic seizures. One of 51 dogs had episodes of undetermined type. Phenotype was unknown for 6 of 51 dogs, and 28 dogs were unaffected. Segregation was consistent with an autosomal recessive trait. CONCLUSIONS AND CLINICAL IMPORTANCE: This movement disorder is prevalent in the Chinook breed, and consistent with a partially penetrant autosomal recessive or polygenic trait. Insufficient evidence exists for definitive localization; episodes may be of basal nuclear origin, but atypical seizures and muscle membrane disorders remain possible etiologies. The generalized seizures may be a variant phenotype of the same mutation that results in dyskinesia, or the 2 syndromes may be independent.

A genome scan for quantitative trait loci influencing carcass, post-natal growth and reproductive traits in commercial Angus cattle.

To gain insight into the number of loci of large effect that underlie variation in cattle, a quantitative trait locus (QTL) scan for 14 economically important traits was performed in two commercial Angus populations using 390 microsatellites, 11 single nucleotide polymorphisms (SNPs) and one duplication loci. The first population comprised 1769 registered Angus bulls born between 1955 and 2003, with Expected Progeny Differences computed by the American Angus Association. The second comprised 38 half-sib families containing 1622 steers with six post-natal growth and carcass phenotypes. Linkage analysis was performed by half-sib least squares regression with gridqtl or Bayesian Markov chain Monte Carlo analysis of complex pedigrees with loki. Of the 673 detected QTL, only 118 have previously been reported, reflecting both the conservative approach to QTL reporting in the literature, and the more liberal approach taken in this study. From 33 to 71% of the genetic variance and 35 to 56% of the phenotypic variance in each trait was explained by the detected QTL. To analyse the effects of 11 SNPs and one duplication locus within candidate genes on each trait, a single marker analysis was performed by fitting an additive allele substitution model in both mapping populations. There were 53 associations detected between the SNP/duplication loci and traits with -log(10) P(nominal) ≥ 4.0, where each association explained 0.92% to 4.4% of the genetic variance and 0.01% to 1.86% of the phenotypic variance. Of these associations, only six SNP/duplication loci were located within 8 cM of a QTL peak for the trait, with two being located at the QTL peak: SST_DG156121:c.362A>G for ribeye muscle area and TG_X05380:c.422C>T for calving ease. Strong associations between several SNP/duplication loci and trait variation were obtained in the absence of any detected linked QTL. However, we reject the causality of several commercialized DNA tests, including an association between TG_X05380:c.422C>T and marbling in Angus cattle.

Genomic tools for characterizing monogenic and polygenic traits in ruminants--using the bovine as an example.

Next generation sequencing platforms have democratized genome sequencing. Large genome centers are no longer required to produce genome sequences costing millions. A few lanes of paired-end sequence on an Illumina Genome Analyzer, costing < $10,000, will produce more sequence than generated only a few years ago to produce the human and cow assemblies. The de novo assembly of large numbers of short reads into a high-quality whole-genome sequence is now technically feasible and will allow the whole genome sequencing and assembly of a broad spectrum of ruminant species. Next-generation sequencing instruments are also proving very useful for transcriptome or resequencing projects in which the entire RNA population produced by a tissue, or the entire genomes of individual animals are sequenced, and the produced reads are aligned to a reference assembly. We have used this strategy to examine gene expression differences in tissues from cattle differing in feed efficiency, to perform genome-wide single nucleotide polymorphism discovery for the construction of ultrahigh-density genotyping assays, and in combination with genome-wide association analysis, for the identification of mutations responsible for Mendelian diseases. The new 800K SNP bovine genotyping assays possess the resolution to map trait associations to the locations of individual genes and the 45 million polymorphisms identified in > 180X genome sequence coverage on over 200 animals can be queried to identify the polymorphisms present within positional candidate genes. These new tools should rapidly allow the identification of genes and mutations underlying variation in cattle production and reproductive traits.

Selection of single-nucleotide polymorphisms and quality of genotypes used in genomic evaluation of dairy cattle in the United States and Canada.

Nearly 57,000 single-nucleotide polymorphisms (SNP) genotyped with the Illumina BovineSNP50 BeadChip (Illumina Inc., San Diego, CA) were investigated to determine usefulness of the associated SNP for genomic prediction. Genotypes were obtained for 12,591 bulls and cows, and SNP were selected based on 5,503 bulls with genotypes from a larger set of SNP. The following SNP were deleted: 6,572 that were monomorphic, 3,213 with scoring problems (primarily because of poor definition of clusters and excess number of clusters), and 3,649 with a minor allele frequency of <2%. Number of SNP for each minor allele frequency class (> or =2%) was fairly uniform (777 to 1,004). For 5 contiguous SNP assigned to chromosome 7, no bulls were heterozygous, which indicated that those SNP are actually on the nonpseudoautosomal portion of the X chromosome. Another 178 SNP that were not assigned to a chromosome but that had many fewer heterozygotes than expected were also assigned to the X chromosome. Existence of Hardy-Weinberg equilibrium was investigated by comparing observed with expected heterozygosity. For 11 SNP, the observed percentage of heterozygous individuals differed from the expected by >15%; therefore, those SNP were deleted. For 2,628 SNP, the genotype at another SNP was highly correlated (i.e., genotypes were identical for >99.5% of bulls), and those were deleted. After edits, 40,874 SNP remained. A parent-progeny conflict was declared when the genotypes were alternate homozygotes. Mean number of conflicts was 2.3 when pedigree was correct and 2,411 when it was incorrect. The sire was genotyped for >93% of animals. Maternal grandsire genotype was similarly checked; however, because alternate homozygotes could be valid, a conflict threshold of 16% was used to indicate a need for further investigation. Genotyping consistency was investigated for 21 bulls genotyped twice with differences primarily from SNP that were not scored in one of the genotypes. Concordance for readable SNP was extremely high (99.96-100%). Thousands of SNP that were polymorphic in Holsteins were monomorphic in Jerseys or Brown Swiss, which indicated that breed-specific SNP sets are required or that all breeds need to be considered in the SNP selection process. Genotypes from the Illumina BovineSNP50 BeadChip are of high accuracy and provide the basis for genomic evaluations in the United States and Canada.

A whole genome association analysis identifies loci associated with Mycobacterium avium subsp. paratuberculosis infection status in US holstein cattle.

The purpose of this study was to identify loci associated with Mycobacterium avium subspecies paratuberculosis (Map) infection status in US Holsteins using the Illumina BovineSNP50 BeadChip whole genome single nucleotide polymorphism (SNP) assay. Two hundred forty-five cows from dairies in New York, Pennsylvania and Vermont enrolled in longitudinal herd studies between January 1999 and November 2007 were assessed for the presence of Map in both faecal and tissue samples. An animal was considered tissue infected if any sample contained at least one colony forming unit of Map per gram of tissue (CFU/g) and the same definition was employed for faecal samples. Each animal was genotyped with the Illumina BovineSNP50 BeadChip and after quality assurance filtering, 218 animals and 45 683 SNPs remained. We sought to identify loci associated with four different case/control classifications: presence of Map in the tissue, presence of Map in faeces, presence of Map in both tissue and faeces and presence of Map in tissue but not faeces. A case-control genome wide association study was conducted to test the four different classifications of Map infection status (cases) when compared with a Map-negative control group (control). Regions on chromosomes 1, 5, 7, 8, 16, 21 and 23 were identified with moderate significance (P < 5 x 10(-5)). Two regions, one on chromosome 3 (near EDN2) and another on chromosome 9 (no positional gene candidates), were identified with a high level of association to the presence of Map in tissue and both tissue and faeces respectively (P < 5 x 10(-7), genome-wide Bonferonni P < 0.05).