Multi-omics Datasets Research Articles

SoyOD: An Integrated Soybean Multi-omics Database for Mining Genes and Biological Research.

Soybean is a globally important crop for food, feed, oil, and nitrogen fixation. A variety of multi-omics studies has been carried out, generating datasets ranging from genotype to phenotype. In order to efficiently utilize these data for basic and applied research, a soybean multi-omics database with extensive data coverage and comprehensive data analysis tools was established. The Soybean Omics Database (SoyOD) integrates important new datasets with existing public datasets to form the most comprehensive collection of soybean multi-omics information. Compared to existing soybean databases, SoyOD incorporates an extensive collection of novel data derived from the deep-sequencing of 984 germplasms, 162 novel transcriptome datasets from seeds at different developmental stages, 53 phenotypic datasets, and more than 2500 phenotypic images. In addition, SoyOD integrates existing data resources, including 59 assembled genomes, genetic variation data from 3904 soybean accessions, 225 sets of phenotypic data, and 1097 transcriptomic sequences covering 507 different tissues and treatment conditions. Moreover, SoyOD can be used to mine candidate genes for important agronomic traits, as shown in a case study on plant height. Additionally, powerful analytical and easy-to-use toolkits enable users to easily access the available multi-omics datasets, and to rapidly search genotypic and phenotypic data in a particular germplasm. The novelty, comprehensiveness, and user-friendly features of SoyOD make it a valuable resource for soybean molecular breeding and biological research. SoyOD is publicly accessible at https://bis.zju.edu.cn/soyod.

Abstract 4120480: HIV-Nef extracellular vesicles utilize the novel Btk-NFκB-MerTK signaling axis to impair macrophage efferocytosis and promote atherosclerosis: A multiomics study

Background: People living with HIV (PLWH) on anti-retroviral therapy remain at risk for cardiovascular diseases, including atherosclerosis. We hypothesized that persistent viral protein (HIV-Nef) in extracellular vesicles (EVs) modulate macrophage heterogeneity to impair atheroprotective efferocytosis to accelerate cardiovascular disease. Methods and Results: Macrophage heterogeneity was characterized in human primary macrophages (50,931 cells; 4 donors) stimulated with EVs engineered to contain HIV-Nef by simultaneous scRNAseq and scATACseq. Among 16 clusters (Fig.1A), two inflammatory Nef dominant clusters were characterized using gene set enrichment analysis. Gene regulatory network analysis of scRNAseq (pySCENIC) and transcription factor footprinting analysis of sc-ATACseq (ChromVar) suggest elevated NFκB activation in these clusters. Whole cell and sub-cellular compartmentalized proteomics of nucleus, cytosol, cell surface and secreted EVs indicate changes in immune response related biological pathways. Network analysis of our multi-omics data predicts Bruton Tyrosine Kinase (Btk) signaling as a potential contributor to Nef induced macrophage dysregulation. Pathway enrichment analysis of these multiomics dataset suggest Nef influences atheroprotective efferocytosis through the Btk-NFκB signaling axis. Spectral flow cytometry, high content imaging and multiplexed qPCR showed reduction in key effector of efferocytosis, MerTK. Btk inhibition using siRNA, reversible and irreversible Btk inhibitors restored MerTK expression and rescued efferocytosis. Importantly, CRISPRa based overexpression of MerTK in human primary macrophages rescued Nef impaired efferocytosis. Injection of HIV-Nef EVs into male and female C57BL6/J mice impaired efferocytosis of peritoneal and bone marrow derived macrophages which was rescued with Btk inhibition in vivo . Critically, injection of HIV-Nef EV into male and female Ldlr -/- enhanced atherogenesis with larger aortic wall thickness and necrotic core (Fig.1B). Conclusion: Persistent Nef in EV in PLWH may modulate macrophage heterogeneity to impair efferocytosis. These findings may help develop novel atheroprotective therapies by restoring macrophage efferocytosis.

TMET-39. DUAL TARGETING OF MAP KINASE SIGNALING AND METABOLISM IS A THERAPEUTIC VULNERABILITY IN PEDIATRIC DIFFUSE MIDLINE GLIOMAS

