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

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

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

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

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

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

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

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

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

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

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

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

Proposed model for the high rate of rearrangement and rapid migration observed in some IncA/C plasmid lineages.

IncA/C plasmids are a class of plasmids from the Enterobacteriaceae that are relatively large (49 to >180 kbp), that are readily transferred by conjugation, and that carry multiple antimicrobial resistance genes. Reconstruction of the phylogeny of these plasmids has been difficult because of the high rate of remodeling by recombination-mediated horizontal gene transfer (HGT). We hypothesized that evaluation of nucleotide polymorphisms relative to the rate of HGT would help to develop a clock to show whether anthropic practices have had significant influences on the lineages of the plasmid. A system was developed to rapidly sequence up to 191 known open reading frames from each of 39 recently isolated IncA/C plasmids from a diverse panel of Salmonella enterica and Escherichia coli strains. With these data plus sequences from GenBank, we were able to distinguish six distinct lineages that had extremely low numbers of polymorphisms within each lineage, especially among the largest group designated as group 1. Two regions, each about half the plasmid in size, could be distinguished with a separate lineal pattern. The distribution of group 1 showed that it has migrated extremely rapidly with fewer polymorphisms than can be expected in 2,000 years. Remodeling by frequent HGT was evident, with a pattern that appeared to have the highest rate just upstream of the putative conjugation origin of transfer (oriT). It seems likely that when an IncA/C plasmid is transferred by conjugation there is an opportunity for plasmid remodeling adjacent to the oriT, which was also adjacent to a multiple antimicrobial resistance gene cassette.

The Role of Aphid Abundance, Species Diversity, and Virus Titer in the Spread of Sweetpotato Potyviruses in Louisiana and Mississippi.

Sweet potato feathery mottle virus (SPFMV), Sweet potato virus G (SPVG), and Sweet potato virus 2 (SPV2) are sweetpotato (Ipomoea batatas) potyviruses nonpersistently transmitted by aphids. Our objective was to determine how aphid abundance, aphid species diversity, and virus titers relate to the spread of SPFMV, SPVG, and SPV2 in Louisiana and Mississippi sweetpotato fields. The most abundant aphid species were Aphis gossypii, Myzus persicae, Rhopalosiphum padi, and Therioaphis trifolii. Aphids were captured during the entire crop cycle but virus infection of sentinel plants occurred mainly during the months of June to August. SPFMV was more commonly detected than SPVG or SPV2 in sentinel plants. Virus titers for SPFMV were higher in samples beginning in late June. Because significant aphid populations were present during April to June when virus titers were low in sweetpotato and there was very little virus infection of sentinel plants, low virus titers may have limited aphid acquisition and transmission opportunities. This is the first study to comprehensively examine aphid transmission of potyviruses in sweetpotato crops in the United States and includes the first report of R. maidis and R. padi as vectors of SPFMV, though they were less efficient than A. gossypii or M. persicae.

A region on BTA14 that includes the positional candidate genes LYPLA1, XKR4 and TMEM68 is associated with feed intake and growth phenotypes in cattle(1).

Feed cost for beef cattle is the largest expense incurred by cattle producers. The development of genetic markers to enhance selection of more efficient animals that require less feed while still achieving acceptable levels of production has the potential to substantially reduce production costs. A genome-wide marker association approach based on the Illumina BovineSNP50 BeadChip™ was used to identify genomic regions affecting average daily feed intake (ADFI), average daily gain (ADG) and residual feed intake traits in a population of 1159 crossbred steers. This approach identified a region on BTA14 from 22.02 to 23.92 Mb containing several single-nucleotide polymorphisms (SNPs) that have significant association with at least one of the traits. Two genes in this region, lysophospholipase 1 (LYPLA1) and transmembrane protein 68 (TMEM68), appeared to be logical positional and functional candidate genes. LYPLA1 deacylates ghrelin, a hormone involved in the regulation of appetite in the rat stomach, while TMEM68 is expressed in bovine rumen, abomasum, intestine and adipose tissue in cattle, and likely affects lipid biosynthetic processes. SNPs lying in or near these two genes were identified by sequencing a subset of animals with extreme phenotypes. A total of 55 SNPs were genotyped and tested for association with the same population of steers. After correction for multiple testing, five markers within 22.79-22.84 Mb, located downstream of TMEM68, and between TMEM68 and the neighbouring gene XKR4, were significant for both ADFI and ADG. Genetic markers predictive of feed intake and weight gain phenotypes in this population of cattle may be useful for the identification and selection of animals that consume less feed, although further evaluation of these markers for effects on other production traits and validation in additional populations will be required.

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

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

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.

Association of single nucleotide polymorphisms in the ANKRA2 and CD180 genes with bovine respiratory disease and presence of Mycobacterium avium subsp. paratuberculosis(1).

