Oliver P. Couture

Gene expression profiling of the short-term adaptive response to acute caloric restriction in liver and adipose tissues of pigs differing in feed efficiency

Residual feed intake (RFI) is a measure of feed efficiency, in which low RFI denotes improved feed efficiency. Caloric restriction (CR) is associated with feed efficiency in livestock species and to human health benefits, such as longevity and cancer prevention. We have developed pig lines that differ in RFI, and we are interested in identifying the genes and pathways that underlie feed efficiency. Prepubertal Yorkshire gilts with low RFI (n = 10) or high RFI (n = 10) were fed ad libitum or fed at restricted intake of 80% of maintenance energy requirements for 8 days. We measured serum metabolites and hormones and generated transcriptional profiles of liver and subcutaneous adipose tissue on these animals. Overall, 6,114 genes in fat and 305 genes in liver were differentially expressed (DE) in response to CR, and 311 genes in fat and 147 genes in liver were DE due to RFI differences. Pathway analyses of CR-induced DE genes indicated a dramatic switch to a conservation mode of energy usage by down-regulating lipogenesis and steroidogenesis in both liver and fat. Interestingly, CR altered expression of genes in immune and cell cycle/apoptotic pathways in fat, which may explain part of the CR-driven lifespan enhancement. In silico analysis of transcription factors revealed ESR1 as a putative regulator of the adaptive response to CR, as several targets of ESR1 in our DE fat genes were annotated as cell cycle/apoptosis genes. The lipid metabolic pathway was overrepresented by down-regulated genes due to both CR and low RFI. We propose a common energy conservation mechanism, which may be controlled by PPARA, PPARG, and/or CREB in both CR and feed-efficient pigs.

Identification of differential gene expression during porcine conceptus rapid trophoblastic elongation and attachment to uterine luminal epithelium

Early embryonic development in the pig is characterized by a rapid elongation of the conceptus trophectoderm on days 11-12 of gestation. Initially, the conceptus trophoblast is morphologically rearranged from a 10-mm sphere into a tubular shape, transitioning into a thin filamentous form >150 mm in length in 2-3 h, followed by continued expansion within the uterine lumen for several days. Conceptus elongation is critical for establishing adequate placental surface area needed for embryo and fetal survival throughout gestation. The objective of this study was to characterize conceptus gene expression during trophoblastic elongation and the early attachment to the uterine endometrium on days 11-14 of gestation with the GeneChip Porcine Genome Array. In all, 3,759 different probe sets were statistically different in at least one comparison [spherical vs. tubular, spherical vs. day 12 filamentous (D12F), spherical vs. day 14 filamentous (D14F), tubular vs. D12F, tubular vs. D14F, and D12F vs. D14F]. When restricted to the spherical vs. D12F and D12F vs. D14F comparisons, 482 and 232 genes, respectively, were statistically different with greater than twofold change in expression. Utilization of k-means clustering, in addition to the Database for Annotation, Visualization, and Integrated Discovery (DAVID), identified genes of interest. Quantitative RT-PCR expression profiles for interferon-gamma (IFNG), heat shock protein 27 kDa (HSPB1), angiomotin, B-cell linker (BLNK), chemokine ligand 14 (CXCL14), parathyroid hormone-like hormone (PTHLH), and maspin were supportive of the GeneChip Porcine Genome Array data.

Microarray gene expression profiles of fasting induced changes in liver and adipose tissues of pigs expressing the melanocortin-4 receptor D298N variant

