Association of ARRDC3 and NFIA variants with bovine congestive heart failure in feedlot cattle

Background

Bovine congestive heart failure (BCHF) has become increasingly prevalent among feedlot cattle in the Western Great Plains of North America with up to 7% mortality in affected herds. BCHF is an untreatable complex condition involving pulmonary hypertension that culminates in right ventricular failure and death. Genes associated with BCHF in feedlot cattle have not been previously identified. Our aim was to search for genomic regions associated with this disease.

Methods

A retrospective, matched case-control design with 102 clinical BCHF cases and their unaffected pen mates was used in a genome-wide association study. Paired nominal data from approximately 560,000 filtered single nucleotide polymorphisms (SNPs) were analyzed with McNemar's test.

Results

Two independent genomic regions were identified as having the most significant association with BCHF: the arrestin domain-containing protein 3 gene ( ARRDC3), and the nuclear factor IA gene ( NFIA, mid- p-values, 1x10 -8 and 2x10 -7, respectively). Animals with two copies of risk alleles at either gene were approximately eight-fold more likely to have BCHF than their matched pen mates with either one or zero risk alleles at both genes (CI 95 = 3-17). Further, animals with two copies of risk alleles at both genes were 28-fold more likely to have BCHF than all others ( p-value = 1×10 -7, CI 95 = 4-206). A missense variant in ARRDC3 (C182Y) represents a potential functional variant since the C182 codon is conserved among all other jawed vertebrate species observed. A two-SNP test with markers in both genes showed 29% of 273 BCHF cases had homozygous risk genotypes in both genes, compared to 2.5% in 198 similar unaffected feedlot cattle. This and other DNA tests may be useful for identifying feedlot animals with the highest risk for BCHF in the environments described here.

Conclusions

Although pathogenic roles for variants in the ARRDC3 and NFIA genes are unknown, their discovery facilitates classifying animals by genetic risk and allows cattle producers to make informed decisions for selective breeding and animal health management.

Two bovine hepacivirus genome sequences from U.S. cattle

Bovine hepacivirus (BoHV) is closely related to the hepatitis C virus (HCV) in humans and can cause both acute and chronic liver infections in cattle. BoHV was first identified in Ghana and Germany in 2015 and since then it has been detected and characterized in other countries around the world, but no strains have been sequenced from U.S. cattle. To date, BoHV has been classified into 2 genotypes (1 and 2), with genotype 1 being further divided into 11 subtypes (A-K). However, the true genetic diversity of BoHV is likely underestimated given limited surveillance and a lack of published genome sequences. Here, we sequenced 2 nearly complete BoHV genomes from serum samples collected in 2019 from beef cattle in Missouri. Sequence comparisons and phylogenetic analysis showed that isolate MARC/2019/60 had high sequence homology with genotype 1, subtype E isolates from China. In contrast, isolate MARC/2019/50 represented a novel BoHV subtype within genotype 2. Thus, we report the first genomic characterization of BoHV isolates from U.S. cattle, and the second complete BoHV2 genome worldwide. This work increases our knowledge of the global genetic diversity of BoHV and demonstrates the co-circulation of divergent BoHV strains in U.S. cattle.

Escherichia coli O157:H7 tir 255 T > A allele strains differ in chromosomal and plasmid composition

Shiga toxin-producing Escherichia coli (STEC) O157:H7 strains with the T allele in the translocated intimin receptor polymorphism (tir) 255 A > T gene associate with human disease more than strains with an A allele; however, the allele is not thought to be the direct cause of this difference. We sequenced a diverse set of STEC O157:H7 strains (26% A allele, 74% T allele) to identify linked differences that might underlie disease association. The average chromosome and pO157 plasmid size and gene content were significantly greater within the tir 255 A allele strains. Eighteen coding sequences were unique to tir 255 A allele chromosomes, and three were unique to tir 255 T allele chromosomes. There also were non-pO157 plasmids that were unique to each tir 255 allele variant. The overall average number of prophages did not differ between tir 255 allele strains; however, there were different types between the strains. Genomic and mobile element variation linked to the tir 255 polymorphism may account for the increased frequency of the T allele isolates in human disease.

First gene-edited calf with reduced susceptibility to a major viral pathogen

Bovine viral diarrhea virus (BVDV) is one of the most important viruses affecting the health and well-being of bovine species throughout the world. Here, we used CRISPR-mediated homology-directed repair and somatic cell nuclear transfer to produce a live calf with a six amino acid substitution in the BVDV binding domain of bovine CD46. The result was a gene-edited calf with dramatically reduced susceptibility to infection as measured by reduced clinical signs and the lack of viral infection in white blood cells. The edited calf has no off-target edits and appears normal and healthy at 20 months of age without obvious adverse effects from the on-target edit. This precision bred, proof-of-concept animal provides the first evidence that intentional genome alterations in the CD46 gene may reduce the burden of BVDV-associated diseases in cattle and is consistent with our stepwise, in vitro and ex vivo experiments with cell lines and matched fetal clones.

Recent Emergence of Bovine Coronavirus Variants with Mutations in the Hemagglutinin-Esterase Receptor Binding Domain in U.S. Cattle

Bovine coronavirus (BCoV) has spilled over to many species, including humans, where the host range variant coronavirus OC43 is endemic. The balance of the opposing activities of the surface spike (S) and hemagglutinin-esterase (HE) glycoproteins controls BCoV avidity, which is critical for interspecies transmission and host adaptation. Here, 78 genomes were sequenced directly from clinical samples collected between 2013 and 2022 from cattle in 12 states, primarily in the Midwestern U.S. Relatively little genetic diversity was observed, with genomes having >98% nucleotide identity. Eleven isolates collected between 2020 and 2022 from four states (Nebraska, Colorado, California, and Wisconsin) contained a 12 nucleotide insertion in the receptor-binding domain (RBD) of the HE gene similar to one recently reported in China, and a single genome from Nebraska collected in 2020 contained a novel 12 nucleotide deletion in the HE gene RBD. Isogenic HE proteins containing either the insertion or deletion in the HE RBD maintained esterase activity and could bind bovine submaxillary mucin, a substrate enriched in the receptor 9-O-acetylated-sialic acid, despite modeling that predicted structural changes in the HE R3 loop critical for receptor binding. The emergence of BCoV with structural variants in the RBD raises the possibility of further interspecies transmission.

