ScalaFlux: A scalable approach to quantify fluxes in metabolic subnetworks.

13C-metabolic flux analysis (13C-MFA) allows metabolic fluxes to be quantified in living organisms and is a major tool in biotechnology and systems biology. Current 13C-MFA approaches model label propagation starting from the extracellular 13C-labeled nutrient(s), which limits their applicability to the analysis of pathways close to this metabolic entry point. Here, we propose a new approach to quantify fluxes through any metabolic subnetwork of interest by modeling label propagation directly from the metabolic precursor(s) of this subnetwork. The flux calculations are thus purely based on information from within the subnetwork of interest, and no additional knowledge about the surrounding network (such as atom transitions in upstream reactions or the labeling of the extracellular nutrient) is required. This approach, termed ScalaFlux for SCALAble metabolic FLUX analysis, can be scaled up from individual reactions to pathways to sets of pathways. ScalaFlux has several benefits compared with current 13C-MFA approaches: greater network coverage, lower data requirements, independence from cell physiology, robustness to gaps in data and network information, better computational efficiency, applicability to rich media, and enhanced flux identifiability. We validated ScalaFlux using a theoretical network and simulated data. We also used the approach to quantify fluxes through the prenyl pyrophosphate pathway of Saccharomyces cerevisiae mutants engineered to produce phytoene, using a dataset for which fluxes could not be calculated using existing approaches. A broad range of metabolic systems can be targeted with minimal cost and effort, making ScalaFlux a valuable tool for the analysis of metabolic fluxes.

Statistical control of peptide and protein error rates in large-scale targeted data-independent acquisition analyses.

Liquid chromatography coupled to tandem mass spectrometry (LC-MS/MS) is the main method for high-throughput identification and quantification of peptides and inferred proteins. Within this field, data-independent acquisition (DIA) combined with peptide-centric scoring, as exemplified by the technique SWATH-MS, has emerged as a scalable method to achieve deep and consistent proteome coverage across large-scale data sets. We demonstrate that statistical concepts developed for discovery proteomics based on spectrum-centric scoring can be adapted to large-scale DIA experiments that have been analyzed with peptide-centric scoring strategies, and we provide guidance on their application. We show that optimal tradeoffs between sensitivity and specificity require careful considerations of the relationship between proteins in the samples and proteins represented in the spectral library. We propose the application of a global analyte constraint to prevent the accumulation of false positives across large-scale data sets. Furthermore, to increase the quality and reproducibility of published proteomic results, well-established confidence criteria should be reported for the detected peptide queries, peptides and inferred proteins.

OpenMS: a flexible open-source software platform for mass spectrometry data analysis.

