Untargeted and stable isotope-assisted metabolomic analysis of MDA-MB-231 cells under hypoxia

Introduction

Hypoxia commonly occurs in cancers and is highly related with the occurrence, development and metastasis of cancer. Treatment of triple negative breast cancer remains challenge. Knowledge about the metabolic status of triple negative breast cancer cell lines in hypoxia is valuable for the understanding of molecular mechanisms of this tumor subtype to develop effective therapeutics.

Objectives

Comprehensively characterize the metabolic profiles of triple negative breast cancer cell line MDA-MB-231 in normoxia and hypoxia and the pathways involved in metabolic changes in hypoxia.

Methods

Differences in metabolic profiles affected pathways of MDA-MB-231 cells in normoxia and hypoxia were characterized using GC–MS based untargeted and stable isotope assisted metabolomic techniques.

Results

Thirty-three metabolites were significantly changed in hypoxia and nine pathways were involved. Hypoxia increased glycolysis, inhibited TCA cycle, pentose phosphate pathway and pyruvate carboxylation, while increased glutaminolysis in MDA-MB-231 cells.

Conclusion

The current results provide metabolic differences of MDA-MB-231 cells in normoxia and hypoxia conditions as well as the involved metabolic pathways, demonstrating the power of combined use of untargeted and stable isotope-assisted metabolomic methods in comprehensive metabolomic analysis.

     

A methodology for elucidating regulatory mechanisms leading to changes in lipid profiles

Introduction

It is difficult to elucidate the metabolic and regulatory factors causing lipidome perturbations.

Objectives

This work simplifies this process.

Methods

A method has been developed to query an online holistic lipid metabolic network (of 7923 metabolites) to extract the pathways that connect the input list of lipids.

Results

The output enables pathway visualisation and the querying of other databases to identify potential regulators. When used to a study a plasma lipidome dataset of polycystic ovary syndrome, 14 enzymes were identified, of which 3 are linked to ELAVL1—an mRNA stabiliser.

Conclusion

This method provides a simplified approach to identifying potential regulators causing lipid-profile perturbations.

     

Dynamic cross-talk analysis among TNF-R,TLR-4 and IL-1R signalings in TNFα-induced inflammatory responses

Background

Development in systems biology research has accelerated in recent years, and the reconstructions for molecular networks can provide a global view to enable in-depth investigation on numerous system properties in biology. However, we still lack a systematic approach to reconstruct the dynamic protein-protein association networks at different time stages from high-throughput data to further analyze the possible cross-talks among different signaling/regulatory pathways.

Methods

In this study we integrated protein-protein interactions from different databases to construct the rough protein-protein association networks (PPANs) during TNFα-induced inflammation. Next, the gene expression profiles of TNFα-induced HUVEC and a stochastic dynamic model were used to rebuild the significant PPANs at different time stages, reflecting the development and progression of endothelium inflammatory responses. A new cross-talk ranking method was used to evaluate the potential core elements in the related signaling pathways of toll-like receptor 4 (TLR-4) as well as receptors for tumor necrosis factor (TNF-R) and interleukin-1 (IL-1R).

Results

The highly ranked cross-talks which are functionally relevant to the TNFα pathway were identified. A bow-tie structure was extracted from these cross-talk pathways, suggesting the robustness of network structure, the coordination of signal transduction and feedback control for efficient inflammatory responses to different stimuli. Further, several characteristics of signal transduction and feedback control were analyzed.

Conclusions

A systematic approach based on a stochastic dynamic model is proposed to generate insight into the underlying defense mechanisms of inflammation via the construction of corresponding signaling networks upon specific stimuli. In addition, this systematic approach can be applied to other signaling networks under different conditions in different species. The algorithm and method proposed in this study could expedite prospective systems biology research when better experimental techniques for protein expression detection and microarray data with multiple sampling points become available in the future.

     

From correlation to causation: analysis of metabolomics data using systems biology approaches

Introduction

Metabolomics is a well-established tool in systems biology, especially in the top–down approach. Metabolomics experiments often results in discovery studies that provide intriguing biological hypotheses but rarely offer mechanistic explanation of such findings. In this light, the interpretation of metabolomics data can be boosted by deploying systems biology approaches.

Objectives

This review aims to provide an overview of systems biology approaches that are relevant to metabolomics and to discuss some successful applications of these methods.

Methods

