Bjoern Usadel

Trimmomatic: a flexible trimmer for Illumina sequence data

Motivation: Although many next-generation sequencing (NGS) read preprocessing tools already existed, we could not find any tool or combination of tools that met our requirements in terms of flexibility, correct handling of paired-end data and high performance. We have developed Trimmomatic as a more flexible and efficient preprocessing tool, which could correctly handle paired-end data. Results: The value of NGS read preprocessing is demonstrated for both reference-based and reference-free tasks. Trimmomatic is shown to produce output that is at least competitive with, and in many cases superior to, that produced by other tools, in all scenarios tested. Availability and implementation: Trimmomatic is licensed under GPL V3. It is cross-platform (Java 1.5+ required) and available at http://www.usadellab.org/cms/index.php?page=trimmomatic Contact: usadel@bio1.rwth-aachen.de Supplementary information: Supplementary data are available at Bioinformatics online.

Starch as a major integrator in the regulation of plant growth

Rising demand for food and bioenergy makes it imperative to breed for increased crop yield. Vegetative plant growth could be driven by resource acquisition or developmental programs. Metabolite profiling in 94 Arabidopsis accessions revealed that biomass correlates negatively with many metabolites, especially starch. Starch accumulates in the light and is degraded at night to provide a sustained supply of carbon for growth. Multivariate analysis revealed that starch is an integrator of the overall metabolic response. We hypothesized that this reflects variation in a regulatory network that balances growth with the carbon supply. Transcript profiling in 21 accessions revealed coordinated changes of transcripts of more than 70 carbon-regulated genes and identified 2 genes (myo-inositol-1-phosphate synthase, a Kelch-domain protein) whose transcripts correlate with biomass. The impact of allelic variation at these 2 loci was shown by association mapping, identifying them as candidate lead genes with the potential to increase biomass production.

Integration of metabolite with transcript and enzyme activity profiling during diurnal cycles in Arabidopsis rosettes

BackgroundGenome-wide transcript profiling and analyses of enzyme activities from central carbon and nitrogen metabolism show that transcript levels undergo marked and rapid changes during diurnal cycles and after transfer to darkness, whereas changes in activities are smaller and delayed. In the starchless pgm mutant, where sugars are depleted every night, there are accentuated diurnal changes in transcript levels. Enzyme activities in this mutant do not show larger diurnal changes; instead, they shift towards the levels found in the wild type after several days of darkness. This indicates that enzyme activities change slowly, integrating the changes in transcript levels over several diurnal cycles.ResultsTo generalize this conclusion, 137 metabolites were profiled using gas and liquid chromatography coupled to mass spectroscopy. The amplitudes of the diurnal changes in metabolite levels in pgm were (with the exception of sugars) similar or smaller than in the wild type. The average levels shifted towards those found after several days of darkness in the wild type. Examples include increased levels of amino acids due to protein degradation, decreased levels of fatty acids, increased tocopherol and decreased myo-inositol. Many metabolite-transcript correlations were found and the proportion of transcripts correlated with sugars increased dramatically in the starchless mutant.ConclusionRapid diurnal changes in transcript levels are integrated over time to generate quasi-stable changes across large sectors of metabolism. This implies that correlations between metabolites and transcripts are due to regulation of gene expression by metabolites, rather than metabolites being changed as a consequence of a change in gene expression.

Global Transcript Levels Respond to Small Changes of the Carbon Status during Progressive Exhaustion of Carbohydrates in Arabidopsis Rosettes

The balance between the supply and utilization of carbon (C) changes continually. It has been proposed that plants respond in an acclimatory manner, modifying C utilization to minimize harmful periods of C depletion. This hypothesis predicts that signaling events are initiated by small changes in C status. We analyzed the global transcriptional response to a gradual depletion of C during the night and an extension of the night, where C becomes severely limiting from 4 h onward. The response was interpreted using published datasets for sugar, light, and circadian responses. Hundreds of C-responsive genes respond during the night and others very early in the extended night. Pathway analysis reveals that biosynthesis and cellular growth genes are repressed during the night and genes involved in catabolism are induced during the first hours of the extended night. The C response is amplified by an antagonistic interaction with the clock. Light signaling is attenuated during the 24-h light/dark cycle. A model was developed that uses the response of 22K genes during a circadian cycle and their responses to C and light to predict global transcriptional responses during diurnal cycles of wild-type and starchless pgm mutant plants and an extended night in wild-type plants. By identifying sets of genes that respond at different speeds and times during C depletion, our extended dataset and model aid the analysis of candidates for C signaling. This is illustrated for AKIN10 and four bZIP transcription factors, and sets of genes involved in trehalose signaling, protein turnover, and starch breakdown.

