The gut microbiota elicits a profound metabolic reorientation in the mouse jejunal mucosa during conventionalisation

Objective: Proper interactions between the intestinal mucosa, gut microbiota and nutrient flow are required to establish homoeostasis of the host. Since the proximal part of the small intestine is the first region where these interactions occur, and since most of the nutrient absorption occurs in the jejunum, it is important to understand the dynamics of metabolic responses of the mucosa in this intestinal region.Design: Germ-free mice aged 8-10 weeks were conventionalised with faecal microbiota, and responses of the jejunal mucosa to bacterial colonisation were followed over a 30-day time course. Combined transcriptome, histology, (1)H NMR metabonomics and microbiota phylogenetic profiling analyses were used.Results: The jejunal mucosa showed a two-phase response to the colonising microbiota. The acute-phase response, which had already started 1 day after conventionalisation, involved repression of the cell cycle and parts of the basal metabolism. The secondary-phase response, which was consolidated during conventionalisation (days 4-30), was characterised by a metabolic shift from an oxidative energy supply to anabolic metabolism, as inferred from the tissue transcriptome and metabonome changes. Detailed transcriptome analysis identified tissue transcriptional signatures for the dynamic control of the metabolic reorientation in the jejunum. The molecular components identified in the response signatures have known roles in human metabolic disorders, including insulin sensitivity and type 2 diabetes mellitus.Conclusion: This study elucidates the dynamic jejunal response to the microbiota and supports a prominent role for the jejunum in metabolic control, including glucose and energy homoeostasis. The molecular signatures of this process may help to find risk markers in the declining insulin sensitivity seen in human type 2 diabetes mellitus, for instance.

p90 ribosomal S6 kinases play a significant role in early gene regulation in the cardiomyocyte response to Gq protein-coupled receptor stimuli, endothelin-1 and α1-adrenergic receptor agonists

Extracellular signal-regulated kinases 1/2 (ERK1/2) and their substrates, p90 ribosomal S6 kinases (RSKs), phosphorylate different transcription factors, contributing differentially to transcriptomic profiles. In cardiomyocytes, ERK1/2 are required for >70% of the transcriptomic response to endothelin-1. Here, we investigated the role of RSKs in the transcriptomic responses to Gq protein-coupled receptor agonists, endothelin-1, phenylephrine (generic α1-adrenergic receptor agonist) and A61603 (α1A-adrenergic receptor selective). Phospho-ERK1/2 and phospho-RSKs appeared in cardiomyocyte nuclei within 2-3 min of stimulation (endothelin-1>a61603≈phenylephrine). All agonists increased nuclear RSK2, but only endothelin-1 increased nuclear RSK1 content. PD184352 (inhibits ERK1/2 activation) and BI-D1870 (inhibits RSKs) were used to dissect the contribution of RSKs to the endothelin-1-responsive transcriptome. Of 213 RNAs upregulated at 1 h, 51% required RSKs for upregulation whereas 29% required ERK1/2 but not RSKs. The transcriptomic response to phenylephrine overlapped with, but was not identical to, endothelin-1. As with endothelin-1, PD184352 inhibited upregulation of most phenylephrine-responsive transcripts, but the greater variation in effects of BI-D1870 suggests that differential RSK signalling influences global gene expression. A61603 induced similar changes in RNA expression in cardiomyocytes as phenylephrine, indicating that the signal was mediated largely through α1A-adrenergic receptors. A61603 also increased expression of immediate early genes in perfused adult rat hearts and, as in cardiomyocytes, upregulation of the majority of genes was inhibited by PD184352. PD184352 or BI-D1870 prevented the increased surface area induced by endothelin-1 in cardiomyocytes. Thus, RSKs play a significant role in regulating cardiomyocyte gene expression and hypertrophy in response to Gq protein-coupled receptor stimulation.

The early transcriptomic response to interleukin 1β and interleukin 33 in rat neonatal cardiomyocytes

