Roman L. Hruska U.S. Meat Animal Research Center

For 10,000 years pigs and humans have shared a close and complex relationship. From domestication to modern breeding practices, humans have shaped the genomes of domestic pigs. Here we present the assembly and analysis of the genome sequence of a female domestic Duroc pig (Sus scrofa) and a comparison with the genomes of wild and domestic pigs from Europe and Asia. Wild pigs emerged in South East Asia and subsequently spread across Eurasia. Our results reveal a deep phylogenetic split between European and Asian wild boars ∼1 million years ago, and a selective sweep analysis indicates selection on genes involved in RNA processing and regulation. Genes associated with immune response and olfaction exhibit fast evolution. Pigs have the largest repertoire of functional olfactory receptor genes, reflecting the importance of smell in this scavenging animal. The pig genome sequence provides an important resource for further improvements of this important livestock species, and our identification of many putative disease-causing variants extends the potential of the pig as a biomedical model. This study presents the assembly and analysis of the genome sequence of a female domestic Duroc pig and a comparison with the genomes of wild and domestic pigs from Europe and Asia; the results shed light on the evolutionary relationship between European and Asian wild boars. The domestic pig (Sus scrofa) is an important livestock species, its genome shaped by thousands of years of domestication and, latterly, sophisticated breeding practices. A high-quality draft genome sequence for a female domestic Duroc pig is published in this issue of Nature, under the auspices of the Swine Genome Sequencing Consortium. Comparisons of the genomes of wild and domestic pigs shed light on the evolutionary relationship between European and Asian wild boars, and reveal the rapid evolution of genes involved in the immune response and in olfaction. The authors identify many possible disease-causing gene variants, increasing the potential of the pig as a biomedical model, and present a detailed analysis of endogenous porcine retroviruses, knowledge of which is important for the possible use of pigs in xenotransplantation.

A visibly distinct muscular hypertrophy (mh), commonly known as double muscling, occurs with high frequency in the Belgian Blue and Piedmontese cattle breeds. The autosomal recessive mh locus causing double-muscling condition in these cattle maps to bovine chromosome 2 within the same interval as myostatin, a member of the TGF-beta superfamily of genes. Because targeted disruption of myostatin in mice results in a muscular phenotype very similar to that seen in double-muscled cattle, we have evaluated this gene as a candidate gene for double-muscling condition by cloning the bovine myostatin cDNA and examining the expression pattern and sequence of the gene in normal and double-muscled cattle. The analysis demonstrates that the levels and timing of expression do not appear to differ between Belgian Blue and normal animals, as both classes show expression initiating during fetal development and being maintained in adult muscle. Moreover, sequence analysis reveals mutations in heavy-muscled cattle of both breeds. Belgian Blue cattle are homozygous for an 11-bp deletion in the coding region that is not detected in cDNA of any normal animals examined. This deletion results in a frame-shift mutation that removes the portion of the Myostatin protein that is most highly conserved among TGF-beta family members and that is the portion targeted for disruption in the mouse study. Piedmontese animals tested have a G-A transition in the same region that changes a cysteine residue to a tyrosine. This mutation alters one of the residues that are hallmarks of the TGF-beta family and are highly conserved during evolution and among members of the gene family. It therefore appears likely that the mh allele in these breeds involves mutation within the myostatin gene and that myostatin is a negative regulator of muscle growth in cattle as well as mice.

The success of genome-wide association (GWA) studies for the detection of sequence variation affecting complex traits in human has spurred interest in the use of large-scale high-density single nucleotide polymorphism (SNP) genotyping for the identification of quantitative trait loci (QTL) and for marker-assisted selection in model and agricultural species. A costeffective and efficient approach for the development of a custom genotyping assay interrogating 54,001 SNP loci to support GWA applications in cattle is described. A novel algorithm for achieving a compressed inter-marker interval distribution proved remarkably successful, with median interval of 37 kb and maximum predicted gap of ,350 kb. The assay was tested on a panel of 576 animals from 21 cattle breeds and six outgroup species and revealed that from 39,765 to 46,492 SNP are polymorphic within individual breeds (average minor allele frequency (MAF) ranging from 0.24 to 0.27). The assay also identified 79 putative copy number variants in cattle. Utility for GWA was demonstrated by localizing known variation for coat color and the presence/absence of horns to their correct genomic locations. The combination of SNP selection and the novel spacing algorithm allows an efficient approach for the development of high-density genotyping platforms in species having full or even moderate quality draft sequence. Aspects of the approach can be exploited in species which lack an available genome sequence. The BovineSNP50 assay described here is commercially available from Illumina and provides a robust platform for mapping disease genes and QTL in cattle.

