Successful development of methodology for detection of hapten-specific contact hypersensitivity (CHS) memory in swine.

Hapten contact hypersensitivity (CHS) elicits a well-documented inflammation response that can be used to illustrate training of immune cells through hapten-specific CHS memory. The education of hapten-specific memory T cells has been well-established, recent research in mice has expanded the "adaptive" characteristic of a memory response from solely a function of the adaptive immune system, to innate cells as well. To test whether similar responses are seen in a non-rodent model, we used hapten-specific CHS to measure the ear inflammation response of outbred pigs to dinitrofluorobenzene (DNFB), oxazolone (OXA), or vehicle controls. We adapted mouse innate memory literature protocols to the domestic pig model. Animals were challenged up to 32 days post initial sensitization exposure to the hapten, and specific ear swelling responses to this challenge were significant for 7, 21, and 32 days post-sensitization. We established hapten-specific CHS memory exists in a non-rodent model. We also developed a successful protocol for demonstrating these CHS responses in a porcine system.

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.

Genomewide association of piglet responses to infection with one of two porcine reproductive and respiratory syndrome virus isolates.

Porcine reproductive and respiratory syndrome (PRRS) is a devastating disease in the swine industry. Identification of host genetic factors that enable selection for improved performance during PRRS virus (PRRSV) infection would reduce the impact of this disease on animal welfare and production efficiency. We conducted genomewide association study (GWAS) analyses of data from 13 trials of approximately 200 commercial crossbred nursery-age piglets that were experimentally infected with 1 of 2 type 2 isolates of PRRSV (NVSL 97-7985 [NVSL] and KS2006-72109 [KS06]). Phenotypes analyzed were viral load (VL) in blood during the first 21 d after infection (dpi) and weight gain (WG) from 0 to 42 dpi. We accounted for the previously identified QTL in the region on SSC4 in our models to increase power to identify additional regions. Many regions identified by single-SNP analyses were not identified using Bayes-B, but both analyses identified the same regions on SSC3 and SSC5 to be associated with VL in the KS06 trials and on SSC6 in the NVSL trials ( < 5 × 10); for WG, regions on SSC5 and SSC17 were associated in the NVSL trials ( < 3 × 10). No regions were identified with either method for WG in the KS06 trials. Except for the region on SSC4, which was associated with VL for both isolates (but only with WG for NVSL), identified regions did not overlap between the 2 PRRSV isolate data sets, despite high estimates of the genetic correlation between isolates for traits based on these data. We also identified genomic regions whose associations with VL or WG interacted with either PRRSV isolate or with genotype at the SSC4 QTL. Gene ontology (GO) annotation terms for genes located near moderately associated SNP ( < 0.003) were enriched for multiple immunologically (VL) and metabolism- (WG) related GO terms. The biological relevance of these regions suggests that, although it may increase the number of false positives, the use of single-SNP analyses and a relaxed threshold also increased the identification of true positives. In conclusion, although only the SSC4 QTL was associated with response to both PRRSV isolates, genes near associated SNP were enriched for the same GO terms across PRRSV isolates, suggesting that host responses to these 2 isolates are affected by the actions of many genes that function together in similar biological processes.

Identification of potential serum biomarkers to predict feed efficiency in young pigs.

Identification of biomarkers for feed efficiency in livestock will aid in the efficient production of high-quality protein to meet the demands of a growing population. The overall objective of this research was to identify biomarkers in serum for swine feed efficiency and to discover pathways affected by divergent selection for residual feed intake (RFI). Serum was collected from young pigs (between 35 and 42 d of age) from 2 lines of pigs that have been genetically selected to be either more efficient (low-RFI) or less efficient (high-RFI). After blood collection, during finishing, pigs from each line were placed on either a low-energy/high-fiber diet or a traditional high-energy/low-fiber diet to test for any diet effects on RFI. Subsets of 6 pigs per line within each diet were used in 3 independent experiments. Pigs with extreme RFI phenotypes from the low-energy/high-fiber diet were used to confirm the results from the first 2 comparisons. Two-dimensional difference in gel electrophoresis and mass spectrometry were used to identify proteins with different abundances between RFI line or finishing diet. Three proteins had consistent and significant ( < 0.05) RFI line differences for both diets: gelsolin, vitronectin, and serine protease inhibitor A3 (serpinA3). Abundance of gelsolin, a protein with roles in actin filament assembly and immune response, was greater in the more efficient low-RFI pigs (9 to 39%). Vitronectin was also more abundant in the low-RFI pigs (39 to 56%) and has known roles in blood homeostasis and may regulate adiposity. SerpinA3 is a member of a very large family of proteins referred to as serine protease inhibitors. A total of 14 spots that were more abundant in the low-RFI line, some at least twice as abundant, were identified as serpinA3. Multiple isoforms of serpinA3 have been reported (serpinA3-1 to serpinA3-4 in pigs and serpinA3-1 to serpinA3-8 in cattle) with serpinA3 having many different functions dependent on isoform. Gelsolin, vitronectin, and serpinA3 are 3 proteins that may play direct and important biological roles in the pathways that control RFI and, ultimately, feed efficiency through energy utilization and homeostasis. These data demonstrate that serum proteins can be a useful source of potential biomarkers for feed efficiency and provide information on pathways with distinct expression patterns between animals that differ in feed efficiency.

