Search

Objective: Evaluate a precision medicine approach to confirm a tentative diagnosis of fatal malignant hyperthermia (MH) in isoflurane-anesthetized pet dogs by identifying novel risk variants in known MH susceptibility genes., Study Design: Retrospective case series., Animals: A male Pit Bull mix aged 7 years (case #1), a male Golden Retriever aged 12 months (case #2) and the dam and sire of case #2., Methods: Available case histories and medical records were reviewed. Missense variants in MH susceptibility genes RYR-1, CACNA1S and STAC3 (case #2 only) were identified by next-generation sequencing of DNA from each case and the parents of case #2 with confirmation by Sanger sequencing. The pathogenicity of variants was evaluated by multiple in silico approaches., Results: Both cases demonstrated clinical signs during isoflurane anesthesia consistent with volatile anesthetic-induced MH, including tachypnea, tachycardia, severe hyperthermia and muscle rigidity. Despite whole body cooling and other treatments, both dogs died after cardiac arrest within 15 minutes of detecting hyperthermia. Gene sequencing identified novel missense RYR-1 variants in case #1 (p.Gly2375Arg) and case #2 (p.Pro152Leu). Both variants were likely pathogenic based on multiple criteria, including gene location, amino acid alteration and population allele frequency. The case #1 variant was identical to a known human diagnostic MH variant (p.Gly2375Arg). Neither parent of case #2 had the case #2 variant, indicating this variant was not inherited, but arose de novo in a germ cell of either parent or early in embryogenesis. Whole genome sequence analysis confirmed parentage. Two missense variants were identified in CACNA1S. Both variants were considered nonpathogenic. No variants were identified in STAC3., Conclusions and Clinical Relevance: Like humans, MH susceptibility in dogs is associated with different rare variants located in pathogenic hotspots in the RYR-1 gene. Next-generation sequencing is a useful tool to assist in the definitive diagnosis of MH in dogs., (Copyright © 2024 The Author(s). Published by Elsevier Ltd.. All rights reserved.)

Background/Objectives: The dairy industry relies on reproductive efficiency to maintain efficient milk production. Spontaneous abortion (SA), defined as pregnancy loss between gestation days 42 and 260, occurred in 4.5% of the artificially inseminated (AI) Holstein heifers and 31.6% of the embryo transfer (ET) recipient Holstein heifers that received in vitro-produced frozen embryos on a single dairy farm in Idaho. Methods: A genome-wide association analysis (GWAA) was performed to identify the associations (FDR p < 0.05) with SA in heifers that were bred by AI (1351 controls that delivered at term and 63 cases that aborted) that conceived following the first insemination, as well as in 59 controls and 273 cases of ET recipient heifers pregnant from the first ET. Results: There were 216 loci and 413 positional candidate genes associated (FDR p < 0.05) with SA in the heifers bred by AI in a recessive model and no loci associated with SA in the ET recipients. Conclusions: The identification of loci associated with SA in the heifers bred by AI may be used to reduce fetal loss through genomic selection.

Reference genomes of cattle and sheep have lacked contiguous assemblies of the sex-determining Y chromosome. Here, we assemble complete and gapless telomere to telomere (T2T) Y chromosomes for these species. We find that the pseudo-autosomal regions are similar in length, but the total chromosome size is substantially different, with the cattle Y more than twice the length of the sheep Y. The length disparity is accounted for by expanded ampliconic region in cattle. The genic amplification in cattle contrasts with pseudogenization in sheep suggesting opposite evolutionary mechanisms since their divergence 19MYA. The centromeres also differ dramatically despite the close relationship between these species at the overall genome sequence level. These Y chromosomes have been added to the current reference assemblies in GenBank opening new opportunities for the study of evolution and variation while supporting efforts to improve sustainability in these important livestock species that generally use sire-driven genetic improvement strategies., (© 2024. This is a U.S. Government work and not under copyright protection in the US; foreign copyright protection may apply.)

Breastfeeding provides many health benefits, but its impact on respiratory health remains unclear. This study addresses the complex and dynamic nature of the mother-milk-infant triad by investigating maternal genomic factors regulating human milk oligosaccharides (HMOs), and their associations with respiratory health among human milk-fed infants. Nineteen HMOs are quantified from 980 mothers of the CHILD Cohort Study. Genome-wide association studies identify HMO-associated loci on chromosome 19p13.3 and 19q13.33 (lowest P = 2.4e-118), spanning several fucosyltransferase (FUT) genes. We identify novel associations on chromosome 3q27.3 for 6'-sialyllactose (P = 2.2e-9) in the sialyltransferase (ST6GAL1) gene. These, plus additional associations on chromosomes 7q21.32, 7q31.32 and 13q33.3, are replicated in the independent INSPIRE Cohort. Moreover, gene-environment interaction analyses suggest that fucosylated HMOs may modulate overall risk of recurrent wheeze among preschoolers with variable genetic risk scores (P < 0.01). Thus, we report novel genetic factors associated with HMOs, some of which may protect the respiratory health of children., (© 2024. The Author(s).)

Heifer conception rate to the first service (HCR1) is defined as the number of heifers that become pregnant to the first breeding service compared to the heifers bred. This study aimed to identify loci associated and gene sets enriched for HCR1 for heifers that were bred by artificial insemination (AI, n = 2829) or were embryo transfer (ET, n = 2086) recipients, by completing a genome-wide association analysis and gene set enrichment analysis using SNP data (GSEA-SNP). Three unique loci, containing four positional candidate genes, were associated ( p < 1 × 10 -5 ) with HCR1 for ET recipients, while the GSEA-SNP identified four gene sets (NES ≥ 3) and sixty-two leading edge genes (LEGs) enriched for HCR1. While no loci were associated with HCR1 bred by AI, one gene set and twelve LEGs were enriched (NES ≥ 3) for HCR1 with the GSEA-SNP. This included one gene ( PKD2 ) shared between HCR1 AI and ET services. Identifying loci associated or enriched for HCR1 provides an opportunity to use them as genomic selection tools to facilitate the selection of cattle with higher reproductive efficiency, and to better understand embryonic loss.

Ovine footrot is an infectious disease with important contributions from Dichelobacter nodosus and Fusobacterium necrophorum . Footrot is characterized by separation of the hoof from underlying tissue, and this causes severe lameness that negatively impacts animal wellbeing, growth, and profitability. Large economic losses result from lost production as well as treatment costs, and improved genetic tools to address footrot are a valuable long-term goal. Prior genetic studies had examined European wool sheep, but hair sheep breeds such as Katahdin and Blackbelly have been reported to have increased resistance to footrot, as well as to intestinal parasites. Thus, footrot condition scores were collected from 251 U.S. sheep including Katahdin, Blackbelly, and European-influenced crossbred sheep with direct and imputed genotypes at OvineHD array (>500,000 single nucleotide polymorphism) density. Genome-wide association was performed with a mixed model accounting for farm and principal components derived from animal genotypes, as well as a random term for the genomic relationship matrix. We identified three genome-wide significant associations, including SNPs in or near GBP6 and TCHH . We also identified 33 additional associated SNPs with genome-wide suggestive evidence, including a cluster of 6 SNPs in a peak near the genome-wide significance threshold located near the glutamine transporter gene SLC38A1 . These findings suggest genetic susceptibility to footrot may be influenced by genes involved in divergent biological processes such as immune responses, nutrient availability, and hoof growth and integrity. This is the first genome-wide study to investigate susceptibility to footrot by including hair sheep and also the first study of any kind to identify multiple genome-wide significant associations with ovine footrot. These results provide a foundation for developing genetic tests for marker-assisted selection to improve resistance to ovine footrot once additional steps like fine mapping and validation are complete., Competing Interests: The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest., (Copyright © 2024 Cinar, Oliveira, Hadfield, Lichtenwalner, Brzozowski, Settlemire, Schoenian, Parker, Neibergs, Cockett and White.)

Pseudorabies virus (PRV)-the causative agent of Aujeszky's disease-was eliminated from commercial pig production herds in the United States (US) in 2004; however, PRV remains endemic among invasive feral swine ( Sus scrofa ). The circulation of PRV among abundant, widespread feral swine populations poses a sustained risk for disease spillover to production herds. Risk-based surveillance has been successfully implemented for PRV in feral swine populations in the US. However, understanding the role of host genetics in infection status may offer new insights into the epidemiology and disease dynamics of PRV that can be applied to management strategies. Genetic mechanisms underlying host susceptibility to PRV are relatively unknown; therefore, we sought to identify genomic regions associated with PRV infection status among naturally infected feral swine using genome-wide association studies (GWAS) and gene set enrichment analysis of single nucleotide polymorphism data (GSEA-SNP). Paired serological and genotypic data were collected from 6,081 feral swine distributed across the invaded range within the contiguous US. Three complementary study populations were developed for GWAS: 1) comprehensive population consisting of feral swine throughout the invaded range within the contiguous US; 2) population of feral swine under high, but temporally variable PRV infection pressure; and 3) population of feral swine under temporally stable, high PRV infection pressure. We identified one intronic SNP associated with PRV infection status within candidate gene AKAP6 on autosome 7. Various gene sets linked to metabolic pathways were enriched in the GSEA-SNP. Ultimately, improving disease surveillance efforts in feral swine will be critical to further understanding of the role host genetics play in PRV infection status, helping secure the health of commercial pork production., Competing Interests: The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest., (Copyright © 2023 Bowden, Kiser, Miller, Buckley, Boggiatto, Giglio, Brown, Garrick, Neibergs, Piaggio, Speidel and Smyser.)

Background: Delayed postoperative hemorrhage (DEPOH) is an important health concern for Scottish deerhounds., Hypothesis/objectives: Identify genes associated with DEPOH in Scottish deerhounds., Animals: Two hundred sixty-nine privately owned Scottish deerhounds., Methods: Retrospective case-control study. DEPOH cases and controls were identified through an owner health survey. Genome-wide association analysis was performed using whole genome sequences from 8 cases and 17 controls. All cases and controls were genotyped for selected variants., Results: Of 269 dogs, 10 met inclusion and exclusion criteria for DEPOH, while 62 controls had undergone similar surgical procedures without DEPOH. Genome-wide association analysis identified a single locus on chromosome 9 spanning 40 genes. One of these genes (SERPINF2 encoding alpha-2 antiplasmin) was directly linked to the pathophysiology of DEPOH. The entire cohort was genotyped for a missense SERPINF2 variant (c.605 C>T; p.A202V). Compared to dogs with the reference C/C genotype, the likelihood of DEPOH was significantly higher for dogs with the T/T genotype (odds ratio [OR] = 1235; 95% confidence interval [CI] = 23-6752; P = 0.0005) and with the C/T genotype (OR = 28; 95% CI = 1.4-542; P = 0.03)., Conclusions and Clinical Importance: SERPINF2 is associated with DEPOH in Scottish deerhounds. Genetic testing might be able to identify dogs that are susceptible to DEPOH., (© 2023 The Authors. Journal of Veterinary Internal Medicine published by Wiley Periodicals LLC on behalf of American College of Veterinary Internal Medicine.)

Monocytes are a core component of the immune system that arise from bone marrow and differentiate into cells responsible for phagocytosis and antigen presentation. Their derivatives are often responsible for the initiation of the adaptive immune response. Monocytes and macrophages are central in both controlling and propagating infectious diseases such as infection by Coxiella burnetii and small ruminant lentivirus in sheep. Genotypes from 513 Rambouillet, Polypay, and Columbia sheep (Ovis aries) were generated using the Ovine SNP50 BeadChip. Of these sheep, 222 animals were subsequently genotyped with the Ovine Infinium® HD SNP BeadChip to increase SNP coverage. Data from the 222 HD genotyped sheep were combined with the data from an additional 258 unique sheep to form a 480-sheep reference panel; this panel was used to impute the low-density genotypes to the HD genotyping density. Then, a genome-wide association analysis was conducted to identify loci associated with absolute monocyte counts from blood. The analysis used a single-locus mixed linear model implementing EMMAX with age and ten principal components as fixed effects. Two genome-wide significant peaks (p < 5x10-7) were identified on chromosomes 9 and 1, and ten genome-wide suggestive peaks (p < 1x10-5) were identified on chromosomes 1, 2, 3, 4, 9, 10, 15, and 16. The identified loci were within or near genes including KCNK9, involved into cytokine production, LY6D, a member of a superfamily of genes, some of which subset monocyte lineages, and HMGN1, which encodes a chromatin regulator associated with myeloid cell differentiation. Further investigation of these loci is being conducted to understand their contributions to monocyte counts. Investigating the genetic basis of monocyte lineages and numbers may in turn provide information about pathogens of veterinary importance and elucidate fundamental immunology., Competing Interests: The authors have declared that no competing interests exist.

Bovine coronavirus (BCoV) is associated with respiratory and enteric infections in both dairy and beef cattle worldwide. It is also one of a complex of pathogens associated with bovine respiratory disease (BRD), which affects millions of cattle annually. The objectives of this study were to identify loci and heritability estimates associated with BCoV infection and BRD in dairy calves and feedlot cattle. Dairy calves from California ( n = 1,938) and New Mexico ( n = 647) and feedlot cattle from Colorado ( n = 915) and Washington ( n = 934) were tested for the presence of BCoV when classified as BRD cases or controls following the McGuirk scoring system. Two comparisons associated with BCoV were investigated: (1) cattle positive for BCoV (BCoV + ) were compared to cattle negative for BCoV (BCoV - ) and (2) cattle positive for BCoV and affected with BRD (BCoV + BRD + ) were compared to cattle negative for BCoV and BRD (BCoV - BRD - ). The Illumina BovineHD BeadChip was used for genotyping, and genome-wide association analyses (GWAA) were performed using EMMAX (efficient mixed-model association eXpedited). The GWAA for BCoV + identified 51 loci ( p < 1 × 10 -5 ; 24 feedlot, 16 dairy, 11 combined) associated with infection with BCoV. Three loci were associated with BCoV + across populations. Heritability estimates for BCoV + were 0.01 for dairy, 0.11 for feedlot cattle, and 0.03 for the combined population. For BCoV + BRD + , 80 loci ( p < 1 × 10 -5 ; 26 feedlot, 25 dairy, 29 combined) were associated including 14 loci across populations. Heritability estimates for BCoV + BRD + were 0.003 for dairy, 0.44 for feedlot cattle, and 0.07 for the combined population. Several positional candidate genes associated with BCoV and BRD in this study have been associated with other coronaviruses and respiratory infections in humans and mice. These results suggest that selection may reduce susceptibility to BCoV infection and BRD in cattle., Competing Interests: The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest., (Copyright © 2021 Kiser and Neibergs.)

Human milk oligosaccharides (HMO), the third most abundant component of human milk, are thought to be important contributors to infant health. Studies have provided evidence that geography, stage of lactation, and Lewis and secretor blood groups are associated with HMO profile. However, little is known about how variation across the genome may influence HMO composition among women in various populations. In this study, we performed genome-wide association analyses of 395 women from 8 countries to identify genetic regions associated with 19 different HMO. Our data support FUT2 as the most significantly associated (P < 4.23 -9 to P < 4.5 -70 ) gene with seven HMO and provide evidence of balancing selection for FUT2. Although polymorphisms in FUT3 were also associated with variation in lacto-N-fucopentaose II and difucosyllacto-N-tetrose, we found little evidence of selection on FUT3. To our knowledge, this is the first report of the use of genome-wide association analyses on HMO., (Copyright © 2021 The Authors. Published by Elsevier Inc. All rights reserved.)

Small ruminant lentivirus (SRLV) causes Maedi-Visna or Ovine Progressive Pneumonia in sheep and creates insidious livestock production losses. This retrovirus is closely related to human immunodeficiency virus and currently has no vaccines or cure. Genetic marker assisted selection for sheep disease resiliency presents an attractive management solution. Previously, we identified a region containing a cluster of zinc finger genes that had association with ovine SRLV proviral concentration. Trait-association analysis validated a small insertion/deletion variant near ZNF389 (rs397514112) in multiple sheep breeds. In the current study, 543 sheep from two distinct populations were genotyped at 34 additional variants for fine mapping of the regulatory elements within this locus. Variants were selected based on ChIP-seq annotation data from sheep alveolar macrophages that defined active cis-regulatory elements predicted to influence zinc finger gene expression. We present a haplotype block of variants within regulatory elements that have improved associations and larger effect sizes (up to 4.7-fold genotypic difference in proviral concentration) than the previously validated ZNF389 deletion marker. Hypotheses for the underlying causal mutation or mutations are presented based on changes to in silico transcription factor binding sites. These variants offer alternative markers for selective breeding and are targets for future functional mutation assays.

Bovine paratuberculosis, caused by Mycobacterium avium subspecies paratuberculosis ( MAP ), continues to impact the dairy industry through increased morbidity, mortality, and lost production. Although genome-wide association analyses (GWAAs) have identified loci associated with susceptibility to MAP , limited progress has been made in identifying mutations that cause disease susceptibility. A 235-kb region on Bos taurus chromosome 3 (BTA3), containing a 70-kb haplotype block surrounding endothelin 2 ( EDN2 ), has previously been associated with the risk of MAP infection. EDN2 is highly expressed in the gut and is involved in intracellular calcium signaling and a wide array of biological processes. The objective of this study was to identify putative causal mutations for disease susceptibility in the region surrounding EDN2 in Holstein and Jersey cattle. Using sequence data from 10 Holstein and 10 Jersey cattle, common variants within the 70-kb region containing EDN2 were identified. A custom SNP genotyping array fine-mapped the region using 221 Holstein and 51 Jersey cattle and identified 17 putative causal variants ( P < 0.01) located in the 5' region of EDN2 and a SNP in the 3' UTR ( P = 0.00009) associated with MAP infection. MicroRNA interference assays, mRNA stability assays, and electrophoretic mobility shift assays were performed to determine if allelic changes at each SNP resulted in differences in EDN2 stability or expression. Two SNPs [ rs109651404 (G/A) and rs110287192 (G/T)] located within the promoter region of EDN2 displayed differential binding affinity for transcription factors in binding sequences harboring the alternate SNP alleles. The luciferase reporter assay revealed that the transcriptional activity of the EDN2 promoter was increased ( P < 0.05) with the A allele for rs109651404 and the G allele for rs110287192 . These results suggest that the variants rs109651404 and rs110287192 are mutations that alter transcription and thus may alter susceptibility to MAP infection in Holstein and Jersey cattle., Competing Interests: The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest., (Copyright © 2021 Kiser, Wang, Zanella, Scraggs, Neupane, Cantrell, Van Tassell, White, Taylor and Neibergs.)

Bovine respiratory disease (BRD) is the costliest disease affecting the cattle industry globally. Orthomyxoviruses, influenza C virus (ICV) and influenza D virus (IDV) have recently been implicated to play a role in BRD. However, there are contradicting reports about the association of IDV and ICV to BRD. Using the largest cohort study (cattle, n = 599) to date we investigated the association of influenza viruses in cattle with BRD. Cattle were scored for respiratory symptoms and pooled nasal and pharyngeal swabs were tested for bovine viral diarrhea virus, bovine herpesvirus 1, bovine respiratory syncytial virus, bovine coronavirus, ICV and IDV by real-time PCR. Cattle that have higher viral loads of IDV and ICV also have greater numbers of co-infecting viruses than controls. More strikingly, 2 logs higher IDV viral RNA in BRD-symptomatic cattle that are co-infected animals than those infected with IDV alone. Our results strongly suggest that ICV and IDV may be significant contributors to BRD., (Copyright © 2020 Elsevier Inc. All rights reserved.)

Paratuberculosis (pTB), also known as Johne's disease (JD), is a contagious, chronic, and granulomatous inflammatory disease of the intestines of ruminants which is caused by Mycobacterium avium subsp. paratuberculosis (MAP) infection, resulting in billions of dollars in economic losses worldwide. Since, currently, no effective cure is available for MAP infection, it is important to explore the genetic variants that affect the host MAP susceptibility. The aim of this study was to analyze a potential association between EDN2 synonymous gene mutations (rs110287192, rs109651404 and rs136707411), that modifies susceptibility to pTB. EDN2 rs110287192, rs109651404 and rs136707411 mutations were genotyped in 68 infected and 753 healthy animals from East Anatolian Red crossbred, Anatolian Black crossbred and Holstein breed cattle by using Custom TaqMan SNP Genotyping Assays. For pTB status, serum antibody levels S/P ≥ 1.0 were assessed in carriers of the different EDN2 genotypes. EDN2 rs110287192 mutation showed a significant association with bovine pTB (adj. p < 0.05). For rs110287192 locus, the odd ratios for GG and TG genotypes versus TT genotypes were 1.73; (95% CI = 0.34-8.59) and 0.53 (95% CI = 0.12-2.37) respectively, which indicated that proportion of TG heterozygotes were significantly higher in control animals as compared to pTB animals. On the other hand, while rs136707411 mutation showed a suggestive association with pTB status in the examined cattle population (nominal p < 0.05); no association was detected between rs109651404 genotypes and pTB status. Selecting animals against rs110287192-GG genotype may decrease the risk of pTB in cattle of the Bos taurus taurus subspecies., Competing Interests: The authors have declared that no competing interests exist.

Heifer conception rate (HCR) is defined as the percentage of inseminated heifers that become pregnant at each service. The genome-wide association analyses in this study focused on identifying the loci associated with Holstein heifer ( n = 2013) conception rate at first service (HCR1) and the number of times bred (TBRD) to achieve a pregnancy. There were 348 unique loci associated ( p < 5 × 10 -8 ) with HCR1 and 615 unique loci associated ( p < 5 × 10 -8 ) with TBRD. The two phenotypes shared 302 loci, and 56 loci were validated in independent cattle populations. There were 52 transcription factor binding sites (TFBS) and 552 positional candidate genes identified in the HCR1- and TBRD-associated loci. The positional candidate genes and the TFBS associated with HCR1 and TBRD were used in the ingenuity pathway analysis (IPA). In the IPA, 11 pathways, 207 master regulators and 11 upstream regulators were associated ( p < 1.23 × 10 -5 ) with HCR1 and TBRD. The validated loci associated with both HCR1 and TBRD make good candidates for genomic selection and further investigations to elucidate the mechanisms associated with subfertility and infertility.

Background: The objectives of this study were to identify loci, positional candidate genes, gene-sets, and pathways associated with spontaneous abortion (SA) in cattle and compare these results with previous human SA studies to determine if cattle are a good SA model for humans. Pregnancy was determined at gestation day 35 for Holstein heifers and cows. Genotypes from 43,984 SNPs of 499 pregnant heifers and 498 pregnant cows that calved at full term (FT) were compared to 62 heifers and 28 cows experiencing SA. A genome-wide association analysis, gene-set enrichment analysis-single nucleotide polymorphism, and ingenuity pathway analysis were used to identify regions, pathways, and master regulators associated with SA in heifers, cows, and a combined population., Results: Twenty-three loci and 21 positional candidate genes were associated ( p < 1 × 10 -5 ) with SA and one of these (KIR3DS1) has been associated with SA in humans. Eight gene-sets (NES > 3.0) were enriched in SA and one was previously reported as enriched in human SA. Four master regulators ( p < 0.01) were associated with SA within two populations., Conclusions: One locus associated with SA was validated and 39 positional candidate and leading-edge genes and 2 gene-sets were enriched in SA in cattle and in humans.

Background: Subfertility is a major issue facing the dairy industry as the average US Holstein cow conception rate (CCR) is approximately 35%. The genetics underlying the physiological processes responsible for CCR, the proportion of cows able to conceive and maintain a pregnancy at each breeding, are not well characterized. The objectives of this study were to identify loci, positional candidate genes, and transcription factor binding sites (TFBS) associated with CCR and determine if there was a genetic correlation between CCR and milk production in primiparous Holstein cows. Cows were bred via artificial insemination (AI) at either observed estrus or timed AI and pregnancy status was determined at day 35 post-insemination. Additive, dominant, and recessive efficient mixed model association expedited (EMMAX) models were used in two genome-wide association analyses (GWAA). One GWAA focused on CCR at first service (CCR1) comparing cows that conceived and maintained pregnancy to day 35 after the first AI (n = 494) to those that were open after the first AI (n = 538). The second GWAA investigated loci associated with the number of times bred (TBRD) required for conception in cows that either conceived after the first AI (n = 494) or repeated services (n = 472)., Results: The CCR1 GWAA identified 123, 198, and 76 loci associated (P < 5 × 10 - 08 ) in additive, dominant, and recessive models, respectively. The TBRD GWAA identified 66, 95, and 33 loci associated (P < 5 × 10 - 08 ) in additive, dominant, and recessive models, respectively. Four of the top five loci were shared in CCR1 and TBRD for each GWAA model. Many of the associated loci harbored positional candidate genes and TFBS with putative functional relevance to fertility. Thirty-six of the loci were validated in previous GWAA studies across multiple breeds. None of the CCR1 or TBRD associated loci were associated with milk production, nor was their significance with phenotypic and genetic correlations to 305-day milk production., Conclusions: The identification and validation of loci, positional candidate genes, and TFBS associated with CCR1 and TBRD can be utilized to improve, and further characterize the processes involved in cattle fertility.

Background: Subfertility is one challenge facing the dairy industry as the average Holstein heifer conception rate (HCR), the proportion of heifers that conceive and maintain a pregnancy per breeding, is estimated at 55-60%. Of the loci associated with HCR, few have been validated in an independent cattle population, limiting their usefulness for selection or furthering our understanding of the mechanisms involved in successful pregnancy. Therefore, the objectives here were to identify loci associated with HCR: 1) to the first artificial insemination (AI) service (HCR1), 2) to repeated AI services required for a heifer to conceive (TBRD) and 3) to validate loci previously associated with fertility. Breeding and health records from 3359 Holstein heifers were obtained after heifers were bred by AI at observed estrus, with pregnancy determined at day 35 via palpation. Heifer DNA was genotyped using the Illumina BovineHD BeadChip, and genome-wide association analyses (GWAA) were performed with additive, dominant and recessive models using the Efficient Mixed Model Association eXpedited (EMMAX) method with a relationship matrix for two phenotypes. The HCR1 GWAA compared heifers that were pregnant after the first AI service (n = 497) to heifers that were open following the first AI service (n = 405), which included those that never conceived. The TBRD GWAA compared only those heifers which did conceive, across variable numbers of AI service (n = 712). Comparison of loci previously associated with fertility, HCR1 or TBRD were considered the same locus for validation when in linkage disequilibrium (D' > 0.7)., Results: The HCR1 GWAA identified 116, 187 and 28 loci associated (P < 5 × 10 - 8 ) in additive, dominant and recessive models, respectively. The TBRD GWAA identified 235, 362, and 69 QTL associated (P < 5 × 10 - 8 ) with additive, dominant and recessive models, respectively. Loci previously associated with fertility were in linkage disequilibrium with 22 loci shared with HCR1 and TBRD, 5 HCR1 and 19 TBRD loci., Conclusions: Loci associated with HCR1 and TBRD that have been identified and validated can be used to improve HCR through genomic selection, and to better understand possible mechanisms associated with subfertility.

Background: National genetic evaluations for disease resistance do not exist, precluding the genetic improvement of cattle for these traits. We imputed BovineHD genotypes to whole genome sequence for 2703 Holsteins that were cases or controls for Bovine Respiratory Disease and sampled from either California or New Mexico to construct and compare genomic prediction models. The sequence variation reference dataset comprised variants called for 1578 animals from Run 5 of the 1000 Bull Genomes Project, including 450 Holsteins and 29 animals sequenced from this study population. Genotypes for 9,282,726 variants with minor allele frequencies ≥5% were imputed and used to obtain genomic predictions in GEMMA using a Bayesian Sparse Linear Mixed Model., Results: Variation explained by markers increased from 13.6% using BovineHD data to 14.4% using imputed whole genome sequence data and the resolution of genomic regions detected as harbouring QTL substantially increased. Explained variation in the analysis of the combined California and New Mexico data was less than when data for each state were separately analysed and the estimated genetic correlation between risk of Bovine Respiratory Disease in California and New Mexico Holsteins was - 0.36. Consequently, genomic predictions trained using the data from one state did not accurately predict disease risk in the other state. To determine if a prediction model could be developed with utility in both states, we selected variants within genomic regions harbouring: 1) genes involved in the normal immune response to infection by pathogens responsible for Bovine Respiratory Disease detected by RNA-Seq analysis, and/or 2) QTL identified in the association analysis of the imputed sequence variants. The model based on QTL selected variants is biased but when trained in one state generated BRD risk predictions with positive accuracies in the other state., Conclusions: We demonstrate the utility of sequence-based and biology-driven model development for genomic selection. Disease phenotypes cannot be routinely recorded in most livestock species and the observed phenotypes may vary in their genomic architecture due to variation in the pathogen composition across environments. Elucidation of trait biology and genetic architecture may guide the development of prediction models with utility across breeds and environments.