Data from 3 summer feedlot studies were utilized to determine the environmental factors that influence heat stress in cattle and also to determine wind speed (WSPD; m.s(-1)) and solar radiation (RAD; W.m(-2)) adjustments to the temperature-humidity index (THI). Visual assessments of heat stress, based on panting scores (0 = no panting to 4 = severe panting), were collected from 1400 to 1700. Mean daily WSPD, black globe temperature at 1500, and minimums for nighttime WSPD, nighttime black globe THI, and daily relative humidity were found to have the greatest influence on panting score from 1400 to 1700 (R2 = 0.61). From hourly values for THI, WSPD, and RAD, panting score was determined to equal -7.563 + (0.121 x THI) - (0.241 x WSPD) + (0.00082 x RAD) (R2 = 0.49). Using the ratio of WSPD to THI and RAD to THI (- 1.992 and 0.0068 for WSPD and RAD, respectively), adjustments to the THI were derived for WSPD and RAD. On the basis of these ratios and the average hourly data for 1400 to 1700, the THI, adjusted for WSPD and RAD, equals [4.51 + THI - (1.992 x WSPD) + (0.0068 x RAD)]. Four separate cattle studies, comparable in size, type of cattle, and number of observations to the 3 original studies, were utilized to evaluate the accuracy of the THI equation adjusted for WSPD and RAD, and the relationship between the adjusted THI and panting score. Mean panting score derived from individual observations of black-hided cattle in these 4 studies were 1.22, 0.94, 1.32, and 2.00 vs. the predicted panting scores of 1.15, 1.17, 1.30, and 1.96, respectively. Correlations between THI and panting score in these studies ranged from r = 0.47 to 0.87. Correlations between the adjusted THI and mean panting score ranged from r = 0.64 to 0.80. These adjustments would be most appropriate to use, within a day, to predict THI during the afternoon hours using hourly data or current conditions. In addition to afternoon conditions, nighttime conditions, including minimum WSPD, minimum black globe THI, and minimum THI, were also found to influence heat stress experienced by cattle. Although knowledge of THI alone is beneficial in determining the potential for heat stress, WSPD and RAD adjustments to the THI more accurately assess animal discomfort.

We report the most extensive physically anchored linkage map for cattle produced to date. Three-hundred thirteen genetic markers ordered in 30 linkage groups, anchored to 24 autosomal chromosomes (n = 29), the X and Y chromosomes, four unanchored syntenic groups and two unassigned linkage groups spanning 2464 cM of the bovine genome are summarized. The map also assigns 19 type I loci to specific chromosomes and/or syntenic groups and four cosmid clones containing informative microsatellites to chromosomes 13, 25 and 29 anchoring syntenic groups U11, U7 and U8, respectively. This map provides the skeletal framework prerequisite to development of a comprehensive genetic map for cattle and analysis of economic trait loci (ETL).

BACKGROUND: Major advances in selection progress for cattle have been made following the introduction of genomic tools over the past 10-12 years. These tools depend upon the Bos taurus reference genome (UMD3.1.1), which was created using now-outdated technologies and is hindered by a variety of deficiencies and inaccuracies. RESULTS: We present the new reference genome for cattle, ARS-UCD1.2, based on the same animal as the original to facilitate transfer and interpretation of results obtained from the earlier version, but applying a combination of modern technologies in a de novo assembly to increase continuity, accuracy, and completeness. The assembly includes 2.7 Gb and is >250× more continuous than the original assembly, with contig N50 >25 Mb and L50 of 32. We also greatly expanded supporting RNA-based data for annotation that identifies 30,396 total genes (21,039 protein coding). The new reference assembly is accessible in annotated form for public use. CONCLUSIONS: We demonstrate that improved continuity of assembled sequence warrants the adoption of ARS-UCD1.2 as the new cattle reference genome and that increased assembly accuracy will benefit future research on this species.