Abstract Metabolic reprogramming in pediatric diffuse midline glioma is driven by gene expression changes induced by the hallmark histone mutation H3K27M, which results in aberrantly permissive activation of oncogenic signaling pathways. Previous studies of diffuse midline glioma with altered H3K27 (DMG-H3K27a) have shown that the RAS pathway, specifically through its downstream kinase, extracellular-signal-related kinase 5 (ERK5), is critical for tumor growth. However, there remains a knowledge gap in identifying effectors of ERK signaling and their roles in DMG-H3K27a. Leveraging several patient multi-omic datasets, we identified a unique relationship between ERK5 and brain tumor glycolysis. We establish that ERK5 is a critical regulator of cell proliferation and glycolysis in DMG-H3K27a, where ERK5 knockdown leads to decreases in glycolysis and increases in mitochondrial metabolism. We demonstrate through loss of function and overexpression studies that ERK5 mediates glycolysis through regulation of glycolytic enzyme Phosphofructo-2-Kinase/Fructose-2,6-Biphosphatase 3 (PFKFB3). PFKFB3 is a protein that catalyzes the synthesis or degradation of fructose-2,6-bisphosphate (F2,6BP). Furthermore, we establish a novel mechanism, in which ERK5 mediates PFKFB3 expression via activation of MEF2A, an oncogenic transcription factor. DMG-H3K27a cells are sensitive to the loss or inhibition of PFKFB3 both in vitro and in orthotopic mouse models. However, monotherapy is ineffective in heterogenous diseases such as gliomas, including DMG-H3K27a. To combat this, we show multi-targeted therapy against the ERK5-PFKFB3 axis is synergistic in vitro and in vivo. Multi-targeted drug therapy against the ERK5-PFKFB3 axis, such as with small-molecule inhibitors, may represent a promising therapeutic approach in patients with pediatric diffuse midline glioma. Moreover, we have identified resistance to our dual therapy using several systems biology approaches, which we are currently investigating and validating. Collectively, our study defines a novel, targetable axis in DMG-H3K27a and highlights the importance of tumor cell metabolism in defining therapeutic vulnerabilities.

N6-methyladenosine regulators in hepatocellular carcinoma: investigating the precise definition and clinical applications of biomarkers

BackgroundAccurately identifying effective biomarkers and translating them into clinical practice have significant implications for improving clinical outcomes in hepatocellular carcinoma (HCC). In this study, our objective is to explore appropriate methods to improve the accuracy of biomarker identification and investigate their clinical value.MethodsConcentrating on the N6-methyladenosine (m6A) modification regulators, we utilized dozens of multi-omics HCC datasets to analyze the expression patterns and genetic features of m6A regulators. Through the integration of big data analysis with function experiments, we have redefined the biological roles of m6A regulators in HCC. Based on the key regulators, we constructed m6A risk models and explored their clinical value in estimating prognosis and guiding personalized therapy for HCC.ResultsMost m6A regulators exhibit abnormal expression in HCC, and their expression is influenced by copy number variations (CNV) and DNA methylation. Large-scale data analysis has revealed the biological roles of many key m6A regulators, and these findings are well consistent with experimental results. The m6A risk models offer significant prognostic value. Moreover, they assist in reassessing the therapeutic potential of drugs such as sorafenib, gemcitabine, CTLA4 and PD1 blockers in HCC.ConclusionsOur findings suggest that the mutual validation of big data analysis and functional experiments may facilitate the precise identification and definition of biomarkers, and our m6A risk models may have the potential to guide personalized chemotherapy, targeted treatment, and immunotherapy decisions in HCC.

Single-Cell Multiomics Profiling Reveals Heterogeneity of Müller Cells in the Oxygen-Induced Retinopathy Model.

