Purpose: In this preliminary study, recently developed unsupervised nonlinear dimension reduction ͑DR͒ and data representation techniques were applied to computer-extracted breast lesion feature spaces across three separate imaging modalities: Ultrasound ͑U.S.͒ with 1126 cases, dynamic contrast enhanced magnetic resonance imaging with 356 cases, and full-field digital mammography with 245 cases. Two methods for nonlinear DR were explored: Laplacian eigenmaps ͓M. Belkin and P. Niyogi, "Laplacian eigenmaps for dimensionality reduction and data representation," Neural Comput. 15, 1373-1396 ͑2003͔͒ and t-distributed stochastic neighbor embedding ͑t-SNE͒ ͓L. van der Maaten and G. Hinton, "Visualizing data using t-SNE," J. Mach. Learn. Res. 9, 2579-2605 ͑2008͔͒. Methods: These methods attempt to map originally high dimensional feature spaces to more human interpretable lower dimensional spaces while preserving both local and global information. The properties of these methods as applied to breast computer-aided diagnosis ͑CADx͒ were evaluated in the context of malignancy classification performance as well as in the visual inspection of the sparseness within the two-dimensional and three-dimensional mappings. Classification performance was estimated by using the reduced dimension mapped feature output as input into both linear and nonlinear classifiers: Markov chain Monte Carlo based Bayesian artificial neural network ͑MCMC-BANN͒ and linear discriminant analysis. The new techniques were compared to previously developed breast CADx methodologies, including automatic relevance determination and linear stepwise ͑LSW͒ feature selection, as well as a linear DR method based on principal component analysis. Using ROC analysis and 0.632+ bootstrap validation, 95% empirical confidence intervals were computed for the each classifier's AUC performance. Results: In the large U.S. data set, sample high performance results include, AUC 0.632+ = 0.88 with 95% empirical bootstrap interval ͓0.787;0.895͔ for 13 ARD selected features and AUC 0.632+ = 0.87 with interval ͓0.817;0.906͔ for four LSW selected features compared to 4D t-SNE mapping ͑from the original 81D feature space͒ giving AUC 0.632+ = 0.90 with interval ͓0.847;0.919͔, all using the MCMC-BANN. Conclusions: Preliminary results appear to indicate capability for the new methods to match or exceed classification performance of current advanced breast lesion CADx algorithms. While not appropriate as a complete replacement of feature selection in CADx problems, DR techniques offer a complementary approach, which can aid elucidation of additional properties associated with the data. Specifically, the new techniques were shown to possess the added benefit of delivering sparse lower dimensional representations for visual interpretation, revealing intricate data structure of the feature space.