sigReannot: an oligo-set re-annotation pipeline based on similarities with the Ensembl transcripts and Unigene clusters

Identification of the Efficacy of Ex Situ Conservation of Ammopiptanthus nanus Based on Its ETS-SSR Markers

Ammopiptanthus possesses ancestral traits and, as a tertiary relict, is one of the surviving remnants of the ancient Mediterranean retreat and climate drought. It is also the only genus of super xerophytic, evergreen, broad-leaved shrubs. Ammopiptanthus nanus, one of the two species in this genus, is predominantly found in extremely arid and frigid environments, and is increasingly threatened with extinction. Study of the species' genetic diversity is thus beneficial for its survival and the efficacy of ex situ conservation efforts. Based on transcriptome data, 15 pairs of effective EST-SSR were screened to evaluate A. nanus genetic diversity. In all, 87 samples from three populations were evaluated, the results of which show that ex situ conservation in the Wuqia region needs to be supplemented. Conservation and breeding of individual A. nanus offspring should be strengthened in the future to ensure their progeny continue to exhibit high genetic diversity and variation.

Molecular signatures of muscle growth and composition deciphered by the meta-analysis of age-related public transcriptomics data

The lean-to-fat ratio is a major issue in the beef meat industry from both carcass and meat production perspectives. This industrial perspective has motivated meat physiologists to use transcriptomics technologies to decipher mechanisms behind fat deposition within muscle during the time course of muscle growth. However, synthetic biological information from this volume of data remains to be produced to identify mechanisms found in various breeds and rearing practices. We conducted a meta-analysis on 10 transcriptomic data sets stored in public databases, from the longissimus thoracis of five different bovine breeds divergent by age. We updated gene identifiers on the last version of the bovine genome (UCD1.2), and the 715 genes common to the 10 studies were subjected to the meta-analysis. Of the 238 genes differentially expressed (DEG), we identified a transcriptional signature of the dynamic regulation of glycolytic and oxidative metabolisms that agrees with a known shift between those two pathways from the animal puberty. We proposed some master genes of the myogenesis, namely MYOG and MAPK14, as probable regulators of the glycolytic and oxidative metabolisms. We also identified overexpressed genes related to lipid metabolism (APOE, LDLR, MXRA8, and HSP90AA1) that may contribute to the expected enhanced marbling as age increases. Lastly, we proposed a transcriptional signature related to the induction (YBX1) or repression (MAPK14, YWAH, ERBB2) of the commitment of myogenic progenitors into the adipogenic lineage. The relationships between the abundance of the identified mRNA and marbling values remain to be analyzed in a marbling biomarkers discovery perspectives.

Transcriptome analyses of liver in newly-hatched chicks during the metabolic perturbation of fasting and re-feeding reveals THRSPA as the key lipogenic transcription factor

Background

The fasting-refeeding perturbation has been used extensively to reveal specific genes and metabolic pathways that control energy metabolism in the chicken. Most global transcriptional scans of the fasting-refeeding response in liver have focused on juvenile chickens that were 1, 2 or 4 weeks old. The present study was aimed at the immediate post-hatch period, in which newly-hatched chicks were subjected to fasting for 4, 24 or 48 h, then refed for 4, 24 or 48 h, and compared with a fully-fed control group at each age (D1-D4).

Results

Visual analysis of hepatic gene expression profiles using hierarchical and K-means clustering showed two distinct patterns, genes with higher expression during fasting and depressed expression upon refeeding and those with an opposing pattern of expression, which exhibit very low expression during fasting and more abundant expression with refeeding. Differentially-expressed genes (DEGs), identified from five prominent pair-wise contrasts of fed, fasted and refed conditions, were subjected to Ingenuity Pathway Analysis. This enabled mapping of analysis-ready (AR)-DEGs to canonical and metabolic pathways controlled by distinct gene interaction networks. The largest number of hepatic DEGs was identified by two contrasts: D2FED48h/D2FAST48h (968 genes) and D2FAST48h/D3REFED24h (1198 genes). The major genes acutely depressed by fasting and elevated upon refeeding included ANGTPL, ATPCL, DIO2, FASN, ME1, SCD, PPARG, SREBP2 and THRSPA—a primary lipogenic transcription factor. In contrast, major lipolytic genes were up-regulated by fasting or down-regulated after refeeding, including ALDOB, IL-15, LDHB, LPIN2, NFE2L2, NR3C1, NR0B1, PANK1, PPARA, SERTAD2 and UPP2.

Conclusions

Transcriptional profiling of liver during fasting/re-feeding of newly-hatched chicks revealed several highly-expressed upstream regulators, which enable the metabolic switch from fasted (lipolytic/gluconeogenic) to fed or refed (lipogenic/thermogenic) states. This rapid homeorhetic shift of whole-body metabolism from a catabolic-fasting state to an anabolic-fed state appears precisely orchestrated by a small number of ligand-activated transcription factors that provide either a fasting-lipolytic state (PPARA, NR3C1, NFE2L2, SERTAD2, FOX01, NR0B1, RXR) or a fully-fed and refed lipogenic/thermogenic state (THRSPA, SREBF2, PPARG, PPARD, JUN, ATF3, CTNNB1). THRSPA has emerged as the key transcriptional regulator that drives lipogenesis and thermogenesis in hatchling chicks, as shown here in fed and re-fed states.

