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.

Inhibition of discoidin domain receptor (DDR)-1 with nilotinib alters CSF miRNAs and is associated with reduced inflammation and vascular fibrosis in Alzheimer's disease.

Discoidin Domain Receptor (DDR)-1 is activated by collagen. Nilotinib is a tyrosine kinase inhibitor that is FDA-approved for leukemia and potently inhibits DDR-1. Individuals diagnosed with mild-moderate Alzheimer's disease (AD) treated with nilotinib (versus placebo) for 12 months showed reduction of amyloid plaque and cerebrospinal fluid (CSF) amyloid, and attenuation of hippocampal volume loss. However, the mechanisms are unclear. Here, we explored unbiased next generation whole genome miRNA sequencing from AD patients CSF and miRNAs were matched with their corresponding mRNAs using gene ontology. Changes in CSF miRNAs were confirmed via measurement of CSF DDR1 activity and plasma levels of AD biomarkers. Approximately 1050 miRNAs are detected in the CSF but only 17 miRNAs are specifically altered between baseline and 12-month treatment with nilotinib versus placebo. Treatment with nilotinib significantly reduces collagen and DDR1 gene expression (upregulated in AD brain), in association with inhibition of CSF DDR1. Pro-inflammatory cytokines, including interleukins and chemokines are reduced along with caspase-3 gene expression. Specific genes that indicate vascular fibrosis, e.g., collagen, Transforming Growth Factors (TGFs) and Tissue Inhibitors of Metalloproteases (TIMPs) are altered by DDR1 inhibition with nilotinib. Specific changes in vesicular transport, including the neurotransmitters dopamine and acetylcholine, and autophagy genes, including ATGs, indicate facilitation of autophagic flux and cellular trafficking. Inhibition of DDR1 with nilotinib may be a safe and effective adjunct treatment strategy involving an oral drug that enters the CNS and adequately engages its target. DDR1 inhibition with nilotinib exhibits multi-modal effects not only on amyloid and tau clearance but also on anti-inflammatory markers that may reduce cerebrovascular fibrosis.

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.

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.

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.

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).

Hepatitis C virus Genotype 1a core gene nucleotide patterns associated with hepatocellular carcinoma risk.

