De novo assembly and annotation of the North American bison (Bison bison) reference genome and subsequent variant identification.

Genomic tools have improved the ability to manage bison populations and enhanced efforts to conserve this iconic species. These tools have been particularly useful for detecting introgression of cattle genome within bison herds but are limited by the need to use the cattle genome as a surrogate for mapping reads. This complicates efforts to distinguish the species of origin of chromosomal segments in individual bison at the genomic level. An assembly (Bison_UMD1.0) based on 75X genome coverage by Illumina and 454 reads was generated using the MaSuRCA assembler, generating a 2.81 Gigbases de novo reference genome from American bison. Comparison of bison and domestic cattle references identified 28 443 364 single nucleotide variants and 2 627 645 insertions/deletions distinguishing the species. Sequence alignment of an additional 12 modern bison samples and two historic bison samples to domestic cattle and bison references provides a dataset of genomic variants defining the different species and within-species variation. This first annotated draft assembly represents a resource for the management and conservation of bison, as well as a means to study the effects on the genome of interspecies hybridization. The comparisons of historical bison sequences with the new bison reference identified genomic differences between modern and pre-population bottleneck bison. The results support the application of genomics to enhance future research on disease, the establishment of satellite conservation herds and insight into bison and cattle speciation. The first genome assembly for bison and dataset provides a foundation that can be built upon as genetic technologies improve over the years.

Development of polymorphic markers in the immune gene complex loci of cattle.

The addition of cattle health and immunity traits to genomic selection indices holds promise to increase individual animal longevity and productivity, and decrease economic losses from disease. However, highly variable genomic loci that contain multiple immune-related genes were poorly assembled in the first iterations of the cattle reference genome assembly and underrepresented during the development of most commercial genotyping platforms. As a consequence, there is a paucity of genetic markers within these loci that may track haplotypes related to disease susceptibility. By using hierarchical assembly of bacterial artificial chromosome inserts spanning 3 of these immune-related gene regions, we were able to assemble multiple full-length haplotypes of the major histocompatibility complex, the leukocyte receptor complex, and the natural killer cell complex. Using these new assemblies and the recently released ARS-UCD1.2 reference, we aligned whole-genome shotgun reads from 125 sequenced Holstein bulls to discover candidate variants for genetic marker development. We selected 124 SNPs, using heuristic and statistical models to develop a custom genotyping panel. In a proof-of-principle study, we used this custom panel to genotype 1,797 Holstein cows exposed to bovine tuberculosis (bTB) that were the subject of a previous GWAS study using the Illumina BovineHD array. Although we did not identify any significant association of bTB phenotypes with these new genetic markers, 2 markers exhibited substantial effects on bTB phenotypic prediction. The models and parameters trained in this study serve as a guide for future marker discovery surveys particularly in previously unassembled regions of the cattle genome.

Symposium review: Advances in sequencing technology herald a new frontier in cattle genomics and genome-enabled selection.

The cattle reference genome assembly has underpinned major innovations in beef and dairy genetics through genome-enabled selection, including removal of deleterious recessive variants and selection for favorable alleles affecting quantitative production traits. The initial reference assemblies, up to and including UMD3.1 and Btau4.1, were based on a combination of clone-by-clone sequencing of bacterial artificial chromosome clones generated from blood DNA of a Hereford bull and whole-genome shotgun sequencing of blood DNA from his inbred daughter/granddaughter named L1 Dominette 01449 (Dominette). The approach introduced assembly gaps, misassemblies, and errors, and it limited the ability to assemble regions that undergo rearrangement in blood cells, such as immune gene clusters. Nonetheless, the reference supported the creation of genotyping tools and provided a basis for many studies of gene expression. Recently, long-read sequencing technologies have emerged that facilitated a re-assembly of the reference genome, using lung tissue from Dominette to resolve many of the problems and providing a bridge to place historical studies in common context. The new reference, ARS-UCD1.2, successfully assembled germline immune gene clusters and improved overall continuity (i.e., reduction of gaps and inversions) by over 250-fold. This reference properly places nearly all of the legacy genetic markers used for over a decade in the industry. In this review, we discuss the improvements made to the cattle reference; remaining issues present in the assembly; tools developed to support genome-based studies in beef and dairy cattle; and the emergence of newer genome assembly methods that are producing even higher-quality assemblies for other breeds of cattle at a fraction of the cost. The new frontier for cattle genomics research will likely include a transition from the individual Hereford reference genome, to a "pan-genome" reference, representing all the DNA segments existing in commonly used cattle breeds, bringing the cattle reference into line with the current direction of human genome research.

