Hierarchical time-series analysis of dynamic bioprocess systems.

BACKGROUND: Monoclonal antibodies (mAbs) are leading types of 'blockbuster' biotherapeutics worldwide; they have been successfully used to treat various cancers and chronic inflammatory and autoimmune diseases. Biotherapeutics process development and manufacturing are complicated due to lack of understanding the factors that impact cell productivity and product quality attributes. Understanding complex interactions between cells, media, and process parameters on the molecular level is essential to bring biomanufacturing to the next level. This can be achieved by analyzing cell culture metabolic levels connected to vital process parameters like viable cell density (VCD). However, VCD and metabolic profiles are dynamic parameters and inherently correlated with time, leading to a significant correlation without actual causality. Many time-series methods deal with such issues. However, with metabolic profiling, the number of measured variables vastly exceeds the number of experiments, making most of existing methods ill-suited and hard to interpret. METHODS AND MAJOR RESULTS: Here we propose an alternative workflow using hierarchical dimension reduction to visualize and interpret the relation between evolution of metabolic profiles and dynamic process parameters. The first step of proposed method is focused on finding predictive relation between metabolic profiles and process parameter at all time points using OPLS regression. For each time point, the p(corr) obtained from OPLS model is considered as a differential metabogram and is further assessed using principal components analysis (PCA). CONCLUSIONS: Compared to traditional batch modeling, applying proposed methodology on metabolic data from Chinese Hamster Ovary (CHO) antibody production characterized the dynamic relation between metabolic profiles and critical process parameters.

A whole-genome sequenced control population in northern Sweden reveals subregional genetic differences.

The number of national reference populations that are whole-genome sequenced are rapidly increasing. Partly driving this development is the fact that genetic disease studies benefit from knowing the genetic variation typical for the geographical area of interest. A whole-genome sequenced Swedish national reference population (n = 1000) has been recently published but with few samples from northern Sweden. In the present study we have whole-genome sequenced a control population (n = 300) (ACpop) from Västerbotten County, a sparsely populated region in northern Sweden previously shown to be genetically different from southern Sweden. The aggregated variant frequencies within ACpop are publicly available (DOI 10.17044/NBIS/G000005) to function as a basic resource in clinical genetics and for genetic studies. Our analysis of ACpop, representing approximately 0.11% of the population in Västerbotten, indicates the presence of a genetic substructure within the county. Furthermore, a demographic analysis showed that the population from which samples were drawn was to a large extent geographically stationary, a finding that was corroborated in the genetic analysis down to the level of municipalities. Including ACpop in the reference population when imputing unknown variants in a Västerbotten cohort resulted in a strong increase in the number of high-confidence imputed variants (up to 81% for variants with minor allele frequency < 5%). ACpop was initially designed for cancer disease studies, but the genetic structure within the cohort will be of general interest for all genetic disease studies in northern Sweden.

Overexpression of VviPGIP1 and NtCAD14 in Tobacco Screened Using Glycan Microarrays Reveals Cell Wall Reorganisation in the Absence of Fungal Infection.

The expression of Vitis vinifera polygalacturonase inhibiting protein 1 (VviPGIP1) in Nicotiana tabacum has been linked to modifications at the cell wall level. Previous investigations have shown an upregulation of the lignin biosynthesis pathway and reorganisation of arabinoxyloglucan composition. This suggests cell wall tightening occurs, which may be linked to defence priming responses. The present study used a screening approach to test four VviPGIP1 and four NtCAD14 overexpressing transgenic lines for cell wall alterations. Overexpressing the tobacco-derived cinnamyl alcohol dehydrogenase (NtCAD14) gene is known to increase lignin biosynthesis and deposition. These lines, particularly PGIP1 expressing plants, have been shown to lead to a decrease in susceptibility towards grey rot fungus Botrytis cinerea. In this study the aim was to investigate the cell wall modulations that occurred prior to infection, which should highlight potential priming phenomena and phenotypes. Leaf lignin composition and relative concentration of constituent monolignols were evaluated using pyrolysis gas chromatography. Significant concentrations of lignin were deposited in the stems but not the leaves of NtCAD14 overexpressing plants. Furthermore, no significant changes in monolignol composition were found between transgenic and wild type plants. The polysaccharide modifications were quantified using gas chromatography (GC-MS) of constituent monosaccharides. The major leaf polysaccharide and cell wall protein components were evaluated using comprehensive microarray polymer profiling (CoMPP). The most significant changes appeared at the polysaccharide and protein level. The pectin fraction of the transgenic lines had subtle variations in patterning for methylesterification epitopes for both VviPGIP1 and NtCAD14 transgenic lines versus wild type. Pectin esterification levels have been linked to pathogen defence in the past. The most marked changes occurred in glycoprotein abundance for both the VviPGIP1 and NtCAD14 lines. Epitopes for arabinogalactan proteins (AGPs) and extensins were notably altered in transgenic NtCAD14 tobacco.

