ANEXdb: an integrated animal ANnotation and microarray EXpression database

Comparison of organelle genomes between endangered mangrove plant Dolichandrone spathacea to terrestrial relative provides insights into its origin and adaptative evolution

Dolichandrone spathacea is a mangrove associate with high medicinal and ecological values. However, due to the dual-pressure of climate change and human activities, D. spathacea has become endangered in China. Moreover, misidentification between D. spathacea and its terrestrial relative D. cauda-felina poses further challenges to field protection and proper medicinal usage of D. spathacea. Thus, to address these problems, we sequenced and assembled mitochondrial (mt) and chloroplast (cp) genomes for both D. spathacea and D. cauda-felina. Comparative analysis revealed apparently different size and scaffold number between the two mt genomes, but a high similarity between the cp genomes. Eight regions with high sequence divergence were identified between the two cp genomes, which might be used for developing candidate DNA markers for distinguishing the two species. The splitting between D. spathacea and D. cauda-felina was inferred to occur at ~6.8 - 7.7 million years ago (Mya), which may be driven by the environment fluctuations in late Miocene. In the cp genome, 12 genes related to the expression of photosynthesis-associated proteins were detected with signatures of positive selection, which may contribute to the origin and evolutionary adaptation of Dolichandrone mangrove species. These new findings do not only enrich organelle genomic resources of Dolichandrone species, but also provide important genetic clues for improving the conservation and proper usage of endangered mangrove associate D. spathacea.

Re-Expression of Tafazzin Isoforms in TAZ-Deficient C6 Glioma Cells Restores Cardiolipin Composition but Not Proliferation Rate and Alterations in Gene Expression

Tafazzin—an acyltransferase—is involved in cardiolipin (CL) remodeling. CL is associated with mitochondrial function, structure and more recently with cell proliferation. Various tafazzin isoforms exist in humans. The role of these isoforms in cardiolipin remodeling is unknown. Aim of this study was to investigate if specific isoforms like Δ5 can restore the wild type phenotype with respect to CL composition, cellular proliferation and gene expression profile. In addition, we aimed to determine the molecular mechanism by which tafazzin can modulate gene expression by applying promoter analysis and (Ingenuity Pathway Analyis) IPA to genes regulated by TAZ-deficiency. Expression of Δ5 and rat full length TAZ in C6-TAZ- cells could fully restore CL composition and—as proven for Δ5—this is naturally associated with restoration of mitochondrial respiration. A similar restoration of CL-composition could not be observed after re-expression of an enzymatically dead full-length rat TAZ (H69L; TAZMut). Re-expression of only rat full length TAZ could restore proliferation rate. Surprisingly, the Δ5 variant failed to restore wild-type proliferation. Further, as expected, re-expression of the TAZMut variant completely failed to reverse the gene expression changes, whereas re-expression of the TAZ-FL variant largely did so and the Δ5 variant to somewhat less extent. Very likely TAZ-deficiency provokes substantial long-lasting changes in cellular lipid metabolism which contribute to changes in proliferation and gene expression, and are not or only very slowly reversible.

Integration of "omics" techniques: Dronedarone affects cardiac remodeling in the infarction border zone