Improved annotation of the domestic pig genome through integration of Iso-Seq and RNA-seq data.

BACKGROUND: Our understanding of the pig transcriptome is limited. RNA transcript diversity among nine tissues was assessed using poly(A) selected single-molecule long-read isoform sequencing (Iso-seq) and Illumina RNA sequencing (RNA-seq) from a single White cross-bred pig. RESULTS: Across tissues, a total of 67,746 unique transcripts were observed, including 60.5% predicted protein-coding, 36.2% long non-coding RNA and 3.3% nonsense-mediated decay transcripts. On average, 90% of the splice junctions were supported by RNA-seq within tissue. A large proportion (80%) represented novel transcripts, mostly produced by known protein-coding genes (70%), while 17% corresponded to novel genes. On average, four transcripts per known gene (tpg) were identified; an increase over current EBI (1.9 tpg) and NCBI (2.9 tpg) annotations and closer to the number reported in human genome (4.2 tpg). Our new pig genome annotation extended more than 6000 known gene borders (5' end extension, 3' end extension, or both) compared to EBI or NCBI annotations. We validated a large proportion of these extensions by independent pig poly(A) selected 3'-RNA-seq data, or human FANTOM5 Cap Analysis of Gene Expression data. Further, we detected 10,465 novel genes (81% non-coding) not reported in current pig genome annotations. More than 80% of these novel genes had transcripts detected in > 1 tissue. In addition, more than 80% of novel intergenic genes with at least one transcript detected in liver tissue had H3K4me3 or H3K36me3 peaks mapping to their promoter and gene body, respectively, in independent liver chromatin immunoprecipitation data. CONCLUSIONS: These validated results show significant improvement over current pig genome annotations.

Short communication: Identification of the pseudoautosomal region in the Hereford bovine reference genome assembly ARS-UCD1.2.

In cattle, the X chromosome accounts for approximately 3 and 6% of the genome in bulls and cows, respectively. In spite of the large size of this chromosome, very few studies report analysis of the X chromosome in genome-wide association studies and genomic selection. This lack of genetic interrogation is likely due to the complexities of undertaking these studies given the hemizygous state of some, but not all, of the X chromosome in males. The first step in facilitating analysis of this gene-rich chromosome is to accurately identify coordinates for the pseudoautosomal boundary (PAB) to split the chromosome into a region that may be treated as autosomal sequence (pseudoautosomal region) and a region that requires more complex statistical models. With the recent release of ARS-UCD1.2, a more complete and accurate assembly of the cattle genome than was previously available, it is timely to fine map the PAB for the first time. Here we report the use of SNP chip genotypes, short-read sequences, and long-read sequences to fine map the PAB (X chromosome:133,300,518) and simultaneously determine the neighboring regions of reduced homology and true pseudoautosomal region. These results greatly facilitate the inclusion of the X chromosome in genome-wide association studies, genomic selection, and other genetic analysis undertaken on this reference genome.

Analysis of the gut bacterial communities in beef cattle and their association with feed intake, growth, and efficiency.

The impetus behind the global food security challenge is direct, with the necessity to feed almost 10 billion people by 2050. Developing a food-secure world, where people have access to a safe and sustainable food supply, is the principal goal of this challenge. To achieve this end, beef production enterprises must develop methods to produce more pounds of animal protein with less. Selection for feed-efficient beef cattle using genetic improvement technologies has helped to understand and improve the stayability and longevity of such traits within the herd. Yet genetic contributions to feed efficiency have been difficult to identify, and differing genetics, feed regimens, and environments among studies contribute to great variation and interpretation of results. With increasing evidence that hosts and their microbiomes interact in complex associations and networks, examining the gut microbial population variation in feed efficiency may lead to partially clarifying the considerable variation in the efficiency of feed utilization. The use of metagenomics and high-throughput sequencing has greatly impacted the study of the ruminant gut. The ability to interrogate these systems at great depth has permitted a greater understanding of the microbiological and molecular mechanisms involved in ruminant nutrition and health. Although the microbial communities of the reticulorumen have been well documented to date, our understanding of the populations within the gastrointestinal tract as a whole is limited. The composition and phylogenetic diversity of the gut microbial community are critical to the overall well-being of the host and must be determined to fully understand the relationship between the microbiomes within segments of the cattle gastrointestinal tract and feed efficiency, ADG, and ADFI. This review addresses recent research regarding the bacterial communities along the gastrointestinal tract of beef cattle; their association with ADG, ADFI, and feed efficiency; and the potential implications for beef production.

