A look under the hood of genomic-estimated breed compositions for brangus cattle: What have we learned?

The Brangus cattle were developed to utilize the superior traits of Angus and Brahman cattle. Their genetic compositions are expected to be stabilized at 3/8 Brahman and 5/8 Angus. Previous studies have shown more than expected Angus lineage with Brangus cattle, and the reasons are yet to be investigated. In this study, we revisited the breed compositions for 3,605 Brangus cattle from three perspectives: genome-wise (GBC), per chromosomes (CBC), and per chromosome segments (SBC). The former (GBC) depicted an overall picture of the "mosaic" genome of the Brangus attributable to their ancestors, whereas the latter two criteria (CBC and SBC) corresponded to local ancestral contributions. The average GBC for the 3,605 Brangus cattle were 70.2% Angus and 29.8% Brahman. The K-means clustering supported the postulation of the mixture of 1/2 Ultrablack (UB) animals in Brangus. For the non-UB Brangus animals, the average GBC were estimated to be 67.4% Angus and 32.6% Brahman. The 95% confidence intervals of their overall GBC were 60.4%-73.5% Angus and 26.5%-39.6% Brahman. Possibly, genetic selection and drifting have resulted in an approximately 5% average deviation toward Angus lineage. The estimated ancestral contributions by chromosomes were heavily distributed toward Angus, with 27 chromosomes having an average Angus CBC greater than 62.5% but only two chromosomes (5 and 20) having Brahman CBC greater than 37.5%. The chromosomal regions with high Angus breed proportions were prevalent, tending to form larger blocks on most chromosomes. In contrast, chromosome segments with high Brahman breed proportion were relatively few and isolated, presenting only on seven chromosomes. Hence, genomic hitchhiking effects were strong where Angus favorable alleles resided but weak where Brahman favorable alleles were present. The functions of genes identified in the chromosomal regions with high ( ≥ 75 % ) Angus compositions were diverse yet may were related to growth and body development. In contrast, the genes identified in the regions with high ( ≥ 37.5 % ) Brahman compositions were primarily responsible for disease resistance. In conclusion, we have addressed the questions concerning the Brangus genetic make-ups. The results can help form a dynamic picture of the Brangus breed formation and the genomic reshaping.

Relationship of molecular breeding value for beef tenderness with heifer traits through weaning of their first calf.

Polymorphisms in µ-calpain (CAPN1) that beneficially associate with beef tenderness are reported to antagonistically associate with calving day in beef heifers and post-partum interval to estrus in beef cows. We, therefore, hypothesized that a molecular breeding value for slice shear force, calculated based on CAPN1 and calpastatin (CAST) genotypes, would demonstrate an antagonistic relationship between genomically predicted slice shear force and ordinal calving date in replacement beef heifers. A secondary objective of this study was to evaluate the association of a polymorphism in diacylglycerol O-acyltransferase (DGAT1) with reproductive traits in beef heifers. One hundred eighty-seven MARC III heifers (» Angus, » Hereford, » Red Poll, and » Pinzgauer) that had been selectively bred to increase the frequency of these polymorphisms were submitted for monthly ultrasound exams beginning at 333 d of age and continuing until the start of breeding to determine pubertal status. At the last exam before breeding, all antral follicles were counted, and the length and height of each ovary was measured to determine if genomic selection for slice shear force associated with ovarian follicle number. Calving date, calf gender, and calf birth weight were recorded at parturition. Regression analysis of the molecular breeding value for slice shear force of the heifers on ordinal calving date indicated no association between genomic prediction of tenderness and calving date (P = 0.16); however, there was a tendency for age at puberty to be delayed in heifers as genetic merit for tenderness improved (P = 0.09). The results of the present study indicate that within experimental precision, selecting for tenderness using genomic predictions had minimal or no antagonistic association with reproductive performance in heifers. Further analysis of reproductive performance as cows is needed within this population but applying these genetic markers to select for tenderness in steers does not antagonize reproductive traits influencing conception or first calf birth date and birth weight in replacement beef heifers.