We review the most recent applications of systems biology tools in the field of metabolomics, such as network inference and analysis, metabolic modelling and pathways analysis.

Results

We offer an ample overview of systems biology tools that can be applied to address metabolomics problems. The characteristics and application results of these tools are discussed also in a comparative manner.

Conclusions

Systems biology-enhanced analysis of metabolomics data can provide insights into the molecular mechanisms originating the observed metabolic profiles and enhance the scientific impact of metabolomics studies.

     

Integrating proteomic and phosphoproteomic data for pathway analysis in breast cancer

Background

As protein is the basic unit of cell function and biological pathway, shotgun proteomics, the large-scale analysis of proteins, is contributing greatly to our understanding of disease mechanisms. Proteomics study could detect the changes of both protein expression and modification. With the releases of large-scale cancer proteome studies, how to integrate acquired proteomic and phosphoproteomic data in more comprehensive pathway analysis becomes implemented, but remains challenging. Integrative pathway analysis at proteome level provides a systematic insight into the signaling network adaptations in the development of cancer.

Results

Here we integrated proteomic and phosphoproteomic data to perform pathway prioritization in breast cancer. We manually collected and curated breast cancer well-known related pathways from the literature as target pathways (TPs) or positive control in method evaluation. Three different strategies including Hypergeometric test based over-representation analysis, Kolmogorov-Smirnov (K-S) test based gene set analysis and topology-based pathway analysis, were applied and evaluated in integrating protein expression and phosphorylation. In comparison, we also assessed the ranking performance of the strategy using information of protein expression or protein phosphorylation individually. Target pathways were ranked more top with the data integration than using the information from proteomic or phosphoproteomic data individually. In the comparisons of pathway analysis strategies, topology-based method outperformed than the others. The subtypes of breast cancer, which consist of Luminal A, Luminal B, Basal and HER2-enriched, vary greatly in prognosis and require distinct treatment. Therefore we applied topology-based pathway analysis with integrating protein expression and phosphorylation profiles on four subtypes of breast cancer. The results showed that TPs were enriched in all subtypes but their ranks were significantly different among the subtypes. For instance, p53 pathway ranked top in the Basal-like breast cancer subtype, but not in HER2-enriched type. The rank of Focal adhesion pathway was more top in HER2- subtypes than in HER2+ subtypes. The results were consistent with some previous researches.

Conclusions

The results demonstrate that the network topology-based method is more powerful by integrating proteomic and phosphoproteomic in pathway analysis of proteomics study. This integrative strategy can also be used to rank the specific pathways for the disease subtypes.

     

Integrative pathway-based survival prediction utilizing the interaction between gene expression and DNA methylation in breast cancer

Background

Integrative analysis on multi-omics data has gained much attention recently. To investigate the interactive effect of gene expression and DNA methylation on cancer, we propose a directed random walk-based approach on an integrated gene-gene graph that is guided by pathway information.

Methods

Our approach first extracts a single pathway profile matrix out of the gene expression and DNA methylation data by performing the random walk over the integrated graph. We then apply a denoising autoencoder to the pathway profile to further identify important pathway features and genes. The extracted features are validated in the survival prediction task for breast cancer patients.

Results

The results show that the proposed method substantially improves the survival prediction performance compared to that of other pathway-based prediction methods, revealing that the combined effect of gene expression and methylation data is well reflected in the integrated gene-gene graph combined with pathway information. Furthermore, we show that our joint analysis on the methylation features and gene expression profile identifies cancer-specific pathways with genes related to breast cancer.

Conclusions

In this study, we proposed a DRW-based method on an integrated gene-gene graph with expression and methylation profiles in order to utilize the interactions between them. The results showed that the constructed integrated gene-gene graph can successfully reflect the combined effect of methylation features on gene expression profiles. We also found that the selected features by DA can effectively extract topologically important pathways and genes specifically related to breast cancer.

     

Proteomic profiling of HBV infected liver biopsies with different fibrotic stages

Background

Hepatitis B virus (HBV) is a global health problem, and infected patients if left untreated may develop cirrhosis and eventually hepatocellular carcinoma. This study aims to enlighten pathways associated with HBV related liver fibrosis for delineation of potential new therapeutic targets and biomarkers.

Methods