Analysis of Plant Cell Walls Using High-Throughput Profiling Techniques with Multivariate Methods.

Plant cell walls are composed of a number of coextensive polysaccharide-rich networks (i.e., pectin, hemicellulose, protein). Polysaccharide-rich cell walls are important in a number of biological processes including fruit ripening, plant-pathogen interactions (e.g., pathogenic fungi), fermentations (e.g., winemaking), and tissue differentiation (e.g., secondary cell walls). Applying appropriate methods is necessary to assess biological roles as for example in putative plant gene functional characterization (e.g., experimental evaluation of transgenic plants). Obtaining datasets is relatively easy, using for example gas chromatography-mass spectrometry (GC-MS) methods for monosaccharide composition, Fourier transform infrared spectroscopy (FT-IR) and comprehensive microarray polymer profiling (CoMPP); however, analyzing the data requires implementing statistical tools for large-scale datasets. We have validated and implemented a range of multivariate data analysis methods on datasets from tobacco, grapevine, and wine polysaccharide studies. Here we present the workflow from processing samples to acquiring data to performing data analysis (particularly principal component analysis (PCA) and orthogonal projection to latent structure (OPLS) methods).

A geographically matched control population efficiently limits the number of candidate disease-causing variants in an unbiased whole-genome analysis.

Whole-genome sequencing is a promising approach for human autosomal dominant disease studies. However, the vast number of genetic variants observed by this method constitutes a challenge when trying to identify the causal variants. This is often handled by restricting disease studies to the most damaging variants, e.g. those found in coding regions, and overlooking the remaining genetic variation. Such a biased approach explains in part why the genetic causes of many families with dominantly inherited diseases, in spite of being included in whole-genome sequencing studies, are left unsolved today. Here we explore the use of a geographically matched control population to minimize the number of candidate disease-causing variants without excluding variants based on assumptions on genomic position or functional predictions. To exemplify the benefit of the geographically matched control population we apply a typical disease variant filtering strategy in a family with an autosomal dominant form of colorectal cancer. With the use of the geographically matched control population we end up with 26 candidate variants genome wide. This is in contrast to the tens of thousands of candidates left when only making use of available public variant datasets. The effect of the local control population is dual, it (1) reduces the total number of candidate variants shared between affected individuals, and more importantly (2) increases the rate by which the number of candidate variants are reduced as additional affected family members are included in the filtering strategy. We demonstrate that the application of a geographically matched control population effectively limits the number of candidate disease-causing variants and may provide the means by which variants suitable for functional studies are identified genome wide.

Cytokine Profiles in Autoantibody Defined Subgroups of Systemic Lupus Erythematosus.