Regulatory Features Underlying Pollination-Dependent and -Independent Tomato Fruit Set Revealed by Transcript and Primary Metabolite Profiling

Indole Acetic Acid 9 (IAA9) is a negative auxin response regulator belonging to the Aux/IAA transcription factor gene family whose downregulation triggers fruit set before pollination, thus giving rise to parthenocarpy. In situ hybridization experiments revealed that a tissue-specific gradient of IAA9 expression is established during flower development, the release of which upon pollination triggers the initiation of fruit development. Comparative transcriptome and targeted metabolome analysis uncovered important features of the molecular events underlying pollination-induced and pollination-independent fruit set. Comprehensive transcriptomic profiling identified a high number of genes common to both types of fruit set, among which only a small subset are dependent on IAA9 regulation. The fine-tuning of Aux/IAA and ARF genes and the downregulation of TAG1 and TAGL6 MADS box genes are instrumental in triggering the fruit set program. Auxin and ethylene emerged as the most active signaling hormones involved in the flower-to-fruit transition. However, while these hormones affected only a small number of transcriptional events, dramatic shifts were observed at the metabolic and developmental levels. The activation of photosynthesis and sucrose metabolism-related genes is an integral regulatory component of fruit set process. The combined results allow a far greater comprehension of the regulatory and metabolic events controlling early fruit development both in the presence and absence of pollination/fertilization.

The Lipopolysaccharide of Sinorhizobium meliloti Suppresses Defense-Associated Gene Expression in Cell Cultures of the Host Plant Medicago truncatula

In the establishment of symbiosis between Medicago truncatula and the nitrogen-fixing bacterium Sinorhizobium meliloti, the lipopolysaccharide (LPS) of the microsymbiont plays an important role as a signal molecule. It has been shown in cell cultures that the LPS is able to suppress an elicitor-induced oxidative burst. To investigate the effect of S. meliloti LPS on defense-associated gene expression, a microarray experiment was performed. For evaluation of the M. truncatula microarray datasets, the software tool MapMan, which was initially developed for the visualization of Arabidopsis (Arabidopsis thaliana) datasets, was adapted by assigning Medicago genes to the ontology originally created for Arabidopsis. This allowed functional visualization of gene expression of M. truncatula suspension-cultured cells treated with invertase as an elicitor. A gene expression pattern characteristic of a defense response was observed. Concomitant treatment of M. truncatula suspension-cultured cells with invertase and S. meliloti LPS leads to a lower level of induction of defense-associated genes compared to induction rates in cells treated with invertase alone. This suppression of defense-associated transcriptional rearrangement affects genes induced as well as repressed by elicitation and acts on transcripts connected to virtually all kinds of cellular processes. This indicates that LPS of the symbiont not only suppresses fast defense responses as the oxidative burst, but also exerts long-term influences, including transcriptional adjustment to pathogen attack. These data indicate a role for LPS during infection of the plant by its symbiotic partner.

Diurnal changes of polysome loading track sucrose content in the rosette of wildtype Arabidopsis and the starchless pgm mutant

Ribosome loading is closely coordinated with the sucrose supply during diurnal cycles in Arabidopsis. Growth is driven by newly fixed carbon in the light, but at night it depends on reserves, like starch, that are laid down in the light. Unless plants coordinate their growth with diurnal changes in the carbon supply, they will experience acute carbon starvation during the night. Protein synthesis represents a major component of cellular growth. Polysome loading was investigated during the diurnal cycle, an extended night, and low CO2 in Arabidopsis (Arabidopsis thaliana) Columbia (Col-0) and in the starchless phosphoglucomutase (pgm) mutant. In Col-0, polysome loading was 60% to 70% in the light, 40% to 45% for much of the night, and less than 20% in an extended night, while in pgm, it fell to less than 25% early in the night. Quantification of ribosomal RNA species using quantitative reverse transcription-polymerase chain reaction revealed that polysome loading remained high for much of the night in the cytosol, was strongly light dependent in the plastid, and was always high in mitochondria. The rosette sucrose content correlated with overall and with cytosolic polysome loading. Ribosome abundance did not show significant diurnal changes. However, compared with Col-0, pgm had decreased and increased abundance of plastidic and mitochondrial ribosomes, respectively. Incorporation of label from 13CO2 into protein confirmed that protein synthesis continues at a diminished rate in the dark. Modeling revealed that a decrease in polysome loading at night is required to balance protein synthesis with the availability of carbon from starch breakdown. Costs are also reduced by using amino acids that accumulated in the previous light period. These results uncover a tight coordination of protein synthesis with the momentary supply of carbon.