A Causality Perspective of Genomic Breed Composition for Composite Animals.

Genomic breed composition (GBC) of an individual animal refers to the partition of its genome according to the inheritance from its ancestors or ancestral breeds. For crossbred or composite animals, knowing their GBC is useful to estimate heterosis, to characterize their actual inheritance from foundation breeds, and to make management decisions for crossbreeding programs. Various statistical approaches have been proposed to estimate GBC in animals, but the interpretations of estimates have varied with these methods. In the present study, we proposed a causality interpretation of GBC based on path analysis. We applied this method to estimating GBC in two composite breeds of beef cattle, namely Brangus and Beefmaster. Three SNP panels were used to estimate GBC: a 10K SNP panel consisting of 10,226 common SNPs across three GeneSeek Genomic Profiler (GGP) bovine SNP arrays (GGP 30K, GGP 40K, and GGP 50K), and two subsets (1K and 5K) of uniformly distributed SNPs. The path analysis decomposed the relationships between the ancestors and the composite animals into direct and indirect path effects, and GBC was measured by the relative ratio of the coefficients of direct (D-GBC) and combined (C-GBC) effects from each ancestral breed to the progeny, respectively. Estimated GBC varied only slightly between different genotyping platforms and between the three SNP panels. In the Brangus cattle, because the two ancestral breeds had a very distant relationship, the estimated D-GBC and C-GBC were comparable to each other in the path analysis, and they corresponded roughly to the estimated GBC from the linear regression and the admixture model. In the Beefmaster, however, the strong relationship in allelic frequencies between Hereford and Shorthorn imposed a challenge for the linear regression and the admixture model to estimated GBC reliably. Instead, D-GBC by the path analysis included only direct ancestral effects, and it was robust to bias due to high genomic correlations between reference (ancestral) breeds.

Genomic mating as sustainable breeding for Chinese indigenous Ningxiang pigs.

An important economic reason for the loss of local breeds is that they tend to be less productive, and hence having less market value than commercial breeds. Nevertheless, local breeds often have irreplaceable values, genetically and sociologically. In the breeding programs with local breeds, it is crucial to balance the selection for genetic gain and the maintaining of genetic diversity. These two objectives are often conflicting, and finding the optimal point of the trade-off has been a challenge for breeders. Genomic selection (GS) provides a revolutionary tool for the genetic improvement of farm animals. At the same time, it can increase inbreeding and produce a more rapid depletion of genetic variability of the selected traits in future generations. Optimum-contribution selection (OCS) represents an approach to maximize genetic gain while constraining inbreeding within a targeted range. In the present study, 515 Ningxiang pigs were genotyped with the Illumina Porcine SNP60 array or the GeneSeek Genomic Profiler Porcine 50K array. The Ningxiang pigs were found to be highly inbred at the genomic level. Average locus-wise inbreeding coefficients were 0.41 and 0.37 for the two SNP arrays used, whereas genomic inbreeding coefficients based on runs of homozygosity were 0.24 and 0.25, respectively. Simulated phenotypic data were used to assess the utility of genomic OCS (GOCS) in comparison with GS without inbreeding control. GOCS was conducted under two scenarios, selecting sires only (GOCS_S) or selecting sires and dams (GOCS_SD), while kinships were constrained on selected parents. The genetic gain for average daily body weight gain (ADG) per generation was between 18.99 and 20.55 g with GOCS_S, and between 23.20 and 28.92 with GOCS_SD, and it varied from 25.38 to 48.38 g under GS without controlling inbreeding. While the rate of genetic gain per generation obtained using GS was substantially larger than that obtained by the two scenarios of genomic OCS in the beginning generations of selection, the difference in the genetic gain of ADG between GS and GOCS reduced quickly in latter generations. At generation ten, the difference in the realized rates of genetic gain between GS and GOCS_SD diminished and ended up with even a slightly higher genetic gain with GOCS_SD, due to the rapid loss of genetic variance with GS and fixation of causative genes. The rate of inbreeding was mostly maintained below 5% per generation with genomic OCS, whereas it increased to between 10.5% and 15.3% per generation with GS. Therefore, genomic OCS appears to be a sustainable strategy for the genetic improvement of local breeds such as Ningxiang pigs, but keeping mind that a variety of GOCS methods exist and the optimal forms remain to be exploited further.

