Search

Beef × dairy crossbred cattle (n = 615) were used to evaluate the effect of preharvest indicator traits and genotypes on the accuracy of estimated breeding values (EBVs) of seedstock candidates for selection. Genotypes for 100,000 single nucleotide polymorphisms were provided by the American Simmental Association of purebred and crossbred seedstock animals (n = 2,632). Five hundred and ninety-five of the 615 beef × dairy cattle had carcass camera and ultrasound data. Phenotypes were not used for any of the seedstock animals even though some may have had performance and ultrasound data. We estimated the genomic relationship matrix among 3,247 animals including both phenotyped and unphenotyped animals. We computed genetic parameters among 37 traits using 666 bivariate restricted maximum likelihood analyses. The accuracy of EBV depends on heritability. For the sake of brevity, we report accuracy for marbling as a proxy for other traits with similar heritability. We focus on accuracy for marbling because marbling is the primary determinant of carcass value. We computed EBV for all 3,247 animals for marbling based on camera data postharvest using best linear unbiased prediction. We report evidence of overlap in causative genes among postharvest carcass traits; marbling, ribeye area, yield grade, fat thickness, and hot carcass weight (HCW) based on genetic correlations. Genetic correlations range from -0.73 to 0.89. Several live animal traits (frame size, body weight and ultrasound fat thickness and ribeye area) were genetically correlated with postharvest traits; including HCW, ribeye area, yield grade, fat thickness, and marbling. Genetic correlations between pre- and postharvest traits ranged from -0.53 to 0.95. Accuracy for marbling ranged from 0.64 to 0.80 for animals with marbling recorded, and from 0.09 to 0.60 for animals without marbling recorded. The accuracy of animals without phenotypes was related to the genomic relationship between animals with phenotype and those without. Live animal traits were useful for predicting economically important carcass traits based on genetic correlations. The accuracy of EBV for seedstock animals that were not phenotyped was low, but this is consistent with theory, and accuracy is expected to increase with the addition of genotypes and carcass data from beef × dairy animals., (Published by Oxford University Press on behalf of the American Society of Animal Science 2024.)

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)., (Published by Oxford University Press on behalf of the American Society of Animal Science 2023.)

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., (Published by Oxford University Press on behalf of Genetics Society of America 2021.)

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., Competing Interests: The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest., (Copyright © 2021 Abrams, McDaneld, Keele, Chitko-McKown, Kuehn and Gonda.)

One approach to reducing calving difficulty is to select heifers with higher breeding value for calving ease. Calving ease is often associated with lower birth weight and that may result in other possible effects on lifetime productivity. Females from experimental select and control calving ease lines within each of the seven populations were compared. Random samples of 720 heifers from lines selected for better calving ease breeding values and 190 heifers from control lines selected for average birth weights were followed through four parities. Select and control lines within the same population were selected to achieve similar yearling weight breeding values. Weights of sampled heifers in select lines were 2.6 kg ( P < 0.01) lighter at birth but not different from control lines at weaning. Select lines had significantly shorter hip height, lighter mature weight, and greater calving success at second parity. Their calves were born significantly earlier with lighter weights and less assistance. Significant interactions with parity showed fewer calves assisted and greater calf survival to weaning as heifers but negligible differences with control lines in later parities. Steer progeny sampled from these dams in select lines ( n = 204) were not different from steers in control lines ( n = 91) for hot carcass weight but had significantly greater fat depth. Two production systems were compared considering the seven populations as replicates. The systems differed in selection history of females (select and control lines) and the use of bulls within their lines as young cows, but used the same bulls in both lines as older cows. Cows were culled after single unsuccessful breeding and kept for up to four parities. Select line cows tended ( P ≤ 0.10) to wean more calves and stay in the herd longer. They were assisted significantly fewer times at calving and had greater calf weight gain to weaning when evaluated over their herd life. Mature weights were lighter in select lines, but marketable cow weight from the systems was nearly identical. Control lines did have more marketable young cow weight and select lines older cow weight. Weaned calf weight per heifer starting the system was significantly greater for the select heifer system due to greater survival of calves from heifers and greater calving success at second parity. No important unfavorable effects of genetic differences in calving ease were identified in this experiment., (Published by Oxford University Press on behalf of the American Society of Animal Science 2021.)

The U.S. Meat Animal Research Center was the first entity in the United States to import the Romanov breed and it has been maintained as a closed flock for over 30 yr. Incorporating this super-prolific breed into crossbred and composite populations has resulted in large improvements in ewe productivity. However, few have quantified factors contributing to genetic and nongenetic variation in ewe reproduction and lamb growth within purebred Romanov populations, which were the objectives of this study. The pedigree contained a total of 8,683 lambs born to 218 and 1,600 unique sires and dams, respectively. Number of lambs born on a per ewe exposed (NLBE) and lambing (NLBL) basis were analyzed in univariate repeatability animal models. As expected, the proportion of phenotypic variance (σP2) in litter size attributable to additive genetic (0.06 to 0.08) and permanent environmental (0.05 to 0.07) effects of the ewe was low. The service sire permanent environmental effect contributed to a small but significant amount of σP2 in NLBE (0.03) but not NLBL. However, the service sire additive genetic effect did not influence σP2 in NLBE or NLBL. Lamb body weight was recorded at birth (BWB) and upon weaning from either milk replacer (~30 d; BWW-N) or their dam (~60 d; BWW-D) and were analyzed in a three-trait model with random additive direct and maternal effects. Estimated direct heritabilities were low for all body weight (BW) traits (0.07 to 0.10). Maternal heritability was moderate for BWB (0.34) but low for weaning BW (0.11 to 0.18). This was the first to report direct and maternal genetic correlations between BW of nursery- and dam-reared lambs, and both were estimated to be moderate (0.43 to 0.47). Additionally, the direct and maternal effects of BWB were more strongly correlated with BWW-N (0.74 to 0.82) than BWW-D (0.17 to 0.33). Despite inbreeding coefficients having increased at a rate of 0.33% per birth year (1986 to 2019) in this flock, they were not consistently associated with reductions in ewe or lamb performance. Parameter estimates generally agreed with those from less-prolific breeds, and results indicate that selection can be an effective means of improving subcomponents of ewe productivity., (Published by Oxford University Press on behalf of the American Society of Animal Science 2020.)

Surveys of microbial populations in environmental niches of interest often utilize sequence variation in the gene encoding the ribosomal small subunit (the 16S rRNA gene). Generally, these surveys target the 16S genes using semi-degenerate primers to amplify portions of a subset of bacterial species, sequence the amplicons in bulk, and assign to putative taxonomic categories by comparison to databases purporting to connect specific sequences in the main variable regions of the gene to specific organisms. Due to sequence length constraints of the most popular bulk sequencing platforms, the primers selected amplify one to three of the nine variable regions, and taxonomic assignment is based on relatively short stretches of sequence (150-500 bases). We demonstrate that taxonomic assignment is improved through reduced unassigned reads by including a survey of near-full-length sequences specific to the target environment, using a niche of interest represented by the upper respiratory tract (URT) of cattle. We created a custom Bovine URT database from these longer sequences for assignment of shorter, less expensive reads in comparisons of the upper respiratory tract among individual animals. This process improves the ability to detect changes in the microbial populations of a given environment, and the accuracy of defining the content of that environment at increasingly higher taxonomic resolution., Competing Interests: The co-author KD Dedonder is currently employed by Veterinary and Biomedical Research Center, Inc, Manhattan. However, the past work was completed at Kansas State University during their education program. This is reflected in the affiliation statement. The current company provided no support in the form of salary for author KDD, and did not have any role in the study design, data collection and analysis, decision to publish, or preparation of the manuscript. There are no relevant patents or products to declare. This does not affect our adherence to PLOS ONE policies on data or materials sharing.