Assessing peripheral blood cell profile of Yorkshire pigs divergently selected for residual feed intake.

The cost of feed is a serious issue in the pork industry, contributing about 65 to 75% of the total production cost. To prevent economic losses and decreased productivity of the herd, it is important to select for animals that eat less for the same lean gain, or more efficient animals. Residual feed intake (RFI) is the difference between observed feed intake and expected feed intake based on estimated maintenance and production requirements. Selection for decreased RFI, or more efficient animals, is a potential solution to higher feed costs in pig production. However, animals that are highly selected for decreased RFI may have reduced energy input to the immune system and fail to withstand diseases and stressors after infection that negatively impact profitability. The objective of this study was to evaluate differences in circulating blood cell profiles at a young age between 2 lines of Yorkshire pigs that were divergently selected for RFI as well as the heritability of these traits, to investigate effects of selection for RFI on immune system parameters, and to identify potential biomarkers for feed efficiency. Previous work has shown that the 2 lines had diverged for IGF-1 in serum in young pigs and, therefore, this stage was investigated for other potential physiological differences. Blood samples were drawn for a complete blood count (CBC) analysis from 517 gilts and barrows, ages 35 to 42 d, across the 2 lines. In general, the low-RFI line had lower numbers of specific types of white blood cells but higher hemoglobin concentration and red blood cell volume compared to the high-RFI line. No significant correlations were found between CBC traits and RFI across and within the lines (0.05 < < 0.1). Of the 15 CBC traits that were measured, 3 were highly heritable (0.56 < < 0.62), 9 were moderately heritable (0.12 < < 0.47), and 3 were lowly heritable ( < 0.12), suggesting a substantial genetic component for CBC traits and that selection for CBC traits could be effective. Our results also show that selection for RFI has significantly impacted the number of circulating blood cells. In this experiment, we studied only healthy animals that were not under known pathogen challenge; therefore, our results cannot be directly applied to a disease challenge situation. Future work will be to challenge the animals and determine the effect of challenge on CBC levels.

Investigation of the efficacy of albumin removal procedures on porcine serum proteome profile.

Improving the ability to predict livestock performance using biomarkers will provide a benefit for livestock genetic evaluation and improvement. The most practical biological sample to screen for development of biomarkers is serum due to the ease of collection. However, protein profiles in serum are complex and dynamic. Strategies are needed to manage variation in serum proteins used for biomarker identification. Albumin is the most abundant protein in serum, comprising over 50% of the overall protein content, and has historically been depleted from serum before biomarker identification. The objective of this study was to investigate the use of gel-based proteomic techniques to evaluate the need for porcine albumin depletion in biomarker identification. Albumin is known to bind many proteins in the blood, thus potential biomarkers could be removed during albumin depletion. Using two-dimensional difference in gel electrophoresis (2D-DIGE), we show whole serum can be used for biomarker discovery. The data obtained show that albumin removal methods are effective for porcine sera. Over 85% of the protein spots resolved on at least half of the gels were changed in abundance between whole and albumin depleted sera. Of the 204 protein spots significantly altered in abundance, 59 were changed over 400%. However, albumin removal also altered the serum proteome in an unpredictable manner; in the depleted sera, 86 protein spots were increased in abundance and 118 were decreased. Furthermore, the abundance of 59.4% of the protein spots in the albumin depleted samples had a larger standard error than whole sera. However, the resolution of albumin in 2D-DIGE analysis of whole sera permitted the detection and quantification of substantial numbers of proteins. Thus, it is proposed that whole serum can be used in a gel-based proteomics system for the identification of porcine biomarkers.

