G. L. Bennett

Linkage mapping bovine EST-based SNP

BackgroundExisting linkage maps of the bovine genome primarily contain anonymous microsatellite markers. These maps have proved valuable for mapping quantitative trait loci (QTL) to broad regions of the genome, but more closely spaced markers are needed to fine-map QTL, and markers associated with genes and annotated sequence are needed to identify genes and sequence variation that may explain QTL.ResultsBovine expressed sequence tag (EST) and bacterial artificial chromosome (BAC)sequence data were used to develop 918 single nucleotide polymorphism (SNP) markers to map genes on the bovine linkage map. DNA of sires from the MARC reference population was used to detect SNPs, and progeny and mates of heterozygous sires were genotyped. Chromosome assignments for 861 SNPs were determined by twopoint analysis, and positions for 735 SNPs were established by multipoint analyses. Linkage maps of bovine autosomes with these SNPs represent 4585 markers in 2475 positions spanning 3058 cM . Markers include 3612 microsatellites, 913 SNPs and 60 other markers. Mean separation between marker positions is 1.2 cM. New SNP markers appear in 511 positions, with mean separation of 4.7 cM. Multi-allelic markers, mostly microsatellites, had a mean (maximum) of 216 (366) informative meioses, and a mean 3-lod confidence interval of 3.6 cM Bi-allelic markers, including SNP and other marker types, had a mean (maximum) of 55 (191) informative meioses, and were placed within a mean 8.5 cM 3-lod confidence interval. Homologous human sequences were identified for 1159 markers, including 582 newly developed and mapped SNP.ConclusionAddition of these EST- and BAC-based SNPs to the bovine linkage map not only increases marker density, but provides connections to gene-rich physical maps, including annotated human sequence. The map provides a resource for fine-mapping quantitative trait loci and identification of positional candidate genes, and can be integrated with other data to guide and refine assembly of bovine genome sequence. Even after the bovine genome is completely sequenced, the map will continue to be a useful tool to link observable phenotypes and animal genotypes to underlying genes and molecular mechanisms influencing economically important beef and dairy traits.

Partial-genome evaluation of postweaning feed intake and efficiency of crossbred beef cattle

The effects of individual SNP and the variation explained by sets of SNP associated with DMI, metabolic midtest BW, BW gain, and feed efficiency, expressed as phenotypic and genetic residual feed intake, were estimated from BW and the individual feed intake of 1,159 steers on dry lot offered a 3.0 Mcal/kg ration for at least 119 d before slaughter. Parents of these F(1) × F(1) (F(1)(2)) steers were AI-sired F(1) progeny of Angus, Charolais, Gelbvieh, Hereford, Limousin, Red Angus, and Simmental bulls mated to US Meat Animal Research Center Angus, Hereford, and MARC III composite females. Steers were genotyped with the BovineSNP50 BeadChip assay (Illumina Inc., San Diego, CA). Effects of 44,163 SNP having minor allele frequencies >0.05 in the F(1)(2) generation were estimated with a mixed model that included genotype, breed composition, heterosis, age of dam, and slaughter date contemporary groups as fixed effects, and a random additive genetic effect with recorded pedigree relationships among animals. Variance in this population attributable to sets of SNP was estimated with models that partitioned the additive genetic effect into a polygenic component attributable to pedigree relationships and a genotypic component attributable to genotypic relationships. The sets of SNP evaluated were the full set of 44,163 SNP and subsets containing 6 to 40,000 SNP selected according to association with phenotype. Ninety SNP were strongly associated (P < 0.0001) with at least 1 efficiency or component trait; these 90 accounted for 28 to 46% of the total additive genetic variance of each trait. Trait-specific sets containing 96 SNP having the strongest associations with each trait explained 50 to 87% of additive variance for that trait. Expected accuracy of steer breeding values predicted with pedigree and genotypic relationships exceeded the accuracy of their sires predicted without genotypic information, although gains in accuracy were not sufficient to encourage that performance testing be replaced by genotyping and genomic evaluations.

Use of bovine EST data and human genomic sequences to map 100 gene-specific bovine markers

