Genomewide association study of lung lesions in cattle using sample pooling

Genome-Wide Association Study Reveals Quantitative Trait Loci and Candidate Genes Associated with High Interferon-gamma Production in Holstein Cattle Naturally Infected with Mycobacterium Bovis

Mycobacterium bovis (Mb) is the causative agent of bovine tuberculosis (bTb). Genetic selection aiming to identify less susceptible animals has been proposed as a complementary measure in ongoing programs toward controlling Mb infection. However, individual animal phenotypes for bTb based on interferon-gamma (IFNɣ) and its use in bovine selective breeding programs have not been explored. In the current study, IFNɣ production was measured using a specific IFNɣ ELISA kit in bovine purified protein derivative (bPPD)-stimulated blood samples collected from Holstein cattle. DNA isolated from the peripheral blood samples collected from the animals included in the study was genotyped with the EuroG Medium Density bead Chip, and the genotypes were imputed to whole-genome sequences. A genome-wide association analysis (GWAS) revealed that the IFNɣ in response to bPPD was associated with a specific genetic profile (heritability = 0.23) and allowed the identification of 163 SNPs, 72 quantitative trait loci (QTLs), 197 candidate genes, and 8 microRNAs (miRNAs) associated with this phenotype. No negative correlations between this phenotype and other phenotypes and traits included in the Spanish breeding program were observed. Taken together, our results define a heritable and distinct immunogenetic profile associated with strong production of IFNɣ in response to Mb.

Identification of Susceptibility Genes Underlying Bovine Respiratory Disease in Xinjiang Brown Cattle Based on DNA Methylation

DNA methylation is a form of epigenetic regulation, having pivotal parts in controlling cellular expansion and expression levels within genes. Although blood DNA methylation has been studied in humans and other species, its prominence in cattle is largely unknown. This study aimed to methodically probe the genomic methylation map of Xinjiang brown (XJB) cattle suffering from bovine respiratory disease (BRD), consequently widening cattle blood methylome ranges. Genome-wide DNA methylation profiling of the XJB blood was investigated through whole-genome bisulfite sequencing (WGBS). Many differentially methylated regions (DMRs) obtained by comparing the cases and controls groups were found within the CG, CHG, and CHH (where H is A, T, or C) sequences (16,765, 7502, and 2656, respectively), encompassing 4334 differentially methylated genes (DMGs). Furthermore, GO/KEGG analyses showed that some DMGs were involved within immune response pathways. Combining WGBS-Seq data and existing RNA-Seq data, we identified 71 significantly differentially methylated (DMGs) and expressed (DEGs) genes (p < 0.05). Next, complementary analyses identified nine DMGs (LTA, STAT3, IKBKG, IRAK1, NOD2, TLR2, TNFRSF1A, and IKBKB) that might be involved in the immune response of XJB cattle infected with respiratory diseases. Although further investigations are needed to confirm their exact implication in the involved immune processes, these genes could potentially be used for a marker-assisted selection of animals resistant to BRD. This study also provides new knowledge regarding epigenetic control for the bovine respiratory immune process.

Estimation of genetic parameters for visceral diseases of fattening Japanese Black cattle in Shimane Prefecture, Japan

In this study, the genetic parameters of major visceral diseases were estimated using the postmortem inspection records of 9057 fattening Japanese Black cattle in Shimane Prefecture, Japan, and the genetic correlation between visceral diseases and carcass traits was analyzed. There were six visceral diseases with a prevalence of 5% or higher, namely, pleurisy, pneumonia, bovine abdominal fat necrosis (BFN), rumenitis, hemorrhagic hepatitis, and perihepatitis. Variance components were estimated using the Gibbs sampling method, and the heritability of the visceral disease ranged from 0.07 to 0.49 for perihepatitis and BFN, respectively. Significant negative genetic correlations were identified between pleurisy and rib thickness (-0.32), BFN and carcass weight (-0.29), and BFN and rib eye area (-0.22). No significant genetic correlation was observed among the visceral diseases. The least squares analysis of variance suggested that some visceral diseases decrease the value of carcass traits. In particular, carcass weight and rib eye area in individuals with BFN were 11.7 kg and 1.87 cm2 lower than those of healthy cattle, respectively. Thus, it was inferred that genetic factors were involved in the visceral diseases of fattening Japanese Black cattle in Shimane Prefecture.