In the heart, inflammatory cytokines including interleukin (IL) 1β are implicated in regulating adaptive and maladaptive changes, whereas IL33 negatively regulates cardiomyocyte hypertrophy and promotes cardioprotection. These agonists signal through a common co-receptor but, in cardiomyocytes, IL1β more potently activates mitogen-activated protein kinases and NFκB, pathways that regulate gene expression. We compared the effects of external application of IL1β and IL33 on the cardiomyocyte transcriptome. Neonatal rat cardiomyocytes were exposed to IL1β or IL33 (0.5, 1 or 2h). Transcriptomic profiles were determined using Affymetrix rat genome 230 2.0 microarrays and data were validated by quantitative PCR. IL1β induced significant changes in more RNAs than IL33 and, generally, to a greater degree. It also had a significantly greater effect in downregulating mRNAs and in regulating mRNAs associated with selected pathways. IL33 had a greater effect on a small, select group of specific transcripts. Thus, differences in intensity of intracellular signals can deliver qualitatively different responses. Quantitatively different responses in production of receptor agonists and transcription factors may contribute to qualitative differences at later times resulting in different phenotypic cellular responses.

Gut bacteria-host metabolic interplay during conventionalisation of the mouse germfree colon

The interplay between dietary nutrients, gut microbiota and mammalian host tissues of the gastrointestinal tract is recognised as highly relevant for host health. Combined transcriptome, metabonome and microbial profiling tools were employed to analyse the dynamic responses of germfree mouse colonic mucosa to colonisation by normal mouse microbiota (conventionalisation) at different time-points during 16 days. The colonising microbiota showed a shift from early (days 1 and 2) to later colonisers (days 8 and 16). The dynamic changes in the microbial community were rapidly reflected by the urine metabolic profiles (day 1) and at later stages (day 4 onward) by the colon mucosa transcriptome and metabolic profiles. Correlations of host transcriptomes, metabolite patterns and microbiota composition revealed associations between Bacilli and Proteobacteria, and differential expression of host genes involved in energy and anabolic metabolism. Differential gene expression correlated with scyllo- and myo-inositol, glutamine, glycine and alanine levels in colonic tissues during the time span of conventionalisation. Our combined time-resolved analyses may help to expand the understanding of host-microbe molecular interactions during the microbial establishment.

Error, reproducibility and sensitivity: a pipeline for data processing of Agilent oligonucleotide expression arrays

Background: Expression microarrays are increasingly used to obtain large scale transcriptomic information on a wide range of biological samples. Nevertheless, there is still much debate on the best ways to process data, to design experiments and analyse the output. Furthermore, many of the more sophisticated mathematical approaches to data analysis in the literature remain inaccessible to much of the biological research community. In this study we examine ways of extracting and analysing a large data set obtained using the Agilent long oligonucleotide transcriptomics platform, applied to a set of human macrophage and dendritic cell samples. Results: We describe and validate a series of data extraction, transformation and normalisation steps which are implemented via a new R function. Analysis of replicate normalised reference data demonstrate that intrarray variability is small (only around 2 of the mean log signal), while interarray variability from replicate array measurements has a standard deviation (SD) of around 0.5 log(2) units (6 of mean). The common practise of working with ratios of Cy5/Cy3 signal offers little further improvement in terms of reducing error. Comparison to expression data obtained using Arabidopsis samples demonstrates that the large number of genes in each sample showing a low level of transcription reflect the real complexity of the cellular transcriptome. Multidimensional scaling is used to show that the processed data identifies an underlying structure which reflect some of the key biological variables which define the data set. This structure is robust, allowing reliable comparison of samples collected over a number of years and collected by a variety of operators. Conclusions: This study outlines a robust and easily implemented pipeline for extracting, transforming normalising and visualising transcriptomic array data from Agilent expression platform. The analysis is used to obtain quantitative estimates of the SD arising from experimental (non biological) intra- and interarray variability, and for a lower threshold for determining whether an individual gene is expressed. The study provides a reliable basis for further more extensive studies of the systems biology of eukaryotic cells.

A comparison of the Thlaspi caerulescens and Thlaspi arvense shoot transcriptomes

Whole-genome transcriptome profiling is revealing how biological systems are regulated at the transcriptional level. This study reports the development of a robust method to profile and compare the transcriptomes of two nonmodel plant species, Thlaspi caerulescens, a zinc (Zn) hyperaccumulator, and Thlaspi arvense, a nonhyperaccumulator, using Affymetrix Arabidopsis thaliana ATH1-121501 GeneChip (R) arrays (Affymetrix, Santa Clara, CA, USA). Transcript abundance was quantified in the shoots of agar- and compost-grown plants of both species. Analyses were optimized using a genomic DNA (gDNA)-based probe-selection strategy based on the hybridization efficiency of Thlaspi gDNA with corresponding A. thaliana probes. In silico alignments of GeneChip (R) probes with Thlaspi gene sequences, and quantitative real-time PCR, confirmed the validity of this approach. Approximately 5000 genes were differentially expressed in the shoots of T. caerulescens compared with T. arvense, including genes involved in Zn transport and compartmentalization. Future functional analyses of genes identified as differentially expressed in the shoots of these closely related species will improve our understanding of the molecular mechanisms of Zn hyperaccumulation.

