Loci associated with conception rate in crossbred beef heifers

The inability of beef cattle to maintain full term pregnancies has become an economic concern for the beef industry. Herd management and nutritional improvements have alleviated environmental impacts on embryonic and fetal loss, yet additional gains can be made through genomic selection. The objectives of this study were to identify loci and gene-sets in crossbred beef heifers associated with the number of services required to become pregnant (TBRD) and heifer conception rate at first service (HCR1). Heifers (n = 709) from a commercial beef operation underwent one round of artificial insemination, before exposure to bulls for natural service for 50 days. Pregnancy and time of conception was determined by ultrasound 35 days after the breeding season. Heifers were genotyped using the GeneSeek (Lincoln, NE) Bovine GGP50K BeadChip prior to genome-wide association analyses (GWAA) conducted using an EIGENSTRAT-like model to identify loci associated (P < 1 × 10⁻⁵) with TBRD and HCR1. One locus was associated (P = 8.97 × 10⁻⁶) with TBRD on BTA19 and included the positional candidate gene ASIC2, which is differentially expressed in the endometrium of fertility classified heifers, and the positional candidate gene, SPACA3. Gene-set enrichment analyses using SNP (GSEA-SNP) data, was performed and identified one gene-set, oxidoreductase activity, acting on paired donors, with incorporation or reduction of molecular oxygen as enriched (NES = 3.15) with TBRD and contained nine leading edge genes that contributed to the enrichment of the gene set. The enriched gene-set is involved in catalyzing oxidation-reduction reactions, which have been associated with oxidative stressors impacting pregnancy success. No loci were associated nor gene-sets enriched with HCR1. Identification of loci, positional candidate genes, gene-sets and leading edge genes enriched for fertility facilitate genomic selection that allows producers to select for reproductively superior cattle, reduce costs associated with infertility, and increase percent calf crop.

Validating loci associated with bovine respiratory disease complex in pre-weaned Holstein calves

Bovine respiratory disease (BRD) is considered one of the most economically important diseases in the cattle industry. Ultimately, the selection of cattle that are less susceptible to disease will allow producers to reduce the prevalence of BRD and lessen its economic impact. The objective of this study was to validate previously identified loci associated with susceptibility to BRD in an independent population of 140 pre-weaned Holstein calves from Wisconsin (WI). Using the McGuirk health scoring system, calves were classified as either clinically affected with BRD (n = 35) or healthy (n = 105). Additive genotypic tests were performed for genomic regions previously associated with susceptibility to BRD in calves from California (CA) and New Mexico (NM). Using this method, 4 loci (P < 0.01) consisting of 10 SNP were validated in the WI population, including 2 loci from CA, 1 locus from NM, and 1 locus from a combined CA + NM population. Most of the positional candidate genes and transcription factor binding site motifs associated with these loci have functions related to innate and adaptive immune responses. The validation of loci associated with susceptibility to BRD in independent populations allows producers to more reliably select cattle that are less susceptible to BRD, improving animal welfare, decreasing the annual revenue losses, and lowering the prevalence of the disease.

Genome-wide association analysis and gene set enrichment analysis with SNP data identify genes associated with 305-day milk yield in Holstein dairy cows

Milk production traits, such as 305-day milk yield (305MY), have been under direct selection to improve production in dairy cows. Over the past 50 years, the average milk yield has nearly doubled, and over 56% of the increase is attributable to genetic improvement. As such, additional improvements in milk yield are still possible as new loci are identified. The objectives of this study were to detect SNPs and gene sets associated with 305MY in order to identify new candidate genes contributing to variation in milk production. A population of 781 primiparous Holstein cows from six central Washington dairies with records of 305MY and energy corrected milk were used to perform a genome-wide association analysis (GWAA) using the Illumina BovineHD BeadChip (777 962 SNPs) to identify QTL associated with 305MY (P < 1.0 × 10^-5 ). A gene set enrichment analysis with SNP data (GSEA-SNP) was performed to identify gene sets (normalized enrichment score > 3.0) and leading edge genes (LEGs) influencing 305MY. The GWAA identified three QTL comprising 34 SNPs and 30 positional candidate genes. In the GSEA-SNP, five gene sets with 58 unique and 24 shared LEGs contributed to 305MY. Identification of QTL and LEGs associated with 305MY can provide additional targets for genomic selection to continue to improve 305MY in dairy cattle.