A survey was performed to estimate the frequency of enterohemorrhagic Escherichia coli O157:H7 or O157:nonmotile (EHEC O157) in feces and on hides within groups of fed cattle from single sources (lots) presented for slaughter at meat processing plants in the Midwestern United States, as well as frequency of carcass contamination during processing from cattle within the same lots. Of 29 lots sampled, 72% had at least one EHEC O157-positive fecal sample and 38% had positive hide samples. Overall, EHEC O157 prevalence in feces and on hides was 28% (91 of 327) and 11% (38 of 355), respectively. Carcass samples were taken at three points during processing: preevisceration, postevisceration before antimicrobial intervention, and postprocessing after carcasses entered the cooler. Of 30 lots sampled, 87% had at least one EHEC O157-positive preevisceration sample, 57% of lots were positive postevisceration, and 17% had positive postprocessing samples. Prevalence of EHEC O157 in the three postprocessing samples was 43% (148 of 341), 18% (59 of 332) and 2% (6 of 330), respectively. Reduction in carcass prevalence from preevisceration to postprocessing suggests that sanitary procedures were effective within the processing plants. Fecal and hide prevalence were significantly correlated with carcass contamination (P = 0.001), indicating a role for control of EHEC O157 in live cattle.

BACKGROUND: The dissection of complex traits of economic importance to the pig industry requires the availability of a significant number of genetic markers, such as single nucleotide polymorphisms (SNPs). This study was conducted to discover several hundreds of thousands of porcine SNPs using next generation sequencing technologies and use these SNPs, as well as others from different public sources, to design a high-density SNP genotyping assay. METHODOLOGY/PRINCIPAL FINDINGS: A total of 19 reduced representation libraries derived from four swine breeds (Duroc, Landrace, Large White, Pietrain) and a Wild Boar population and three restriction enzymes (AluI, HaeIII and MspI) were sequenced using Illumina's Genome Analyzer (GA). The SNP discovery effort resulted in the de novo identification of over 372K SNPs. More than 549K SNPs were used to design the Illumina Porcine 60K+SNP iSelect Beadchip, now commercially available as the PorcineSNP60. A total of 64,232 SNPs were included on the Beadchip. Results from genotyping the 158 individuals used for sequencing showed a high overall SNP call rate (97.5%). Of the 62,621 loci that could be reliably scored, 58,994 were polymorphic yielding a SNP conversion success rate of 94%. The average minor allele frequency (MAF) for all scorable SNPs was 0.274. CONCLUSIONS/SIGNIFICANCE: Overall, the results of this study indicate the utility of using next generation sequencing technologies to identify large numbers of reliable SNPs. In addition, the validation of the PorcineSNP60 Beadchip demonstrated that the assay is an excellent tool that will likely be used in a variety of future studies in pigs.

Adam Phillippy, Curtis Van Tassell, Timothy Smith and colleagues present a new reference genome assembly for the domestic goat using a pipeline that improves contiguity of the assembly by more than 250-fold. The pipeline uses a combination of short- and long-read sequencing, optical mapping, and chromatin interaction mapping. The decrease in sequencing cost and increased sophistication of assembly algorithms for short-read platforms has resulted in a sharp increase in the number of species with genome assemblies. However, these assemblies are highly fragmented, with many gaps, ambiguities, and errors, impeding downstream applications. We demonstrate current state of the art for de novo assembly using the domestic goat (Capra hircus) based on long reads for contig formation, short reads for consensus validation, and scaffolding by optical and chromatin interaction mapping. These combined technologies produced what is, to our knowledge, the most continuous de novo mammalian assembly to date, with chromosome-length scaffolds and only 649 gaps. Our assembly represents a ∼400-fold improvement in continuity due to properly assembled gaps, compared to the previously published C. hircus assembly, and better resolves repetitive structures longer than 1 kb, representing the largest repeat family and immune gene complex yet produced for an individual of a ruminant species.