Retinal neovascularization poses heightened risks of vision loss and blindness. Despite its clinical significance, the molecular mechanisms underlying the pathogenesis of retinal neovascularization remain elusive. This study utilized single-cell multiomics profiling in an oxygen-induced retinopathy (OIR) model to comprehensively investigate the intricate molecular landscape of retinal neovascularization. Mice were exposed to hyperoxia to induce the OIR model, and retinas were isolated for nucleus isolation. The cellular landscape of the single-nucleus suspensions was extensively characterized through single-cell multiomics sequencing. Single-cell data were integrated with genome-wide association study (GWAS) data to identify correlations between ocular cell types and diabetic retinopathy. Cell communication analysis among cells was conducted to unravel crucial ligand-receptor signals. Trajectory analysis and dynamic characterization of Müller cells were performed, followed by integration with human retinal data for pathway analysis. The multiomics dataset revealed six major ocular cell classes, with Müller cells/astrocytes showing significant associations with proliferative diabetic retinopathy (PDR). Cell communication analysis highlighted pathways that are associated with vascular proliferation and neurodevelopment, such as Vegfa-Vegfr2, Igf1-Igf1r, Nrxn3-Nlgn1, and Efna5-Epha4. Trajectory analysis identified a subset of Müller cells expressing genes linked to photoreceptor degeneration. Multiomics data integration further unveiled positively regulated genes in OIR Müller cells/astrocytes associated with axon development and neurotransmitter transmission. This study significantly advances our understanding of the intricate cellular and molecular mechanisms underlying retinal neovascularization, emphasizing the pivotal role of Müller cells. The identified pathways provide valuable insights into potential therapeutic targets for PDR, offering promising directions for further research and clinical interventions.

Decoding the epigenetics and chromatin loop dynamics of androgen receptor-mediated transcription

Androgen receptor (AR)-mediated transcription plays a critical role in development and prostate cancer growth. AR drives gene expression by binding to thousands of cis-regulatory elements (CRE) that loop to hundreds of target promoters. With multiple CREs interacting with a single promoter, it remains unclear how individual AR bound CREs contribute to gene expression. To characterize the involvement of these CREs, we investigate the AR-driven epigenetic and chromosomal chromatin looping changes by generating a kinetic multi-omic dataset comprised of steady-state mRNA, chromatin accessibility, transcription factor binding, histone modifications, chromatin looping, and nascent RNA. Using an integrated regulatory network, we find that AR binding induces sequential changes in the epigenetic features at CREs, independent of gene expression. Further, we show that binding of AR does not result in a substantial rewiring of chromatin loops, but instead increases the contact frequency of pre-existing loops to target promoters. Our results show that gene expression strongly correlates to the changes in contact frequency. We then propose and experimentally validate an unbalanced multi-enhancer model where the impact on gene expression of AR-bound enhancers is heterogeneous, and is proportional to their contact frequency with target gene promoters. Overall, these findings provide insights into AR-mediated gene expression upon acute androgen simulation and develop a mechanistic framework to investigate nuclear receptor mediated perturbations.

Multimodal AI/ML for discovering novel biomarkers and predicting disease using multi-omics profiles of patients with cardiovascular diseases

Cardiovascular diseases (CVDs) are complex, multifactorial conditions that require personalized assessment and treatment. Advancements in multi-omics technologies, namely RNA sequencing and whole-genome sequencing, have provided translational researchers with a comprehensive view of the human genome. The efficient synthesis and analysis of this data through integrated approach that characterizes genetic variants alongside expression patterns linked to emerging phenotypes, can reveal novel biomarkers and enable the segmentation of patient populations based on personalized risk factors. In this study, we present a cutting-edge methodology rooted in the integration of traditional bioinformatics, classical statistics, and multimodal machine learning techniques. Our approach has the potential to uncover the intricate mechanisms underlying CVD, enabling patient-specific risk and response profiling. We sourced transcriptomic expression data and single nucleotide polymorphisms (SNPs) from both CVD patients and healthy controls. By integrating these multi-omics datasets with clinical demographic information, we generated patient-specific profiles. Utilizing a robust feature selection approach, we identified a signature of 27 transcriptomic features and SNPs that are effective predictors of CVD. Differential expression analysis, combined with minimum redundancy maximum relevance feature selection, highlighted biomarkers that explain the disease phenotype. This approach prioritizes both biological relevance and efficiency in machine learning. We employed Combination Annotation Dependent Depletion scores and allele frequencies to identify variants with pathogenic characteristics in CVD patients. Classification models trained on this signature demonstrated high-accuracy predictions for CVD. The best performing of these models was an XGBoost classifier optimized via Bayesian hyperparameter tuning, which was able to correctly classify all patients in our test dataset. Using SHapley Additive exPlanations, we created risk assessments for patients, offering further contextualization of these predictions in a clinical setting. Across the cohort, RPL36AP37 and HBA1 were scored as the most important biomarkers for predicting CVDs. A comprehensive literature review revealed that a substantial portion of the diagnostic biomarkers identified have previously been associated with CVD. The framework we propose in this study is unbiased and generalizable to other diseases and disorders.