The aim of this study was to evaluate how the cytokine profiles differed between autoantibody based subgroups of systemic lupus erythematosus (SLE). SLE is a systemic autoimmune disease, characterized by periods of flares (active disease) and remission (inactive disease). The disease can affect many organ systems, e.g., skin, joints, kidneys, heart, and the central nervous system (CNS). SLE patients often have an overproduction of cytokines, e.g., interferons, chemokines, and interleukins. The high cytokine levels are part of the systemic inflammation, which can lead to tissue injury. In the present study, SLE patients were divided into five groups based on their autoantibody profiles. We thus defined these five groups: ANA negative, antiphospholipid (aPL) positive, anti-Sm/anti-RNP positive, Sjögren's syndrome (SS) antigen A and B positive, and patients positive for more than one type of autoantibodies (other SLE). Cytokines were measured using Mesoscale Discovery (MSD) multiplex analysis. On the basis of the cytokine data, ANA negative patients were the most deviating subgroup, with lower levels of interferon (IFN)-γ, tumor necrosis factor (TNF)-α, interleukin (IL)-12/IL-23p40, and interferon gamma-induced protein (IP)-10. Despite low cytokine levels in the ANA negative group, autoantibody profiles did not discriminate between different cytokine patterns.

Association between pathology and texture features of multi parametric MRI of the prostate.

The role of multi-parametric (mp)MRI in the diagnosis and treatment of prostate cancer has increased considerably. An alternative to visual inspection of mpMRI is the evaluation using histogram-based (first order statistics) parameters and textural features (second order statistics). The aims of the present work were to investigate the relationship between benign and malignant sub-volumes of the prostate and textures obtained from mpMR images. The performance of tumor prediction was investigated based on the combination of histogram-based and textural parameters. Subsequently, the relative importance of mpMR images was assessed and the benefit of additional imaging analyzed. Finally, sub-structures based on the PI-RADS classification were investigated as potential regions to automatically detect maligned lesions. Twenty-five patients who received mpMRI prior to radical prostatectomy were included in the study. The imaging protocol included T2, DWI, and DCE. Delineation of tumor regions was performed based on pathological information. First and second order statistics were derived from each structure and for all image modalities. The resulting data were processed with multivariate analysis, using PCA (principal component analysis) and OPLS-DA (orthogonal partial least squares discriminant analysis) for separation of malignant and healthy tissue. PCA showed a clear difference between tumor and healthy regions in the peripheral zone for all investigated images. The predictive ability of the OPLS-DA models increased for all image modalities when first and second order statistics were combined. The predictive value reached a plateau after adding ADC and T2, and did not increase further with the addition of other image information. The present study indicates a distinct difference in the signatures between malign and benign prostate tissue. This is an absolute prerequisite for automatic tumor segmentation, but only the first step in that direction. For the specific identified signature, DCE did not add complementary information to T2 and ADC maps.

Haralick texture features from apparent diffusion coefficient (ADC) MRI images depend on imaging and pre-processing parameters.

In recent years, texture analysis of medical images has become increasingly popular in studies investigating diagnosis, classification and treatment response assessment of cancerous disease. Despite numerous applications in oncology and medical imaging in general, there is no consensus regarding texture analysis workflow, or reporting of parameter settings crucial for replication of results. The aim of this study was to assess how sensitive Haralick texture features of apparent diffusion coefficient (ADC) MR images are to changes in five parameters related to image acquisition and pre-processing: noise, resolution, how the ADC map is constructed, the choice of quantization method, and the number of gray levels in the quantized image. We found that noise, resolution, choice of quantization method and the number of gray levels in the quantized images had a significant influence on most texture features, and that the effect size varied between different features. Different methods for constructing the ADC maps did not have an impact on any texture feature. Based on our results, we recommend using images with similar resolutions and noise levels, using one quantization method, and the same number of gray levels in all quantized images, to make meaningful comparisons of texture feature results between different subjects.

Heterozygous colon cancer-associated mutations of SAMHD1 have functional significance.