We hypothesized cattle that differed in BW gain had different digestive tract microbiota. Two experiments were conducted. In both experiments, steers received a diet that consisted of 8.0% chopped alfalfa hay, 20% wet distillers grain with solubles, 67.75% dry-rolled corn, and 4.25% vitamin/mineral mix (including monensin) on a dry matter basis. Steers had ad libitum access to feed and water. In experiment 1, 144 steers (age = 310 ± 1.5 d; BW = 503 ± 37.2 kg) were individually fed for 105 d. Ruminal digesta samples were collected from eight steers with the greatest (1.96 ± 0.02 kg/d) and eight steers with the least ADG (1.57 ± 0.02 kg/d) that were within ±0.32 SD of the mean (10.1 ± 0.05 kg/d) dry matter. In experiment 2, 66 steers (age = 396 ± 1 d; BW = 456 ± 5 kg) were individually fed for 84 d. Rumen, duodenum, jejunum, ileum, cecum, and colon digesta samples were collected from eight steers with the greatest (2.39 ± 0.06 kg/d) and eight steers with the least ADG (1.85 ± 0.06 kg/d) that were within ±0.55 SD of the mean dry matter intake (11.9 ± 0.1 kg/d). In both studies, DNA was isolated and the V1 to V3 regions of the 16S rRNA gene were sequenced. Operational taxonomic units were classified using 0.03 dissimilarity and identified using the Greengenes 16S rRNA gene database. In experiment 1, there were no differences in the Chao1, Shannon, Simpson, and InvSimpson diversity indexes or the permutation multivariate analysis of variance (PERMANOVA; P = 0.57). The hierarchical test returned six clades as being differentially abundant between steer classifications (P < 0.05). In experiment 2, Chao1, Shannon, Simpson, and InvSimpson diversity indexes and PERMANOVA between steer classified as less or greater ADG did not differ (P > 0.05) for the rumen, duodenum, ileum, cecum, and colon. In the jejunum, there tended to be a difference in the Chao1 (P = 0.09) and Simpson diversity (P = 0.09) indexes between steer classifications, but there was no difference in the Shannon (P = 0.14) and InvSimpson (P = 0.14) diversity indexes. Classification groups for the jejunum differed (P = 0.006) in the PERMANOVA. The hierarchical dependence false discovery rate procedure returned 11 clades as being differentially abundant between steer classifications in the jejunum (P < 0.05). The majority of the OTU were in the Families Corynebacteriaceae and Coriobacteriaceae. This study suggests that intestinal differences in the microbiota of ruminants may be associated with animal performance., (© Published by Oxford University Press on behalf of the American Society of Animal Science 2020.)

Larkspur (Delphinium spp.) poisoning is a long-term problem for cattle grazing on rangelands of western North America. Recent research has shown that both plant and animal-based factors are critical in understanding and mitigating larkspur poisoning in cattle. Non-toxicological factors including sex, age, cattle breed, and plant chemotype affect cattle responses to larkspur. For example, Angus heifers are more susceptible to larkspur intoxication than are steers or bulls. Young cattle appear to be more susceptible to larkspur poisoning than mature animals. Beef breeds of cattle are more susceptible to larkspur intoxication than dairy breeds. In addition to animal factors, plant alkaloid composition (chemotype) affects the potential toxicity for cattle because of differences in the ratios and concentrations of highly toxic N-(methylsuccinimido) anthranoyllycoctonine (MSAL)-type alkaloids compared to less lethal non-MSAL-type alkaloids. Animal- and plant-based factors can provide substantial information to inform livestock producers on management to reduce risk and cattle losses to various larkspur species in western North America., (Published by Elsevier Ltd.)

Disease incidence is intimately associated with an animal's commensal bacteria populations (microbiome), as microbes that are involved with morbidity and mortality are commonly found in animals with no sign of disease. An understanding of the animal's resident respiratory pathogens, in the upper nasal cavity prior to weaning, may help us to understand the impact of these pathogens on incidence of respiratory disease. For this research, the overall goal was to characterize bacterial populations associated with calves at an early age and through time periods prior to weaning in 3 herds at the U.S. Meat Animal Research Center. Nasal swabs from the upper nasal cavity were collected at initial vaccination (approximately 40 d of age), preconditioning (approximately 130 d of age), and weaning (approximately 150 d of age) in 2015 and 2016. DNA was extracted from nasal swabs and combined into 2 pools of 10 animals for each sampling time point, in each herd, for a total of 6 pools at each sampling time point and 18 pools for all sampling time points within each year. To evaluate and compare the microbiome of each pooled sample, hypervariable regions 1 through 3 along the 16S ribosomal RNA (rRNA) gene were amplified by PCR and sequenced using next-generation sequencing (Illumina MiSeq) for identification of the bacterial taxa present. Alpha and beta diversity were also measured. Overall, microbial communities were different between combinations of sampling year, herd location, and sampling time prior to weaning as shown by beta diversity. Analysis of these specific respiratory pathogens prior to weaning will present a clearer picture of the distribution of microbial populations in animals prior to weaning and not exhibiting clinical signs of respiratory disease. Therefore, evaluation of the animal's resident bacterial populations in the upper nasal cavity during different phases of the beef production system may help us to understand the impact of the microbiome on incidence of respiratory disease in cattle., (Published by Oxford University Press on behalf of the American Society of Animal Science 2019.)

Larkspur (Delphinium spp.) poisoning is a long-term problem for cattle grazing on rangelands of western North America. Results from preliminary experiments have suggested that differences in larkspur toxicity may exist between heifers and bulls. The objective of this study was to compare the physiological responses of yearling Angus heifers, steers, and bulls with a standardized dose of Delphinium barbeyi and to test the hypothesis that the response is sex dependent. Clinical signs of intoxication, including muscle coordination and function, were measured 24 h after oral dosing with larkspur by walking the cattle at a pace of 5 to 6 km h-1 for up to 40 min on an oval dirt track. Due to the experimental methods used, the variation in susceptibility to larkspur was not quantifiable for walking times of 0 or 40 min or more. Larkspur susceptible animals that were not able to walk (0 min; 36% of the animals) or larkspur resistant animals that walked the entire test period of 40 min (9% of the animals) resulted in censored or truncated data. The statistical methods (censReg and lmec) were used to adjust for data truncation or censoring. The heifers were only able to walk -8.9 ± 3.9 min (65.5% censored on the left) compared with 13.2 ± 3.7 min for bulls and 15.9 ± 2.7 min for steers. When heifers were compared with bulls and steers together, heifers walked 23.4 ± 4.5 min less (P < 0.0001). Serum alkaloid concentrations were measured immediately before walking, and deltaline concentrations averaged 266 ± 28, 131 ± 20, and 219 ± 28 ng mL-1 for all heifers, steers, and bulls, respectively, and serum methyllycaconitine concentrations averaged 660 ± 46, 397 ± 32, and 612 ± 34 ng mL-1 for all heifers, steers, and bulls, respectively. The relative risk of a zero walk time for yearling heifers is 330% that of yearling bulls (P = 0.0008). These results suggest that yearling Angus heifers are more susceptible to larkspur intoxication and, when possible, heifers should be kept from grazing larkspur-infested rangelands as a simple management tool to reduce the risk of fatal poisoning., (Published by Oxford University Press on behalf of the American Society of Animal Science 2019.)

Objective: The purpose of this study was to evaluate potential relationships between cytokine gene expression, complete blood counts (CBC) and animals that were sick or would become sick. The CBC and the transcript abundance of cytokines and their receptors expressed in leukocytes were measured from calves at two early timepoints, and again after diagnosis with bovine respiratory disease (BRD)., Results: Blood was collected from calves at pre-conditioning (n = 796) and weaning (n = 791) for CBC. Blood counts were also measured for the calves with BRD (n = 13), and asymptomatic calves (n = 75) after weaning. The CBC were compared for these animals at 3 time points. At diagnosis, neutrophils were higher and basophils lower in sick animals (P < 0.05). To further characterize BRD responses, transcript abundance of 84 cytokine genes were evaluated in 5 calves with BRD and 9 asymptomatic animals at all time points. There was more data for CBC than transcript abundance; hence, animal and temporary environmental correlations between CBC and transcript abundance were exploited to improve the power of the transcript abundance data. Expression of CCL16, CXCR1, CCR1 was increased in BRD positive animals compared to controls (P-corrected < 0.1). Cytokine expression data may help to provide insight into an animal's health.

The effect of age on larkspur poisoning of cattle is unknown. An experiment consisting of oral dosing of dried, ground, Delphinium barbeyi to ten Angus steers as yearlings, and again at two years was performed. There was a significant difference between the responses of yearling and two year old steers (P = 0.0015), with yearling steers being more susceptible. These results suggest that the adverse response of Angus cattle to larkspur is age-dependent., (Published by Elsevier Ltd.)