Priority-Elastic net for binary disease outcome prediction based on multi-omics data

BackgroundHigh-dimensional omics data integration has emerged as a prominent avenue within the healthcare industry, presenting substantial potential to improve predictive models. However, the data integration process faces several challenges, including data heterogeneity, priority sequence in which data blocks are prioritized for rendering predictive information contained in multiple blocks, assessing the flow of information from one omics level to the other and multicollinearity.MethodsWe propose the Priority-Elastic net algorithm, a hierarchical regression method extending Priority-Lasso for the binary logistic regression model by incorporating a priority order for blocks of variables while fitting Elastic-net models sequentially for each block. The fitted values from each step are then used as an offset in the subsequent step. Additionally, we considered the adaptive elastic-net penalty within our priority framework to compare the results.ResultsThe Priority-Elastic net and Priority-Adaptive Elastic net algorithms were evaluated on a brain tumor dataset available from The Cancer Genome Atlas (TCGA), accounting for transcriptomics, proteomics, and clinical information measured over two glioma types: Lower-grade glioma (LGG) and glioblastoma (GBM).ConclusionOur findings suggest that the Priority-Elastic net is a highly advantageous choice for a wide range of applications. It offers moderate computational complexity, flexibility in integrating prior knowledge while introducing a hierarchical modeling perspective, and, importantly, improved stability and accuracy in predictions, making it superior to the other methods discussed. This evolution marks a significant step forward in predictive modeling, offering a sophisticated tool for navigating the complexities of multi-omics datasets in pursuit of precision medicine’s ultimate goal: personalized treatment optimization based on a comprehensive array of patient-specific data. This framework can be generalized to time-to-event, Cox proportional hazards regression and multicategorical outcomes. A practical implementation of this method is available upon request in R script, complete with an example to facilitate its application.

Evaluating the performance of multi-omics integration: a thyroid toxicity case study.