Transcriptional profiling coupled with blood metabolite analyses were used to identify porcine genes and pathways that respond to a fasting treatment or to a D298N missense mutation in the melanocortin-4 receptor (MC4R) gene. Gilts (12 homozygous for D298 and 12 homozygous for N298) were either fed ad libitum or fasted for 3 days. Fasting decreased body weight, backfat, and serum urea concentration and increased serum nonesterified fatty acid. In response to fasting, 7,029 genes in fat and 1,831 genes in liver were differentially expressed (DE). MC4R genotype did not significantly affect gene expression, body weight, backfat depth, or any measured serum metabolite concentration. Pathway analyses of fasting-induced DE genes indicated that lipid and steroid synthesis was downregulated in both liver and fat. Fasting increased expression of genes involved in glucose sparing pathways, such as oxidation of amino acids and fatty acids in liver, and in extracellular matrix pathways, such as cell adhesion and adherens junction in fat. Additionally, we identified DE transcription factors (TF) that regulate many DE genes. This confirms the involvement of TF, such as PPARG, SREBF1, and CEBPA, which are known to regulate the fasting response, and implicates additional TF, such as ESR1. Interestingly, ESR1 controls several fasting induced genes in fat that are involved in cell matrix morphogenesis. Our findings indicate a transcriptional response to fasting in two key metabolic tissues of pigs, which was corroborated by changes in blood metabolites, and the involvement of novel putative transcriptional regulators in the immediate adaptive response to fasting.

ANEXdb: an integrated animal ANnotation and microarray EXpression database

To determine annotations of the sequence elements on microarrays used for transcriptional profiling experiments in livestock species, currently researchers must either use the sparse direct annotations available for these species or create their own annotations. ANEXdb (http://www.anexdb.org) is an open-source web application that supports integrated access of two databases that house microarray expression (ExpressDB) and EST annotation (AnnotDB) data. The expression database currently supports storage and querying of Affymetrix-based expression data as well as retrieval of experiments in a form ready for NCBI-GEO submission; these services are available online. AnnotDB currently houses a novel assembly of approximately 1.6 million unique porcine-expressed sequence reads called the Iowa Porcine Assembly (IPA), which consists of 140,087 consensus sequences, the Iowa Tentative Consensus (ITC) sequences, and 103,888 singletons. The IPA has been annotated via transfer of information from homologs identified through sequence alignment to NCBI RefSeq. These annotated sequences have been mapped to the Affymetrix porcine array elements, providing annotation for 22,569 of the 23,937 (94%) porcine-specific probe sets, of which 19,253 (80%) are linked to an NCBI RefSeq entry. The ITC has also been mined for sequence variation, providing evidence for up to 202,383 SNPs, 62,048 deletions, and 958 insertions in porcine-expressed sequence. These results create a single location to obtain porcine annotation of and sequence variation in differently expressed genes in expression experiments, thus permitting possible identification of causal variants in such genes of interest. The ANEXdb application is open source and available from SourceForge.net.

Expression of collagen genes in the cones of skin in the Duroc/Yorkshire porcine model of fibroproliferative scarring.

During the past decades there has been minimal improvement in prevention and treatment of hypertrophic scarring. Reasons include the lack of a validated animal model, imprecise techniques to dissect scar into the histologic components, and limited methodology for measurement of gene expression. These problems have been addressed with the Duroc/Yorkshire model of healing, laser capture microdissection, and the Affymetrix Porcine GeneChip®. Here we compared collagen gene expression in fibroproliferative healing in the Duroc breed to nonfibroproliferative healing in the Yorkshires. We made shallow and deep dorsal wounds, biopsied at 1, 2, 3, 12, and 20 weeks. We sampled the dermal cones by laser capture microdissection, extracted and amplified the RNA, and hybridized Affymetrix Porcine GeneChips®. We also obtained samples of human hypertrophic scar approximately 20 weeks postinjury. Data were normalized and statistical analysis performed with mixed linear regression using the Bioconductor R/maanova package. Genes for further analysis were also restricted with four biologic criteria, including that the 20-week deep Duroc expression match the human samples. Eleven ollagen genes and seven collagen types were differentially over expressed in deep Duroc wounds including 1a1, 1a2, 3a1, 4a1, 4a2, 5a1, 5a2, 5a3, 6a3 (transcript variant 5), 14a1 and 15a1. COL7a1 gene was differentially under expressed in deep Duroc wounds. The results suggest that collagens I, III, IV, V, VI, VII, XIV, and XV1 are involved in the process of fibroproliferative scarring. With these clues, we will attempt to construct the regulatory pathway(s) of fibroproliferative healing.

Functional Genomics Unique to Week 20 Post Wounding in the Deep Cone/Fat Dome of the Duroc/Yorkshire Porcine Model of Fibroproliferative Scarring

Background Hypertrophic scar was first described over 100 years ago; PubMed has more than 1,000 references on the topic. Nevertheless prevention and treatment remains poor, because 1) there has been no validated animal model; 2) human scar tissue, which is impossible to obtain in a controlled manner, has been the only source for study; 3) tissues typically have been homogenized, mixing cell populations; and 4) gene-by-gene studies are incomplete. Methodology/Principal Findings We have assembled a system that overcomes these barriers and permits the study of genome-wide gene expression in microanatomical locations, in shallow and deep partial-thickness wounds, and pigmented and non-pigmented skin, using the Duroc(pigmented fibroproliferative)/Yorkshire(non-pigmented non-fibroproliferative) porcine model. We used this system to obtain the differential transcriptome at 1, 2, 3, 12 and 20 weeks post wounding. It is not clear when fibroproliferation begins, but it is fully developed in humans and the Duroc breed at 20 weeks. Therefore we obtained the derivative functional genomics unique to 20 weeks post wounding. We also obtained long-term, forty-six week follow-up with the model. Conclusions/Significance 1) The scars are still thick at forty-six weeks post wounding further validating the model. 2) The differential transcriptome provides new insights into the fibroproliferative process as several genes thought fundamental to fibroproliferation are absent and others differentially expressed are newly implicated. 3) The findings in the derivative functional genomics support old concepts, which further validates the model, and suggests new avenues for reductionist exploration. In the future, these findings will be searched for directed networks likely involved in cutaneous fibroproliferation. These clues may lead to a better understanding of the systems biology of cutaneous fibroproliferation, and ultimately prevention and treatment of hypertrophic scarring.