Background: Feed intake and body weight gain are economically important inputs and outputs of beef production systems. The purpose of this study was to discover differentially expressed genes that will be robust for feed intake and gain across a large segment of the cattle industry. Transcriptomic studies often suffer from issues with reproducibility and cross-validation. One way to improve reproducibility is by integrating multiple datasets via meta-analysis. RNA sequencing (RNA-Seq) was performed on longissimus dorsi muscle from 80 steers (5 cohorts, each with 16 animals) selected from the outside fringe of a bivariate gain and feed intake distribution to understand the genes and pathways involved in feed efficiency. In each cohort, 16 steers were selected from one of four gain and feed intake phenotypes (n = 4 per phenotype) in a 2 × 2 factorial arrangement with gain and feed intake as main effect variables. Each cohort was analyzed as a single experiment using a generalized linear model and results from the 5 cohort analyses were combined in a meta-analysis to identify differentially expressed genes (DEG) across the cohorts., Results: A total of 51 genes were differentially expressed for the main effect of gain, 109 genes for the intake main effect, and 11 genes for the gain x intake interaction (P corrected  < 0.05). A jackknife sensitivity analysis showed that, in general, the meta-analysis produced robust DEGs for the two main effects and their interaction. Pathways identified from over-represented genes included mitochondrial energy production and oxidative stress pathways for the main effect of gain due to DEG including GPD1, NDUFA6, UQCRQ, ACTC1, and MGST3. For intake, metabolic pathways including amino acid biosynthesis and degradation were identified, and for the interaction analysis the pathways identified included GADD45, pyridoxal 5'phosphate salvage, and caveolar mediated endocytosis signaling., Conclusions: Variation among DEG identified by cohort suggests that environment and breed may play large roles in the expression of genes associated with feed efficiency in the muscle of beef cattle. Meta-analyses of transcriptome data from groups of animals over multiple cohorts may be necessary to elucidate the genetics contributing these types of biological phenotypes.

Individuals often respond differently to the same vaccine; some of this variation may be caused by genetic differences among animals. Our objective was to estimate heritability and identify genomic regions associated with humoral response to an Escherichia coli O157:H7 vaccine in beef cattle. Crossbred beef cattle (n = 651) were vaccinated with a commercially available E. coli O157:H7 vaccine. Serum was collected at time of initial vaccination (d 0), booster (d 21), and d 56 after initial vaccination. Total antibodies specific to siderophore receptor and porin proteins in the vaccine were quantified by enzyme-linked immunosorbent assay. Genomic DNA was isolated from whole blood and genotyped with the bovine GeneSeek Genomic Profiler-High Density 78K or 26K Single Nucleotide Polymorphism BeadChip and imputed to 777,000 SNP genotypes. Heritability was estimated by restricted maximum likelihood (REML) using both 1) pedigree and 2) genomic relationships among individuals. Fixed effects were contemporary group, calf age, sex, principal components from SNP genotype data, and pedigree-derived heterozygosity effects. Additive and dominance effects of SNPs were estimated individually while accounting for contemporary group, sex, and the top 20 principal components calculated from the genomic relationship matrix. Heritability of initial response to vaccination (d 21 -d 0) was 0.10 ± 0.175 using pedigree relationships and 0.14 ± 0.149 using genomic relationships, but neither estimate was statistically different from zero. Heritability of booster (d 56 -d 21) and overall (d 56 -d 0) responses were low and not statistically significant from zero. There were no clusters of linked SNP associated with vaccine response, but eight regionally isolated SNPs were significantly associated with initial or overall response to vaccination. Regional genetic variation for initial response to an E. coli O157:H7 vaccine was observed, although overall heritability of this response was not statistically significant from zero.

Bovine respiratory disease complex (BRDC) is a multifactor disease, and disease incidence may be associated with an animal's commensal bacterial populations (microbiome) in the upper nasal cavity. Identifying these commensal bacterial populations in the upper nasal cavity may help us to understand the impact of the microbiome on incidence of BRDC in cattle. Various sampling techniques have previously been utilized to evaluate the microbiome of different locations of the upper nasal cavity in cattle. Therefore, our objective was to determine whether bacterial populations of the nasal cavity vary based on these sampling locations. Two common sampling techniques were evaluated, including 6-inch nasal swabs and deep nasopharyngeal swabs. Nasal swabs from calves were collected when the animal was diagnosed with BRDC after weaning in the feedlot in addition to collection of samples from asymptomatic cohorts. Samples were pooled in groups based on year the animal was in the feedlot (2015 or 2016), when the animal was diagnosed with BRDC (1 to 5 weeks after weaning), type of sample (6-inch nasal swab or deep nasopharyngeal swab), and health status (diagnosis with BRDC or control). Variable regions 1 through 3 along the 16S rRNA gene were amplified by PCR and sequenced using next-generation sequencing (Illumina MiSeq) for identification of the bacterial taxa present. Overall, sampling site did not consistently influence diversity of the bacterial populations of the upper nasal cavity. However, the effect of disease incidence on the microbiome was depended on sampling time after weaning (P = 0.0462) for 2015, while the main effects of sampling time after weaning (P = 0.00992) and disease phenotype (P = 0.012) were significant for 2016. These data for 2016 demonstrate that in addition to bacterial profiles changing throughout weaning, calves diagnosed with BRDC have different bacterial profiles compared to their control cohorts. In addition, evaluation of the microbiome identified predominant bacteria genera in the upper nasal cavity included those previously reported to be associated with cattle diagnosed with BRDC including Mycoplasma sp., Psychrobacter sp., and Mannheimia sp. In summary, these results demonstrate that shorter, less invasive 6-inch nasal swabs produce similar results to deep nasopharyngeal swabs.

Four chemically similar alkaloids, anabasine, anabaseine, epibatidine and dimethylphenylpiperazinium (DMPP), are potent nicotinic acetylcholine receptor agonists of fetal muscle nicotinic acetylcholine receptors in human TE-671 cells. Based on results with these cells, we hypothesized that the alkaloids would completely inhibit ultrasound-monitored fetal movement in a goat model. Different, single doses of anabasine, anabaseine, epibatidine, DMPP, or saline control were administered I.V. to pregnant goats on day 40 of gestation and the number of fetal movements per 5 min sample was measured by ultrasound at times 0, 0.5, 1, 2, 4 and 8 h. The differences among does in fetal movements were more consistent at dosing and following recovery for doses of anabasine above 0.125 mg/kg compared to the other compounds and dosages. Anabasine actions were dose-dependent with an IC 50 value of ∼0.1 mg/kg, and, at a dose of 0.8 mg/kg, completely inhibited fetal movement for 1.5 h after dosing. Anabaseine, epibatidine, and DMPP failed to completely inhibit fetal movement in day 40 pregnant goats at doses predicted to be effective. These results suggest that while experiments with TE-671 cells provide valuable information and predictions of the actions of plant alkaloids on fetal movement, in vivo experiments are still required in order to determine the ability of an alkaloid to inhibit fetal movement in livestock species. Moreover, other pharmacological properties such as receptor differences between mammalian species and differences in the pharmacokinetic properties of the alkaloids also are likely to weaken teratologic predictions based solely on the in vitro data., (Published by Elsevier Ltd.)

Mesenteric fat, a depot within the visceral fat, accumulates in cattle during maturation and finishing and may be a potential source of production inefficiency. The aim of this study was to determine whether the genes expressed in the mesenteric fat of steers were associated with body weight gain and feed intake. Sixteen steers chosen by their rank of distance from the bivariate mean for gain and feed intake were used for this study. Mesenteric fat was obtained and evaluated for differences in gene expression. A total of 1831 genes were identified as differentially expressed among steers with variation in feed intake and gain. Many of these genes were involved with metabolic processes such as proteolysis, transcription and translation. In addition, the Gene Ontology annotations including transport and localization were both over-represented among the differentially expressed genes. Pathway analysis was also performed on the differentially expressed genes. The superoxide radical degradation pathway was identified as over-represented based on the differential expression of the genes GPX7, SOD2 and TYRP1, suggesting a potential role for oxidative stress or inflammatory pathways among low gain-high intake animals. GPX7 and SOD2 were in lower transcript abundance, and TYRP1 was higher in transcript abundance among the low gain-high feed intake animals. The retinoate biosynthesis pathway was also enriched due to the differential expression of the genes AKR1C3, ALDH8A1, RDH8, RDH13 and SDR9C7. These genes were all more highly expressed in the low gain-high intake animals. The glycerol degradation and granzyme A signaling pathways were both associated with gain. Three glycerol kinase genes and the GZMA gene were differentially expressed among high vs. low gain animals. Mesenteric fat is a metabolically active tissue, and in this study, genes involved in proteolysis, transcription, translation, transport immune function, glycerol degradation and oxidative stress were differentially expressed among beef steers with variation in body weight gain and feed intake., (© 2017 Stichting International Foundation for Animal Genetics.)