Consumption of an adequate volume of high-quality colostrum soon after birth is critical for a calf's health. Few studies have focused on the genetics associated with colostrum production, even though several dairy herds in the United States have reported incidents of low to no colostrum production during the fall and winter seasons. The objectives of this study were to identify loci associated with quantity and quality of colostrum production in a herd of Jersey cattle (n = 345) and to identify potential positional candidate genes and/or transcription factor binding site motifs located near associated loci. Cattle that freshened between the months of October and December of 2016 at a single dairy were enrolled in the study and produced on average 3.03 kg of colostrum at their first milking. This study included 112 cattle genotyped with the GeneSeek GGP50k BeadChip and another 233 cattle previously genotyped with various other arrays. The 233 cattle genotyped at lower densities were imputed to the GGP50k BeadChip density using BEAGLE 4.1.1, and 2 genome-wide association analyses (GWAA) were conducted using an additive efficient mixed-model association expedited method with a genomic relationship matrix (EMMAX-GRM). The first GWAA investigated loci associated with colostrum quantity and identified 7 loci: 6 that were moderately associated (5 × 10-07 > P < 1 × 10-05) and 1 that was strongly associated (P < 5 × 10-07). The second GWAA investigated colostrum quality and identified 1 moderately (5 × 10-07 > P < 1 × 10-05) associated locus. Five loci harbored positional candidate genes which had functional relevance to colostrum production, and 1 locus located on BTA10 contained a transcription factor binding site motif for TFAP2A which has previously been linked to mammary gland development. Pseudoheritability estimates were moderate for colostrum quality (0.19 ± 0.06) and high for colostrum quantity (0.76 ± 0.11), suggesting that genomic selection for these traits would be possible. Diminished colostrum quantity or quality can have a significant impact on herd health and herd economics. The identification of loci, positional candidate genes, and transcription factor binding site motifs associated with colostrum production could be used in genomic selection to allow producers to select for cattle with good colostrum production, improving calf health, and reducing economic losses to the herd., (© The Author(s) 2018. Published by Oxford University Press on behalf of the American Society of Animal Science. All rights reserved. For permissions, please e-mail: journals.permissions@oup.com.)