Association of ARRDC3 and NFIA variants with bovine congestive heart failure in feedlot cattle

Background:   Bovine congestive heart failure (BCHF) has become increasingly prevalent among feedlot cattle in the Western Great Plains of North America with up to 7% mortality in affected herds.  BCHF is an untreatable complex condition involving pulmonary hypertension that culminates in right ventricular failure and death. Genes associated with BCHF in feedlot cattle have not been previously identified.  Our aim was to search for genomic regions associated with this disease. Methods:  A retrospective, matched case-control design with 102 clinical BCHF cases and their unaffected pen mates was used in a genome-wide association study.  Paired nominal data from approximately 560,000 filtered single nucleotide polymorphisms (SNPs) were analyzed with McNemar’s test. Results:  The most significant genome-wide association was in the arrestin domain-containing protein 3 gene (ARRDC3), followed by the nuclear factor IA gene (NFIA, mid-p-values, 1x10-8 and 2x10-7, respectively).  Animals with homozygous risk alleles at either gene were approximately eight-fold more likely to have BCHF than their matched pen mates without those risk alleles (CI95 = 3-17).  Animals with homozygous risk alleles at both genes were 28-fold more likely to have BCHF than all others (p-value = 1x10-7, CI95 = 4-206).  A linked missense variant in ARRDC3 (C182Y) represents a potential functional variant as the C182 codon is conserved among all other jawed vertebrate species observed.  A DNA test with two markers showed 29% of 273 BCHF cases had homozygous risk alleles in both genes, compared to 2.5% in 198 similar unaffected feedlot cattle. This DNA test may be useful for identifying feedlot animals with the highest risk for BCHF in the environments described here. Conclusions: Although pathogenic roles for ARRDC3 and NFIA variants associated with BCHF are unknown, their discovery facilitates classifying animals by genetic risk and allows cattle producers to make informed decisions for selective breeding and animal health management.

Detection of bovine inflammatory cytokines IL-1β, IL-6, and TNF-α with a multiplex electrochemiluminescent assay platform

Commercially available bovine-specific assays are limited in number, and multiplex assays for this species are rare. Our objective was to develop a multiplex assay for the bovine inflammatory cytokines IL-1β, IL-6, and TNF-α using the Meso Scale Discovery U-PLEX platform. "Do-It-Yourself" ELISA kits that contained polyclonal antibodies, both unlabeled and biotinylated, and the specific recombinant bovine cytokine standard, were purchased for each of these three cytokines. The biotinylated antibodies were coupled to linkers that bind to specific locations within each well of the U-PLEX plate. Unique linkers were used for each of the cytokines. The unlabeled antibodies were conjugated with electrochemiluminescent labels to serve as detection antibodies. Each cytokine assay was optimized individually prior to performing an optimization on the multiplex assay containing reagents for all three cytokines. To calculate cytokine concentrations, standard curves were developed using the recombinant cytokines and were run concurrently on each plate. Standard curves for IL-1β and TNF-α were run at concentrations ranging from 0 to 50,000 pg/mL, and for IL-6 from 0 to 10,000 pg/mL. The average lowest level of detection concentration measured by the standard curves were 5.3 pg/mL, 0.92 pg/mL, and 22.34 pg/mL for IL-1β, IL-6, and TNF-α respectively, as determined by data from seven plates containing bovine plasma samples from a combination of healthy and diseased cattle. The U-PLEX platform was a viable means to develop custom analyte- and species-specific multiplex assays using privately developed or purchased sets of commercially available reagents.

Cytokine and Haptoglobin Profiles From Shipping Through Sickness and Recovery in Metaphylaxis- or Un-Treated Cattle

Fifty-six head of cattle, 28 animals with bovine respiratory disease complex (BRDC), and 28 healthy animals that were matched by treatment, sale barn of origin, day, and interactions among these variables, were identified from a population of 180 animals (60 each purchased at three sale barns located in Missouri, Tennessee, and Kentucky) enrolled in a study comparing animals receiving metaphylaxis to saline-treated controls. Cattle were transported to a feedlot in KS and assigned to treatment group. Blood samples were collected at Day 0 (at sale barn), Day 1, Day 9, and Day 28 (at KS feedlot), and transported to the US Meat Animal Research Center in Clay Center, NE where plasma was harvested and stored at −80°C until assayed for the cytokines IFN-γ, IL-1β, IL-6, and TNF-α, and the acute stress protein haptoglobin (HPT). Our objectives were to determine if cytokine and haptoglobin profiles differed between control and metaphylaxis treatment groups over time, and if profiles differed between animals presenting with BRDC and those that remained healthy. There was no difference between the treated animals and their non-treated counterparts for any of the analytes measured. Sale barn of origin tended to affect TNF-α concentration. Differences for all analytes changed over days, and on specific days was associated with state of origin and treatment. The Treatment by Day by Case interaction was significant for HPT. The analyte most associated with BRDC was HPT on D9, possibly indicating that many of the cattle were not exposed to respiratory pathogens prior to entering the feedlot.

Rapid estimation of Salmonella enterica contamination level in ground beef - Application of the time-to-positivity method using a combination of molecular detection and direct plating

Little progress has been made in decreasing the incidence rate of salmonellosis in the US over the past decade. Mitigating the contribution of contaminated raw meat to the salmonellosis incidence rate requires rapid methods for quantifying Salmonella, so that highly contaminated products can be removed before entering the food chain. Here we evaluated the use of Time-to-Positivity (TTP) as a rapid, semi-quantitative approach for estimating Salmonella contamination levels in ground beef. Growth rates of 14 Salmonella strains (inoculated at log 1 to -2 CFU/g) were characterized in lean ground beef mTSB enrichments and time-to-detection was determined using culture and molecular detection methods. Enrichments were sampled at five timepoints and results were used to construct a prediction model of estimated contamination level by TTP (superscript indicates time in hours) defined as TTP⁴: ≥5 CFU/g; TTP⁶: ≤5, ≥1 CFU/g; TTP⁸: ≤1, ≥0.01 CFU/g; with samples negative at 8 h estimated ≤0.01 CFU/g. Model performance measures showed high sensitivity (100%) and specificity (83% and 93% for two detection methods) for samples with a TTP⁴, with false negative rates of 0%.