The focal point of this limited review is bioenergetic research conducted in the Biological Engineering Research Unit at the U.S. Meat Animal Research Center (MARC), using recently developed instrumentation and analytical techniques. The dynamics of observed thermoregulatory responses in cattle to thermal heat load challenges are explored, with an emphasis on physiological and behavioral parameters of body temperature, respiration rate, and feed intake. Observations of body temperature, especially tympanic temperature, have shown hot environments to cause phase shifts, increased amplitude, and increased means for diurnal rhythms. Fractal analysis of body temperature records obtained at 2- to 10-min intervals has been found to be robust for objectively differentiating among responses of cattle in cool to hot environments, and it indicates a stress threshold of approximately 25 degrees C (coincident with declining feed intake). Other analyses determined a 21 degrees C threshold for increased respiration rate. The reported observations and analyses provide further understanding of how and why the animals respond to environmental challenges, an understanding that is necessary for refining performance models and developing energetic and thermoregulatory models. The dynamic responses are discussed in the context of establishing criteria for proactive environmental management for cattle during hot weather, using heat waves as an example.

BACKGROUND: Centromeres are essential for chromosome segregation, yet their DNA sequences evolve rapidly. In most animals and plants that have been studied, centromeres contain megabase-scale arrays of tandem repeats. Despite their importance, very little is known about the degree to which centromere tandem repeats share common properties between different species across different phyla. We used bioinformatic methods to identify high-copy tandem repeats from 282 species using publicly available genomic sequence and our own data. RESULTS: Our methods are compatible with all current sequencing technologies. Long Pacific Biosciences sequence reads allowed us to find tandem repeat monomers up to 1,419 bp. We assumed that the most abundant tandem repeat is the centromere DNA, which was true for most species whose centromeres have been previously characterized, suggesting this is a general property of genomes. High-copy centromere tandem repeats were found in almost all animal and plant genomes, but repeat monomers were highly variable in sequence composition and length. Furthermore, phylogenetic analysis of sequence homology showed little evidence of sequence conservation beyond approximately 50 million years of divergence. We find that despite an overall lack of sequence conservation, centromere tandem repeats from diverse species showed similar modes of evolution. CONCLUSIONS: While centromere position in most eukaryotes is epigenetically determined, our results indicate that tandem repeats are highly prevalent at centromeres of both animal and plant genomes. This suggests a functional role for such repeats, perhaps in promoting concerted evolution of centromere DNA across chromosomes.

The domestication and development of cattle has considerably impacted human societies, but the histories of cattle breeds and populations have been poorly understood especially for African, Asian, and American breeds. Using genotypes from 43,043 autosomal single nucleotide polymorphism markers scored in 1,543 animals, we evaluate the population structure of 134 domesticated bovid breeds. Regardless of the analytical method or sample subset, the three major groups of Asian indicine, Eurasian taurine, and African taurine were consistently observed. Patterns of geographic dispersal resulting from co-migration with humans and exportation are recognizable in phylogenetic networks. All analytical methods reveal patterns of hybridization which occurred after divergence. Using 19 breeds, we map the cline of indicine introgression into Africa. We infer that African taurine possess a large portion of wild African auroch ancestry, causing their divergence from Eurasian taurine. We detect exportation patterns in Asia and identify a cline of Eurasian taurine/indicine hybridization in Asia. We also identify the influence of species other than Bos taurus taurus and B. t. indicus in the formation of Asian breeds. We detect the pronounced influence of Shorthorn cattle in the formation of European breeds. Iberian and Italian cattle possess introgression from African taurine. American Criollo cattle originate from Iberia, and not directly from Africa with African ancestry inherited via Iberian ancestors. Indicine introgression into American cattle occurred in the Americas, and not Europe. We argue that cattle migration, movement and trading followed by admixture have been important forces in shaping modern bovine genomic variation.

We report a bovine linkage map constructed with 1236 polymorphic DNA markers and 14 erythrocyte antigens and serum proteins. The 2990-cM map consists of a sex-specific, X chromosome linkage group and 29 sex-averaged, autosomal linkage groups with an average interval size of 2.5 cM. The map contains 627 new markers and 623 previously linked markers, providing a basis for integrating the four published bovine maps. Orientation and chromosomal assignment of all the linkage groups, except BTA20 and BTA22, was provided by 88 markers that were assigned previously to chromosomes. This map provides sufficient marker density for genomic scans of populations segregating quantitative trait loci (QTL) and subsequent implementation of marker-assisted selection (MAS) mating schemes.

