Oligonucleotide microarray analysis of age-related gene expression profiles in miniature pigs

Increasing gene discovery and coverage using RNA-seq of globin RNA reduced porcine blood samples

Background

Transcriptome analysis of porcine whole blood has several applications, which include deciphering genetic mechanisms for host responses to viral infection and vaccination. The abundance of alpha- and beta-globin transcripts in blood, however, impedes the ability to cost-effectively detect transcripts of low abundance. Although protocols exist for reduction of globin transcripts from human and mouse/rat blood, preliminary work demonstrated these are not useful for porcine blood Globin Reduction (GR). Our objectives were to develop a porcine specific GR protocol and to evaluate the GR effects on gene discovery and sequence read coverage in RNA-sequencing (RNA-seq) experiments.

Results

A GR protocol for porcine blood samples was developed using RNase H with antisense oligonucleotides specifically targeting porcine hemoglobin alpha (HBA) and beta (HBB) mRNAs. Whole blood samples (n = 12) collected in Tempus tubes were used for evaluating the efficacy and effects of GR on RNA-seq. The HBA and HBB mRNA transcripts comprised an average of 46.1% of the mapped reads in pre-GR samples, but those reads reduced to an average of 8.9% in post-GR samples. Differential gene expression analysis showed that the expression level of 11,046 genes were increased, whereas 34 genes, excluding HBA and HBB, showed decreased expression after GR (FDR <0.05). An additional 815 genes were detected only in post-GR samples.

Conclusions

Our porcine specific GR primers and protocol minimize the number of reads of globin transcripts in whole blood samples and provides increased coverage as well as accuracy and reproducibility of transcriptome analysis. Increased detection of low abundance mRNAs will ensure that studies relying on transcriptome analyses do not miss information that may be vital to the success of the study.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-954) contains supplementary material, which is available to authorized users.

Genome-wide DNA methylation changes in skeletal muscle between young and middle-aged pigs

BACKGROUND

Age-related physiological, biochemical and functional changes in mammalian skeletal muscle have been shown to begin at the mid-point of the lifespan. However, the underlying changes in DNA methylation that occur during this turning point of the muscle aging process have not been clarified. To explore age-related genomic methylation changes in skeletal muscle, we employed young (0.5 years old) and middle-aged (7 years old) pigs as models to survey genome-wide DNA methylation in the longissimus dorsi muscle using a methylated DNA immunoprecipitation sequencing approach.

RESULTS

We observed a tendency toward a global loss of DNA methylation in the gene-body region of the skeletal muscle of the middle-aged pigs compared with the young group. We determined the genome-wide gene expression pattern in the longissimus dorsi muscle using microarray analysis and performed a correlation analysis using DMR (differentially methylated region)-mRNA pairs, and we found a significant negative correlation between the changes in methylation levels within gene bodies and gene expression. Furthermore, we identified numerous genes that show age-related methylation changes that are potentially involved in the aging process. The methylation status of these genes was confirmed using bisulfite sequencing PCR. The genes that exhibited a hypomethylated gene body in middle-aged pigs were over-represented in various proteolysis and protein catabolic processes, suggesting an important role for these genes in age-related muscle atrophy. In addition, genes associated with tumorigenesis exhibited aged-related differences in methylation and expression levels, suggesting an increased risk of disease associated with increased age.

CONCLUSIONS

This study provides a comprehensive analysis of genome-wide DNA methylation patterns in aging pig skeletal muscle. Our findings will serve as a valuable resource in aging studies, promoting the pig as a model organism for human aging research and accelerating the development of comparative animal models in aging research.

Differential gene expression profiling of porcine epithelial cells infected with three enterotoxigenic Escherichia coli strains

Background

Enterotoxigenic Escherichia coli (ETEC) is one of the most important pathogenic bacteria causing severe diarrhoea in human and pigs. In ETEC strains, the fimbrial types F4 and F18 are commonly found differently colonized within the small intestine and cause huge economic losses in the swine industry annually worldwide. To address the underlying mechanism, we performed a transcriptome study of porcine intestinal epithelial cells (IPEC-J2) with and without infection of three representative ETEC strains.

Results

A total 2443, 3493 and 867 differentially expressed genes were found in IPEC-J2 cells infected with F4ab ETEC (C_F4ab), with F4ac ETEC (C_F4ac) and with F18ac ETEC (C_F18ac) compared to the cells without infection (control), respectively. The number of differentially expressed genes between C_F4ab and C_F4ac, C_F4ab and C_F18ac, and C_F4ac and C_F18ac were 77, 1446 and 1629, respectively. The gene ontology and pathway analysis showed that the differentially expressed genes in C_F4abvs control are significantly involved in cell-cycle progress and amino acid metabolism, while the clustered terms of the differentially expressed genes in C_F4acvs control comprise immune, inflammation and wounding response and apoptosis as well as cell cycle progress and proteolysis. Differentially expressed genes between C_F18acvs control are mainly involved in cell-cycle progression and immune response. Furthermore, fundamental differences were observed in expression levels of immune-related genes among the three ETEC treatments, especially for the important pro-inflammatory molecules, including IL-6, IL-8, TNF-α, CCL20, CXCL2 etc.

Conclusions

The discovery in this study provides insights into the interaction of porcine intestinal epithelial cells with F4 ETECs and F18 ETEC, respectively. The genes induced by ETECs with F4 versus F18 fimbriae suggest why ETEC with F4 may be more virulent compared to F18 which seems to elicit milder effects.

Oligonucleotide microarray analysis of dietary-induced hyperlipidemia gene expression profiles in miniature pigs

Background

Hyperlipidemia animal models have been established, but complete gene expression profiles of the transition from normal lipid levels have not been obtained. Miniature pigs are useful model animals for gene expression studies on dietary-induced hyperlipidemia because they have a similar anatomy and digestive physiology to humans, and blood samples can be obtained from them repeatedly.

Methodology

Two typical dietary treatments were used for dietary-induced hyperlipidemia models, by using specific pathogen-free (SPF) Clawn miniature pigs. One was a high-fat and high-cholesterol diet (HFCD) and the other was a high-fat, high-cholesterol, and high-sucrose diet (HFCSD). Microarray analyses were conducted from whole blood samples during the dietary period and from white blood cells at the end of the dietary period to evaluate the transition of expression profiles of the two dietary models.

Principal Findings

Variations in whole blood gene expression intensity within the HFCD or the HFCSD group were in the same range as the controls provide with normal diet at all periods. This indicates uniformity of dietary-induced hyperlipidemia for our dietary protocols. Gene ontology- (GO) based functional analyses revealed that characteristics of the common changes between HFCD and HFCSD were involved in inflammatory responses and reproduction. The correlation coefficient between whole blood and white blood cell expression profiles at 27 weeks with the HFCSD diet was significantly lower than that of the control and HFCD diet groups. This may be due to the effects of RNA originating from the tissues and/or organs.

Conclusions

No statistically significant differences in fasting plasma lipids and glucose levels between the HFCD and HFCSD groups were observed. However, blood RNA analyses revealed different characteristics corresponding to the dietary protocols. In this study, whole blood RNA analyses proved to be a useful tool to evaluate transitions in dietary-induced hyperlipidemia gene expression profiles in miniature pigs.