Long Qu

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.

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.

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.

A Fast EM Algorithm for BayesA-Like Prediction of Genomic Breeding Values

Prediction accuracies of estimated breeding values for economically important traits are expected to benefit from genomic information. Single nucleotide polymorphism (SNP) panels used in genomic prediction are increasing in density, but the Markov Chain Monte Carlo (MCMC) estimation of SNP effects can be quite time consuming or slow to converge when a large number of SNPs are fitted simultaneously in a linear mixed model. Here we present an EM algorithm (termed “fastBayesA”) without MCMC. This fastBayesA approach treats the variances of SNP effects as missing data and uses a joint posterior mode of effects compared to the commonly used BayesA which bases predictions on posterior means of effects. In each EM iteration, SNP effects are predicted as a linear combination of best linear unbiased predictions of breeding values from a mixed linear animal model that incorporates a weighted marker-based realized relationship matrix. Method fastBayesA converges after a few iterations to a joint posterior mode of SNP effects under the BayesA model. When applied to simulated quantitative traits with a range of genetic architectures, fastBayesA is shown to predict GEBV as accurately as BayesA but with less computing effort per SNP than BayesA. Method fastBayesA can be used as a computationally efficient substitute for BayesA, especially when an increasing number of markers bring unreasonable computational burden or slow convergence to MCMC approaches.

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.

A fast and efficient Gibbs sampler for BayesB in whole-genome analyses

Background In whole-genome analyses, the number p of marker covariates is often much larger than the number n of observations. Bayesian multiple regression models are widely used in genomic selection to address this problem of

Improved Estimation of the Noncentrality Parameter Distribution from a Large Number of t-Statistics, with Applications to False Discovery Rate Estimation in Microarray Data Analysis

p\gg n.