Dronedarone improves microvascular flow during atrial fibrillation and reduces the infarct size in acute models of myocardial infarction. However, dronedarone might be harmful in patients with recent decompensated heart failure and increases mortality in patients with permanent atrial fibrillation. A pathophysiological explanation for these discrepant data is lacking. This study investigated the effects of dronedarone on gene and protein expression in the infarcted area and border zone in pigs subjected to anterior ischemia/reperfusion myocardial infarction. The ischemia/reperfusion myocardial infarction was induced in 16 pigs. Eight pigs were treated with dronedarone for 28 days after myocardial infarction, the remaining pigs served as control. Microarray-based transcriptome profiling and 2D-DIGE-based proteome analysis were used to assess the effects of dronedarone on left ventricular gene expression in healthy (LV), infarcted (MI), and border zone tissue. Selected targets were validated by RT-qPCR or immunoblot analyses, with special emphasize given to the transcriptome/proteome overlap. Combined "omics" analysis was performed to identify most significant disease and function charts affected by dronedarone and to establish an integrated network. The levels of 879 (BZ) or 7 (MI) transcripts and 51 (LV) or 15 (BZ) proteins were significantly altered by dronedarone, pointing to a substantial efficacy of dronedarone in the border zone. Transcriptome and proteome data indicate that dronedarone influences post-infarction remodeling processes and identify matricellular proteins as major targets of dronedarone in this setting. This finding is fully supported by the disease and function charts as well as by the integrated network established by combined "omics". Dronedarone therapy alters myocardial gene expression after acute myocardial infarction with pronounced effects in the border zone. Dronedarone promotes infarct healing via regulation of periostin and might contribute to the limitation of its expansion as well as cardiac rupture. Thus, there are no experimental hints that dronedarone per se has direct harmful effects after MI in ventricular tissue. Impact statement Dronedarone reduced the infarct size in models of acute myocardial infarction (MI). Here, we show that dronedarone attenuates many of the substantial changes in gene expression that are provoked by acute myocardial infarction (AMI) in pigs. Dronedarone modifies the expression of gene panels related to post-infarction cardiac healing and remodeling processes and, most remarkably, this occurs predominantly in the infarction border-zone and much less so in the vital or infarcted myocardium. Combined "omics" identified matricellular proteins and ECM as major dronedarone-regulated targets and emphasizes their relevance for Disease Charts and Tox Function Charts associated with tissue remodeling and cellular movement. The results demonstrate dronedarone's capability of regulating cardiac repair and remodeling processes specifically in the infarction border zone and identify underlying mechanisms and pathways that might be employed in future therapeutic strategies to improve long-term cardiac tissue function and stability.

Atrial fibrillation and rapid acute pacing regulate adipocyte/adipositas-related gene expression in the atria

PURPOSE

Atrial fibrillation (AF) has been associated with increased volumes of epicardial fat and atrial adipocyte accumulation. Underlying mechanisms are not well understood. This study aims to identify rapid atrial pacing (RAP)/AF-dependent changes in atrial adipocyte/adipositas-related gene expression (AARE).

METHODS

Right atrial (RA) and adjacent epicardial adipose tissue (EAT) samples were obtained from 26 patients; 13 with AF, 13 in sinus rhythm (SR). Left atrial (LA) samples were obtained from 9 pigs (5 RAP, 4 sham-operated controls). AARE was analyzed using microarrays and RT-qPCR. The impact of diabetes/obesity on gene expression was additionally determined in RA samples (RAP ex vivo and controls) from 3 vs. 6 months old ZDF rats.

RESULTS

RAP in vivo of pigs resulted in substantial changes of AARE, with 66 genes being up- and 53 down-regulated on the mRNA level. Differential expression during adipocyte differentiation was confirmed using 3T3-L1 cells. In patients with AF (compared to SR), a comparable change in RA mRNA levels concerned a fraction of genes only (RETN, IGF1, HK2, PYGM, LOX, and NR4A3). RA and EAT were affected by AF to a different extent. In patients, concomitant disease contributes to AARE changes.

CONCLUSIONS

RAP, and to lesser extent AF, provoke significant changes in atrial AARE. In chronic AF, activation of this gene panel is very likely mediated by AF itself, AF risk factors and concomitant diseases. This may facilitate the development of an AF substrate by increasing atrial ectopic fat and fat infiltration of the atrial myocardium.

Endometrial gene expression profiling in pregnant Meishan and Yorkshire pigs on day 12 of gestation

Background

Litter size in pigs is a major factor affecting the profitability in the pig industry. The peri-implantation window in pigs is characterized by the coordinated interactions between the maternal uterine endometrium and the rapidly elongating conceptuses and represents a period of time during which a large percentage of the developing conceptuses are lost. However, the gene expression and regulatory networks in the endometrium contributing to the establishment of the maternal: placental interface remain poorly understood.

Results

