Multi-omics Data Integration Research Articles

Metabolomics of head and neck cancer in biofluids: an integrative systematic review.

Head and neck cancer (HNC) is the fifth most common cancer globally. Diagnosis at early stages are critical to reduce mortality and improve functional and esthetic outcomes associated with HNC. Metabolomics is a promising approach for discovery of biomarkers and metabolic pathways for risk assessment and early detection of HNC. To summarize and consolidate the available evidence on metabolomics and HNC in plasma/serum, saliva, and urine. A systematic search of experimental research was executed using PubMed and Web of Science. Available data on areas under the curve was extracted. Metabolic pathway enrichment analysis were performed to identify metabolic pathways altered in HNC. Fifty-four studies were eligible for data extraction (33 performed in plasma/serum, 15 in saliva and 6 in urine). Metabolites with high discriminatory performance for detection of HNC included single metabolites and combination panels of several lysoPCs, pyroglutamate, glutamic acid, glucose, tartronic acid, arachidonic acid, norvaline, linoleic acid, propionate, acetone, acetate, choline, glutamate and others. The glucose-alanine cycle and the urea cycle were the most altered pathways in HNC, among other pathways (i.e. gluconeogenesis, glycine and serine metabolism, alanine metabolism, etc.). Specific metabolites that can potentially serve as complementary less- or non-invasive biomarkers, as well as metabolic pathways integrating the data from the available studies, are presented. The present work highlights utility of metabolite-based biomarkers for risk assessment, early detection, and prognostication of HNC, as well as facilitates incorporation of available metabolomics studies into multi-omics data integration and big data analytics for personalized health.

Multiomic Investigations into Lung Health and Disease.

Diseases of the lung account for more than 5 million deaths worldwide and are a healthcare burden. Improving clinical outcomes, including mortality and quality of life, involves a holistic understanding of the disease, which can be provided by the integration of lung multi-omics data. An enhanced understanding of comprehensive multiomic datasets provides opportunities to leverage those datasets to inform the treatment and prevention of lung diseases by classifying severity, prognostication, and discovery of biomarkers. The main objective of this review is to summarize the use of multiomics investigations in lung disease, including multiomics integration and the use of machine learning computational methods. This review also discusses lung disease models, including animal models, organoids, and single-cell lines, to study multiomics in lung health and disease. We provide examples of lung diseases where multi-omics investigations have provided deeper insight into etiopathogenesis and have resulted in improved preventative and therapeutic interventions.

Integration of multi-omics data using adaptive graph learning and attention mechanism for patient classification and biomarker identification

With the rapid development and accumulation of high-throughput sequencing technology and omics data, many studies have conducted a more comprehensive understanding of human diseases from a multi-omics perspective. Meanwhile, graph-based methods have been widely used to process multi-omics data due to its powerful expressive ability. However, most existing graph-based methods utilize fixed graphs to learn sample embedding representations, which often leads to sub-optimal results. Furthermore, treating embedding representations of different omics equally usually cannot obtain more reasonable integrated information. In addition, the complex correlation between omics is not fully taken into account. To this end, we propose an end-to-end interpretable multi-omics integration method, named MOGLAM, for disease classification prediction. Dynamic graph convolutional network with feature selection is first utilized to obtain higher quality omic-specific embedding information by adaptively learning the graph structure and discover important biomarkers. Then, multi-omics attention mechanism is applied to adaptively weight the embedding representations of different omics, thereby obtaining more reasonable integrated information. Finally, we propose omic-integrated representation learning to capture complex common and complementary information between omics while performing multi-omics integration. Experimental results on three datasets show that MOGLAM achieves superior performance than other state-of-the-art multi-omics integration methods. Moreover, MOGLAM can identify important biomarkers from different omics data types in an end-to-end manner.

Multiomics: A user-friendly multi-omics data harmonisation R pipeline

Data from multiple omics layers of a biological system is growing in quantity, heterogeneity and dimensionality. Simultaneous multi-omics data integration is of immense interest to researchers as it has potential to unlock previously hidden biomolecular relationships leading to early diagnosis, prognosis, and expedited treatments. Many tools for multi-omics data integration are developed. However, these tools are often restricted to highly specific experimental designs, types of omics data, and specific data formats. A major limitation of the field is the lack of a pipeline that can accept data in unrefined form to preserve maximum biology in an individual dataset prior to integration. We fill this gap by developing a flexible, generic multi-omics pipeline called multiomics, to facilitate general-purpose data exploration and analysis of heterogeneous data. The pipeline takes unrefined multi-omics data as input, sample information and user-specified parameters to generate a list of output plots and data tables for quality control and downstream analysis. We have demonstrated its application on a sepsis case study. We enabled limited checkpointing functionality where intermediate output is staged to allow continuation after errors or interruptions in the pipeline and generate a script for reproducing the analysis to improve reproducibility. Our pipeline can be installed as an R package or manually from the git repository, and is accompanied by detailed documentation with walkthroughs on three case studies.