Use of overlapping DNA pools to discern genetic differences despite pooling error

Genotyping pools of commercial cattle and individual seedstock animals may reveal hidden relationships between sectors enabling use of commercial data for genetic evaluation. However, commercial data capture may be compromised by inexact pool formation. We aimed to estimate the concordance between distances or genomic covariance among pooling allele frequencies (PAFs) of DNA pools comprised of 100 animals with 0% or 50% overlap of animals in common between pools. Cattle lung samples were collected from a commercial beef processing plant on a single day. Six pools of 100 animals each were constructed so that overlap between pools was 0% or 50%. Two pools of all 200 animals were constructed to estimate PAFs for all 200 animals. Frozen lung tissue (0.01 g) from each animal was weighed into a tube containing a pool; there were two pools of 200 animals each and six pools of 100 animals each. Every contribution of an individual animal was an independent measurement to insure independence of pooling errors. Lung samples were kept on dried ice during the pooling process to keep them from thawing. The eight pools were then assayed for approximately 100,000 single nucleotide polymorphisms (SNP). PAF for each SNP and pool was based on the relative intensity of the two dyes used to detect the alleles rather than genotype calls which are not tractable from pooling data. Euclidean distances and genomic relationships among the PAFs for the eight pools were estimated and distances were tested for concordance with pool overlap using permutation-based analysis of distance. Distances among pools were concordant with the planned overlap of animals shared between pools (P = 0.0024); pool overlap accounted for 70% of the variation and pooling error accounted for 30%. Pools containing 100 animals with no overlap were the most distant from one another and pools with 50% overlap were the least distant. This work shows that we can discern differences in distance between pairs of overlapping DNA pools sharing 0% and 50% of the animals. Genomic correlations among nonoverlapping pools indicated that nonoverlapping pool pairs did not share many related animals because genomic correlations were near zero for these pairs. On the other hand, one pair of nonoverlapping pools likely contained related animals between pools because the correlation was 0.21. Pools sharing 50% overlap ranged in genomic relationship between 0.21 and 0.39 (N = 12).

Transcript Characteristics on the Susceptibility Difference of Bovine Respiratory Disease

Bovine respiratory disease (BRD) is one of the major health issues in the cattle industry, resulting in significant financial crises globally. There is currently no good treatment, and cattle are made resistant to pneumonia through disease-resistant breeding. The serial blood samples from six Xinjiang brown (XJB) calves were collected for the RNA sequencing (RNA-seq). The obtained six samples were grouped into two groups, in each group as infected with BRD and healthy calves, respectively. In our study, the differential expression mRNAs were detected by using RNA-seq and constructed a protein-protein interaction (PPI) network related to the immunity in cattle. The key genes were identified by protein interaction network analysis, and the results from RNA-seq were verified using reverse transcription-quantitative polymerase chain reaction (RT-qPCR). A total of 488 differentially expressed (DE) mRNAs were identified. Importantly, the enrichment analysis of these identified DEGs classified them as mainly enriched in the regulation and immune response processes. The 16 hub genes were found to be related to immune pathways categorized by PPIs analysis. Results revealed that many hub genes were related to the immune response to respiratory disease. These results will provide the basis for a better understanding of the molecular mechanism of bovine resistance to BRD.

Accuracy of GEBV of sires based on pooled allele frequency of their progeny