A Topological Map of the Compartmentalized Arabidopsis thaliana Leaf Metabolome

Background The extensive subcellular compartmentalization of metabolites and metabolism in eukaryotic cells is widely acknowledged and represents a key factor of metabolic activity and functionality. In striking contrast, the knowledge of actual compartmental distribution of metabolites from experimental studies is surprisingly low. However, a precise knowledge of, possibly all, metabolites and their subcellular distributions remains a key prerequisite for the understanding of any cellular function. Methodology/Principal Findings Here we describe results for the subcellular distribution of 1,117 polar and 2,804 lipophilic mass spectrometric features associated to known and unknown compounds from leaves of the model plant Arabidopsis thaliana. Using an optimized non-aqueous fractionation protocol in conjunction with GC/MS- and LC/MS-based metabolite profiling, 81.5% of the metabolic data could be associated to one of three subcellular compartments: the cytosol (including the mitochondria), vacuole, or plastids. Statistical analysis using a marker-‘free’ approach revealed that 18.5% of these metabolites show intermediate distributions, which can either be explained by transport processes or by additional subcellular compartments. Conclusion/Significance Next to a functional and conceptual workflow for the efficient, highly resolved metabolite analysis of the fractionated Arabidopsis thaliana leaf metabolome, a detailed survey of the subcellular distribution of several compounds, in the graphical format of a topological map, is provided. This complex data set therefore does not only contain a rich repository of metabolic information, but due to thorough validation and testing by statistical methods, represents an initial step in the analysis of metabolite dynamics and fluxes within and between subcellular compartments.

Demethylesterification of Cell Wall Pectins in Arabidopsis Plays a Role in Seed Germination

The methylesterification status of cell wall homogalacturonans, mediated through the action of pectin methylesterases (PMEs), influences the biophysical properties of plant cell walls such as elasticity and porosity, important parameters for cell elongation and water uptake. The completion of seed germination requires cell wall extensibility changes in both the radicle itself and in the micropylar tissues surrounding the radicle. In wild-type seeds of Arabidopsis (Arabidopsis thaliana), PME activities peaked around the time of testa rupture but declined just before the completion of germination (endosperm weakening and rupture). We overexpressed an Arabidopsis PME inhibitor to investigate PME involvement in seed germination. Seeds of the resultant lines showed a denser methylesterification status of their cell wall homogalacturonans, but there were no changes in the neutral sugar and uronic acid composition of the cell walls. As compared with wild-type seeds, the PME activities of the overexpressing lines were greatly reduced throughout germination, and the low steady-state levels neither increased nor decreased. The most striking phenotype was a significantly faster rate of germination, which was not connected to altered testa rupture morphology but to alterations of the micropylar endosperm cells, evident by environmental scanning electron microscopy. The transgenic seeds also exhibited an apparent reduced sensitivity to abscisic acid with respect to its inhibitory effects on germination. We speculate that PME activity contributes to the temporal regulation of radicle emergence in endospermic seeds by altering the mechanical properties of the cell walls and thereby the balance between the two opposing forces of radicle elongation and mechanical resistance of the endosperm.

Algorithm-driven Artifacts in median polish summarization of Microarray data

BackgroundHigh-throughput measurement of transcript intensities using Affymetrix type oligonucleotide microarrays has produced a massive quantity of data during the last decade. Different preprocessing techniques exist to convert the raw signal intensities measured by these chips into gene expression estimates. Although these techniques have been widely benchmarked in the context of differential gene expression analysis, there are only few examples where their performance has been assessed in respect to coexpression-based studies such as sample classification.ResultsIn the present paper we benchmark the three most used normalization procedures (MAS5, RMA and GCRMA) in the context of inter-array correlation analysis, confirming and extending the finding that RMA and GCRMA consistently overestimate sample similarity upon normalization. We determine that median polish summarization is responsible for generating a large proportion of these over-similarity artifacts. Furthermore, we show that most affected probesets show also internal signal disagreement, and tend to be composed by individual probes hitting different gene transcripts. We finally provide a correction to the RMA/GCRMA summarization procedure that massively reduces inter-array correlation artifacts, without affecting the detection of differentially expressed genes.ConclusionsWe propose tRMA as a modification of RMA to normalize microarray experiments for correlation-based analysis.