Copy number variations (CNVs) are large insertions, deletions or duplications in the genome that vary between members of a species and are known to affect a wide variety of phenotypic traits. In this study, we identified CNVs in a population of bulls using low coverage next-generation sequence data. First, in order to determine a suitable strategy for CNV detection in our data, we compared the performance of three distinct CNV detection algorithms on benchmark CNV datasets and concluded that using the multiple sample read depth approach was the best method for identifying CNVs in our sequences. Using this technique, we identified a total of 1341 copy number variable regions (CNVRs) from genome sequences of 154 purebred sires used in Cycle VII of the USMARC Germplasm Evaluation Project. These bulls represented the seven most popular beef breeds in the United States: Hereford, Charolais, Angus, Red Angus, Simmental, Gelbvieh and Limousin. The CNVRs covered 6.7% of the bovine genome and spanned 2465 protein-coding genes and many known quantitative trait loci (QTL). Genes harbored in the CNVRs were further analyzed to determine their function as well as to find any breed-specific differences that may shed light on breed differences in adaptation, health and production., (Published 2016. This article is a U.S. Government work and is in the public domain in the USA.)

Background: Feed intake and gain are economically important traits in beef production. The rumen wall interacts with feed, microbial populations, and fermentation products important to cattle nutrition. As such, it is likely to be a critical component in the beef steer's ability to utilize feedstuffs efficiently. To identify genes associated with steer feed intake and body weight gain traits, and to gain an understanding of molecules and pathways involved in feed intake and utilization, RNA sequencing (RNA-Seq) was performed on rumen papillae from 16 steers with variation in gain and feed intake. Four steers were chosen from each of the four Cartesian quadrants for gain×feed intake and used to generate individual RNA-Seq libraries., Results: Normalized read counts from all of the mapped reads from each of the four groups of animals were individually compared to the other three groups. In addition, differentially expressed genes (DEGs) between animals with high and low gain, as well as high and low intake were also evaluated. A total of 931 genes were differentially expressed in the analyses of the individual groups. Eighty-nine genes were differentially expressed between high and low gain animals; and sixty-nine were differentially expressed in high versus low intake animals. Several of the genes identified in this study have been previously associated with feed efficiency. Among those are KLK10, IRX3, COL1A1, CRELD2, HDAC10, IFITM3, and VIM., Conclusions: Many of the genes identified in this study are involved with immune function, inflammation, apoptosis, cell growth/proliferation, nutrient transport, and metabolic pathways and may be important predictors of feed intake and gain in beef cattle., (Published by Elsevier B.V.)

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.

The genome sequence was obtained from 270 sires used in the Germplasm Evaluation (GPE) project. These bulls included 154 purebred AI sires from GPE Cycle VII breeds (Hereford, Angus, Simmental, Limousin, Charolais, Gelbvieh, and Red Angus), 83 F crosses of those breeds, and 33 AI sires from 8 other breeds. The exome capture sequence targeting coding regions of the genome was obtained from 176 of these bulls. Sequence reads were mapped to the UMD 3.1 bovine genome assembly; a mean of 2.5-fold (x) coverage per bull was obtained from the genomic sequence, and the targeted exons were covered at a mean of 20.0x. Over 28.8 million biallelic sequence variants were detected where each allele was present in at least 3 different bulls. These included 22.0 million previously reported variants and 94.1% of the 774,660 autosomal and BTA X SNP on the BovineHD BeadChip assay (HD). More than 92% of the variants detected in targeted exons were also detected from the low-coverage genome sequence. Less than 1% of the variants detected from the combined genome and exome sequence occurred in annotated protein-coding sequences and 5' and 3' untranslated regions (UTR) surrounding the 19,994 annotated protein coding regions. Variation was detected in the coding sequence or UTR of 96.8% of the genes: loss-of-function variants were predicted for 3,298 genes, 14,973 contained nonsynonymous variants, 11,276 had variation in UTR, and 17,721 genes contained synonymous variants. Minor allele frequencies (MAF) were <0.05 for 47.8% of the coding sequence and UTR variants, and MAF distributions were skewed toward low MAF. In contrast, 11.1% of the HD SNP detected in these bulls had MAF < 0.05, and the distribution was skewed toward higher MAF. Genes involved in immune system processes and immune response were overrepresented among those genes containing high MAF loss-of-function and nonsynonymous polymorphisms. Detected variants were submitted to the National Center for Biotechnology Information genetic variation database (dbSNP) under the handle MARC, batch GPE&#95;Bull&#95;GenEx.

Bovine respiratory disease complex (BRDC) is the most expensive disease in beef cattle in the United States costing the industry at least US$1 billion annually. Bovine respiratory disease complex causes damage to lung tissue resulting in persistent lung lesions observable at slaughter. Severe lung lesions at harvest have been associated with decreased preharvest ADG and increased clinical BRDC in the feedlot. Our objective was to identify SNP that are associated with severe lung lesions observed at harvest in feedlot cattle. We conducted a genomewide association study (GWAS) using a case-control design for severe lung lesions in fed cattle at slaughter using the Illumina Bovine HD array (approximately 770,000 SNP) and sample pooling. Lung samples were collected from 11,520 young cattle, a portion of which had not been treated with antibiotics (participating in a "natural" marketing program), at a large, commercial beef processing plant in central Nebraska. Lung samples with lesions (cases) and healthy lungs (controls) were collected when both phenotypes were in close proximity on the viscera (offal) table. We constructed 60 case and 60 control pools with 96 animals per pool. Pools were constructed by sampling sequence to ensure that case and control pool pairs were matched by proximity on the processing line. The Bovine HD array (770,000 SNP) was run on all pools. Fourteen SNP on BTA 2, 3, 4, 9, 11, 14, 15, 22, 24, and 25 were significant at the genomewide experiment-wise error rate of 5% ( ≤ 1.49 × 10). Eighty-five SNP on 28 chromosomes achieved a false discovery rate of 5% ( ≤ 5.38 × 10). Significant SNP were near (±100 kb) genes involved in tissue repair and regeneration, tumor suppression, cell proliferation, apoptosis, control of organ size, and immunity. Based on 85 significantly associated SNP in or near a collection of genes with diverse function on 28 chromosomes, we conclude that the genomic footprint of lung lesions is complex. A complex genomic footprint (genes and regulatory elements that affect the trait) is consistent with what is known about the cause of the disease: complex interactions among multiple viral and bacterial pathogens along with several environmental factors including dust, commingling, transportation, and stress. Characterization of sequence variation near significant SNP will enable accurate and cost effective genome-enhanced genetic evaluations for BRDC resistance in AI bulls and seed stock populations.

Reproductive success is an important component of commercial beef cattle production, and identification of DNA markers with predictive merit for reproductive success would facilitate accurate prediction of mean daughter pregnancy rate, enabling effective selection of bulls to improve female fertility. A previous study identified SNP associated with beef cattle reproductive efficiency based on a genomewide association analysis approach using genotyping multiple-animal pools of DNA to increase the number of animals that could be genotyped with available resources. For the current study, we expand on this previous study by individually genotyping cattle from the pooling study for 89 SNP that were associated with female pregnancy rate. The aims of the study were to confirm the results of the pooling study and, more specifically, identify modes of gene action and DNA variations such as haplotypes that would not be possible with pooled genotyping. Eighty-nine SNP selected from the pooling study were evaluated using the Sequenom MassARRAY system to individually genotype animals from populations evaluated in the pooling study, including both and breeds. From this research, regions on chromosomes 5 (26.3-48.1 Mb; UMD3.1 assembly) and 9 (37,436,575 bp; UMD3.1 assembly), first identified in the previous pooling study, were shown through individual genotyping to harbor genetic variation ( < 0.05 genomewide significance) affecting reproductive efficiency in interspecific crosses ( and ) of cattle. Each of these markers exhibited additive (vs. dominant) gene action. Additionally, a haplotype block harboring an allele of origin with negative effects on reproduction was identified on chromosome 5 in interspecific composite breeds of × composites.

Bovine viral diarrhea viruses (BVDV) comprise a diverse group of viruses that cause disease in cattle. BVDV may establish both transient and persistent infections depending on the developmental stage of the animal at exposure. The objective was to determine whether genomic regions harboring single nucleotide polymorphisms (SNPs) could be associated with the presence or absence of persistent BVDV infection. A genome-wide association approach based on 777,000 SNP markers was used. Samples of animals identified as positive (n = 1200) or negative (n = 1200) for the presence of BVDV in skin samples (n = 1200) were used. DNA samples were combined in 24 pools (100 animals per pool). One SNP, significant at the 5 percent genome-wide level (P = 9.41 × 10(-8) ), was detected on chromosome 14, located at position 80 675 176 bp. Fifteen SNPs, residing on chromosomes 1, 2, 6, 8, 10, 15 and 18, were moderately associated (P < 1 × 10(-5) ) with persistent BVDV infection. Results show that genes harboring or neighboring significant SNPs are involved in leucopenia, signal transduction, RNA splicing and DNA methylation processes., (Published 2014. This article is a U.S. Government work and is in the public domain in the USA.)