We characterized the endometrial gene expression profile during the peri-implantation stage of development by comparing two breeds that demonstrate very different reproductive efficiencies. We employed the porcine Affymetrix GeneChip® to assay the transcriptomic profiles of genes expressed in the uterine endometrium obtained from Meishan and Yorkshire gilts (n = 4 for each breed) on day 12 of gestation (M12 and Y12, respectively). Total of 17,076 probesets were identified as "present" in at least two arrays. A mixed model-based statistical analysis predicted a total of 2,656 (q < 0.1) transcripts as differentially expressed between Meishan and Yorkshire pigs. Eighteen differentially expressed transcripts of interest were validated by quantitative real-time PCR. Gene ontology (GO) annotation revealed that the known functions of the differentially expressed genes were involved in a series of important biological processes relevant to early pregnancy establishment in the pig.

Conclusions

The results identified endometrial gene expression profiles of two breeds differing in litter size and identified candidate genes that are related to known physiological pathways related to reproductive prolificacy. These findings provide a deeper understanding of molecular pathways differing between two breeds at the critical peri-implantation stage of pregnancy, which can be utilized to better understand the events contributing to pregnancy establishment in the pig.

Transcriptome analysis of porcine M. semimembranosus divergent in intramuscular fat as a consequence of dietary protein restriction

Background

Intramuscular fat (IMF) content is positively correlated with aspects of pork palatability, including flavour, juiciness and overall acceptability. The ratio of energy to protein in the finishing diet of growing pigs can impact on IMF content with consequences for pork quality. The objective of this study was to compare gene expression profiles of Musculus semimembranosus (SM) of animals divergent for IMF as a consequence of protein dietary restriction in an isocaloric diet. The animal model was derived through the imposition of low or high protein diets during the finisher stage in Duroc gilts. RNA was extracted from post mortem SM tissue, processed and hybridised to Affymetrix porcine GeneChip® arrays.

Results

IMF content of SM muscle was increased on the low protein diet (3.60 ± 0.38% versus 1.92 ± 0.35%). Backfat depth was also greater in animals on the low protein diet, and average daily gain and feed conversion ratio were lower, but muscle depth, protein content and moisture content were not affected. A total of 542 annotated genes were differentially expressed (DE) between animals on low and high protein diets, with 351 down-regulated and 191 up-regulated on the low protein diet. Transcript differences were validated for a subset of DE genes by qPCR. Alterations in functions related to cell cycle, muscle growth, extracellular matrix organisation, collagen development, lipogenesis and lipolysis, were observed. Expression of adipokines including LEP, TNFα and HIF1α were increased and the hypoxic stress response was induced. Many of the identified transcriptomic responses have also been observed in genetic and fetal programming models of differential IMF accumulation, indicating they may be robust biological indicators of IMF content.

Conclusion

An extensive perturbation of overall energy metabolism in muscle occurs in response to protein restriction. A low protein diet can modulate IMF content of the SM by altering gene pathways involved in lipid biosynthesis and degradation; however this nutritional challenge negatively impacts protein synthesis pathways, with potential consequences for growth.

Vascular transcriptional alterations produced by juvenile obesity in Ossabaw swine

We adopted a transcriptome-wide microarray analysis approach to determine the extent to which vascular gene expression is altered as a result of juvenile obesity and identify obesity-responsive mRNAs. We examined transcriptional profiles in the left anterior descending coronary artery (LAD), perivascular fat adjacent to the LAD, and descending thoracic aorta between obese (n = 5) and lean (n = 6) juvenile Ossabaw pigs (age = 22 wk). Obesity was experimentally induced by feeding the animals a high-fat/high-fructose corn syrup/high-cholesterol diet for 16 wk. We found that expression of 189 vascular cell genes in the LAD and expression of 165 genes in the thoracic aorta were altered with juvenile obesity (false discovery rate ≤ 10%) with an overlap of only 28 genes between both arteries. Notably, a number of genes found to be markedly upregulated in the LAD of obese pigs are implicated in atherosclerosis, including ACP5, LYZ, CXCL14, APOE, PLA2G7, LGALS3, SPP1, ITGB2, CYBB, and P2RY12. Furthermore, pathway analysis revealed the induction of proinflammatory and pro-oxidant pathways with obesity primarily in the LAD. Gene expression in the LAD perivascular fat was minimally altered with juvenile obesity. Together, we provide new evidence that obesity produces artery-specific changes in pretranslational regulation with a clear upregulation of proatherogenic genes in the LAD. Our data may offer potential viable drug targets and mechanistic insights regarding the molecular precursors involved in the origins of overnutrition and obesity-associated vascular disease. In particular, our results suggest that the oxidized LDL/LOX-1/NF-κB signaling axis may be involved in the early initiation of a juvenile obesity-induced proatherogenic coronary artery phenotype.

