Han Rauwerda

Potential new antiepileptogenic targets indicated by microarray analysis in a rat model for temporal lobe epilepsy.

To get insight into the mechanisms that may lead to progression of temporal lobe epilepsy, we investigated gene expression during epileptogenesis in the rat. RNA was obtained from three different brain regions [CA3, entorhinal cortex (EC), and cerebellum (CB)] at three different time points after electrically induced status epilepticus (SE): acute phase [group D (1 d)], latent period [group W (1 week)], and chronic epileptic period [group M (3–4 months)]. A group that was stimulated but that had not experienced SE and later epilepsy was also included (group nS). Gene expression analysis was performed using the Affymetrix Gene Chip System (RAE230A). We used GENMAPP and Gene Ontology to identify global biological trends in gene expression data. The immune response was the most prominent process changed during all three phases of epileptogenesis. Synaptic transmission was a downregulated process during the acute and latent phases. GABA receptor subunits involved in tonic inhibition were persistently downregulated. These changes were observed mostly in both CA3 and EC but not in CB. Rats that were stimulated but that did not develop spontaneous seizures later on had also some changes in gene expression, but this was not reflected in a significant change of a biological process. These data suggest that the targeting of specific genes that are involved in these biological processes may be a promising strategy to slow down or prevent the progression of epilepsy. Especially genes related to the immune response, such as complement factors, interleukins, and genes related to prostaglandin synthesis and coagulation pathway may be interesting targets.

Analysis of temporal gene expression during Bacillus subtilis spore germination and outgrowth.

Bacillus subtilis forms dormant spores upon nutrient depletion. Under favorable environmental conditions, the spore breaks its dormancy and resumes growth in a process called spore germination and outgrowth. To elucidate the physiological processes that occur during the transition of the dormant spore to an actively growing vegetative cell, we studied this process in a time-dependent manner by a combination of microscopy, analysis of extracellular metabolites, and a genome-wide analysis of transcription. The results indicate the presence of abundant levels of late sporulation transcripts in dormant spores. In addition, the results suggest the existence of a complex and well-regulated spore outgrowth program, involving the temporal expression of at least 30% of the B. subtilis genome.

Dynamic changes of proteases and protease inhibitors revealed by microarray analysis in CA3 and entorhinal cortex during epileptogenesis in the rat.

Summary:  We investigated expression of genes involved in the proteolytic process during epileptogenesis in a rat model of temporal lobe epilepsy (TLE). In a previous microarray study we found prominent activation of this process, which reached highest expression during the acute and latent phase (1 week after SE) in CA3 and entorhinal cortex (EC). Detailed analysis shows differences in dynamics of the changes of several protease genes such as cathepsins, caspases, matrix metalloproteinases, and plasminogen activators. Most genes were acutely upregulated while others were mainly activated during the latent phase. Interestingly several proteolytic genes were still elevated in the chronic epileptic phase. Various protease inhibitors followed a similar time course. The identification of changes in the activation of genes involved in proteolysis at critical phases during epileptogenesis could point to potential time specific targets for intervention. The fact that several proteolytic genes were still activated in the chronic epileptic phase makes them interesting candidates to modify and slow down seizure progression.

Analysing Scientific Workflows: Why Workflows Not Only Connect Web Services

Life science workflow systems are developed to help life scientists to conveniently connect various programs and web services. In practice however, much time is spent on data conversion, because web services provided by different organisations use different data formats. We have analysed all the Taverna workflows available at the my Experiment web site on December 11, 2008. Our analysis of the tasks in these workflows shows several noticeable aspects: their number ranges from 1 to 70 tasks per workflow; 18% of the workflows consist of a single task.Of the tasks used are 22% web services; local services, i.e. tasks executed by the workflow system itself, are very popular and cover 57% of tasks; tasks implemented by the workflow designer, scripting tasks, are is also used often (14%). Our analysis shows that over 30\% of tasks are related to data conversion.

Salvaging Affymetrix probes after probe-level re-annotation