A major unresolved issue is how the uterus influences infertility and subfertility in cattle. Serial embryo transfer was previously used to classify heifers as high-fertile (HF), subfertile (SF), or infertile (IF). To assess pregnancy loss, two in vivo-produced embryos were transferred into HF, SF, and IF heifers on day 7, and pregnancy outcome was assessed on day 17. Pregnancy rate was substantially higher in HF (71%) and SF (90%) than IF (20%) heifers. Elongating conceptuses were about twofold longer in HF than SF heifers. Transcriptional profiling detected relatively few differences in the endometrium of nonpregnant HF, SF, and IF heifers. In contrast, there was a substantial difference in the transcriptome response of the endometrium to pregnancy between HF and SF heifers. Considerable deficiencies in pregnancy-dependent biological pathways associated with extracellular matrix structure and organization as well as cell adhesion were found in the endometrium of SF animals. Distinct gene expression differences were also observed in conceptuses from HF and SF animals, with many of the genes decreased in SF conceptuses known to be embryonic lethal in mice due to defects in embryo and/or placental development. Analyses of biological pathways, key players, and ligand-receptor interactions based on transcriptome data divulged substantial evidence for dysregulation of conceptus-endometrial interactions in SF animals. These results support the ideas that the uterus impacts conceptus survival and programs conceptus development, and ripple effects of dysregulated conceptus-endometrial interactions elicit loss of the postelongation conceptus in SF cattle during the implantation period of pregnancy., Competing Interests: Conflict of interest statement: Editor R.M.R. has not published any research papers with any of the authors within the last four years and does not have collaborative research grants with them.