Estimation of Genomic Breed Composition for Purebred and Crossbred Animals Using Sparsely Regularized Admixture Models.

A variety of statistical methods, such as admixture models, have been used to estimate genomic breed composition (GBC). These methods, however, tend to produce non-zero components to reference breeds that shared some genomic similarity with a test animal. These non-essential GBC components, in turn, offset the estimated GBC for the breed to which it belongs. As a result, not all purebred animals have 100% GBC of their respective breeds, which statistically indicates an elevated false-negative rate in the identification of purebred animals with 100% GBC as the cutoff. Otherwise, a lower cutoff of estimated GBC will have to be used, which is arbitrary, and the results are less interpretable. In the present study, three admixture models with regularization were proposed, which produced sparse solutions through suppressing the noise in the estimated GBC due to genomic similarities. The regularization or penalty forms included the L1 norm penalty, minimax concave penalty (MCP), and smooth clipped absolute deviation (SCAD). The performances of these regularized admixture models on the estimation of GBC were examined in purebred and composite animals, respectively, and compared to that of the non-regularized admixture model as the baseline model. The results showed that, given optimal values for λ, the three sparsely regularized admixture models had higher power and thus reduced the false-negative rate for the breed identification of purebred animals than the non-regularized admixture model. Of the three regularized admixture models, the two with a non-convex penalty outperformed the one with L1 norm penalty. In the Brangus, a composite cattle breed, estimated GBC were roughly comparable among the four admixture models, but all the four models underestimated the GBC for these composite animals when non-ancestral breeds were included as the reference. In conclusion, the admixture models with sparse regularization gave more parsimonious, consistent and interpretable results of estimated GBC for purebred animals than the non-regularized admixture model. Nevertheless, the utility of regularized admixture models for estimating GBC in crossbred or composite animals needs to be taken with caution.

Heat stress-induced deficits in growth, metabolic efficiency, and cardiovascular function coincided with chronic systemic inflammation and hypercatecholaminemia in ractopamine-supplemented feedlot lambs.

Heat stress hinders growth and well-being in livestock, an effect that is perhaps exacerbated by the ß1 agonist ractopamine. Heat stress deficits are mediated in part by reduced feed intake, but other mechanisms involved are less understood. Our objective was to determine the direct impact of heat stress on growth and well-being in ractopamine-supplemented feedlot lambs. Commercial wethers were fed under heat stress (40 °C) for 30 d, and controls (18 °C) were pair-fed. In a 2 × 2 factorial, lambs were also given a daily gavage of 0 or 60 mg ractopamine. Growth, metabolic, cardiovascular, and stress indicators were assessed throughout the study. At necropsy, 9th to 12th rib sections (four-rib), internal organs, and feet were assessed, and sartorius muscles were collected for ex vivo glucose metabolic studies. Heat stress increased (P < 0.05) rectal temperatures and respiration rates throughout the study and reduced (P < 0.05) weight gain and feed efficiency over the first week, ultrasonic loin-eye area and loin depth near the end of the study, and four-rib weight at necropsy. Fat content of the four-rib and loin were also reduced (P < 0.05) by heat stress. Ractopamine increased (P < 0.05) loin weight and fat content and partially moderated the impact of heat stress on rectal temperature and four-rib weight. Heat stress reduced (P < 0.05) spleen weight, increased (P < 0.05) adrenal and lung weights, and was associated with hoof wall overgrowth but not organ lesions. Ractopamine did not affect any measured indicators of well-being. Heat stress reduced (P < 0.05) blood urea nitrogen and increased (P < 0.05) circulating monocytes, granulocytes, and total white blood cells as well as epinephrine, TNFα, cholesterol, and triglycerides. Cortisol and insulin were not affected. Heat stress reduced (P < 0.05) blood pressure and heart rates in all lambs and increased (P < 0.05) left ventricular wall thickness in unsupplemented but not ractopamine-supplemented lambs. No cardiac arrhythmias were observed. Muscle glucose uptake did not differ among groups, but insulin-stimulated glucose oxidation was reduced (P < 0.05) in muscle from heat-stressed lambs. These findings demonstrate that heat stress impairs growth, metabolism, and well-being even when the impact of feed intake is eliminated by pair-feeding and that systemic inflammation and hypercatecholaminemia likely contribute to these deficits. Moreover, ractopamine improved muscle growth indicators without worsening the effects of heat stress.