Specific sequence changes in codons 70 and 91 of the hepatitis C virus genotype 1b (HCV GT1b) core gene have been associated with increased risk of hepatocellular carcinoma (HCC). Essentially all previous studies were conducted in Asian populations with a wide range of liver disease, and none were conducted specifically in GT1a-infected individuals. We conducted a pilot study in a multiethnic population in the USA with HCV-related cirrhosis to determine if this association extended to GT1a-infected individuals and to determine if other sequence changes in the HCV core gene were associated with HCC risk. HCV core gene sequences from sera of 90 GT1 HCV carriers with cirrhosis (42 with HCC) were analysed using standard RT-PCR-based procedures. Nucleotide sequence data were compared with reference sequences available from GenBank. The frequency of sequence changes in codon 91 was not statistically different between HCC (7/19) and non-HCC (11/22) GT1b carriers. In GT1a carriers, sequence changes in codon 91 were observed less often than in GT1b carriers but were not observed in non-HCC subjects (4/23 vs 0/26, P = 0.03, Fisher's exact test). Sequence changes in codon 70 were not distributed differently between HCC and non-HCC GT1a and 1b carriers. Most importantly, for GT1a carriers, a panel of specific nucleotide changes in other codons was collectively present in all subjects with HCC, but not in any of the non-HCC patients. The utility of this test panel for early detection of HCC in GT1a-infected individuals needs to be assessed in larger populations, including longitudinal studies.

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.

RNA-Seq Data Analysis.

RNA-Seq data analysis stands as a vital part of genomics research, turning vast and complex datasets into meaningful biological insights. It is a field marked by rapid evolution and ongoing innovation, necessitating a thorough understanding for anyone seeking to unlock the potential of RNA-Seq data. In this chapter, we describe the intricate landscape of RNA-seq data analysis, elucidating a comprehensive pipeline that navigates through the entirety of this complex process. Beginning with quality control, the chapter underscores the paramount importance of ensuring the integrity of RNA-seq data, as it lays the groundwork for subsequent analyses. Preprocessing is then addressed, where the raw sequence data undergoes necessary modifications and enhancements, setting the stage for the alignment phase. This phase involves mapping the processed sequences to a reference genome, a step pivotal for decoding the origins and functions of these sequences.Venturing into the heart of RNA-seq analysis, the chapter then explores differential expression analysis-the process of identifying genes that exhibit varying expression levels across different conditions or sample groups. Recognizing the biological context of these differentially expressed genes is pivotal; hence, the chapter transitions into functional analysis. Here, methods and tools like Gene Ontology and pathway analyses help contextualize the roles and interactions of the identified genes within broader biological frameworks. However, the chapter does not stop at conventional analysis methods. Embracing the evolving paradigms of data science, it delves into machine learning applications for RNA-seq data, introducing advanced techniques in dimension reduction and both unsupervised and supervised learning. These approaches allow for patterns and relationships to be discerned in the data that might be imperceptible through traditional methods.

Short-Read RNA-Seq.

RNA sequencing (RNA-Seq) has emerged as a powerful and versatile tool for the comprehensive analysis of transcriptomes and has been widely used to investigate gene expression, copy number variation, alternative splicing, and novel transcript discovery. This chapter outlines the methodology for conducting short-read RNA-Seq, starting from RNA enrichment to library preparation and sequencing. Throughout the chapter, practical tips and best practices are provided to guide researchers in order to optimize each step of the RNA-Seq workflow. Multiple quality control steps throughout the workflow that are critical to obtain high-quality RNA-Seq data are also discussed.

Non-fatal outcomes of COVID-19 disease in pediatric organ transplantation associates with down-regulation of senescence pathways.

This is a cross-sectional study examining kinetics and durability of immune response in children with solid organ transplants (SOTs) who had COVID-19 disease between November 2020 through June 2022, who were followed for 60-days at a single transplant center. Blood was collected between 1-14 (acute infection), and 15-60 days of a positive PCR (convalescence). SOT children with peripheral blood mononuclear cells (PBMC) cryopreserved before 2019 were non-infected controls (ctrls). PBMCs stimulated with 15-mer peptides from spike protein and anti-CD49d/anti-CD28. Testing done included mass cytometry, mi-RNA sequencing with confirmatory qPCR. 38 children formed the study cohort, 10 in the acute phase and 8 in the convalescence phase. 20 subjects were non-infected controls. Two subjects had severe disease. Subjects in the acute and convalescent phases were different subjects. The median age and tacrolimus level at blood draw was not significantly different. There was no death, and no subject was lost to follow-up. During acute infection CD57 expression was low in NKT, Th17 effector memory, memory Treg, CD4-CD8-, and γδT cells (p = 0.01, p = 0.04, p = 0.03, p = 0.03, p = 0.004 respectively). The frequencies of NK and Th2 effector memory cells increased (p = 0.01, p = 0.02) during acute infection. Non-switched memory B and CD8 central memory cell frequencies were decreased during acute infection (p = 0.02; p = 0.02), but the decrease in CD8 central memory cells did not persist. CD4-CD8- and CD14 monocyte frequencies increased during recovery (p = 0.03; p = 0.007). Our observations suggest down regulation of CD57 with absence of NK cell contraction protect against death from COVID-19 disease in children with SOTs.

Suppressor T helper type 17 cell responses in intestinal transplant recipients with allograft rejection.

BACKGROUND: Intestinal transplant (ITx) rejection is associated with memory T helper type 17 cell (Th17) infiltration of grafted tissues. Modulation of Th17 effector cell response is facilitated by T regulatory (Treg) cells, but a phenotypic characterization of this process is lacking in the context of allograft rejection. METHODS: Flow cytometry was performed to examine the expression of surface receptors, cytokines, and transcription factors in Th17 and Treg cells in ITx control (n = 34) and rejection patients (n = 23). To elucidate key pathways guiding the rejection biology, we utilized RNA sequencing (RNAseq) and assessed epigenetic stability through pyrosequencing of the Treg-specific demethylated region (TSDR). RESULTS: We found that intestinal allograft rejection is characterized by Treg cellular infiltrates, which are polarized toward Th17-type chemokine receptor, ROR-γt transcription factor expression, and cytokine production. These Treg cell subsets have maintained epigenetic stability, as defined by FoxP3-TSDR methylation status, but displayed upregulation of functional Treg and purinergic signaling genes by RNAseq analysis such as CD39, in keeping with suppressor Th17 properties. CONCLUSION: We show that ITx rejection is associated with increased polarized cells that express a Th17-like phenotype concurrent with regulatory purinergic markers.

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.

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.

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.

MicroRNA-483-5p Inhibits Hepatocellular Carcinoma Cell Proliferation, Cell Steatosis, and Fibrosis by Targeting PPARα and TIMP2.

MicroRNAs (miRNAs) are small non-coding RNA molecules that bind with the 3' untranslated regions (UTRs) of genes to regulate expression. Downregulation of miR-483-5p (miR-483) is associated with the progression of hepatocellular carcinoma (HCC). However, the significant roles of miR-483 in nonalcoholic fatty liver disease (NAFLD), alcoholic fatty liver diseases (AFLD), and HCC remain elusive. In the current study, we investigated the biological significance of miR-483 in NAFLD, AFLD, and HCC in vitro and in vivo. The downregulation of miR-483 expression in HCC patients' tumor samples was associated with Notch 3 upregulation. Overexpression of miR-483 in a human bipotent progenitor liver cell line HepaRG and HCC cells dysregulated Notch signaling, inhibited cell proliferation/migration, induced apoptosis, and increased sensitivity towards antineoplastic agents sorafenib/regorafenib. Interestingly, the inactivation of miR-483 upregulated cell steatosis and fibrosis signaling by modulation of lipogenic and fibrosis gene expression. Mechanistically, miR-483 targets PPARα and TIMP2 gene expression, which leads to the suppression of cell steatosis and fibrosis. The downregulation of miR-483 was observed in mice liver fed with a high-fat diet (HFD) or a standard Lieber-Decarli liquid diet containing 5% alcohol, leading to increased hepatic steatosis/fibrosis. Our data suggest that miR-483 inhibits cell steatosis and fibrogenic signaling and functions as a tumor suppressor in HCC. Therefore, miR-483 may be a novel therapeutic target for NAFLD/AFLD/HCC management in patients with fatty liver diseases and HCC.

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.

Using a spike-in experiment to evaluate analysis of LC-MS data.

BACKGROUND: Recent advances in liquid chromatography-mass spectrometry (LC-MS) technology have led to more effective approaches for measuring changes in peptide/protein abundances in biological samples. Label-free LC-MS methods have been used for extraction of quantitative information and for detection of differentially abundant peptides/proteins. However, difference detection by analysis of data derived from label-free LC-MS methods requires various preprocessing steps including filtering, baseline correction, peak detection, alignment, and normalization. Although several specialized tools have been developed to analyze LC-MS data, determining the most appropriate computational pipeline remains challenging partly due to lack of established gold standards. RESULTS: The work in this paper is an initial study to develop a simple model with "presence" or "absence" condition using spike-in experiments and to be able to identify these "true differences" using available software tools. In addition to the preprocessing pipelines, choosing appropriate statistical tests and determining critical values are important. We observe that individual statistical tests could lead to different results due to different assumptions and employed metrics. It is therefore preferable to incorporate several statistical tests for either exploration or confirmation purpose. CONCLUSIONS: The LC-MS data from our spike-in experiment can be used for developing and optimizing LC-MS data preprocessing algorithms and to evaluate workflows implemented in existing software tools. Our current work is a stepping stone towards optimizing LC-MS data acquisition and testing the accuracy and validity of computational tools for difference detection in future studies that will be focused on spiking peptides of diverse physicochemical properties in different concentrations to better represent biomarker discovery of differentially abundant peptides/proteins.

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.

Convolutional Neural Network-Based Compound Fingerprint Prediction for Metabolite Annotation.

Metabolite annotation has been a challenging issue especially in untargeted metabolomics studies by liquid chromatography coupled with mass spectrometry (LC-MS). This is in part due to the limitations of publicly available spectral libraries, which consist of tandem mass spectrometry (MS/MS) data acquired from just a fraction of known metabolites. Machine learning provides the opportunity to predict molecular fingerprints based on MS/MS data. The predicted molecular fingerprints can then be used to help rank putative metabolite IDs obtained by using either the precursor mass or the formula of the unknown metabolite. This method is particularly useful to help annotate metabolites whose corresponding MS/MS spectra are missing or cannot be matched with those in accessible spectral libraries. We investigated a convolutional neural network (CNN) for molecular fingerprint prediction based on data acquired by MS/MS. We used more than 680,000 MS/MS spectra obtained from the MoNA repository and NIST 20, representing about 36,000 compounds for training and testing our CNN model. The trained CNN model is implemented as a python package, MetFID. The package is available on GitHub for users to enter their MS/MS spectra and corresponding putative metabolite IDs to obtain ranked lists of metabolites. Better performance is achieved by MetFID in ranking putative metabolite IDs using the CASMI 2016 benchmark dataset compared to two other machine learning-based tools (CSI:FingerID and ChemDistiller).