Integration of Imaging Genomics Data for the Study of Alzheimer's Disease Using Joint-Connectivity-Based Sparse Nonnegative Matrix Factorization

Singular Value Decomposition-Driven Non-negative Matrix Factorization with Application to Identify the Association Patterns of Sarcoma Recurrence

Sarcomas are malignant tumors from mesenchymal tissue and are characterized by their complexity and diversity. The high recurrence rate making it important to understand the mechanisms behind their recurrence and to develop personalized treatments and drugs. However, previous studies on the association patterns of multi-modal data on sarcoma recurrence have overlooked the fact that genes do not act independently, but rather function within signaling pathways. Therefore, this study collected 290 whole solid images, 869 gene and 1387 pathway data of over 260 sarcoma samples from UCSC and TCGA to identify the association patterns of gene-pathway-cell related to sarcoma recurrences. Meanwhile, considering that most multi-modal data fusion methods based on the joint non-negative matrix factorization (NMF) model led to poor experimental repeatability due to random initialization of factorization parameters, the study proposed the singular value decomposition (SVD)-driven joint NMF model by applying the SVD method to calculate initialized weight and coefficient matrices to achieve the reproducibility of the results. The results of the experimental comparison indicated that the SVD algorithm enhances the performance of the joint NMF algorithm. Furthermore, the representative module indicated a significant relationship between genes in pathways and image features. Multi-level analysis provided valuable insights into the connections between biological processes, cellular features, and sarcoma recurrence. In addition, potential biomarkers were uncovered, while various mechanisms of sarcoma recurrence were identified from an imaging genetic perspective. Overall, the SVD-NMF model affords a novel perspective on combining multi-omics data to explore the association related to sarcoma recurrence.

Deep self-reconstruction driven joint nonnegative matrix factorization model for identifying multiple genomic imaging associations in complex diseases

OBJECTIVE

Comprehensive analysis of histopathology images and transcriptomics data enables the identification of candidate biomarkers and multimodal association patterns. Most existing multimodal data association studies are derived from extensions of the joint nonnegative matrix factorization model for identifying complex data associations, which can make full use of clinical prior information. However, the raw data were usually taken as the input without considering the underlying complex multi-subspace structure, influencing the subsequent integration analysis results.

METHODS

This study proposed a deep-self reconstructed joint nonnegative matrix factorization (DSRJNMF) model to use self-expressive properties to reconstruct the raw data to characterize the similarity structure associated with clinical labels. Then, the sparsity, orthogonality, and regularization constraints constructed from prior information are added to the DSRJNMF model to determine the sparse set of biologically relevant features across modalities.

RESULTS

The algorithm has been applied to identify the imaging genetic association of triple negative breast cancer (TNBC). Multilevel experimental results demonstrate that the proposed algorithm better estimates potential associations between pathological image features and miRNA-gene and identifies consistent multimodal imaging genetic biomarkers to guide the interpretation of TNBC.

CONCLUSION

The propose method provides a novel idea of data association analysis oriented to complex diseases.

Integrating multi-omics data of childhood asthma using a deep association model

Childhood asthma is one of the most common respiratory diseases with rising mortality and morbidity. The multi-omics data is providing a new chance to explore collaborative biomarkers and corresponding diagnostic models of childhood asthma. To capture the nonlinear association of multi-omics data and improve interpretability of diagnostic model, we proposed a novel deep association model (DAM) and corresponding efficient analysis framework. First, the Deep Subspace Reconstruction was used to fuse the omics data and diagnostic information, thereby correcting the distribution of the original omics data and reducing the influence of unnecessary data noises. Second, the Joint Deep Semi-Negative Matrix Factorization was applied to identify different latent sample patterns and extract biomarkers from different omics data levels. Third, our newly proposed Deep Orthogonal Canonical Correlation Analysis can rank features in the collaborative module, which are able to construct the diagnostic model considering nonlinear correlation between different omics data levels. Using DAM, we deeply analyzed the transcriptome and methylation data of childhood asthma. The effectiveness of DAM is verified from the perspectives of algorithm performance and biological significance on the independent test dataset, by ablation experiment and comparison with many baseline methods from clinical and biological studies. The DAM-induced diagnostic model can achieve a prediction AUC of 0.912, which is higher than that of many other alternative methods. Meanwhile, relevant pathways and biomarkers of childhood asthma are also recognized to be collectively altered on the gene expression and methylation levels. As an interpretable machine learning approach, DAM simultaneously considers the non-linear associations among samples and those among biological features, which should help explore interpretative biomarker candidates and efficient diagnostic models from multi-omics data analysis for human complex diseases.