Classification of 16S rRNA reads is improved using a niche-specific database constructed by near-full length sequencing

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.

A Computational Method to Quantify the Effects of Slipped Strand Mispairing on Bacterial Tetranucleotide Repeats

The virulence and pathogenicity of bacterial pathogens are related to their adaptability to changing environments. One process enabling adaptation is based on minor changes in genome sequence, as small as a few base pairs, within segments of genome called simple sequence repeats (SSRs) that consist of multiple copies of a short sequence (from one to several nucleotides), repeated in series. SSRs are found in eukaryotes as well as prokaryotes, and length variation in them occurs at frequencies up to a million-fold higher than bacterial point mutations through the process of slipped strand mispairing (SSM) by DNA polymerase during replication. The characterization of SSR length by standard sequencing methods is complicated by the appearance of length variation introduced during the sequencing process that obscures the lower abundance repeat number variants in a population. Here we report a computational approach to correct for sequencing process-induced artifacts, validated for tetranucleotide repeats by use of synthetic constructs of fixed, known length. We apply this method to a laboratory culture of Histophilus somni, prepared from a single colony, and demonstrate that the culture consists of populations of distinct sequence phase and length variants at individual tetranucleotide SSR loci.

Comparative genomics of Salmonella enterica serovar Montevideo reveals lineage-specific gene differences that may influence ecological niche association

Salmonella enterica serovar Montevideo has been linked to recent foodborne illness outbreaks resulting from contamination of products such as fruits, vegetables, seeds and spices. Studies have shown that Montevideo also is frequently associated with healthy cattle and can be isolated from ground beef, yet human salmonellosis outbreaks of Montevideo associated with ground beef contamination are rare. This disparity fuelled our interest in characterizing the genomic differences between Montevideo strains isolated from healthy cattle and beef products, and those isolated from human patients and outbreak sources. To that end, we sequenced 13 Montevideo strains to completion, producing high-quality genome assemblies of isolates from human patients (n=8) or from healthy cattle at slaughter (n=5). Comparative analysis of sequence data from this study and publicly available sequences (n=72) shows that Montevideo falls into four previously established clades, differentially occupied by cattle and human strains. The results of these analyses reveal differences in metabolic islands, environmental adhesion determinants and virulence factors within each clade, and suggest explanations for the infrequent association between bovine isolates and human illnesses.

Transcriptome Responses to Respiratory Virus Infection of Pigs within the Tracheobronchial Lymphnode Following Infection with PRRSV, PCV2 or IAV

Porcine reproductive and respiratory syndromevirus (PRRSV) is a major respiratory pathogen of swine that has become extremely costly to the swine industry worldwide, often causing losses in production and animal life due to their ease of spread. However, the intracellular changes that occur in pigs following viral respiratory infections are still scantily understood for PRRSV, as well as, other viral respiratory infections. The aim of this study was to acquire a better understanding of PRRS disease by comparing gene expression changes that occur in tracheobronchial lymph nodes (TBLN) of pigs infected with either porcine reproductive and respiratory syndrome virus (PRRSV), porcine circovirus type 2 (PCV2), or swine influenza virus (IAV) infections. The study identified and compared gene expression changes in the TBLN of 16 pigs following infection by PRRSV, PCV2, IAV, or sham inoculation. Total RNA was pooled for each group and time-point (1, 3, 6, and14 DPI) to make 16 libraries, for analysis by Digital Gene Expression Tag Profiling (DGETP). The data underwent standard filtering to generate a list of sequence tag raw counts that were then analyzed using multidimensional and differential expression statistical tests. The results showed that PRRSV, IAV and PCV-2 infections followed a clinical course typical of experimental infection of young pigs with these viruses. Gene expression results echoed this course, as well as, uncovered genes related to shared and unique host immune responses to the 3 viruses. By testing and observing the host response to other respiratory viruses, our study has elucidated similarities and differences that can assist in development of vaccine and therapeutics that shorten or prevent a chronic PRRSV infection.

Genomic signatures of Mannheimia haemolytica that associate with the lungs of cattle with respiratory disease, an integrative conjugative element, and antibiotic resistance genes

BACKGROUND

Mannheimia haemolytica typically resides in cattle as a commensal member of the upper respiratory tract microbiome. However, some strains can invade their lungs and cause respiratory disease and death, including those with multi-drug resistance. A nucleotide polymorphism typing system was developed for M. haemolytica from the genome sequences of 1133 North American isolates, and used to identify genetic differences between isolates from the lungs and upper respiratory tract of cattle with and without clinical signs of respiratory disease.

RESULTS

A total of 26,081 nucleotide polymorphisms were characterized after quality control filtering of 48,403 putative polymorphisms. Phylogenetic analyses of nucleotide polymorphism genotypes split M. haemolytica into two major genotypes (1 and 2) that each were further divided into multiple subtypes. Multiple polymorphisms were identified with alleles that tagged genotypes 1 or 2, and their respective subtypes. Only genotype 2 M. haemolytica associated with the lungs of diseased cattle and the sequence of a particular integrative and conjugative element (ICE). Additionally, isolates belonging to one subtype of genotype 2 (2b), had the majority of antibiotic resistance genes detected in this study, which were assorted into seven combinations that ranged from 1 to 12 resistance genes.

CONCLUSIONS

Typing of diverse M. haemolytica by nucleotide polymorphism genotypes successfully identified associations with diseased cattle lungs, ICE sequence, and antibiotic resistance genes. Management of cattle by their carriage of M. haemolytica could be an effective intervention strategy to reduce the prevalence of respiratory disease and supplemental needs for antibiotic treatments in North American herds.

Observations on macrolide resistance and susceptibility testing performance in field isolates collected from clinical bovine respiratory disease cases