Even small variations in dNTP concentrations decrease DNA replication fidelity, and this observation prompted us to analyze genomic cancer data for mutations in enzymes involved in dNTP metabolism. We found that sterile alpha motif and histidine-aspartate domain-containing protein 1 (SAMHD1), a deoxyribonucleoside triphosphate triphosphohydrolase that decreases dNTP pools, is frequently mutated in colon cancers, that these mutations negatively affect SAMHD1 activity, and that several SAMHD1 mutations are found in tumors with defective mismatch repair. We show that minor changes in dNTP pools in combination with inactivated mismatch repair dramatically increase mutation rates. Determination of dNTP pools in mouse embryos revealed that inactivation of one SAMHD1 allele is sufficient to elevate dNTP pools. These observations suggest that heterozygous cancer-associated SAMHD1 mutations increase mutation rates in cancer cells.

Quantitative proteomics reveals protein profiles underlying major transitions in aspen wood development.

BACKGROUND: Wood development is of outstanding interest both to basic research and industry due to the associated cellulose and lignin biomass production. Efforts to elucidate wood formation (which is essential for numerous aspects of both pure and applied plant science) have been made using transcriptomic analyses and/or low-resolution sampling. However, transcriptomic data do not correlate perfectly with levels of expressed proteins due to effects of post-translational modifications and variations in turnover rates. In addition, high-resolution analysis is needed to characterize key transitions. In order to identify protein profiles across the developmental region of wood formation, an in-depth and tissue specific sampling was performed. RESULTS: We examined protein profiles, using an ultra-performance liquid chromatography/quadrupole time of flight mass spectrometry system, in high-resolution tangential sections spanning all wood development zones in Populus tremula from undifferentiated cambium to mature phloem and xylem, including cell expansion and cell death zones. In total, we analyzed 482 sections, 20-160 µm thick, from four 47-year-old trees growing wild in Sweden. We obtained high quality expression profiles for 3,082 proteins exhibiting consistency across the replicates, considering that the trees were growing in an uncontrolled environment. A combination of Principal Component Analysis (PCA), Orthogonal Projections to Latent Structures (OPLS) modeling and an enhanced stepwise linear modeling approach identified several major transitions in global protein expression profiles, pinpointing (for example) locations of the cambial division leading to phloem and xylem cells, and secondary cell wall formation zones. We also identified key proteins and associated pathways underlying these developmental landmarks. For example, many of the lignocellulosic related proteins were upregulated in the expansion to the early developmental xylem zone, and for laccases with a rapid decrease in early xylem zones. We observed upregulation of two forms of xylem cysteine protease (Potri.002G005700.1 and Potri.005G256000.2; Pt-XCP2.1) in early xylem and their downregulation in late maturing xylem. Our data also show that Pt-KOR1.3 (Potri.003G151700.2) exhibits an expression pattern that supports the hypothesis put forward in previous studies that this is a key xyloglucanase involved in cellulose biosynthesis in primary cell walls and reduction of cellulose crystallinity in secondary walls. CONCLUSION: Our novel multivariate approach highlights important processes and provides confirmatory insights into the molecular foundations of wood development.

Metabolomics study of fatigue in patients with rheumatoid arthritis naïve to biological treatment.

Fatigue occurs in all chronic inflammatory diseases, in cancer, and in some neurological conditions. Patients often regard fatigue as one of their most debilitating problems, but currently there is no established treatment and the mechanisms that lead to and regulate fatigue are incompletely understood. Our objective was to more completely understand the physiology of this phenomenon. Twenty-four patients with rheumatoid arthritis (RA) naïve to treatment with biological drugs were enrolled for the study. Fatigue was measured with a fatigue visual analogue scale (fVAS). Ethylenediaminetetraacetic acid (EDTA) plasma samples were subjected to gas chromatography-time-of-flight mass spectrometry (GC/MS-TOF)-based metabolite profiling. Obtained metabolite data were evaluated by multivariate data analysis with orthogonal projections to latent structures (OPLS) method to pinpoint metabolic changes related to fatigue severity. A significant multivariate OPLS model was obtained between the fVAS scores and the measured metabolic levels. Increasing fatigue scores were associated with a metabolic pattern characterized by down-regulation of metabolites from the urea cycle, fatty acids, tocopherols, aromatic amino acids, and hypoxanthine. Uric acid levels were increased. Apart from fatigue, we found no other disease-related variables that might be responsible for these changes. Our MS-based metabolomic approach demonstrated strong associations between fatigue and several biochemical patterns related to oxidative stress.