Enhanced estimates of carcass and meat quality effects for polymorphisms in myostatin and µ-calpain genes.

The objective of this study was to enhance estimates of additive, dominance, and epistatic effects of marker polymorphisms on beef carcass and quality traits. Myostatin (MSTN) F94L SNP and the µ-calpain (CAPN1) 316 and 4751 SNP haplotype have previously been associated with fat and muscle traits in beef cattle. Multiyear selection in a composite population segregating these polymorphisms increased minor allele (F94L L) and chosen haplotype (CAPN1 CC and GT) frequencies to intermediate levels resulting in more precise estimates of additive and nonadditive genetic effects. During the 3 yr after selection, 176 steers were evaluated for growth, carcass, meat quality, tenderness (n = 103), and meat color traits. The statistical model included year, age of dam, age of the steer, and genotype in a random animal model. The 9 genotypes (3 CAPN1 diplotypes × 3 F94L genotypes) affected marbling score, ribeye area, adjusted fat thickness, vision yield grade (all P < 0.001), slice shear force (P = 0.03), and CIE L* reflectance (P = 0.01). Linear contrasts of the 9 genotypes estimated additive, recessive, and epistatic genetic effects. Significant additive effects of the F94L L allele decreased marbling score, adjusted fat thickness, vision yield grade, and slice shear force; and increased ribeye area and CIE L* reflectance. The homozygous F94L FF and LL genotypes differed by 1.3 to 1.9 phenotypic SD for most carcass traits and by 0.8 to 0.9 SD for slice shear force and CIE L* reflectance but carcass weight differed by only 3 kg (0.1 SD). The L allele was partially recessive to F for ribeye area (P = 0.02) and the heterozygous FL means tended to be closer to the FF genotype than the LL genotype for other carcass traits but differences from additive were not significant. The CAPN1 additive × F94L additive effect on slice shear force was the only significant epistatic estimate. The F94L L allele is prevalent in Limousin but nearly absent in other U.S. purebreds. This allele had about half of the effects on birth weight, muscle, and fat traits reported for severe MSTN mutations in Belgian Blue and Piedmontese breeds. The interaction between MSTN and CAPN1 genotypes may reflect the strong additive effects of MSTN F94L L allele on fat and muscle traits interfering with the phenotypic effect of CAPN1 genotype on meat tenderness.

Comparing SNP panels and statistical methods for estimating genomic breed composition of individual animals in ten cattle breeds.