Tissue samples from 47 HBV infected patients with different fibrotic stages (F1 to F6) were enrolled for 2D-DIGE proteomic screening. Differentially expressed proteins were identified by mass spectrometry and verified by western blotting. Functional proteomic associations were analyzed by EnrichNet application.

Results

Fibrotic stage variations were observed for apolipoprotein A1 (APOA1), pyruvate kinase PKM (KPYM), glyceraldehyde 3-phospahate dehydrogenase (GAPDH), glutamate dehydrogenase (DHE3), aldehyde dehydrogenase (ALDH2), alcohol dehydrogenase (ALDH1A1), transferrin (TRFE), peroxiredoxin 3 (PRDX3), phenazine biosynthesis-like domain-containing protein (PBLD), immuglobulin kappa chain C region (IGKC), annexin A4 (ANXA4), keratin 5 (KRT5). Enrichment analysis with Reactome and Kegg databases highlighted the possible involvement of platelet release, glycolysis and HDL mediated lipid transport pathways. Moreover, string analysis revealed that HIF-1α (Hypoxia-inducible factor 1-alpha), one of the interacting partners of HBx (Hepatitis B X protein), may play a role in the altered glycolytic response and oxidative stress observed in liver fibrosis.

Conclusions

To our knowledge, this is the first protomic research that studies HBV infected fibrotic human liver tissues to investigate alterations in protein levels and affected pathways among different fibrotic stages. Observed changes in the glycolytic pathway caused by HBx presence and therefore its interactions with HIF-1α can be a target pathway for novel therapeutic purposes.

     

Effect of everolimus on the glucose metabolic pathway in mouse skeletal muscle cells (C2C12)

Introduction

Everolimus selectively inhibits mammalian target of rapamycin complex 1 (mTORC1) and exerts an antineoplastic effect. Metabolic disturbance has emerged as a common and unique side effect of everolimus.

Objectives

We used targeted metabolomic analysis to investigate the effects of everolimus on the intracellular glycometabolic pathway.

Methods

Mouse skeletal muscle cells (C2C12) were exposed to everolimus for 48 h, and changes in intracellular metabolites were determined by capillary electrophoresis time-of-flight mass spectrometry. mRNA abundance, protein expression and activity were measured for enzymes involved in glycometabolism and related pathways.

Results

Both extracellular and intracellular glucose levels increased with exposure to everolimus. Most intracellular glycometabolites were decreased by everolimus, including those involved in glycolysis and the pentose phosphate pathway, whereas no changes were observed in the tricarboxylic acid cycle. Everolimus suppressed mRNA expression of enzymes related to glycolysis, downstream of mTOR signaling enzymes and adenosine 5′-monophosphate protein kinases. The activity of key enzymes involved in glycolysis and the pentose phosphate pathway were decreased by everolimus. These results show that everolimus impairs glucose utilization in intracellular metabolism.

Conclusions

The present metabolomic analysis indicates that everolimus impairs glucose metabolism in muscle cells by lowering the activities of glycolysis and the pentose phosphate pathway.

     

Pathway-based analyses

Background

New technologies for acquisition of genomic data, while offering unprecedented opportunities for genetic discovery, also impose severe burdens of interpretation andpenalties for multiple testing.

Methods

The Pathway-based Analyses Group of the Genetic Analysis Workshop 19 (GAW19) sought reduction of multiple-testing burden through various approaches to aggregation of highdimensional data in pathways informed by prior biological knowledge.

Results

Experimental methods testedincluded the use of "synthetic pathways" (random sets of genes) to estimate power and false-positive error rate of methods applied to simulated data; data reduction via independent components analysis, single-nucleotide polymorphism (SNP)-SNP interaction, and use of gene sets to estimate genetic similarity; and general assessment of the efficacy of prior biological knowledge to reduce the dimensionality of complex genomic data.

Conclusions

The work of this group explored several promising approaches to managing high-dimensional data, with the caveat that these methods are necessarily constrained by the quality of external bioinformatic annotation.

     

Characterizing biomarkers in osteosarcoma metastasis based on an ego-network

Objectives

To characterize biomarkers that underlie osteosarcoma (OS) metastasis based on an ego-network.

Results

From the microarray data, we obtained 13,326 genes. By combining PPI data and microarray data, 10,520 shared genes were found and constructed into ego-networks. 17 significant ego-networks were identified with p < 0.05. In the pathway enrichment analysis, seven ego-networks were identified with the most significant pathway.

Conclusions