Transcriptional profiling of liver in riboflavin-deficient chicken embryos explains impaired lipid utilization, energy depletion, massive hemorrhaging, and delayed feathering

BACKGROUND

A strain of Leghorn chickens (rd/rd), unable to produce a functional riboflavin-binding protein, lays riboflavin-deficient eggs, in which all embryos suddenly die at mid-incubation (days 13-15). This malady, caused by riboflavin deficiency, leads to excessive lipid accumulation in liver, impaired β-oxidation of lipid, and severe hypoglycemia prior to death. We have used high-density chicken microarrays for time-course transcriptional scans of liver in chicken embryos between days 9-15 during this riboflavin-deficiency-induced metabolic catastrophe. For comparison, half of rd/rd embryos (n = 16) were rescued from this calamity by injection of riboflavin just prior to incubation of fertile eggs from rd/rd hens.

RESULTS

No significant differences were found between hepatic transcriptomes of riboflavin-deficient and riboflavin-rescued embryos at the first two ages (days 9 and 11). Overall, we found a 3.2-fold increase in the number of differentially expressed hepatic genes between day 13 (231 genes) and day 15 (734 genes). Higher expression of genes encoding the chicken flavoproteome was more evident in rescued- (15 genes) than in deficient-embryos (4 genes) at day 15. Diminished activity of flavin-dependent enzymes in riboflavin-deficient embryos blocks catabolism of yolk lipids, which normally serves as the predominant source of energy required for embryonic development.

CONCLUSIONS

Riboflavin deficiency in mid-stage embryos leads to reduced expression of numerous genes controlling critical functions, including β-oxidation of lipids, blood coagulation and feathering. Surprisingly, reduced expression of feather keratin 1 was found in liver of riboflavin-deficient embryos at e15, which could be related to their delayed feathering and sparse clubbed down. A large number of genes are expressed at higher levels in liver of riboflavin-deficient embryos; these up-regulated genes control lipid storage/transport, gluconeogenesis, ketogenesis, protein catabolism/ubiquitination and cell death.

Divergent selection for residual feed intake affects the transcriptomic and proteomic profiles of pig skeletal muscle

Improving feed efficiency is a relevant strategy to reduce feed cost and environmental waste in livestock production. Selection experiments on residual feed intake (RFI), a measure of feed efficiency, previously indicated that low RFI was associated with lower feed intake, similar growth rate, and greater lean meat content compared with high RFI. To gain insights into the molecular mechanisms underlying these differences, 24 Large White females from 2 lines divergently selected for RFI were examined. Pigs from a low-RFI ("efficient") and high-RFI ("inefficient") line were individually fed ad libitum from 67 d of age (27 kg BW) to slaughter at 115 kg BW (n = 8 per group). Additional pigs of the high-RFI line were feed restricted to the daily feed intake of the ad libitum low-RFI pigs (n = 8) to investigate the impact of selection independently of feed intake. Global gene and protein expression profiles were assessed in the LM collected at slaughter. The analyses involved a porcine commercial microarray and 2-dimensional gel electrophoresis. About 1,000 probes were differentially expressed (P < 0.01) between RFI lines. Only 10% of those probes were also affected by feed restriction. Gene functional classification indicated a greater expression of genes involved in protein synthesis and a lower expression of genes associated with mitochondrial energy metabolism in the low-RFI pigs compared with the high-RFI pigs. At the protein level, 11 unique identified proteins exhibited a differential abundance (P < 0.05) between RFI lines. Differentially expressed proteins were generally not significantly affected by feed restriction. Mitochondrial oxidative proteins such as aconitase hydratase, ATP synthase subunit α, and creatine kinase S-type had a lower abundance in the low-RFI pigs, whereas fructose-biphosphate aldolase A and glyceraldehyde-3-phosphate dehydrogenase, 2 proteins involved in glycolysis, had a greater abundance in those pigs compared with high-RFI pigs. Antioxidant proteins such as superoxide dismutase and glutathione peroxidase 3 at the mRNA level and peroxiredoxin-6 at the protein level were also less expressed in LM of the most efficient pigs, likely related to lower oxidative molecule production. Collectively, both the transcriptomic and proteomic approaches revealed a lower oxidative metabolism in muscle of the low-RFI pigs and all these modifications were largely independent of differences in feed intake.

The fetal ovary exhibits temporal sensitivity to a 'real-life' mixture of environmental chemicals

The development of fetal ovarian follicles is a critical determinant of adult female reproductive competence. Prolonged exposure to environmental chemicals (ECs) can perturb this process with detrimental consequences for offspring. Here we report on the exposure of pregnant ewes to an environmental mixture of ECs derived from pastures fertilized with sewage sludge (biosolids): a common global agricultural practice. Exposure of pregnant ewes to ECs over 80 day periods during early, mid or late gestation reduced the proportion of healthy early stage fetal follicles comprising the ovarian reserve. Mid and late gestation EC exposures had the most marked effects, disturbing maternal and fetal liver chemical profiles, masculinising fetal anogenital distance and greatly increasing the number of altered fetal ovarian genes and proteins. In conclusion, differential temporal sensitivity of the fetus and its ovaries to EC mixtures has implications for adult ovarian function following adverse exposures during pregnancy.