Multi-omics data integration has been repeatedly discussed as the way forward to more comprehensively cover the molecular responses of cells or organisms to chemical exposure in systems toxicology and regulatory risk assessment. In Canzler etal. (Arch Toxicol 94(2):371-388. https://doi.org/10.1007/s00204-020-02656-y ), we reviewed the state of the art in applying multi-omics approaches in toxicological research and chemical risk assessment. We developed best practices for the experimental design of multi-omics studies, omics data acquisition, and subsequent omics data integration. We found that multi-omics data sets for toxicological research questions were generally rare, with no datasets comprising more than two omics layers adhering to these best practices. Due to these limitations, we could not fully assess the benefits of different data integration approaches or quantitatively evaluate the contribution of various omics layers for toxicological research questions. Here, we report on a multi-omics study on thyroid toxicity that we conducted in compliance with these best practices. We induced direct and indirect thyroid toxicity through Propylthiouracil (PTU) and Phenytoin, respectively, in a 28-day plus 14-day recovery oral rat toxicity study. We collected clinical and histopathological data and six omics layers, including the long and short transcriptome, proteome, phosphoproteome, and metabolome from plasma, thyroid, and liver. We demonstrate that the multi-omics approach is superior to single-omics in detecting responses at the regulatory pathway level. We also show how combining omics data with clinical and histopathological parameters facilitates the interpretation of the data. Furthermore, we illustrate how multi-omics integration can hint at the involvement of non-coding RNAs in post-transcriptional regulation. Also, we show that multi-omics facilitates grouping, and we assess how much information individual and combinations of omics layers contribute to this approach.

A Low Transformed Tubal Rank Tensor Model Using a Spatial-Tubal Constraint for Sample Clustering with Cancer Multi-omics Data

Background: Since each dimension of a tensor can store different types of genomics data, compared to matrix methods, utilizing tensor structure can provide a deeper understanding of multi-dimensional data while also facilitating the discovery of more useful information related to cancer. However, in reality, there are issues such as insufficient utilization of prior knowledge in multi-omics data and limitations in the recovery of low-tubal-rank tensors. Therefore, the method proposed in this article was developed. Objective: In this paper, we proposed a low transformed tubal rank tensor model (LTTRT) using a spatial-tubal constraint to accurately partition different types of cancer samples and provide reliable theoretical support for the identification, diagnosis, and treatment of cancer. Method: In the LTTRT method, the transformed tensor nuclear norm based on the transformed tensor singular value decomposition is characterized by the low-rank tensor, which can explore the global low-rank property of the tensor, resolving the challenge of the tensor nuclear norm-based method not achieving the lowest tubal rank. Additionally, the introduction of weighted total variation regularization is conducive to extracting more information from sequencing data in both spatial and tubal dimensions, exploring cross-correlation features of multiple genomic data, and addressing the problem of overlooking prior knowledge from various perspectives. In addition, the L1-norm is used to improve sparsity. A symmetric Gauss‒Seidel-based alternating direction method of multipliers (sGS-ADMM) is used to update the LTTRT model iteratively. Results: The experiments of sample clustering on multiple integrated cancer multi-omics datasets show that the proposed LTTRT method is better than existing methods. Experimental results validate the effectiveness of LTTRT in accurately partitioning different types of cancer samples. result: The experiments of sample clustering on multiple integrated cancer multi-omics datasets show that the proposed LTTRT method is better than existing methods. Experimental results validate the effectiveness of LTTRT in accurately partitioning different types of cancer samples. Conclusion: The LTTRT method achieves precise segmentation of different types of cancer samples.

FOXN3-AS1: A Candidate Prognostic Marker and Epigenetic Target with Immunotherapeutic Implications in Acute Myeloid Leukemia.

We focused on the FOXN3 gene and selected its antisense transcripts (FOXN3-AS1) to investigate its potential involvement in acute myeloid leukemia (AML). Several integrated multi-omics datasets have expanded the horizons of the cancer landscape. With the emergence of new high-throughput technologies, a large number of non-coding RNAs have been confirmed to be involved in the pathogenesis of different types of hematological malignancies. We conducted experimental validation using quantitative polymerase chain reaction (qPCR) with bone marrow specimens from AML patients. Then, Kaplan-Meier (KM) and Receiver Operating Characteristic (ROC) curves were used to substantiate the prognostic association between FOXN3-AS1 and AML patients within the TCGA database. Correlation between FOXN3-AS1 expression and gene mutation, immune, and immune function using Spearman correlation analysis. To explore the physical and functional interaction between FOXN3-AS1 and the DNMT1 protein, we utilized the RPISeq web tool from Iowa State University. Subsequently, we performed qPCR experiments to test the effect of 5AzaC (DNMT1 inhibitor) on FOXN3-AS1 expression AML cell lines (THP1 and OCI-AML3). We leveraged the "OncoPredict" R package in conjunction with the Genomics of Drug Sensitivity (GDSC) database to predict drug response in AML patients expressing FOXN3-AS1. We observed a significant upregulation of FOXN3-AS1 expression in AML patients compared to healthy controls using clinical samples. The TCGA database revealed an association between high FOXN3-AS1 expression and adverse prognosis. In our subsequent analysis, genes with poor prognostic implications in AML patients were exclusively identified in the FOXN3-AS1 high-expression group, further corroborating this relationship. AML patients with higher FOXN3-AS1 expression levels may respond less optimally to immunotherapy than patients with lower levels. Besides, we computationally predicted the interaction of FOXN3- AS1 and DNMT1 protein and experimentally confirmed that DNMT1i (GSK-3484862) affects the expression level of FOXN3-AS1. We also found that the chemotherapy drugs (5-Fluorouralic, Cisplatin, Dactolisib, Sapitinib, Temozolomide, Ulixertinib, Vinorelbine, Ruxolitinib, Osimertinib and Cisplatin) showed favorable responses in AML patients with high FOXN3-AS1 expression levels. Our candidate approach identifies FOXN3-AS1 as a prognostic indicator of survival in AML with a potential immune-related role. The preliminary observations we made on FOXN3-AS1/DNMT1 crosstalk warrant more in-depth invested immunotherapeutic approaches in AML.

Transcriptome-Wide Profiling of Nascent RNA in Neurons with Enriched H3K27ac Signal Elevates eRNA Identification Efficiency.

Growing evidence suggests that activity-dependent gene expression is crucial for neuronal plasticity and behavioral experience. Enhancer RNAs (eRNAs), a class of long noncoding RNAs, play a key role in these processes. However, eRNAs are highly dynamic and are often present at lower levels than their corresponding mRNAs, making them difficult to detect using total RNA-seq techniques. Nascent RNA sequencing, which separates nascent RNAs from the steady-state RNA population, has been shown to increase the ability to detect activity-induced eRNAs with a higher signal-to-noise ratio. However, there is a lack of bioinformatic tools or pipelines for detecting eRNAs utilizing nascent RNA-seq and other multiomics data sets. In this study, we addressed this gap by developing a novel bioinformatic framework, e-finder, for finding eRNAs and have made it available to the scientific community. Additionally, we reanalyzed our previous nascent RNA sequencing data and compared them with total RNA-seq data to identify activity-regulated RNAs in neuronal cell populations. Using H3K27 acetylome data, we characterized activity-dependent eRNAs that drive the transcriptional activity of the target genes. Our analysis identified a subset of eRNAs involved in mediating synapse organization, which showed increased activity-dependent transcription after the potassium chloride stimulation. Notably, our data suggest that nascent RNA-seq with an enriched H3K27ac signal exhibits high resolution to identify potential eRNAs in response to membrane depolarization. Our findings uncover the role of the eRNA-mediated gene activation network in neuronal systems, providing new insights into the molecular processes characterizing neurological diseases.

Integrative multi-omics analysis reveals genetic and heterotic contributions to male fertility and yield in potato.

The genetic analysis of potato is hampered by the complexity of tetrasomic inheritance. An ongoing effort aims to transform the clonally propagated tetraploid potato into a seed-propagated diploid crop, which would make genetic analyses much easier owing to disomic inheritance. Here, we construct and report the large-scale genetic and heterotic characteristics of a diploid F2 potato population derived from the cross of two highly homozygous inbred lines. We investigate 20,382 traits generated from multi-omics dataset and identify 25,770 quantitative trait loci (QTLs). Coupled with gene expression data, we construct a systems-genetics network for gene discovery in potatoes. Importantly, we explore the genetic basis of heterosis in this population, especially for yield and male fertility heterosis. We find that positive heterotic effects of yield-related QTLs and negative heterotic effects of metabolite QTLs (mQTLs) contribute to yield heterosis. Additionally, we identify a PME gene with a dominance heterotic effect that plays an important role in male fertility heterosis. This study provides genetic resources for the potato community and will facilitate the application of heterosis in diploid potato breeding.

Interplay between Lipids and Complement Proteins—How Multiomics Data Integration Can Help Unravel Age-related Macular Degeneration Pathophysiology: A Proof-of-concept Study

RationaleWe hypothesized that the integration of multi-omics data sets can help unravel some complex aspects of AMD pathophysiology. Our objectives were to identify correlation patterns between genetic variants, complement proteins and lipids using a multi-omics data integration approach and to determine how these interconnections affect AMD. Materials and methodsThe analyses were performed in a subset of the Singapore Indian Eye Study (SINDI) study. We randomly selected 155 AMD cases and age- and sex-matched them with 155 controls. We measured in serum samples 35 complement proteins and 66 lipids, and used nine genetic variants. Firstly, a multi-omics data integration method was used to identify correlation patterns between the omics data. Then, we tested possible interactions between the lipids and complement proteins using logistic regression models. ResultsAmong the 155 AMD cases, 93 (60.0%) had early and 62 (40.0%) intermediate AMD. Firstly, we identified two clusters between complement proteins and lipids involving (1) MASP1 and several different HDL particles, and (2) FHR1, CPN2 and CO8G, and sphingomyelin and different cholesterol. Secondly, we identified one interaction between C1R and sphingomyelin with an odds of AMD 2 times higher for individuals with low levels of sphingomyelin and C1R (OR=2.13 [1.09, 4.17]). ConclusionWe report here, using a cutting-edge multi-omics integration approach, the complex interconnections between genetic, metabolomics and proteomic data. This method permitted us to obtain a holistic picture and identify multi-omics signature of AMD pathophysiology. These results advocate for a personalized therapeutic approach that accounts for multiple pathways. However, these results need to be validated in larger studies with different ethnic groups.

ScDRMAE: integrating masked autoencoder with residual attention networks to leverage omics feature dependencies for accurate cell clustering.

Cell clustering is foundational for analyzing the heterogeneity of biological tissues using single-cell sequencing data. With the maturation of single-cell multi-omics sequencing technologies, we can integrate multiple omics data to perform cell clustering, thereby overcoming the limitations of insufficient information from single omics data. Existing methods for cell clustering often only consider the differences in data patterns during the analysis of multi-omics data, but the dependencies between omics features of different cell types also significantly influence cell clustering. Moreover, the high dropout rates in scRNA-seq and scATAC-seq data can impact the performance of cell clustering. We propose a cell clustering model based on a masked autoencoder, scDRMAE. Utilizing a masking mechanism, scDRMAE effectively learns the relationships between different features and imputes false zeros caused by dropout events. To differentiate the importance of various omics data in cell clustering, we dynamically adjust the weights of different omics data through an attention mechanism. Finally, we use the K-means algorithm for cluster analysis of the fused multi-omics data. On commonly used sets of 15 multi-omics datasets, our method demonstrates superior cell clustering performance on multiple metrics compared to other computational methods. In addition, when datasets exhibit varying degrees of dropout noise, our method shows better performance and stronger stability on multiple metrics compared to other methods. Moreover, by analyzing the cell clusters classified by scDRMAE, we identified several biologically significant biomarkers that have been validated, further confirming the effectiveness of scDRMAE in cell clustering from a biological perspective.

Prioritizing disease-related rare variants by integrating gene expression data.

Rare variants, comprising the vast majority of human genetic variations, are likely to have more deleterious impact in the context of human diseases compared to common variants. Here we present carrier statistic, a statistical framework to prioritize disease-related rare variants by integrating gene expression data. By quantifying the impact of rare variants on gene expression, carrier statistic can prioritize those rare variants that have large functional consequence in the patients. Through simulation studies and analyzing real multi-omics dataset, we demonstrated that carrier statistic is applicable in studies with limited sample size (a few hundreds) and achieves substantially higher sensitivity than existing rare variants association methods. Application to Alzheimer's disease reveals 16 rare variants within 15 genes with extreme carrier statistics. We also found strong excess of rare variants among the top prioritized genes in patients compared to that in healthy individuals. The carrier statistic method can be applied to various rare variant types and is adaptable to other omics data modalities, offering a powerful tool for investigating the molecular mechanisms underlying complex diseases.

Harnessing Multi-Omics Strategies and Bioinformatics Innovations for Advancing Soybean Improvement: A Comprehensive Review.

Soybean improvement has entered a new era with the advent of multi-omics strategies and bioinformatics innovations, enabling more precise and efficient breeding practices. This comprehensive review examines the application of multi-omics approaches in soybean-encompassing genomics, transcriptomics, proteomics, metabolomics, epigenomics, and phenomics. We first explore pre-breeding and genomic selection as tools that have laid the groundwork for advanced trait improvement. Subsequently, we dig into the specific contributions of each -omics field, highlighting how bioinformatics tools and resources have facilitated the generation and integration of multifaceted data. The review emphasizes the power of integrating multi-omics datasets to elucidate complex traits and drive the development of superior soybean cultivars. Emerging trends, including novel computational techniques and high-throughput technologies, are discussed in the context of their potential to revolutionize soybean breeding. Finally, we address the challenges associated with multi-omics integration and propose future directions to overcome these hurdles, aiming to accelerate the pace of soybean improvement. This review serves as a crucial resource for researchers and breeders seeking to leverage multi-omics strategies for enhanced soybean productivity and resilience.

Systemic dysregulation and molecular insights into poor influenza vaccine response in the aging population.

Vaccination-induced protection against influenza is greatly diminished and increasingly heterogeneous with age. We investigated longitudinally (up to five time points) a cohort of 234 vaccinated >65-year-old vaccinees with adjuvanted vaccine FluAd across two independent seasons. System-level analyses of multiomics datasets measuring six modalities and serological data revealed that poor responders lacked time-dependent changes in response to vaccination as observed in responders, suggestive of systemic dysregulation in poor responders. Multiomics integration revealed key molecules and their likely role in vaccination response. High prevaccination plasma interleukin-15 (IL-15) concentrations negatively associated with antibody production, further supported by experimental validation in mice revealing an IL-15-driven natural killer cell axis explaining the suppressive role in vaccine-induced antibody production as observed in poor responders. We propose a subset of long-chain fatty acids as modulators of persistent inflammation in poor responders. Our findings provide a potential link between low-grade chronic inflammation and poor vaccination response and open avenues for possible pharmacological interventions to enhance vaccine responses.

An explainable graph neural network approach for effectively integrating multi-omics with prior knowledge to identify biomarkers from interacting biological domains.

The rapid growth of multi-omics datasets, in addition to the wealth of existing biological prior knowledge, necessitates the development of effective methods for their integration. Such methods are essential for building predictive models and identifying disease-related molecular markers. We propose a framework for supervised integration of multi-omics data with biological priors represented as knowledge graphs. Our framework leverages graph neural networks (GNNs) to model the relationships among features from high-dimensional 'omics data and set transformers to integrate low-dimensional representations of 'omics features. Furthermore, our framework incorporates explainability methods to elucidate important biomarkers and extract interaction relationships between biological quantities of interest. We demonstrate the effectiveness of our approach by applying it to Alzheimer's disease (AD) multi-omics data from the ROSMAP cohort, showing that the integration of transcriptomics and proteomics data with AD biological domain network priors improves the prediction accuracy of AD status and highlights functional AD biomarkers.

G-Protein Signaling in Alzheimer's Disease: Spatial Expression Validation of Semi-supervised Deep Learning-Based Computational Framework.

Systemic study of pathogenic pathways and interrelationships underlying genes associated with Alzheimer's disease (AD) facilitates the identification of new targets for effective treatments. Recently available large-scale multiomics datasets provide opportunities to use computational approaches for such studies. Here, we devised a novel disease gene identification (digID) computational framework that consists of a semi-supervised deep learning classifier to predict AD-associated genes and a protein-protein interaction (PPI) network-based analysis to prioritize the importance of these predicted genes in AD. digID predicted 1,529 AD-associated genes and revealed potentially new AD molecular mechanisms and therapeutic targets including GNAI1 and GNB1, two G-protein subunits that regulate cell signaling, and KNG1, an upstream modulator of CDC42 small G-protein signaling and mediator of inflammation and candidate coregulator of amyloid precursor protein (APP). Analysis of mRNA expression validated their dysregulation in AD brains but further revealed the significant spatial patterns in different brain regions as well as among different subregions of the frontal cortex and hippocampi. Super-resolution STochastic Optical Reconstruction Microscopy (STORM) further demonstrated their subcellular colocalization and molecular interactions with APP in a transgenic mouse model of both sexes with AD-like mutations. These studies support the predictions made by digID while highlighting the importance of concurrent biological validation of computationally identified gene clusters as potential new AD therapeutic targets.