Functional Variants Surrounding Endothelin 2 Are Associated With Mycobacterium avium Subspecies paratuberculosis Infection

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.

A Population Genomics Analysis of the Native Irish Galway Sheep Breed

The Galway sheep population is the only native Irish sheep breed and this livestock genetic resource is currently categorised as 'at-risk'. In the present study, comparative population genomics analyses of Galway sheep and other sheep populations of European origin were used to investigate the microevolution and recent genetic history of the breed. These analyses support the hypothesis that British Leicester sheep were used in the formation of the Galway. When compared to conventional and endangered breeds, the Galway breed was intermediate in effective population size, genomic inbreeding and runs of homozygosity. This indicates that, although the Galway breed is declining, it is still relatively genetically diverse and that conservation and management plans informed by genomic information may aid its recovery. The Galway breed also exhibited distinct genomic signatures of artificial or natural selection when compared to other breeds, which highlighted candidate genes that may be involved in production and health traits.

Susceptibility loci revealed for bovine respiratory disease complex in pre-weaned holstein calves

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.

Estimation of inbreeding and effective population size of full-blood Wagyu cattle registered with the American Wagyu Cattle Association

The objective of this research was to examine the population structure of full-blood (100%) Wagyu cattle registered in the United States with the American Wagyu Association, with the aim of estimating and comparing the levels of inbreeding from both pedigree and genotypic data. A total of 4132 full-blood Wagyu cattle pedigrees were assessed and used to compute the inbreeding coefficients (FIT and FST ) and the effective population size (Ne ) from pedigree data for the period 1994 to 2011. In addition to pedigree analysis, 47 full-blood Wagyu cattle representing eight prominent sire lines in the American Wagyu cattle population were genotyped using the Illumina BovineSNP50 BeadChip. Genotypic data were then used to estimate genomic inbreeding coefficients (FROH ) by calculating runs of homozygosity. The mean inbreeding coefficient based on the pedigree data was estimated at 4.80%. The effective population size averaged 17 between the years 1994 and 2011 with an increase of 42.9 in 2000 and a drop of 1.8 in 2011. Examination of the runs of homozygosity revealed that the 47 Wagyu cattle from the eight prominent sire lines had a mean genomic inbreeding coefficient (FROH ) estimated at 9.08% compared to a mean inbreeding coefficient based on pedigree data of 4.8%. These data suggest that the mean genotype inbreeding coefficient of full-blood Wagyu cattle exceeds the inbreeding coefficient identified by pedigree. Inbreeding has increased slowly at a rate of 0.03% per year over the past 17 years. Wagyu breeders should continue to utilize many sires from divergent lines and consider outcrossing to other breeds to enhance genetic diversity and minimize the adverse effects of inbreeding in Wagyu.

Technical note: High fidelity of whole-genome amplified sheep (Ovis aries) deoxyribonucleic acid using a high-density single nucleotide polymorphism array-based genotyping platform

Advances in high-throughput genotyping technologies have afforded researchers the opportunity to study ever-increasing numbers of SNP in animal genomes. However, many studies encounter difficulties in obtaining sufficient quantities of high-quality DNA for such analyses, particularly when the source biological material is limited or degraded. The recent development of in vitro whole-genome amplification approaches has permitted researchers to circumvent these challenges by increasing the amount of usable DNA in normally small-quantity samples. Here, we assess the performance of whole-genome amplification products generated from ovine genomic DNA using a high-throughput SNP genotyping platform, the newly developed Illumina ovineSNP50 BeadChip. Our results demonstrate a high genotype call rate for conventional genomic DNA and whole-genome amplified genomic DNA. The data also reveal an exceptionally high concordance rate ( > or = 99%) between the genotypes generated from whole-genome amplified products and their conventional genomic DNA counterparts. This study supports the use of whole-genome amplification as a viable solution for the analysis of high-density SNP genotypic data using compromised or limited starting material.

Linkage mapping of the locus for inherited ovine arthrogryposis (IOA) to sheep chromosome 5

Arthrogryposis is a congenital malformation affecting the limbs of newborn animals and infants. Previous work has demonstrated that inherited ovine arthrogryposis (IOA) has an autosomal recessive mode of inheritance. Two affected homozygous recessive (art/art) Suffolk rams were used as founders for a backcross pedigree of half-sib families segregating the IOA trait. A genome scan was performed using 187 microsatellite genetic markers and all backcross animals were phenotyped at birth for the presence and severity of arthrogryposis. Pairwise LOD scores of 1.86, 1.35, and 1.32 were detected for three microsatellites, BM741, JAZ, and RM006, that are located on sheep Chr 5 (OAR5). Additional markers in the region were identified from the genetic linkage map of BTA7 and by in silico analyses of the draft bovine genome sequence, three of which were informative. Interval mapping of all autosomes produced an F value of 21.97 (p < 0.01) for a causative locus in the region of OAR5 previously flagged by pairwise linkage analysis. Inspection of the orthologous region of HSA5 highlighted a previously fine-mapped locus for human arthrogryposis multiplex congenita neurogenic type (AMCN). A survey of the HSA5 genome sequence identified plausible candidate genes for both IOA and human AMCN.

48 PREGNANCY AND FETAL CHARACTERISTICS AFTER TRANSFER OF VITRIFIED IN VIVO AND CLONED BOVINE EMBRYOS

This study was conducted to examine pregnancy progression and fetal characteristics following transfer of vitrified bovine nuclear transfer (NT) vs. in vivo-derived (Vivo) embryos. Nuclear transfer was conducted using skin fibroblast cells derived from cultured ear explants taken from an elite Holstein-Friesian dairy cow. Expanding and hatching blastocysts on Day 7 were vitrified using liquid nitrogen surface vitrification (Shen et al. 2004 Reprod. Fertil. Dev. 16, 158). Day 7 in vivo embryos, produced using standard superovulation procedures applied to Holstein-Friesian heifers (n = 6), were vitrified in the same way. Following warming, embryos were transferred to synchronized recipients (2 per recipient; NT: n = 81 recipients; Vivo: n = 20 recipients). Pregnancies were monitored by ultrasound scanning on Days 25, 45, and 75, and a sample of animals was slaughtered at each time point to recover the fetus/placenta for further analyses. On Day 25, the entire conceptus was weighed together; on Day 45, the fetus and placenta were weighed separately; on Day 75, the fetus was dissected and the major organs were weighed. Significantly more animals remained pregnant after transfer of in vivo-derived embryos than NT embryos at all time points: Day 25 (95.0 vs. 61.7%, P < 0.01), Day 45 (93.3 vs. 39.7%, P < 0.001), and Day 75 (90.9 vs. 19.8%, P < 0.001). There was no significant difference (P = 0.10) in the weight of the conceptus on Day 25 from NT transfers (1.14 ± 0.23 g, n = 8) vs. Vivo transfers (0.75 ± 0.19 g, n = 8). On Day 45, there was no significant difference in the weight of either fetus (P = 0.393) or membranes (P = 0.167) between NT embryos (fetus: 2.76 ± 0.40 g, n = 12; membranes 59.0 ± 10.0 g, n = 11) or in vivo-derived embryos (fetus: 2.60 ± 0.15 g, n = 6; membranes 41.8 ± 5.2 g, n = 4). However, on Day 75, the weights of the fetus and several of the major organs were heavier from NT embryos (Table 1). These data suggest that the large calf syndrome is manifested after Day 45. Table 1. Fetal characteristics on Day 75 of gestation following transfer of vitrified in vivo-derived vs. NT blastocysts This work was supported by the U.S. Department of Agriculture.