Proteogenomic data and resources for pan-cancer analysis.

The National Cancer Institute's Clinical Proteomic Tumor Analysis Consortium (CPTAC) investigates tumors from a proteogenomic perspective, creating rich multi-omics datasets connecting genomic aberrations to cancer phenotypes. To facilitate pan-cancer investigations, we have generated harmonized genomic, transcriptomic, proteomic, and clinical data for >1000 tumors in 10 cohorts to create a cohesive and powerful dataset for scientific discovery. We outline efforts by the CPTAC pan-cancer working group in data harmonization, data dissemination, and computational resources for aiding biological discoveries. We also discuss challenges for multi-omics data integration and analysis, specifically the unique challenges of working with both nucleotide sequencing and mass spectrometry proteomics data.

Applications of multi-omics analysis in human diseases.

Multi-omics usually refers to the crossover application of multiple high-throughput screening technologies represented by genomics, transcriptomics, single-cell transcriptomics, proteomics and metabolomics, spatial transcriptomics, and so on, which play a great role in promoting the study of human diseases. Most of the current reviews focus on describing the development of multi-omics technologies, data integration, and application to a particular disease; however, few of them provide a comprehensive and systematic introduction of multi-omics. This review outlines the existing technical categories of multi-omics, cautions for experimental design, focuses on the integrated analysis methods of multi-omics, especially the approach of machine learning and deep learning in multi-omics data integration and the corresponding tools, and the application of multi-omics in medical researches (e.g., cancer, neurodegenerative diseases, aging, and drug target discovery) as well as the corresponding open-source analysis tools and databases, and finally, discusses the challenges and future directions of multi-omics integration and application in precision medicine. With the development of high-throughput technologies and data integration algorithms, as important directions of multi-omics for future disease research, single-cell multi-omics and spatial multi-omics also provided a detailed introduction. This review will provide important guidance for researchers, especially who are just entering into multi-omics medical research.

Integration of multi-omics data to elucidate keystone unknown taxa within microbialite-forming ecosystems.

Microbes continually shape Earth's biochemical and physical landscapes by inhabiting diverse metabolic niches. Despite the important role microbes play in ecosystem functioning, most microbial species remain unknown highlighting a gap in our understanding of structured complex ecosystems. To elucidate the relevance of these unknown taxa, often referred to as "microbial dark matter," the integration of multiple high throughput sequencing technologies was used to evaluate the co-occurrence and connectivity of all microbes within the community. Since there are no standard methodologies for multi-omics integration of microbiome data, we evaluated the abundance of "microbial dark matter" in microbialite-forming communities using different types meta-omic datasets: amplicon, metagenomic, and metatranscriptomic sequencing previously generated for this ecosystem. Our goal was to compare the community structure and abundances of unknown taxa within the different data types rather than to perform a functional characterization of the data. Metagenomic and metatranscriptomic data were input into SortMeRNA to extract 16S rRNA gene reads. The output, as well as amplicon sequences, were processed through QIIME2 for taxonomy analysis. The R package mdmnets was utilized to build co-occurrence networks. Most hubs presented unknown classifications, even at the phyla level. Comparisons of the highest scoring hubs of each data type using sequence similarity networks allowed the identification of the most relevant hubs within the microbialite-forming communities. This work highlights the importance of unknown taxa in community structure and proposes that ecosystem network construction can be used on several types of data to identify keystone taxa and their potential function within microbial ecosystems.

AFEI: adaptive optimized vertical federated learning for heterogeneous multi-omics data integration.