Despite decreasing genotyping costs, in some cases individually genotyping animals is not economically feasible (e.g., in small ruminants). An alternative is to pool DNA, using the pooled allele frequency (PAF) to garner information on performance. Still, the use of PAF for prediction (estimation of genomic breeding values; GEBVs) has been limited. Two potential sources of error on accuracy of GEBV of sires, obtained from PAF of their progeny themselves lacking pedigree information, were tested: (i) pool construction error (unequal contribution of DNA from animals in pools), and (ii) technical error (variability when reading the array). Pooling design (random, extremes, K-means), pool size (5, 10, 25, 50, and 100 individuals), and selection scenario (random, phenotypic) also were considered. These factors were tested by simulating a sheep population. Accuracy of GEBV-the correlation between true and estimated values-was not substantially affected by pool construction or technical error, or selection scenario. A significant interaction, however, between pool size and design was found. Still, regardless of design, mean accuracy was higher for pools of 10 or less individuals. Mean accuracy of GEBV was 0.174 (SE 0.001) for random pooling, and 0.704 (SE 0.004) and 0.696 (SE 0.004) for extreme and K-means pooling, respectively. Non-random pooling resulted in moderate accuracy of GEBV. Overall, pooled genotypes can be used in conjunction with individual genotypes of sires for moderately accurate predictions of their genetic merit with little effect of pool construction or technical error.

Evaluating Accuracy of DNA Pool Construction Based on White Blood Cell Counts

Pooling individual samples prior to DNA extraction can mitigate the cost of DNA extraction and genotyping; however, these methods need to accurately generate equal representation of individuals within pools. The objective of this study was to determine accuracy of pool construction of blood samples based on white blood cell counts compared to two common DNA quantification methods. Fifty individual bovine blood samples were collected, and then pooled with all individuals represented in each pool. Pools were constructed with the target of equal representation of each individual animal based on number of white blood cells, spectrophotometric readings, spectrofluorometric readings, and whole blood volume with 9 pools per method and a total of 36 pools. Pools and individual samples that comprised the pools were genotyped using a commercially available genotyping array. ASReml was used to estimate variance components for individual animal contribution to pools. The correlation between animal contributions between two pools was estimated using bivariate analysis with starting values set to the result of a univariate analysis. Adonis test on distance matrix from the animal correlation showed clustering with method, and higher correlations between methods than within (P < 1 × 10^-6). White blood cell count was predictive of sample representation when compared to pooling based on DNA concentration. Therefore, constructing pools using white blood cell counts prior to DNA extraction may reduce cost associated with DNA extraction and genotyping and improve representation of individuals in a pool.

ATAC-Seq identifies regions of open chromatin in the bronchial lymph nodes of dairy calves experimentally challenged with bovine respiratory syncytial virus

Background

Bovine Respiratory Syncytial Virus (BRSV) is a cause of Bovine Respiratory Disease (BRD). DNA-based biomarkers contributing to BRD resistance are potentially present in non-protein-coding regulatory regions of the genome, which can be determined using ATAC-Seq. The objectives of this study were to: (i) identify regions of open chromatin in DNA extracted from bronchial lymph nodes (BLN) of healthy dairy calves experimentally challenged with BRSV and compare them with those from non-challenged healthy control calves, (ii) elucidate the chromatin regions that were differentially or uniquely open in the BRSV challenged relative to control calves, and (iii) compare the genes found in regions proximal to the differentially open regions to the genes previously found to be differentially expressed in the BLN in response to BRSV and to previously identified BRD susceptibility loci. This was achieved by challenging clinically healthy Holstein-Friesian calves (mean age 143 ± 14 days) with either BRSV inoculum (n = 12) or with sterile phosphate buffered saline (PBS) (n = 6) and preparing and sequencing ATAC-Seq libraries from fresh BLN tissues.

Results

Using Diffbind, 9,144 and 5,096 differentially accessible regions (P < 0.05, FDR < 0.05) were identified between BRSV challenged and control calves employing DeSeq2 and EdgeR, respectively. Additionally, 8,791 chromatin regions were found to be uniquely open in BRSV challenged calves. Seventy-six and 150 of the genes that were previously found to be differentially expressed using RNA-Seq, were located within 2 kb downstream of the differentially accessible regions, and of the regions uniquely open in BRSV challenged calves, respectively. Pathway analyses within ClusterProfiler indicated that these genes were involved in immune responses to infection and participated in the Th1 and Th2 pathways, pathogen recognition and the anti-viral response. There were 237 differentially accessible regions positioned within 40 previously identified BRD susceptibility loci.