ADC texture--an imaging biomarker for high-grade glioma?

PURPOSE: Survival for high-grade gliomas is poor, at least partly explained by intratumoral heterogeneity contributing to treatment resistance. Radiological evaluation of treatment response is in most cases limited to assessment of tumor size months after the initiation of therapy. Diffusion-weighted magnetic resonance imaging (MRI) and its estimate of the apparent diffusion coefficient (ADC) has been widely investigated, as it reflects tumor cellularity and proliferation. The aim of this study was to investigate texture analysis of ADC images in conjunction with multivariate image analysis as a means for identification of pretreatment imaging biomarkers. METHODS: Twenty-three consecutive high-grade glioma patients were treated with radiotherapy (2 Gy/60 Gy) with concomitant and adjuvant temozolomide. ADC maps and T1-weighted anatomical images with and without contrast enhancement were collected prior to treatment, and (residual) tumor contrast enhancement was delineated. A gray-level co-occurrence matrix analysis was performed on the ADC maps in a cuboid encapsulating the tumor in coronal, sagittal, and transversal planes, giving a total of 60 textural descriptors for each tumor. In addition, similar examinations and analyses were performed at day 1, week 2, and week 6 into treatment. Principal component analysis (PCA) was applied to reduce dimensionality of the data, and the five largest components (scores) were used in subsequent analyses. MRI assessment three months after completion of radiochemotherapy was used for classifying tumor progression or regression. RESULTS: The score scatter plots revealed that the first, third, and fifth components of the pretreatment examinations exhibited a pattern that strongly correlated to survival. Two groups could be identified: one with a median survival after diagnosis of 1099 days and one with 345 days, p = 0.0001. CONCLUSIONS: By combining PCA and texture analysis, ADC texture characteristics were identified, which seems to hold pretreatment prognostic information, independent of known prognostic factors such as age, stage, and surgical procedure. These findings encourage further studies with a larger patient cohort.

Lipid mediator profiles differ between lung compartments in asthmatic and healthy humans.

Oxylipins are oxidised fatty acids that can exert lipid mediator functions in inflammation, and several oxylipins derived from arachidonic acid are linked to asthma. This study quantified oxylipin profiles in different regions of the lung to obtain a broad-scale characterisation of the allergic asthmatic inflammation in relation to healthy individuals. Bronchoalveolar lavage fluid (BALF), bronchial wash fluid and endobronchial mucosal biopsies were collected from 16 healthy and 16 mildly allergic asthmatic individuals. Inflammatory cell counts, immunohistochemical staining and oxylipin profiling were performed. Univariate and multivariate statistics were employed to evaluate compartment-dependent and diagnosis-dependent oxylipin profiles in relation to other measured parameters. Multivariate modelling showed significantly different bronchial wash fluid and BALF oxylipin profiles in both groups (R(2)Y[cum]=0.822 and Q(2)[cum]=0.759). Total oxylipin concentrations and five individual oxylipins, primarily from the lipoxygenase (LOX) pathway of arachidonic and linoleic acid, were elevated in bronchial wash fluid from asthmatics compared to that from healthy controls, supported by immunohistochemical staining of 15-LOX-1 in the bronchial epithelium. No difference between the groups was found among BALF oxylipins. In conclusion, bronchial wash fluid and BALF contain distinct oxylipin profiles, which may have ramifications for the study of respiratory diseases. Specific protocols for sampling proximal and distal airways separately should be employed for lipid mediator studies.

EGFR gene variants are associated with specific somatic aberrations in glioma.