Using genomic DNA-based probe-selection to improve the sensitivity of high-density oligonucleotide arrays when applied to heterologous species

High-density oligonucleotide (oligo) arrays are a powerful tool for transcript profiling. Arrays based on GeneChip® technology are amongst the most widely used, although GeneChip® arrays are currently available for only a small number of plant and animal species. Thus, we have developed a method to improve the sensitivity of high-density oligonucleotide arrays when applied to heterologous species and tested the method by analysing the transcriptome of Brassica oleracea L., a species for which no GeneChip® array is available, using a GeneChip® array designed for Arabidopsis thaliana (L.) Heynh. Genomic DNA from B. oleracea was labelled and hybridised to the ATH1-121501 GeneChip® array. Arabidopsis thaliana probe-pairs that hybridised to the B. oleracea genomic DNA on the basis of the perfect-match (PM) probe signal were then selected for subsequent B. oleracea transcriptome analysis using a .cel file parser script to generate probe mask files. The transcriptional response of B. oleracea to a mineral nutrient (phosphorus; P) stress was quantified using probe mask files generated for a wide range of gDNA hybridisation intensity thresholds. An example probe mask file generated with a gDNA hybridisation intensity threshold of 400 removed > 68 % of the available PM probes from the analysis but retained >96 % of available A. thaliana probe-sets. Ninety-nine of these genes were then identified as significantly regulated under P stress in B. oleracea, including the homologues of P stress responsive genes in A. thaliana. Increasing the gDNA hybridisation intensity thresholds up to 500 for probe-selection increased the sensitivity of the GeneChip® array to detect regulation of gene expression in B. oleracea under P stress by up to 13-fold. Our open-source software to create probe mask files is freely available http://affymetrix.arabidopsis.info/xspecies/ webcite and may be used to facilitate transcriptomic analyses of a wide range of plant and animal species in the absence of custom arrays.

Genes driving potato tuber initiation and growth: identification based on transcriptional changes using the POCI array

