Biomarker Discovery for Hepatocellular Carcinoma in Patients with Liver Cirrhosis Using Untargeted Metabolomics and Lipidomics Studies

Hepatocellular carcinoma (HCC), the most prevalent form of liver cancer, is the third leading cause of mortality globally. Patients with HCC have a poor prognosis due to the fact that the emergence of symptoms typically occurs at a late stage of the disease. In addition, conventional biomarkers perform suboptimally when identifying HCC in its early stages, heightening the need for the identification of new and more effective biomarkers. Using metabolomics and lipidomics approaches, this study aims to identify serum biomarkers for identification of HCC in patients with liver cirrhosis (LC). Serum samples from 20 HCC cases and 20 patients with LC were analyzed using ultra-high-performance liquid chromatography-Q Exactive mass spectrometry (UHPLC-Q-Exactive-MS). Metabolites and lipids that are significantly altered between HCC cases and patients with LC were identified. These include organic acids, amino acids, TCA cycle intermediates, fatty acids, bile acids, glycerophospholipids, sphingolipids, and glycerolipids. The most significant variability was observed in the concentrations of bile acids, fatty acids, and glycerophospholipids. In the context of HCC cases, there was a notable increase in the levels of phosphatidylethanolamine and triglycerides, but the levels of fatty acids and phosphatidylcholine exhibited a substantial decrease. In addition, it was observed that all of the identified metabolites exhibited a superior area under the receiver operating characteristic (ROC) curve in comparison to alpha-fetoprotein (AFP). The pathway analysis of these metabolites revealed fatty acid, lipid, and energy metabolism as the most impacted pathways. Putative biomarkers identified in this study will be validated in future studies via targeted quantification.

Deconvolution of immune cell composition and biological age of hepatocellular carcinoma using DNA methylation

Hepatocellular carcinoma (HCC) has been an approved indication for the administration of immunotherapy since 2017, but biomarkers that predict therapeutic response have remained limited. Understanding and characterizing the tumor immune microenvironment enables better classification of these tumors and may reveal biomarkers that predict immunotherapeutic efficacy. In this paper, we applied a cell-type deconvolution algorithm using DNA methylation array data to investigate the composition of the tumor microenvironment in HCC. Using publicly available and in-house datasets with a total cohort size of 57 patients, each with tumor and matched normal tissue samples, we identified key differences in immune cell composition. We found that NK cell abundance was significantly decreased in HCC tumors compared to adjacent normal tissue. We also applied DNA methylation "clocks" which estimate phenotypic aging and compared these findings to expression-based determinations of cellular senescence. Senescence and epigenetic aging were significantly increased in HCC tumors, and the degree of age acceleration and senescence was strongly associated with decreased NK cell abundance. In summary, we found that NK cell infiltration in the tumor microenvironment is significantly diminished, and that this loss of NK abundance is strongly associated with increased senescence and age-related phenotype. These findings point to key interactions between NK cells and the senescent tumor microenvironment and offer insights into the pathogenesis of HCC as well as potential biomarkers of therapeutic efficacy.

Formalin Fixation, Delay to Fixation, and Time in Fixative Adversely Impact Copy Number Variation Analysis by aCGH

Although molecular profiling of DNA isolated from formalin-fixed, paraffin-embedded (FFPE) tumor specimens has become more common in recent years, it remains unclear how discrete FFPE processing variables may affect detection of copy number variation (CNV). To better understand such effects, array comparative genomic hybridization (aCGH) profiles of FFPE renal cell carcinoma specimens that experienced different delays to fixation (DTFs; 1, 2, 3, and 12 hours) and times in fixative (TIFs; 6, 12, 23, and 72 hours) were compared to snap-frozen tumor and blood specimens from the same patients. A greater number of regions containing CNVs relative to commercial reference DNA were detected in DNA from FFPE tumor specimens than snap-frozen tumor specimens even though they originated from the same tumor blocks. Extended DTF and TIF affected the number of DNA segments with a copy number status that differed between FFPE and frozen tumor specimens; a DTF ≥3 hours led to more segments, while a TIF of 72 hours led to fewer segments. Importantly, effects were not random as a higher guanine-cytosine (GC) content and/or a higher percentage of repeats were observed among stable regions. While limiting aCGH analysis to FFPE specimens with a DTF <3 hours and a TIF <72 hours may circumvent some effects, results from FFPE specimens should be validated against fresh or frozen specimens whenever possible.

Deep Learning Based Metabolite Annotation

Metabolite annotation is a major bottleneck in untargeted metabolomics studies by liquid chromatography coupled with mass spectrometry (LC-MS). This is in part due to the limited publicly available spectral libraries, which consist of tandem mass spectrometry (MS/MS) data acquired from just a fraction of known compounds. Machine learning and deep learning methods provide the opportunity to predict molecular fingerprints based on MS/MS data. The predicted molecular fingerprints can then be used to help rank candidate metabolite IDs obtained based on predicted formula or measured precursor m/z of the unknown metabolite. This approach is particularly useful to help annotate metabolites whose corresponding MS/MS spectra cannot be matched with those in spectral libraries. We previously reported application of a convolutional neural network (CNN) for molecular fingerprint prediction using MS/MS spectra obtained from the MoNA repository and NIST 20. In this paper, we investigate high-dimensional representation of the spectral data and molecular fingerprints to improve accuracy in molecular fingerprint prediction.

Cell-type Deconvolution and Age Estimation Using DNA Methylation Reveals NK Cell Deficiency in the Hepatocellular Carcinoma Microenvironment

Hepatocellular carcinoma (HCC) has been an approved indication for the administration of immunotherapy since 2017, but biomarkers that predict therapeutic response have remained limited. Understanding and characterizing the tumor immune microenvironment enables better classification of these tumors and may reveal biomarkers that predict immunotherapeutic efficacy. In this paper, we applied a cell-type deconvolution algorithm using DNA methylation array data to investigate the composition of the tumor microenvironment in HCC. Using two publicly available datasets with a total cohort size of 57 patients, each with tumor and matched normal tissue samples, we identified key differences in immune cell composition. We found that NK cell abundance was significantly decreased in HCC tumors compared to adjacent normal tissue. We also applied DNA methylation "clocks" which estimate phenotypic aging and compared these findings to expression-based determinations of cellular senescence. Senescence and epigenetic aging was significantly increased in HCC tumors, and the degree of age acceleration and senescence was strongly associated with decreased NK cell abundance. In summary, we found that NK cell infiltration in the tumor microenvironment is significantly diminished, and that this loss of NK abundance is strongly associated with increased senescence and age-related phenotype. These findings point to key interactions between NK cells and the senescent tumor microenvironment and offer insights into the pathogenesis of HCC as well as potential biomarkers of therapeutic efficacy.

A Bayesian two-step integrative procedure incorporating prior knowledge for the identification of miRNA-mRNAs involved in hepatocellular carcinoma

Recent studies have confirmed the role of miRNA regulation of gene expression in oncogenesis for various cancers. In parallel, prior knowledge about relationships between miRNA and mRNA have been accumulated from biological experiments or statistical analyses. Improved identification of disease-associated miRNA-mRNA pairs may be achieved by incorporating prior knowledge into integrative genomic analyses. In this study we focus on 39 patients with hepatocellular carcinoma (HCC) and 25 patients with liver cirrhosis and use a flexible Bayesian two-step integrative method. We found 66 significant miRNA-mRNA pairs, several of which contain molecules that have previously been identified as potential biomarkers. These results demonstrate the utility of the proposed approach in providing a better understanding of relationships between different biological levels, thereby giving insights into the biological mechanisms underlying the diseases, while providing a better selection of biomarkers that may serve as diagnostic, prognostic, or therapeutic biomarker candidates.

Metabolomic Analysis of Plasma from Breast Cancer Patients Using Ultra-High-Performance Liquid Chromatography Coupled with Mass Spectrometry: An Untargeted Study