Estimates of epistatic and pleiotropic effects of () and () genetic markers on beef heifer performance traits enhanced by selection.

Genetic marker effects and type of inheritance are estimated with poor precision when minor marker allele frequencies are low. A stable composite population (MARC II) was subjected to marker assisted selection for 2 yr to equalize and genetic marker frequencies to evaluate the epistatic and pleiotropic effects of these markers on BW, reproduction, and first calf performance traits in replacement beef females ( = 171) managed under 2 postweaning development protocols. Traits evaluated on the heifers were birth BW, weaning BW, 11-mo BW, 12-mo BW, 13-mo BW, first breeding season pregnancy evaluation BW, first calving season BW, 11-mo puberty, 12-mo puberty, 13-mo puberty, first breeding season pregnancy, and first calf weaning rate. Additionally, heifer's first calf performance traits of ordinal calving date, first calf birth BW, and first calf weaning BW (with and without age adjustment) were analyzed. Selection to increase minor allele frequencies and balanced sampling across genotype classes enhanced the ability to detect all genetic effects except dominance × dominance epistasis. The × genotype effect was significant ( < 0.05) for 11-mo BW and 12-mo BW and tended to be significant ( = 0.08) for 13-mo BW. Consistently, for all 3 traits, the most significant effect among epistatic × genotype effects was the additive effect, with the G allele decreasing BW. There were no associations between × genotype and fertility related traits ( ≥ 0.46) in this study. Additionally, there were no × genotype associations with first progeny performance traits ( ≥ 0.14). The large effect of the additive × additive interaction on first calf weaning BW was imprecisely estimated, which may warrant further investigation.

Molecular evolution of epizootic hemorrhagic disease viruses in North America based on historical isolates using motif fingerprints.

Epizootic hemorrhagic disease virus (EHDV) is an orbivirus of the Reoviridae family that has significant impact on wild and captive white-tailed deer. Although closely related to bluetongue virus that can cause disease in sheep and cattle, North American EHDV historically has not been associated with disease in cattle or sheep. Severe disease in cattle has been reported with other EHDV strains from East Asia and the Middle East. To understand the potential role of viral genetics in the epidemiology of epizootic hemorrhagic disease, a molecular characterization of North American EHDV strains from 1955 to 2012 was conducted via conventional phylogenetic analysis and a new classification approach using motif fingerprint patterns. Overall, this study indicates that the genetic make-up of EHDV populations in North America have slowly evolved over time. The data also suggested limited reassortment events between serotypes 1 and 2 and introduces a new analysis tool for more detailed sequence pattern analysis.

Microbial community profiles of the jejunum from steers differing in feed efficiency.

Research regarding the association between the microbial community and host feed efficiency in cattle has primarily focused on the rumen. However, the various microbial populations within the gastrointestinal tract as a whole are critical to the overall well-being of the host and need to be examined when determining the interplay between host and nonhost factors affecting feed efficiency. The objective of this study was to characterize the microbial communities of the jejunum among steers differing in feed efficiency. Within 2 contemporary groups of steers, individual ADFI and ADG were determined from animals fed the same diet. At the end of each feeding period, steers were ranked based on their standardized distance from the bivariate mean (ADG and ADFI). Four steers with the greatest deviation within each Cartesian quadrant were sampled ( = 16/group; 2 groups). Bacterial 16S rRNA gene amplicons were sequenced from the jejunum content using next-generation sequencing technology. The phylum Firmicutes accounted for up to 90% of the populations within all samples and was dominated by the families Clostridiaceae and Ruminococcaceae. UniFrac principal coordinate analyses did not indicate any separation of microbial communities within the jejunum based on feed efficiency phenotype, and no significant changes were indicated by bacterial diversity or richness metrics. The relative abundances of microbial populations and operational taxonomic units did reveal significant differences between feed efficiency groups ( < 0.05), including the phylum Proteobacteria ( = 0.030); the families Lachnospiraceae ( = 0.035), Coriobacteriaceae ( = 0.012), and Sphingomonadaceae ( = 0.035); and the genera ( = 0.019), ( = 0.018), and ( = 0.022). The study identified jejunal microbial associations with feed efficiency, ADG, and ADFI. This study suggests the association of the jejunum microbial community as a factor influencing feed efficiency at the 16S level.

Cecum microbial communities from steers differing in feed efficiency.