Conclusions

The identified open chromatin regions are likely to be involved in the regulatory response of gene transcription induced by infection with BRSV. Consequently, they may contain variants which impact resistance to BRD that could be used in breeding programmes to select healthier, more robust cattle.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-020-07268-5.

Genomics: the host's genotype and its relevance to bovine respiratory disease

Genomic variation exists in cattle that affects their susceptibility to the complex of pathogens responsible for bovine respiratory disease (BRD). Heritability estimates and genome-wide association analyses (GWAA) support the role of host genomic variation in BRD susceptibility. Heritability estimates for BRD susceptibility range from 0.02 to 0.29 depending on the population, the definition of the disease, and the accuracy of diagnosis. GWAA have identified genomic regions (loci) associated with BRD in beef and dairy cattle based on a variety of BRD diagnostic criteria. National standards need to be developed for BRD diagnostics and reporting to facilitate selection. Commercial genotyping is available to predict BRD susceptibility in dairy cattle and for the selection of replacement animals. Disease pathogen profiles vary by region and can result in genetic heterogeneity where different loci are important for susceptibility to different BRD pathogens. Although the identification of the BRD pathogens may not be critical for treatment, it is of paramount importance in identifying loci that render cattle susceptible to the disease. Identification of loci associated with host susceptibility to BRD provides a foundation for genomic selection to reduce disease and opens the possibilities to a better understanding of how the host defends itself.

Pooled genotyping strategies for the rapid construction of genomic reference populations1

The growing concern with the environment is making important for livestock producers to focus on selection for efficiency-related traits, which is a challenge for commercial cattle herds due to the lack of pedigree information. To explore a cost-effective opportunity for genomic evaluations of commercial herds, this study compared the accuracy of bulls' genomic estimated breeding values (GEBV) using different pooled genotype strategies. We used ten replicates of previously simulated genomic and phenotypic data for one low (t1) and one moderate (t2) heritability trait of 200 sires and 2,200 progeny. Sire's GEBV were calculated using a univariate mixed model, with a hybrid genomic relationship matrix (h-GRM) relating sires to: 1) 1,100 pools of 2 animals; 2) 440 pools of 5 animals; 3) 220 pools of 10 animals; 4) 110 pools of 20 animals; 5) 88 pools of 25 animals; 6) 44 pools of 50 animals; and 7) 22 pools of 100 animals. Pooling criteria were: at random, grouped sorting by t1, grouped sorting by t2, and grouped sorting by a combination of t1 and t2. The same criteria were used to select 110, 220, 440, and 1,100 individual genotypes for GEBV calculation to compare GEBV accuracy using the same number of individual genotypes and pools. Although the best accuracy was achieved for a given trait when pools were grouped based on that same trait (t1: 0.50-0.56, t2: 0.66-0.77), pooling by one trait impacted negatively on the accuracy of GEBV for the other trait (t1: 0.25-0.46, t2: 0.29-0.71). Therefore, the combined measure may be a feasible alternative to use the same pools to calculate GEBVs for both traits (t1: 0.45-0.57, t2: 0.62-0.76). Pools of 10 individuals were identified as representing a good compromise between loss of accuracy (~10%-15%) and cost savings (~90%) from genotype assays. In addition, we demonstrated that in more than 90% of the simulations, pools present higher sires' GEBV accuracy than individual genotypes when the number of genotype assays is limited (i.e., 110 or 220) and animals are assigned to pools based on phenotype. Pools assigned at random presented the poorest results (t1: 0.07-0.45, t2: 0.14-0.70). In conclusion, pooling by phenotype is the best approach to implementing genomic evaluation using commercial herd data, particularly when pools of 10 individuals are evaluated. While combining phenotypes seems a promising strategy to allow more flexibility to the estimates made using pools, more studies are necessary in this regard.