Defining Mononuclear Phagocyte Subset Homology Across Several Distant Warm-Blooded Vertebrates Through Comparative Transcriptomics

Mononuclear phagocytes are organized in a complex system of ontogenetically and functionally distinct subsets, that has been best described in mouse and to some extent in human. Identification of homologous mononuclear phagocyte subsets in other vertebrate species of biomedical, economic, and environmental interest is needed to improve our knowledge in physiologic and physio-pathologic processes, and to design intervention strategies against a variety of diseases, including zoonotic infections. We developed a streamlined approach combining refined cell sorting and integrated comparative transcriptomics analyses which revealed conservation of the mononuclear phagocyte organization across human, mouse, sheep, pigs and, in some respect, chicken. This strategy should help democratizing the use of omics analyses for the identification and study of cell types across tissues and species. Moreover, we identified conserved gene signatures that enable robust identification and universal definition of these cell types. We identified new evolutionarily conserved gene candidates and gene interaction networks for the molecular regulation of the development or functions of these cell types, as well as conserved surface candidates for refined subset phenotyping throughout species. A phylogenetic analysis revealed that orthologous genes of the conserved signatures exist in teleost fishes and apparently not in Lamprey.

The Longissimus and Semimembranosus muscles display marked differences in their gene expression profiles in pig

BACKGROUND

Meat quality depends on skeletal muscle structure and metabolic properties. While most studies carried on pigs focus on the Longissimus muscle (LM) for fresh meat consumption, Semimembranosus (SM) is also of interest because of its importance for cooked ham production. Even if both muscles are classified as glycolytic muscles, they exhibit dissimilar myofiber composition and metabolic characteristics. The comparison of LM and SM transcriptome profiles undertaken in this study may thus clarify the biological events underlying their phenotypic differences which might influence several meat quality traits.

METHODOLOGY/PRINCIPAL FINDINGS

Muscular transcriptome analyses were performed using a custom pig muscle microarray: the 15 K Genmascqchip. A total of 3823 genes were differentially expressed between the two muscles (Benjamini-Hochberg adjusted P value ≤0.05), out of which 1690 and 2133 were overrepresented in LM and SM respectively. The microarray data were validated using the expression level of seven differentially expressed genes quantified by real-time RT-PCR. A set of 1047 differentially expressed genes with a muscle fold change ratio above 1.5 was used for functional characterization. Functional annotation emphasized five main clusters associated to transcriptome muscle differences. These five clusters were related to energy metabolism, cell cycle, gene expression, anatomical structure development and signal transduction/immune response.

CONCLUSIONS/SIGNIFICANCE

This study revealed strong transcriptome differences between LM and SM. These results suggest that skeletal muscle discrepancies might arise essentially from different post-natal myogenic activities.

Dendritic cell subtypes from lymph nodes and blood show contrasted gene expression programs upon Bluetongue virus infection

Human and animal hemorrhagic viruses initially target dendritic cells (DCs). It has been proposed, but not documented, that both plasmacytoid DCs (pDCs) and conventional DCs (cDCs) may participate in the cytokine storm encountered in these infections. In order to evaluate the contribution of DCs in hemorrhagic virus pathogenesis, we performed a genome-wide expression analysis during infection by Bluetongue virus (BTV), a double-stranded RNA virus that induces hemorrhagic fever in sheep and initially infects cDCs. Both pDCs and cDCs accumulated in regional lymph nodes and spleen during BTV infection. The gene response profiles were performed at the onset of the disease and markedly differed with the DC subtypes and their lymphoid organ location. An integrative knowledge-based analysis revealed that blood pDCs displayed a gene signature related to activation of systemic inflammation and permeability of vasculature. In contrast, the gene profile of pDCs and cDCs in lymph nodes was oriented to inhibition of inflammation, whereas spleen cDCs did not show a clear functional orientation. These analyses indicate that tissue location and DC subtype affect the functional gene expression program induced by BTV and suggest the involvement of blood pDCs in the inflammation and plasma leakage/hemorrhage during BTV infection in the real natural host of the virus. These findings open the avenue to target DCs for therapeutic interventions in viral hemorrhagic diseases.

Gene coexpression network analysis identifies genes and biological processes shared among anterior pituitary and brain areas that affect estrous behavior in dairy cows