The objectives of this study were; first, to describe gamithromycin susceptibility of Mannheimia haemolytica, Pasteurella multocida, and Histophilus somni isolated from cattle diagnosed with bovine respiratory disease (BRD) and previously treated with either gamithromycin for control of BRD (mass medication=MM) or sham-saline injected (control=CON); second, to describe the macrolide resistance genes present in genetically typed M. haemolytica isolates; third, use whole-genome sequencing (WGS) to correlate the phenotypic resistance and genetic determinants for resistance among M. haemolytica isolates. M. haemolytica (n=276), P. multocida (n=253), and H. somni (n=78) were isolated from feedlot cattle diagnosed with BRD. Gamithromycin susceptibility was determined by broth microdilution. Whole-genome sequencing was utilized to determine the presence/absence of macrolide resistance genes and to genetically type M. haemolytica. Generalized linear mixed models were built for analysis. There was not a significant difference between MM and CON groups in regards to the likelihood of culturing a resistant isolate of M. haemolytica or P. multocida. The likelihood of culturing a resistant isolate of M. haemolytica differed significantly by state of origin in this study. A single M. haemolytica genetic subtype was associated with an over whelming majority of the observed resistance. H. somni isolation counts were low and statistical models would not converge. Phenotypic resistance was predicted with high sensitivity and specificity by WGS. Additional studies to elucidate the relationships between phenotypic expression of resistance/genetic determinants for resistance and clinical response to antimicrobials are necessary to inform judicious use of antimicrobials in the context of relieving animal disease and suffering.

Resolving Bovine viral diarrhea virus subtypes from persistently infected U.S. beef calves with complete genome sequence

Bovine viral diarrhea virus (BVDV) is classified into 2 genotypes, BVDV-1 and BVDV-2, each of which contains distinct subtypes with genetic and antigenic variation. To effectively control BVDV by vaccination, it is important to know which subtypes of the virus are circulating and how their prevalence is changing over time. Accordingly, the purpose of our study was to estimate the current prevalence and diversity of BVDV subtypes from persistently infected (PI) beef calves in the central United States. Phylogenetic analysis of the 5'-UTR (5' untranslated region) for 119 virus strains revealed that a majority (82%) belonged to genotype 1b, and the remaining strains were distributed between genotypes 1a (9%) and 2 (8%); however, BVDV-2 subtypes could not be confidently resolved. Therefore, to better define the variability of U.S. BVDV isolates and further investigate the division of BVDV-2 isolates into subtypes, complete genome sequences were obtained for these isolates as well as representatives of BVDV-1a and -1b. Phylogenetic analyses of the complete coding sequence provided more conclusive genetic classification and revealed that U.S. BVDV-2 isolates belong to at least 3 distinct genetic groups that are statistically supported by both complete and individual coding gene analyses. These results show that a more complex set of BVDV-2 subtypes has been circulating in this region than was previously thought.

Reducing assembly complexity of microbial genomes with single-molecule sequencing

BACKGROUND

The short reads output by first- and second-generation DNA sequencing instruments cannot completely reconstruct microbial chromosomes. Therefore, most genomes have been left unfinished due to the significant resources required to manually close gaps in draft assemblies. Third-generation, single-molecule sequencing addresses this problem by greatly increasing sequencing read length, which simplifies the assembly problem.

RESULTS

To measure the benefit of single-molecule sequencing on microbial genome assembly, we sequenced and assembled the genomes of six bacteria and analyzed the repeat complexity of 2,267 complete bacteria and archaea. Our results indicate that the majority of known bacterial and archaeal genomes can be assembled without gaps, at finished-grade quality, using a single PacBio RS sequencing library. These single-library assemblies are also more accurate than typical short-read assemblies and hybrid assemblies of short and long reads.

CONCLUSIONS

Automated assembly of long, single-molecule sequencing data reduces the cost of microbial finishing to $1,000 for most genomes, and future advances in this technology are expected to drive the cost lower. This is expected to increase the number of completed genomes, improve the quality of microbial genome databases, and enable high-fidelity, population-scale studies of pan-genomes and chromosomal organization.

Evolutionary characterization of pig interferon-inducible transmembrane gene family and member expression dynamics in tracheobronchial lymph nodes of pigs infected with swine respiratory disease viruses

Studies have found that a cluster of duplicated gene loci encoding the interferon-inducible transmembrane proteins (IFITMs) family have antiviral activity against several viruses, including influenza A virus. The gene family has 5 and 7 members in humans and mice, respectively. Here, we confirm the current annotation of pig IFITM1, IFITM2, IFITM3, IFITM5, IFITM1L1 and IFITM1L4, manually annotated IFITM1L2, IFITM1L3, IFITM5L, IFITM3L1 and IFITM3L2, and provide expressed sequence tag (EST) and/or mRNA evidence, not contained with the NCBI Reference Sequence database (RefSeq), for the existence of IFITM6, IFITM7 and a new IFITM1-like (IFITM1LN) gene in pigs. Phylogenic analyses showed seven porcine IFITM genes with highly conserved human/mouse orthologs known to have anti-viral activity. Digital Gene Expression Tag Profiling (DGETP) of swine tracheobronchial lymph nodes (TBLN) of pigs infected with swine influenza virus (SIV), porcine pseudorabies virus, porcine reproductive and respiratory syndrome virus or porcine circovirus type 2 over 14 days post-inoculation (dpi) showed that gene expression abundance differs dramatically among pig IFITM family members, ranging from 0 to over 3000 tags per million. In particular, SIV up-regulated IFITM1 by 5.9 fold at 3 dpi. Bayesian framework further identified pig IFITM1 and IFITM3 as differentially expressed genes in the overall transcriptome analysis. In addition to being a component of protein complexes involved in homotypic adhesion, the IFITM1 is also associated with pathways related to regulation of cell proliferation and IFITM3 is involved in immune responses.

Effects of TMEM154 haplotypes 1 and 3 on susceptibility to ovine progressive pneumonia virus following natural exposure in sheep