These significant ego-modules were potential biomarkers that reveal the potential mechanisms in OS metastasis, which may contribute to understanding cancer prognoses and providing new perspectives in the treatment of cancer.

     

The Black Yeasts: an Update on Species Identification and Diagnosis

Purpose of review

Black yeast-like fungi are capable of causing a wide range of infections, including invasive disease. The diagnosis of infections caused by these species can be problematic. We review the changes in the nomenclature and taxonomy of these fungi, and methods used for detection and species identification that aid in diagnosis.

Recent findings

Molecular assays, including DNA barcode analysis and rolling circle amplification, have improved our ability to correctly identify these species. A proteomic approach using matrix-assisted laser desorption/ionization time-of-flight mass spectrometry (MALDI-TOF MS) has also shown promising results. While progress has been made with molecular techniques using direct specimens, data are currently limited.

Summary

Molecular and proteomic assays have improved the identification of black yeast-like fungi. However, improved molecular and proteomic databases and better assays for the detection and identification in direct specimens are needed to improve the diagnosis of disease caused by black yeast-like fungi.

     

bFGF-mediated pluripotency maintenance in human induced pluripotent stem cells is associated with NRAS-MAPK signaling

Background

Human pluripotent stem cells (PSCs) open new windows for basic research and regenerative medicine due to their remarkable properties, i.e. their ability to self-renew indefinitely and being pluripotent. There are different, conflicting data related to the role of basic fibroblast growth factor (bFGF) in intracellular signal transduction and the regulation of pluripotency of PSCs. Here, we investigated the effect of bFGF and its downstream pathways in pluripotent vs. differentiated human induced (hi) PSCs.

Methods

bFGF downstream signaling pathways were investigated in long-term culture of hiPSCs from pluripotent to differentiated state (withdrawing bFGF) using immunoblotting, immunocytochemistry and qPCR. Subcellular distribution of signaling components were investigated by simple fractionation and immunoblotting upon bFGF stimulation. Finally, RAS activity and RAS isoforms were studied using RAS assays both after short- and long-term culture in response to bFGF stimulation.

Results

Our results revealed that hiPSCs were differentiated into the ectoderm lineage upon withdrawing bFGF as an essential pluripotency mediator. Pluripotency markers OCT4, SOX2 and NANOG were downregulated, following a drastic decrease in MAPK pathway activity levels. Notably, a remarkable increase in phosphorylation levels of p38 and JAK/STAT3 was observed in differentiated hiPSCs, while the PI3K/AKT and JNK pathways remained active during differentiation. Our data further indicate that among the RAS paralogs, NRAS predominantly activates the MAPK pathway in hiPSCs.

Conclusion

Collectively, the MAPK pathway appears to be the prime signaling pathway downstream of bFGF for maintaining pluripotency in hiPSCs and among the MAPK pathways, the activity of NRAS-RAF-MEK-ERK is decreased during differentiation, whereas p38 is activated and JNK remains constant.

     

A lost opportunity for science: journals promote data sharing in metabolomics but do not enforce it

Introduction

Data sharing is being increasingly required by journals and has been heralded as a solution to the ‘replication crisis’.

Objectives

(i) Review data sharing policies of journals publishing the most metabolomics papers associated with open data and (ii) compare these journals’ policies to those that publish the most metabolomics papers.

Methods

A PubMed search was used to identify metabolomics papers. Metabolomics data repositories were manually searched for linked publications.

Results

Journals that support data sharing are not necessarily those with the most papers associated to open metabolomics data.

Conclusion

Further efforts are required to improve data sharing in metabolomics.

     

Generalized logical model based on network topology to capture the dynamical trends of cellular signaling pathways

Background

Cellular responses to extracellular perturbations require signaling pathways to capture and transmit the signals. However, the underlying molecular mechanisms of signal transduction are not yet fully understood, thus detailed and comprehensive models may not be available for all the signaling pathways. In particular, insufficient knowledge of parameters, which is a long-standing hindrance for quantitative kinetic modeling necessitates the use of parameter-free methods for modeling and simulation to capture dynamic properties of signaling pathways.

Results

We present a computational model that is able to simulate the graded responses to degradations, the sigmoidal biological relationships between signaling molecules and the effects of scheduled perturbations to the cells. The simulation results are validated using experimental data of protein phosphorylation, demonstrating that the proposed model is capable of capturing the main trend of protein activities during the process of signal transduction. Compared with existing simulators, our model has better performance on predicting the state transitions of signaling networks.