The expression of estrous (sexually receptive) behavior (EB), a key fertility trait in dairy cows, has been declining over the past few decades both in intensity and duration. Improved knowledge of the genomic factors underlying EB, which is currently lacking, may lead to novel applications to enhance fertility. Our objective was to identify genes and biological processes shared among the bovine anterior pituitary (AP) and four brain areas that act together to regulate EB by investigating networks of coexpressed genes between these tissues. We used a systems biology approach called weighted gene coexpression network analysis for defining gene coexpression networks using gene expression data from the following tissues collected from 14 cows at estrus: AP, dorsal hypothalamus (DH), ventral hypothalamus (VH), amygdala (AM), and hippocampus (HC). Consensus modules of coexpressed genes were identified between the networks for the AM-DH, HC-DH, VH-DH, AP-DH, and AM-HC tissue pairs. The correlation between the module's eigengene (weighted average gene expression profile) and levels of EB exhibited by the experimental cows were tested. Estrous behavior-correlated modules were found enriched for gene ontology terms like glial cell development and regulation of neural projection development as well as for Kyoto Encyclopedia of Genes and Genomes pathway terms related to brain degenerative diseases. General cellular processes like oxidative phosphorylation and ribosome and biosynthetic processes were found enriched in several correlated modules, indicating increased transcription and protein synthesis. Stimulation of ribosomal RNA synthesis is known from rodent studies to be a primary event in the activation of neuronal cells and pathways involved in female reproductive behavior and this precedes the estrogen-driven expansion of dendrites and synapses. Similar processes also operate in cows to affect EB. Hub genes within EB-correlated modules (e.g. NEFL, NDRG2, GAP43, THY1, and TCF7L2, among others) are strong candidates among genes regulating EB expression. The study improved our understanding of the genomic regulation of EB in dairy cows by providing new insights into genes and biological processes shared among the bovine AP and brain areas acting together to regulate EB. The new knowledge could lead to the development of novel management strategies to monitor and improve reproductive performance in dairy cows (for example, biomarkers for estrus detection).

Canine recombinant adenovirus vector induces an immunogenicity-related gene expression profile in skin-migrated CD11b⁺ -type DCs

Gene expression profiling of the blood cell response induced early after vaccination has previously been demonstrated to predict the immunogenicity of vaccines. In this study, we evaluated whether the analysis of the gene expression profile of skin-migrated dendritic cells (DCs) could be informative for the in vitro prediction of immunogenicity of vaccine, using canine adenovirus serotype 2 (CAV2) as vaccine vector. CAV2 has been shown to induce immunity to transgenes in several species including sheep and is an interesting alternative to human adenovirus-based vectors, based on the safety records of the parental strain in dogs and the lack of pre-existing immunity in non-host species. Skin-migrated DCs were collected from pseudo-afferent lymph in sheep. Both the CD11b(+) -type and CD103(+) -type skin-migrated DCs were transduced by CAV2. An analysis of the global gene response to CAV2 in the two skin DC subsets showed that the gene response in CD11b(+) -type DCs was far higher and broader than in the CD103(+) -type DCs. A newly released integrative analytic tool from Ingenuity systems revealed that the CAV2-modulated genes in the CD11b(+) -type DCs clustered in several activated immunogenicity-related functions, such as immune response, immune cell trafficking and inflammation. Thus gene profiling in skin-migrated DC in vitro indicates that the CD11b(+) DC type is more responsive to CAV2 than the CD103(+) DC type, and provides valuable information to help in evaluating and possibly improving viral vector vaccine effectiveness.

Mitochondrial function, fatty acid metabolism, and immune system are relevant features of pig adipose tissue development

The molecular mechanisms underlying the genetic control of fat development in humans and livestock species still require characterization. To gain insights on gene expression patterns associated with genetic propensity for adiposity, we compared subcutaneous adipose tissue (SCAT) transcriptomics profiles from two contrasted pig breeds for body fatness. Samples were obtained from Large White (LW; lean phenotype) and Basque pigs (B; low growth and high fat content) at 35 kg (n = 5 per breed) or 145 kg body weight (n = 10 per breed). Using a custom adipose tissue microarray, we found 271 genes to be differentially expressed between the two breeds at both stages, out of which 123 were highly expressed in LW pigs and 148 genes were highly expressed in B pigs. Functional enrichment analysis based on gene ontology (GO) terms highlighted gene groups corresponding to the mitochondrial energy metabolism in LW pigs, whereas immune response was found significantly enriched in B pigs. Genes associated with lipid metabolism, such as ELOVL6, a gene involved in fatty acid elongation, had a lower expression in B compared with LW pigs. Furthermore, despite enlarged adipocyte diameters and higher plasma leptin concentration, B pigs displayed reduced lipogenic enzyme activities compared with LW pigs at 145 kg. Altogether, our results suggest that the development of adiposity was associated with a progressive worsening of the metabolic status, leading to a low-grade inflammatory state, and may thus be of significant interest for both livestock production and human health.

Transcriptomic profiling of proteases and antiproteases in the liver of sexually mature hens in relation to vitellogenesis

Background

Most egg yolk precursors are synthesized by the liver, secreted into the blood and transferred into oocytes, to provide nutrients and bioactive molecules for the avian embryo. Three hundred and sixteen distinct proteins have been identified in egg yolk. These include 37 proteases and antiproteases, which are likely to play a role in the formation of the yolk (vitellogenesis), as regulators of protein metabolism. We used a transcriptomic approach to define the protease and antiprotease genes specifically expressed in the hen liver in relation to vitellogenesis by comparing sexually mature and pre-laying chickens showing different steroid milieu.