Pig immune response to general stimulus and to porcine reproductive and respiratory syndrome virus infection: a meta-analysis approach

BACKGROUND

The availability of gene expression data that corresponds to pig immune response challenges provides compelling material for the understanding of the host immune system. Meta-analysis offers the opportunity to confirm and expand our knowledge by combining and studying at one time a vast set of independent studies creating large datasets with increased statistical power. In this study, we performed two meta-analyses of porcine transcriptomic data: i) scrutinized the global immune response to different challenges, and ii) determined the specific response to Porcine Reproductive and Respiratory Syndrome Virus (PRRSV) infection. To gain an in-depth knowledge of the pig response to PRRSV infection, we used an original approach comparing and eliminating the common genes from both meta-analyses in order to identify genes and pathways specifically involved in the PRRSV immune response. The software Pointillist was used to cope with the highly disparate data, circumventing the biases generated by the specific responses linked to single studies. Next, we used the Ingenuity Pathways Analysis (IPA) software to survey the canonical pathways, biological functions and transcription factors found to be significantly involved in the pig immune response. We used 779 chips corresponding to 29 datasets for the pig global immune response and 279 chips obtained from 6 datasets for the pig response to PRRSV infection, respectively.

RESULTS

The pig global immune response analysis showed interconnected canonical pathways involved in the regulation of translation and mitochondrial energy metabolism. Biological functions revealed in this meta-analysis were centred around translation regulation, which included protein synthesis, RNA-post transcriptional gene expression and cellular growth and proliferation. Furthermore, the oxidative phosphorylation and mitochondria dysfunctions, associated with stress signalling, were highly regulated. Transcription factors such as MYCN, MYC and NFE2L2 were found in this analysis to be potentially involved in the regulation of the immune response. The host specific response to PRRSV infection engendered the activation of well-defined canonical pathways in response to pathogen challenge such as TREM1, toll-like receptor and hyper-cytokinemia/ hyper-chemokinemia signalling. Furthermore, this analysis brought forth the central role of the crosstalk between innate and adaptive immune response and the regulation of anti-inflammatory response. The most significant transcription factor potentially involved in this analysis was HMGB1, which is required for the innate recognition of viral nucleic acids. Other transcription factors like interferon regulatory factors IRF1, IRF3, IRF5 and IRF8 were also involved in the pig specific response to PRRSV infection.

CONCLUSIONS

This work reveals key genes, canonical pathways and biological functions involved in the pig global immune response to diverse challenges, including PRRSV infection. The powerful statistical approach led us to consolidate previous findings as well as to gain new insights into the pig immune response either to common stimuli or specifically to PRRSV infection.

Effects of irbesartan on gene expression revealed by transcriptome analysis of left atrial tissue in a porcine model of acute rapid pacing in vivo

BACKGROUND

Atrial fibrillation (AF) is characterized by electrical and structural remodeling of the atria with atrial fibrosis being one hallmark. Angiotensin II (AngII) is a major contributing factor and blockage of its type I receptor (AT1R) prevents remodeling to some extent. Here we explored the effects of the AT1R antagonist irbesartan on global gene expression and profibrotic signaling pathways after induction of rapid atrial pacing (RAP) in vivo in pigs.

METHODS AND RESULTS

Microarray-based RNA profiling was used to screen left atrial (LA) tissue specimens for differences in atrial gene expression in a model of acute RAP. RAP caused an overall expression profile that reflected AngII-induced ROS production, tissue remodeling, and energy depletion. Of special note, the mRNA levels of EDN1, SGK1, and CTGF encoding pro-endothelin, stress- and glucocorticoid activated kinase-1, and of connective tissue growth factor were identified to be significantly increased after 7h of rapid pacing. These specific expression changes were additionally validated by RT-qPCR or immunoblot analyses in LA, RA, and partly in LV samples. All RAP-induced differential gene expression patterns were partially attenuated in the presence of irbesartan. Similar results were obtained after RAP of HL-1 cardiomyocytes in vitro. Furthermore, exogenously added endothelin-1 (ET1) induced CTGF expression concomitant to the transcriptional activation of SGK1 in HL-1 cells.