Vertical federated learning has gained popularity as a means of enabling collaboration and information sharing between different entities while maintaining data privacy and security. This approach has potential applications in disease healthcare, cancer prognosis prediction, and other industries where data privacy is a major concern. Although using multi-omics data for cancer prognosis prediction provides more information for treatment selection, collecting different types of omics data can be challenging due to their production in various medical institutions. Data owners must comply with strict data protection regulations such as European Union (EU) General Data Protection Regulation. To share patient data across multiple institutions, privacy and security issues must be addressed. Therefore, we propose an adaptive optimized vertical federated-learning-based framework adaptive optimized vertical federated learning for heterogeneous multi-omics data integration (AFEI) to integrate multi-omics data collected from multiple institutions for cancer prognosis prediction. AFEI enables participating parties to build an accurate joint evaluation model for learning more information related to cancer patients from different perspectives, based on the distributed and encrypted multi-omics features shared by multiple institutions. The experimental results demonstrate that AFEI achieves higher prediction accuracy (6.5% on average) than using single omics data by utilizing the encrypted multi-omics data from different institutions, and it performs almost as well as prognosis prediction by directly integrating multi-omics data. Overall, AFEI can be seen as an efficient solution for breaking down barriers to multi-institutional collaboration and promoting the development of cancer prognosis prediction.

Integrating omics databases for enhanced crop breeding.

Crop plant breeding involves selecting and developing new plant varieties with desirable traits such as increased yield, improved disease resistance, and enhanced nutritional value. With the development of high-throughput technologies, such as genomics, transcriptomics, and metabolomics, crop breeding has entered a new era. However, to effectively use these technologies, integration of multi-omics data from different databases is required. Integration of omics data provides a comprehensive understanding of the biological processes underlying plant traits and their interactions. This review highlights the importance of integrating omics databases in crop plant breeding, discusses available omics data and databases, describes integration challenges, and highlights recent developments and potential benefits. Taken together, the integration of omics databases is a critical step towards enhancing crop plant breeding and improving global food security.

Integration of single-cell regulon atlas and multi-omics data for prognostic stratification and personalized treatment prediction in human lung adenocarcinoma

Transcriptional programs are often dysregulated in cancers. A comprehensive investigation of potential regulons is critical to the understanding of tumorigeneses. We first constructed the regulatory networks from single-cell RNA sequencing data in human lung adenocarcinoma (LUAD). We next introduce LPRI (Lung Cancer Prognostic Regulon Index), a precision oncology framework to identify new biomarkers associated with prognosis by leveraging the single cell regulon atlas and bulk RNA sequencing or microarray datasets. We confirmed that LPRI could be a robust biomarker to guide prognosis stratification across lung adenocarcinoma cohorts. Finally, a multi-omics data analysis to characterize molecular alterations associated with LPRI was performed from The Cancer Genome Atlas (TCGA) dataset. Our study provides a comprehensive chart of regulons in LUAD. Additionally, LPRI will be used to help prognostic prediction and developing personalized treatment for future studies.

A review of multi-omics data integration through deep learning approaches for disease diagnosis, prognosis, and treatment.

Accurate diagnosis is the key to providing prompt and explicit treatment and disease management. The recognized biological method for the molecular diagnosis of infectious pathogens is polymerase chain reaction (PCR). Recently, deep learning approaches are playing a vital role in accurately identifying disease-related genes for diagnosis, prognosis, and treatment. The models reduce the time and cost used by wet-lab experimental procedures. Consequently, sophisticated computational approaches have been developed to facilitate the detection of cancer, a leading cause of death globally, and other complex diseases. In this review, we systematically evaluate the recent trends in multi-omics data analysis based on deep learning techniques and their application in disease prediction. We highlight the current challenges in the field and discuss how advances in deep learning methods and their optimization for application is vital in overcoming them. Ultimately, this review promotes the development of novel deep-learning methodologies for data integration, which is essential for disease detection and treatment.

Progress in single-cell multimodal sequencing and multi-omics data integration.

With the rapid advance of single-cell sequencing technology, cell heterogeneity in various biological processes was dissected at different omics levels. However, single-cell mono-omics results in fragmentation of information and could not provide complete cell states. In the past several years, a variety of single-cell multimodal omics technologies have been developed to jointly profile multiple molecular modalities, including genome, transcriptome, epigenome, and proteome, from the same single cell. With the availability of single-cell multimodal omics data, we can simultaneously investigate the effects of genomic mutation or epigenetic modification on transcription and translation, and reveal the potential mechanisms underlying disease pathogenesis. Driven by the massive single-cell omics data, the integration method of single-cell multi-omics data has rapidly developed. Integration of the massive multi-omics single-cell data in public databases in the future will make it possible to construct a cell atlas of multi-omics, enabling us to comprehensively understand cell state and gene regulation at single-cell resolution. In this review, we summarized the experimental methods for single-cell multimodal omics data and computational methods for multi-omics data integration. We also discussed the future development of this field.