Abstract. A system to use bovine EST data in conjunction with human genomic sequence to improve the bovine linkage map over the entire genome or on specific chromosomes was evaluated. Bovine EST sequence was used to provide primer sequences corresponding to bovine genes, while human genomic sequence directed primer design to flank introns and produce amplicons of appropriate size for efficient direct sequencing. The sequence tagged sites (STS) produced in this way from the four sires of the MARC reference families were examined for single nucleotide polymorphisms (SNPs) that could be used to map the corresponding genes. With this approach, along with a primer/extension mass spectrometry SNP genotyping assay, 100 ESTs were placed on the bovine genetic linkage map. The first 70 were chosen at random from bovine EST–human genomic comparisons. An additional 30 ESTs were successfully mapped to bovine Chromosome 19 (BTA19), and comparison of the resulting BTA19 map to the position of the corresponding human orthologs on the HSA17 draft sequences revealed differences in the spacing and order of genes. Over 80% of successful amplicons contained SNPs, indicating that this is an efficient approach to generating EST-associated genetic markers. We have demonstrated the feasibility of constructing a linkage map based on SNPs associated with ESTs and the plausibility of utilizing EST, comparative mapping information, and human sequence data to target regions of the bovine genome for SNP marker development.

Association of bovine neonatal Fc receptor a-chain gene (FCGRT) haplotypes with serum IgG concentration in newborn calves

This report describes allelic variation in FCGRT (which encodes the a-chain of FcRn) and its association with variation of IgG concentration in neonatal calves. Five SNPs were identified by sequencing 1305 bp of FCGRT genomic DNA from a multi-breed panel of 96 cattle and 27 founders of a reference population. These SNPs defined five FCGRT haplotypes that were verified by segregation and used to test association of FCGRT with neonatal IgG concentration in a case-control study. This study established that dams with FCGRT haplotype 3 had a significantly greater risk of failure of passive transfer in their calves (odds ratio [OR] = 3.80, CI95% 1.10–13.18, p = 0.035). Calves with FCGRT haplotype 2 were less likely to have high levels of passively acquired immunoglobulin (OR = 0.18, CI95% 0.05–0.68, p = 0.011). These results indicate that the bovine FCGRT haplotype markers are in linkage disequilibrium with genetic risk factors affecting passive transfer of IgG in beef cattle, an important determinant of neonatal calf morbidity and mortality.

Comparative map alignment of BTA27 and HSA4 and 8 to identify conserved segments of genome containing fat deposition QTL

Abstract. Quantitative trait loci (QTL) associated with fat deposition have been identified on bovine Chromosome 27 (BTA27) in two different cattle populations. To generate more informative markers for verification and refinement of these QTL-containing intervals, we initiated construction of a BTA27 comparative map. Fourteen genes were selected for mapping based on previously identified regions of conservation between the cattle and human genomes. Markers were developed from the bovine orthologs of genes found on human Chromosomes 1 (HSA1), 4, 8, and 14. Twelve genes were mapped on the bovine linkage map by using markers associated with single nucleotide polymorphisms or microsatellites. Seven of these genes were also anchored to the physical map by assignment of fluorescence in situ hybridization probes. The remaining two genes not associated with an identifiable polymorphism were assigned only to the physical map. In all, seven genes were mapped to BTA27. Map information generated from the other seven genes not syntenic with BTA27 refined the breakpoint locations of conserved segments between species and revealed three areas of disagreement with the previous comparative map. Consequently, portions of HSA1 and 14 are not conserved on BTA27, and a previously undefined conserved segment corresponding to HSA8p22 was identified near the pericentromeric region of BTA8. These results show that BTA27 contains two conserved segments corresponding to HSA8p, which are separated by a segment corresponding to HSA4q. Comparative map alignment strongly suggests the conserved segment orthologous to HSA8p21-q11 contains QTL for fat deposition in cattle.

The incidence of Escherichia coli on beef carcasses and its association with aerobic mesophilic plate count categories during the slaughter process.