Apart from the rumen, limited knowledge exists regarding the structure and function of bacterial communities within the gastrointestinal tract and their association with beef cattle feed efficiency. The objective of this study was to characterize the microbial communities of the cecum among steers differing in feed efficiency. Within 2 contemporary groups of steers, individual feed intake and BW gain were determined from animals fed the same diet. Within both of 2 contemporary groups, BW was regressed on feed intake and 4 steers within each Cartesian quadrant were sampled ( = 16/group). Bacterial 16S rRNA gene amplicons were sequenced from the cecal content using next-generation sequencing technology. No significant changes in diversity or richness were detected among quadrants, and UniFrac principal coordinate analysis did not show any differences among quadrants for microbial communities within the cecum. The relative abundances of microbial populations and operational taxonomic units revealed significant differences among feed efficiency groups ( < 0.05). Firmicutes was the dominant cecal phylum in all groups and accounted for up to 81% of the populations among samples. Populations were also dominated by families Ruminococcaceae, Lachnospiraceae, and Clostridiaceae, with significant shifts in the relative abundance of taxa among feed efficiency groups, including families Ruminococcaceae ( = 0.040), Lachnospiraceae ( = 0.020), Erysipelotrichaceae ( = 0.046), and Clostridiaceae ( = 0.043) and genera ( = 0.049), ( = 0.044), ( = 0.042), ( = 0.040), ( = 0.042), and ( = 0.042). The study identified cecal microbial associations with feed efficiency, ADG, and ADFI. This study suggests an association of the cecum microbial community with bovine feed efficiency at the 16S level.

Genomewide association study of lung lesions in cattle using sample pooling.

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.

A polymorphism in myostatin influences puberty but not fertility in beef heifers, whereas µ-calpain affects first calf birth weight.

The use of genetic markers to aid in selection decisions to improve carcass and growth characteristics is of great interest to the beef industry. However, it is important to examine potential antagonistic interactions with fertility in cows before widespread application of marker-assisted selection. The objective of the current experiment was to examine the influence of 2 commercially available markers currently in use for improving carcass traits, the myostatin (MSTN) F94L and µ-calpain (CAPN1) 316 and 4751 polymorphisms, on heifer development and reproductive performance. In Exp. 1, beef heifers (n = 146) were evaluated for growth and reproductive traits over a 3-yr period to determine if these polymorphisms influenced reproductive performance. In Exp. 2, heifers representing the 2 homozygous genotypes for the MSTN F94L polymorphism were slaughtered on d 4 of the estrous cycle and reproductive tracts were collected for morphological examination. In Exp. 1, there was a tendency (P = 0.06) for birth BW to be affected by MSTN with the Leu allele increasing birth BW in an additive fashion. Additionally, MSTN significantly affected the proportion of pubertal heifers by the start of the breeding season (P < 0.05) with the Leu allele additively decreasing the proportion pubertal; however, this did not result in a delay in conception or a decrease in pregnancy rates during the first breeding season (P > 0.15). The GT haplotype of CAPN1, which was previously associated with decreased meat tenderness, was associated with an additive decrease in birth BW of the first calf born to these heifers (P < 0.05). In Exp. 2, there were no differences between the MSTN genotypes for gross or histological morphology of the anterior pituitary, uterus, or ovaries (P > 0.05). From these results, we concluded that the MSTN F94L and CAPN1 polymorphisms can be used to improve carcass traits without compromising fertility in beef heifers. The influence of these markers on cow performance and herd life remains to be determined. While the delay in puberty associated with the MSTN F94L polymorphism did not negatively impact reproductive performance in heifers, caution should be used when combining this marker with other markers for growth or carcass traits until the potential interactions are more clearly understood.

CAPN1, CAST, and DGAT1 genetic effects on preweaning performance, carcass quality traits, and residual variance of tenderness in a beef cattle population selected for haplotype and allele equalization.

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.

Genomewide association study of reproductive efficiency in female cattle.

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.

µ-Calpain, calpastatin, and growth hormone receptor genetic effects on preweaning performance, carcass quality traits, and residual variance of tenderness in Angus cattle selected to increase minor haplotype and allele frequencies.