Integration of gene expression and DNA methylation data across different experiments.

Integrative analysis of multi-omic datasets has proven to be extremely valuable in cancer research and precision medicine. However, obtaining multimodal data from the same samples is often difficult. Integrating multiple datasets of different omics remains a challenge, with only a few available algorithms developed to solve it. Here, we present INTEND (IntegratioN of Transcriptomic and EpigeNomic Data), a novel algorithm for integrating gene expression and DNA methylation datasets covering disjoint sets of samples. To enable integration, INTEND learns a predictive model between the two omics by training on multi-omic data measured on the same set of samples. In comprehensive testing on 11 TCGA (The Cancer Genome Atlas) cancer datasets spanning 4329 patients, INTEND achieves significantly superior results compared with four state-of-the-art integration algorithms. We also demonstrate INTEND's ability to uncover connections between DNA methylation and the regulation of gene expression in the joint analysis of two lung adenocarcinoma single-omic datasets from different sources. INTEND's data-driven approach makes it a valuable multi-omic data integration tool. The code for INTEND is available at https://github.com/Shamir-Lab/INTEND.

Idiopathic Pulmonary Arterial Hypertension: Network-Based Integration of Multi-Omics Data Reveals New Molecular Signatures and Candidate Drugs.

Idiopathic pulmonary arterial hypertension (IPAH) is a progressive disease that affects the pulmonary arteries, resulting in increased pulmonary vascular resistance and right ventricular dysfunction, which can ultimately lead to heart failure and death. The molecular substrates of IPAH are poorly understood while diagnostics and therapeutics innovation remain as unmet needs for this debilitating disease. In this study, a network-based methodology was used to uncover the salient molecular mechanisms of IPAH to inform drug and diagnostic discovery, and personalized medicine. Expression profiling datasets associated with IPAH were obtained from the Gene Expression Omnibus database: GSE15197, GSE113439, GSE53408, and GSE67597. The comparative analysis of mRNA and miRNA expression data and the modular analysis of a transcriptome-based weighted gene coexpression network unraveled disease-specific gene and miRNA signatures. DEAD-box helicase 52 (DDx52), ESF1 nucleolar pre-RNA processing protein (ESF1), heterogeneous nuclear ribonuclearprotein A3 (MNRNPA3), Myosin VA (MYO5A), replication factor C subunit 1 (RFC1), and arginine and serine rich coiled coil 1 (RSRC1) were detected as the salient genes for IPAH. In addition, the salient gene-based drug repositioning analysis identified alvespimycin, tanespimycin, geldanamycin, LY294002, cephaeline, digoxigenin, lanatoside C, helveticoside, trichostatin A, phenoxybenzamine, genistein, pioglitazone, and rosiglitazone as potential drug candidates for IPAH. In conclusion, this study provides new molecular signatures in relation to IPAH and attendant potential drug candidates for further experimental and translational clinical research for patients with IPAH.

MMDAE-HGSOC: A novel method for high-grade serous ovarian cancer molecular subtypes classification based on multi-modal deep autoencoder

High-grade serous ovarian cancer (HGSOC) is a type of ovarian cancer developed from serous tubal intraepithelial carcinoma. The intrinsic differences among molecular subtypes are closely associated with prognosis and pathological characteristics. At present, multi-omics data integration methods include early integration and late integration. Most existing HGSOC molecular subtypes classification methods are based on the early integration of multi-omics data. The mutual interference among multi-omics data is ignored, which affects the effectiveness of feature learning. High-dimensional multi-omics data contains genes unassociated with HGSOC molecular subtypes, resulting in redundant information, which is not conducive to model training. In this paper, we propose a multi-modal deep autoencoder learning method, MMDAE-HGSOC. MiRNA expression, DNA methylation, and copy number variation (CNV) are integrated with mRNA expression data to construct a multi-omics feature space. The multi-modal deep autoencoder network is used to learn the high-level feature representation of multi-omics data. The superposition LASSO (S-LASSO) regression algorithm is proposed to fully obtain the associated genes of HGSOC molecular subtypes. The experimental results show that MMDAE-HGSOC is superior to the existing classification methods. Finally, we analyze the enrichment gene ontology (GO) terms and biological pathways of these significant genes, which are discovered during the gene selection process.

Lung cancer multi-omics digital human avatars for integrating precision medicine into clinical practice: the LANTERN study