Small ruminant lentiviruses (SRLV) adversely affect production and well-being of sheep and goats throughout much of the world. The SRLV, including ovine progressive pneumonia virus (OPPV) in North America, cause lifetime infections, and management procedures to eradicate or reduce disease prevalence are costly. Variants of ovine transmembrane protein 154 gene (TMEM154) affect susceptibility to OPPV. The primary experimental objective was to estimate additive and dominance effects of TMEM154 haplotypes 1 and 3 on susceptibility to OPPV infection following natural exposure. A group of 187 trial lambs was born and raised by mature, infected ewes to ensure natural exposure to OPPV. Parents of trial lambs were heterozygous for haplotypes 1 and 3, producing lambs with diplotypes "1 1," "1 3," and "3 3." A group of 20 sentinel lambs was born and raised by mature, uninfected ewes that were diplotype "1 1." Sentinel lambs had diplotypes "1 1" and "1 3," being sired by the same set of rams as trial lambs. Trial and sentinel lambs were comingled during the experiment. Lambs were weaned at 60 d of age, bled 1 wk after weaning, and thereafter at intervals of 4 or 5 wk until 9 mo of age when OPPV infection status was determined by use of a competitive enzyme-linked immunosorbent assay. Only 1 sentinel lamb became infected. Infection status of trial lambs was analyzed using logistic regression procedures to account for the binary nature of infection status and random effects of sires. Effects of sex, type of birth, type of rearing, age of dam, breed type of dam, and sires were not detected (P>0.20). Infection status was affected by diplotype of lamb (P=0.005), with additive (P=0.002) and dominance (P=0.052) effects identified. Predicted probabilities of infection for lambs with diplotypes "1 1," "1 3," and "3 3" were 0.094, 0.323, and 0.346, respectively. Confidence intervals for probabilities of infection for diplotypes "1 3" and "3 3" were similar, but distinct from diplotype "1 1." These results are consistent with complete dominance of haplotype 3 relative to haplotype 1. The probability of infection at 9 mo of age for lambs with either diplotype "1 3" or "3 3" averaged 3.56 times that of lambs with diplotype "1 1." Genetic susceptibility to OPPV infection can be reduced by selection to increase the frequency of haplotype 1, resulting in a greater proportion of lambs with diplotype "1 1."

Reactomes of porcine alveolar macrophages infected with porcine reproductive and respiratory syndrome virus

Porcine reproductive and respiratory syndrome (PRRS) has devastated pig industries worldwide for many years. It is caused by a small RNA virus (PRRSV), which targets almost exclusively pig monocytes or macrophages. In the present study, five SAGE (serial analysis of gene expression) libraries derived from 0 hour mock-infected and 6, 12, 16 and 24 hours PRRSV-infected porcine alveolar macrophages (PAMs) produced a total 643,255 sequenced tags with 91,807 unique tags. Differentially expressed (DE) tags were then detected using the Bayesian framework followed by gene/mRNA assignment, arbitrary selection and manual annotation, which determined 699 DE genes for reactome analysis. The DAVID, KEGG and REACTOME databases assigned 573 of the DE genes into six biological systems, 60 functional categories and 504 pathways. The six systems are: cellular processes, genetic information processing, environmental information processing, metabolism, organismal systems and human diseases as defined by KEGG with modification. Self-organizing map (SOM) analysis further grouped these 699 DE genes into ten clusters, reflecting their expression trends along these five time points. Based on the number one functional category in each system, cell growth and death, transcription processes, signal transductions, energy metabolism, immune system and infectious diseases formed the major reactomes of PAMs responding to PRRSV infection. Our investigation also focused on dominant pathways that had at least 20 DE genes identified, multi-pathway genes that were involved in 10 or more pathways and exclusively-expressed genes that were included in one system. Overall, our present study reported a large set of DE genes, compiled a comprehensive coverage of pathways, and revealed system-based reactomes of PAMs infected with PRRSV. We believe that our reactome data provides new insight into molecular mechanisms involved in host genetic complexity of antiviral activities against PRRSV and lays a strong foundation for vaccine development to control PRRS incidence in pigs.

Genetic testing for TMEM154 mutations associated with lentivirus susceptibility in sheep

In sheep, small ruminant lentiviruses cause an incurable, progressive, lymphoproliferative disease that affects millions of animals worldwide. Known as ovine progressive pneumonia virus (OPPV) in the U.S., and Visna/Maedi virus (VMV) elsewhere, these viruses reduce an animal’s health, productivity, and lifespan. Genetic variation in the ovine transmembrane protein 154 gene (TMEM154) has been previously associated with OPPV infection in U.S. sheep. Sheep with the ancestral TMEM154 haplotype encoding glutamate (E) at position 35, and either form of an N70I variant, were highly-susceptible compared to sheep homozygous for the K35 missense mutation. Our current overall aim was to characterize TMEM154 in sheep from around the world to develop an efficient genetic test for reduced susceptibility. The average frequency of TMEM154 E35 among 74 breeds was 0.51 and indicated that highly-susceptible alleles were present in most breeds around the world. Analysis of whole genome sequences from an international panel of 75 sheep revealed more than 1,300 previously unreported polymorphisms in a 62 kb region containing TMEM154 and confirmed that the most susceptible haplotypes were distributed worldwide. Novel missense mutations were discovered in the signal peptide (A13V) and the extracellular domains (E31Q, I74F, and I102T) of TMEM154. A matrix-assisted laser desorption/ionization–time-of flight mass spectrometry (MALDI-TOF MS) assay was developed to detect these and six previously reported missense and two deletion mutations in TMEM154. In blinded trials, the call rate for the eight most common coding polymorphisms was 99.4% for 499 sheep tested and 96.0% of the animals were assigned paired TMEM154 haplotypes (i.e., diplotypes). The widespread distribution of highly-susceptible TMEM154 alleles suggests that genetic testing and selection may improve the health and productivity of infected flocks.

Next-generation sequencing of the porcine skeletal muscle transcriptome for computational prediction of microRNA gene targets

BACKGROUND

MicroRNA are a class of small RNAs that regulate gene expression by inhibiting translation of protein encoding transcripts through targeting of a microRNA-protein complex by base-pairing of the microRNA sequence to cognate recognition sequences in the 3' untranslated region (UTR) of the mRNA. Target identification for a given microRNA sequence is generally accomplished by informatics analysis of predicted mRNA sequences present in the genome or in databases of transcript sequence for the tissue of interest. However, gene models for porcine skeletal muscle transcripts in current databases, specifically complete sequence of the 3' UTR, are inadequate for this exercise.