Bovine respiratory disease (BRD) is the most common infectious disease of beef and dairy cattle and is characterized by a complex infectious etiology that includes a variety of viral and bacterial pathogens. We examined the global changes in mRNA abundance in healthy lung and lung lesions and in the lymphoid tissues bronchial lymph node, retropharyngeal lymph node, nasopharyngeal lymph node and pharyngeal tonsil collected at the peak of clinical disease from beef cattle experimentally challenged with either bovine respiratory syncytial virus, infectious bovine rhinotracheitis, bovine viral diarrhea virus, Mannheimia haemolytica or Mycoplasma bovis. We identified signatures of tissue-specific transcriptional responses indicative of tropism in the coordination of host's immune tissue responses to infection by viral or bacterial infections. Furthermore, our study shows that this tissue tropism in host transcriptional response to BRD pathogens results in the activation of different networks of response genes. The differential crosstalk among genes expressed in lymphoid tissues was predicted to be orchestrated by specific immune genes that act as 'key players' within expression networks. The results of this study serve as a basis for the development of innovative therapeutic strategies and for the selection of cattle with enhanced resistance to BRD.

Infertility and subfertility negatively impact the economics and reproductive performance of cattle. Of note, significant pregnancy loss occurs in cattle during the first month of pregnancy, yet little is known about the genetic loci influencing pregnancy success and loss in cattle. To identify quantitative trait loci (QTL) with large effects associated with early pregnancy loss, Angus crossbred heifers were classified based on day 28 pregnancy outcomes to serial embryo transfer. A genome wide association analysis (GWAA) was conducted comparing 30 high fertility heifers with 100% success in establishing pregnancy to 55 subfertile heifers with 25% or less success. A gene set enrichment analysis SNP (GSEA-SNP) was performed to identify gene sets and leading edge genes influencing pregnancy loss. The GWAA identified 22 QTL (p < 1 x 10-5), and GSEA-SNP identified 9 gene sets (normalized enrichment score > 3.0) with 253 leading edge genes. Network analysis identified TNF (tumor necrosis factor), estrogen, and TP53 (tumor protein 53) as the top of 671 upstream regulators (p < 0.001), whereas the SOX2 (SRY [sex determining region Y]-box 2) and OCT4 (octamer-binding transcription factor 4) complex was the top master regulator out of 773 master regulators associated with fertility (p < 0.001). Identification of QTL and genes in pathways that improve early pregnancy success provides critical information for genomic selection to increase fertility in cattle. The identified genes and regulators also provide insight into the complex biological mechanisms underlying pregnancy establishment in cattle.