The increasing amount of available expressed gene sequence data makes whole-transcriptome analysis of certain crop species possible. Potato currently has the second largest number of publicly available expressed sequence tag (EST) sequences among the Solanaceae. Most of these ESTs, plus other proprietary sequences, were combined and used to generate a unigene assembly. The set of 246,182 sequences produced 46,345 unigenes, which were used to design a 44K 60-mer oligo array (Potato Oligo Chip Initiative: POCI). In this study, we attempt to identify genes controlling and driving the process of tuber initiation and growth by implementing large-scale transcriptional changes using the newly developed POCI array. Major gene expression profiles could be identified exhibiting differential expression at key developmental stages. These profiles were associated with functional roles in cell division and growth. A subset of genes involved in the regulation of the cell cycle, based on their Gene Ontology classification, exhibit a clear transient upregulation at tuber onset indicating increased cell division during these stages. The POCI array allows the study of potato gene expression on a much broader level than previously possible and will greatly enhance analysis of transcriptional control mechanisms in a wide range of potato research areas. POCI sequence and annotation data are publicly available through the POCI database (http://pgrc.ipk-gatersleben.de/poci).

A Brassica exon array for whole-transcript gene expression profiling

Affymetrix GeneChip (R) arrays are used widely to study transcriptional changes in response to developmental and environmental stimuli. GeneChip (R) arrays comprise multiple 25-mer oligonucleotide probes per gene and retain certain advantages over direct sequencing. For plants, there are several public GeneChip (R) arrays whose probes are localised primarily in 39 exons. Plant whole-transcript (WT) GeneChip (R) arrays are not yet publicly available, although WT resolution is needed to study complex crop genomes such as Brassica, which are typified by segmental duplications containing paralogous genes and/or allopolyploidy. Available sequence data were sampled from the Brassica A and C genomes, and 142,997 gene models identified. The assembled gene models were then used to establish a comprehensive public WT exon array for transcriptomics studies. The Affymetrix GeneChip (R) Brassica Exon 1.0 ST Array is a 5 mu M feature size array, containing 2.4 million 25-base oligonucleotide probes representing 135,201 gene models, with 15 probes per gene distributed among exons. Discrimination of the gene models was based on an E-value cut-off of 1E(-5), with <= 98 sequence identity. The 135 k Brassica Exon Array was validated by quantifying transcriptome differences between leaf and root tissue from a reference Brassica rapa line (R-o-18), and categorisation by Gene Ontologies (GO) based on gene orthology with Arabidopsis thaliana. Technical validation involved comparison of the exon array with a 60-mer array platform using the same starting RNA samples. The 135 k Brassica Exon Array is a robust platform. All data relating to the array design and probe identities are available in the public domain and are curated within the BrassEnsembl genome viewer at http://www.brassica.info/BrassEnsembl/index.html.

A SEPALLATA gene is involved in the development and ripening of strawberry (Fragaria X ananassa Duch.) fruit, a non-climacteric tissue

Climacteric and non-climacteric fruits have traditionally been viewed as representing two distinct programmes of ripening associated with differential respiration and ethylene hormone effects. In climacteric fruits, such as tomato and banana, the ripening process is marked by increased respiration and is induced and co-ordinated by ethylene, while in non-climacteric fruits, such as strawberry and grape, it is controlled by an ethylene-independent process with little change in respiration rate. The two contrasting mechanisms, however, both lead to texture, colour, and flavour changes that probably reflect some common programmes of regulatory control. It has been shown that a SEPALLATA(SEP)4-like gene is necessary for normal ripening in tomato. It has been demonstrated here that silencing a fruit-related SEP1/2-like (FaMADS9) gene in strawberry leads to the inhibition of normal development and ripening in the petal, achene, and receptacle tissues. In addition, analysis of transcriptome profiles reveals pleiotropic effects of FaMADS9 on fruit development and ripening-related gene expression. It is concluded that SEP genes play a central role in the developmental regulation of ripening in both climacteric and non-climacteric fruits. These findings provide important information to extend the molecular control of ripening in a non-climacteric fruit beyond the limited genetic and cultural options currently available.

Dysregulation of REST-regulated coding and non-coding RNAs in a cellular model of Huntington's disease

Huntingtin (Htt) protein interacts with many transcriptional regulators, with widespread disruption to the transcriptome in Huntington's disease (HD) brought about by altered interactions with the mutant Htt (muHtt) protein. Repressor Element-1 Silencing Transcription Factor (REST) is a repressor whose association with Htt in the cytoplasm is disrupted in HD, leading to increased nuclear REST and concomitant repression of several neuronal-specific genes, including brain-derived neurotrophic factor (Bdnf). Here, we explored a wide set of HD dysregulated genes to identify direct REST targets whose expression is altered in a cellular model of HD but that can be rescued by knock-down of REST activity. We found many direct REST target genes encoding proteins important for nervous system development, including a cohort involved in synaptic transmission, at least two of which can be rescued at the protein level by REST knock-down. We also identified several microRNAs (miRNAs) whose aberrant repression is directly mediated by REST, including miR-137, which has not previously been shown to be a direct REST target in mouse. These data provide evidence of the contribution of inappropriate REST-mediated transcriptional repression to the widespread changes in coding and non-coding gene expression in a cellular model of HD that may affect normal neuronal function and survival.

Transcriptional dysregulation of coding and non-coding genes in cellular models of Huntington's disease

HD (Huntington's disease) is a late onset heritable neurodegenerative disorder that is characterized by neuronal dysfunction and death, particularly in the cerebral cortex and medium spiny neurons of the striatum. This is followed by progressive chorea, dementia and emotional dysfunction, eventually resulting in death. HD is caused by an expanded CAG repeat in the first exon of the HD gene that results in an abnormally elongated polyQ (polyglutamine) tract in its protein product, Htt (Huntingtin). Wild-type Htt is largely cytoplasmic; however, in HD, proteolytic N-terminal fragments of Htt form insoluble deposits in both the cytoplasm and nucleus, provoking the idea that mutHtt (mutant Htt) causes transcriptional dysfunction. While a number of specific transcription factors and co-factors have been proposed as mediators of mutHtt toxicity, the causal relationship between these Htt/transcription factor interactions and HD pathology remains unknown. Previous work has highlighted REST [RE1 (repressor element 1)-silencing transcription factor] as one such transcription factor. REST is a master regulator of neuronal genes, repressing their expression. Many of its direct target genes are known or suspected to have a role in HD pathogenesis, including BDNF (brain-derived neurotrophic factor). Recent evidence has also shown that REST regulates transcription of regulatory miRNAs (microRNAs), many of which are known to regulate neuronal gene expression and are dysregulated in HD. Thus repression of miRNAs constitutes a second, indirect mechanism by which REST can alter the neuronal transcriptome in HD. We will describe the evidence that disruption to the REST regulon brought about by a loss of interaction between REST and mutHtt may be a key contributory factor in the widespread dysregulation of gene expression in HD.

The global effect of follicle-stimulating hormone and tumour necrosis factor α on gene expression in cultured bovine ovarian granulosa cells

Background Oocytes mature in ovarian follicles surrounded by granulosa cells. During follicle growth, granulosa cells replicate and secrete hormones, particularly steroids close to ovulation. However, most follicles cease growing and undergo atresia or regression instead of ovulating. To investigate the effects of stimulatory (follicle-stimulating hormone; FSH) and inhibitory (tumour necrosis factor alpha; TNFα) factors on the granulosa cell transcriptome, bovine ovaries were obtained from a local abattoir and pools of granulosa cells were cultured in vitro for six days under defined serum-free conditions with treatments present on days 3–6. Initially dose–response experiments (n = 4) were performed to determine the optimal concentrations of FSH (0.33 ng/ml) and TNFα (10 ng/ml) to be used for the microarray experiments. For array experiments cells were cultured under control conditions, with FSH, with TNFα, or with FSH plus TNFα (n = 4 per group) and RNA was harvested for microarray analyses. Results Statistical analysis showed primary clustering of the arrays into two groups, control/FSH and TNFα/TNFα plus FSH. The effect of TNFα on gene expression dominated that of FSH, with substantially more genes differentially regulated, and the pathways and genes regulated by TNFα being similar to those of FSH plus TNFα treatment. TNFα treatment reduced the endocrine activity of granulosa cells with reductions in expression of FST, INHA, INBA and AMH. The top-ranked canonical pathways and GO biological terms for the TNFα treatments included antigen presentation, inflammatory response and other pathways indicative of innate immune function and fibrosis. The two most significant networks also reflect this, containing molecules which are present in the canonical pathways of hepatic fibrosis/hepatic stellate cell activation and transforming growth factor β signalling, and these were up regulated. Upstream regulator analyses also predicted TNF, interferons γ and β1 and interleukin 1β. Conclusions In vitro, the transcriptome of granulosa cells responded minimally to FSH compared with the response to TNFα. The response to TNFα indicated an active process akin to tissue remodelling as would occur upon atresia. Additionally there was reduction in endocrine function and induction of an inflammatory response to TNFα that displays features similar to immune cells.

Lunasin in cereal seeds: what is the origin?

Lunasin is a peptide from soybean seeds which has been demonstrated to have anticancer properties. It has also been reported in cereal seeds: wheat, rye, barley and Triticale. However, extensive searches of transcriptome and DNA sequence databases for wheat and other cereals have failed to identify sequences encoding either the lunasin peptide or a precursor protein. This raises the question of the origin of the lunasin reported in cereal grain.

Characterization of the lipoxygenase (LOX) gene family in the Chinese white pear (Pyrus bretschneideri) and comparison with other members of the Rosaceae

Background Lipoxygenases (LOXs), a type of non-haem iron-containing dioxygenase, are ubiquitous enzymes in plants and participate in the formation of fruit aroma which is a very important aspect of fruit quality. Amongst the various aroma volatiles, saturated and unsaturated alcohols and aldehydes provide the characteristic aroma of the fruit. These compounds are formed from unsaturated fatty acids through oxidation, pyrolysis and reduction steps. This biosynthetic pathway involves at least four enzymes, including LOX, the enzyme responsible for lipid oxidation. Although some studies have been conducted on the LOX gene family in several species including Arabidopsis, soybean, cucumber and apple, there is no information from pear; and the evolutionary history of this gene family in the Rosaceae is still not resolved. Results In this study we identified 107 LOX homologous genes from five Rosaceous species (Pyrus bretschneideri, Malus × domestica, Fragaria vesca, Prunus mume and Prunus persica); 23 of these sequences were from pear. By using structure analysis, phylogenic analysis and collinearity analysis, we identified variation in gene structure and revealed the phylogenetic evolutionary relationship of this gene family. Expression of certain pear LOX genes during fruit development was verified by analysis of transcriptome data. Conclusions 23 LOX genes were identified in pear and these genes were found to have undergone a duplication 30–45 MYA; most of these 23 genes are functional. Specific gene duplication was found on chromosome4 in the pear genome. Useful information was provided for future research on the evolutionary history and transgenic research on LOX genes.