High-resolution mass spectrometry (MS) has become an important tool in the life sciences, contributing to the diagnosis and understanding of human diseases, elucidating biomolecular structural information and characterizing cellular signaling networks. However, the rapid growth in the volume and complexity of MS data makes transparent, accurate and reproducible analysis difficult. We present OpenMS 2.0 (http://www.openms.de), a robust, open-source, cross-platform software specifically designed for the flexible and reproducible analysis of high-throughput MS data. The extensible OpenMS software implements common mass spectrometric data processing tasks through a well-defined application programming interface in C++ and Python and through standardized open data formats. OpenMS additionally provides a set of 185 tools and ready-made workflows for common mass spectrometric data processing tasks, which enable users to perform complex quantitative mass spectrometric analyses with ease.

DIANA--algorithmic improvements for analysis of data-independent acquisition MS data.

MOTIVATION: Data independent acquisition mass spectrometry has emerged as a reproducible and sensitive alternative in quantitative proteomics, where parsing the highly complex tandem mass spectra requires dedicated algorithms. Recently, targeted data extraction was proposed as a novel analysis strategy for this type of data, but it is important to further develop these concepts to provide quality-controlled, interference-adjusted and sensitive peptide quantification. RESULTS: We here present the algorithm DIANA and the classifier PyProphet, which are based on new probabilistic sub-scores to classify the chromatographic peaks in targeted data-independent acquisition data analysis. The algorithm is capable of providing accurate quantitative values and increased recall at a controlled false discovery rate, in a complex gold standard dataset. Importantly, we further demonstrate increased confidence gained by the use of two complementary data-independent acquisition targeted analysis algorithms, as well as increased numbers of quantified peptide precursors in complex biological samples. AVAILABILITY AND IMPLEMENTATION: DIANA is implemented in scala and python and available as open source (Apache 2.0 license) or pre-compiled binaries from http://quantitativeproteomics.org/diana. PyProphet can be installed from PyPi (https://pypi.python.org/pypi/pyprophet). SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

DynaMet: a fully automated pipeline for dynamic LC-MS data.

Dynamic isotope labeling data provides crucial information about the operation of metabolic pathways and are commonly generated via liquid chromatography-mass spectrometry (LC-MS). Metabolome-wide analysis is challenging as it requires grouping of metabolite features over different samples. We developed DynaMet for fully automated investigations of isotope labeling experiments from LC-high-resolution MS raw data. DynaMet enables untargeted extraction of metabolite labeling profiles and provides integrated tools for expressive data visualization. To validate DynaMet we first used time course labeling data of the model strain Bacillus methanolicus from (13)C methanol resulting in complex spectra in multicarbon compounds. Analysis of two biological replicates revealed high robustness and reproducibility of the pipeline. In total, DynaMet extracted 386 features showing dynamic labeling within 10 min. Of these features, 357 could be fitted by implemented kinetic models. Feature identification against KEGG database resulted in 215 matches covering multiple pathways of core metabolism and major biosynthetic routes. Moreover, we performed time course labeling experiment with Escherichia coli on uniformly labeled (13)C glucose resulting in a comparable number of detected features with labeling profiles of high quality. The distinct labeling patterns of common central metabolites generated from both model bacteria can readily be explained by one versus multicarbon compound metabolism. DynaMet is freely available as an extension package for Python based eMZed2, an open source framework built for rapid development of LC-MS data analysis workflows.

pyOpenMS: a Python-based interface to the OpenMS mass-spectrometry algorithm library.

pyOpenMS is an open-source, Python-based interface to the C++ OpenMS library, providing facile access to a feature-rich, open-source algorithm library for MS-based proteomics analysis. It contains Python bindings that allow raw access to the data structures and algorithms implemented in OpenMS, specifically those for file access (mzXML, mzML, TraML, mzIdentML among others), basic signal processing (smoothing, filtering, de-isotoping, and peak-picking) and complex data analysis (including label-free, SILAC, iTRAQ, and SWATH analysis tools). pyOpenMS thus allows fast prototyping and efficient workflow development in a fully interactive manner (using the interactive Python interpreter) and is also ideally suited for researchers not proficient in C++. In addition, our code to wrap a complex C++ library is completely open-source, allowing other projects to create similar bindings with ease. The pyOpenMS framework is freely available at https://pypi.python.org/pypi/pyopenms while the autowrap tool to create Cython code automatically is available at https://pypi.python.org/pypi/autowrap (both released under the 3-clause BSD licence).

eMZed: an open source framework in Python for rapid and interactive development of LC/MS data analysis workflows.

SUMMARY: The Python-based, open-source eMZed framework was developed for mass spectrometry (MS) users to create tailored workflows for liquid chromatography (LC)/MS data analysis. The goal was to establish a unique framework with comprehensive basic functionalities that are easy to apply and allow for the extension and modification of the framework in a straightforward manner. eMZed supports the iterative development and prototyping of individual evaluation strategies by providing a computing environment and tools for inspecting and modifying underlying LC/MS data. The framework specifically addresses non-expert programmers, as it requires only basic knowledge of Python and relies largely on existing successful open-source software, e.g. OpenMS. AVAILABILITY: The framework eMZed and its documentation are freely available at http://emzed.biol.ethz.ch/. eMZed is published under the GPL 3.0 license, and an online discussion group is available at https://groups.google.com/group/emzed-users. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Three-dimensional imaging of xylem at cell wall level through near field nano holotomography.

Detailed imaging of the three-dimensionally complex architecture of xylary plants is important for studying biological and mechanical functions of woody plants. Apart from common two-dimensional microscopy, X-ray micro-computed tomography has been established as a three-dimensional (3D) imaging method for studying the hydraulic function of wooden plants. However, this X-ray imaging method can barely reach the resolution needed to see the minute structures (e.g. pit membrane). To complement the xylem structure with 3D views at the nanoscale level, X-ray near-field nano-holotomography (NFH) was applied to analyze the wood species Pinus sylvestris and Fagus sylvatica. The demanded small specimens required focused ion beam (FIB) application. The FIB milling, however, influenced the image quality through gallium implantation on the cell-wall surfaces. The measurements indicated that NFH is appropriate for imaging wood at nanometric resolution. With a 26 nm voxel pitch, the structure of the cell-wall surface in Pinus sylvestris could be visualized in genuine detail. In wood of Fagus sylvatica, the structure of a pit pair, including the pit membrane, between two neighboring fibrous cells could be traced tomographically.

Quantitative Proteome Landscape of the NCI-60 Cancer Cell Lines.

Here we describe a proteomic data resource for the NCI-60 cell lines generated by pressure cycling technology and SWATH mass spectrometry. We developed the DIA-expert software to curate and visualize the SWATH data, leading to reproducible detection of over 3,100 SwissProt proteotypic proteins and systematic quantification of pathway activities. Stoichiometric relationships of interacting proteins for DNA replication, repair, the chromatin remodeling NuRD complex, ß-catenin, RNA metabolism, and prefoldins are more evident than that at the mRNA level. The data are available in CellMiner (discover.nci.nih.gov/cellminercdb and discover.nci.nih.gov/cellminer), allowing casual users to test hypotheses and perform integrative, cross-database analyses of multi-omic drug response correlations for over 20,000 drugs. We demonstrate the value of proteome data in predicting drug response for over 240 clinically relevant chemotherapeutic and targeted therapies. In summary, we present a novel proteome resource for the NCI-60, together with relevant software tools, and demonstrate the benefit of proteome analyses.

Estimating the onset of cambial activity in Scots pine in northern Finland by means of the heat-sum approach.

We estimated the date of onset (Date(est)) of cambial activity by the pinning method in Scots pine (Pinus sylvestris L.) trees at Vanttauskoski (Site 1) and Laanila (Site 2) near the latitudinal limit of Scots pine in northern Finland. In each year and at each site, observations were made on a different set of five trees. The estimated dates of onset of cambial activity were compared with the corresponding heat sums, calculated in degree-days according to two models. Within years, Date(est) varied among trees by up to 15 days at Site 1 and up to 13 days at Site 2. Among years, mean Date(est) varied by 15.3 days at Site 1 and 12.0 days at Site 2. The overall mean Date(est) differed between sites by 6 days (June 5 at Site 1 and June 11 at Site 2). Among all trees in all years, the mean number of degree days (d.d.) calculated from mean daily temperature above a threshold of 5 degrees C before Date(est) ranged from 68.7 to 135 d.d. at Site 1 and from 37.4 to 154.7 d.d. at Site 2. Among years, the mean heat sum before Date(est )ranged from 94 to 112.5 d.d. at Site 1 and from 61.4 to 136 d.d. at Site 2. Variation among years in heat sum before Date(est) at Site 2 was highly significant, indicating that one or more factors other than, or in addition to, heat sum determines the onset of cambial activity in Scots pine. Similar results were obtained when heat sum was computed from the area between the sine wave generated by daily maximum and minimum temperature and the threshold temperature.

Lignification in cell cultures of Pinus radiata: activities of enzymes and lignin topochemistry.

Enzymatic and topochemical aspects of lignification were studied in a Pinus radiata D. Don cell culture system that was induced to differentiate tracheary elements and sclereids with lignified secondary cell walls. The activities of the lignin-related enzymes phenylalanine ammonia lyase (PAL; EC 4.3.1.5) and cinnamyl alcohol dehydrogenase (CAD; EC 1.1.1.195) increased concomitantly with cell differentiation, indicating that the increase in enzyme activity was related to lignification of the cell walls and was not induced by stress. This result also indicates that PAL and CAD are suitable markers for tracheary element differentiation in coniferous gymnosperms. To further characterize lignification in this cell culture system, cellular UV-microspectrophotometry and thioacidolysis were employed. Typical UV-absorption spectra of lignin were obtained from the secondary cell walls of the tracheary elements and sclereids and from the compound middle lamella connecting differentiated cells, and the presence of lignin was confirmed by thioacidolysis. Certain aspects of lignin topochemistry in the cell walls of the tracheary elements were similar to cell walls of P. radiata wood, such as the high lignin concentration in the compound middle lamella connecting adjacent cells and the lower lignin concentration in the secondary cell walls. Therefore, the P. radiata cell culture system appears to be well suited to study the formation of lignified secondary cell walls in coniferous gymnosperms.

Annual Cambial Rhythm in Pinus halepensis and Pinus sylvestris as Indicator for Climate Adaptation.

To understand better the adaptation strategies of intra-annual radial growth in Pinus halepensis and Pinus sylvestris to local environmental conditions, we examined the seasonal rhythm of cambial activity and cell differentiation at tissue and cellular levels. Two contrasting sites differing in temperature and amount of precipitation were selected for each species, one typical for their growth and the other represented border climatic conditions, where the two species coexisted. Mature P. halepensis trees from Mediterranean (Spain) and sub-Mediterranean (Slovenia) sites, and P. sylvestris from sub-Mediterranean (Slovenia) and temperate (Slovenia) sites were selected. Repeated sampling was performed throughout the year and samples were prepared for examination with light and transmission electron microscopes. We hypothesized that cambial rhythm in trees growing at the sub-Mediterranean site where the two species co-exist will be similar as at typical sites for their growth. Cambium in P. halepensis at the Mediterranean site was active throughout the year and was never truly dormant, whereas at the sub-Mediterranean site it appeared to be dormant during the winter months. In contrast, cambium in P. sylvestris was clearly dormant at both sub-Mediterranean and temperate sites, although the dormant period seemed to be significantly longer at the temperate site. Thus, the hypothesis was only partly confirmed. Different cambial and cell differentiation rhythms of the two species at the site where both species co-exist and typical sites for their growth indicate their high but different adaptation strategies in terms of adjustment of radial growth to environmental heterogeneity, crucial for long-term tree performance and survival.

Fast and Efficient XML Data Access for Next-Generation Mass Spectrometry.

MOTIVATION: In mass spectrometry-based proteomics, XML formats such as mzML and mzXML provide an open and standardized way to store and exchange the raw data (spectra and chromatograms) of mass spectrometric experiments. These file formats are being used by a multitude of open-source and cross-platform tools which allow the proteomics community to access algorithms in a vendor-independent fashion and perform transparent and reproducible data analysis. Recent improvements in mass spectrometry instrumentation have increased the data size produced in a single LC-MS/MS measurement and put substantial strain on open-source tools, particularly those that are not equipped to deal with XML data files that reach dozens of gigabytes in size. RESULTS: Here we present a fast and versatile parsing library for mass spectrometric XML formats available in C++ and Python, based on the mature OpenMS software framework. Our library implements an API for obtaining spectra and chromatograms under memory constraints using random access or sequential access functions, allowing users to process datasets that are much larger than system memory. For fast access to the raw data structures, small XML files can also be completely loaded into memory. In addition, we have improved the parsing speed of the core mzML module by over 4-fold (compared to OpenMS 1.11), making our library suitable for a wide variety of algorithms that need fast access to dozens of gigabytes of raw mass spectrometric data. AVAILABILITY: Our C++ and Python implementations are available for the Linux, Mac, and Windows operating systems. All proposed modifications to the OpenMS code have been merged into the OpenMS mainline codebase and are available to the community at https://github.com/OpenMS/OpenMS.

Identification of a hydroxy substituted calamenene--a sesquiterpene associated with wound reactions in non-infected xylem of Tilia spp.

Xylem of lime trees (Tilia spp.) with wound reactions was structurally investigated by scanning (SEM) and transmission electron microscopy (TEM) as well as chemically analyzed by direct thermal desorption-gas chromatography-mass spectrometry (DTD-GC-MS). Wound reactions in the outer xylem lead to distinct discolorations around the wound. Within a 4-week response no fungal infection occurred in discoloured xylem. At the fine structural level, wound reactions become primarily visible as the secretion of dark-staining substances from parenchyma cells into lumens of vessels and fibres. With increasing reaction time vessels aggregate large amounts of secretion products, whereas in fibres wall-associated linings are formed and the inner secondary wall appears incrusted. After 2-3 months a narrow, greenish-brown boundary developed at the transition between the discoloured outer and the unchanged inner xylem. This green-brown boundary layer remained non-infected also in older wounds. DTD-GC-MS analyses revealed that the sesquiterpene Hydroxycalamenene represents a key substance of wound reactions in non-infected lime trees. Other substances such as fatty acids or their esters and coniferyl aldehydes or their derivatives were also found. TEM investigations of the samples after DTD-GC-MS showed less pronounced cell wall-attached linings in fibres as well as reduced incrustation of inner secondary walls. The massive deposits in the vessel lumens remained unchanged. The role of these wound reaction products and their ways of synthesis are discussed.

Topochemical and transmission electron microscopic studies of bacterial decay in pine (Pinus sylvestris L.) harbour foundation piles.

Bacterial degradation of Pinus sylvestris harbour foundation piles was studied topochemically by scanning UV-microspectrophotometry. This analytical technique enables direct imaging of lignin distribution within individual cell wall layers. Additionally, light and transmission electron microscopy (TEM) was used to characterise structural changes of the cell walls. Various decay stages were found in the samples. TEM revealed that the bacterial degradation occurred mainly in the S2, leaving granular remnants in degraded wall portions with lower as well as higher electron density than the surrounding unmodified wall. In the initial stages, topochemical investigations revealed that lignin modification starts in the innermost parts of the secondary wall, most clearly observed in latewood tracheids. During advanced degradation, lignin modification occurs more or less severe in walls of all cell types. However, even in cell portions with intensive decay, the compound middle lamellae and ray tracheids were undegraded. The knowledge about lignin modification at initial stages of wood degradation by bacteria is of fundamental importance to provide more information on the process of cell wall decay.

Suppression of 4-coumarate-CoA ligase in the coniferous gymnosperm Pinus radiata.

Severe suppression of 4-coumarate-coenzyme A ligase (4CL) in the coniferous gymnosperm Pinus radiata substantially affected plant phenotype and resulted in dwarfed plants with a "bonsai tree-like" appearance. Microscopic analyses of stem sections from 2-year-old plants revealed substantial morphological changes in both wood and bark tissues. This included the formation of weakly lignified tracheids that displayed signs of collapse and the development of circumferential bands of axial parenchyma. Acetyl bromide-soluble lignin assays and proton nuclear magnetic resonance studies revealed lignin reductions of 36% to 50% in the most severely affected transgenic plants. Two-dimensional nuclear magnetic resonance and pyrolysis-gas chromatography-mass spectrometry studies indicated that lignin reductions were mainly due to depletion of guaiacyl but not p-hydroxyphenyl lignin. 4CL silencing also caused modifications in the lignin interunit linkage distribution, including elevated beta-aryl ether (beta-O-4 unit) and spirodienone (beta-1) levels, accompanied by lower phenylcoumaran (beta-5), resinol (beta-beta), and dibenzodioxocin (5-5/beta-O-4) levels. A sharp depletion in the level of saturated (dihydroconiferyl alcohol) end groups was also observed. Severe suppression of 4CL also affected carbohydrate metabolism. Most obvious was an up to approximately 2-fold increase in galactose content in wood from transgenic plants due to increased compression wood formation. The molecular, anatomical, and analytical data verified that the isolated 4CL clone is associated with lignin biosynthesis and illustrated that 4CL silencing leads to complex, often surprising, physiological and morphological changes in P. radiata.

Intervessel pit structure and histochemistry of two mangrove species as revealed by cellular UV microspectrophotometry and electron microscopy: intraspecific variation and functional significance.

Intervessel pits play a key role in trees' water transport, lying at the base of drought-induced embolism, and in the regulation of hydraulic conductivity via hydrogels bordering pit canals. Recently, their microstructure has been the focus of numerous studies, but the considerable variation, even within species and the histochemistry of pit membranes, remains largely unexplained. In the present study, intervessel pits of the outermost wood were examined for Avicennia marina, of dry and rainy season wood separately for Rhizophora mucronata. The thickness of the pit membranes was measured on transmission electron micrographs while their topochemical nature was also analyzed via cellular UV microspectrophotometry. Pit membranes of R. mucronata were slightly thicker in dry season wood than in rainy season wood, but their spectra showed for both seasons a lignin and a yet unidentified higher wavelength absorbing component. It was suggested to be a derivative of the deposits, regularly filling pit canals. The vestures of A. marina chemically resembled pit membranes rather than cell walls.

Effect of local heating and cooling on cambial activity and cell differentiation in the stem of Norway spruce (Picea abies).

UNLABELLED: BACKGROUND AND AIMS The effect of heating and cooling on cambial activity and cell differentiation in part of the stem of Norway spruce (Picea abies) was investigated. METHODS: A heating experiment (23-25 degrees C) was carried out in spring, before normal reactivation of the cambium, and cooling (9-11 degrees C) at the height of cambial activity in summer. The cambium, xylem and phloem were investigated by means of light- and transmission electron microscopy and UV-microspectrophotometry in tissues sampled from living trees. KEY RESULTS: Localized heating for 10 d initiated cambial divisions on the phloem side and after 20 d also on the xylem side. In a control tree, regular cambial activity started after 30 d. In the heat-treated sample, up to 15 earlywood cells undergoing differentiation were found to be present. The response of the cambium to stem cooling was less pronounced, and no anatomical differences were detected between the control and cool-treated samples after 10 or 20 d. After 30 d, latewood started to form in the sample exposed to cooling. In addition, almost no radially expanding tracheids were observed and the cambium consisted of only five layers of cells. Low temperatures reduced cambial activity, as indicated by the decreased proportion of latewood. On the phloem side, no alterations were observed among cool-treated and non-treated samples. CONCLUSIONS: Heating and cooling can influence cambial activity and cell differentiation in Norway spruce. However, at the ultrastructural and topochemical levels, no changes were observed in the pattern of secondary cell-wall formation and lignification or in lignin structure, respectively.

Differentiation of terminal latewood tracheids in silver fir trees during autumn.

BACKGROUND AND AIMS: The differentiation of terminal latewood tracheids of silver fir (Abies alba) trees grown in Slovenia was investigated in autumn/winter 2001/2002. METHODS: The experimental trees were divided into three groups: one with narrow annual rings, width less than 1 mm; one with annual ring widths between 1 and 4 mm; and one group with broad rings larger than 4 mm. The differentiation of terminal latewood tracheids was investigated by light-, electron- and UV-microscopy in tissues sampled in October and November 2001 and March 2002. KEY RESULTS: In the middle of October, cambial divisions did not occur any more in any of the trees. In trees with narrow annual rings, cell wall deposition as well as lignification were completed in terminal latewood tracheids at this date, whereas in trees with annual ring widths of more than 1 mm these processes still continued. Electron microscopy as well as UV microscopy revealed an unlignified inner S(2) layer and the absence of S(3) and warty layers. With increasing distance from the cambium, wall formation and lignification gradually appeared to be completed. Samples of all trees taken in the middle of November only contained differentiated terminal latewood tracheids. At the structural and lignin topochemical level, November and March samples showed completed differentiation of walls of terminal latewood tracheids. CONCLUSIONS: In trees with broader annual rings, the final steps of differentiation of the youngest latewood tracheids near the cambium still continued during autumn, but were finished prior to winter. It was concluded from structural observations that duration of cambial activity is longer in trees with broad annual rings than in trees with narrow rings.