Background: Single nucleotide polymorphism (SNP) arrays for domestic cattle have catalyzed the identification of genetic markers associated with complex traits for inclusion in modern breeding and selection programs. Using actual and imputed Illumina 778K genotypes for 3887 U.S. beef cattle from 3 populations (Angus, Hereford, SimAngus), we performed genome-wide association analyses for feed efficiency and growth traits including average daily gain (ADG), dry matter intake (DMI), mid-test metabolic weight (MMWT), and residual feed intake (RFI), with marker-based heritability estimates produced for all traits and populations., Results: Moderate and/or large-effect QTL were detected for all traits in all populations, as jointly defined by the estimated proportion of variance explained (PVE) by marker effects (PVE ≥ 1.0%) and a nominal P-value threshold (P ≤ 5e-05). Lead SNPs with PVE ≥ 2.0% were considered putative evidence of large-effect QTL (n = 52), whereas those with PVE ≥ 1.0% but < 2.0% were considered putative evidence for moderate-effect QTL (n = 35). Identical or proximal lead SNPs associated with ADG, DMI, MMWT, and RFI collectively supported the potential for either pleiotropic QTL, or independent but proximal causal mutations for multiple traits within and between the analyzed populations. Marker-based heritability estimates for all investigated traits ranged from 0.18 to 0.60 using 778K genotypes, or from 0.17 to 0.57 using 50K genotypes (reduced from Illumina 778K HD to Illumina Bovine SNP50). An investigation to determine if QTL detected by 778K analysis could also be detected using 50K genotypes produced variable results, suggesting that 50K analyses were generally insufficient for QTL detection in these populations, and that relevant breeding or selection programs should be based on higher density analyses (imputed or directly ascertained)., Conclusions: Fourteen moderate to large-effect QTL regions which ranged from being physically proximal (lead SNPs ≤ 3Mb) to fully overlapping for RFI, DMI, ADG, and MMWT were detected within and between populations, and included evidence for pleiotropy, proximal but independent causal mutations, and multi-breed QTL. Bovine positional candidate genes for these traits were functionally conserved across vertebrate species.

