Identifying Regions of the Genome Associated with Conception Rate to the First Service in Holstein Heifers Bred by Artificial Insemination and as Embryo Transfer Recipients.

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.

Genomic regions associated with pseudorabies virus infection status in naturally infected feral swine (Sus scrofa).

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.

Identification by whole genome sequencing of genes associated with delayed postoperative hemorrhage in Scottish deerhounds.

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.

Genome-wide association with footrot in hair and wool sheep.

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.

A high-density genome-wide association with absolute blood monocyte count in domestic sheep identifies novel loci.

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.

Identifying Loci Associated With Bovine Corona Virus Infection and Bovine Respiratory Disease in Dairy and Feedlot Cattle.

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.

A DNA Regulatory Element Haplotype at Zinc Finger Genes Is Associated with Host Resilience to Small Ruminant Lentivirus in Two Sheep Populations.

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.

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.

Key genetic variants associated with variation of milk oligosaccharides from diverse human populations.

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.

Influenza C and D viral load in cattle correlates with bovine respiratory disease (BRD): Emerging role of orthomyxoviruses in the pathogenesis of BRD.

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.

The EDN2 rs110287192 gene polymorphism is associated with paratuberculosis susceptibility in multibreed cattle population.

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.

Identification of Loci and Pathways Associated with Heifer Conception Rate in U.S. Holsteins.

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.

Genomic Analysis of Spontaneous Abortion in Holstein Heifers and Primiparous Cows.

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.

Identification of loci associated with conception rate in primiparous Holstein cows.

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.

Validation of 46 loci associated with female fertility traits in cattle.

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.

QTL-mapping and genomic prediction for bovine respiratory disease in U.S. Holsteins using sequence imputation and feature selection.

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.

Rapid Communication: Genome-wide association analyses identify loci associated with colostrum production in Jersey cattle1.

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.

Uterine influences on conceptus development in fertility-classified animals.

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.

Loci and pathways associated with uterine capacity for pregnancy and fertility in beef cattle.

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.

Tissue Tropism in Host Transcriptional Response to Members of the Bovine Respiratory Disease Complex.

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.