CONCLUSIONS

RAP provokes substantial changes in atrial and ventricular myocardial gene expression that could be partly reversed by irbesartan. ET1 contributes to AF-dependent atrial fibrosis by synergistic activity with AngII to stimulate SGK1 expression and enhance phosphorylation of the SGK1 protein which, in turn, induces CTGF. The latter has been consistently associated with tissue fibrosis. These findings suggest ETR antagonists as being beneficial in AF treatment.

A hierarchical semiparametric model for incorporating intergene information for analysis of genomic data

For analysis of genomic data, e.g., microarray data from gene expression profiling experiments, the two-component mixture model has been widely used in practice to detect differentially expressed genes. However, it naïvely imposes strong exchangeability assumptions across genes and does not make active use of a priori information about intergene relationships that is currently available, e.g., gene annotations through the Gene Ontology (GO) project. We propose a general strategy that first generates a set of covariates that summarizes the intergene information and then extends the two-component mixture model into a hierarchical semiparametric model utilizing the generated covariates through latent nonparametric regression. Simulations and analysis of real microarray data show that our method can outperform the naïve two-component mixture model.

Dronedarone prevents microcirculatory abnormalities in the left ventricle during atrial tachypacing in pigs

BACKGROUND AND PURPOSE

Atrial fibrillation induces ischaemic microcirculatory flow abnormalities in the ventricle, contributing to the risk for acute coronary syndromes. We evaluated the effect of dronedarone on ventricular perfusion during rapid atrial pacing (RAP).

EXPERIMENTAL APPROACH

Coronary and fractional flow reserve (CFR/FFR) were measured in the left anterior descending artery in 29 pigs. Six received RAP, six received RAP with dronedarone (RAP/D), seven received dronedarone alone, four received RAP with amiodarone (RAP/A), and six received neither (sham). In ventricular tissue, oxidative stress/ischaemia-related gene and protein expression was evaluated by RT-PCR and Western blotting; Isoprostanes were measured by GC-MS procedures.

KEY RESULTS

CFR was decreased in the RAP group, compared with other groups. FFR was not different between groups. Effective refractory period was reduced in RAP compared with RAP/D. RAP-activated PKC phosphorylation tended to be decreased by dronedarone (P= 0.055) RAP induced NOX-1 and NOX-2 protein and the mRNA for hypoxia-inducible factor-1α (HIF-1α). Dronedarone reduced the pacing-dependent increase in the expression of NOX-2 protein and of HIF-1α mRNA. The oxidative stress marker, F(2)-isoprostane, was increased by RAP and this increase was attenuated by dronedarone. Other oxidative stress/ischaemia-related genes were induced by RAP compared with sham and were decreased by dronedarone treatment. In HL1 cells, dronedarone significantly inhibited the increased phosphorylation of PKCα after oxidative stress, with an almost significant effect (P= 0.059) on that after RAP.

CONCLUSIONS AND IMPLICATIONS

Dronedarone abolished RAP-induced ventricular microcirculatory abnormalities by decreasing oxidative stress/ischaemia-related gene and protein expression in the ventricle.

Pig bone marrow-derived macrophages resemble human macrophages in their response to bacterial lipopolysaccharide

Mouse bone marrow-derived macrophages (BMDM) grown in M-CSF (CSF-1) have been used widely in studies of macrophage biology and the response to TLR agonists. We investigated whether similar cells could be derived from the domestic pig using human rCSF-1 and whether porcine macrophages might represent a better model of human macrophage biology. Cultivation of pig bone marrow cells for 5-7 d in presence of human rCSF-1 generated a pure population of BMDM that expressed the usual macrophage markers (CD14, CD16, and CD172a), were potent phagocytic cells, and produced TNF in response to LPS. Pig BMDM could be generated from bone marrow cells that had been stored frozen and thawed so that multiple experiments can be performed on samples from a single animal. Gene expression in pig BMDM from outbred animals responding to LPS was profiled using Affymetrix microarrays. The temporal cascade of inducible and repressible genes more closely resembled the known responses of human than mouse macrophages, sharing with humans the regulation of genes involved in tryptophan metabolism (IDO, KYN), lymphoattractant chemokines (CCL20, CXCL9, CXCL11, CXCL13), and the vitamin D3-converting enzyme, Cyp27B1. Conversely, in common with published studies of human macrophages, pig BMDM did not strongly induce genes involved in arginine metabolism, nor did they produce NO. These results establish pig BMDM as an alternative tractable model for the study of macrophage transcriptional control.