Results

Using a 20 K chicken oligoarray, a total of 582 genes were shown to be over-expressed in the liver of sexually mature hens (1.2 to 67 fold-differences). Eight of the top ten over-expressed genes are known components of the egg yolk or perivitelline membrane. This list of 582 genes contains 12 proteases and 3 antiproteases. We found that “uncharacterized protein LOC419301/similar to porin” (GeneID:419301), an antiprotease and “cathepsin E-A-like/similar to nothepsin” (GeneID:417848), a protease, were the only over-expressed candidates (21-fold and 35-fold difference, respectively) that are present in the egg yolk. Additionally, we showed the 4-fold over-expression of “ovochymase-2/similar to oviductin” (GeneID:769290), a vitelline membrane-specific protease.

Conclusions

Our approach revealed that three proteases and antiproteases are likely to participate in the formation of the yolk. The role of the other 12 proteases and antiproteases which are over-expressed in our model remains unclear. At least 1/3 of proteases and antiproteases identified in egg yolk and vitelline membrane proteomes are expressed similarly in the liver regardless of the maturity of hens, and have been initially identified as regulators of haemostasis and inflammatory events. The lack of effect of sex steroids on these genes expressed in the liver but the products of which are found in the yolk suggests that these may be passively incorporated into the yolk rather than actively produced for that purpose. These results raise the question of the biological significance of egg yolk proteases and antiproteases, and more generally of all minor proteins that have been identified in egg yolk.

Transcriptional responses of resistant and susceptible fish clones to the bacterial pathogen Flavobacterium psychrophilum

Flavobacterium psychrophilum is a bacterial species that represents one of the most important pathogens for aquaculture worldwide, especially for salmonids. To gain insights into the genetic basis of the natural resistance to F. psychrophilum, we selected homozygous clones of rainbow trout with contrasted susceptibility to the infection. We compared the transcriptional response to the bacteria in the pronephros of a susceptible and a resistant line by micro-array analysis five days after infection. While the basal transcriptome of healthy fish was significantly different in the resistant and susceptible lines, the transcriptome modifications induced by the bacteria involved essentially the same genes and pathways. The response to F. psychrophilum involved antimicrobial peptides, complement, and a number of enzymes and chemokines. The matrix metalloproteases mmp9 and mmp13 were among the most highly induced genes in both genetic backgrounds. Key genes of both pro- and anti-inflammatory response such as IL1 and IL10, were up-regulated with a greater magnitude in susceptible animals where the bacterial load was also much higher. While higher resistance to F. psychrophilum does not seem to be based on extensive differences in the orientation of the immune response, several genes including complement C3 showed stronger induction in the resistant fish. They may be important for the variation of susceptibility to the infection.

Genetic susceptibility to S. aureus mastitis in sheep: differential expression of mammary epithelial cells in response to live bacteria or supernatant

Staphylococcus aureus is a prevalent pathogen for mastitis in dairy ruminants and is responsible for both clinical and subclinical mastitis. Mammary epithelial cells (MEC) represent not only a physical barrier against bacterial invasion but are also active players of the innate immune response permitting infection clearance. To decipher their functions in general and in animals showing different levels of genetic predisposition to Staphylococcus in particular, MEC from ewes undergoing a divergent selection on milk somatic cell count were stimulated by S. aureus. MEC response was also studied according to the stimulation condition with live bacteria or culture supernatant. The early MEC response was studied during a 5 h time course by microarray to identify differentially expressed genes with regard to the host genetic background and as a function of the conditions of stimulation. In both conditions of stimulation, metabolic processes were altered, the apoptosis-associated pathways were considerably modified, and inflammatory and immune responses were enhanced with the upregulation of il1a, il1b, and tnfa and several chemokines known to enhance neutrophil (cxcl8) or mononuclear leukocyte (ccl20) recruitment. Genes associated with oxidative stress were increased after live bacteria stimulation, whereas immune response-related genes were higher after supernatant stimulation in the early phase. Only 20 genes were differentially expressed between Staphylococcus spp-mastitis resistant and susceptible animals without any clearly defined role on the control of infection. To conclude, this suggests that MEC may not represent the cell type at the origin of the difference of mastitis susceptibility, at least as demonstrated in our genetic model. Supernatant or heat-killed S. aureus produce biological effects that are essentially different from those induced by live bacteria.

SigReannot-mart: a query environment for expression microarray probe re-annotations

Expression microarrays are commonly used to study transcriptomes. Most of the arrays are now based on oligo-nucleotide probes. Probe design being a tedious task, it often takes place once at the beginning of the project. The oligo set is then used for several years. During this time period, the knowledge gathered by the community on the genome and the transcriptome increases and gets more precise. Therefore re-annotating the set is essential to supply the biologists with up-to-date annotations. SigReannot-mart is a query environment populated with regularly updated annotations for different oligo sets. It stores the results of the SigReannot pipeline that has mainly been used on farm and aquaculture species. It permits easy extraction in different formats using filters. It is used to compare probe sets on different criteria, to choose the set for a given experiment to mix probe sets in order to create a new one.