METHODOLOGY/PRINCIPAL FINDINGS

To provide data necessary to identify gene targets for microRNA in porcine skeletal muscle, normalized cDNA libraries were sequenced using Roche 454 GS-FLX pyrosequencing and de novo assembly of transcripts enriched in the 3' UTR was performed using the MIRA sequence assembly program. Over 725 million bases of sequence were generated, which assembled into 18,202 contigs. Sequence reads were mapped to a 3' UTR database containing porcine sequences. The 3' UTR that mapped to the database were examined to predict targets for previously identified microRNA that had been separately sequenced from the same porcine muscle sample used to generate the cDNA libraries. For genes with microRNA-targeted 3' UTR, KEGG pathways were computationally determined in order to identify potential functional effects of these microRNA-targeted transcripts.

CONCLUSIONS

Through next-generation sequencing of transcripts expressed in skeletal muscle, mapping reads to a 3' UTR database, and prediction of microRNA target sites in the 3' UTR, our results identified genes expressed in porcine skeletal muscle and predicted the microRNA that target these genes. Additionally, identification of pathways regulated by these microRNA-targeted genes provides us with a set of genes that can be further evaluated for their potential role in skeletal muscle development and growth.

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.

Phylogeny of Shiga toxin-producing Escherichia coli O157 isolated from cattle and clinically ill humans

Cattle are a major reservoir for Shiga toxin-producing Escherichia coli O157 (STEC O157) and harbor multiple genetic subtypes that do not all associate with human disease. STEC O157 evolved from an E. coli O55:H7 progenitor; however, a lack of genome sequence has hindered investigations on the divergence of human- and/or cattle-associated subtypes. Our goals were to 1) identify nucleotide polymorphisms for STEC O157 genetic subtype detection, 2) determine the phylogeny of STEC O157 genetic subtypes using polymorphism-derived genotypes and a phage insertion typing system, and 3) compare polymorphism-derived genotypes identified in this study with pulsed field gel electrophoresis (PFGE), the current gold standard for evaluating STEC O157 diversity. Using 762 nucleotide polymorphisms that were originally identified through whole-genome sequencing of 189 STEC O157 human- and cattle-isolated strains, we genotyped a collection of 426 STEC O157 strains. Concatenated polymorphism alleles defined 175 genotypes that were tagged by a minimal set of 138 polymorphisms. Eight major lineages of STEC O157 were identified, of which cattle are a reservoir for seven. Two lineages regularly harbored by cattle accounted for the majority of human disease in this study, whereas another was rarely represented in humans and may have evolved toward reduced human virulence. Notably, cattle are not a known reservoir for E. coli O55:H7 or STEC O157:H⁻ (the first lineage to diverge within the STEC O157 serogroup), which both cause human disease. This result calls into question how cattle may have originally acquired STEC O157. The polymorphism-derived genotypes identified in this study did not surpass PFGE diversity assessed by BlnI and XbaI digestions in a subset of 93 strains. However, our results show that they are highly effective in assessing the evolutionary relatedness of epidemiologically unrelated STEC O157 genetic subtypes, including those associated with the cattle reservoir and human disease.

Bacterial community analysis of beef cattle feedlots reveals that pen surface is distinct from feces

The surface of beef cattle feedlot pens is commonly conceptualized as being packed uncomposted manure. Despite the important role that the feedlot pen may play in the transmission of veterinary and zoonotic pathogens, the bacterial ecology of feedlot surface material is not well understood. Our present study characterized the bacterial communities of the beef cattle feedlot pen surface material using 3647 full-length 16S rDNA sequences, and we compared the community composition of feedlot pens to the fecal source material. The feedlot surface composite was represented by members of the phylum Actinobacteria (42%), followed by Firmicutes (24%), Bacteroidetes (24%), and Proteobacteria (9%). The feedlot pen surface material bacterial communities were clearly distinct from those of the feces from animals in the same pen. Comparisons with previously published results of feces from the animals in the same pen reveal that, of 139 genera identified, only 25 were present in both habitats. These results indicate that, microbiologically, the feedlot pen surface material is separate and distinct from the fecal source material, suggesting that bacteria that originate in cattle feces face different selection pressures and survival challenges during their tenure in the feedlot pen, as compared to their residence in the gastrointestinal tract.

Virulence-associated and antibiotic resistance genes of microbial populations in cattle feces analyzed using a metagenomic approach

The bovine fecal microbiota impacts human food safety as well as animal health. Although the bacteria of cattle feces have been well characterized using culture-based and culture-independent methods, techniques have been lacking to correlate total community composition with community function. We used high throughput sequencing of total DNA extracted from fecal material to characterize general community composition and examine the repertoire of microbial genes present in beef cattle feces, including genes associated with antibiotic resistance and bacterial virulence. Results suggest that traditional 16S sequencing using "universal" primers to generate full-length sequence may under represent Acitinobacteria and Proteobacteria. Over eight percent (8.4%) of the sequences from our beef cattle fecal pool sample could be categorized as virulence genes, including a suite of genes associated with resistance to antibiotic and toxic compounds (RATC). This is a higher proportion of virulence genes found in Sargasso sea, chicken cecum, and cow rumen samples, but comparable to the proportion found in Antarctic marine derived lake, human fecal, and farm soil samples. The quantitative nature of metagenomic data, combined with the large number of RATC classes represented in samples from widely different habitats indicates that metagenomic data can be used to track relative amounts of antibiotic resistance genes in individual animals over time. Consequently, these data can be used to generate sample-specific and temporal antibiotic resistance gene profiles to facilitate an understanding of the ecology of the microbial communities in each habitat as well as the epidemiology of antibiotic resistant gene transport between and among habitats.

An atlas of bovine gene expression reveals novel distinctive tissue characteristics and evidence for improving genome annotation

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.

The genome sequence of taurine cattle: a window to ruminant biology and evolution

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.