Distinct peripheral blood RNA responses to Salmonella in pigs differing in Salmonella shedding levels: intersection of IFNG, TLR and miRNA pathways

Transcriptomic analysis of the response to bacterial pathogens has been reported for several species, yet few studies have investigated the transcriptional differences in whole blood in subjects that differ in their disease response phenotypes. Salmonella species infect many vertebrate species, and pigs colonized with Salmonella enterica serovar Typhimurium (ST) are usually asymptomatic, making detection of these Salmonella-carrier pigs difficult. The variable fecal shedding of Salmonella is an important cause of foodborne illness and zoonotic disease. To investigate gene pathways and biomarkers associated with the variance in Salmonella shedding following experimental inoculation, we initiated the first analysis of the whole blood transcriptional response induced by Salmonella. A population of pigs (n = 40) was inoculated with ST and peripheral blood and fecal Salmonella counts were collected between 2 and 20 days post-inoculation (dpi). Two groups of pigs with either low shedding (LS) or persistent shedding (PS) phenotypes were identified. Global transcriptional changes in response to ST inoculation were identified by Affymetrix Genechip® analysis of peripheral blood RNA at day 0 and 2 dpi. ST inoculation triggered substantial gene expression changes in the pigs and there was differential expression of many genes between LS and PS pigs. Analysis of the differential profiles of gene expression within and between PS and LS phenotypic classes identified distinct regulatory pathways mediated by IFN-γ, TNF, NF-κB, or one of several miRNAs. We confirmed the activation of two regulatory factors, SPI1 and CEBPB, and demonstrated that expression of miR-155 was decreased specifically in the PS animals. These data provide insight into specific pathways associated with extremes in Salmonella fecal shedding that can be targeted for further exploration on why some animals develop a carrier state. This knowledge can also be used to develop rational manipulations of genetics, pharmaceuticals, nutrition or husbandry methods to decrease Salmonella colonization, shedding and spread.

Left ventricular remodeling in swine after myocardial infarction: a transcriptional genomics approach

Despite the apparent appropriateness of left ventricular (LV) remodeling following myocardial infarction (MI), it poses an independent risk factor for development of heart failure. There is a paucity of studies into the molecular mechanisms of LV remodeling in large animal species. We took an unbiased molecular approach to identify candidate transcription factors (TFs) mediating the genetic reprogramming involved in post-MI LV remodeling in swine. Left ventricular tissue was collected from remote, non-infarcted myocardium, 3 weeks after MI-induction or sham-surgery. Microarray analysis identified 285 upregulated and 278 downregulated genes (FDR < 0.05). Of these differentially expressed genes, the promoter regions of the human homologs were searched for common TF binding sites (TFBS). Eighteen TFBS were overrepresented >two-fold (p < 0.01) in upregulated and 13 in downregulated genes. Left ventricular nuclear protein extracts were assayed for DNA-binding activity by protein/DNA array. Out of 345 DNA probes, 30 showed signal intensity changes >two-fold. Five TFs were identified in both TFBS and protein/DNA array analyses, which showed matching changes for COUP-TFII and glucocorticoid receptor (GR) only. Treatment of swine with the GR antagonist mifepristone after MI reduced the post-MI increase in LV mass, but LV dilation remained unaffected. Thus, using an unbiased approach to study post-MI LV remodeling in a physiologically relevant large animal model, we identified COUP-TFII and GR as potential key mediators of post-MI remodeling.

Impact of exercise training on endothelial transcriptional profiles in healthy swine: a genome-wide microarray analysis