Bovine respiratory disease complex (BRDC) is an important cause of mortality and morbidity in cattle; costing the dairy and beef industries millions of dollars annually, despite the use of vaccines and antibiotics. BRDC is caused by one or more of several viruses (bovine respiratory syncytial virus, bovine herpes type 1 also known as infectious bovine rhinotracheitis, and bovine viral diarrhea virus), which predispose animals to infection with one or more bacteria. These include: Pasteurella multocida, Mannheimia haemolytica, Mycoplasma bovis, and Histophilus somni. Some cattle appear to be more resistant to BRDC than others. We hypothesize that appropriate immune responses to these pathogens are subject to genetic control. To determine which genes are involved in the immune response to each of these pathogens it was first necessary to experimentally induce infection separately with each pathogen to document clinical and pathological responses in animals from which tissues were harvested for subsequent RNA sequencing. Herein these infections and animal responses are described.

Susceptibility to bovine respiratory disease (BRD) is multi-factorial and is influenced by stress in conjunction with infection by both bacterial and viral pathogens. While vaccination is broadly used in an effort to prevent BRD, it is far from being fully protective and cases diagnosed from a combination of observed clinical signs without any attempt at identifying the causal pathogens are usually treated with antibiotics. Dairy and beef cattle losses from BRD are profound worldwide and genetic studies have now been initiated to elucidate host loci which underlie susceptibility with the objective of enabling molecular breeding to reduce disease prevalence. In this study, we employed RNA sequencing to examine the bronchial lymph node transcriptomes of controls and beef cattle which had individually been experimentally challenged with bovine respiratory syncytial virus, infectious bovine rhinotracheitis, bovine viral diarrhea virus, Pasteurella multocida, Mannheimia haemolytica or Mycoplasma bovis to identify the genes that are involved in the bovine immune response to infection. We found that 142 differentially expressed genes were located in previously described quantitative trait locus regions associated with risk of BRD. Mutations affecting the expression or amino acid composition of these genes may affect disease susceptibility and could be incorporated into molecular breeding programs. Genes involved in innate immunity were generally found to be differentially expressed between the control and pathogen-challenged animals suggesting that variation in these genes may lead to a heritability of susceptibility that is pathogen independent. However, we also found pathogen-specific expression profiles which suggest that host genetic variation for BRD susceptibility is pathogen dependent.

Unlabelled: Bovine respiratory disease (BRD) is a common health problem for both dairy and beef cattle, resulting in significant economic loses. In order to identify viruses associated with BRD, we used a metagenomics approach to enrich and sequence viral nucleic acids in the nasal swabs of 50 young dairy cattle with symptoms of BRD. Following deep sequencing, de novo assembly, and translated protein sequence similarity searches, numerous known and previously uncharacterized viruses were identified. Bovine adenovirus 3, bovine adeno-associated virus, bovine influenza D virus, bovine parvovirus 2, bovine herpesvirus 6, bovine rhinitis A virus, and multiple genotypes of bovine rhinitis B virus were identified. The genomes of a previously uncharacterized astrovirus and picobirnaviruses were also partially or fully sequenced. Using real-time PCR, the rates of detection of the eight viruses that generated the most reads were compared for the nasal secretions of 50 animals with BRD versus 50 location-matched healthy control animals. Viruses were detected in 68% of BRD-affected animals versus 16% of healthy control animals. Thirty-eight percent of sick animals versus 8% of controls were infected with multiple respiratory viruses. Significantly associated with BRD were bovine adenovirus 3 (P < 0.0001), bovine rhinitis A virus (P = 0.005), and the recently described bovine influenza D virus (P = 0.006), which were detected either alone or in combination in 62% of animals with BRD. A metagenomics and real-time PCR detection approach in carefully matched cases and controls can provide a rapid means to identify viruses associated with a complex disease, paving the way for further confirmatory tests and ultimately to effective intervention strategies., Importance: Bovine respiratory disease is the most economically important disease affecting the cattle industry, whose complex root causes include environmental, genetics, and infectious factors. Using an unbiased metagenomics approach, we characterized the viruses in respiratory secretions from BRD cases and identified known and previously uncharacterized viruses belonging to seven viral families. Using a case-control format with location-matched animals, we compared the rates of viral detection and identified 3 viruses associated with severe BRD signs. Combining a metagenomics and case-control format can provide candidate pathogens associated with complex infectious diseases and inform further studies aimed at reducing their impact., (Copyright © 2015, American Society for Microbiology. All Rights Reserved.)