Transcript profiles in longissimus dorsi muscle and subcutaneous adipose tissue: a comparison of pigs with different postweaning growth rates.

Although most pigs recover rapidly from stresses associated with the transition of weaning, a portion of the population lags behind their contemporaries in growth performance. The underlying biological and molecular mechanisms involved in postweaning differences in growth performance are poorly understood. The objective of this experiment was to use transcriptional profiling of skeletal muscle and adipose tissue to develop a better understanding of the metabolic basis for poor weaned-pig transition. A total of 1,054 pigs was reared in commercial conditions and weighed at birth, weaning, and 3 wk postweaning. Transition ADG (tADG) was calculated as the ADG for the 3-wk period postweaning. Nine pigs from both the lowest 10th percentile (low tADG) and the 60th to 70th percentile (high tADG) were harvested at 3 wk postweaning. Differential expression analysis was conducted in longissimus dorsi muscle (LM) and subcutaneous adipose tissue using RNA-Seq methodology. In LM, 768 transcripts were differentially expressed (DE), 327 with higher expression in low tADG and 441 with higher expression in high tADG pigs (q < 0.10). Expression patterns measured in LM by RNA-Seq were verified in 30 of 32 transcripts using quantitative PCR. No DE transcripts were identified in adipose tissue. To identify biological functions potentially underlying the effects of tADG on skeletal muscle metabolism and physiology, functional annotation analysis of the DE transcripts was conducted using DAVID and Pathway Studio analytic tools. The group of DE genes with lower expression in LM of low tADG pigs was enriched in genes with functions related to muscle contraction, glucose metabolism, cytoskeleton organization, muscle development, and response to hormone stimulus (enrichment score > 1.3). The list of DE genes with higher expression in low tADG LM was enriched in genes with functions related to protein catabolism (enrichment score > 1.3). Analysis of known gene-gene interactions identified possible regulators of these differences in gene expression in LM of high and low tADG pigs; these include forkhead box O1 (FOXO1), growth hormone (GH1), and the glucocorticoid receptor (NR3C1). Differences in gene expression between poor transitioning pigs and their contemporaries indicate a shift to decreased protein synthesis, increased protein degradation, and reduced glucose metabolism in the LM of low tADG pigs.

Analysis of blood leukocytes in a naturally occurring immunodeficiency of pigs shows the defect is localized to B and T cells.

Severe combined immunodeficiency (SCID) is the result of a set of inherited genetic defects which render components of the immune response nonfunctional. In Arabian horses, Jack Russell terriers, and mice, the disorder is a consequence of the absence of T and B lymphocytes, while natural killer (NK) cell and other leukocyte populations remain intact. Preliminary analysis of a naturally acquired form of inherited SCID in a line of pigs showed several defects in the architecture and composition of secondary lymphoid organs. In this study, a quantitative assessment of lymphocyte populations in affected and normal littermates showed depleted T or B lymphocyte populations in affected pigs; however, NK cells and neutrophils were present in numbers comparable to unaffected littermates. The results indicate that the immune defect in pigs shares the same features as other SCID-affected species.

Integrating comparative expression profiling data and association of SNPs with Salmonella shedding for improved food safety and porcine disease resistance.

Salmonella in swine is a major food safety problem, as the majority of US swine herds are Salmonella-positive. Salmonella can be shed from colonized swine and contaminate (i) neighbouring pigs; (ii) slaughter plants and pork products; (iii) edible crops when swine manure is used as a fertilizer; and (iv) water supplies if manure used as crop fertilizer runs off into streams and waterways. A potentially powerful method of addressing pre-harvest food safety at the farm level is through genetic improvement of disease resistance in animals. In this research, we describe a successful strategy for discovering genetic variation at candidate genes associated with disease resistance in pigs. This involves integrating our recent global gene expression analysis of the porcine response to Salmonella with information from the literature about important candidate genes. We identified single-nucleotide polymorphisms (SNPs) in these functional candidate genes and genotyped three independent pig populations that had data on Salmonella faecal shedding or internal burden (total n = 377) at these loci. Of 31 SNPs genotyped, 21 SNPs segregated in at least two populations with a minor allele frequency of 15% or greater. Statistical analysis revealed thirteen SNPs associated with Salmonella faecal shedding or tissue colonization, with an estimated proportion of false positives (PFP) ≤0.2. The genes with associated SNPs included GNG3, NCF2, TAP1, VCL, AMT, CCR1, CD163, CCT7, EMP1 and ACP2. These associations provide new information about the mechanisms of porcine host response to Salmonella and may be useful in improving genetic resistance to this bacterium.