BackgroundAffymetrix GeneChips can be re-annotated at the probe-level by breaking up the original probe-sets and recomposing new probe-sets based on up-to-date genomic knowledge, such as available in Entrez Gene. This results in custom Chip Description Files (CDF). Using these custom CDFs improves the quality of the data and thus the results of related gene expression studies. However, 44–71% of the probes on a GeneChip are lost in this re-annotation process. Although generally aimed at less known genes, losing these probes obviously means a substantial loss of expensive experiment data. Biologists are therefore very reluctant to adopt this approach.FindingsWe aimed to re-introduce the non-affected Affymetrix probe-sets after these re-annotation procedures. For this, we developed an algorithm (CDF-Merger) and applied it to standard Affymetrix CDFs and custom Brainarray CDFs to obtain Hybrid CDFs. Thus, salvaging lost Affymetrix probes with our CDF-Merger restored probe content up to 94%. Because the salvaged probes (up to 54% of the probe content on the arrays) represent less-reliable probe-sets, we made the origin of all probe-set definitions traceable, so biologists can choose at any time in their analyses, which subset of probe-sets they want to use.ConclusionThe availability of up-to-date Hybrid CDFs plus R environment allows for easy implementation of our approach.

Transcriptome analysis of Traf6 function in the innate immune response of zebrafish embryos

TRAF6 is a key player at the cross-roads of development and immunity. The analysis of its in vivo molecular function is a great challenge since severe developmental defects and early lethality caused by Traf6 deficiency in knock-out mice interfere with analyses of the immune response. In this study we have used a new strategy to analyze the function of Traf6 in a zebrafish-Salmonella infectious disease model. In our approach the effect of a Traf6 translation-blocking morpholino was titrated down to avoid developmental defects and the response to infection under these conditions was studied using the combination of microarray analysis and whole transcriptome deep sequencing. Transcriptome profiling of the traf6 knock-down allowed the identification of a gene set whose responsiveness during infection is highly dependent on Traf6. Expression trend analysis based on the resulting datasets identified nine clusters of genes with characteristic transcription response profiles, demonstrating Traf6 has a dynamic role as a positive and negative regulator. Among the Traf6-dependent genes was a large set of well known anti-microbial and inflammatory genes. Additionally, we identified several genes which were not previously linked to a response to microbial infection, such as the fertility hormone gene gnrh2 and the DNA-damage regulated autophagy modulator 1 gene dram1. With the use of the zebrafish embryo model we have now analyzed the in vivo function of Traf6 in the innate immune response without interference of adaptive immunity.

RNA isolation method for single embryo transcriptome analysis in zebrafish

BackgroundTranscriptome analysis during embryogenesis usually requires pooling of embryos to obtain sufficient RNA. Hence, the measured levels of gene-expression represent the average mRNA levels of pooled samples and the biological variation among individuals is confounded. This can irreversibly reduce the robustness, resolution, or expressiveness of the experiment. Therefore, we developed a robust method to isolate abundant high-quality RNA from individual embryos to perform single embryo transcriptome analyses using zebrafish as a model organism. Available methods for embryonic zebrafish RNA isolation minimally utilize ten embryos. Further downscaling of these methods to one embryo is practically not feasible.FindingsWe developed a single embryo RNA extraction method based on sample homogenization in liquid nitrogen, RNA extraction with phenol and column purification. Evaluation of this method showed that: the quality of the RNA was very good with an average RIN value of 8.3-8.9; the yield was always ≥ 200 ng RNA per embryo; the method was applicable to all stages of zebrafish embryogenesis; the success rate was almost 100%; and the extracted RNA performed excellent in microarray experiments in that the technical variation was much lower than the biological variation.ConclusionsPresented is a high-quality, robust RNA isolation method. Obtaining sufficient RNA from single embryos eliminates the necessity of sample pooling and its associated drawbacks. Although our RNA isolation method has been setup for transcriptome analysis in zebrafish, it can also be used for other model systems and other applications like (q)PCR and transcriptome sequencing.

The absence of Ser389 phosphorylation in p53 affects the basal gene expression level of many p53-dependent genes and alters the biphasic response to UV exposure in mouse embryonic fibroblasts