Breast cancer (BC) is one of the leading causes of cancer mortality in women worldwide, and therefore, novel biomarkers for early disease detection are critically needed. We performed herein an untargeted plasma metabolomic profiling of 55 BC patients and 55 healthy controls (HC) using ultra-high performance liquid chromatography coupled with quadrupole time-of-flight mass spectrometry (UHPLC/Q-TOF-MS). Pre-processed data revealed 2494 ions in total. Data matrices’ paired t-tests revealed 792 ions (both positive and negative) which presented statistically significant changes (FDR < 0.05) in intensity levels between cases versus controls. Metabolites identified with putative names via MetaboQuest using MS/MS and mass-based approaches included amino acid esters (i.e., N-stearoyl tryptophan, L-arginine ethyl ester), dipeptides (ile-ser, met-his), nitrogenous bases (i.e., uracil derivatives), lipid metabolism-derived molecules (caproleic acid), and exogenous compounds from plants, drugs, or dietary supplements. LASSO regression selected 16 metabolites after several variables (TNM Stage, Grade, smoking status, menopausal status, and race) were adjusted. A predictive conditional logistic regression model on the 16 LASSO selected ions provided a high diagnostic performance with an area-under-the-curve (AUC) value of 0.9729 (95% CI 0.96−0.98) on all 55 samples. This study proves that BC possesses a specific metabolic signature that could be exploited as a novel metabolomics-based approach for BC detection and characterization. Future studies of large-scale cohorts are needed to validate these findings.

Integrative Analysis of DNA Methylation and microRNA Expression Reveals Mechanisms of Racial Heterogeneity in Hepatocellular Carcinoma

Pathologic alterations in epigenetic regulation have long been considered a hallmark of many cancers, including hepatocellular carcinoma (HCC). In a healthy individual, the relationship between DNA methylation and microRNA (miRNA) expression maintains a fine balance; however, disruptions in this harmony can aid in the genesis of cancer or the propagation of existing cancers. The balance between DNA methylation and microRNA expression and its potential disturbance in HCC can vary by race. There is emerging evidence linking epigenetic events including DNA methylation and miRNA expression to cancer disparities. In this paper, we evaluate the epigenetic mechanisms of racial heterogenity in HCC through an integrated analysis of DNA methylation, miRNA, and combined regulation of gene expression. Specifically, we generated DNA methylation, mRNA-seq, and miRNA-seq data through the analysis of tumor and adjacent non-tumor liver tissues from African Americans (AA) and European Americans (EA) with HCC. Using mixed ANOVA, we identified cytosine-phosphate-guanine (CpG) sites, mRNAs, and miRNAs that are significantly altered in HCC vs. adjacent non-tumor tissue in a race-specific manner. We observed that the methylome was drastically changed in EA with a significantly larger number of differentially methylated and differentially expressed genes than in AA. On the other hand, the miRNA expression was altered to a larger extent in AA than in EA. Pathway analysis functionally linked epigenetic regulation in EA to processes involved in immune cell maturation, inflammation, and vascular remodeling. In contrast, cellular proliferation, metabolism, and growth pathways are found to predominate in AA as a result of this epigenetic analysis. Furthermore, through integrative analysis, we identified significantly differentially expressed genes in HCC with disparate epigenetic regulation, associated with changes in miRNA expression for AA and DNA methylation for EA.

Integrative Analysis to Identify Race-Associated Metabolite Biomarkers for Hepatocellular Carcinoma

Compared to European-Americans (EAs), the incidence of hepatocellular carcinoma (HCC) is higher in African-Americans (AAs) and is associated with more advanced tumor stage at diagnosis and lower survival rates. The increasing burden makes discovery of novel diagnostic, prognostic, and therapeutic biomarkers distinguishing HCC from underlying cirrhosis a significant focus. In this study, we analyzed tissue and serum samples from 40 HCC cases and 25 patients with liver cirrhosis to identify candidate biomarkers that distinguish HCC from cirrhotic patients in a race specific manner. Through integrative analysis of transcriptomic and metabolomic data, we investigated candidate metabolite biomarkers that are specific to AAs and EAs. The results from this demonstrate the utility of integrating transcriptomic and metabolomic data to prioritize clinically and biologically relevant metabolite biomarkers that can increase understanding of molecular mechanisms driving HCC in different racial groups.

Epigenetic changes associated with mechanisms of disparities in hepatocellular carcinoma

In addition to socioeconomic influences, biological factors are believed to play a role in health disparities. In this paper, we investigate miRNA, mRNA, and DNA methylation patterns that contribute to disparities in hepatocellular carcinoma (HCC). This is accomplished by integration of mRNA-Seq, miRNA-Seq, and DNA methylation data we acquired by analysis of liver tissues from 30 HCC patients consisting of European Americans (EAs), African Americans (AAs), and Asian Americans (Asians). Mixed-ANOVA models are applied to identify miRNAs, mRNAs, and DNA methylation sites that are significantly altered in tumor vs. adjacent normal tissues in a race-specific manner. Through integrated analysis, a refined list of differentially expressed mRNAs is obtained by selecting those that are targets of differentially expressed miRNAs and consist of promoter regions that are differentially methylated.

Identification of miRNA-mRNA associations in hepatocellular carcinoma using hierarchical integrative model

Background

The established role miRNA-mRNA regulation of gene expression has in oncogenesis highlights the importance of integrating miRNA with downstream mRNA targets. These findings call for investigations aimed at identifying disease-associated miRNA-mRNA pairs. Hierarchical integrative models (HIM) offer the opportunity to uncover the relationships between disease and the levels of different molecules measured in multiple omic studies.

Methods

The HIM model we formulated for analysis of mRNA-seq and miRNA-seq data can be specified with two levels: (1) a mechanistic submodel relating mRNAs to miRNAs, and (2) a clinical submodel relating disease status to mRNA and miRNA, while accounting for the mechanistic relationships in the first level.

Results

mRNA-seq and miRNA-seq data were acquired by analysis of tumor and normal liver tissues from 30 patients with hepatocellular carcinoma (HCC). We analyzed the data using HIM and identified 157 significant miRNA-mRNA pairs in HCC. The majority of these molecules have already been independently identified as being either diagnostic, prognostic, or therapeutic biomarker candidates for HCC. These pairs appear to be involved in processes contributing to the pathogenesis of HCC involving inflammation, regulation of cell cycle, apoptosis, and metabolism. For further evaluation of our method, we analyzed miRNA-seq and mRNA-seq data from TCGA network. While some of the miRNA-mRNA pairs we identified by analyzing both our and TCGA data are previously reported in the literature and overlap in regulation and function, new pairs have been identified that may contribute to the discovery of novel targets.

Conclusion

The results strongly support the hypothesis that miRNAs are important regulators of mRNAs in HCC. Furthermore, these results emphasize the biological relevance of studying miRNA-mRNA pairs.

Multi-omic Pathway and Network Analysis to Identify Biomarkers for Hepatocellular Carcinoma

The threat of Hepatocellular Carcinoma (HCC) is a growing problem, with incidence rates anticipated to near double over the next two decades. The increasing burden makes discovery of novel diagnostic, prognostic, and therapeutic biomarkers distinguishing HCC from underlying cirrhosis a significant focus. In this study, we analyzed tissue and serum samples from 40 HCC cases and 25 patients with liver cirrhosis (CIRR) to better understand the mechanistic differences between HCC and CIRR. Through pathway and network analysis, we are able to take a systems biology approach to conduct multi-omic analysis of transcriptomic, glycoproteomic, and metabolomic data acquired through various platforms. As a result, we are able to identify the FXR/RXR Activation pathway as being represented by molecules spanning multiple molecular compartments in these samples. Specifically, serum metabolites deoxycholate and chenodeoxycholic acid and serum glycoproteins C4A/C4B, KNG1, and HPX are biomarker candidates identified from this analysis that are of interest for future targeted studies. These results demonstrate the integrative power of multi-omic analysis to prioritize clinically and biologically relevant biomarker candidates that can increase understanding of molecular mechanisms driving HCC and make an impact in patient care.

