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

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⁻⁸) with HCR1 and 615 unique loci associated (p < 5 × 10⁻⁸) 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⁻⁵) 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.

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

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.