Gene set enrichment analysis of SNP data in dairy and beef cattle with bovine respiratory disease

Bovine respiratory disease (BRD) is a complex disease that is associated with infection by bacterial and viral pathogens when cattle fail to adequately respond to stress. The objective of this study was to use gene set enrichment analysis of SNP data (GSEA-SNP) and a network analysis (ingenuity pathway analysis) to identify gene sets, genes within gene sets (leading-edge genes) and upstream regulators associated with BRD in pre-weaned dairy calves and beef feedlot cattle. BRD cases and controls were diagnosed using the McGuirk health scoring system. Holstein calves were sampled from commercial calf-raising facilities in California (1003 cases and 1011 controls) and New Mexico (376 cases and 372 controls). Commercial feedlot cattle were sampled from Colorado (500 cases and 499 controls) and Washington (504 cases and 497 controls). There were 102 and 237 unique leading-edge genes identified in the dairy calf and beef cattle populations respectively. Six leading-edge genes (ADIPOQ, HTR2A, MIF, PDE6G, PRDX3 and SNCA) were associated with BRD in both dairy and beef cattle. Network analysis identified glucose as the most influential upstream regulator in dairy cattle, whereas in beef cattle, TNF was the most influential upstream regulator. The genes, gene sets and upstream regulators associated with BRD have common functions associated with immunity, inflammation and pulmonary disease and provide insights into the mechanisms that are critical to BRD susceptibility in cattle.

Rapid Communication: Subclinical bovine respiratory disease - loci and pathogens associated with lung lesions in feedlot cattle

Bovine respiratory disease (BRD) is an economically important disease of feedlot cattle that is caused by viral and bacterial pathogen members of the BRD complex. Many cases of subclinical BRD go untreated and are not detected until slaughter, when lung lesions are identified. The objectives of this study were to identify which BRD pathogens were associated with the presence of lung lesions at harvest and to identify genomic loci that were associated with susceptibility to lung lesions as defined by consolidation of the lung and/or the presence of fibrin tissue. Steers from a Colorado feedlot (n = 920) were tested for the presence of viral and bacterial pathogens using deep pharyngeal and mid-nasal swabs collected on entry into the study. Pathogen profiles were compared between cattle with or without lung consolidation (LC), fibrin tissue in the lung (FT), a combination of LC and FT in the same lung (lung lesions [LL]), and hyperinflated lungs (HIF) at harvest. Genotyping was conducted using the Illumina BovineHD BeadChip. Genomewide association analyses (GWAA) were conducted using EMMAX (efficient mixed-model association eXpedited), and pseudoheritabilities were estimated. The pathogen profile comparisons revealed that LC (P = 0.01, odds ratio [OR] = 3.37) and LL cattle (P = 0.04, OR = 4.58) were more likely to be infected with bovine herpes virus-1 and that HIF cattle were more likely to be infected with Mycoplasma spp. (P = 0.04, OR = 4.33). Pseudoheritability estimates were 0.25 for LC, 0.00 for FT, 0.28 for LL, and 0.13 for HIF. Because pseudoheritability for FT was estimated to be 0, GWAA results for FT were not reported. There were 4 QTL that were moderately associated (P < 1 × 10⁻⁵) with only LC, 2 that were associated with only LL, and 1 that was associated with LC and LL. Loci associated with HIF included 12 that were moderately associated and 3 that were strongly associated (uncorrected P < 5 × 10−7). A 24-kb region surrounding significant lead SNP was investigated to identify positional candidate genes. Many positional candidate genes underlying or flanking the detected QTL have been associated with signal transduction, cell adhesion, or gap junctions, which have functional relevance to the maintenance of lung health. The identification of pathogens and QTL associated with the presence of lung abnormalities in cattle exhibiting subclinical BRD allows the identification of loci that may not be detected through manifestation of clinical disease alone.