BACKGROUND: SNPs are informative to estimate genomic breed composition (GBC) of individual animals, but selected SNPs for this purpose were not made available in the commercial bovine SNP chips prior to the present study. The primary objective of the present study was to select five common SNP panels for estimating GBC of individual animals initially involving 10 cattle breeds (two dairy breeds and eight beef breeds). The performance of the five common SNP panels was evaluated based on admixture model and linear regression model, respectively. Finally, the downstream implication of GBC on genomic prediction accuracies was investigated and discussed in a Santa Gertrudis cattle population. RESULTS: There were 15,708 common SNPs across five currently-available commercial bovine SNP chips. From this set, four subsets (1,000, 3,000, 5,000, and 10,000 SNPs) were selected by maximizing average Euclidean distance (AED) of SNP allelic frequencies among the ten cattle breeds. For 198 animals presented as Akaushi, estimated GBC of the Akaushi breed (GBCA) based on the admixture model agreed very well among the five SNP panels, identifying 166 animals with GBCA = 1. Using the same SNP panels, the linear regression approach reported fewer animals with GBCA = 1. Nevertheless, estimated GBCA using both models were highly correlated (r = 0.953 to 0.992). In the genomic prediction of a Santa Gertrudis population (and crosses), the results showed that the predictability of molecular breeding values using SNP effects obtained from 1,225 animals with no less than 0.90 GBC of Santa Gertrudis (GBCSG) decreased on crossbred animals with lower GBCSG. CONCLUSIONS: Of the two statistical models used to compute GBC, the admixture model gave more consistent results among the five selected SNP panels than the linear regression model. The availability of these common SNP panels facilitates identification and estimation of breed compositions using currently-available bovine SNP chips. In view of utility, the 1 K panel is the most cost effective and it is convenient to be included as add-on content in future development of bovine SNP chips, whereas the 10 K and 16 K SNP panels can be more resourceful if used independently for imputation to intermediate or high-density genotypes.

µ-Calpain (CAPN1), calpastatin (CAST), and growth hormone receptor (GHR) genetic effects on Angus beef heifer performance traits and reproduction.

Genetic marker effects and type of inheritance are estimated with poor precision when minor marker allele frequencies are low. An Angus population was subjected to marker assisted selection for multiple years to equalize CAPN1 haplotypes, CAST, and GHR genetic marker frequencies. The objective was to estimate the pleiotropic effects of these carcass quality oriented markers for body weight, reproduction, and first calf performance traits in 174 replacement beef females which were managed under 2 post-weaning development protocols. Heifers were weighed at 11-, 12-, and 13-mo, at first breeding season pregnancy evaluation, and prior to first calving season. Pubertal status was determined at 11-, 12-, and 13-mo of age. Antral follicles were counted, reproductive tracts were scored, and tract dimensions were measured at 13-mo. Body condition and hip height were scored and measured at pregnancy evaluation and prior to calving season. Heifer pregnancy and weaning rates and ordinal birth date were recorded. Calf body weights at birth and weaning were analyzed. Single df linear contrasts for recessive effects of the GHR heterozygous genotype showed significant decreases of 2.5-3.6% in 11-, 12-, and 13-mo heifer body weights and heifer weight prior to calving. The additive differences between GHR homozygotes were small and not significant for all body weights measured but a 1 wk difference in calf birth date was significant. For all 13-mo uterine measurements, scores, and antral follicle counts, only the CAST dominance contrast for medium antral follicle count was significant. The CAPN1 haplotype with a strong additive effect for increased beef tenderness also had a significant additive effect on calving date. Heifers homozygous for the tender haplotype calved 7.9 days later than heifers homozygous for the tough haplotype. Most heifer reproductive traits were not significantly affected by CAST and CAPN1 markers that are widely used to improve beef tenderness by selection and breeders should not be concerned with how these markers affect reproduction and other heifer traits with the possible exception of CAPN1 effects on calving date.

Associations between infectious bovine keratoconjunctivitis at weaning and ultrasongraphically measured body composition traits in yearling cattle.