The objective was to determine whether single nucleotide polymorphisms (SNPs) in the ANKRA2 and CD180 genes are associated with incidence of bovine respiratory disease (BRD) and presence of Mycobacterium avium subsp. paratuberculosis (MAP) in cattle. Two independent populations were used. The first population (BRD-affected; N = 90) was composed of 31 half-sib progeny, from a Brahman × Angus sire, that were treated for BRD. Untreated offspring from the sire were selected to serve as controls. The second population (MAP-infected) of 330 animals of unknown parentage was evaluated for the presence of MAP in ileocecal lymph node and classified as positive or negative. Markers in both genes were assessed for association in these two populations. In the BRD-affected population, five SNPs in the ANKRA2 gene were significantly associated (P < 0.05), and two SNPs were highly associated (P < 0.01) with incidence of BRD. In addition, two SNPs in the CD180 gene were found to be associated with this trait. In the MAP-infected population, one SNP in the ANKRA2 gene was significantly associated (P < 0.05) with the presence or absence of MAP, and a SNP in the CD180 gene was highly associated (P < 0.01) with the trait. Haplotypes, using significant markers, showed a positive association with both incidence of BRD (P = 0.0001) and with the presence of MAP (P = 0.0032). Markers in the ANKRA2 and CD180 genes are associated with the ability of the animal to cope with pathogens.

The glucose lowering effect of an oral insulin (Capsulin) during an isoglycaemic clamp study in persons with type 2 diabetes.

AIM: Randomized, open, single-centre, two-way crossover study comparing the pharmacokinetic (PK) and pharmacodynamic (PD) properties of subcutaneous (sc) regular human insulin (Actrapid) and oral insulin in a capsule form (Capsulin). METHODS: Sixteen persons (12 males) with type 2 diabetes on oral hypoglycaemic agents (OHAs) participated. Mean (s.d.) age 60.2 (5.5) years, BMI 28.3 (3.4) kg/m(2), haemoglobin A(1c) (HbA(1c)) 7.4% (1.1). Two 6-h isoglycaemic glucose clamp studies were conducted 11 days apart. All subjects received in random order 12U sc Actrapid on one clamp study day and either 150U or 300U Capsulin (Cap) on the other day. Glucose infusion rates (GIRs), plasma insulin and C-peptide concentrations were determined throughout each 6-h isoglycaemic clamp. Between the clamp study days, all patients received 150U Capsulin twice daily, dropping all their standard OHAs apart from metformin. Self-monitored blood glucose (SMBG) levels were taken four times a day between the clamp study days. RESULTS: Administration of either Actrapid or Capsulin (150 and 300U) increased GIRs reaching a maximum values at approximately 280-330 min. Overall values for maximum GIR values were higher for Actrapid than either dose of Capsulin (p < 0.05). The significantly greater systemic insulin concentrations following Actrapid were reflected in the AUC(0-6 h) (910 +/- 270 vs. 472 +/- 245 pmol h/L; 950 +/- 446 vs. 433 +/- 218 pmol h/L; both p < 0.05 for Actrapid vs. 150U Capsulin and 300U Capsulin respectively). No difference was observed between 150U and 300U Capsulin. During the repeat-dosing period, good safety and tolerability were observed with Capsulin, and SMBG levels remained stable. At the poststudy visit, significant falls in HbA(1c), weight and triglycerides were observed. CONCLUSIONS: Administration of the oral insulin Capsulin preparation demonstrated a significant hypoglycaemic action over a period of 6 h associated with only a small increase in circulating plasma insulin concentrations.

Markers on bovine chromosome 20 associated with carcass quality and composition traits and incidence of contracting infectious bovine keratoconjunctivitis.

The objective of this study was to use single nucleotide polymorphisms (SNP) located on bovine chromosome 20 to fine map a previously identified QTL associated with the incidence of infectious bovine keratoconjunctivitis (IBK). Crossbred steers (GPE 7; n = 539) derived from sires of 7 Bos taurus breeds and having veterinary records related to IBK were used to test the association of a total of 105 SNP located under the most relevant region of the QTL. Five SNP were significantly associated with IBK (P < 0.05), as animals inheriting differing genotypes from individual SNP exhibited significantly different incidence rates of IBK. The population also had numerous other phenotypes, supporting evaluation of association of the 105 markers with carcass traits to identify potential antagonistic effects of implementing a marker-assisted selection program for IBK susceptibility. An association of 2 SNP for marbling and tenderness was identified, along with 3 SNP associated with the percentage of carcasses classified as choice. Four SNP were significantly associated with fat yield, 2 SNP with longissimus muscle area, and 2 additional SNP with dressing percentage. The association of these markers indicates that the evaluated QTL region may, in fact, harbor the causative mutations responsible for the variation observed in IBK susceptibility and carcass quality and composition traits. Thus, further evaluation of SNP in this region is necessary in order to identify mutations accounting for the largest degree of variation for IBK and carcass traits.

Vibrational modes and the dynamic solvent effect in electron and proton transfer.

This article primarily reviews recent work on ultrafast experiments on excited state intramolecular electron and proton transfer, with an emphasis on experiments on chemical systems that have been analyzed theoretically. In particular, those systems that have been quantitatively characterized by static spectroscopy, which provides detailed information about the reaction potential energy surface and about other parameters that are necessary to make a direct comparison to theoretical predictions, are described.