Genetic marker effects and interactions are estimated with poor precision when minor marker allele frequencies are low. An Angus population was subjected to marker assisted selection for multiple years to increase divergent haplotype and minor marker allele frequencies to 1) estimate effect size and mode of inheritance for previously reported SNP on targeted beef carcass quality traits; 2) estimate effects of previously reported SNP on nontarget performance traits; and 3) evaluate tenderness SNP specific residual variance models compared to a single residual variance model for tenderness. Divergent haplotypes within µ-calpain (CAPN1), and SNP within calpastatin (CAST) and growth hormone receptor (GHR) were successfully selected to increase their frequencies. Traits evaluated were birth BW, weaning BW, final BW, fat thickness, LM area, USDA marbling score, yield grade, slice shear force (SSF), and visible and near infrared predicted slice shear force. Both CAPN1 and CAST exhibited additive (P < 0.001) modes of inheritance for SSF and neither exhibited dominance (P ≥ 0.19). Furthermore, the interaction between CAPN1 and CAST for SSF was not significant (P = 0.55). Estimated additive effects of CAPN1 (1.049 kg) and CAST (1.257 kg) on SSF were large in this study. Animals homozygous for tender alleles at both CAPN1 and CAST would have 4.61 kg lower SSF (38.6% of the mean) than animals homozygous tough for both markers. There was also an effect of CAST on yield grade (P < 0.02). The tender CAST allele was associated with more red meat yield and less trimmable fat. There were no significant effects (P ≥ 0.23) for GHR on any of the traits evaluated in this study. Furthermore, CAST specific residual variance models were found to fit significantly better (P < 0.001) than single residual variance models for SSF, with the tougher genotypes having larger residual variance. Thus, the 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 confirms the importance of CAPN1 and CAST for tenderness in beef, provides a new effect of CAST on beef tenderness, and questions the utility of GHR as a selection marker for beef quality.

Selection for genetic markers in beef cattle reveals complex associations of thyroglobulin and casein1-s1 with carcass and meat traits.

Genetic markers in casein (CSN1S1) and thyroglobulin (TG) genes have previously been associated with fat distribution in cattle. Determining the nature of these genetic associations (additive, recessive, or dominant) has been difficult, because both markers have small minor allele frequencies in most beef cattle populations. This results in few animals homozygous for the minor alleles. selection to increase the frequencies of the minor alleles for 2 SNP markers in these genes was undertaken in a composite population. The objective was to obtain better estimates of genetic effects associated with these markers and determine if there were epistatic interactions. Selection increased the frequencies of minor alleles for both SNP from <0.30 to 0.45. Bulls (n = 24) heterozygous for both SNP were used in 3 yr to produce 204 steer progeny harvested at an average age of 474 d. The combined effect of the 9 CSN1S1 × TG genotypes was associated with carcass-adjusted fat thickness (P < 0.06) and meat tenderness predicted at the abattoir by visible and near-infrared reflectance spectroscopy (P < 0.04). Genotype did not affect BW from birth through harvest, ribeye area, marbling score, slice shear force, or image-based yield grade (P > 0.10). Additive, dominance, and epistatic SNP association effects were estimated from genotypic effects for adjusted fat thickness and predicted meat tenderness. Adjusted fat thickness showed a dominance association with TG SNP (P < 0.06) and an epistatic additive CSN1S1 × additive TG association (P < 0.03). For predicted meat tenderness, heterozygous TG meat was more tender than meat from either homozygote (P < 0.002). Dominance and epistatic associations can result in different SNP allele substitution effects in populations where SNP have the same linkage disequilibrium with causal mutations but have different frequencies. Although the complex associations estimated in this study would contribute little to within-population selection response, they could be important for marker-assisted management or reciprocal selection schemes.

Accuracy of genomic breeding values in multibreed beef cattle populations derived from deregressed breeding values and phenotypes.

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_GPE) and 1,834 high accuracy bull sires of the 2,000 Bull Project (2000_BULL) representing influential breeds in the U.S. beef cattle industry. The 2000_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_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_GPE) relative to 2000_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.

Associations of genetic markers in cattle receiving differing implant protocols.