Genetic marker effects and type of inheritance are estimated with poor precision when minor marker allele frequencies are low. A stable composite population (MARC III) was subjected to marker-assisted selection for multiple years to equalize specific marker frequencies to 1) estimate effect size and mode of inheritance for previously reported SNP on targeted beef carcass quality traits (n=254), 2) estimate pleiotropic effects of previously reported SNP on nontarget performance traits (n=542 or 254), and 3) evaluate tenderness SNP specific residual variance for LM tenderness. Three haplotypes within μ-calpain (CAPN1), a SNP in calpastatin (CAST), and a dinucleotide substitution in diacylglycerol O-acyltransferase 1 (DGAT1) were successfully selected to equalize their frequencies. Traits evaluated were birth BW, weaning BW, yearling BW, final BW, dressing percent, HCW, fat thickness, LM area, USDA marbling score, yield grade, LM slice shear force (SSF), and visible and near-infrared (VISNIR)-predicted SSF. While the CAPN1 genotype effect on SSF was not significant (P=0.12), the direction and size of CAPN1 contrasts were consistent with previous research. Effects on SSF between divergent CAPN1 haplotypes (1.153 kg) and the additive effect of CAST (0.902 kg) were large, and animals homozygous for tender alleles at both CAPN1 and CAST would have 4.11 kg lower SSF (27.5% of the mean) than animals homozygous tough for both markers. Furthermore, the interaction between CAPN1 and CAST for SSF was not significant (P=0.40). There were significant effects for DGAT1 on adjusted fat thickness (P=0.02) and VISNIR-predicted SSF (P<0.001) with additive and dominance modes of inheritance (P<0.05) for both traits. Furthermore, CAST genotype specific residual variance models fit significantly better (P<0.001) than single residual variance models for SSF, with the tougher genotypes having progressively larger residual (and hence phenotypic) variances. Therefore, risk of a tough steak from the undesired CAST genotype is increased through both an increase in mean and an increase in variation. This work is supportive of the importance of CAPN1 and CAST for mean tenderness in beef, confirms an effect of CAST on beef LM tenderness variation, and identifies an effect of DGAT1 on subcutaneous fat thickness.

Colostrum intake is critical to a piglet's survival and can be measured by precipitating out the γ-immunoglobulins from serum with ammonium sulfate (immunocrit). Genetic analysis of immunocrits on 5312 piglets indicated that the heritabilities (se) for direct and maternal effects were 0.13 (0.06) and 0.53 (0.08) respectively. To identify QTL for direct genetic effects, piglets with the highest and lowest immunocrits from 470 litters were selected. Six sets of DNA pools were created based on sire of the litter. These 12 DNA pools were applied to Illumina Porcine SNP60 BeadChips. Normalized X and Y values were analyzed. Three different SNP selection methods were used: deviation of the mean from high vs. low pools, the deviation adjusted for variance based on binomial theory and ANOVA. The 25 highest ranking SNPs were selected from each evaluation for further study along with 12 regions selected based on a five-SNP window approach. Selected SNPs were individually genotyped in the 988 piglets included in pools as well as in 524 piglets that had intermediate immunocrits. Association analyses were conducted fitting an animal model using the estimated genetic parameters. Nineteen SNPs were nominally associated (P < 0.01) with immunocrit values, of which nine remained significant (P < 0.05) after Bonferroni correction, located in 16 genomic regions on 13 chromosomes. In conclusion, the pooling strategy reduced the cost to scan the genome by more than 80% and identified genomic regions associated with a piglet's ability to acquire γ-immunoglobulin from colostrum. Each method to rank SNPs from the pooled analyses contributed unique validated markers, suggesting that multiple analyses will reveal more QTL than a single analysis., (Published 2014. This article is a U.S. Government work and is in the public domain in the USA.)

Reproductive efficiency is of economic importance in commercial beef cattle production, as failure to achieve pregnancy reduces the number of calves marketed per cow exposed. Identification of genetic markers with predictive merit for reproductive success would facilitate early selection of sires with daughters having improved reproductive rate without increasing generation intervals. To identify regions of the genome harboring variation affecting reproductive success, we applied a genomewide association study (GWAS) approach based on the >700,000 SNP marker assay, using a procedure based on genotyping multianimal pools of DNA to increase the number of animals that could be genotyped with available resources. Cows from several populations were classified according to reproductive efficiency, and DNA was pooled within population and phenotype prior to genotyping. Populations evaluated included a research population at the U.S. Meat Animal Research Center, 2 large commercial ranch populations, and a number of smaller populations (<100 head) across the United States. We detected 2 SNP with significant genomewide association (P ≤ 1.49 × 10(-7)), on BTA21 and BTA29, 3 SNP with suggestive associations (P ≤ 2.91 × 10(-6)) on BTA5, and 1 SNP with suggestive association each on BTA1 and BTA25. In addition to our novel findings, we confirmed previously published associations for SNP on BTA-X and all autosomes except 3 (BTA21, BTA22, and BTA28) encompassing substantial breed diversity including Bos indicus and Bos taurus breeds. The study identified regions of the genome associated with reproductive efficiency, which are being targeted for further analysis to develop robust marker systems, and demonstrated that DNA pooling can be used to substantially reduce the cost of GWAS in cattle.

Reproductive efficiency is arguably the most economically important trait in commercial beef cattle production, as failure to achieve pregnancy reduces the number of calves marketed per cow exposed to breeding. Identification of variation in the genome with predictive merit for reproductive success would facilitate accurate prediction of daughter pregnancy rate in sires enabling effective selection of bulls producing daughters with improved fertility. Copy number variation (CNV) is one form of structural variation that results in abnormal copies of DNA in the genome. Our lab previously reported a region between 25 and 70 Mb on chromosome 5 associated with reproductive efficiency. To further evaluate this region for genomic variations such as CNV, we initially applied a genomewide association approach based on the >700,000 SNP marker assay and PennCNV analysis to 68 Bos indicus cross females from a large commercial ranch in central Florida. A genomic segment located on chromosome 5, spanning the region of approximately 58.37 to 58.44 Mb (Bovine UMD3.1 assembly) was identified as containing a deletion associated with decreased reproductive efficiency. To verify this deletion, a quantitative real-time PCR test was developed to evaluate additional females in the central Florida and U.S. Meat Animal Research Center (USMARC) populations. The frequency of the homozygous deletion was 7% (23/319) in the central Florida population (Bos indicus cross) for females with 2 consecutive failed breeding seasons and 0% in the USMARC population (predominantly Bos taurus) of low reproductive females. To date, we have not identified the deletion in Bos taurus cattle, suggesting that the deletion is specific to Bos indicus cattle. From these data, we have identified a deletion on chromosome 5 associated with reproductive efficiency in Bos indicus-influenced cattle.

To assist cattle producers transition from microsatellite (MS) to single nucleotide polymorphism (SNP) genotyping for parental verification we previously devised an effective and inexpensive method to impute MS alleles from SNP haplotypes. While the reported method was verified with only a limited data set (N = 479) from Brown Swiss, Guernsey, Holstein, and Jersey cattle, some of the MS-SNP haplotype associations were concordant across these phylogenetically diverse breeds. This implied that some haplotypes predate modern breed formation and remain in strong linkage disequilibrium. To expand the utility of MS allele imputation across breeds, MS and SNP data from more than 8000 animals representing 39 breeds (Bos taurus and B. indicus) were used to predict 9410 SNP haplotypes, incorporating an average of 73 SNPs per haplotype, for which alleles from 12 MS markers could be accurately be imputed. Approximately 25% of the MS-SNP haplotypes were present in multiple breeds (N = 2 to 36 breeds). These shared haplotypes allowed for MS imputation in breeds that were not represented in the reference population with only a small increase in Mendelian inheritance inconsistancies. Our reported reference haplotypes can be used for any cattle breed and the reported methods can be applied to any species to aid the transition from MS to SNP genetic markers. While ~91% of the animals with imputed alleles for 12 MS markers had ≤1 Mendelian inheritance conflicts with their parents' reported MS genotypes, this figure was 96% for our reference animals, indicating potential errors in the reported MS genotypes. The workflow we suggest autocorrects for genotyping errors and rare haplotypes, by MS genotyping animals whose imputed MS alleles fail parentage verification, and then incorporating those animals into the reference dataset.

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.