Background: Bovine respiratory disease complex (BRDC) is an infectious disease of cattle that is caused by a combination of viral and/or bacterial pathogens. Selection for cattle with reduced susceptibility to respiratory disease would provide a permanent tool for reducing the prevalence of BRDC. The objective of this study was to identify BRDC susceptibility loci in pre-weaned Holstein calves as a prerequisite to using genetic improvement as a tool for decreasing the prevalence of BRDC. High density SNP genotyping with the Illumina BovineHD BeadChip was conducted on 1257 male and 757 female Holstein calves from California (CA), and 767 calves identified as female from New Mexico (NM). Of these, 1382 were classified as BRDC cases, and 1396 were classified as controls, with all phenotypes assigned using the McGuirk health scoring system. During the acquisition of blood for DNA isolation, two deep pharyngeal and one mid-nasal diagnostic swab were obtained from each calf for the identification of bacterial and viral pathogens. Genome-wide association analyses were conducted using four analytical approaches (EIGENSTRAT, EMMAX-GRM, GBLUP and FvR). The most strongly associated SNPs from each individual analysis were ranked and evaluated for concordance. The heritability of susceptibility to BRDC in pre-weaned Holstein calves was estimated., Results: The four statistical approaches produced highly concordant results for 373 top ranked SNPs that defined 126 chromosomal regions for the CA population. Similarly, in NM, 370 SNPs defined 138 genomic regions that were identified by all four approaches. When the two populations were combined (i.e., CA + NM) and analyzed, 324 SNPs defined 116 genomic regions that were associated with BRDC across all analytical methods. Heritability estimates for BRDC were 21% for both CA and NM as individual populations, but declined to 13% when the populations were combined., Conclusions: Four analytical approaches utilizing both single and multi-marker association methods revealed common genomic regions associated with BRDC susceptibility that can be further characterized and used for genomic selection. Moderate heritability estimates were observed for BRDC susceptibility in pre-weaned Holstein calves, thereby supporting the application of genomic selection to reduce the prevalence of BRDC in U.S. Holsteins.

Background: The identification of genetic markers associated with complex traits that are expensive to record such as feed intake or feed efficiency would allow these traits to be included in selection programs. To identify large-effect QTL, we performed a series of genome-wide association studies and functional analyses using 50 K and 770 K SNP genotypes scored in 5,133 animals from 4 independent beef cattle populations (Cycle VII, Angus, Hereford and Simmental×Angus) with phenotypes for average daily gain, dry matter intake, metabolic mid-test body weight and residual feed intake., Results: A total of 5, 6, 11 and 10 significant QTL (defined as 1-Mb genome windows with Bonferroni-corrected P-value<0.05) were identified for average daily gain, dry matter intake, metabolic mid-test body weight and residual feed intake, respectively. The identified QTL were population-specific and had little overlap across the 4 populations. The pleiotropic or closely linked QTL on BTA 7 at 23 Mb identified in the Angus population harbours a promising candidate gene ACSL6 (acyl-CoA synthetase long-chain family member 6), and was the largest effect QTL associated with dry matter intake and mid-test body weight explaining 10.39% and 14.25% of the additive genetic variance, respectively. Pleiotropic or closely linked QTL associated with average daily gain and mid-test body weight were detected on BTA 6 at 38 Mb and BTA 7 at 93 Mb confirming previous reports. No QTL for residual feed intake explained more than 2.5% of the additive genetic variance in any population. Marker-based estimates of heritability ranged from 0.21 to 0.49 for residual feed intake across the 4 populations., Conclusions: This GWAS study, which is the largest performed for feed efficiency and its component traits in beef cattle to date, identified several large-effect QTL that cumulatively explained a significant percentage of additive genetic variance within each population. Differences in the QTL identified among the different populations may be due to differences in power to detect QTL, environmental variation, or differences in the genetic architecture of trait variation among breeds. These results enhance our understanding of the biology of growth, feed intake and utilisation in beef cattle.

Infertility and subfertility are important and pervasive reproductive problems in both domestic animals and humans. The majority of embryonic loss occurs during the first three weeks of pregnancy in cattle and women due, in part, to inadequate endometrial receptivity for support of embryo implantation. To identify heifers of contrasting fertility, serial rounds of artificial insemination (AI) were conducted in 201 synchronized crossbred beef heifers. The heifers were then fertility classified based on number of pregnancies detected on day 35 in four AI opportunities. Heifers, classified as having high fertility, subfertility or infertility, were selected for further study. The fertility-classified heifers were superovulated and flushed, and the recovered embryos were graded and then transferred to synchronized recipients. Quantity of embryos recovered per flush, embryo quality, and subsequent recipient pregnancy rates did not differ by fertility classification. Two in vivo-produced bovine embryos (stage 4 or 5, grade 1 or 2) were then transferred into each heifer on day 7 post-estrus. Pregnancy rates were greater in high fertility than lower fertility heifers when heifers were used as embryo recipients. The reproductive tracts of the classified heifers were obtained on day 14 of the estrous cycle. No obvious morphological differences in reproductive tract structures and histology of the uterus were observed in the heifers. Microarray analysis revealed differences in the endometrial transcriptome based on fertility classification. A genome-wide association study, based on SNP genotyping, detected 7 moderate associations with fertility across 6 different chromosomes. Collectively, these studies support the idea that innate differences in uterine function underlie fertility and early pregnancy loss in ruminants. Cattle with defined early pregnancy success or loss is useful to elucidate the complex biological and genetic mechanisms governing endometrial receptivity and uterine competency for pregnancy.

Background: Bovine paratuberculosis (ParaTB) also known as Johne's disease, is a contagious fatal disease resulting from infection by Mycobacterium avium subspecies paratuberculosis (MAP). Previous studies have identified loci associated with ParaTB using different measurements to define cases and controls. The objective of this study was to combine the data from two recent studies to identify genetic loci associated with MAP tissue infection and humoral immune response, defined by MAP ELISA-positive cattle, by comparing cases and control animals for one or both measures of infection., Methodology/principal Findings: The two populations used for the association analyses were a cohort of MAP tissue infected animals and control Holstein cows from the USA and the second cohort composed of ELISA-positive and ELISA-negative Holstein cows from Italy. Altogether 1190 cattle were genotyped with the Illumina BovineSNP50 BeadChip. SNP markers were removed if the minor allele frequency <0.01 or genotyping failure was >5%. Animals were removed with >5% genotyping failure. Whole genome association analyses were conducted with the GRAMMAR-CG method using two different definitions of control populations., Conclusion/significance: The analyses identified several loci (P<5 e-05) associated with ParaTB, defined by positive ELISA and presence of bacteria in tissue compared to ELISA and tissue negative animals, on chromosomes 1, 12 and 15 and one unassigned SNP. These results confirmed associations on chromosome 12 and the unassigned SNP with ParaTB which had been found in the Italian population alone. Furthermore, several additional genomic regions were found associated with ParaTB when ELISA and tissue positive animals were compared with tissue negative samples. These loci were on chromosomes 1, 6, 7, 13, 16, 21,23 and 25 (P<5 e-05). The results clearly indicate the importance of the phenotype definition when seeking to identify markers associated with different disease responses.

