The origin of life and native tissue development are dependent on the heterogeneity of pluripotent stem cells. Bone marrow mesenchymal stem cells (BMMSCs) are located in a complicated niche with variable matrix stiffnesses, resulting in divergent stem cell fates. However, how stiffness drives stem cell fate remains unknown. For this study, we performed whole-gene transcriptomics and precise untargeted metabolomics sequencing to elucidate the complex interaction network of stem cell transcriptional and metabolic signals in extracellular matrices (ECMs) with different stiffnesses, and we propose a potential mechanism involved in stem cell fate decision. In a stiff (39~45 kPa) ECM, biosynthesis of aminoacyl-tRNA was up-regulated, and increased osteogenesis was also observed. In a soft (7~10 kPa) ECM, biosynthesis of unsaturated fatty acids and deposition of glycosaminoglycans were increased, accompanied by enhanced adipogenic/chondrogenic differentiation of BMMSCs. In addition, a panel of genes responding to the stiffness of the ECM were validated in vitro, mapping out the key signaling network that regulates stem cells’ fate decisions. This finding of “stiffness-dependent manipulation of stem cell fate” provides a novel molecular biological basis for development of potential therapeutic targets within tissue engineering, from both a cellular metabolic and a biomechanical perspective.