Many traits affecting profitability and sustainability of meat, milk, and fiber production are polygenic, with no single gene having an overwhelming influence on observed variation. No knowledge of the specific genes controlling these traits has been needed to make substantial improvement through selection. Significant gains have been made through phenotypic selection enhanced by pedigree relationships and continually improving statistical methodology. Genomic selection, recently enabled by assays for dense SNP located throughout the genome, promises to increase selection accuracy and accelerate genetic improvement by emphasizing the SNP most strongly correlated to phenotype although the genes and sequence variants affecting phenotype remain largely unknown. These genomic predictions theoretically rely on linkage disequilibrium (LD) between genotyped SNP and unknown functional variants, but familial linkage may increase effectiveness when predicting individuals related to those in the training data. Genomic selection with functional SNP genotypes should be less reliant on LD patterns shared by training and target populations, possibly allowing robust prediction across unrelated populations. Although the specific variants causing polygenic variation may never be known with certainty, a number of tools and resources can be used to identify those most likely to affect phenotype. Associations of dense SNP genotypes with phenotype provide a 1-dimensional approach for identifying genes affecting specific traits; in contrast, associations with multiple traits allow defining networks of genes interacting to affect correlated traits. Such networks are especially compelling when corroborated by existing functional annotation and established molecular pathways. The SNP occurring within network genes, obtained from public databases or derived from genome and transcriptome sequences, may be classified according to expected effects on gene products. As illustrated by functionally informed genomic predictions being more accurate than naive whole-genome predictions of beef tenderness, coupling evidence from livestock genotypes, phenotypes, gene expression, and genomic variants with existing knowledge of gene functions and interactions may provide greater insight into the genes and genomic mechanisms affecting polygenic traits and facilitate functional genomic selection for economically important traits.

