Peter Sykacek

The transcriptome of syncytia induced by the cyst nematode Heterodera schachtii in Arabidopsis roots

Arabidopsis thaliana is a host for the sugar beet cyst nematode Heterodera schachtii. Juvenile nematodes invade the roots and induce the development of a syncytium, which functions as a feeding site for the nematode. Here, we report on the transcriptome of syncytia induced in the roots of Arabidopsis. Microaspiration was employed to harvest pure syncytium material, which was then used to prepare RNA for hybridization to Affymetrix GeneChips. Initial data analysis showed that the gene expression in syncytia at 5 and 15 days post-infection did not differ greatly, and so both time points were compared together with control roots. Out of a total of 21 138 genes, 18.4% (3893) had a higher expression level and 15.8% (3338) had a lower expression level in syncytia, as compared with control roots, using a multiple-testing corrected false discovery rate of below 5%. A gene ontology (GO) analysis of up- and downregulated genes showed that categories related to high metabolic activity were preferentially upregulated. A principal component analysis was applied to compare the transcriptome of syncytia with the transcriptome of different Arabidopsis organs (obtained by the AtGenExpress project), and with specific root tissues. This analysis revealed that syncytia are transcriptionally clearly different from roots (and all other organs), as well as from other root tissues.

Detecting and correcting systematic variation in large-scale RNA sequencing data

High-throughput RNA sequencing (RNA-seq) enables comprehensive scans of entire transcriptomes, but best practices for analyzing RNA-seq data have not been fully defined, particularly for data collected with multiple sequencing platforms or at multiple sites. Here we used standardized RNA samples with built-in controls to examine sources of error in large-scale RNA-seq studies and their impact on the detection of differentially expressed genes (DEGs). Analysis of variations in guanine-cytosine content, gene coverage, sequencing error rate and insert size allowed identification of decreased reproducibility across sites. Moreover, commonly used methods for normalization (cqn, EDASeq, RUV2, sva, PEER) varied in their ability to remove these systematic biases, depending on sample complexity and initial data quality. Normalization methods that combine data from genes across sites are strongly recommended to identify and remove site-specific effects and can substantially improve RNA-seq studies.

Bacterial protein signals are associated with Crohn’s disease

Objective No Crohn’s disease (CD) molecular maker has advanced to clinical use, and independent lines of evidence support a central role of the gut microbial community in CD. Here we explore the feasibility of extracting bacterial protein signals relevant to CD, by interrogating myriads of intestinal bacterial proteomes from a small number of patients and healthy controls. Design We first developed and validated a workflow—including extraction of microbial communities, two-dimensional difference gel electrophoresis (2D-DIGE), and LC-MS/MS—to discover protein signals from CD-associated gut microbial communities. Then we used selected reaction monitoring (SRM) to confirm a set of candidates. In parallel, we used 16S rRNA gene sequencing for an integrated analysis of gut ecosystem structure and functions. Results Our 2D-DIGE-based discovery approach revealed an imbalance of intestinal bacterial functions in CD. Many proteins, largely derived from Bacteroides species, were over-represented, while under-represented proteins were mostly from Firmicutes and some Prevotella members. Most overabundant proteins could be confirmed using SRM. They correspond to functions allowing opportunistic pathogens to colonise the mucus layers, breach the host barriers and invade the mucosae, which could still be aggravated by decreased host-derived pancreatic zymogen granule membrane protein GP2 in CD patients. Moreover, although the abundance of most protein groups reflected that of related bacterial populations, we found a specific independent regulation of bacteria-derived cell envelope proteins. Conclusions This study provides the first evidence that quantifiable bacterial protein signals are associated with CD, which can have a profound impact on future molecular diagnosis.

Model-based probe set optimization for high-performance microarrays