Metabolomic Analysis of Liver Tissues for Characterization of Hepatocellular Carcinoma

Hepatocellular carcinoma (HCC) causes more than half a million annual deaths worldwide. Understanding the mechanisms contributing to HCC development is highly desirable for improved surveillance, diagnosis, and treatment. Liver tissue metabolomics has the potential to reflect the physiological changes behind HCC development. Also, it allows identification of biomarker candidates for future evaluation in biofluids and investigation of racial disparities in HCC. Tumor and nontumor tissues from 40 patients were analyzed by both gas chromatography-mass spectrometry (GC-MS) and liquid chromatography-mass spectrometry (LC-MS) platforms to increase the metabolome coverage. The levels of the metabolites extracted from solid liver tissue of the HCC area and adjacent non-HCC area were compared. Among the analytes detected by GC-MS and LC-MS with significant alterations, 18 were selected based on biological relevance and confirmed metabolite identification. These metabolites belong to TCA cycle, glycolysis, purines, and lipid metabolism and have been previously reported in liver metabolomic studies where high correlation with HCC progression is implied. We demonstrated that metabolites related to HCC pathogenesis can be identified through liver tissue metabolomic analysis. Additionally, this study has enabled us to identify race-specific metabolites associated with HCC.

Discovery of potential urine-accessible metabolite biomarkers associated with muscle disease and corticosteroid response in the mdx mouse model for Duchenne

Urine is increasingly being considered as a source of biomarker development in Duchenne Muscular Dystrophy (DMD), a severe, life-limiting disorder that affects approximately 1 in 4500 boys. In this study, we considered the mdx mice-a murine model of DMD-to discover biomarkers of disease, as well as pharmacodynamic biomarkers responsive to prednisolone, a corticosteroid commonly used to treat DMD. Longitudinal urine samples were analyzed from male age-matched mdx and wild-type mice randomized to prednisolone or vehicle control via liquid chromatography tandem mass spectrometry. A large number of metabolites (869 out of 6,334) were found to be significantly different between mdx and wild-type mice at baseline (Bonferroni-adjusted p-value < 0.05), thus being associated with disease status. These included a metabolite with m/z = 357 and creatine, which were also reported in a previous human study looking at serum. Novel observations in this study included peaks identified as biliverdin and hypusine. These four metabolites were significantly higher at baseline in the urine of mdx mice compared to wild-type, and significantly changed their levels over time after baseline. Creatine and biliverdin levels were also different between treated and control groups, but for creatine this may have been driven by an imbalance at baseline. In conclusion, our study reports a number of biomarkers, both known and novel, which may be related to either the mechanisms of muscle injury in DMD or prednisolone treatment.

INDEED: R package for network based differential expression analysis