Database URL:http://sigreannot-mart.toulouse.inra.fr/

Infection of nonhost species dendritic cells in vitro with an attenuated myxoma virus induces gene expression that predicts its efficacy as a vaccine vector

Recombinant myxoma virus (MYXV) can be produced without a loss of infectivity, and its highly specific host range makes it an ideal vaccine vector candidate, although careful examination of its interaction with the immune system is necessary. Similar to rabbit bone marrow-derived dendritic cells (BM-DCs), ovine dendritic cells can be infected by SG33, a MYXV vaccine strain, and support recombinant antigen expression. The frequency of infected cells in the nonhost was lower and the virus cycle was abortive in these cell types. Among BM-DC subpopulations, Langerhans cell-like DCs were preferentially infected at low multiplicities of infection. Interestingly, ovine BM-DCs remained susceptible to MYXV after maturation, although apoptosis occurred shortly after infection as a function of the virus titer. When gene expression was assessed in infected BM-DC cultures, type I interferon (IFN)-related and inflammatory genes were strongly upregulated. DC gene expression profiles were compared with the profiles produced by other poxviruses in interaction with DCs, but very few commonalities were found, although genes that were previously shown to predict vaccine efficacy were present. Collectively, these data support the idea that MYXV permits efficient priming of adaptive immune responses and should be considered a promising vaccine vector along with other poxviruses.

Regional regulation of transcription in the bovine genome

Eukaryotic genes are distributed along chromosomes as clusters of highly expressed genes termed RIDGEs (Regions of IncreaseD Gene Expression) and lowly expressed genes termed anti-RIDGEs, interspersed among genes expressed at intermediate levels or not expressed. Previous studies based on this observation suggested a dual mechanism of gene regulation, where, in addition to transcription factors, the chromosomal domain influences the expression level of their embedded genes. The objectives here were to provide evidence for the existence of chromosomal regional regulation of transcription in the bovine genome, to analyse the genomic features of genes located within RIDGEs versus anti-RIDGEs and tissue-specific genes versus housekeeping and to examine the genomic distribution of genes subject to positive selection in bovines. Gene expression analysis of four brain tissues and the anterior pituitary of 28 cows identified 70 RIDGEs and 41 anti-RIDGEs (harbouring 3735 and 1793 bovine genes respectively) across the bovine genome which are significantly higher than expected by chance. Housekeeping genes (defined here as genes expressed in all five tissues) were over-represented within RIDGEs but tissue-specific genes (genes expressed in only one of the five tissues) were not. Housekeeping genes and genes within RIDGEs had, in general, higher expression levels and GC content but shorter gene lengths and intron lengths than tissue-specific genes and genes within anti-RIDGES. Our findings suggest the existence of chromosomal regional regulation of transcription in the bovine genome. The genomic features observed for genes within RIDGEs and housekeeping genes in bovines agree with previous studies in several other species further strengthening the hypothesis of selective pressure to keep the highly and widely expressed genes short and compact for transcriptional efficiency. Further, positively selected genes were found non-randomly distributed on the genome with a preference for RIDGEs and regions of intermediate gene expression compared to anti-RIDGEs.

Strengthening insights into host responses to mastitis infection in ruminants by combining heterogeneous microarray data sources

Background

Gene expression profiling studies of mastitis in ruminants have provided key but fragmented knowledge for the understanding of the disease. A systematic combination of different expression profiling studies via meta-analysis techniques has the potential to test the extensibility of conclusions based on single studies. Using the program Pointillist, we performed meta-analysis of transcription-profiling data from six independent studies of infections with mammary gland pathogens, including samples from cattle challenged in vivo with S. aureus, E. coli, and S. uberis, samples from goats challenged in vivo with S. aureus, as well as cattle macrophages and ovine dendritic cells infected in vitro with S. aureus. We combined different time points from those studies, testing different responses to mastitis infection: overall (common signature), early stage, late stage, and cattle-specific.

Results

Ingenuity Pathway Analysis of affected genes showed that the four meta-analysis combinations share biological functions and pathways (e.g. protein ubiquitination and polyamine regulation) which are intrinsic to the general disease response. In the overall response, pathways related to immune response and inflammation, as well as biological functions related to lipid metabolism were altered. This latter observation is consistent with the milk fat content depression commonly observed during mastitis infection. Complementarities between early and late stage responses were found, with a prominence of metabolic and stress signals in the early stage and of the immune response related to the lipid metabolism in the late stage; both mechanisms apparently modulated by few genes, including XBP1 and SREBF1.

The cattle-specific response was characterized by alteration of the immune response and by modification of lipid metabolism. Comparison of E. coli and S. aureus infections in cattle in vivo revealed that affected genes showing opposite regulation had the same altered biological functions and provided evidence that E. coli caused a stronger host response.

Conclusions

This meta-analysis approach reinforces previous findings but also reveals several novel themes, including the involvement of genes, biological functions, and pathways that were not identified in individual studies. As such, it provides an interesting proof of principle for future studies combining information from diverse heterogeneous sources.