Use of bioinformatic SNP predictions in differentially expressed genes to find SNPs associated with Salmonella colonization in swine.

Asymptomatic Salmonella-carrier pigs present a major problem in preharvest food safety, with a recent survey indicating >50% of swine herds in the United States have Salmonella-positive animals. Salmonella-carrier pigs serve as a reservoir for contamination of neighbouring pigs, abattoir pens and pork products. In addition, fresh produce as well as water can be contaminated with Salmonella from manure used as fertilizer. Control of Salmonella at the farm level could be through genetic improvement of porcine disease resistance, a potentially powerful method of addressing preharvest pork safety. In this research, we integrate gene expression profiling data and sequence alignment-based prediction of single nucleotide polymorphisms (SNPs) to successfully identify SNPs in functional candidate genes to test for the associations with swine response to Salmonella. A list of 2527 genes that were differentially regulated in porcine whole blood in response to infection with Salmonella enterica serovar Typhimurium were selected. In those genes, SNPs were predicted using ANEXdb alignments based on stringent clustering of all publically available porcine cDNA and expressed sequence tag (EST) sequences. A set of 30 mostly non-synonymous SNPs were selected for genotype analysis of four independent populations (n = 750) with Salmonella faecal shedding or tissue colonization phenotypes. Nine SNPs segregated with minor allele frequency ≥15% in at least two populations. Statistical analysis revealed SNPs associated with Salmonella shedding, such as haptoglobin (HP, p = 0.001, q = 0.01), neutrophil cytosolic factor 2 (NCF2 #2, p = 0.04, q = 0.21) and phosphogluconate dehydrogenase (p = 0.066, q = 0.21). These associations may be useful in identifying and selecting pigs with improved resistance to this bacterium.

Methods for transcriptomic analyses of the porcine host immune response: application to Salmonella infection using microarrays.

Technological developments in both the collection and analysis of molecular genetic data over the past few years have provided new opportunities for an improved understanding of the global response to pathogen exposure. Such developments are particularly dramatic for scientists studying the pig, where tools to measure the expression of tens of thousands of transcripts, as well as unprecedented data on the porcine genome sequence, have combined to expand our abilities to elucidate the porcine immune system. In this review, we describe these recent developments in the context of our work using primarily microarrays to explore gene expression changes during infection of pigs by Salmonella. Thus while the focus is not a comprehensive review of all possible approaches, we provide links and information on both the tools we use as well as alternatives commonly available for transcriptomic data collection and analysis of porcine immune responses. Through this review, we expect readers will gain an appreciation for the necessary steps to plan, conduct, analyze and interpret the data from transcriptomic analyses directly applicable to their research interests.

Correlating blood immune parameters and a CCT7 genetic variant with the shedding of Salmonella enterica serovar Typhimurium in swine.

The porcine response to Salmonella infection is critical for control of Salmonella fecal shedding and the establishment of Salmonella carrier status. In this study, 40 crossbred pigs were intranasally inoculated with Salmonella enterica serovar Typhimurium (Salmonella Typhimurium) and monitored for Salmonella fecal shedding and blood immune parameters at 2, 7, 14 and 20 days post-inoculation (dpi). Using a multivariate permutation test, a positive correlation was observed between Salmonella Typhimurium shedding levels at 2 and 7dpi and serum interferon-gamma (IFNgamma) levels at 2dpi (p<0.05), with Salmonella being shed in greater numbers from animals with higher IFNgamma levels. A positive correlation was also observed between IFNgamma levels and the number of banded neutrophils (2dpi), circulating neutrophils (7 and 14dpi), monocytes (7dpi), and white blood cells (WBCs) (7, 14 and 20dpi). We have further performed association studies on these immune response parameters as well as shedding status of the Salmonella-infected pigs with a single nucleotide polymorphism (SNP) in the porcine gene CCT7, previously shown by our group to be transcriptionally up-regulated in swine experimentally inoculated with Salmonella Typhimurium. Our analyses with the 40 pigs suggest a positive association (p=0.0012) of SNP genotype A/G at position AK240296.c1153G>A of the CCT7 gene with Salmonella shedding at 7dpi compared to the G/G homozygote genotype. Linking specific genes and genetic polymorphisms with the porcine immune response to Salmonella infection and shedding may identify potential markers for carrier pigs as well as targets for disease diagnosis, intervention and prevention.

Computational integration of structural and functional genomics data across species to develop information on the porcine inflammatory gene regulatory pathway.

We are investigating the porcine gut immune response to infection through gene expression profiling. Porcine Affymetrix GeneChip data was obtained from RNA prepared from mesenteric lymph node of swine infected with either Salmonella enterica serovar Typhimurium (ST) or S. Choleraesuis (SC) for 0, 8, 24, 48 or 504 hours post-inoculation (hpi). In total, 2365 genes with statistical evidence for differential expression (DE; p < 0.01, q < 0.26, fold-change > 2) between at least two time-points were identified. Comparative Gene Ontology analyses revealed that a high proportion of annotated DE genes in both infections are involved in immune and defence responses. Hierarchical clustering of expression patterns and annotations showed that 22 of the 83 genes upregulated from 8-24 hpi in the SC infection are known NF-kappaB targets. The promoter sequences of human genes orthologous to the DE genes were collected and TFM-Explorer was used to identify a set of 72 gene promoters with significant over-representation of NF-kappaB DNA-binding motifs. All 22 known NF-kappaB target genes are in this list; we hypothesize that the remaining 51 genes are un-recognized NF-kappaB targets. Integration of these results and verification of putative target genes will increase our understanding of the porcine response pathways responding to bacterial infection.

Integration of structural and functional genomics.

This paper introduces a special issue of Animal Genetics, which is devoted to the recent symposium held at Iowa State University entitled 'Integration of Structural and Functional Genomics'. We describe issues and needs that confront the animal genomics community, and describe how this symposium was structured to address these issues by improving communication and collaboration across species and disciplines. The session topics and oral presentations are briefly described for all invited speakers.

Mapping of 443 porcine EST improves the comparative maps for SSC1 and SSC7 with the human genome.

Numerous mapping studies of complex traits in the pig have resulted in quantitative trait loci (QTL) intervals of 10-20 cM. To improve the chances to identify the genes located in such intervals, increased expressed sequence tags (EST)-based marker density, coupled with comparative mapping with species whose genomes have been sequenced such as human and mouse, is the most efficient tool. In this study, we mapped 443 porcine EST with a radiation hybrid (RH) panel (384 had LOD > 6.0) and a somatic cell hybrid panel. Requiring no discrepancy between two-point and multipoint RH data allowed robust assignment of 309 EST, of which most were located on porcine chromosomes (SSC) 1, 4, 7, 8 and X. Moreover, we built framework maps for two chromosomes, SSC1 and SSC7, with mapped QTL in regions with known rearrangement between pig and human genomes. Using the Blast tool, we found orthologies between 407 of the 443 pig cDNA sequences and human genes, or to existing pig genes. Our porcine/human comparative mapping results reveal possible new homologies for SSC1, SSC3, SSC5, SSC6, SSC12 and SSC14 and add markers in synteny breakpoints for chromosome 7.

Physical mapping of genes in the porcine ovarian transcriptome.

Reproductive efficiency and associated traits are of major economic importance to the swine industry and have been more difficult to improve genetically than other production traits. Integration of phenotypical data with gene mapping and expression studies provides a powerful approach for dissection of the genetic basis regulating complex traits. We developed a total of 101 polymerase chain reaction-based markers, representing 91 unique genes, for expressed sequence tags previously reported to be putatively differentially expressed in the porcine ovarian transcriptome of a swine line selected on an index of high ovulation rate and embryonic survival. These were subsequently used in physical mapping experiments with a porcine radiation hybrid and somatic cell hybrid panels. Our results increased the information content of the porcine physical map useful for comparative mapping by c. 10%. Moreover, the mapped genes are likely to be biologically relevant to the molecular mechanisms that control ovulation rate in the pig. A total of 12 differentially expressed genes were mapped to regions previously reported to contain quantitative trait loci affecting swine ovulation rate.