A sequencing strategy for identifying variation throughout the prion gene of BSE-affected cattle

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.

Prion gene haplotypes of U.S. cattle

BACKGROUND

Bovine spongiform encephalopathy (BSE) is a fatal neurological disorder characterized by abnormal deposits of a protease-resistant isoform of the prion protein. Characterizing linkage disequilibrium (LD) and haplotype networks within the bovine prion gene (PRNP) is important for 1) testing rare or common PRNP variation for an association with BSE and 2) interpreting any association of PRNP alleles with BSE susceptibility. The objective of this study was to identify polymorphisms and haplotypes within PRNP from the promoter region through the 3'UTR in a diverse sample of U.S. cattle genomes.

RESULTS

A 25.2-kb genomic region containing PRNP was sequenced from 192 diverse U.S. beef and dairy cattle. Sequence analyses identified 388 total polymorphisms, of which 287 have not previously been reported. The polymorphism alleles define PRNP by regions of high and low LD. High LD is present between alleles in the promoter region through exon 2 (6.7 kb). PRNP alleles within the majority of intron 2, the entire coding sequence and the untranslated region of exon 3 are in low LD (18.0 kb). Two haplotype networks, one representing the region of high LD and the other the region of low LD yielded nineteen different combinations that represent haplotypes spanning PRNP. The haplotype combinations are tagged by 19 polymorphisms (htSNPS) which characterize variation within and across PRNP.

CONCLUSION

The number of polymorphisms in the prion gene region of U.S. cattle is nearly four times greater than previously described. These polymorphisms define PRNP haplotypes that may influence BSE susceptibility in cattle.

Characterization of 954 bovine full-CDS cDNA sequences

Background

Genome assemblies rely on the existence of transcript sequence to stitch together contigs, verify assembly of whole genome shotgun reads, and annotate genes. Functional genomics studies also rely on transcript sequence to create expression microarrays or interpret digital tag data produced by methods such as Serial Analysis of Gene Expression (SAGE). Transcript sequence can be predicted based on reconstruction from overlapping expressed sequence tags (EST) that are obtained by single-pass sequencing of random cDNA clones, but these reconstructions are prone to errors caused by alternative splice forms, transcripts from gene families with related sequences, and expressed pseudogenes. These errors confound genome assembly and annotation. The most useful transcript sequences are derived by complete insert sequencing of clones containing the entire length, or at least the full protein coding sequence (CDS) portion, of the source mRNA. While the bovine genome sequencing initiative is nearing completion, there is currently a paucity of bovine full-CDS mRNA and protein sequence data to support bovine genome assembly and functional genomics studies. Consequently, the production of high-quality bovine full-CDS cDNA sequences will enhance the bovine genome assembly and functional studies of bovine genes and gene products. The goal of this investigation was to identify and characterize the full-CDS sequences of bovine transcripts from clones identified in non-full-length enriched cDNA libraries. In contrast to several recent full-length cDNA investigations, these full-CDS cDNAs were selected, sequenced, and annotated without the benefit of the target organism's genomic sequence, by using comparison of bovine EST sequence to existing human mRNA to identify likely full-CDS clones for full-length insert cDNA (FLIC) sequencing.

Results

The predicted bovine protein lengths, 5' UTR lengths, and Kozak consensus sequences from 954 bovine FLIC sequences (bFLICs; average length 1713 nt, representing 762 distinct loci) are all consistent with previously sequenced mammalian full-length transcripts.

Conclusion

In most cases, the bFLICs span the entire CDS of the genes, providing the basis for creating predicted bovine protein sequences to support proteomics and comparative evolutionary research as well as functional genomics and genome annotation. The results demonstrate the utility of the comparative approach in obtaining predicted protein sequences in other species.

Use of bovine single nucleotide polymorphism markers to verify sample tracking in beef processing

OBJECTIVE

To determine whether a selected set of 20 single nucleotide polymorphism (SNP) markers derived from beef cattle populations can be used to verify sample tracking in a commercial slaughter facility that processes primarily market (ie, culled) dairy cows.

DESIGN

Prospective, blinded validation study.

ANIMALS

165 cows and 3 bulls from 18 states (82% Holstein, 8% other dairy breeds, and 10% beef breeds).

PROCEDURE

Blood was collected by venipuncture from randomly chosen animals just prior to slaughter. The purported corresponding liver samples were collected during beef processing, and genotype profiles were obtained for each sample.

RESULTS

On the basis of SNP allele frequencies in these cattle, the mean probability that 2 randomly selected individuals would possess identical genotypes at all 20 loci was 4.3 x 10(-8). Thus, the chance of a coincidental genotype match between 2 animals was 1 in 23 million. Genotype profiles confirmed appropriate matching for 152 of the 168 (90.5%) purported blood-liver sample pairs and revealed mismatching for 16 (9.5%) pairs. For the 16 mismatched sample pairs, 33% to 76% of the 20 SNP genotypes did not match (mean, 52%). Discordance that could be attributed to genotyping error was estimated to be < 1% on the basis of results for split samples.

CONCLUSIONS AND CLINICAL RELEVANCE

Results suggest that this selected set of 20 bovine SNP markers is sufficiently informative to verify accuracy of sample tracking in slaughter plants that process beef or dairy cattle. These or similar SNP markers may facilitate high-throughput, DNA-based, traceback programs designed to detect drug residues in tissues, control of animal diseases, and enhance food safety.

Linkage mapping bovine EST-based SNP

Background

Existing linkage maps of the bovine genome primarily contain anonymous microsatellite markers. These maps have proved valuable for mapping quantitative trait loci (QTL) to broad regions of the genome, but more closely spaced markers are needed to fine-map QTL, and markers associated with genes and annotated sequence are needed to identify genes and sequence variation that may explain QTL.

Results