Transcriptomic analysis of milk somatic cells in mastitis resistant and susceptible sheep upon challenge with Staphylococcus epidermidis and Staphylococcus aureus

Background

The existence of a genetic basis for host responses to bacterial intramammary infections has been widely documented, but the underlying mechanisms and the genes are still largely unknown. Previously, two divergent lines of sheep selected for high/low milk somatic cell scores have been shown to be respectively susceptible and resistant to intramammary infections by Staphylococcus spp. Transcriptional profiling with an 15K ovine-specific microarray of the milk somatic cells of susceptible and resistant sheep infected successively by S. epidermidis and S. aureus was performed in order to enhance our understanding of the molecular and cellular events associated with mastitis resistance.

Results

The bacteriological titre was lower in the resistant than in the susceptible animals in the 48 hours following inoculation, although milk somatic cell concentration was similar. Gene expression was analysed in milk somatic cells, mainly represented by neutrophils, collected 12 hours post-challenge. A high number of differentially expressed genes between the two challenges indicated that more T cells are recruited upon inoculation by S. aureus than S. epidermidis. A total of 52 genes were significantly differentially expressed between the resistant and susceptible animals. Further Gene Ontology analysis indicated that differentially expressed genes were associated with immune and inflammatory responses, leukocyte adhesion, cell migration, and signal transduction. Close biological relationships could be established between most genes using gene network analysis. Furthermore, gene expression suggests that the cell turn-over, as a consequence of apoptosis/granulopoiesis, may be enhanced in the resistant line when compared to the susceptible line.

Conclusions

Gene profiling in resistant and susceptible lines has provided good candidates for mapping the biological pathways and genes underlying genetically determined resistance and susceptibility towards Staphylococcus infections, and opens new fields for further investigation.

Gene expression patterns in four brain areas associate with quantitative measure of estrous behavior in dairy cows

BACKGROUND

The decline noticed in several fertility traits of dairy cattle over the past few decades is of major concern. Understanding of the genomic factors underlying fertility, which could have potential applications to improve fertility, is very limited. Here, we aimed to identify and study those genes that associated with a key fertility trait namely estrous behavior, among genes expressed in four bovine brain areas (hippocampus, amygdala, dorsal hypothalamus and ventral hypothalamus), either at the start of estrous cycle, or at mid cycle, or regardless of the phase of cycle.

RESULTS

An average heat score was calculated for each of 28 primiparous cows in which estrous behavior was recorded for at least two consecutive estrous cycles starting from 30 days post-partum. Gene expression was then measured in brain tissue samples collected from these cows, 14 of which were sacrificed at the start of estrus and 14 around mid cycle. For each brain area, gene expression was modeled as a function of the orthogonally transformed average heat score values using a Bayesian hierarchical mixed model. Genes whose expression patterns showed significant linear or quadratic relationships with heat scores were identified. These included genes expected to be related to estrous behavior as they influence states like socio-sexual behavior, anxiety, stress and feeding motivation (OXT, AVP, POMC, MCHR1), but also genes whose association with estrous behavior is novel and warrants further investigation.

CONCLUSIONS

Several genes were identified whose expression levels in the bovine brain associated with the level of expression of estrous behavior. The genes OXT and AVP play major roles in regulating estrous behavior in dairy cows. Genes related to neurotransmission and neuronal plasticity are also involved in estrous regulation, with several genes and processes expressed in mid-cycle probably contributing to proper expression of estrous behavior in the next estrus. Studying these genes and the processes they control improves our understanding of the genomic regulation of estrous behavior expression.

A factor model to analyze heterogeneity in gene expression

BACKGROUND

Microarray technology allows the simultaneous analysis of thousands of genes within a single experiment. Significance analyses of transcriptomic data ignore the gene dependence structure. This leads to correlation among test statistics which affects a strong control of the false discovery proportion. A recent method called FAMT allows capturing the gene dependence into factors in order to improve high-dimensional multiple testing procedures. In the subsequent analyses aiming at a functional characterization of the differentially expressed genes, our study shows how these factors can be used both to identify the components of expression heterogeneity and to give more insight into the underlying biological processes.

RESULTS

The use of factors to characterize simple patterns of heterogeneity is first demonstrated on illustrative gene expression data sets. An expression data set primarily generated to map QTL for fatness in chickens is then analyzed. Contrarily to the analysis based on the raw data, a relevant functional information about a QTL region is revealed by factor-adjustment of the gene expressions. Additionally, the interpretation of the independent factors regarding known information about both experimental design and genes shows that some factors may have different and complex origins.

CONCLUSIONS

As biological information and technological biases are identified in what was before simply considered as statistical noise, analyzing heterogeneity in gene expression yields a new point of view on transcriptomic data.

Using transcriptome profiling to characterize QTL regions on chicken chromosome 5

Background