Identification of genes associated with susceptibility to Mycobacterium avium ssp. paratuberculosis (Map) tissue infection in Holstein cattle using gene set enrichment analysis-SNP

Multiple genome-wide association analyses have investigated susceptibility to bovine paratuberculosis, but few loci have been identified across independent cattle populations. A SNP-based gene set enrichment analysis (GSEA-SNP) allows expanded identification of genes with moderate effects on a trait through the enrichment of gene sets instead of identifying only few loci with large effects. Therefore, the objective of this study was to identify genes that were moderately associated with Mycobacterium avium ssp. paratuberculosis (Map) tissue infection using GSEA-SNP in Holstein cattle from the Pacific Northwest (PNW; n = 205) and from the PNW and Northeast (PNW+NE; n = 245) which were previously genotyped with the Illumina BovineSNP50 BeadChip. The GSEA-SNP utilized 4389 gene sets from five databases. For each annotated gene in the UMD3.1 assembly (n = 19,723), the most significant SNP within each gene and its surrounding region (10 kb up- and downstream) was selected as a proxy for that gene. Any gene set with a normalized enrichment score > 2.5 was considered enriched. Thirteen gene sets (8 PNW GSEA-SNP; 5 PNW+NE) were enriched in these analyses and all have functions that relate to nuclear factor kappa beta. Nuclear factor kappa beta is critical to gut immune responses, implicated in host immune responses to other mycobacterial diseases, and has established roles in inflammation as well as cancer. Gene sets and genes moderately associated with Map infection could be used in genomic selection to allow producers to select for less susceptible cattle, lower the prevalence of the disease, and reduce economic losses.

Identification of loci associated with susceptibility to subspecies () tissue infection in cattle

Johne's disease is a contagious bacterial infection of cattle caused by ssp. (). A previous genome-wide association analysis (GWAA) in Holstein cattle identified QTL on BTA3 and BTA9 that were highly associated (P < 5 × 10) and on BTA1, BTA16, and BTA21 that were moderately associated (P < 5 × 10) with Map tissue infection. The objectives of this study were to validate previous GWAA results in Jersey cattle ( = 57), Holstein cattle from the Pacific Northwest (PNW, = 205) and a combined Holstein population from the PNW and the Northeast (PNW + NE, = 423), and also identify new loci associated with tissue infection. DNA was genotyped using the Illumina BovineSNP50 BeadChip, and the PNW + NE data was also imputed to whole genome sequence level using Run4 of the 1000 Bull Genomes project with Beagle v 4.1 and FImpute. Cases were ileocecal node positive and controls were negative for by quantitative PCR (qPCR). Individuals were removed for SNP call rate < 90%, and SNP were removed for genotype call rate < 90% or minor allele frequency < 1%. For the Jersey, PNW, and PNW + NE, GWAA were conducted using an allelic dosage model. For the PNW and the PNW + NE, an additional efficient mixed-model association eXpedited (EMMAX) analysis was performed using additive, dominance and recessive models. Seven QTL on BTA22 were identified in the Jersey population with the most significant ( = 4.45 × 10) located at 21.7 megabases (Mb). Six QTL were associated in the PNW and the PNW + NE analyses, including a QTL previously identified on BTA16 in the NE population. The most significant locus for the PNW was located on BTA21 at 61 Mb ( = 8.61 × 10) while the most significant locus for the PNW + NE was on BTA12 at 90 Mb ( = 2.33 × 10). No additional QTL were identified with the imputed GWAA. Putative positional candidate genes were identified within 50 kb 5' and 3' of each QTL. Two positional candidate genes were identified in Jersey cattle, 1 identified in the PNW and 8 in the PNW + NE populations. Many identified positional candidate genes are involved in signal transduction, have immunological functions, or have putative functional relevance in entry into host cells. This study supported 2 previously identified SNP within a QTL on BTA16 and identified 16 new QTL, including 2 found in the PNW and the PNW+NE, associated with tissue infection.