Bovine expressed sequence tag (EST) and bacterial artificial chromosome (BAC)sequence data were used to develop 918 single nucleotide polymorphism (SNP) markers to map genes on the bovine linkage map. DNA of sires from the MARC reference population was used to detect SNPs, and progeny and mates of heterozygous sires were genotyped. Chromosome assignments for 861 SNPs were determined by twopoint analysis, and positions for 735 SNPs were established by multipoint analyses. Linkage maps of bovine autosomes with these SNPs represent 4585 markers in 2475 positions spanning 3058 cM . Markers include 3612 microsatellites, 913 SNPs and 60 other markers. Mean separation between marker positions is 1.2 cM. New SNP markers appear in 511 positions, with mean separation of 4.7 cM. Multi-allelic markers, mostly microsatellites, had a mean (maximum) of 216 (366) informative meioses, and a mean 3-lod confidence interval of 3.6 cM Bi-allelic markers, including SNP and other marker types, had a mean (maximum) of 55 (191) informative meioses, and were placed within a mean 8.5 cM 3-lod confidence interval. Homologous human sequences were identified for 1159 markers, including 582 newly developed and mapped SNP.

Conclusion

Addition of these EST- and BAC-based SNPs to the bovine linkage map not only increases marker density, but provides connections to gene-rich physical maps, including annotated human sequence. The map provides a resource for fine-mapping quantitative trait loci and identification of positional candidate genes, and can be integrated with other data to guide and refine assembly of bovine genome sequence. Even after the bovine genome is completely sequenced, the map will continue to be a useful tool to link observable phenotypes and animal genotypes to underlying genes and molecular mechanisms influencing economically important beef and dairy traits.

Integrating linkage and radiation hybrid mapping data for bovine chromosome 15

BACKGROUND

Bovine chromosome (BTA) 15 contains a quantitative trait loci (QTL) for meat tenderness, as well as several breaks in synteny with human chromosome (HSA) 11. Both linkage and radiation hybrid (RH) maps of BTA 15 are available, but the linkage map lacks gene-specific markers needed to identify genes underlying the QTL, and the gene-rich RH map lacks associations with marker genotypes needed to define the QTL. Integrating the maps will provide information to further explore the QTL as well as refine the comparative map between BTA 15 and HSA 11. A recently developed approach to integrating linkage and RH maps uses both linkage and RH data to resolve a consensus marker order, rather than aligning independently constructed maps. Automated map construction procedures employing this maximum-likelihood approach were developed to integrate BTA RH and linkage data, and establish comparative positions of BTA 15 markers with HSA 11 homologs.

RESULTS

The integrated BTA 15 map represents 145 markers; 42 shared by both data sets, 36 unique to the linkage data and 67 unique to RH data. Sequence alignment yielded comparative positions for 77 bovine markers with homologs on HSA 11. The map covers approximately 32% of HSA 11 sequence in five segments of conserved synteny, another 15% of HSA 11 is shared with BTA 29. Bovine and human order are consistent in portions of the syntenic segments, but some rearrangement is apparent. Comparative positions of gene markers near the meat tenderness QTL indicate the region includes separate segments of HSA 11. The two microsatellite markers flanking the QTL peak are between defined syntenic segments.

CONCLUSIONS

Combining data to construct an integrated map not only consolidates information from different sources onto a single map, but information contributed from each data set increases the accuracy of the map. Comparison of bovine maps with well annotated human sequence can provide useful information about genes near mapped bovine markers, but bovine gene order may be different than human. Procedures to connect genetic and physical mapping data, build integrated maps for livestock species, and connect those maps to more fully annotated sequence can be automated, facilitating the maintenance of up-to-date maps, and providing a valuable tool to further explore genetic variation in livestock.

Positional candidate gene selection from livestock EST databases using Gene Ontology

MOTIVATION

The number of expressed sequence tags (ESTs) in GenBank has now surpassed 200,000 for cattle and 100,000 for swine. The Institute of Genome Research (TIGR) has organized these sequences into approximately 60,000 non-redundant consensus sequences (identified by TIGR Gene Indices) for cattle and 40,000 for swine. Anonymous ESTs are of limited value unless they are connected to function. Functional information is difficult to manage electronically because of heterogeneity of meaning and form among databases. The Gene Ontology (GO) Consortium has produced ontologies for gene function with consistent meaning and form across species. Linking livestock EST to gene function through similarity with sequences from other annotation-rich mammals could accelerate: (1) the discovery of positional candidate genes underlying a livestock quantitative trait locus (QTL) and (2) comparative mapping between livestock and other mammals (e.g. humans, mouse and rat). We initiated this investigation to determine if incorporation of the GO into the annotation process could accelerate livestock positional candidate gene discovery.

RESULTS

We have associated livestock ESTs with GO nodes through sequence similarity to the NCBI Reference Sequences (RefSeq). Positional candidate genes are identified within minutes that otherwise required days. The schema described here accommodates queries that return GO nodes from terms familiar to biologists, such as gene name, alternate/alias symbol, and OMIM phenotype.

AVAILABILITY

Scripts and schema are available on request from the authors.

Selection and use of SNP markers for animal identification and paternity analysis in U.S. beef cattle

DNA marker technology represents a promising means for determining the genetic identity and kinship of an animal. Compared with other types of DNA markers, single nucleotide polymorphisms (SNPs) are attractive because they are abundant, genetically stable, and amenable to high-throughput automated analysis. In cattle, the challenge has been to identify a minimal set of SNPs with sufficient power for use in a variety of popular breeds and crossbred populations. This report describes a set of 32 highly informative SNP markers distributed among 18 autosomes and both sex chromosomes. Informativity of these SNPs in U.S. beef cattle populations was estimated from the distribution of allele and genotype frequencies in two panels: one consisting of 96 purebred sires representing 17 popular breeds, and another with 154 purebred American Angus from six herds in four Midwestern states. Based on frequency data from these panels, the estimated probability that two randomly selected, unrelated individuals will possess identical genotypes for all 32 loci was 2.0 x 10(-13) for multi-breed composite populations and 1.9 x 10(-10) for purebred Angus populations. The probability that a randomly chosen candidate sire will be excluded from paternity was estimated to be 99.9% and 99.4% for the same respective populations. The DNA immediately surrounding the 32 target SNPs was sequenced in the 96 sires of the multi-breed panel and found to contain an additional 183 polymorphic sites. Knowledge of these additional sites, together with the 32 target SNPs, allows the design of robust, accurate genotype assays on a variety of high-throughput SNP genotyping platforms.