Gene expression in hypothalamus, liver, and adipose tissues and food intake response to melanocortin-4 receptor agonist in pigs expressing melanocortin-4 receptor mutations.

Transcriptional profiling was used to identify genes and pathways that responded to intracerebroventricular injection of melanocortin-4 receptor (MC4R) agonist [Nle(4), d-Phe(7)]-α-melanocyte stimulating hormone (NDP-MSH) in pigs homozygous for the missense mutation in the MC4R, D298 allele (n = 12), N298 allele (n = 12), or heterozygous (n = 12). Food intake (FI) was measured at 12 and 24 h after treatment. All pigs were killed at 24 h after treatment, and hypothalamus, liver, and back-fat tissue was collected. NDP-MSH suppressed (P < 0.004) FI at 12 and 24 h in all animals after treatment. In response to NDP-MSH, 278 genes in hypothalamus (q ≤ 0.07, P ≤ 0.001), 249 genes in liver (q ≤ 0.07, P ≤ 0.001), and 5,066 genes in fat (q ≤ 0.07, P ≤ 0.015) were differentially expressed. Pathway analysis of NDP-MSH-induced differentially expressed genes indicated that genes involved in cell communication, nucleotide metabolism, and signal transduction were prominently downregulated in the hypothalamus. In both liver and adipose tissue, energy-intensive biosynthetic and catabolic processes were downregulated in response to NDP-MSH. This included genes encoding for biosynthetic pathways such as steroid and lipid biosynthesis, fatty acid synthesis, and amino acid synthesis. Genes involved in direct energy-generating processes, such as oxidative phosphorylation, electron transport, and ATP synthesis, were upregulated, whereas TCA-associated genes were prominently downregulated in NDP-MSH-treated pigs. Our data also indicate a metabolic switch toward energy conservation since genes involved in energy-intensive biosynthetic and catabolic processes were downregulated in NDP-MSH-treated pigs.

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.

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.