Background: Like human immunodeficiency virus (HIV), ovine lentivirus (OvLV) is macrophage-tropic and causes lifelong infection. OvLV infects one quarter of U.S. sheep and induces pneumonia and body condition wasting. There is no vaccine to prevent OvLV infection and no cost-effective treatment for infected animals. However, breed differences in prevalence and proviral concentration have indicated a genetic basis for susceptibility to OvLV. A recent study identified TMEM154 variants in OvLV susceptibility. The objective here was to identify additional loci associated with odds and/or control of OvLV infection., Methodology/principal Findings: This genome-wide association study (GWAS) included 964 sheep from Rambouillet, Polypay, and Columbia breeds with serological status and proviral concentration phenotypes. Analytic models accounted for breed and age, as well as genotype. This approach identified TMEM154 (nominal P=9.2×10(-7); empirical P=0.13), provided 12 additional genomic regions associated with odds of infection, and provided 13 regions associated with control of infection (all nominal P<1 × 10(-5)). Rapid decline of linkage disequilibrium with distance suggested many regions included few genes each. Genes in regions associated with odds of infection included DPPA2/DPPA4 (empirical P=0.006), and SYTL3 (P=0.051). Genes in regions associated with control of infection included a zinc finger cluster (ZNF192, ZSCAN16, ZNF389, and ZNF165; P=0.001), C19orf42/TMEM38A (P=0.047), and DLGAP1 (P=0.092)., Conclusions/significance: These associations provide targets for mutation discovery in sheep susceptibility to OvLV. Aside from TMEM154, these genes have not been associated previously with lentiviral infection in any species, to our knowledge. Further, data from other species suggest functional hypotheses for future testing of these genes in OvLV and other lentiviral infections. Specifically, SYTL3 binds and may regulate RAB27A, which is required for enveloped virus assembly of human cytomegalovirus. Zinc finger transcription factors have been associated with positive selection for repression of retroviral replication. DLGAP1 binds and may regulate DLG1, a known regulator of HIV infectivity.

Bovine respiratory disease (BRD) is considered to be the most costly infectious disease in the cattle industry. Bovine viral diarrhea virus (BVDV) is one of the pathogens involved with the BRD complex of disease. BVDV infection also negatively impacts cow reproduction and calf performance. Loci associated with persistently infected animals (BVD-PI) and linked with BRD have previously been identified near 14 Mb on bovine chromosome 2 (BTA2) and 15.3 Mb on bovine chromosome 26 (BTA26). The objective of this study was to refine the loci associated with BVD-PI and linked with BRD. Association testing for BVD-PI was performed on a population of 65 BVD-PI calves, 51 of their dams, and 60 unaffected calves (controls) with 142 single nucleotide polymorphisms (SNPs) on BTA2 and 173 SNPs on BTA26. Comparisons were made between BVD-PI calves and controls calves and the dams of BVD-PI calves and controls calves. For the linkage analysis of BRD, the same markers were used to genotype two half-sib families consisting of the sires and 72 BRD positive and 148 BRD negative offspring. Using an allelic chi-square test, 11 loci on BTA2 and 8 loci on BTA26 were associated with the dams of the BVD-PI calves (P < 0.05) and 4 loci on BTA2 and 11 loci on BTA26 were associated with BVD-PI calves. This demonstrates that although some of the loci on BTA2 and BTA26 are jointly involved in the fetal and dam response to BVD-PI infection, there are loci that are solely associated with the maternal or fetal susceptibility to disease. One locus on BTA2 and two loci on BTA26 were found to be linked (P < 0.05) with BRD. The regions linked with BRD were also associated with BVD-PI demonstrating that both the broad (BRD) and narrow (BVD-PI) definition of disease identified shared genomic regions as important in disease susceptibility. These results further refined the loci associated with BVD-PI and linked with BRD.

SNP-based gene-set enrichment analysis from single nucleotide polymorphisms, or GSEA-SNP, is a tool to identify candidate genes based on enrichment analysis of sets of genes rather than single SNP associations. The objective of this study was to identify modest-effect genes associated with Mycobacterium avium subsp. paratuberculosis (Map) tissue infection or fecal shedding using GSEA-SNP applied to KEGG pathways or Gene Ontology (GO) gene sets. The Illumina Bovine SNP50 BeadChip was used to genotype 209 Holstein cows for the GSEA-SNP analyses. For each of 13,744 annotated genes genome-wide located within 50 kb of a Bovine SNP50 SNP, the single SNP with the highest Cochran-Armitage Max statistic was used as a proxy statistic for that gene's strength of affiliation with Map. Gene-set enrichment was tested using a weighted Kolmogorov-Smirnov-like running sum statistic with data permutation to adjust for multiple testing. For tissue infection and fecal shedding, no gene sets in KEGG pathways or in GO sets for molecular function or cellular component were enriched for signal. The GO biological process gene set for positive regulation of cell motion (GO:0051272, q = 0.039, 5/11 genes contributing to the core enrichment) was enriched for Map tissue infection, while no GO biological process gene sets were enriched for fecal shedding. GSEA-SNP complements traditional SNP association approaches to identify genes of modest effects as well as genes with larger effects as demonstrated by the identification of one locus that we previously found to be associated with Map tissue infection using a SNP-by-SNP genome-wide association study.

The Pecorans (higher ruminants) are believed to have rapidly speciated in the Mid-Eocene, resulting in five distinct extant families: Antilocapridae, Giraffidae, Moschidae, Cervidae, and Bovidae. Due to the rapid radiation, the Pecoran phylogeny has proven difficult to resolve, and 11 of the 15 possible rooted phylogenies describing ancestral relationships among the Antilocapridae, Giraffidae, Cervidae, and Bovidae have each been argued as representations of the true phylogeny. Here we demonstrate that a genome-wide single nucleotide polymorphism (SNP) genotyping platform designed for one species can be used to genotype ancient DNA from an extinct species and DNA from species diverged up to 29 million years ago and that the produced genotypes can be used to resolve the phylogeny for this rapidly radiated infraorder. We used a high-throughput assay with 54,693 SNP loci developed for Bos taurus taurus to rapidly genotype 678 individuals representing 61 Pecoran species. We produced a highly resolved phylogeny for this diverse group based upon 40,843 genome-wide SNP, which is five times as many informative characters as have previously been analyzed. We also establish a method to amplify and screen genomic information from extinct species, and place Bison priscus within the Bovidae. The quality of genotype calls and the placement of samples within a well-supported phylogeny may provide an important test for validating the fidelity and integrity of ancient samples. Finally, we constructed a phylogenomic network to accurately describe the relationships between 48 cattle breeds and facilitate inferences concerning the history of domestication and breed formation.

Background: In 2006, an atypical U.S. case of bovine spongiform encephalopathy (BSE) was discovered in Alabama and later reported to be polymorphic for glutamate (E) and lysine (K) codons at position 211 in the bovine prion protein gene (Prnp) coding sequence. A bovine E211K mutation is important because it is analogous to the most common pathogenic mutation in humans (E200K) which causes hereditary Creutzfeldt - Jakob disease, an autosomal dominant form of prion disease. The present report describes a high-throughput matrix-associated laser desorption/ionization-time-of-flight mass spectrometry assay for scoring the Prnp E211K variant and its use to determine an upper limit for the K211 allele frequency in U.S. cattle., Results: The K211 allele was not detected in 6062 cattle, including those from five commercial beef processing plants (3892 carcasses) and 2170 registered cattle from 42 breeds. Multiple nearby polymorphisms in Prnp coding sequence of 1456 diverse purebred cattle (42 breeds) did not interfere with scoring E211 or K211 alleles. Based on these results, the upper bounds for prevalence of the E211K variant was estimated to be extremely low, less than 1 in 2000 cattle (Bayesian analysis based on 95% quantile of the posterior distribution with a uniform prior)., Conclusion: No groups or breeds of U.S. cattle are presently known to harbor the Prnp K211 allele. Because a carrier was not detected, the number of additional atypical BSE cases with K211 will also be vanishingly low.