BackgroundThe current management of lung cancer patients has reached a high level of complexity. Indeed, besides the traditional clinical variables (e.g., age, sex, TNM stage), new omics data have recently been introduced in clinical practice, thereby making more complex the decision-making process. With the advent of Artificial intelligence (AI) techniques, various omics datasets may be used to create more accurate predictive models paving the way for a better care in lung cancer patients.MethodsThe LANTERN study is a multi-center observational clinical trial involving a multidisciplinary consortium of five institutions from different European countries. The aim of this trial is to develop accurate several predictive models for lung cancer patients, through the creation of Digital Human Avatars (DHA), defined as digital representations of patients using various omics-based variables and integrating well-established clinical factors with genomic data, quantitative imaging data etc. A total of 600 lung cancer patients will be prospectively enrolled by the recruiting centers and multi-omics data will be collected. Data will then be modelled and parameterized in an experimental context of cutting-edge big data analysis. All data variables will be recorded according to a shared common ontology based on variable-specific domains in order to enhance their direct actionability. An exploratory analysis will then initiate the biomarker identification process. The second phase of the project will focus on creating multiple multivariate models trained though advanced machine learning (ML) and AI techniques for the specific areas of interest. Finally, the developed models will be validated in order to test their robustness, transferability and generalizability, leading to the development of the DHA. All the potential clinical and scientific stakeholders will be involved in the DHA development process. The main goals aim of LANTERN project are: i) To develop predictive models for lung cancer diagnosis and histological characterization; (ii) to set up personalized predictive models for individual-specific treatments; iii) to enable feedback data loops for preventive healthcare strategies and quality of life management.DiscussionThe LANTERN project will develop a predictive platform based on integration of multi-omics data. This will enhance the generation of important and valuable information assets, in order to identify new biomarkers that can be used for early detection, improved tumor diagnosis and personalization of treatment protocols.Ethics Committee approval number5420 − 0002485/23 from Fondazione Policlinico Universitario Agostino Gemelli IRCCS – Università Cattolica del Sacro Cuore Ethics Committee.Trial registrationclinicaltrial.gov - NCT05802771.

Integrative multi-omics analysis reveals novel idiopathic pulmonary fibrosis endotypes associated with disease progression

BackgroundIdiopathic pulmonary fibrosis (IPF) is characterized by the accumulation of extracellular matrix in the pulmonary interstitium and progressive functional decline. We hypothesized that integration of multi-omics data would identify clinically meaningful molecular endotypes of IPF.MethodsThe IPF-PRO Registry is a prospective registry of patients with IPF. Proteomic and transcriptomic (including total RNA [toRNA] and microRNA [miRNA]) analyses were performed using blood collected at enrollment. Molecular data were integrated using Similarity Network Fusion, followed by unsupervised spectral clustering to identify molecular subtypes. Cox proportional hazards models tested the relationship between these subtypes and progression-free and transplant-free survival. The molecular subtypes were compared to risk groups based on a previously described 52-gene (toRNA expression) signature. Biological characteristics of the molecular subtypes were evaluated via linear regression differential expression and canonical pathways (Ingenuity Pathway Analysis [IPA]) over-representation analyses.ResultsAmong 232 subjects, two molecular subtypes were identified. Subtype 1 (n = 105, 45.3%) and Subtype 2 (n = 127, 54.7%) had similar distributions of age (70.1 +/- 8.1 vs. 69.3 +/- 7.6 years; p = 0.31) and sex (79.1% vs. 70.1% males, p = 0.16). Subtype 1 had more severe disease based on composite physiologic index (CPI) (55.8 vs. 51.2; p = 0.002). After adjusting for CPI and antifibrotic treatment at enrollment, subtype 1 experienced shorter progression-free survival (HR 1.79, 95% CI 1.28,2.56; p = 0.0008) and similar transplant-free survival (HR 1.30, 95% CI 0.87,1.96; p = 0.20) as subtype 2. There was little agreement in the distribution of subjects to the molecular subtypes and the risk groups based on 52-gene signature (kappa = 0.04, 95% CI= -0.08, 0.17), and the 52-gene signature risk groups were associated with differences in transplant-free but not progression-free survival. Based on heatmaps and differential expression analyses, proteins and miRNAs (but not toRNA) contributed to classification of subjects to the molecular subtypes. The IPA showed enrichment in pulmonary fibrosis-relevant pathways, including mTOR, VEGF, PDGF, and B-cell receptor signaling.ConclusionsIntegration of transcriptomic and proteomic data from blood enabled identification of clinically meaningful molecular endotypes of IPF. If validated, these endotypes could facilitate identification of individuals likely to experience disease progression and enrichment of clinical trials.Trial registrationNCT01915511