An analysis of 535 prefabricated beef carcass samples taken in three processing plants demonstrated an association between the mesophilic aerobic plate count (APC) class and the incidence of obtaining an Escherichia coli-positive sample. Beef carcasses were sampled from three separate plants; one was a fed-beef processing plant and the other two were cow/bull plants. Samples were obtained by sponging and were analyzed for APC and E. coli. When samples were classified into four APC levels or classes (class 1: < 2, class 2: > or = 2 and < 3, class 3: > or = 3 and < 4, and class 4: > or = 4 log CFU/cm2), a trend indicating that samples from higher APC classes were more likely to be positive for E. coli biotype 1 was observed. Of the APC class 4 samples (> or = 4 log CFU/cm2), 88% were positive for the presence of E. coli, as opposed to 21% in APC class 1 (< 2 log CFU/cm2). Univariate chi-square analysis of the resulting contingency tables from reclassified data (class 1: < 2, class 2: > or = 2 and < 3, and class 3: > or = 3 log CFU/cm2) indicated a strong association between APC class and the incidence (presence or absence) of an E. coli-positive sample. Using multivariate analysis to account for influences of plant and within plant processing site, the data indicated a strong positive linear trend between the presence of E. coli and the APC class.

The change in differing leukocyte populations during vaccination to bovine respiratory disease and their correlations with lung scores, health records, and average daily gain

Bovine respiratory disease (BRD) is the most economically important disease in U.S. feedlots. Infection can result in morbidity, mortality, and reduced average daily gain. Cheap and reliable genetic methods of prediction and protection from BRD would be highly advantageous to the industry. The immune response may correlate with BRD incidence. Cattle (n = 2,182) were vaccinated against common viral and bacterial pathogens of BRD. Two blood samples were collected, one during booster vaccination and one 21d later, enabling 3 phenotypes for each trait [prebooster (pre), postbooster (post), and delta (post minus pre)]. From the blood samples innate and adaptive responses [counts of white blood cells (WBC), neutrophils, lymphocytes, monocytes, eosinophils, and basophils] were measured. In addition, feedlot ADG and binary traits [health records (HR; 0 = healthy, 1 = ill) and lung scores (LS; collected at harvest; 0 = no lesions, 1 = lesions)] were also recorded. Traits ADG, HR, and LS have all been significantly correlated with infection to BRD. In this investigation we aimed to find correlations between the immune response and ADG, HR, and LS to find an easily measurable trait that would be a good predictor of BRD resistance after vaccination. The results showed an average positive delta for the innate immune response (eosinophils, basophils, neutrophils), whereas the adaptive immune response had an average negative delta (lymphocytes). Overall, we discovered that the immune responses had moderately high heritabilities (h(2); lowest: delta monocytes, 0.21 ± 0.05; greatest: pre lymphocytes: 0.5 ± 0.05), with lymphocytes having the greatest h(2) throughout the study (h(2) ≥ 0.41). All genetic correlations were calculated using bivariate REML models. Although LS did not significantly correlate with any of the immune phenotypes, both ADG (post lymphocytes, -0.24 ± 0.12) and HR (pre eosinophils, -0.67 ± 0.29; delta WBC, -0.5 ± 0.24, and delta lymphocytes, -0.67 ± 0.21) did. All the significant genetic correlations with HR were negative; resistance to BRD appears to be a function of greater delta lymphocytes and WBC. The increase in eosinophils may potentially link its role in decreasing lymphocytes. These results may enable producers to predict if revaccination, quarantine, and breeding of animals is required to reduce the incidence of BRD postvaccination. In addition, immunological phenotypes maybe used to aid genomic selection indices to select animals with greater rates of protection after BRD vaccination.

Physical and linkage mapping of mammary-derived expressed sequence tags in cattle

This study describes the physical and linkage mapping of 42 gene-associated markers developed from mammary gland-derived expressed sequence tags to the cattle genome. Of the markers, 25 were placed on the USDA reference linkage map and 37 were positioned on the Roslin 3000-rad radiation hybrid (RH) map, with 20 assignments shared between the maps. Although no novel regions of conserved synteny between the cattle and the human genomes were identified, the coverage was extended for bovine chromosomes 3, 7, 15, and 29 compared with previously published comparative maps between human and bovine genomes. Overall, these data improve the resolution of the human-bovine comparative maps and will assist future efforts to integrate bovine RH and linkage map data.