Multi-omics integrative analysis is pivotal for elucidating complex molecular mechanisms and biological processes, yet remains challenging due to the high dimensionality and heterogeneity of multi-omics data. Here we describe MIA, a machine learning framework for multi-omics integrative analysis. Unlike existing algorithms that rely on two-dimensional representations, MIA employs a three-dimensional tensor representation coupled with tensor decomposition, Fuzzy C-Means, and an enhanced random forest model to jointly enable accurate sample stratification and feature discovery. Benchmarking on simulated datasets demonstrates that MIA achieves higher accuracy in both clustering and feature identification by comparison with extant algorithms. Application to TCGA datasets further shows the ability of MIA to stratify samples and identify features associated with significantly clinical outcomes. Notably, as applied to glioblastoma, MIA is capable to detect three previously uncharacterized subtypes with distinct prognostic profiles and uncover critical feature genes linked to subtyping. These feature genes are associated with therapeutic response and retain discriminative capacity for major GBM cellular subpopulations at the single-cell level. Collectively, these results establish MIA as a generalizable computational framework for multi-omics integrative analysis, enabling systematic molecular subtyping and significant feature discovery across complex biological systems.