Genetic structure and admixture in sheep from terminal breeds in the United States.

Selection for performance in diverse production settings has resulted in variation across sheep breeds worldwide. Although sheep are an important species to the United States, the current genetic relationship among many terminal sire breeds is not well characterized. Suffolk, Hampshire, Shropshire and Oxford (terminal) and Rambouillet (dual purpose) sheep (n = 248) sampled from different flocks were genotyped using the Applied Biosystems Axiom Ovine Genotyping Array (50K), and additional Shropshire sheep (n = 26) using the Illumina Ovine SNP50 BeadChip. Relationships were investigated by calculating observed heterozygosity, inbreeding coefficients, eigenvalues, pairwise Wright's FST estimates and an identity by state matrix. The mean observed heterozygosity for each breed ranged from 0.30 to 0.35 and was consistent with data reported in other US and Australian sheep. Suffolk from two different regions of the United States (Midwest and West) clustered separately in eigenvalue plots and the rectangular cladogram. Further, divergence was detected between Suffolk from different regions with Wright's FST estimate. Shropshire animals showed the greatest divergence from other terminal breeds in this study. Admixture between breeds was examined using admixture, and based on cross-validation estimates, the best fit number of populations (clusters) was K = 6. The greatest admixture was observed within Hampshire, Suffolk, and Shropshire breeds. When plotting eigenvalues, US terminal breeds clustered separately in comparison with sheep from other locations of the world. Understanding the genetic relationships between terminal sire breeds in sheep will inform us about the potential applicability of markers derived in one breed to other breeds based on relatedness.

Neuropathic pain in a primary care electronic health record database.

BACKGROUND: Neuropathic pain (NP) is common in the adult population but is difficult to study in electronic health record (EHR) databases because it is a symptom rather than a pathologic diagnosis. The first step in studying NP in EHR databases is to develop methods for identifying patients with NP. The objectives of this study were to develop estimates of the prevalence of NP among patients in a primary care EHR database and describe these patients' demographic characteristics and health-care utilization. METHODS: This was a retrospective cohort study of de-identified data from a 5-year period (2005-2010) from 23 general practitioners (GPs) in 10 primary care practices in southwestern Ontario, Canada. International Classification of Diseases version 9 (ICD-9) diagnostic codes and medication prescriptions were used to identify patients with certain and probable NP. RESULTS: Different methods produced prevalence estimates ranging from 1.5% (for certain NP in the epidemiologically rigorous period cohort) to 11.2% (for certain NP + probable NP in the more inclusive database cohort). Patients in the NP groups had more GP visits, specialist referrals and analgesic prescriptions than patients without NP. CONCLUSION: This study represents a step towards being able to utilize EHR databases to study NP by proposing methods to identify patients with certain and probable NP in a primary care EHR database. Validation against a gold standard is the next step.

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.

Identification of loci associated with tolerance to Johne's disease in Holstein cattle.

Johne's disease, caused by Mycobacterium avium subspecies paratuberculosis (Map), is a fatal disease in cattle. The objective of this study was to identify loci associated with tolerance in cows infected with Map. Tolerance was defined as a cow's fitness at a given level of Map infection intensity. Fitness was measured by Map faecal cultures, and Map infection intensity was measured by culturing four gut tissues. The quantitative phenotype of tolerance was defined by numerical indexes of cultures of peak (peak tolerance, PT) and average (average tolerance, AT) faecal and tissue Map from 245 Holstein cows. The categorical phenotype was defined as: ≥ 100 cfu Map tissue infection, and faecal shedding ≥ 75 cfu (intolerant) or <10 cfu (tolerant cows). In 94 cows, Map was identified in ≥ 1 tissue, including 44 cows with ≥ 100 Map tissue cfu and 36 with ≥ 1 faecal cfu. A genome-wide association analysis was performed after filtering, leaving genotypes for 45,789 SNPs in 90 animals for the quantitative phenotype and 16 cases and 25 controls for the categorical analysis of tolerance. rs41748405:A>C (BTA15) was associated with PT (P = 1.12 × 10(-7)) and AT (P = 2.17 × 10(-6)). Associations were identified with PT and adjacent SNPs ss61512613:A>G and ss61530518:A>G (BTA6) (P < 3.0 × 10(-5)), and with AT for ss61469568:A>G (BTA 2) (P = 3.3 × 10(-5)) and ss86284768:A>G (BTA1) (P = 3.31 × 10(-5)). For the categorical phenotype, an association was found with ss8632653:A>G (BTA6) (P < 5.0 × 10(-5)). This is the first study to identify loci associated with tolerance to Johne's disease.

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.

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).

Linkage disequilibrium and signatures of selection on chromosomes 19 and 29 in beef and dairy cattle.

The objective of this study was to quantify the extent of linkage disequilibrium (LD) on bovine chromosomes 19 and 29 and to study the pattern of selection signatures in beef and dairy breeds (Angus and Holstein) of Bos taurus. The extent of LD was estimated for 370 and 186 single nucleotide polymorphism markers on BTA19 and 29 respectively using the square of the correlation coefficient (r(2)) among alleles at pairs of loci. A comparison of the extent of LD found that the decline of LD followed a similar pattern in both breeds. We observed long-range LD and found that LD dissipates to background levels at a locus separation of about 20 Mb on both chromosomes. Along each chromosome, patterns of LD were variable in both breeds. We find that a minimum of 30 000 informative and evenly spaced markers would be required for whole-genome association studies in cattle. In addition, we have identified chromosomal regions that show some evidence of selection for economically important traits in Angus and Holstein cattle. The results of this study are of importance for the design and application of association studies.

Construction of bovine whole-genome radiation hybrid and linkage maps using high-throughput genotyping.

High-density whole-genome maps are essential for ordering genes or markers and aid in the assembly of genome sequence. To increase the density of markers on the bovine radiation hybrid map, and hence contribute to the assembly of the bovine genome sequence, an Illumina BeadStation was used to simultaneously type large numbers of markers on the Roslin-Cambridge 3000 rad bovine-hamster whole-genome radiation hybrid panel (WGRH3000). In five multiplex reactions, 6738 sequence tagged site (STS) markers were successfully typed on the WGRH3000 panel DNA. These STSs harboured SNPs that were developed as a result of the bovine genome sequencing initiative. Typically, the most time consuming and expensive part of creating high-density radiation hybrid (RH) maps is genotyping the markers on the RH panel with conventional approaches. Using the method described in this article, we have developed a high-density whole-genome RH map with 4690 loci and a linkage map with 2701 loci, with direct comparison to the bovine whole-genome sequence assembly (Btau_2.0) in a fraction of the time it would have taken with conventional typing and genotyping methods.

Assessment of complications of distal radius fractures and development of a complication checklist.

The purposes of this study were to determine the overall incidence of distal radius fracture (DRF) complications, determine the incidence and types of DRF complications in a consecutive cohort of 250 patients with DRFs, describe DRF complications reported by patients compared with those reported by physicians, and formulate a DRF complication checklist to improve recording of DRF complications. We found that the overall complication rates vary widely (6% to 80%). Physician-reported complication data were collected for 236 patients, and a physician-reported complication rate of 27% was determined. A patient-reported complication rate of 21% was found for 207 patients whose patient-reported data were collected. We also noted that patients and physicians assess DRF complications differently: patients are more focused on symptoms than diagnoses. A DRF complication checklist was developed to improve prospective data collection. The checklist includes a classification for all DRF complications and allows for assessment of severity of each complication.