Genomic analysis of bovine respiratory disease and lung consolidation in preweaned Holstein calves using clinical scoring and lung ultrasound

Bovine respiratory disease (BRD) is a leading cause of morbidity and mortality in dairy calves, with detrimental long-term effects that include stunted growth, increased age at first calving, and decreased milk production in first lactation. The objectives of this study were to establish a protocol for objective and efficient assessment of BRD phenotypes in preweaned dairy calves, develop a genomic reference population with well-defined clinical and subclinical phenotypes, identify chromosomal regions associated with BRD in a genome-wide association study, estimate genetic parameters of BRD, and predict genomic breeding values of dairy calves. A total of 1,107 Holstein calves from 6 dairy farms in southern Wisconsin were examined using clinical respiratory scoring and lung ultrasound at 3 and 6 wk of age. The clinical respiratory score was based on visual appraisal of eyes, nose, ears, cough, and temperature. Lung ultrasound scores were assigned based on the amount of consolidation present. Calves were genotyped with a commercially available SNP array and after quality control and imputation to higher density, 690,291 SNP markers and 1,014 individuals remained. Single-step genome-wide association study and single-step genomic best linear unbiased prediction were applied to binary 3- and 6-wk phenotypes considered as overall respiratory healthy versus affected (RESP3, RESP6) or as presence or absence of lung consolidation (CON3, CON6). Lung ultrasound combined with a clinical scoring system allowed for efficient and objective assessment for the prevalence of BRD. Proportions of variance attributed to 1-Mb non-overlapping windows suggested genomic regions that may contain putative candidate genes, most notably regions on Bos taurus autosomes 1, 6, 7, 10, 11, 12, 15, 17, 18, 27, and 28 that explained 0.70 to 1.45% of the genetic variance. Heritability estimates were higher at 3 wk (0.214 and 0.241 for CON3 and RESP3, respectively) than 6 wk (0.084 and 0.111 for CON6 and RESP6, respectively), and mean reliabilities of genomic estimated breeding vales for calves with genotypes and phenotypes ranged from 0.12 for CON6 to 0.30 for RESP3.

Genome-wide association study of lung lesions and pleurisy in New Zealand lambs

Pneumonia is an important issue for sheep production, leading to reduced growth rate and a predisposition to pleurisy. The objective of this study was to identify loci associated with pneumonic lesions and pleurisy in New Zealand progeny test lambs. The lungs from 3,572 progeny-test lambs were scored for presence and severity of pneumonic lesions and pleurisy at slaughter. Animals were genotyped using the Illumina Ovine Infinium HD SNP BeadChip (606,006 markers). The heritability of lung lesion score and pleurisy were calculated using the genomic relationship matrix, and genome-wide association analyses were conducted using EMMAX and haplotype trend regression. At slaughter, 35% of lambs had pneumonic lesions, with 9% showing lesions on more than half of any individual lobe. The number of lambs recorded as having pleurisy by the processing plants was 9%. Heritability estimates for pneumonic lesions and pleurisy scores adjusted for heteroscedasticity (CPSa and PLEURa) were 0.16 (± 0.03) and 0.05 (± 0.02), respectively. Five single-nucleotide polymorphisms (SNPs) were significantly associated with pneumonic lesions at the genome-wide level, and additional 37 SNPs were suggestively significant. Four SNPs were significantly associated with pleurisy, with an additional 11 SNPs reaching the suggestive level of significance. There were no regions that overlapped between the 2 traits. Multiple SNPs were in regions that contained genes involved in either the DNA damage response or the innate immune response, including several that had previously been reported to have associations with respiratory disease. Both EMMAX and HTR analyses of pleurisy data showed a significant peak on chromosome 2, located downstream from the transcription factor SP3. SP3 activates or suppresses the expression of numerous genes, including several genes with known functions in the immune system. This study identified several SNPs associated with genes involved in both the innate immune response and the response to DNA damage that are associated with pneumonic lesions and pleurisy in lambs at slaughter. Additionally, the identification in sheep of several SNPs within genes that have previously been associated with the respiratory system in cattle, pigs, rats, and mice indicates that there may be common pathways that underlie the response to invasion by respiratory pathogens in multiple species.

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.