Integration of Multi-omic Data in a Molecular Tumor Board Reveals EGFR-Associated ALK-Inhibitor Resistance in a Patient With Inflammatory Myofibroblastic Cancer.

Inflammatory myofibroblastic tumors (IMTs) are intermediate-grade mesenchymal neoplasms commonly characterized by chromosomal rearrangements causing constitutive activation of anaplastic lymphoma kinase (ALK) and/or ALK mutations causing reduced sensitivity to ALK tyrosine kinase inhibitors (TKI). We present a patient with an IMT who initially responded to first-line alectinib, but who later suffered disease relapse and presently survives with moderate residual disease after receiving second-line lorlatinib. Biopsy specimens were analyzed using next generation sequencing (DNA-seq and RNA-seq) and reverse phase protein microarray (RPPA) as part of an institutional Molecular Tumor Board (MTB) study. An EML4-ALK rearrangement and EGFR activation (pEGFRY1068) were present in both the primary and recurrent tumors, while a secondary ALK I1171N mutation was exclusive to the latter. EGFR signaling in the background of a secondary ALK mutation is correlated with reduced ALK TKI sensitivity in vitro, implicating an important mechanism of drug resistance development in this patient. The RPPA results also critically demonstrate that ALK signaling (ALKY1604) was not activated in the recurrent tumor, thereby indicating that standard-of-care use of third- or fourth-line ALK TKI would not likely be efficacious or durable. These results underscore the importance of real-time clinical integration of functional protein drug target activation data with NGS in the MTB setting for improving selection of patient-tailored therapy.

MRGCN: cancer subtyping with multi-reconstruction graph convolutional network using full and partial multi-omics dataset.

Cancer is a molecular complex and heterogeneous disease. Each type of cancer is usually composed of several subtypes with different treatment responses and clinical outcomes. Therefore, subtyping is a crucial step in cancer diagnosis and therapy. The rapid advances in high-throughput sequencing technologies provide an increasing amount of multi-omics data, which benefits our understanding of cancer genetic architecture, and yet poses new challenges in multi-omics data integration. We propose a graph convolutional network model, called MRGCN for multi-omics data integrative representation. MRGCN simultaneously encodes and reconstructs multiple omics expression and similarity relationships into a shared latent embedding space. In addition, MRGCN adopts an indicator matrix to denote the situation of missing values in partial omics, so that the full and partial multi-omics processing procedures are combined in a unified framework. Experimental results on 11 multi-omics datasets show that cancer subtypes obtained by MRGCN with superior enriched clinical parameters and log-rank test P-values in survival analysis over many typical integrative methods. https://github.com/Polytech-bioinf/MRGCN.git https://figshare.com/articles/software/MRGCN/23058503.

Multi-omic assessment shows dysregulation of pulmonary and systemic immunity to e-cigarette exposure

Electronic cigarette (Ecig) use has become more common, gaining increasing acceptance as a safer alternative to tobacco smoking. However, the 2019 outbreak of Ecig and Vaping-Associated Lung Injury (EVALI) alerted the community to the potential for incorporation of deleterious ingredients such as vitamin E acetate into products without adequate safety testing. Understanding Ecig induced molecular changes in the lung and systemically can provide a path to safety assessment and protect consumers from unsafe formulations. While vitamin E acetate has been largely removed from commercial and illicit products, many Ecig products contain additives that remain largely uncharacterized. In this study, we determined the lung-specific effects as well as systemic immune effects in response to exposure to a common Ecig base, propylene glycol and vegetable glycerin (PGVG), with and without a 1% addition of phytol, a diterpene alcohol that has been found in commercial products. We exposed animals to PGVG with and without phytol and assessed metabolite, lipid, and transcriptional markers in the lung. We found both lung-specific as well as systemic effects in immune parameters, metabolites, and lipids. Phytol drove modest changes in lung function and increased splenic CD4 T cell populations. We also conducted multi-omic data integration to better understand early complex pulmonary responses, highlighting a central enhancement of acetylcholine responses and downregulation of palmitic acid connected with conventional flow cytometric assessments of lung, systemic inflammation, and pulmonary function. Our results demonstrate that Ecig exposure not only leads to changes in pulmonary function but also affects systemic immune and metabolic parameters.