While the salutary effects of exercise training on conduit artery endothelial cells have been reported in animals and humans with cardiovascular risk factors or disease, whether a healthy endothelium is alterable with exercise training is less certain. The purpose of this study was to evaluate the impact of exercise training on transcriptional profiles in normal endothelial cells using a genome-wide microarray analysis. Brachial and internal mammary endothelial gene expression was compared between a group of healthy pigs that exercise trained for 16-20 wk (n = 8) and a group that remained sedentary (n = 8). We found that a total of 130 genes were upregulated and 84 genes downregulated in brachial artery endothelial cells with exercise training (>1.5-fold and false discovery rate <15%). In contrast, a total of 113 genes were upregulated and 31 genes downregulated in internal mammary artery endothelial cells using the same criteria. Although there was an overlap of 66 genes (59 upregulated and 7 downregulated with exercise training) between the brachial and internal mammary arteries, the identified endothelial gene networks and biological processes influenced by exercise training were distinctly different between the brachial and internal mammary arteries. These data indicate that a healthy endothelium is indeed responsive to exercise training and support the concept that the influence of physical activity on endothelial gene expression is not homogenously distributed throughout the vasculature.

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 differences in healthy brachial and femoral arteries of Rapacz familial hypercholesterolemic swine

The mechanisms underlying the unequal distribution of atherosclerotic disease in the peripheral arteries are currently unclear. Gene expression differences in healthy arteries may influence the heterogeneous distribution of atherosclerosis. Therefore, this investigation compares gene expression in healthy atheroprotected brachial and atherosusceptible femoral arteries of young and disease free Rapacz familial hypercholesterolemic (FHC) swine. We hypothesized that transcripts related to atherosusceptibility would be differentially expressed between these arteries prior to the onset of disease. Femoral and brachial arteries were harvested from four 13-day-old Rapacz FHC swine. No atherosclerotic disease was detected using Sudan IV, Verhoeff-van Gieson, and hematoxylin-eosin staining. Gene expression was quantified using Affymetrix GeneChip Porcine Genome Arrays. An average of 15,552 probe sets had detectable transcripts, while 430 probe sets showed a significant differential expression between arteries (false discovery rate < 0.05). The human orthologs of 63 probe sets with differential expression and a 1.5-fold or greater transcript abundance between arteries are associated with Wnt/β-catenin, lysophospholipid, and Ca-signaling, as well as apoptosis. This is the first investigation reporting that differences in relative abundance of gene expression exist between brachial and femoral arteries in young Rapacz FHC swine prior to the development of atherosclerotic lesions.

Porcine S100A8 and S100A9: molecular characterizations and crucial functions in response to Haemophilus parasuis infection

S100 calcium-binding protein A8 (S100A8) and S100 calcium-binding protein A9 (S100A9) are pivotal mediators of inflammatory and protective anti-infection responses for the mammalian host. In this study, we present the molecular cloning of porcine S100A8 (pS100A8) and porcine S100A9 (pS100A9). Both genes comprise 3 exons and 2 introns and are located on pig chromosome 4q21-q23 (closely linked to SW512). Homology comparison to other mammalian species affirmed that critical functional amino acids for post-transcriptional modification, inflammatory regulation, and formation of heterodimeric complexes exist in pS100A8 and pS100A9. Under normal conditions, both genes are preferentially expressed in porcine immune or immune-related organs, e.g., bone marrow, spleen, lymph nodes, and lung. Upon stimulation in porcine whole blood cultures with LPS or Poly(I:C), they are dramatically induced. Interestingly, the maximum increase of mRNA levels in blood cultures of Meishan pigs is significantly greater than that in Duroc pigs. We previously showed that pS100A8 and pS100A9 mRNA were up-regulated following Haemophilus parasuis (HPS) infection. We herein further confirm their up-regulation at the protein level in multiple HPS infected tissues (spleen, lung and liver). Functional cluster and network analysis based on our previous microarray data discovered that CEBPB may be one of the key transcription factors. A pS100A8/pS100A9-CASP3-SLC1A2 pathway regulating lipid metabolism was found. Both of their pro- and anti-inflammatory functions in response to HPS infection are highlighted.

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.