With recent advancement of omics technologies, fueled by decreased cost and increased number of available datasets, computational methods for differential expression analysis are sought to identify disease-associated biomolecules. Conventional differential expression analysis methods (e.g. student's t-test, ANOVA) focus on assessing mean and variance of biomolecules in each biological group. On the other hand, network-based approaches take into account the interactions between biomolecules in choosing differentially expressed ones. These interactions are typically evaluated by correlation methods that tend to generate over-complicated networks due to many seemingly indirect associations. In this paper, we introduce a new R/Bioconductor package INDEED that allows users to construct a sparse network based on partial correlation, and to identify biomolecules that have significant changes both at individual expression and pairwise interaction levels. We applied INDEED for analysis of two omic datasets acquired in a cancer biomarker discovery study to help rank disease-associated biomolecules. We believe biomolecules selected by INDEED lead to improved sensitivity and specificity in detecting disease status compared to those selected by conventional statistical methods. Also, INDEED's framework is amenable to further expansion to integrate networks from multi-omic studies, thereby allowing selection of reliable disease-associated biomolecules or disease biomarkers.

Multi-omic approaches for characterization of hepatocellular carcinoma

Multi-omic approaches offer the opportunity to characterize complex diseases such as cancer at various molecular levels. In this paper, we present transcriptomic, proteomic/glycoproteomic, glycomic, and metabolomic (TPGM) data we acquired by analysis of liver tissues from hepatocellular carcinoma (HCC) cases and patients with liver cirrhosis. We evaluated changes in the levels of transcripts, proteins, glycans, and metabolites between tumor and cirrhotic tissues by statistical methods. We demonstrated the potential of multi-omic approaches and network analysis to investigate the interactions among these biomolecules in the progression of liver cirrhosis to HCC. Also, we showed the significance of multi-omic approaches to identify pathways altered in HCC.

Protein network construction using reverse phase protein array data

In this paper, we introduce a novel computational method for constructing protein networks based on reverse phase protein array (RPPA) data to identify complex patterns in protein signaling. The method is applied to phosphoproteomic profiles of basal expression and activation/phosphorylation of 76 key signaling proteins in three breast cancer cell lines (MCF7, LCC1, and LCC9). Temporal RPPA data are acquired at 48h, 96h, and 144h after knocking down four genes in separate experiments. These genes are selected from a previous study as important determinants for breast cancer survival. Interaction networks are constructed by analyzing the expression levels of protein pairs using a multivariate analysis of variance model. A new scoring criterion is introduced to determine relevant protein pairs. Through a network topology based analysis, we search for wiring patterns to identify key proteins that are associated with significant changes in expression levels across various experimental conditions.

Metabolomic Characterization of Hepatocellular Carcinoma in Patients with Liver Cirrhosis for Biomarker Discovery

Background: Metabolomics plays an important role in providing insight into the etiology and mechanisms of hepatocellular carcinoma (HCC). This is accomplished by a comprehensive analysis of patterns involved in metabolic alterations in human specimens. This study compares the levels of plasma metabolites in HCC cases versus cirrhotic patients and evaluates the ability of candidate metabolites in distinguishing the two groups. Also, it investigates the combined use of metabolites and clinical covariates for detection of HCC in patients with liver cirrhosis.Methods: Untargeted analysis of metabolites in plasma from 128 subjects (63 HCC cases and 65 cirrhotic controls) was conducted using gas chromatography coupled to mass spectrometry (GC-MS). This was followed by targeted evaluation of selected metabolites. LASSO regression was used to select a set of metabolites and clinical covariates that are associated with HCC. The performance of candidate biomarkers in distinguishing HCC from cirrhosis was evaluated through a leave-one-out cross-validation based on area under the receiver operating characteristics (ROC) curve.Results: We identified 11 metabolites and three clinical covariates that differentiated HCC cases from cirrhotic controls. Combining these features in a panel for disease classification using support vector machines (SVM) yielded better area under the ROC curve compared with alpha-fetoprotein (AFP).Conclusions: This study demonstrates the combination of metabolites and clinical covariates as an effective approach for early detection of HCC in patients with liver cirrhosis.Impact: Further investigation of these findings may improve understanding of HCC pathophysiology and possible implication of the metabolites in HCC prevention and diagnosis. Cancer Epidemiol Biomarkers Prev; 26(5); 675-83. ©2016 AACR.

Metabolomic profiling of breast tumors using ductal fluid

Identification of new biomarkers for breast cancer remains critical in order to enhance early detection of the disease and improve its prognosis. Towards this end, we performed an untargeted metabolomic analysis of breast ductal fluid using an ultra-performance liquid chromatography coupled with a quadrupole time-of-light (UPLC-QTOF) mass spectrometer. We investigated the metabolomic profiles of breast tumors using ductal fluid samples collected by ductal lavage (DL). We studied fluid from both the affected breasts and the unaffected contralateral breasts (as controls) from 43 women with confirmed unilateral breast cancer. Using this approach, we identified 1560 ions in the positive mode and 538 ions in the negative mode after preprocessing of the UPLC‑QTOF data. Paired t-tests applied on these data matrices identified 209 ions (positive and negative modes combined) with significant change in intensity level between affected and unaffected control breasts (adjusted p-values <0.05). Among these, 83 ions (39.7%) showed a fold change (FC) >1.2 and 66 ions (31.6%) were identified with putative compound names. The metabolites that we identified included endogenous metabolites such as amino acid derivatives (N-Acetyl-DL-tryptophan) or products of lipid metabolism such as N-linoleoyl taurine, trans-2-dodecenoylcarnitine, lysophosphatidylcholine LysoPC(18:2(9Z,12Z)), glycerophospholipids PG(18:0/0:0), and phosphatidylserine PS(20:4(5Z,8Z,11Z,14Z). Generalized LASSO regression further selected 21 metabolites when race, menopausal status, smoking, grade and TNM stage were adjusted for. A predictive conditional logistic regression model, using the LASSO selected 21 ions, provided diagnostic accuracy with the area under the curve of 0.956 (sensitivity/specificity of 0.907/0.884). This is the first study that shows the feasibility of conducting a comprehensive metabolomic profiling of breast tumors using breast ductal fluid to detect changes in the cellular microenvironment of the tumors and shows the potential for this approach to be used to improve detection of breast cancer.

INDEED: Integrated differential expression and differential network analysis of omic data for biomarker discovery

Differential expression (DE) analysis is commonly used to identify biomarker candidates that have significant changes in their expression levels between distinct biological groups. One drawback of DE analysis is that it only considers the changes on single biomolecule level. Recently, differential network (DN) analysis has become popular due to its capability to measure the changes on biomolecular pair level. In DN analysis, network is typically built based on correlation and biomarker candidates are selected by investigating the network topology. However, correlation tends to generate over-complicated networks and the selection of biomarker candidates purely based on network topology ignores the changes on single biomolecule level. In this paper, we propose a novel approach, INDEED, that builds sparse differential network based on partial correlation and integrates DE and DN analyses for biomarker discovery. We applied this approach on real proteomic and glycomic data generated by liquid chromatography coupled with mass spectrometry for hepatocellular carcinoma (HCC) biomarker discovery study. For each omic data, we used one dataset to select biomarker candidates, built a disease classifier and evaluated the performance of the classifier on an independent dataset. The biomarker candidates, selected by INDEED, were more reproducible across independent datasets, and led to a higher classification accuracy in predicting HCC cases and cirrhotic controls compared with those selected by separate DE and DN analyses. INDEED also identified some candidates previously reported to be relevant to HCC, such as intercellular adhesion molecule 2 (ICAM2) and c4b-binding protein alpha chain (C4BPA), which were missed by both DE and DN analyses. In addition, we applied INDEED for survival time prediction based on transcriptomic data acquired by analysis of samples from breast cancer patients. We selected biomarker candidates and built a regression model for survival time prediction based on a gene expression dataset and patients' survival records. We evaluated the performance of the regression model on an independent dataset. Compared with the biomarker candidates selected by DE and DN analyses, those selected through INDEED led to more accurate survival time prediction.

Abstract 450: Serum biomarkers in patients treated with stereotactic body radiation therapy (SBRT) for prostate cancer

Purpose: The goal of this study was to determine the feasibility for developing prognostic protein biomarkers in serum samples from patients undergoing Stereotactic Body Radiation Therapy (SBRT) for organ confined prostate cancers.

Methods: 130 patients presenting with organ confined prostate cancers were treated with SBRT to doses of 35-36.25 Gy in 5 fractions. Peripheral blood samples were collected prospectively from patients at time 0 prior to radiation and serially after the first treatment (24 hours), follow-up months 1, 3, 6, 9, and 12 during the first year, and every 6 months thereafter, for up to 3 years. Processed study samples were analyzed by SomaLogic, Inc., using the SOMAscan Version 3 proteomic assay. Statistical analysis was performed on log-transformed data, with statistically significant differences identified by evaluating false discovery rate corrected p-values. Protein correlations were discovered with the Jonckheere-Terpstra (JT) test. Ingenuity Pathway Analysis (IPA) software was used to analyze cellular signaling pathways. PSA levels and clinical recurrence information were prospectively obtained at each follow-up visit and maintained in an institutional database.

Results:

Patient stratification by risk assessment identified 27 low-, 71 intermediate- and 32 high-risk patients in the study cohort. Proteins that function in cell proliferation, angiogenesis, protein signaling, gonad development, and cell migration correlated with Gleason's grade. CGA.LHB, KLK3, and CNTFR correlated both with tumor stage and Gleason's grade. With a median follow up of 3 years, ten patients experienced biochemical failures. While no proteins identified as differentially expressed in recurrent patients achieved significance, IPA pathway analysis of the top differentially expressed proteins converged on the IL-6 canonical pathway.

183 proteins were attributed to radiation effects on differential expression in patients by comparing pre-treatment to the 3 months post-treatment specimens. IPA network pathway analyses of these paired samples revealed changes in cell activation and signaling pathways, with the primary regulatory networks associated with ERK1/2, NF-κB, IFNG, ADIPOQ, histone 3, and histone 4. Of particular interest, high adiponectin levels in patients presenting with higher tumor stage decreased after radiation exposure, underscoring the potential for this molecule as a signal for determining response to radiation treatment.

Conclusion:

These hypothesis-generating experiments show differential serum protein levels in prostate cancer patient cohorts that have been stratified by risk factors and in longitudinal studies of patients following treatment with SBRT. Candidate biomarker proteins and molecular pathways have been identified for validation as potential predictors of prostate cancer response to radiation treatment.

Citation Format: Einsley Janowski, Simeng Suy, Rency S. Varghese, Olga Timofeeva, Stuart G. Field, Rachel Ostroff, Steve Williams, Anatoly Dritschilo, Sean P. Collins. Serum biomarkers in patients treated with stereotactic body radiation therapy (SBRT) for prostate cancer. [abstract]. In: Proceedings of the 107th Annual Meeting of the American Association for Cancer Research; 2016 Apr 16-20; New Orleans, LA. Philadelphia (PA): AACR; Cancer Res 2016;76(14 Suppl):Abstract nr 450.

Assessment of GC-MS in Detecting Changes in the Levels of Metabolites Using a Spike-in Experiment in Human Plasma

Gas Chromatography coupled with Mass Spectrometry (GC-MS) has been broadly used for the detection of changes in metabolite levels in complex samples. Internal Standards (IS) spiked into a complex background at different concentrations help assess the capability of GC-MS in detecting changes in metabolite levels. This study uses a Latin square design to evaluate the ability of GC-MS in full scan and Single Ion Monitoring (SIM) modes to detect changes among IS spiked into human plasma samples at varying concentrations. Statistical analysis of the data demonstrates the potential of GC-MS to detect true differences over a wide range of concentration levels.

GC-MS Based Plasma Metabolomics for Identification of Candidate Biomarkers for Hepatocellular Carcinoma in Egyptian Cohort

This study evaluates changes in metabolite levels in hepatocellular carcinoma (HCC) cases vs. patients with liver cirrhosis by analysis of human blood plasma using gas chromatography coupled with mass spectrometry (GC-MS). Untargeted metabolomic analysis of plasma samples from participants recruited in Egypt was performed using two GC-MS platforms: a GC coupled to single quadruple mass spectrometer (GC-qMS) and a GC coupled to a time-of-flight mass spectrometer (GC-TOFMS). Analytes that showed statistically significant changes in ion intensities were selected using ANOVA models. These analytes and other candidates selected from related studies were further evaluated by targeted analysis in plasma samples from the same participants as in the untargeted metabolomic analysis. The targeted analysis was performed using the GC-qMS in selected ion monitoring (SIM) mode. The method confirmed significant changes in the levels of glutamic acid, citric acid, lactic acid, valine, isoleucine, leucine, alpha tocopherol, cholesterol, and sorbose in HCC cases vs. patients with liver cirrhosis. Specifically, our findings indicate up-regulation of metabolites involved in branched-chain amino acid (BCAA) metabolism. Although BCAAs are increasingly used as a treatment for cancer cachexia, others have shown that BCAA supplementation caused significant enhancement of tumor growth via activation of mTOR/AKT pathway, which is consistent with our results that BCAAs are up-regulated in HCC.

LC-MS/MS-based serum proteomics for identification of candidate biomarkers for hepatocellular carcinoma

Associating changes in protein levels with the onset of cancer has been widely investigated to identify clinically relevant diagnostic biomarkers. In the present study, we analyzed sera from 205 patients recruited in the United States and Egypt for biomarker discovery using label-free proteomic analysis by LC-MS/MS. We performed untargeted proteomic analysis of sera to identify candidate proteins with statistically significant differences between hepatocellular carcinoma (HCC) and patients with liver cirrhosis. We further evaluated the significance of 101 proteins in sera from the same 205 patients through targeted quantitation by MRM on a triple quadrupole mass spectrometer. This led to the identification of 21 candidate protein biomarkers that were significantly altered in both the United States and Egyptian cohorts. Among the 21 candidates, ten were previously reported as HCC-associated proteins (eight exhibiting consistent trends with our observation), whereas 11 are new candidates discovered by this study. Pathway analysis based on the significant proteins reveals upregulation of the complement and coagulation cascades pathway and downregulation of the antigen processing and presentation pathway in HCC cases versus patients with liver cirrhosis. The results of this study demonstrate the power of combining untargeted and targeted quantitation methods for a comprehensive serum proteomic analysis, to evaluate changes in protein levels and discover novel diagnostic biomarkers. All MS data have been deposited in the ProteomeXchange with identifier PXD001171 (http://proteomecentral.proteomexchange.org/dataset/PXD001171).

Pathway and network approaches for identification of cancer signature markers from omics data

The advancement of high throughput omic technologies during the past few years has made it possible to perform many complex assays in a much shorter time than the traditional approaches. The rapid accumulation and wide availability of omic data generated by these technologies offer great opportunities to unravel disease mechanisms, but also presents significant challenges to extract knowledge from such massive data and to evaluate the findings. To address these challenges, a number of pathway and network based approaches have been introduced. This review article evaluates these methods and discusses their application in cancer biomarker discovery using hepatocellular carcinoma (HCC) as an example.

LC-MS profiling of N-Glycans derived from human serum samples for biomarker discovery in hepatocellular carcinoma

Defining clinically relevant biomarkers for early stage hepatocellular carcinoma (HCC) in a high-risk population of cirrhotic patients has potentially far-reaching implications for disease management and patient health. Changes in glycan levels have been associated with the onset of numerous diseases including cancer. In the present study, we used liquid chromatography coupled with electrospray ionization mass spectrometry (LC–ESI-MS) to analyze N-glycans in sera from 183 participants recruited in Egypt and the U.S. and identified candidate biomarkers that distinguish HCC cases from cirrhotic controls. N-Glycans were released from serum proteins and permethylated prior to the LC–ESI-MS analysis. Through two complementary LC–ESI-MS quantitation approaches, global profiling and targeted quantitation, we identified 11 N-glycans with statistically significant differences between HCC cases and cirrhotic controls. These glycans can further be categorized into four structurally related clusters, matching closely with the implications of important glycosyltransferases in cancer progression and metastasis. The results of this study illustrate the power of the integrative approach combining complementary LC–ESI-MS based quantitation approaches to investigate changes in N-glycan levels between HCC cases and patients with liver cirrhosis.

Evaluation of metabolite biomarkers for hepatocellular carcinoma through stratified analysis by gender, race, and alcoholic cirrhosis

BACKGROUND

The effects of hepatocellular carcinoma on liver metabolism and circulating metabolites have been subjected to continuing investigation. This study compares the levels of selected metabolites in sera of hepatocellular carcinoma cases versus patients with liver cirrhosis and evaluates the influence of gender, race, and alcoholic cirrhosis on the performance of the metabolites as candidate biomarkers for hepatocellular carcinoma.

METHODS

Targeted quantitation of 15 metabolites is performed by selected research monitoring in sera from 89 Egyptian subjects (40 hepatocellular carcinoma cases and 49 cirrhotic controls) and 110 U.S. subjects (56 hepatocellular carcinoma cases and 54 cirrhotic controls). Logistic regression models are used to evaluate the ability of these metabolites in distinguishing hepatocellular carcinoma cases from cirrhotic controls. The influences of gender, race, and alcoholic cirrhosis on the performance of the metabolites are analyzed by stratified logistic regression.

RESULTS

Two metabolites are selected on the basis of their significance to both cohorts. Although both metabolites discriminate hepatocellular carcinoma cases from cirrhotic controls in males and Caucasians, they are insignificant in females and African Americans. One metabolite is significant in patients with alcoholic cirrhosis and the other in nonalcoholic cirrhosis.

CONCLUSIONS

The study demonstrates the potential of two metabolites as candidate biomarkers for hepatocellular carcinoma by combining them with α-fetoprotein (AFP) and gender. Stratified statistical analyses reveal that gender, race, and alcoholic cirrhosis affect the relative levels of small molecules in serum.

IMPACT

The findings of this study contribute to a better understanding of the influence of gender, race, and alcoholic cirrhosis in investigating small molecules as biomarkers for hepatocellular carcinoma.

Functional proteomics analysis to study ATM dependent signaling in response to ionizing radiation

Ataxia telangiectasia (AT) is a human genetic disease characterized by radiation sensitivity, impaired neuronal development and predisposition to cancer. Using a genetically defined model cell system consisting of cells expressing a kinase dead or a kinase proficient ATM gene product, we previously reported systemic alterations in major metabolic pathways that translate at the gene expression, protein and small molecule metabolite levels. Here, we report ionizing radiation induced stress response signaling arising from perturbations in the ATM gene, by employing a functional proteomics approach. Functional pathway analysis shows robust translational and post-translational responses under ATM proficient conditions, which include enrichment of proteins in the Ephrin receptor and axonal guidance signaling pathways. These molecular networks offer a hypothesis generating function for further investigations of cellular stress responses.

LC-MS based serum metabolomics for identification of hepatocellular carcinoma biomarkers in Egyptian cohort

Although hepatocellular carcinoma (HCC) has been subjected to continuous investigation and its symptoms are well-known, early stage diagnosis of this disease remains difficult and the survival rate after diagnosis is typically very low (3-5%). Early and accurate detection of metabolic changes in the sera of patients with liver cirrhosis can help improve the prognosis of HCC and lead to a better understanding of its mechanism at the molecular level, thus providing patients with in-time treatment of the disease. In this study, we compared metabolite levels in sera of 40 HCC patients and 49 cirrhosis patients from Egypt by using ultraperformance liquid chromatography coupled with quadrupole time-of-flight mass spectrometer (UPLC-QTOF MS). Following data preprocessing, the most relevant ions in distinguishing HCC cases from cirrhotic controls are selected by statistical methods. Putative metabolite identifications for these ions are obtained through mass-based database search. The identities of some of the putative identifications are verified by comparing their MS/MS fragmentation patterns and retention times with those from authentic compounds. Finally, the serum samples are reanalyzed for quantitation of selected metabolites as candidate biomarkers of HCC. This quantitation was performed using isotope dilution by selected reaction monitoring (SRM) on a triple quadrupole linear ion trap (QqQLIT) coupled to UPLC. Statistical analysis of the UPLC-QTOF data identified 274 monoisotopic ion masses with statistically significant differences in ion intensities between HCC cases and cirrhotic controls. Putative identifications were obtained for 158 ions by mass based search against databases. We verified the identities of selected putative identifications including glycholic acid (GCA), glycodeoxycholic acid (GDCA), 3β, 6β-dihydroxy-5β-cholan-24-oic acid, oleoyl carnitine, and Phe-Phe. SRM-based quantitation confirmed significant differences between HCC and cirrhotic controls in metabolite levels of bile acid metabolites, long chain carnitines and small peptide. Our study provides useful insight into appropriate experimental design and computational methods for serum biomarker discovery using LC-MS/MS based metabolomics. This study has led to the identification of candidate biomarkers with significant changes in metabolite levels between HCC cases and cirrhotic controls. This is the first MS-based metabolic biomarker discovery study on Egyptian subjects that led to the identification of candidate metabolites that discriminate early stage HCC from patients with liver cirrhosis.

Utilization of metabolomics to identify serum biomarkers for hepatocellular carcinoma in patients with liver cirrhosis

Characterizing the metabolic changes pertaining to hepatocellular carcinoma (HCC) in patients with liver cirrhosis is believed to contribute towards early detection, treatment, and understanding of the molecular mechanisms of HCC. In this study, we compare metabolite levels in sera of 78 HCC cases with 184 cirrhotic controls by using ultra performance liquid chromatography coupled with a hybrid quadrupole time-of-flight mass spectrometry (UPLC-QTOF MS). Following data preprocessing, the most relevant ions in distinguishing HCC cases from patients with cirrhosis are selected by parametric and non-parametric statistical methods. Putative metabolite identifications for these ions are obtained through mass-based database search. Verification of the identities of selected metabolites is conducted by comparing their MS/MS fragmentation patterns and retention time with those from authentic compounds. Quantitation of these metabolites is performed in a subset of the serum samples (10 HCC and 10 cirrhosis) using isotope dilution by selected reaction monitoring (SRM) on triple quadrupole linear ion trap (QqQLIT) and triple quadrupole (QqQ) mass spectrometers. The results of this analysis confirm that metabolites involved in sphingolipid metabolism and phospholipid catabolism such as sphingosine-1-phosphate (S-1-P) and lysophosphatidylcholine (lysoPC 17:0) are up-regulated in sera of HCC vs. those with liver cirrhosis. Down-regulated metabolites include those involved in bile acid biosynthesis (specifically cholesterol metabolism) such as glycochenodeoxycholic acid 3-sulfate (3-sulfo-GCDCA), glycocholic acid (GCA), glycodeoxycholic acid (GDCA), taurocholic acid (TCA), and taurochenodeoxycholate (TCDCA). These results provide useful insights into HCC biomarker discovery utilizing metabolomics as an efficient and cost-effective platform. Our work shows that metabolomic profiling is a promising tool to identify candidate metabolic biomarkers for early detection of HCC cases in high risk population of cirrhotic patients.

Ion annotation-assisted analysis of LC-MS based metabolomic experiment

Background

Analysis of multiple LC-MS based metabolomic studies is carried out to determine overlaps and differences among various experiments. For example, in large metabolic biomarker discovery studies involving hundreds of samples, it may be necessary to conduct multiple experiments, each involving a subset of the samples due to technical limitations. The ions selected from each experiment are analyzed to determine overlapping ions. One of the challenges in comparing the ion lists is the presence of a large number of derivative ions such as isotopes, adducts, and fragments. These derivative ions and the retention time drifts need to be taken into account during comparison.

Results

We implemented an ion annotation-assisted method to determine overlapping ions in the presence of derivative ions. Following this, each ion is represented by the monoisotopic mass of its cluster. This mass is then used to determine overlaps among the ions selected across multiple experiments.

Conclusion

The resulting ion list provides better coverage and more accurate identification of metabolites compared to the traditional method in which overlapping ions are selected on the basis of individual ion mass.

Abstract 4795: Metabolic biomarkers of hepatocellular carcinoma

Characterizing the metabolic changes pertaining to hepatocellular carcinoma (HCC) in patients with liver cirrhosis (LC) is believed to contribute towards early detection, treatment, and our understanding of the molecular mechanisms of HCC. Prospectively evaluated adult patients with LC were recruited from our hepatology practice. All patients were diagnosed to have LC on the basis of established clinical, laboratory and imaging criteria. Cases were diagnosed to have HCC based on either histology or characteristic imaging features. The LC controls were required to be HCC free for at least 6 months from the time of study entry. We compared metabolite levels in sera of 78 HCC cases with 184 LC controls by using ultra performance liquid chromatography coupled with a hybrid quadrupole time-of-flight mass spectrometry (UPLC-QTOF MS). Following data preprocessing, the most relevant ions in distinguishing HCC cases from LC controls were selected by statistical methods. Putative metabolite identifications for these ions were obtained through mass-based database search. Pathway analysis of these putative identifications reveals significant enrichment of metabolites involved in bile acid biosynthesis, porphyrin and chlorophyll metabolism, and arachidonic acid metabolism. Verification of the identities of selected metabolites was conducted by comparing their MS/MS fragmentation patterns and retention time with those from authentic compounds. Quantitation of these metabolites was performed in a subset of the serum samples (10 HCC and 10 LC) using isotope dilution by selected reaction monitoring (SRM) on triple quadrupole linear ion trap (QqQLIT) mass spectrometer. The results of this analysis confirms that metabolites involved in phospholipid catabolism such as sphingosine-1-phosphate (S-1-P) and lysophosphatidylcholine are up-regulated in sera of HCC vs. those with LC. Down-regulated metabolites include those involved in cholesterol metabolism such as glycochenodeoxycholic acid 3-sulfate (3-sulfo-GCDCA), glycocholic acid (GCA), glycodeoxycholic acid (GDCA), taurocholic acid (TCA), taurochenodeoxycholate (TCDCA), as well as those involved in heme catabolism including billuribin, I-Urobili, and bilirubin glucuronide. The preliminary results obtained in this study provide useful insights into our understanding of metabolic differences between HCC and LC.

Citation Format: {Authors}. {Abstract title} [abstract]. In: Proceedings of the 103rd Annual Meeting of the American Association for Cancer Research; 2012 Mar 31-Apr 4; Chicago, IL. Philadelphia (PA): AACR; Cancer Res 2012;72(8 Suppl):Abstract nr 4795. doi:1538-7445.AM2012-4795

LC-MS data analysis for differential protein expression detection

In proteomic studies, liquid chromatography coupled with mass spectrometry (LC-MS) is a common platform to compare the abundance of various peptides that characterize particular proteins in biological samples. Each LC-MS run generates data consisting of thousands of peak intensities for peptides represented by retention time (RT) and mass-to-charge ratio (m/z) values. In label-free differential protein expression studies, multiple LC-MS runs are compared to identify differentially abundant peptides between distinct biological groups. This approach presents a computational challenge because of the following reasons (i) substantial variation in RT across multiple runs due to the LC instrument conditions and the variable complexity of peptide mixtures, (ii) variation in m/z values due to occasional drift in the calibration of the mass spectrometry instrument, and (iii) variation in peak intensities caused by various factors including noise and variability in sample handling and processing. In this chapter, we present computational methods for quantification and comparison of peptides by label-free LC-MS analysis. We discuss data preprocessing methods for alignment and normalization of LC-MS data. Also, we present multivariate statistical methods and pattern recognition methods for detection of differential protein expression from preprocessed LC-MS data.

Abstract 5569: N-glycosylation in cultured cells and hepatocellular carcinoma

Purpose: Incidence of hepatocellular cancer (HCC) continues to increase primarily due to hepatitis C viral (HCV) infection. Prognosis and survival of patients is highly affected by the disease stage at the time of diagnosis. Our mass spectrometric study evaluated N-glycans during the progression of HCV infection to cancer. Glycosylation of immunoglobulins and other serum proteins was examined in cultured cells and in serum of patients with hepatocellular carcinoma. Methods and Results: Protein associated N-glycans were released with PNGaseF and analyzed by MALDI-TOF/TOF analysis following solid phase permethylation. Analysis of less than 0.04 ml of serum led to relative quantification of 70 N-glycan structures. Immunoaffinity isolation of immunoglobulins and other serum proteins allowed us to study protein-specific glycosylation. Progression of HCV infection to HCC was strongly associated with changes in the glycosylation of immunoglobulins. In a pilot case-control study (25 HCC cases and 35 controls), N-glycosylation of eight glycans was significantly different in HCC compared to chronic liver disease controls. The glycosylation of serum proteins was further compared to the proteins secreted by the Huh-7.5 cell line. Conclusion: This study demonstrates mass spectrometric analysis of 70 N-glycans in serum, cultured cells, and isolated proteins. Our results show that the analysis of permethylated N-glycans helps to define changes associated with the progression of HCV infection to HCC. Evaluation of glycan abundance suggests the potential to use glycans for the detection of liver disease.

Citation Format: {Authors}. {Abstract title} [abstract]. In: Proceedings of the 101st Annual Meeting of the American Association for Cancer Research; 2010 Apr 17-21; Washington, DC. Philadelphia (PA): AACR; Cancer Res 2010;70(8 Suppl):Abstract nr 5569.

Identification of N-glycan serum markers associated with hepatocellular carcinoma from mass spectrometry data

Glycocylation represents the most complex and widespread post-translational modifications in human proteins. The variation of glycosylation is closely related to oncogenic transformation. Therefore, profiling of glycans detached from proteins is a promising strategy to identify biomarkers for cancer detection. This study identified candidate glycan biomarkers associated with hepatocellular carcinoma by mass spectrometry. Specifically, mass spectrometry data were analyzed with a peak selection procedure which incorporates multiple random sampling strategies with recursive feature selection based on support vector machines. Ten peak sets were obtained from different combinations of samples. Seven peaks were shared by each of the 10 peaksets, in which 7-12 peaks were selected, indicating 58-100% of peaks were shared by the 10 peaksets. Support vector machines and hierarchical clustering method were used to evaluate the performance of the peaksets. The predictive performance of the seven peaks was further evaluated by using 19 newly generated MALDI-TOF spectra. Glycan structures for four glycans of the seven peaks were determined. Literature search indicated that the structures of the four glycans could be found in some cancer-related glycoproteins. The method of this study is significant in deriving consistent, accurate, and biological significant glycan marker candidates for hepatocellular carcinoma diagnosis.

Enhanced expression of SOS1 is detected in prostate cancer epithelial cells from African-American men

African-American (AA) men experience an increased risk of developing prostate cancers as well as increased mortality following treatment as compared to European-American (EA) men. The aim of our study was to identify biological factors with the potential to predispose AA men to prostate tumor progression and metastasis. To identify cancer-specific gene expression patterns in AA men, we established primary prostate cancer epithelial cells from 14 AA and 13 EA men. High-throughput microarrays were used to investigate differences in global gene expression comparing the two groups. Quantitative RT-PCR and immunohistochemistry validated mRNA and protein expression levels. RNAi knockdowns provided support for biological significance for the identified genes in prostate cancer cells. Son of sevenless homolog 1 (SOS1) was overexpressed in AA male-derived primary prostate cancer epithelial cells. Depletion of SOS1 in PC3 and DU145 prostate cancer cells resulted in decreased capacities for cell proliferation, migration and invasion, at least partially through inhibition of extracellular signal-regulated kinase 1 and 2. Tissue microarray analyses of SOS1 expression in prostate carcinomas correlated with Gleason's grades of tumors, consistent with a possible role in prostate cancer progression. Investigation of prostate cancer-derived epithelial cells has led to identification of SOS1 as a potential candidate biomarker and molecular therapeutic target in prostate cancer in AA men, consistent with the hypothesis that a biological basis exists for prostate cancer aggressiveness in AA men.

Integrated peptide and glycan biomarker discovery using MALDI-TOF mass spectrometry

Quantitative comparison of peptides and glycans in serum is conducted using matrix-assisted laser desorption/ionization-time of flight mass spectrometry (MALDI-TOF MS) to identify biomarkers. A peak selection algorithm is developed to identify a panel of integrated peptide and glycan peaks to distinguish hepatocellular carcinoma (HCC) cases from high-risk population of patients with chronic liver disease (CLD). Candidate peptide and glycan markers selected frequently in multiple runs of the algorithm are presented. The performance of these markers is evaluated in terms of their ability to distinguish HCC cases from patients with CLD in a blinded validation set.

Detection of hepatocellular carcinoma using glycomic analysis

PURPOSE

Hepatocellular carcinoma (HCC) represents an increasing health problem in the United States. Serum alpha-fetoprotein, the currently used clinical marker, is elevated in only approximately 60% of HCC patients; therefore, the identification of additional markers is expected to have significant public health impact. The objective of our study was to quantitatively assess N-glycans originating from serum glycoproteins as alternative markers for the detection of HCC.

EXPERIMENTAL DESIGN

We used matrix-assisted laser desorption/ionization time-of-flight mass spectrometry for quantitative comparison of 83 N-glycans in serum samples of 202 participants (73 HCC cases, 77 age- and gender-matched cancer-free controls, and 52 patients with chronic liver disease). N-glycans were enzymatically released from serum glycoproteins and permethylated before mass spectrometric quantification.

RESULTS

The abundance of 57 N-glycans was significantly altered in HCC patients compared with controls. The sensitivity of six individual glycans evaluated for separation of HCC cases from population controls ranged from 73% to 90%, and the specificity ranged from 36% to 91%. A combination of three selected N-glycans was sufficient to classify HCC with 90% sensitivity and 89% specificity in an independent validation set of patients with chronic liver disease. The three N-glycans remained associated with HCC after adjustment for chronic viral infection and other known covariates, whereas the other glycans increased significantly at earlier stages of the progression of chronic viral infection to HCC.

CONCLUSION

A set of three identified N-glycans is sufficient for the detection of HCC with 90% prediction accuracy in a population with high rates of hepatitis C viral infection. Further evaluation of a wider clinical utility of these candidate markers is warranted.

Ant colony optimization for biomarker identification from MALDI-TOF mass spectra

We present a novel method that combines ant colony optimization with support vector machines (ACO-SVM) to select candidate biomarkers from MALDI-TOF serum profiles of hepatocellular carcinoma (HCC) patients and matched controls. The method identified relevant mass points that achieve high sensitivity and specificity in distinguishing HCC patients from healthy individuals. The results indicate that the MALDI-TOF technology could provide the means to discover novel biomarkers for HCC.

Analysis of MALDI-TOF mass spectrometry data for discovery of peptide and glycan biomarkers of hepatocellular carcinoma

This paper presents computational methods to analyze MALDI-TOF mass spectrometry data for quantitative comparison of peptides and glycans in serum. The methods are applied to identify candidate biomarkers in serum samples of 203 participants from Egypt; 73 hepatocellular carcinoma (HCC) cases, 52 patients with chronic liver disease (CLD) consisting of cirrhosis and fibrosis cases, and 78 population controls. Two complementary sample preparation methods were applied prior to generating mass spectra: (1) low molecular weight (LMW) enrichment of each serum sample was carried out for MALDI-TOF quantification of peptides, and (2) glycans were enzymatically released from proteins in each serum sample and permethylated for MALDI-TOF quantification of glycans. A peak selection algorithm was applied to identify the most useful peptide and glycan peaks for accurate detection of HCC cases from high-risk population of patients with CLD. In addition to global peaks selected by the whole population based approach, where identically labeled patients are treated as a single group, subgroup-specific peaks were identified by searching for peaks that are differentially abundant in a subgroup of patients only. The peak selection process was preceded by peak screening, where we eliminated peaks that have significant association with covariates such as age, gender, and viral infection based on the peptide and glycan spectra from population controls. The performance of the selected peptide and glycan peaks was evaluated in terms of their ability in detecting HCC cases from patients with CLD in a blinded validation set and through the cross-validation method. Finally, we investigated the possibility of using both peptides and glycans in a panel to enhance the diagnostic capability of these candidate markers. Further evaluation is needed to examine the potential clinical utility of the candidate peptide and glycan markers identified in this study.

Classification algorithms for phenotype prediction in genomics and proteomics

This paper gives an overview of statistical and machine learning-based feature selection and pattern classification algorithms and their application in molecular cancer classification or phenotype prediction. In particular, the paper focuses on the use of these computational methods for gene and peak selection from microarray and mass spectrometry data, respectively. The selected features are presented to a classifier for phenotype prediction.

Analysis of MALDI-TOF mass spectrometry data for detection of glycan biomarkers

We present a computational framework for analysis of MALDI-TOF mass spectrometry data to enable quantitative comparison of glycans in serum. The proposed framework enables a systematic selection of glycan structures that have good generalization capability in distinguishing subjects from two pre-labeled groups. We applied the proposed method for a biomarker discovery study that involves 203 participants from Cairo, Egypt; 73 hepatocellular carcinoma (HCC) cases, 52 patients with chronic liver disease (CLD), and 78 healthy individuals. Glycans were enzymatically released from proteins in serum and permethylated prior to mass spectrometric quantification. A subset of the participants (35 HCC and 35 CLD cases) was used as a training set to select global and subgroup-specific peaks. The peak selection step is preceded by peak screening, where we eliminate peaks that seem to have association with covariates such as age, gender, and viral infection based on the 78 spectra from healthy individuals. To ensure that the global peaks have good generalization capability, we subjected the entire spectral preprocessing and peak selection step to a cross-validation; a randomly selected subset of the training set was used for spectral preprocessing and peak selection in multiple runs with resubstitution. In addition to global peak identification method, we describe a new approach that allows the selection of subgroup-specific glycans by searching for glycans that display differential abundance in a subgroup of patients only. The performance of the global and subgroup-specific peaks is evaluated via a blinded independent set that comprises of 38 HCC and 17 CLD cases. Further evaluation of the potential clinical utility of the selected global and subgroup-specific candidate markers is needed.

Candidate markers for the detection of hepatocellular carcinoma in low-molecular weight fraction of serum

Hepatocellular carcinoma (HCC) represents an important public health problem in Egypt where up to 90% of HCC cases are attributable to hepatitis C viral (HCV) infection. Serum alpha-fetoprotein is elevated in only approximately 60% of HCC patients. The development of effective markers for the detection of HCC could have an impact on cancer mortality and significant public health implications worldwide. The objective of our study was to assess six candidate markers for detection of HCC identified by mass spectrometric analysis of enriched serum. The study examined 78 HCC cases and 72 age- and gender-matched cancer-free controls recruited from the Egyptian population. Matrix-assisted laser desorption-ionization time-of-flight mass spectrometric analysis of enriched low-molecular weight fraction of serum was used for identification of the candidate markers. Our analyses show that all six candidate markers are associated with HCC after adjustment for important covariates including HCV and hepatitis B viral infections. The marker candidates are independently predictive of HCC with areas under the receiver operating characteristic (AuROC) curve ranging from 63-93%. A combination of the six markers improves prediction accuracy to 100% sensitivity, 91% specificity and 98% AuROC curve in an independent test set of 50 patients. Two of the candidate markers were identified by sequencing as fragments of complement C3 and C4. In conclusion, a set of six peptides distinguished with high prediction accuracy HCC from controls in an Egyptian population with a high rate of chronic HCV infection. Further evaluation of these marker candidates for the diagnosis of HCC is needed.

Peak selection from MALDI-TOF mass spectra using ant colony optimization

MOTIVATION

Due to the large number of peaks in mass spectra of low-molecular-weight (LMW) enriched sera, a systematic method is needed to select a parsimonious set of peaks to facilitate biomarker identification. We present computational methods for matrix-assisted laser desorption/ionization time-of-flight (MALDI-TOF) spectral data preprocessing and peak selection. In particular, we propose a novel method that combines ant colony optimization (ACO) with support vector machines (SVM) to select a small set of useful peaks.

RESULTS

The proposed hybrid ACO-SVM algorithm selected a panel of eight peaks out of 228 candidate peaks from MALDI-TOF spectra of LMW enriched sera. An SVM classifier built with these peaks achieved 94% sensitivity and 100% specificity in distinguishing hepatocellular carcinoma from cirrhosis in a blind validation set of 69 samples. Area under the receiver operating characteristic (ROC) curve was 0.996. The classification capability of these peaks is compared with those selected by the SVM-recursive feature elimination method.

AVAILABILITY

Supplementary material and MATLAB scripts to implement the methods described in this article are available at http://microarray.georgetown.edu/web/files/bioinf.htm.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Enrichment of low molecular weight fraction of serum for MS analysis of peptides associated with hepatocellular carcinoma

A challenging aspect of biomarker discovery in serum is the interference of abundant proteins with identification of disease-related proteins and peptides. This study describes enrichment of serum by denaturing ultrafiltration, which enables an efficient profiling and identification of peptides up to 5 kDa. We consistently detect several hundred peptide-peaks in MALDI-TOF and SELDI-TOF spectra of enriched serum. The sample preparation is fast and reproducible with an average CV for all 276 peaks in the MALDI-TOF spectrum of 11%. Compared to unenriched serum, the number of peaks in enriched spectra is 4 times higher at an S/N ratio of 5 and 20 times higher at an S/N ratio of 10. To demonstrate utility of the methods, we compared 20 enriched sera of patients with hepatocellular carcinoma (HCC) and 20 age-matched controls using MALDI-TOF. The comparison of 332 peaks at p < 0.001 identified 45 differentially abundant peaks that classified HCC with 90% accuracy in this small pilot study. Direct TOF/TOF sequencing of the most abundant peptide matches with high probability des-Ala-fibrinopeptide A. This study shows that enrichment of the low molecular weight fraction of serum facilitates an efficient discovery of peptides that could serve as biomarkers for detection of HCC as well as other diseases.

Analysis of mass spectral serum profiles for biomarker selection

MOTIVATION

Mass spectrometric profiles of peptides and proteins obtained by current technologies are characterized by complex spectra, high dimensionality and substantial noise. These characteristics generate challenges in the discovery of proteins and protein-profiles that distinguish disease states, e.g. cancer patients from healthy individuals. We present low-level methods for the processing of mass spectral data and a machine learning method that combines support vector machines, with particle swarm optimization for biomarker selection.

RESULTS

The proposed method identified mass points that achieved high prediction accuracy in distinguishing liver cancer patients from healthy individuals in SELDI-QqTOF profiles of serum.

AVAILABILITY

MATLAB scripts to implement the methods described in this paper are available from the HWR's lab website http://lombardi.georgetown.edu/labpage

Increasing the efficiency of fuzzy logic-based gene expression data analysis

DNA microarray technology can accommodate a multifaceted analysis of the expression of genes in an organism. The wealth of spatiotemporal data generated by this technology allows researchers to potentially reverse engineer a particular genetic network. "Fuzzy logic" has been proposed as a method to analyze the relationships between genes and help decipher a genetic network. This method can identify interacting genes that fit a known "fuzzy" model of gene interaction by testing all combinations of gene expression profiles. This paper introduces improvements made over previous fuzzy gene regulatory models in terms of computation time and robustness to noise. Improvement in computation time is achieved by using a cluster analysis as a preprocessing method to reduce the total number of gene combinations analyzed. This approach speeds up the algorithm by a factor of 50% with minimal effect on the results. The model's sensitivity to noise is reduced by implementing appropriate methods of "fuzzy rule aggregation" and "conjunction" that produce reliable results in the face of minor changes in model input.