Although many QTL for various traits have been mapped in livestock, location confidence intervals remain wide that makes difficult the identification of causative mutations. The aim of this study was to test the contribution of microarray data to QTL detection in livestock species. Three different but complementary approaches are proposed to improve characterization of a chicken QTL region for abdominal fatness (AF) previously detected on chromosome 5 (GGA5).

Results

Hepatic transcriptome profiles for 45 offspring of a sire known to be heterozygous for the distal GGA5 AF QTL were obtained using a 20 K chicken oligochip. mRNA levels of 660 genes were correlated with the AF trait. The first approach was to dissect the AF phenotype by identifying animal subgroups according to their 660 transcript profiles. Linkage analysis using some of these subgroups revealed another QTL in the middle of GGA5 and increased the significance of the distal GGA5 AF QTL, thereby refining its localization. The second approach targeted the genes correlated with the AF trait and regulated by the GGA5 AF QTL region. Five of the 660 genes were considered as being controlled either by the AF QTL mutation itself or by a mutation close to it; one having a function related to lipid metabolism (HMGCS1). In addition, a QTL analysis with a multiple trait model combining this 5 gene-set and AF allowed us to refine the QTL region. The third approach was to use these 5 transcriptome profiles to predict the paternal Q versus q AF QTL mutation for each recombinant offspring and then refine the localization of the QTL from 31 cM (100 genes) at a most probable location confidence interval of 7 cM (12 genes) after determining the recombination breakpoints, an interval consistent with the reductions obtained by the two other approaches.

Conclusion

The results showed the feasibility and efficacy of the three strategies used, the first revealing a QTL undetected using the whole population, the second providing functional information about a QTL region through genes related to the trait and controlled by this region (HMGCS1), the third could drastically refine a QTL region.

Comparison of three microarray probe annotation pipelines: differences in strategies and their effect on downstream analysis

Background

Reliable annotation linking oligonucleotide probes to target genes is essential for functional biological analysis of microarray experiments. We used the IMAD, OligoRAP and sigReannot pipelines to update the annotation for the ARK-Genomics Chicken 20 K array as part of a joined EADGENE/SABRE workshop. In this manuscript we compare their annotation strategies and results. Furthermore, we analyse the effect of differences in updated annotation on functional analysis for an experiment involving Eimeria infected chickens and finally we propose guidelines for optimal annotation strategies.

Results

IMAD, OligoRAP and sigReannot update both annotation and estimated target specificity. The 3 pipelines can assign oligos to target specificity categories although with varying degrees of resolution. Target specificity is judged based on the amount and type of oligo versus target-gene alignments (hits), which are determined by filter thresholds that users can adjust based on their experimental conditions. Linking oligos to annotation on the other hand is based on rigid rules, which differ between pipelines.

For 52.7% of the oligos from a subset selected for in depth comparison all pipelines linked to one or more Ensembl genes with consensus on 44.0%. In 31.0% of the cases none of the pipelines could assign an Ensembl gene to an oligo and for the remaining 16.3% the coverage differed between pipelines. Differences in updated annotation were mainly due to different thresholds for hybridisation potential filtering of oligo versus target-gene alignments and different policies for expanding annotation using indirect links. The differences in updated annotation packages had a significant effect on GO term enrichment analysis with consensus on only 67.2% of the enriched terms.

Conclusion

In addition to flexible thresholds to determine target specificity, annotation tools should provide metadata describing the relationships between oligos and the annotation assigned to them. These relationships can then be used to judge the varying degrees of reliability allowing users to fine-tune the balance between reliability and coverage. This is important as it can have a significant effect on functional microarray analysis as exemplified by the lack of consensus on almost one third of the terms found with GO term enrichment analysis based on updated IMAD, OligoRAP or sigReannot annotation.

Methods for interpreting lists of affected genes obtained in a DNA microarray experiment

BACKGROUND

The aim of this paper was to describe and compare the methods used and the results obtained by the participants in a joint EADGENE (European Animal Disease Genomic Network of Excellence) and SABRE (Cutting Edge Genomics for Sustainable Animal Breeding) workshop focusing on post analysis of microarray data. The participating groups were provided with identical lists of microarray probes, including test statistics for three different contrasts, and the normalised log-ratios for each array, to be used as the starting point for interpreting the affected probes. The data originated from a microarray experiment conducted to study the host reactions in broilers occurring shortly after a secondary challenge with either a homologous or heterologous species of Eimeria.

RESULTS

Several conceptually different analytical approaches, using both commercial and public available software, were applied by the participating groups. The following tools were used: Ingenuity Pathway Analysis, MAPPFinder, LIMMA, GOstats, GOEAST, GOTM, Globaltest, TopGO, ArrayUnlock, Pathway Studio, GIST and AnnotationDbi. The main focus of the approaches was to utilise the relation between probes/genes and their gene ontology and pathways to interpret the affected probes/genes. The lack of a well-annotated chicken genome did though limit the possibilities to fully explore the tools. The main results from these analyses showed that the biological interpretation is highly dependent on the statistical method used but that some common biological conclusions could be reached.

CONCLUSION

It is highly recommended to test different analytical methods on the same data set and compare the results to obtain a reliable biological interpretation of the affected genes in a DNA microarray experiment.