Genomic selection involves the assessment of genetic merit through prediction equations that allocate genetic variation with dense marker genotypes. It has the potential to provide accurate breeding values for selection candidates at an early age and facilitate selection for expensive or difficult to measure traits. Accurate across-breed prediction would allow genomic selection to be applied on a larger scale in the beef industry, but the limited availability of large populations for the development of prediction equations has delayed researchers from providing genomic predictions that are accurate across multiple beef breeds. In this study, the accuracy of genomic predictions for 6 growth and carcass traits were derived and evaluated using 2 multibreed beef cattle populations: 3,358 crossbred cattle of the U.S. Meat Animal Research Center Germplasm Evaluation Program (USMARC&#95;GPE) and 1,834 high accuracy bull sires of the 2,000 Bull Project (2000&#95;BULL) representing influential breeds in the U.S. beef cattle industry. The 2000&#95;BULL EPD were deregressed, scaled, and weighted to adjust for between- and within-breed heterogeneous variance before use in training and validation. Molecular breeding values (MBV) trained in each multibreed population and in Angus and Hereford purebred sires of 2000&#95;BULL were derived using the GenSel BayesCπ function (Fernando and Garrick, 2009) and cross-validated. Less than 10% of large effect loci were shared between prediction equations trained on (USMARC&#95;GPE) relative to 2000&#95;BULL although locus effects were moderately to highly correlated for most traits and the traits themselves were highly correlated between populations. Prediction of MBV accuracy was low and variable between populations. For growth traits, MBV accounted for up to 18% of genetic variation in a pooled, multibreed analysis and up to 28% in single breeds. For carcass traits, MBV explained up to 8% of genetic variation in a pooled, multibreed analysis and up to 42% in single breeds. Prediction equations trained in multibreed populations were more accurate for Angus and Hereford subpopulations because those were the breeds most highly represented in the training populations. Accuracies were less for prediction equations trained in a single breed due to the smaller number of records derived from a single breed in the training populations.

Background: WC1 co-receptors belong to the scavenger receptor cysteine-rich (SRCR) superfamily and are encoded by a multi-gene family. Expression of particular WC1 genes defines functional subpopulations of WC1(+) γδ T cells. We have previously identified partial or complete genomic sequences for thirteen different WC1 genes through annotation of the bovine genome Btau&#95;3.1 build. We also identified two WC1 cDNA sequences from other cattle that did not correspond to sequences in the Btau&#95;3.1 build. Their absence in the Btau&#95;3.1 build may have reflected gaps in the genome assembly or polymorphisms among animals. Since the response of γδ T cells to bacterial challenge is determined by WC1 gene expression, it was critical to understand whether individual cattle or breeds differ in the number of WC1 genes or display polymorphisms., Results: Real-time quantitative PCR using DNA from the animal whose genome was sequenced ("Dominette") and sixteen other animals representing ten breeds of cattle, showed that the number of genes coding for WC1 co-receptors is thirteen. The complete coding sequences of those thirteen WC1 genes is presented, including the correction of an error in the WC1-2 gene due to mis-assembly in the Btau&#95;3.1 build. All other cDNA sequences were found to agree with the previous annotation of complete or partial WC1 genes. PCR amplification and sequencing of the most variable N-terminal SRCR domain (domain 1 which has the SRCR "a" pattern) of each of the thirteen WC1 genes showed that the sequences are highly conserved among individuals and breeds. Of 160 sequences of domain 1 from three breeds of cattle, no additional sequences beyond the thirteen described WC1 genes were found. Analysis of the complete WC1 cDNA sequences indicated that the thirteen WC1 genes code for three distinct WC1 molecular forms., Conclusion: The bovine WC1 multi-gene family is composed of thirteen genes coding for three structural forms whose sequences are highly conserved among individual cattle and breeds. The sequence diversity necessary for WC1 genes to function as a multi-genic pattern recognition receptor array is encoded in the genome, rather than generated by recombinatorial diversity or hypermutation.

Reproductive efficiency is of economic importance in commercial beef cattle production, since failure to achieve pregnancy reduces the number of calves marketed. Identification of genetic markers with predictive merit for reproductive success would facilitate early selection of females and avoid inefficiencies associated with sub-fertile cows. To identify regions of the genome harboring variation affecting reproductive success, we applied a genome-wide association approach based on the >700,000 SNP marker assay. To include the largest number of individuals possible under the available budget, cows from several populations were assigned to extremes for reproductive efficiency, and DNA was pooled within population and phenotype before genotyping. Surprisingly, pools prepared from DNA of low reproductive cattle returned fluorescence intensity data intermediate between fertile females and males for SNP mapped to the Y chromosome (i.e., male sex chromosome). The presence of Y-associated material in low reproductive heifers or cows was confirmed by Y-directed PCR, which revealed that 21 to 29% of females in the low reproductive category were positive by a Y chromosome PCR test normally used to sex embryos. The presence of the Y chromosome anomaly was further confirmed with application of additional Y-specific PCR amplicons, indicating the likelihood of the presence of some portion of male sex chromosome in female cattle in various beef cattle herds across the U.S. Discovery of this Y anomaly in low reproductive females may make an important contribution to management of reproductive failures in beef cattle operations.

The promise of genomic selection is accurate prediction of the genetic potential of animals from their genotypes. Simple DNA tests might replace low-accuracy predictions for expensive or lowly heritable measures of puberty and fertility based on performance and pedigree. Knowing with some certainty which DNA variants (e.g., SNP) affect puberty and fertility is the best way to fulfill the promise. Several SNP from the BovineSNP50 assay have tentatively been associated with reproductive traits including age at puberty, antral follicle count, and pregnancy observed on different sets of heifers. However, sample sizes are too small and SNP density is too sparse to definitively determine genomic regions harboring causal variants affecting reproductive success. Additionally, associations between individual SNP and similar phenotypes are inconsistent across data sets, and genomic predictions do not appear to be globally applicable to cattle of different breeds. Discrepancies may be a result of different QTL segregating in the sampled populations, differences in linkage disequilibrium (LD) patterns such that the same SNP are not correlated with the same QTL, and spurious correlations with phenotype. Several approaches can be used independently or in combination to improve detection of genomic factors affecting heifer puberty and fertility. Larger samples and denser SNP will increase power to detect real associations with SNP having more consistent LD with underlying QTL. Meta-analysis combining results from different studies can also be used to effectively increase sample size. High-density genotyping with heifers pooled by pregnancy status or early and late puberty can be a cost-effective means to sample large numbers. Networks of genes, implicated by associations with multiple traits correlated with puberty and fertility, could provide insight into the complex nature of these traits, especially if corroborated by functional annotation, established gene interaction pathways, and transcript expression. Example analyses are provided to demonstrate how integrating information about gene function and regulation with statistical associations from whole-genome SNP genotyping assays might enhance knowledge of genomic mechanisms affecting puberty and fertility, enabling reliable DNA tests to guide heifer selection decisions.

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.

Knowledge of breed composition can be useful in multiple aspects of cattle production, and can be critical for analyzing the results of whole genome-wide association studies currently being conducted around the world. We examine the feasibility and accuracy of using genotype data from the most prevalent bovine genome-wide association studies platform, the Illumina BovineSNP50 array (Illumina Inc., San Diego, CA), to estimate breed composition for individual breeds of cattle. First, allele frequencies (of Illumina-defined allele B) of SNP on the array for each of 16 beef cattle breeds were defined by genotyping a large set of more than 2,000 bulls selected in cooperation with the respective breed associations to be representative of their breed. With these breed-specific allele frequencies, the breed compositions of approximately 2,000 two-, three-, and four-way cross (of 8 breeds) cattle produced at the US Meat Animal Research Center were predicted by using a simple multiple regression technique or Mendel (http://www.genetics.ucla.edu/software/mendel) and their genotypes from the Illumina BovineSNP50 array, and were then compared with pedigree-based estimates of breed composition. The accuracy of marker-based breed composition estimates was 89% when using either estimation method for all breeds except Angus and Red Angus (averaged 79%), based on comparing estimates with pedigree-based average breed composition. Accuracy increased to approximately 88% when these 2 breeds were combined into an aggregate Angus group. Additionally, we used a subset of these markers, approximately 3,000 that populate the Illumina Bovine3K (Illumina Inc.), to see whether breed composition could be estimated with similar accuracy when using this reduced panel of SNP makers. When breed composition was estimated using only SNP in common with the Bovine 3K array, accuracy was slightly reduced to 83%. These results suggest that SNP data from these arrays could be used to estimate breed composition in most US beef cattle in situations where pedigree is not known (e.g., multiple-sire natural service matings, non-source-verified animals in feedlots or at slaughter). This approach can aid analyses that depend on knowledge of breed composition, including identification and adjustment of breed-based population stratification, when performing genome-wide association studies on populations with incomplete pedigrees. In addition, SNP-based breed composition estimates may facilitate fitting cow germplasm to the environment, managing cattle in the feedlot, and tracing disease cases back to the geographic region or farm of origin.

Background: A back curvature defect similar to kyphosis in humans has been observed in swine herds. The defect ranges from mild to severe curvature of the thoracic vertebrate in split carcasses and has an estimated heritability of 0.3. The objective of this study was to identify genomic regions that affect this trait., Results: Single nucleotide polymorphism (SNP) associations performed with 198 SNPs and microsatellite markers in a Duroc-Landrace-Yorkshire resource population (U.S. Meat Animal Research Center, USMARC resource population) of swine provided regions of association with this trait on 15 chromosomes. Positional candidate genes, especially those involved in human skeletal development pathways, were selected for SNP identification. SNPs in 16 candidate genes were genotyped in an F2 population (n = 371) and the USMARC resource herd (n = 1,257) with kyphosis scores. SNPs in KCNN2 on SSC2, RYR1 and PLOD1 on SSC6 and MYST4 on SSC14 were significantly associated with kyphosis in the resource population of swine (P ≤ 0.05). SNPs in CER1 and CDH7 on SSC1, PSMA5 on SSC4, HOXC6 and HOXC8 on SSC5, ADAMTS18 on SSC6 and SOX9 on SSC12 were significantly associated with the kyphosis trait in the F2 population of swine (P ≤ 0.05)., Conclusions: These data suggest that this kyphosis trait may be affected by several loci and that these may differ by population. Carcass value could be improved by effectively removing this undesirable trait from pig populations.

Genomic structural variation is an important and abundant source of genetic and phenotypic variation. Here, we describe the first systematic and genome-wide analysis of copy number variations (CNVs) in modern domesticated cattle using array comparative genomic hybridization (array CGH), quantitative PCR (qPCR), and fluorescent in situ hybridization (FISH). The array CGH panel included 90 animals from 11 Bos taurus, three Bos indicus, and three composite breeds for beef, dairy, or dual purpose. We identified over 200 candidate CNV regions (CNVRs) in total and 177 within known chromosomes, which harbor or are adjacent to gains or losses. These 177 high-confidence CNVRs cover 28.1 megabases or approximately 1.07% of the genome. Over 50% of the CNVRs (89/177) were found in multiple animals or breeds and analysis revealed breed-specific frequency differences and reflected aspects of the known ancestry of these cattle breeds. Selected CNVs were further validated by independent methods using qPCR and FISH. Approximately 67% of the CNVRs (119/177) completely or partially span cattle genes and 61% of the CNVRs (108/177) directly overlap with segmental duplications. The CNVRs span about 400 annotated cattle genes that are significantly enriched for specific biological functions, such as immunity, lactation, reproduction, and rumination. Multiple gene families, including ULBP, have gone through ruminant lineage-specific gene amplification. We detected and confirmed marked differences in their CNV frequencies across diverse breeds, indicating that some cattle CNVs are likely to arise independently in breeds and contribute to breed differences. Our results provide a valuable resource beyond microsatellites and single nucleotide polymorphisms to explore the full dimension of genetic variability for future cattle genomic research.

Chromosomal regions harboring variation affecting cattle birth weight and BW gain to 1 yr of age were identified by marker association using the highly parallel BovineSNP50 BeadChip (50K) assay composed of 54,001 individual SNP. Genotypes were obtained from progeny (F(1); 590 steers) and 2-, 3-, and 4-breed cross grandprogeny (F(1)(2) = F(1) x F(1); 1,306 steers and 707 females) of 150 AI sires representing 7 breeds (22 sires per breed; Angus, Charolais, Gelbvieh, Hereford, Limousin, Red Angus, and Simmental). Genotypes and birth, weaning, and yearling BW records were used in whole-genome association analyses to estimate effects of individual SNP on growth. Traits analyzed included growth component traits: birth weight (BWT), 205-d adjusted birth to weaning BW gain (WG), 160-d adjusted postweaning BW gain (PWG); cumulative traits: 205-d adjusted weaning weight (WW = BWT + WG) and 365-d adjusted yearling weight (YW = BWT + WG + PWG); and indexes of relative differences between postnatal growth and birth weight. Modeled fixed effects included additive effects of calf and dam SNP genotype, year-sex-management contemporary groups, and covariates for calf and dam breed composition and heterosis. Direct and maternal additive polygenic effects and maternal permanent environment effects were random. Missing genotypes, including 50K genotypes of most dams, were approximated with a single-locus BLUP procedure from pedigree relationships and known 50K genotypes. Various association criteria were applied: stringent tests to account for multiple testing but with limited power to detect associations with small effects, and relaxed nominal P that may detect SNP associated with small effects but include excessive false positive associations. Genomic locations of the 231 SNP meeting stringent criteria generally coincided with described previously QTL affecting growth traits. The 12,425 SNP satisfying relaxed tests were located throughout the genome. Most SNP associated with BWT and postnatal growth affected components in the same direction, although detection of SNP associated with one component independent of others presents a possible opportunity for SNP-assisted selection to increase postnatal growth relative to BWT.

Background: A comprehensive transcriptome survey, or gene atlas, provides information essential for a complete understanding of the genomic biology of an organism. We present an atlas of RNA abundance for 92 adult, juvenile and fetal cattle tissues and three cattle cell lines., Results: The Bovine Gene Atlas was generated from 7.2 million unique digital gene expression tag sequences (300.2 million total raw tag sequences), from which 1.59 million unique tag sequences were identified that mapped to the draft bovine genome accounting for 85% of the total raw tag abundance. Filtering these tags yielded 87,764 unique tag sequences that unambiguously mapped to 16,517 annotated protein-coding loci in the draft genome accounting for 45% of the total raw tag abundance. Clustering of tissues based on tag abundance profiles generally confirmed ontology classification based on anatomy. There were 5,429 constitutively expressed loci and 3,445 constitutively expressed unique tag sequences mapping outside annotated gene boundaries that represent a resource for enhancing current gene models. Physical measures such as inferred transcript length or antisense tag abundance identified tissues with atypical transcriptional tag profiles. We report for the first time the tissue-specific variation in the proportion of mitochondrial transcriptional tag abundance., Conclusions: The Bovine Gene Atlas is the deepest and broadest transcriptome survey of any livestock genome to date. Commonalities and variation in sense and antisense transcript tag profiles identified in different tissues facilitate the examination of the relationship between gene expression, tissue, and gene function.

To understand the biology and evolution of ruminants, the cattle genome was sequenced to about sevenfold coverage. The cattle genome contains a minimum of 22,000 genes, with a core set of 14,345 orthologs shared among seven mammalian species of which 1217 are absent or undetected in noneutherian (marsupial or monotreme) genomes. Cattle-specific evolutionary breakpoint regions in chromosomes have a higher density of segmental duplications, enrichment of repetitive elements, and species-specific variations in genes associated with lactation and immune responsiveness. Genes involved in metabolism are generally highly conserved, although five metabolic genes are deleted or extensively diverged from their human orthologs. The cattle genome sequence thus provides a resource for understanding mammalian evolution and accelerating livestock genetic improvement for milk and meat production.

Background: In 2006, an atypical U.S. case of bovine spongiform encephalopathy (BSE) was discovered in Alabama and later reported to be polymorphic for glutamate (E) and lysine (K) codons at position 211 in the bovine prion protein gene (Prnp) coding sequence. A bovine E211K mutation is important because it is analogous to the most common pathogenic mutation in humans (E200K) which causes hereditary Creutzfeldt - Jakob disease, an autosomal dominant form of prion disease. The present report describes a high-throughput matrix-associated laser desorption/ionization-time-of-flight mass spectrometry assay for scoring the Prnp E211K variant and its use to determine an upper limit for the K211 allele frequency in U.S. cattle., Results: The K211 allele was not detected in 6062 cattle, including those from five commercial beef processing plants (3892 carcasses) and 2170 registered cattle from 42 breeds. Multiple nearby polymorphisms in Prnp coding sequence of 1456 diverse purebred cattle (42 breeds) did not interfere with scoring E211 or K211 alleles. Based on these results, the upper bounds for prevalence of the E211K variant was estimated to be extremely low, less than 1 in 2000 cattle (Bayesian analysis based on 95% quantile of the posterior distribution with a uniform prior)., Conclusion: No groups or breeds of U.S. cattle are presently known to harbor the Prnp K211 allele. Because a carrier was not detected, the number of additional atypical BSE cases with K211 will also be vanishingly low.

Background: Classical and atypical bovine spongiform encephalopathies (BSEs) are cattle prion diseases. Distinct bovine prion gene (PRNP) alleles have been associated with classical and atypical BSE susceptibility. However, the full extent of PRNP allele association with BSE susceptibility is not known. A systematic sequence-based genotyping method that detects variation throughout PRNP would be useful for: 1) detecting rare PRNP alleles that may be present in BSE-affected animals and 2) testing PRNP alleles for an association with either classical or atypical BSE susceptibility., Findings: We improved a Sanger-based sequencing strategy for detecting bovine PRNP variation through all exons, introns, and part of the promoter (25.2 kb). Our current method can detect 389 known and other potentially unknown PRNP polymorphisms that may be present in BSE-affected cattle. We determined PRNP genotypes for the first U.S. BSE case and her sire. Previously unknown PRNP polymorphisms were not detected in either animal and all PRNP genotypes support the sire-daughter relationship., Conclusion: The methodologies described here characterize variation throughout PRNP. Consequently, rare PRNP alleles that may be present in BSE-affected cattle can be detected.