Integrating Multi-omics Data for Alzheimer's Disease to Explore Its Biomarkers Via the Hypergraph-Regularized Joint Deep Semi-Non-Negative Matrix Factorization Algorithm

Alzheimer's disease (AD) is a progressive and irreversible neurodegenerative disorder. Its etiology may be associated with genetic, environmental, and lifestyle factors. With the advancement of technology, the integration of genomics, transcriptomics, and imaging data related to AD allows simultaneous exploration of molecular information at different levels and their interaction within the organism. This paper proposes a hypergraph-regularized joint deep semi-non-negative matrix factorization (HR-JDSNMF) algorithm to integrate positron emission tomography (PET), single-nucleotide polymorphism (SNP), and gene expression data for AD. The method employs matrix factorization techniques to nonlinearly decompose the original data at multiple layers, extracting deep features from different omics data, and utilizes hypergraph mining to uncover high-order correlations among the three types of data. Experimental results demonstrate that this approach outperforms several matrix factorization-based algorithms and effectively identifies multi-omics biomarkers for AD. Additionally, single-cell RNA sequencing (scRNA-seq) data for AD were collected, and genes within significant modules were used to categorize different types of cell clusters into high and low-risk cell groups. Finally, the study extensively explores the differences in differentiation and communication between these two cell types. The multi-omics biomarkers unearthed in this study can serve as valuable references for the clinical diagnosis and drug target discovery for AD. The realization of the algorithm in this paper code is available at https://github.com/ShubingKong/HR-JDSNMF .

Estrogen Receptor Genes, Cognitive Decline, and Alzheimer Disease

BACKGROUND AND OBJECTIVES

Lifetime risk of Alzheimer disease (AD) dementia is twofold higher in women compared with men, and low estrogen levels in postmenopause have been suggested as a possible contributor. We examined 3 ER (GPER1, ER2, and ER1) variants in association with AD traits as an indirect method to test the association between estrogen and AD in women. Although the study focus was on women, in a comparison, we separately examined ER molecular variants in men.

METHODS

Participants were followed for an average of 10 years in one of the 2 longitudinal clinical pathologic studies of aging. Global cognition was assessed using a composite score derived from 19 neuropsychological tests' scores. Postmortem pathologic assessment included examination of 3 AD (amyloid-β and tau tangles determined by immunohistochemistry, and a global AD pathology score derived from diffuse and neurotic plaques and neurofibrillary tangle count) and 8 non-AD pathology indices. ER molecular genomic variants included genotyping and examining ER DNA methylation and RNA expression in brain regions including the dorsolateral prefrontal cortex (DLPFC) that are major players in cognition and often have AD pathology.

RESULTS

The mean age of women (N = 1711) at baseline was 78.0 (SD = 7.7) years. In women, GPER1 molecular variants had the most consistent associations with AD traits. GPER1 DNA methylation was associated with cognitive decline, tau tangle density, and global AD pathology score. GPER1 RNA expression in DLPFC was related to cognitive decline and tau tangle density. Other associations included associations of ER2 and ER1 sequence variants and DNA methylation with cognition. RNA expressions in DLPFC of genes involved in signaling mechanisms of activated ERs were also associated with cognitive decline and tau tangle density in women. In men (N = 651, average age at baseline: 77.4 [SD = 7.3]), there were less robust associations between ER molecular genomic variants and AD cognitive and pathologic traits. No consistent association was seen between ER molecular genomic variations and non-AD pathologies in either of the sexes.

DISCUSSION

ER DNA methylation and RNA expression, and to some extent ER polymorphisms, were associated with AD cognitive and pathologic traits in women, and to a lesser extent in men.

A multi-task SCCA method for brain imaging genetics and its application in neurodegenerative diseases

BACKGROUND AND OBJECTIVES

In brain imaging genetics, multi-task sparse canonical correlation analysis (MTSCCA) is effective to study the bi-multivariate associations between genetic variations such as single nucleotide polymorphisms (SNPs) and multi-modal imaging quantitative traits (QTs). However, most existing MTSCCA methods are neither supervised nor capable of distinguishing the shared patterns of multi-modal imaging QTs from the specific patterns.