A number of gene variants have been associated with an increased risk of developing glioma. We hypothesized that the reported risk variants may be associated with tumor genomic instability. To explore potential correlations between germline risk variants and somatic genetic events, we analyzed matched tumor and blood samples from 95 glioma patients by means of SNP genotyping. The generated genotype data was used to calculate genome-wide allele-specific copy number profiles of the tumor samples. We compared the copy number profiles across samples and found two EGFR gene variants (rs17172430 and rs11979158) that were associated with homozygous deletion at the CDKN2A/B locus. One of the EGFR variants (rs17172430) was also associated with loss of heterozygosity at the EGFR locus. Our findings were confirmed in a separate dataset consisting of matched blood and tumor samples from 300 glioblastoma patients, compiled from publically available TCGA data. These results imply there is a functional effect of germline EGFR variants on tumor progression.

Dynamic ultrasound imaging--a multivariate approach for the analysis and comparison of time-dependent musculoskeletal movements.

BACKGROUND: Muscle functions are generally assumed to affect a wide variety of conditions and activities, including pain, ischemic and neurological disorders, exercise and injury. It is therefore very desirable to obtain more information on musculoskeletal contributions to and activity during clinical processes such as the treatment of muscle injuries, post-surgery evaluations, and the monitoring of progressive degeneration in neuromuscular disorders.The spatial image resolution achievable with ultrasound systems has improved tremendously in the last few years and it is nowadays possible to study skeletal muscles in real-time during activity. However, ultrasound imaging has an inherent problem that makes it difficult to compare different measurement series or image sequences from two or more subjects. Due to physiological differences between different subjects, the ultrasound sequences will be visually different - partly because of variation in probe placement and partly because of the difficulty of perfectly reproducing any given movement. METHODS: Ultrasound images of the biceps and calf of a single subject were transformed to achieve congruence and then efficiently compressed and stacked to facilitate analysis using a multivariate method known as O2PLS. O2PLS identifies related and unrelated variation in and between two sets of data such that different phases of the studied movements can be analysed. The methodology was used to study the dynamics of the Achilles tendon and the calf and also the Biceps brachii and upper arm. The movements of these parts of the body are both of interest in clinical orthopaedic research. RESULTS: This study extends the novel method of multivariate analysis of congruent images (MACI) to facilitate comparisons between two series of ultrasound images. This increases its potential range of medical applications and its utility for detecting, visualising and quantifying the dynamics and functions of skeletal muscle. CONCLUSIONS: The most important results of this study are that MACI with O2PLS is able to consistently extract meaningful variability from pairs of ultrasound sequences. The MACI method with O2PLS is a powerful tool with great potential for visualising and comparing dynamics between movements. It has many potential clinical applications in the study of muscle injuries, post-surgery evaluations and evaluations of rehabilitation, and the assessment of athletic training interventions.

Metabolic responses to change in disease activity during tumor necrosis factor inhibition in patients with rheumatoid arthritis.

Assessment of disease activity in patients with rheumatoid arthritis (RA) is of importance in the evaluation of treatment. The most important measure of disease activity is the Disease Activity Score counted in 28 joints (DAS28). In this study, we evaluated whether metabolic profiling could complement current measures of disease activity. Fifty-six patients, in two separate studies, were followed for two years after commencing anti-TNF therapy. DAS28 was assessed, and metabolic profiles were recorded at defined time points. Correlations between metabolic profile and DAS28 scores were analyzed using multivariate statistics. The metabolic responses to lowering DAS28 scores varied in different patients but could predict DAS28 scores at the individual and subgroup level models. The erythrocyte sedimentation rate (ESR) component in DAS28 was most correlated to the metabolite data, pointing to inflammation as the primary effect driving metabolic profile changes. Patients with RA had differing metabolic response to changes in DAS28 following anti-TNF therapy. This suggests that discovery of new metabolic biomarkers for disease activity will derive from studies at the individual and subgroup level. Increased inflammation, measured as ESR, was the main common effect seen in metabolic profiles from periods associated with high DAS28.

Modelling human musculoskeletal functional movements using ultrasound imaging.