We report the highest density genetic linkage map for a livestock species produced to date. Three published maps for Sus scrofa were merged by genotyping virtually every publicly available microsatellite across a single reference population to yield 1042 linked loci, 536 of which are novel assignments, spanning 2286.2 cM (average interval 2.23 cM) in 19 linkage groups (18 autosomal and X chromosomes, n = 19). Linkage groups were constructed de novo and mapped by locus content to avoid propagation of errors in older genotypes. The physical and genetic maps were integrated with 123 informative loci assigned previously by fluorescence in situ hybridization (FISH). Fourteen linkage groups span the entire length of each chromosome. Coverage of chromosomes 11, 12, 15, and 18 will be evaluated as more markers are physically assigned. Marker-deficient regions were identified only on 11q1.7-qter and 14 cen-q1.2. Recombination rates (cM/Mbp) varied between and within chromosomes. Short chromosomal arms recombined at higher rates than long arms, and recombination was more frequent in telomeric regions than in pericentric regions. The high-resolution comprehensive map has the marker density needed to identify quantitative trait loci (QTL), implement marker-assisted selection or introgression and YAC contig construction or chromosomal microdissection.

Animal production results in conversion of feeds into valuable products such as meat, milk, eggs, and wool as well as into unavoidable and less desirable waste products. Intensification of animal numbers and increasing urbanization has resulted in considerable attention to odorous gases produced from animal wastes. It is clear that animal manure was, and still is, a valuable resource. However, it may be a major obstacle to future development of the animal industry if its impact on the environment is not properly controlled. Poor odor prevention and control from animal wastes is related to a lack of knowledge of the fundamental nature of odor and its production by farm animals. Odor, like noise, is a nuisance or disturbance and there is no universally accepted definition of an objectionable odor. Thus, regulation and control of odors in the environment is difficult because of the technical difficulties of defining odor limits and their measurement and evaluation. A variety of direct (sensory) and indirect (analytical instruments) methods for measuring odor intensity and determination of individual or key odor components are discussed. The biological origins of the four principal classes of odor compounds, namely branched- and straight-chain VFA, ammonia and volatile amines, indoles and phenols, and the volatile sulfur-containing compounds, are reviewed. Because more than 50% of N from animals is excreted as urea, one strategy to conserve N in waste is to inhibit the urease enzyme that converts urea to ammonia. Laboratory studies to evaluate di- and triamide compounds to control urea hydrolysis in slurries of cattle and swine wastes are presented. Finally, a brief overview of various intervention strategies is provided. Multiple combinations of nutritional management, housing systems, treatment options as well as storage and disposal of animal wastes will be required to reduce environmental pollution and provide for long-term sustainable growth.

Cell surface charge and hydrophobicity of Bacillus subtilis, Escherichia coli O157:H7, Listeria monocytogenes, Salmonella typhimurium, Serratia marcescens, Staphylococcus aureus, and Staphylococcus epidermidis were determined by hydrocarbon adherence, hydrophobic interaction, and electrostatic interaction chromatography. Surface charge and hydrophobicity were compared with the initial attachment values and rates of attachment of the bacteria to meat surfaces. There was a linear correlation between the relative negative charge on the bacterial cell surface and initial attachment to lean beef muscle (r2 = 0.885) and fat tissue (r2 = 0.777). Hydrophobicity correlated well with attachment to fat tissue only. The relative hydrophobicity of each bacterium was dependent on the specific method of determination, with wide variations noted between methods.