A major challenge in microarray design is the selection of highly specific oligonucleotide probes for all targeted genes of interest, while maintaining thermodynamic uniformity at the hybridization temperature. We introduce a novel microarray design framework (Thermodynamic Model-based Oligo Design Optimizer, TherMODO) that for the first time incorporates a number of advanced modelling features: (i) A model of position-dependent labelling effects that is quantitatively derived from experiment. (ii) Multi-state thermodynamic hybridization models of probe binding behaviour, including potential cross-hybridization reactions. (iii) A fast calibrated sequence-similarity-based heuristic for cross-hybridization prediction supporting large-scale designs. (iv) A novel compound score formulation for the integrated assessment of multiple probe design objectives. In contrast to a greedy search for probes meeting parameter thresholds, this approach permits an optimization at the probe set level and facilitates the selection of highly specific probe candidates while maintaining probe set uniformity. (v) Lastly, a flexible target grouping structure allows easy adaptation of the pipeline to a variety of microarray application scenarios. The algorithm and features are discussed and demonstrated on actual design runs. Source code is available on request.

A friendly statistics package for microarray analysis

SUMMARYnThe friendly statistics package for microarray analysis (FSPMA) is a tool that aims to fill the gap between simple to use and powerful analysis. FSPMA is a platform-independent R-package that allows efficient exploration of microarray data without the need for computer programming. Analysis is based on a mixed model ANOVA library (YASMA) that was extended to allow more flexible comparisons and other useful operations like k nearest neighbour imputing and spike-based normalization. Processing is controlled by a definition file that specifies all the steps necessary to derive analysis results from quantified microarray data. In addition to providing analysis without programming, the definition file also serves as exact documentation of all the analysis steps.nnnAVAILABILITYnThe library is available under GPL 2 license and, together with additional information, provided at http://www.ccbi.cam.ac.uk/software/psyk/software.html#fspma

The impact of quantitative optimization of hybridization conditions on gene expression analysis

BackgroundWith the growing availability of entire genome sequences, an increasing number of scientists can exploit oligonucleotide microarrays for genome-scale expression studies. While probe-design is a major research area, relatively little work has been reported on the optimization of microarray protocols.ResultsAs shown in this study, suboptimal conditions can have considerable impact on biologically relevant observations. For example, deviation from the optimal temperature by one degree Celsius lead to a loss of up to 44% of differentially expressed genes identified. While genes from thousands of Gene Ontology categories were affected, transcription factors and other low-copy-number regulators were disproportionately lost. Calibrated protocols are thus required in order to take full advantage of the large dynamic range of microarrays.For an objective optimization of protocols we introduce an approach that maximizes the amount of information obtained per experiment. A comparison of two typical samples is sufficient for this calibration. We can ensure, however, that optimization results are independent of the samples and the specific measures used for calibration. Both simulations and spike-in experiments confirmed an unbiased determination of generally optimal experimental conditions.ConclusionsWell calibrated hybridization conditions are thus easily achieved and necessary for the efficient detection of differential expression. They are essential for the sensitive pro filing of low-copy-number molecules. This is particularly critical for studies of transcription factor expression, or the inference and study of regulatory networks.

RNA Interference in Ageing Research – A Mini-Review

Background: The search for genetic mechanisms affecting life-span and ageing represents an important part of ageing research, especially since the discovery of single-gene mutations with dramatic effects on these traits. Due to its relative ease of use and its power to specifically target arbitrary genes, RNA interference (RNAi) has rapidly been adopted as a technique for silencing gene expression. The feasibility of genome-wide RNAi screens potentially much simplifies the identification of novel ageing-related genes. Objective: In a review of applications of RNAi in ageing research with a focus on the model organisms Caenorhabditis elegans and Drosophila melanogaster and discussing recent technical developments, we aim to highlight the current and future impact of this technology in the field. Method: We show how RNAi has successfully been used to complement classic mutant studies. Moreover, we discuss the novel opportunities and challenges of an application of RNAi in genome-wide screens in D. melanogaster, which has become possible with the recent availability of a comprehensive transgenic RNAi library for the fly. We highlight, in particular, how the flexible control of RNAi induction can support the study of dynamic processes like ageing through specific experiments and the development of matching computational methods. In an overview of complementary approaches we discuss the challenge of extracting insight from the high-dimensional measurement datasets that are required for the study of dynamic effects and interaction dependencies. Conclusion: RNAi has emerged as a powerful tool for the study of ageing, allowing the further characterization of the roles of specific genes in the ageing process as well as the efficient identification of new genes implicated. RNAi has contributed to our understanding of age-related diseases especially by making genes amenable to manipulation for which mutants were not easily available. Recent developments enable genome-wide screens with unprecedented temporal and spatial control of RNAi induction. Specific RNAi time-course experiments provide an opportunity for the analysis of high-resolution gene expression profiles capturing the dynamics of ageing-relevant processes and gene interactions. Research exploiting new avenues opened by the growing RNAi toolbox will considerably contribute to the next steps in researching the genetics of ageing and age-related diseases.