OBJECTIVE: To assess associations between infectious bovine keratoconjunctivitis (IBK) diagnosed at weaning and production traits in yearling beef calves. DESIGN: Retrospective population-based cohort study. ANIMALS: 1,882 Angus calves. PROCEDURES: Angus calves from 1 farm were evaluated over 7 years. The association between yearling body production traits and detection of IBK lesions at weaning was evaluated. RESULTS: Yearlings that had evidence of IBK at weaning had less 12th rib fat depth, ribeye area, and body weight than did cohorts without evidence of IBK. Average daily gain was greater in cattle that had IBK lesions at weaning, but this did not offset lower body weight at weaning. CONCLUSIONS AND CLINICAL RELEVANCE: The associations between IBK at weaning and production variables persisted well into the postweaning period, and there appeared to be a relationship between decreased body composition traits at yearling evaluation and IBK infection before weaning.

Body composition and gene expression QTL mapping in mice reveals imprinting and interaction effects.

BACKGROUND: Shifts in body composition, such as accumulation of body fat, can be a symptom of many chronic human diseases; hence, efforts have been made to investigate the genetic mechanisms that underlie body composition. For example, a few quantitative trait loci (QTL) have been discovered using genome-wide association studies, which will eventually lead to the discovery of causal mutations that are associated with tissue traits. Although some body composition QTL have been identified in mice, limited research has been focused on the imprinting and interaction effects that are involved in these traits. Previously, we found that Myostatin genotype, reciprocal cross, and sex interacted with numerous chromosomal regions to affect growth traits. RESULTS: Here, we report on the identification of muscle, adipose, and morphometric phenotypic QTL (pQTL), translation and transcription QTL (tQTL) and expression QTL (eQTL) by applying a QTL model with additive, dominance, imprinting, and interaction effects. Using an F2 population of 1000 mice derived from the Myostatin-null C57BL/6 and M16i mouse lines, six imprinted pQTL were discovered on chromosomes 6, 9, 10, 11, and 18. We also identified two IGF1 and two Atp2a2 eQTL, which could be important trans-regulatory elements. pQTL, tQTL and eQTL that interacted with Myostatin, reciprocal cross, and sex were detected as well. Combining with the additive and dominance effect, these variants accounted for a large amount of phenotypic variation in this study. CONCLUSIONS: Our study indicates that both imprinting and interaction effects are important components of the genetic model of body composition traits. Furthermore, the integration of eQTL and traditional QTL mapping may help to explain more phenotypic variation than either alone, thereby uncovering more molecular details of how tissue traits are regulated.

Genome-wide association and prediction of direct genomic breeding values for composition of fatty acids in Angus beef cattle.

BACKGROUND: As consumers continue to request food products that have health advantages, it will be important for the livestock industry to supply a product that meet these demands. One such nutrient is fatty acids, which have been implicated as playing a role in cardiovascular disease. Therefore, the objective of this study was to determine the extent to which molecular markers could account for variation in fatty acid composition of skeletal muscle and identify genomic regions that harbor genetic variation. RESULTS: Subsets of markers on the Illumina 54K bovine SNPchip were able to account for up to 57% of the variance observed in fatty acid composition. In addition, these markers could be used to calculate a direct genomic breeding values (DGV) for a given fatty acids with an accuracy (measured as simple correlations between DGV and phenotype) ranging from -0.06 to 0.57. Furthermore, 57 1-Mb regions were identified that were associated with at least one fatty acid with a posterior probability of inclusion greater than 0.90. 1-Mb regions on BTA19, BTA26 and BTA29, which harbored fatty acid synthase, Sterol-CoA desaturase and thyroid hormone responsive candidate genes, respectively, explained a high percentage of genetic variance in more than one fatty acid. It was also observed that the correlation between DGV for different fatty acids at a given 1-Mb window ranged from almost 1 to -1. CONCLUSIONS: Further investigations are needed to identify the causal variants harbored within the identified 1-Mb windows. For the first time, Angus breeders have a tool whereby they could select for altered fatty acid composition. Furthermore, these reported results could improve our understanding of the biology of fatty acid metabolism and deposition.

Genome-wide association study of infectious bovine keratoconjunctivitis in Angus cattle.