Conclusion

The proposed simulation tool provides a valuable resource for modeling cellular signaling pathways using a knowledge-based method.

     

Comprehensive annotation of microRNA expression profiles

Background

MicroRNAs (miRNAs) regulate many biological processes by post-translational gene silencing. Analysis of miRNA expression profiles is a reliable method for investigating particular biological processes due to the stability of miRNA and the development of advanced sequencing methods. However, this approach is limited by the broad specificity of miRNAs, which may target several mRNAs.

Result

In this study, we developed a method for comprehensive annotation of miRNA array or deep sequencing data for investigation of cellular biological effects. Using this method, the specific pathways and biological processes involved in Alzheimer’s disease were predicted with high correlation in four independent samples. Furthermore, this method was validated for evaluation of cadmium telluride (CdTe) nanomaterial cytotoxicity. As a result, apoptosis pathways were selected as the top pathways associated with CdTe nanoparticle exposure, which is consistent with previous studies.

Conclusions

Our findings contribute to the validation of miRNA microarray or deep sequencing results for early diagnosis of disease and evaluation of the biological safety of new materials and drugs.

     

Feed-forward neural network model for hunger and satiety related VAS score prediction

Background

An artificial neural network approach was chosen to model the outcome of the complex signaling pathways in the gastro-intestinal tract and other peripheral organs that eventually produce the satiety feeling in the brain upon feeding.

Methods

A multilayer feed-forward neural network was trained with sets of experimental data relating concentration-time courses of plasma satiety hormones to Visual Analog Scales (VAS) scores. The network successfully predicted VAS responses from sets of satiety hormone data obtained in experiments using different food compositions.

Results

The correlation coefficients for the predicted VAS responses for test sets having i) a full set of three satiety hormones, ii) a set of only two satiety hormones, and iii) a set of only one satiety hormone were 0.96, 0.96, and 0.89, respectively. The predicted VAS responses discriminated the satiety effects of high satiating food types from less satiating food types both in orally fed and ileal infused forms.

Conclusions

From this application of artificial neural networks, one may conclude that neural network models are very suitable to describe situations where behavior is complex and incompletely understood. However, training data sets that fit the experimental conditions need to be available.

     

Group-wise ANOVA simultaneous component analysis for designed <Emphasis Type="Italic">omics</Emphasis> experiments

Introduction

Modern omics experiments pertain not only to the measurement of many variables but also follow complex experimental designs where many factors are manipulated at the same time. This data can be conveniently analyzed using multivariate tools like ANOVA-simultaneous component analysis (ASCA) which allows interpretation of the variation induced by the different factors in a principal component analysis fashion. However, while in general only a subset of the measured variables may be related to the problem studied, all variables contribute to the final model and this may hamper interpretation.

Objectives

We introduce here a sparse implementation of ASCA termed group-wise ANOVA-simultaneous component analysis (GASCA) with the aim of obtaining models that are easier to interpret.

Methods

GASCA is based on the concept of group-wise sparsity introduced in group-wise principal components analysis where structure to impose sparsity is defined in terms of groups of correlated variables found in the correlation matrices calculated from the effect matrices.

Results

The GASCA model, containing only selected subsets of the original variables, is easier to interpret and describes relevant biological processes.

Conclusions

GASCA is applicable to any kind of omics data obtained through designed experiments such as, but not limited to, metabolomic, proteomic and gene expression data.

     

Quantitative analysis of tetrahydrofolate metabolites from <Emphasis Type="Italic">clostridium autoethanogenum</Emphasis>

Introduction

Quantification of tetrahydrofolates (THFs), important metabolites in the Wood–Ljungdahl pathway (WLP) of acetogens, is challenging given their sensitivity to oxygen.

Objective

To develop a simple anaerobic protocol to enable reliable THFs quantification from bioreactors.

Methods

Anaerobic cultures were mixed with anaerobic acetonitrile for extraction. Targeted LC–MS/MS was used for quantification.

Results

Tetrahydrofolates can only be quantified if sampled anaerobically. THF levels showed a strong correlation to acetyl-CoA, the end product of the WLP.

Conclusion

Our method is useful for relative quantification of THFs across different growth conditions. Absolute quantification of THFs requires the use of labelled standards.