The potential interaction of growth-promoting implants and genetic markers previously reported to be associated with growth, carcass traits, and tenderness was evaluated. Two implant protocols were applied to subsets of steers (n = 383) and heifers (n = 65) that were also genotyped for 47 SNP reported to be associated with variation in growth, fat thickness, LM area, marbling, or tenderness. The "mild" protocol consisted of a single terminal implant [16 mg estradiol benzoate (EB), 80 mg trenbalone acetate (TBA) or 8 mg EB, 80 mg TBA given to steers and heifers, respectively]. The "aggressive" protocol consisted of both a growing implant (8 mg EB, 40 mg TBA) for the lightest half of the animals on the aggressive protocol and 2 successive implants (28 mg EB, 200 mg TBA) given to all animals assigned to the aggressive treatment. Implant protocol had measurable impact on BW and ADG (P < 0.05), with the aggressive protocol increasing these traits before the terminal implant (relative to the mild protocol), whereas the mild protocol increased ADG after the terminal implant so that the final BW and ADG over the experimental period were similar between protocols. Animals on the aggressive protocol had significantly increased (P < 0.05) LM area (1.9 cm(2)), slice shear force (1.4 kg), and intact desmin (0.05 units), but decreased (P < 0.05) marbling score (49 units) and adjusted fat thickness (0.1 cm), and yield grade (0.15 units). Among both treatments, 8 of 9 growth-related SNP were associated with BW or ADG, and 6 of 17 tenderness-related SNP were associated with slice shear force or intact desmin. Favorable growth alleles generally were associated with increased carcass yield traits but decreased tenderness. Similarly, favorable tenderness genotypes for some markers were associated with decreased BW and ADG. Some interactions of implant protocol and genotype were noted, with some growth SNP alleles increasing the effect of the aggressive protocol. In contrast, putative beneficial effects of favorable tenderness SNP alleles were mitigated by the effects of aggressive implant. These type of antagonisms of management variables and genotypes must be accounted for in marker assisted selection (MAS) programs, and our results suggest that MAS could be used to manage, but likely will not eliminate negative impact of implants on quality.

Comparison of bacterial communities in faeces of beef cattle fed diets containing corn and wet distillers' grain with solubles.

AIM: The mammalian intestinal microflora has been shown to impact host physiology. In cattle, intestinal bacteria are also associated with faecal contamination of environmental sources and human illness via foodborne pathogens. Use of wet distillers' grains with solubles (WDGS) in cattle feed creates a gastrointestinal environment where some bacterial species are enriched. Here, we examine if a diet containing 40% WDGS results in fundamentally different microbial community structures. METHODS AND RESULTS: The 20,002 16S r-RNA gene sequences from 20 beef cattle were analysed using Sanger sequencing methods. At the genus level, Prevotella (Gram negative) and Anaerobacter (Gram positive) were the most frequently occurring bacteria in our beef cattle faecal samples. Diet-associated differences in prevalence were noted for Prevotella but not Anaerobacter. CONCLUSIONS: Diet affects community structure. Faecal communities of co-housed beef cattle are not identical. SIGNIFICANCE AND IMPACT OF THE STUDY: It is known that a diet of 40% corn-based WDGS increases the generic Escherichia coli in the faeces and enriches E. coli O157:H7. The results from the current study suggest that in addition to previously observed changes in E. coli, the entire bacterial community structure is different for animals fed 40% corn-based WDGS compared to a traditional corn-finishing diet.

Genetic markers on BTA14 predictive for residual feed intake in beef steers and their effects on carcass and meat quality traits.

With the high cost of feed for animal production, genetic selection for animals that metabolize feed more efficiently could result in substantial cost savings for cattle producers. The purpose of this study was to identify DNA markers predictive for differences among cattle for traits associated with feed efficiency. Crossbred steers were fed a high-corn diet for 140 days and average daily feed intake (ADFI), average daily gain (ADG), and residual feed intake (RFI) phenotypes were obtained. A region on chromosome 14 was previously associated with RFI in this population of animals. To develop markers with the highest utility for predicting an animal's genetic potential for RFI, we genotyped additional markers within this chromosomal region. These polymorphisms were genotyped on the same animals (n = 1066) and tested for association with ADFI, ADG and RFI. Six markers within this region were associated with RFI (P ≤ 0.05). After conservative correction for multiple testing, one marker at 25.09 Mb remained significant (P = 0.02) and is responsible for 3.6% of the RFI phenotypic variation in this population of animals. Several of these markers were also significant for ADG, although none were significant after correction. Marker alleles with positive effects on ADG corresponded to lower RFI, suggesting an effect increasing growth without increasing feed intake. All markers were also assessed for their effects on meat quality and carcass traits. All of the markers associated with RFI were associated with adjusted fat thickness (AFT, P ≤ 0.009) and three were also associated with hot carcass weight (HCW, P ≤ 0.003). Marker alleles associated with lower RFI were also associated with reduced AFT, and if they were associated for HCW, the effect was an increase in weight. These markers may be useful as prediction tools for animals that utilize feed more efficiently; however, validation with additional populations of cattle is required.