A class III transcription factor composed of RNA.

It is generally assumed that the machinery that transcribes genes is composed entirely of polypeptides. However, in vitro transcription by silkworm RNA polymerase III requires a transcription factor that is not a polypeptide. This component, TFIIIR, is distinct from the previously identified transcription components: RNA polymerase III, and the accessory factors TFIIIA, TFIIIB, TFIIIC, and TFIIID. The newly discovered TFIIIR is a macromolecule that appears to be composed of RNA. It is resistant to heat, detergent, phenol, protease, and deoxyribonuclease, but it is sensitive to alkali and ribonuclease.

Evaluation of biochemical parameters and genetic markers for association with meat tenderness in South African feedlot cattle.

A large proportion of South African feedlot cattle are crossbreds of Brahman (BrX, Bos indicus), and Simmental (SiX, Bos taurus). A sample of 20 grain fed bulls from each of these crossbreeds was used to compare meat quality with that of the small frame indigenous Nguni (NgX, Sanga) by evaluating a variety of biochemical and genetic parameters previously shown to be associated with meat tenderness. Shear force values were generally high (5.6kg average at 14days post mortem), with SiX animals higher than BrX or NgX (P=0.051) despite higher calpastatin:calpain ratio in BrX (P<0.05). Calpain activity and cold shortening were both correlated with tenderness for all classes. The sample size was too small to accurately estimate genotypic effects of previously published markers in the CAST and CAPN1 genes, but the allele frequencies suggest that only modest progress would be possible in these South African crossbreds using these markers.

Recognition and management of dermatologic complications from interventional radiology procedures.

A variety of dermatologic complications can occur after interventional radiology procedures, including fluoroscopy-induced radiation dermatitis, thermal skin injury from tumor ablation, non-target embolization to the skin, allergic reactions related to interventional radiology procedures, and dermatitis and infections at catheter sites. Yet, interventional radiologists typically lack training in dermatology. This review focuses on recognition of dermatologic complications and introduces basic principles for management of these complications. By taking a more active role in the diagnosis, management, and follow-up of dermatologic complications, interventional radiologists can improve the care for patients suffering iatrogenic skin inury.

Cardiovascular exercise intervention improves the primary antibody response to keyhole limpet hemocyanin (KLH) in previously sedentary older adults.

Based upon a prior cross-sectional study, we hypothesized that an aerobic exercise intervention in sedentary older adults would improve a primary T cell-dependent immune response. Participants were a subset of older subjects from a large, ongoing exercise intervention study who were randomly assigned to either an aerobic exercise (Cardio, n=30, 68.9+0.8 years) or flexibility/balance (Flex, n=20, 69.9+1.2 years) intervention. The intervention consisted of either three aerobic sessions for 30-60 min at 55-70% VO(2 max) or two 60 min flexibility/balance sessions weekly for 10 months. Eight months into the intervention, samples were collected before intramuscular administration of KLH (125 microg), followed by sampling at 2, 3, and 6 weeks post-KLH. Serum anti-KLH IgM, IgG1, and IgG2 was measured by ELISA. Physiological and psychosocial measures were also assessed pre- and post-intervention. While there was no difference in the anti-KLH IgG2 response between groups, Cardio displayed significantly (p<0.05) higher anti-KLH IgG1 (at weeks 2, 3, and 6 post) and IgM responses when compared to Flex. Despite cardiovascular intervention-induced improvement in physical fitness (approximately 11% vs. 1% change in VO(2 peak) in Cardio vs. Flex, respectively), we found no relationship between improved fitness and enhanced anti-KLH antibody responses. Optimism, perceived stress, and affect were all associated with enhanced immune response. We have shown for the first time that cardiovascular training in previously sedentary elderly results in significantly higher primary IgG1 and IgM antibody responses, while having no effect on IgG2 production.

Effect of porcine reproductive and respiratory syndrome virus on porcine alveolar macrophage function as determined using serial analysis of gene expression (SAGE).

Porcine reproductive and respiratory syndrome virus (PRRSV) is a major pathogen of swine worldwide and causes considerable economic loss. The main target of infection is the porcine alveolar macrophage (PAM). Infection of PAMs by PRRSV causes significant changes in their function by mechanisms that are not understood. We have employed Serial Analysis of Gene Expression (SAGE) to examine the global expression of genes in PRRSV-infected PAMs. Total cellular RNAwas prepared from in vitro mock-infected and PRRSV strain VR-2332-infected PAMs at 0, 6, 12, 16 and 24 hours after infection, and subjected to SAGE analysis to obtain > 100,000 tags per time point. These sequences were processed to account for sequencing error before generating tag:count lists. These lists were deposited into a modified Identitag database for mapping to porcine and PRRSV genes. Identified unique mRNAtags were analyzed for their identity and relative abundance. Examination of the SAGE data indicated that there were changes in gene expression occurring in the PRRSV-infected PAMs over time post-infection. More than 400 unique tags with significantly altered expression levels were identified (p < 0.01 with Bonferroni correction). The validity and kinetics of expression of SAGE identified genes were evaluated using real-time RT-PCR.