Human-Mediated Introgression of Haplotypes in a Modern Dairy Cattle Breed

Domestic animals can serve as model systems of adaptive introgression and their genomic signatures. In part, their usefulness as model systems is due to their well-known histories. Different breeding strategies such as introgression and artificial selection have generated numerous desirable phenotypes and superior performance in domestic animals. The modern Danish Red Dairy Cattle is studied as an example of an introgressed population. It originates from crossing the traditional Danish Red Dairy Cattle with the Holstein and Brown Swiss breeds, both known for high milk production. This crossing happened, among other things due to changes in the production system, to raise milk production and overall performance. The genomes of modern Danish Red Dairy Cattle are heavily influenced by regions introgressed from the Holstein and Brown Swiss breeds and under subsequent selection in the admixed population. The introgressed proportion of the genome was found to be highly variable across the genome. Haplotypes introgressed from Holstein and Brown Swiss contained or overlapped known genes affecting milk production, as well as protein and fat content (CD14, ZNF215, BCL2L12, and THRSP for Holstein origin and ITPR2, BCAT1, LAP3, and MED28 for Brown Swiss origin). Genomic regions with high introgression signals also contained genes and enriched QTL associated with calving traits, body confirmation, feed efficiency, carcass, and fertility traits. These introgressed signals with relative identity-by-descent scores larger than the median showing Holstein or Brown Swiss introgression are mostly significantly correlated with the corresponding test statistics from signatures of selection analyses in modern Danish Red Dairy Cattle. Meanwhile, the putative significant introgressed signals have a significant dependency with the putative significant signals from signatures of selection analyses. Artificial selection has played an important role in the genomic footprints of introgression in the genome of modern Danish Red Dairy Cattle. Our study on a modern cattle breed contributes to an understanding of genomic consequences of selective introgression by demonstrating the extent to which adaptive effects contribute to shape the specific genomic consequences of introgression.

Identification of polymorphisms in the bovine collagenous lectins and their association with infectious diseases in cattle

Infectious diseases are a significant issue in animal production systems, including both the dairy and beef cattle industries. Understanding and defining the genetics of infectious disease susceptibility in cattle is an important step in the mitigation of their impact. Collagenous lectins are soluble pattern recognition receptors that form an important part of the innate immune system, which serves as the first line of host defense against pathogens. Polymorphisms in the collagenous lectin genes have been shown in previous studies to contribute to infectious disease susceptibility, and in cattle, mutations in two collagenous lectin genes (MBL1 and MBL2) are associated with mastitis. To further characterize the contribution of variation in the bovine collagenous lectins to infectious disease susceptibility, we used a pooled NGS approach to identify short nucleotide variants (SNVs) in the collagenous lectins (and regulatory DNA) of cattle with (n = 80) and without (n = 40) infectious disease. Allele frequency analysis identified 74 variants that were significantly (p < 5 × 10⁻⁶) associated with infectious disease, the majority of which were clustered in a 29-kb segment upstream of the collectin locus on chromosome 28. In silico analysis of the functional effects of all the variants predicted 11 SNVs with a deleterious effect on protein structure and/or function, 148 SNVs that occurred within potential transcription factor binding sites, and 31 SNVs occurring within potential miRNA binding elements. This study provides a detailed look at the genetic variation of the bovine collagenous lectins and identifies potential genetic markers for infectious disease susceptibility.

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.

Feed Intake and Weight Changes in Bos indicus-Bos taurus Crossbred Steers Following Bovine Viral Diarrhea Virus Type 1b Challenge Under Production Conditions

