Search

Ion trajectory simulation in mass spectrometry systems from injection to detection is technically challenging but very important for better understanding the ion dynamics in instrument development. Here, we present SimELIT (

Integration of genome-scale metabolic networks with reactive transport models (RTMs) is an advanced simulation technique that enables predicting the changes of microbial growth and metabolism in space and time. Despite promising demonstrations in the past, computational inefficiency has been pointed out as a critical issue to overcome because it requires repeated implementation of linear programming (LP) to get flux balance analysis (FBA) solutions in every time step and every spatial grid. To address this challenge, we propose a new simulation method where we train/validate artificial neural networks (ANNs) using FBA solutions randomly sampled and incorporate the resulting reduced-order FBA model (represented as algebraic equations) into RTMs as source/sink terms. We demonstrate the efficiency of our method via a case study ofShewanella oneidensisMR-1 strain. During the aerobic growth on lactate,S. oneidensisproduces metabolic byproducts (such as pyruvate and acetate), which are subsequently consumed as alternative carbon sources when the preferred ones are depleted. Simulating such intricate dynamics posed a considerable challenge, which we overcame by adopting the cybernetic approach that describes metabolic switches as the outcome of dynamic competition among multiple growth options. In both zero-dimensional batch and one-dimensional column configurations, the ANN-based reduced-order models achieved substantial reduction of computational time by several orders of magnitude compared to the original LP-based FBA models. Importantly, the ANN models produced robust solutions without any special measures to prevent numerical instability. These developments significantly promote our ability to utilize genome-scale networks in complex, multi-physics, and multi-dimensional ecosystem modeling.

Signal digitization is a commonly overlooked part of ion mobility-mass spectrometry (IMS-MS) workflows, yet it greatly affects signal-to-noise ratio and MS resolution measurements. Here, we report on the integration of a 2 GS/s, 14-bit ADC with structures for lossless ion manipulations (SLIM-IMS-MS) and compare the performance to a commonly used 8-bit ADC. The 14-bit ADC provided a reduction in the digitized noise by a factor of ∼6, owing largely to the use of smaller bit sizes. The low baseline allowed threshold voltage levels to be set very close to the MCP baseline voltage, allowing for as much signal to be acquired as possible without overloading or excessive digitization of MCP baseline noise. Analyses of Agilent tuning mixture ions and a mixture of heavy labeled phosphopeptides showed that the 14-bit ADC provided a ∼1.5-2× signal-to-noise (S/N) increase for high intensity ions, such as the Agilent tuning mixture ions and the 2+ and 3+ charge states of many phosphopeptide constituents. However, signal enhancements were as much as 10-fold for low intensity ions, and the 14-bit ADC enabled discernible signal intensities otherwise lost using an 8-bit digitizer. Additionally, the 14-bit ADC required ∼14-fold fewer mass spectra to be averaged to produce a mass spectrum with a similar S/N as the 8-bit ADC, demonstrating ∼10× higher measurement throughput. The high resolution, low baseline, and fast speed of the new 14-bit ADC enables high performance digitization of MS, IMS-MS, and SLIM-IMS-MS spectra and provides a much better picture of analyte profiles in complex mixtures.

Modern mass spectrometry-based workflows employing hybrid instrumentation and orthogonal separations collect multidimensional data, potentially allowing deeper understanding in omics studies through adoption of artificial intelligence methods. However, the large volume of these rich spectra challenges existing data storage and access technologies, therefore precluding informatics advancements. We present MZA (pronounced m-za), the mass-to-charge (m/z) generic data storage and access tool designed to facilitate software development and artificial intelligence research in multidimensional mass spectrometry measurements. Composed by a data conversion tool and a simple file structure based on the HDF5 format, MZA provides easy, cross-platform and cross-programming language access to raw MS-data, enabling fast development of new tools in data science programming languages such as Python and R. The software executable, example MS-data and example Python and R scripts are freely available at https://github.com/PNNL-m-q/mza

Metaproteomics has been increasingly utilized for high-throughput characterization of proteins in complex environments and has been demonstrated to provide insights into microbial composition and functional roles. However, significant challenges remain in metaproteomic data analysis, including creation of a sample-specific protein sequence database. A well-matched database is a requirement for successful metaproteomics analysis, and the accuracy and sensitivity of PSM identification algorithms suffer when the database is incomplete or contains extraneous sequences. When matched DNA sequencing data of the sample is unavailable or incomplete, creating the proteome database that accurately represents the organisms in the sample is a challenge. Here, we leverage a

Graphical abstract, Microbial communities organize into spatial patterns that are largely governed by interspecies interactions. This phenomenon is an important metric for understanding community functional dynamics, yet the use of spatial patterns for predicting microbial interactions is currently lacking. Here we propose supervised deep learning as a new tool for network inference. An agent-based model was used to simulate the spatiotemporal evolution of two interacting organisms under diverse growth and interaction scenarios, the data of which was subsequently used to train deep neural networks. For small-size domains (100 µm × 100 µm) over which interaction coefficients are assumed to be invariant, we obtained fairly accurate predictions, as indicated by an average R2 value of 0.84. In application to relatively larger domains (450 µm × 450 µm) where interaction coefficients are varying in space, deep learning models correctly predicted spatial distributions of interaction coefficients without any additional training. Lastly, we evaluated our model against real biological data obtained using Pseudomonas fluorescens and Escherichia coli co-cultures treated with polymeric chitin or N-acetylglucosamine, the hydrolysis product of chitin. While P. fluorescens can utilize both substrates for growth, E. coli lacked the ability to degrade chitin. Consistent with our expectations, our model predicted context-dependent interactions across two substrates, i.e., degrader-cheater relationship on chitin polymers and competition on monomers. The combined use of the agent-based model and machine learning algorithm successfully demonstrates how to infer microbial interactions from spatially distributed data, presenting itself as a useful tool for the analysis of more complex microbial community interactions.

Lung cancer is the leading cause of cancer-related deaths in the United States, with estimates of 236,740 new cases and 118,830 deaths in 2022 secondary to the disease. Blood-based liquid biopsies hold promise to reduce morbidity and mortality from lung cancer by enabling early detection to downstage disease at diagnosis, theragnostic identification of patients most likely to be helped or harmed by therapy, monitoring of therapeutic efficacy, and detection of residual disease. PrognomiQ’s multi-omics platform comprehensively profiles proteins, metabolites, lipids, mRNA, and cfDNA in blood samples which can be used for the development of liquid biopsy tests with high sensitivity and specificity for lung cancer. We conducted a case-control study comprising 1031 subjects: 361 subjects with untreated non-small-cell lung cancer (NSCLC) and 670 matched controls which included 340 subjects with salient pulmonary and gastrointestinal co-morbidities. Blood samples from each subject were processed to provide 7 different `omics readouts. LCMS was used to detect and quantify proteins, metabolites, and lipids. In addition, cfDNA and mRNA were assayed using next-generation sequencing. cfDNA reads were analyzed to estimate fragment-lengths, copy-number variation, and CpG site methylation. All molecular data were normalized using standard methods specific to each assay. Univariate analyses of cases vs controls were performed to identify differentially abundant features on all available samples per assay. We detected 9,868 proteins, 605 lipids, 329 metabolites, and 109,070 mRNA transcripts. Of these, 3,098 proteins, 210 lipids, 57 metabolites, and 30,236 mRNA transcripts were significantly different (FWER < 0.05) in cases versus controls. Gene set enrichment analysis on statistically significant transcripts and proteins identified multiple gene-ontology terms associated with cancer including the Wnt signaling process and IgA immunoglobulin complex, respectively. From cfDNA data, we identified 234 non-contiguous genomic regions associated with the fragment-length disorder, 4,790 with copy-number variation, and 74 differentially methylated genomic regions spanning 184 CpG sites (FWER < 0.05). With the premise that deviations from copy number neutrality are more likely to indicate a tumor contribution, we then focused our examination on those differentially expressed proteins that overlap with differentially expressed mRNA transcripts as well as CNV genomic regions. We identified 52 protein coding genes including E-cadherin (associated with EMT) and related binding proteins such as RAB11B, CAPZB, EPS15, FLNB, MYH9, STK24 and YWHAE. Ongoing machine-learning-based classifier training to distinguish between cancer and non-cancer can serve as the basis for the development of high-sensitivity liquid-biopsy tests for lung cancer. Citation Format: Jinlyung Choi, Ajinkya Kokate, Ehdieh Khaledian, Manway Liu, Preethi Prasad, John Blume, Jessica Chan, Rea Cuaresma, Kevin Dai, Manoj Khadka, Thidar Khin, Yuya Kodama, Joon-Yong Lee, Hoda Malekpour, Megan Mora, Nithya Mudaliar, Sara Nouri Golmaei, Madhuvanthi Ramaiah, Saividya Ramaswamy, Peter Spiro, Dijana Vitko, Kavya Swaminathan, James Yee, Brian Young, Chinmay Belthangady, Bruce Wilcox, Brian Koh, Philip Ma. Biomarker discovery in non-small-cell lung cancer enabled by deep multi-omics profiling of proteins, metabolites, transcripts, and genes in blood. [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2023; Part 1 (Regular and Invited Abstracts); 2023 Apr 14-19; Orlando, FL. Philadelphia (PA): AACR; Cancer Res 2023;83(7_Suppl):Abstract nr 6606.

Pancreatic ductal adenocarcinoma (PDAC) is currently the 3rd leading cause of cancer-related deaths in the US. Although the all-stage 5-year survival rate is ~10%, early-stage 5-year survival is markedly superior and in excess of 40%. Hence, early detection of PDAC via blood-based liquid biopsies holds promise to reduce morbidity and mortality. PrognomiQ’s multi-omics platform performs deep and unbiased molecular profiling of blood samples to detect proteins, metabolites, lipids, mRNA, miRNA, cfDNA fragmentation and copy-number, and CpG methylation. Here we report results from training and validation of a classifier on a subset of that multi-omic data with the potential to enable the development of high sensitivity and specificity tests for early detection of PDAC.We conducted a case-control study comprising 146 subjects across 16 clinical sites, including 63 pathology-confirmed, untreated PDAC cases (12 stage I, 8 stage II, 4 stage III, 36 stage IV, and 3 stage unknown) and 83 age- and gender- matched controls without any known cancer. For each subject, venous blood samples including plasma were collected. Unbiased LCMS was used to detect and quantify proteins, and targeted, multiplexed MRM-LCMS assays were used for both metabolites and lipids. After data processing, we detected 54,114 proteomic features, 898 lipids, and 373 metabolites. 445 proteomic features, 170 lipids, and 37 metabolites were found to be significantly different as determined by Bonferroni-corrected Wilcoxon tests with FWER < 0.05. For classification, the dataset was split into training (37 cases and 37 controls) and validation (26 cases and 46 controls) sets, with control for collection site and date, age, and gender. XGBoost models were constructed for each analyte class using ten repeats of 10-fold cross-validation. To improve specificity to PDAC, all proteomic features which mapped to GOBP terms associated with acute-phase response, inflammation, and immune response were excluded prior to training. The best-performing hyperparameters were used for a final model built on the full training set and then used for inference on the validation set. At 99% specificity, the proteomic classifier had sensitivities of 77%, 57%, and 88% for Stages 1-4, Stages 1-2, and Stages 3-4, respectively, estimated by bootstrap re-sampling of the validation results. Metabolomics had sensitivities of 81%, 71%, and 88%. Lipidomics had sensitivities of 65%, 71%, and 65%. A joint, multi-omic model was constructed by averaging the scaled probabilities of all models. This joint model improved performance at 99% specificity with sensitivities of 92%, 86%, and 94%, highlighting the synergy of multi-omics data, particularly phenotypically related omics such as those described here. Multi-omic classifiers such as these can serve as the foundation for blood-based liquid biopsies for the early detection of PDAC. Citation Format: John Blume, Ghristine Bundalian, Jessica Chan, Connie Chao-Shern, Jinlyung Choi, Rea Cuaresma, Kevin Dai, Sara N. Golmaei, Jun Heok Jang, Manoj Khadka, Ehdieh Khaledian, Thidar Khin, Yuya Kodama, Ajinkya Kokate, Joon-Yong Lee, Manway Liu, Hoda Malekpour, Megan Mora, Nithya Mudaliar, Preethi Prasad, Madhuvanthi Ramaiah, Saividya Ramaswamy, Peter Spiro, Kavya Swaminathan, Dijana Vitko, James Yee, Brian Young, Susan Zhang, Chinmay Belthangady, Bruce Wilcox, Brian Koh, Philip Ma. A multi-omics classifier achieves high sensitivity and specificity for pancreatic ductal adenocarcinoma in a case-control study of 146 subjects [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2023; Part 1 (Regular and Invited Abstracts); 2023 Apr 14-19; Orlando, FL. Philadelphia (PA): AACR; Cancer Res 2023;83(7_Suppl):Abstract nr 6597.

Current coronavirus disease 2019 (COVID-19) vaccines effectively reduce overall morbidity and mortality and are vitally important to controlling the pandemic. Individuals who previously recovered from COVID-19 have enhanced immune responses after vaccination (hybrid immunity) compared with their naïve-vaccinated peers; however, the effects of post-vaccination breakthrough infections on humoral immune response remain to be determined. Here, we measure neutralizing antibody responses from 104 vaccinated individuals, including those with breakthrough infections, hybrid immunity, and no infection history. We find that human immune sera after breakthrough infection and vaccination after natural infection broadly neutralize SARS-CoV-2 (severe acute respiratory coronavirus 2) variants to a similar degree. Although age negatively correlates with antibody response after vaccination alone, no correlation with age was found in breakthrough or hybrid immune groups. Together, our data suggest that the additional antigen exposure from natural infection substantially boosts the quantity, quality, and breadth of humoral immune response regardless of whether it occurs before or after vaccination.

Climate change is globally affecting rainfall patterns, necessitating the improvement of drought tolerance in crops. Sorghum bicolor is a drought-tolerant cereal capable of producing high yields under water scarcity conditions. Functional stay-green sorghum genotypes can maintain green leaf area and efficient grain filling in terminal post-flowering water deprivation, a period of ~10 weeks. To obtain molecular insights into these characteristics, two drought-tolerant genotypes, BTx642 and RTx430, were grown in control and terminal post-flowering drought field plots in the Central Valley of California. Photosynthetic, photoprotective, water dynamics, and biomass traits were quantified and correlated with metabolomic data collected from leaves, stems, and roots at multiple timepoints during drought. Physiological and metabolomic data was then compared to longitudinal RNA sequencing data collected from these two genotypes. The metabolic response to drought highlights the uniqueness of the post-flowering drought acclimation relative to pre-flowering drought. The functional stay-green genotype BTx642 specifically induced photoprotective responses in post-flowering drought supporting a putative role for photoprotection in the molecular basis of the functional stay-green trait. Specific genes are highlighted that may contribute to post-flowering drought tolerance and that can be targeted in crops to maximize yields under limited water input conditions.HighlightPathways contributing to the long-term maintenance of photosynthetic activity in terminal post-flowering drought are revealed by a comprehensive approach combining in-field photosynthetic physiological analysis, metabolomics, and transcriptomics.

Pancreatic cancer is the third leading cause of cancer-related deaths in the United States. Disease biomarkers quantified from blood-based assays may help reduce mortality by enabling early detection, treatment selection, or response and resistance assessment. PrognomiQ has developed a multi-omics assay and analysis platform that comprehensively profiles blood samples to detect proteins, metabolites, lipids, mRNA, miRNA, cfDNA fragments, and methylation at CpG sites. This platform can provide deep insights into the biology of pancreatic cancer and could enable the development of high sensitivity and specificity tests for the early detection of pancreatic cancer. We conducted a case-control study comprising 196 subjects: 92 with untreated pancreatic cancer and 104 matched controls without pancreatic cancer. For each subject, blood was collected in assay-specific tubes and processed to provide 7 different `omics readouts. cfDNA and mRNA were isolated from samples and assayed following standard NGS protocols. cfDNA fragments were processed to estimate fragment-length disorder and copy-number variation along with CpG site methylation. In addition, targeted and untargeted LCMS were used to detect and quantify proteins, metabolites, and lipids. After normalization, non-parametric univariate analyses of cases versus controls were performed to identify differentially abundant features on all available samples for each assay. Unsupervised learning was used to investigate the separation of subjects into groups based on disease status for the subset of 157 subjects for which complete data on all 7 readouts were available. We detected 2,812 proteins, 811 lipids, 373 metabolites, and 110,864 mRNA transcripts in all samples where data for each assay was available. Of these, 275 proteins, 232 lipids, 49 metabolites, and 3385 mRNA transcripts were significantly different (FWER < 0.05) in cases versus controls. From cfDNA data, we identified 35 non-contiguous genomic regions associated with fragment-length disorder, 557 with copy-number variation, and 5 with multiple, differentially methylated CpGs (FWER < 0.05) that aggregately span 307 protein-coding genes; of these, the overlap with the differentially expressed proteins included E-cadherin (tumor suppressor) and N-cadherin (involved in epithelial-to-mesenchymal transition). Statistically significant genes and proteins were found to be associated with processes including Wnt signaling, regulation of focal adhesion assembly, and actin cytoskeleton organization. Multi-omics, unsupervised learning showed separation of early- and late-stage cases and controls. High-dimensional bioinformatics analyses systematically decomposed each `omics data type into joint and orthogonal components associated with pancreatic cancer. Ongoing multivariate analyses, including supervised machine learning, will further elucidate the biology of pancreatic cancer development, and serve as the foundation for high-sensitivity blood tests for the early detection and monitoring of pancreatic cancer. Citation Format: Ehdieh Khaledian, Preethi Prasad, John Blume, Ghristine Bundalian, Connie Chao-Shern, Jinlyung Choi, Rea Cuaresma, Jared Deyarmin, Jun Heok Jang, Manoj Khadka, Thidar Khin, Yuya Kodama, Ajinkya Kokate, Joon-Yong Lee, Manway Liu, Nithya Mudaliar, Madhuvanthi Ramaiah, Saividya Ramaswamy, Peter Spiro, Kavya Swaminathan, Preston Williams, Mi Yang, James Yee, Brian Young, Robert Zawada, Susan Zhang, Chinmay Belthangady, Bruce Wilcox, Philip Ma. High-dimensional, multi-omics analyses of proteins, metabolites, transcripts, and genes enable biomarker discovery in early- and late-stage pancreatic cancer [abstract]. In: Proceedings of the AACR Special Conference on Pancreatic Cancer; 2022 Sep 13-16; Boston, MA. Philadelphia (PA): AACR; Cancer Res 2022;82(22 Suppl):Abstract nr A038.