ABSTRACT Phosphorylation is important in p53-mediated DNA damage responses. After UV irradiation, p53 is phosphorylated specifically at murine residue Ser389. Phosphorylation mutant p53.S389A cells and mice show reduced apoptosis and compromised tumor suppression after UV irradiation. We investigated the underlying cellular processes by time-series analysis of UV-induced gene expression responses in wild-type, p53.S389A, and p53−/− mouse embryonic fibroblasts. The absence of p53.S389 phosphorylation already causes small endogenous gene expression changes for 2,253, mostly p53-dependent, genes. These genes showed basal gene expression levels intermediate to the wild type and p53−/−, possibly to readjust the p53 network. Overall, the p53.S389A mutation lifts p53-dependent gene repression to a level similar to that of p53−/− but has lesser effect on p53-dependently induced genes. In the wild type, the response of 6,058 genes to UV irradiation was strictly biphasic. The early stress response, from 0 to 3 h, results in the activation of processes to prevent the accumulation of DNA damage in cells, whereas the late response, from 12 to 24 h, relates more to reentering the cell cycle. Although the p53.S389A UV gene response was only subtly changed, many cellular processes were significantly affected. The early response was affected the most, and many cellular processes were phase-specifically lost, gained, or altered, e.g., induction of apoptosis, cell division, and DNA repair, respectively. Altogether, p53.S389 phosphorylation seems essential for many p53 target genes and p53-dependent processes.

WRKY40 and WRKY6 act downstream of the green leaf volatile E‐2‐hexenal in Arabidopsis

Plants are known to be responsive to volatiles, but knowledge about the molecular players involved in transducing their perception remains scarce. We study the response of Arabidopsis thaliana to E-2-hexenal, one of the green leaf volatiles (GLV) that is produced upon wounding, herbivory or infection with pathogens. We have taken a transcriptomics approach to identify genes that are induced by E-2-hexenal, but not by defence hormones or other GLVs. Furthermore, by studying the promoters of early E-2-hexenal-induced genes we determined that the only statistically enriched cis-element was the W-box motif. Since members of the plant-specific family of WRKY transcription factors act in trans on this cis-element, we focused on WRKY6, 40 and 53 that were most strongly induced by E-2-hexenal. Root elongation of Arabidopsis seedlings of the wrky40 wrky6 double mutant was much less inhibited than in wt plants, similar to the E-2-hexenal-responsive mutant her1, which is perturbed in γ-amino butyric acid (GABA) metabolism. The induction of several of the E-2-hexenal-specific genes was much higher in the wrky40, wrky6 or wrky40 wrky6 mutants, including GAD4, a glutamate decarboxylase that catalyzes the formation of GABA from glutamate. In conclusion, WRKY6 and 40 seem to act as important players transducing E-2-hexenal perception.

Improving small RNA-seq by using a synthetic spike-in set for size-range quality control together with a set for data normalization

There is an increasing interest in complementing RNA-seq experiments with small-RNA (sRNA) expression data to obtain a comprehensive view of a transcriptome. Currently, two main experimental challenges concerning sRNA-seq exist: how to check the size distribution of isolated sRNAs, given the sensitive size-selection steps in the protocol; and how to normalize data between samples, given the low complexity of sRNA types. We here present two separate sets of synthetic RNA spike-ins for monitoring size-selection and for performing data normalization in sRNA-seq. The size-range quality control (SRQC) spike-in set, consisting of 11 oligoribonucleotides (10–70 nucleotides), was tested by intentionally altering the size-selection protocol and verified via several comparative experiments. We demonstrate that the SRQC set is useful to reproducibly track down biases in the size-selection in sRNA-seq. The external reference for data-normalization (ERDN) spike-in set, consisting of 19 oligoribonucleotides, was developed for sample-to-sample normalization in differential-expression analysis of sRNA-seq data. Testing and applying the ERDN set showed that it can reproducibly detect differential expression over a dynamic range of 218. Hence, biological variation in sRNA composition and content between samples is preserved while technical variation is effectively minimized. Together, both spike-in sets can significantly improve the technical reproducibility of sRNA-seq.