Bayesian modelling of shared gene function

MOTIVATIONnBiological assays are often carried out on tissues that contain many cell lineages and active pathways. Microarray data produced using such material therefore reflect superimpositions of biological processes. Analysing such data for shared gene function by means of well-matched assays may help to provide a better focus on specific cell types and processes. The identification of genes that behave similarly in different biological systems also has the potential to reveal new insights into preserved biological mechanisms.nnnRESULTSnIn this article, we propose a hierarchical Bayesian model allowing integrated analysis of several microarray data sets for shared gene function. Each gene is associated with an indicator variable that selects whether binary class labels are predicted from expression values or by a classifier which is common to all genes. Each indicator selects the component models for all involved data sets simultaneously. A quantitative measure of shared gene function is obtained by inferring a probability measure over these indicators. Through experiments on synthetic data, we illustrate potential advantages of this Bayesian approach over a standard method. A shared analysis of matched microarray experiments covering (a) a cycle of mouse mammary gland development and (b) the process of in vitro endothelial cell apoptosis is proposed as a biological gold standard. Several useful sanity checks are introduced during data analysis, and we confirm the prior biological belief that shared apoptosis events occur in both systems. We conclude that a Bayesian analysis for shared gene function has the potential to reveal new biological insights, unobtainable by other means.nnnAVAILABILITYnAn online supplement and MatLab code are available at http://www.sykacek.net/research.html#mcabf

Bayesian assignment of gene ontology terms to gene expression experiments

Motivation: Gene expression assays allow for genome scale analyses of molecular biological mechanisms. State-of-the-art data analysis provides lists of involved genes, either by calculating significance levels of mRNA abundance or by Bayesian assessments of gene activity. A common problem of such approaches is the difficulty of interpreting the biological implication of the resulting gene lists. This lead to an increased interest in methods for inferring high-level biological information. A common approach for representing high level information is by inferring gene ontology (GO) terms which may be attributed to the expression data experiment. Results: This article proposes a probabilistic model for GO term inference. Modelling assumes that gene annotations to GO terms are available and gene involvement in an experiment is represented by a posterior probabilities over gene-specific indicator variables. Such probability measures result from many Bayesian approaches for expression data analysis. The proposed model combines these indicator probabilities in a probabilistic fashion and provides a probabilistic GO term assignment as a result. Experiments on synthetic and microarray data suggest that advantages of the proposed probabilistic GO term inference over statistical test-based approaches are in particular evident for sparsely annotated GO terms and in situations of large uncertainty about gene activity. Provided that appropriate annotations exist, the proposed approach is easily applied to inferring other high level assignments like pathways. Availability: Source code under GPL license is available from the author. Contact: peter.sykacek@boku.ac.at

Root traits of European Vicia faba cultivars-Using machine learning to explore adaptations to agroclimatic conditions

Faba bean (Vicia faba L.) is an important source of protein, but breeding for increased yield stability and stress tolerance is hampered by the scarcity of phenotyping information. Because comparisons of cultivars adapted to different agroclimatic zones improve our understanding of stress tolerance mechanisms, the root architecture and morphology of 16 European faba bean cultivars were studied at maturity. Different machine learning (ML) approaches were tested in their usefulness to analyse trait variations between cultivars. A supervised, that is, hypothesis-driven, ML approach revealed that cultivars from Portugal feature greater and coarser but less frequent lateral roots at the top of the taproot, potentially enhancing water uptake from deeper soil horizons. Unsupervised clustering revealed that trait differences between northern and southern cultivars are not predominant but that two cultivar groups, independently from major and minor types, differ largely in overall root system size. Methodological guidelines on how to use powerful ML methods such as random forest models for enhancing the phenotypical exploration of plants are given.