As a complement to the Bovine HapMap Consortium project, we initiated a systematic study of the copy numbervariation (CNV) within the same cattle population using array comparative genomic hybridization (array CGH). Oligonucleotide CGH arrays were designed and fabricated to cover all chromosomes with an average interval of 6 kb using the latest bovine genome assembly. In the initial screening, three Holstein bulls were selected to represent major paternal lineages of the Holstein breed with some maternal linkages between these lines. Dual-label hybridizations were performed using either Hereford L1 Dominette 01449 or L1 Domino 99375 as reference. The CNVs were represented by gains and losses of normalized fluorescence intensities relative to the reference. The data presented here, for the first time, demonstrated that significant amounts of germline and fewer somatic CNVs exist in cattle, that many CNVs are common both across diverse cattle breeds and among individuals within a breed, and that array CGH is an effective tool to systematically detect bovine CNV. Selected CNVs have been confirmed by independent methods using real-time (RT) PCR. The strategy used in this study, based on genome higher-orderarchitecture variation, is a powerful approach to generating resources for the identification of novel genomic variation and candidate genes for economically important traits.

Background: Cattle are important agriculturally and relevant as a model organism. Previously described genetic and radiation hybrid (RH) maps of the bovine genome have been used to identify genomic regions and genes affecting specific traits. Application of these maps to identify influential genetic polymorphisms will be enhanced by integration with each other and with bacterial artificial chromosome (BAC) libraries. The BAC libraries and clone maps are essential for the hybrid clone-by-clone/whole-genome shotgun sequencing approach taken by the bovine genome sequencing project., Results: A bovine BAC map was constructed with HindIII restriction digest fragments of 290,797 BAC clones from animals of three different breeds. Comparative mapping of 422,522 BAC end sequences assisted with BAC map ordering and assembly. Genotypes and pedigree from two genetic maps and marker scores from three whole-genome RH panels were consolidated on a 17,254-marker composite map. Sequence similarity allowed integrating the BAC and composite maps with the bovine draft assembly (Btau3.1), establishing a comprehensive resource describing the bovine genome. Agreement between the marker and BAC maps and the draft assembly is high, although discrepancies exist. The composite and BAC maps are more similar than either is to the draft assembly., Conclusion: Further refinement of the maps and greater integration into the genome assembly process may contribute to a high quality assembly. The maps provide resources to associate phenotypic variation with underlying genomic variation, and are crucial resources for understanding the biology underpinning this important ruminant species so closely associated with humans.

The objectives of this study were to 1) estimate the allelic frequencies in US beef cattle of 6 DNA markers reported to be associated with variation in dairy production traits; and 2) evaluate the association of these markers with beef production traits. Several genetic markers have been associated with milk yield or composition, including polymorphisms in secreted phosphoprotein 1 (SPP1; also called osteopontin), growth hormone receptor (GHR), casein S1 (CSN1S1), diacylglycerol O-acyltransferase 1 (DGAT1), peroxisome proliferator-activated receptor gamma co-activator-1alpha (PPARGC1A), and ATP-binding cassette subfamily G (white) member 2 (ABCG2). Allelic frequencies for these 6 markers, and their association with 21 phenotypes, were evaluated in 2 crossbred beef cattle populations that sample influential industry sires. Five of 6 markers were segregating in beef cattle populations; the exception was ABCG2. The SPP1 marker was associated with yearling weight (P = 0.025), live weight at slaughter (P = 0.016), postweaning ADG (P = 0.007), and HCW (P = 0.007) in a large, multisire population representing the 7 most populous beef breeds in the United States. Postweaning growth trait associations were confirmed in an independent population of similar construction, including sires from tropically adapted breeds. The SPP1 marker was associated with yearling weight (P = 0.034), live weight at slaughter (P = 0.011), and postweaning ADG (P = 0.015) and showed a trend toward association with HCW (P = 0.083) in this population. Whereas DGAT1, GHR, and CSN1S1 polymorphisms showed association with some traits in individual populations, the lack of consistent predictive merit between populations indicates they may not be suited for beef cattle selection. No significant associations were observed for the PPARGC1A marker and any of 21 recorded traits, indicating this marker had no apparent value in selection for the beef cattle traits tested in these populations. The SPP1 marker had consistent associations and effect sizes (10.5 to 11.5 kg of live weight at slaughter) in both populations, providing strong evidence for utility of the SPP1 marker for postweaning growth in beef cattle.

Background: Repetitive elements comprise approximately 45% of mammalian genomes and are increasingly known to impact genomic function by contributing to the genomic architecture, by direct regulation of gene expression and by affecting genomic size, diversity and evolution. The ubiquity and increasingly understood importance of repetitive elements contribute to the need to identify and annotate them. We set out to identify previously uncharacterized repetitive DNA in the porcine genome. Once found, we characterized the prevalence of these repeats in other mammals., Results: We discovered 27 repetitive elements in 220 BACs covering 1% of the porcine genome (Comparative Vertebrate Sequencing Initiative; CVSI). These repeats varied in length from 55 to 1059 nucleotides. To estimate copy numbers, we went to an independent source of data, the BAC-end sequences (Wellcome Trust Sanger Institute), covering approximately 15% of the porcine genome. Copy numbers in BAC-ends were less than one hundred for 6 repeat elements, between 100 and 1000 for 16 and between 1,000 and 10,000 for 5. Several of the repeat elements were found in the bovine genome and we have identified two with orthologous sites, indicating that these elements were present in their common ancestor. None of the repeat elements were found in primate, rodent or dog genomes. We were unable to identify any of the replication machinery common to active transposable elements in these newly identified repeats., Conclusion: The presence of both orthologous and non-orthologous sites indicates that some sites existed prior to speciation and some were generated later. The identification of low to moderate copy number repetitive DNA that is specific to artiodactyls will be critical in the assembly of livestock genomes and studies of comparative genomics.