BACKGROUND: Infectious Bovine Keratoconjunctivitis (IBK) in beef cattle, commonly known as pinkeye, is a bacterial disease caused by Moraxellabovis. IBK is characterized by excessive tearing and ulceration of the cornea. Perforation of the cornea may also occur in severe cases. IBK is considered the most important ocular disease in cattle production, due to the decreased growth performance of infected individuals and its subsequent economic effects. IBK is an economically important, lowly heritable categorical disease trait. Mass selection of unaffected animals has not been successful at reducing disease incidence. Genome-wide studies can determine chromosomal regions associated with IBK susceptibility. The objective of the study was to detect single-nucleotide polymorphism (SNP) markers in linkage disequilibrium (LD) with genetic variants associated with IBK in American Angus cattle. RESULTS: The proportion of phenotypic variance explained by markers was 0.06 in the whole genome analysis of IBK incidence classified as two, three or nine categories. Whole-genome analysis using any categorisation of (two, three or nine) IBK scores showed that locations on chromosomes 2, 12, 13 and 21 were associated with IBK disease. The genomic locations on chromosomes 13 and 21 overlap with QTLs associated with Bovine spongiform encephalopathy, clinical mastitis or somatic cell count. CONCLUSIONS: Results of these genome-wide analyses indicated that if the underlying genetic factors confer not only IBK susceptibility but also IBK severity, treating IBK phenotypes as a two-categorical trait can cause information loss in the genome-wide analysis. These results help our overall understanding of the genetics of IBK and have the potential to provide information for future use in breeding schemes.

Whole genome analysis of infectious bovine keratoconjunctivitis in Angus cattle using Bayesian threshold models.

Infectious bovine keratoconjunctivitis (IBK), also known as pinkeye, is characterized by damage to the cornea and is an economically important, lowly heritable, categorical disease trait in beef cattle. Scores of eye damage were collected at weaning on 858 Angus cattle. SNP genotypes for each animal were obtained from BovineSNP50 Infinium-beadchips. Simultaneous associations of all SNP with IBK phenotype were determined using Bayes-C that treats SNP effects as random with equal variance for an assumed fraction (π=0.999) of SNP having no effect on IBK scores. Bayes-C threshold models were used to estimate SNP effects by classifying IBK into two, three or nine ordered categories. Magnitudes of genetic variances estimated in localized regions across the genome indicated that SNP within the most informative regions accounted for much of the genetic variance of IBK and pointed out some degree of association to IBK. There are many candidate genes in these regions which could include a gene or group of genes associated with bacterial disease in cattle.

Association of toll-like receptor four single nucleotide polymorphisms with incidence of infectious bovine keratoconjunctivitis (IBK) in cattle.

Toll-like receptor 4 (TLR4) is a receptor protein that binds pathogen ligands, which are mainly associated with Gram-negative bacteria. The objective of this study was to investigate the association of nucleotide polymorphisms in TLR4 with infectious bovine keratoconjunctivitis (IBK), or pinkeye, incidence in American Angus cattle. Animals with previously calculated breeding values for IBK susceptibility were used to identify two SNPs in TLR4; Int1 (A/G) in intron1 (-26 Ex2 position) and Ex3 (C/T) in exon3 (1,678 position). To investigate the possible role of these SNPs in IBK susceptibility, the disease incidence information was collected on 370 calves raised in Iowa at two time points-June or August (disease season) and October (at weaning) and genotyped using PCR-RFLP protocols. In statistical models including year, pasture management group, and SNP, the Int1 SNP had a significant effect on IBK infection rates both in-season (P < 0.05) and at weaning (P < 0.01), whereas the Ex3 SNP was not significant (P > 0.79) at either time point. Furthermore, the Int1 SNP alone could account for 2.1% of phenotypic variation in IBK infection during the disease season and 3.0% of phenotypic variation in IBK infection at the time of weaning. These data indicate that there is a relationship between Int1 genotype and the rate of IBK infection in American Angus cattle.