METHODS

A new diagnosis-guided MTSCCA (DDG-MTSCCA) with parameter decomposition and graph-guided pairwise group lasso penalty was proposed. Specifically, the multi-tasking modeling paradigm enables us to comprehensively identify risk genetic loci by jointly incorporating multi-modal imaging QTs. The regression sub-task was raised to guide the selection of diagnosis-related imaging QTs. To reveal the diverse genetic mechanisms, the parameter decomposition and different constraints were utilized to facilitate the identification of modality-consistent and -specific genotypic variations. Besides, a network constraint was added to find out meaningful brain networks. The proposed method was applied to synthetic data and two real neuroimaging data sets respectively from Alzheimer's disease neuroimaging initiative (ADNI) and Parkinson's progression marker initiative (PPMI) databases.

RESULTS

Compared with the competitive methods, the proposed method exhibited higher or comparable canonical correlation coefficients (CCCs) and better feature selection results. In particular, in the simulation study, DDG-MTSCCA showed the best anti-noise ability and achieved the highest average hit rate, about 25% higher than MTSCCA. On the real data of Alzheimer's disease (AD) and Parkinson's disease (PD), our method obtained the highest average testing CCCs, about 40% ∼ 50% higher than MTSCCA. Especially, our method could select more comprehensive feature subsets, and the top five SNPs and imaging QTs were all disease-related. The ablation experimental results also demonstrated the significance of each component in the model, i.e., the diagnosis guidance, parameter decomposition, and network constraint.

CONCLUSIONS

These results on simulated data, ADNI and PPMI cohorts suggested the effectiveness and generalizability of our method in identifying meaningful disease-related markers. DDG-MTSCCA could be a powerful tool in brain imaging genetics, worthy of in-depth study.

Application of orthogonal sparse joint non-negative matrix factorization based on connectivity in Alzheimer's disease research

Based on the mining of micro- and macro-relationships of genetic variation and brain imaging data, imaging genetics has been widely applied in the early diagnosis of Alzheimer's disease (AD). However, effective integration of prior knowledge remains a barrier to determining the biological mechanism of AD. This paper proposes a new connectivity-based orthogonal sparse joint non-negative matrix factorization (OSJNMF-C) method based on integrating the structural magnetic resonance image, single nucleotide polymorphism and gene expression data of AD patients; the correlation information, sparseness, orthogonal constraint and brain connectivity information between the brain image data and genetic data are designed as constraints in the proposed algorithm, which efficiently improved the accuracy and convergence through multiple iterative experiments. Compared with the competitive algorithm, OSJNMF-C has significantly smaller related errors and objective function values than the competitive algorithm, showing its good anti-noise performance. From the biological point of view, we have identified some biomarkers and statistically significant relationship pairs of AD/mild cognitive impairment (MCI), such as rs75277622 and BCL7A, which may affect the function and structure of multiple brain regions. These findings will promote the prediction of AD/MCI.

Integration of RNA molecules data with prior-knowledge driven Joint Deep Semi-Negative Matrix Factorization for heart failure study

Heart failure (HF) is the main manifestation of cardiovascular disease. Recent studies have shown that various RNA molecules and their complex connections play an essential role in HF’s pathogenesis and pathological progression. This paper aims to mine key RNA molecules associated with HF. We proposed a Prior-knowledge Driven Joint Deep Semi-Negative Matrix Factorization (PD-JDSNMF) model that uses a hierarchical nonlinear feature extraction method that integrates three types of data: mRNA, lncRNA, and miRNA. The PPI information is added to the model as prior knowledge, and the Laplacian constraint is used to help the model resist the noise in the genetic data. We used the PD-JDSNMF algorithm to identify significant co-expression modules. The elements in the module are then subjected to bioinformatics analysis and algorithm performance analysis. The results show that the PD-JDSNMF algorithm can robustly select biomarkers associated with HF. Finally, we built a heart failure diagnostic model based on multiple classifiers and using the Top 13 genes in the significant module, the AUC of the internal test set was up to 0.8714, and the AUC of the external validation set was up to 0.8329, which further confirmed the effectiveness of the PD-JDSNMF algorithm.