The objective of this study was to characterize genetic, environmental, and economic factors related to the incidence of bovine respiratory disease (BRD) in feedlot calves. Records from 18,112 calves representing 9 breeds (Angus, Braunvieh, Charolais, Gelbvieh, Hereford, Limousin, Pinzgauer, Red Poll, and Simmental) and 3 composite types (MARC I, MARC II, and MARC III) over a 15-yr period (1987 to 2001) were evaluated. Disease incidence was observed and recorded by station veterinary and technical staff. The incidence of BRD varied across years, with the annual observed incidence ranging from 5 to 44%. From 1987 to 1992, the annual average incidence generally exceeded 20%. However, in later years the annual incidence did not exceed 14%. The epidemiological pattern indicated that BRD infection increased dramatically after 5 d on feed and remained high until approximately 80 d on feed. Previous BRD infection during the preweaning period did not influence subsequent BRD infection in the feedlot. Steers were more likely to become sick with BRD than heifers; castration before entry in the feedlot may be a predisposing cause. Few significant differences among breeds were detected for BRD incidence. Adjusted solutions from mixed model analyses indicated that Herefords were generally more susceptible to BRD infection (P < 0.05) than MARC I and III composite types. Composite breed types had similar susceptibility compared with other purebred breeds. Mortality associated with BRD was greatest in Red Poll calves (9%) compared with the average over all breeds (4%). Estimates of heritability for resistance to BRD ranged from 0.04 to 0.08 +/- 0.01. When the observed heritability was transformed to an underlying continuous scale, the estimate increased to 0.18. Selection for resistance to BRD could be effective if phenotypes for BRD resistance were known. Thus, development of an inexpensive and humane method of challenging animals with BRD to determine resistance would be an important step in reducing the incidence of BRD. This study also demonstrated that producer-collected field data could be used for selection against this disease. The economic loss associated with lower gains and treatment costs for BRD infection in a 1,000-cattle feedlot was estimated as dollar 13.90 per animal, not including labor and associated handling costs.

BACKGROUND: The domestic pig (Sus scrofa) is important both as a food source and as a biomedical model given its similarity in size, anatomy, physiology, metabolism, pathology, and pharmacology to humans. The draft reference genome (Sscrofa10.2) of a purebred Duroc female pig established using older clone-based sequencing methods was incomplete, and unresolved redundancies, short-range order and orientation errors, and associated misassembled genes limited its utility. RESULTS: We present 2 annotated highly contiguous chromosome-level genome assemblies created with more recent long-read technologies and a whole-genome shotgun strategy, 1 for the same Duroc female (Sscrofa11.1) and 1 for an outbred, composite-breed male (USMARCv1.0). Both assemblies are of substantially higher (>90-fold) continuity and accuracy than Sscrofa10.2. CONCLUSIONS: These highly contiguous assemblies plus annotation of a further 11 short-read assemblies provide an unprecedented view of the genetic make-up of this important agricultural and biomedical model species. We propose that the improved Duroc assembly (Sscrofa11.1) become the reference genome for genomic research in pigs.

The cattle rumen has a diverse microbial ecosystem that is essential for the host to digest plant material. Extremes in body weight (BW) gain in mice and humans have been associated with different intestinal microbial populations. The objective of this study was to characterize the microbiome of the cattle rumen among steers differing in feed efficiency. Two contemporary groups of steers (n=148 and n=197) were fed a ration (dry matter basis) of 57.35% dry-rolled corn, 30% wet distillers grain with solubles, 8% alfalfa hay, 4.25% supplement, and 0.4% urea for 63 days. Individual feed intake (FI) and BW gain were determined. Within contemporary group, the four steers within each Cartesian quadrant were sampled (n=16/group) from the bivariate distribution of average daily BW gain and average daily FI. Bacterial 16S rRNA gene amplicons were sequenced from the harvested bovine rumen fluid samples using next-generation sequencing technology. No significant changes in diversity or richness were indicated, and UniFrac principal coordinate analysis did not show any separation of microbial communities within the rumen. However, the abundances of relative microbial populations and operational taxonomic units did reveal significant differences with reference to feed efficiency groups. Bacteroidetes and Firmicutes were the dominant phyla in all ruminal groups, with significant population shifts in relevant ruminal taxa, including phyla Firmicutes and Lentisphaerae, as well as genera Succiniclasticum, Lactobacillus, Ruminococcus, and Prevotella. This study suggests the involvement of the rumen microbiome as a component influencing the efficiency of weight gain at the 16S level, which can be utilized to better understand variations in microbial ecology as well as host factors that will improve feed efficiency.

We report the most extensive genetic linkage map for a livestock species produced to date. We have linked 376 microsatellite (MS) loci with seven restriction fragment length polymorphic loci in a backcross reference population. The 383 markers were placed into 24 linkage groups which span 1997 cM. Seven additional MS did not fall into a linkage group. Linkage groups are assigned to 13 autosomes and the X chromosome (haploid n = 19). This map provides the basis for genetic analysis of quantitative inheritance of phenotypic and physiologic traits in swine.