BACKGROUND: A widespread and fundamental assumption in the health sciences is that muscle functions are related to a wide variety of conditions, for example pain, ischemic and neurological disorder, exercise and injury. It is therefore highly desirable to study musculoskeletal contributions in clinical applications such as the treatment of muscle injuries, post-surgery evaluations, monitoring of progressive degeneration in neuromuscular disorders, and so on.The spatial image resolution in ultrasound systems has improved tremendously in the last few years and nowadays provides detailed information about tissue characteristics. It is now possible to study skeletal muscles in real-time during activity. METHODS: The ultrasound images are transformed to be congruent and are effectively compressed and stacked in order to be analysed with multivariate techniques. The method is applied to a relevant clinical orthopaedic research field, namely to describe the dynamics in the Achilles tendon and the calf during real-time movements. RESULTS: This study introduces a novel method to medical applications that can be used to examine ultrasound image sequences and to detect, visualise and quantify skeletal muscle dynamics and functions. CONCLUSIONS: This new objective method is a powerful tool to use when visualising tissue activity and dynamics of musculoskeletal ultrasound registrations.

Chemometrics in metabolomics--a review in human disease diagnosis.

Metabolomics is a post genomic research field concerned with developing methods for analysis of low molecular weight compounds in biological systems, such as cells, organs or organisms. Analyzing metabolic differences between unperturbed and perturbed systems, such as healthy volunteers and patients with a disease, can lead to insights into the underlying pathology. In metabolomics analysis, large amounts of data are routinely produced in order to characterize samples. The use of multivariate data analysis techniques and chemometrics is a commonly used strategy for obtaining reliable results. Metabolomics have been applied in different fields such as disease diagnosis, toxicology, plant science and pharmaceutical and environmental research. In this review we take a closer look at the chemometric methods used and the available results within the field of disease diagnosis. We will first present some current strategies for performing metabolomics studies, especially regarding disease diagnosis. The main focus will be on data analysis strategies and validation of multivariate models, since there are many pitfalls in this regard. Further, we highlight the most interesting metabolomics publications and discuss these in detail; additional studies are mentioned as a reference for the interested reader. A general trend is an increased focus on biological interpretation rather than merely the ability to classify samples. In the conclusions, the general trends and some recommendations for improving metabolomics data analysis are provided.

Integrated analysis of transcript, protein and metabolite data to study lignin biosynthesis in hybrid aspen.

Tree biotechnology will soon reach a mature state where it will influence the overall supply of fiber, energy and wood products. We are now ready to make the transition from identifying candidate genes, controlling important biological processes, to discovering the detailed molecular function of these genes on a broader, more holistic, systems biology level. In this paper, a strategy is outlined for informative data generation and integrated modeling of systematic changes in transcript, protein and metabolite profiles measured from hybrid aspen samples. The aim is to study characteristics of common changes in relation to genotype-specific perturbations affecting the lignin biosynthesis and growth. We show that a considerable part of the systematic effects in the system can be tracked across all platforms and that the approach has a high potential value in functional characterization of candidate genes.

Statistically integrated metabonomic-proteomic studies on a human prostate cancer xenograft model in mice.

A novel statistically integrated proteometabonomic method has been developed and applied to a human tumor xenograft mouse model of prostate cancer. Parallel 2D-DIGE proteomic and 1H NMR metabolic profile data were collected on blood plasma from mice implanted with a prostate cancer (PC-3) xenograft and from matched control animals. To interpret the xenograft-induced differences in plasma profiles, multivariate statistical algorithms including orthogonal projection to latent structure (OPLS) were applied to generate models characterizing the disease profile. Two approaches to integrating metabonomic data matrices are presented based on OPLS algorithms to provide a framework for generating models relating to the specific and common sources of variation in the metabolite concentrations and protein abundances that can be directly related to the disease model. Multiple correlations between metabolites and proteins were found, including associations between serotransferrin precursor and both tyrosine and 3-D-hydroxybutyrate. Additionally, a correlation between decreased concentration of tyrosine and increased presence of gelsolin was also observed. This approach can provide enhanced recovery of combination candidate biomarkers across multi-omic platforms, thus, enhancing understanding of in vivo model systems studied by multiple omic technologies.