Bovine viral diarrhea virus (BVDV) has major impacts on beef cattle production worldwide, but the understanding of host animal genetic influence on illness is limited. This study evaluated rectal temperature, weight change and feed intake in Bos indicus crossbred steers (n = 366) that were challenged with BVDV Type 1b, and where family lines were stratified across three vaccine treatments of modified live (MLV), killed, (KV) or no vaccine (NON). Pyrexia classification based on 40.0 °C threshold following challenge and vaccine treatment were investigated for potential interactions with sire for weight change and feed intake following challenge. Pyrexia classification affected daily feed intake (ADFI, p = 0.05), and interacted with day (p < 0.001) for ADFI. Although low incidence of clinical signs was observed, there were marked reductions in average daily gain (ADG) and cumulative feed intake during the first 14 day post-challenge; ADG (CV of 104%) and feed efficiency were highly variable in the 14-day period immediately post-challenge as compared to the subsequent 14-day periods. A sire × vaccine strategy interaction affected ADFI (p < 0.001), and a sire by time period interaction affected ADG (p = 0.03) and total feed intake (p = 0.03). This study demonstrates that different coping responses may exist across genetic lines to the same pathogen, and that subclinical BVDV infection has a measurable impact on cattle production measures.

Genomewide association study of liver abscess in beef cattle

Fourteen percent of U.S. cattle slaughtered in 2011 had liver abscesses, resulting in reduced carcass weight, quality, and value. Liver abscesses can result from a common bacterial cause, , which inhabits rumen lesions caused by acidosis and subsequently escapes into the blood stream, is filtered by the liver, and causes abscesses in the liver. Our aim was to identify SNP associated with liver abscesses in beef cattle. We used lung samples as a DNA source because they have low economic value, they have abundant DNA, and we had unrestricted access to sample them. We collected 2,304 lung samples from a beef processing plant: 1,152 from animals with liver abscess and 1,152 from animals without liver abscess. Lung tissue from pairs of animals, 1 with abscesses and another without, were collected from near one another on the viscera table to ensure that pairs of phenotypically extreme animals came from the same lot. Within each phenotype (abscess or no abscess), cattle were pooled by slaughter sequence into 12 pools of 96 cattle for each phenotype for a total of 24 pools. The pools were constructed by equal volume of frozen lung tissue from each animal. The DNA needed to allelotype each pool was then extracted from pooled lung tissue and the BovineHD Bead Array (777,962 SNP) was run on all 24 pools. Total intensity (TI), an indicator of copy number variants, was the sum of intensities from red and green dyes. Pooling allele frequency (PAF) was red dye intensity divided TI. Total intensity and PAF were weighted by the inverse of their respective genomic covariance matrices computed over all SNP across the genome. A false discovery rate ≤ 5% was achieved for 15 SNP for PAF and 20 SNP for TI. Genes within 50 kbp from significant SNP were in diverse pathways including maintenance of pH homeostasis in the gastrointestinal tract, maintain immune defenses in the liver, migration of leukocytes from the blood into infected tissues, transport of glutamine into the kidney in response to acidosis to facilitate production of bicarbonate to increase pH, aggregate platelets to liver injury to facilitate liver repair, and facilitate axon guidance. Evidence from the 35 detected SNP associations combined with evidence of polygenic variation indicate that there is adequate genetic variation in incidence rate of liver abscesses, which could be exploited to select sires for reduced susceptibility to subacute acidosis and associated liver abscess.

Use of Genomic Tools to Improve Cattle Health in the Context of Infectious Diseases

Although infectious diseases impose a heavy economic burden on the cattle industry, the etiology of many disorders that affect livestock is not fully elucidated, and effective countermeasures are often lacking. The main tools available until now have been vaccines, antibiotics and antiparasitic drugs. Although these have been very successful in some cases, the appearance of parasite and microbial resistance to these treatments is a cause of concern. Next-generation sequencing provides important opportunities to tackle problems associated with pathogenic illnesses. This review describes the rapid gains achieved to track disease progression, identify the pathogens involved, and map pathogen interactions with the host. Use of novel genomic tools subsequently aids in treatment development, as well as successful creation of breeding programs aimed toward less susceptible livestock. These may be important tools for mitigating the long term effects of combating infection and helping reduce the reliance on antibiotic treatment.