The immense and growing repositories of transcriptional data may contain critical insights for developing new therapies. Current approaches to mining these data largely rely on binary classifications of disease vs. control, and are not able to incorporate measures of disease severity. We report an analytical approach to integrate ordinal clinical information with transcriptomics. We apply this method to public data for a large cohort of Huntington's disease patients and controls, identifying and prioritizing phenotype-associated genes. We verify the role of a high-ranked gene in dysregulation of sphingolipid metabolism in the disease and demonstrate that inhibiting the enzyme, sphingosine-1-phosphate lyase 1 (SPL), has neuroprotective effects in Huntington's disease models. Finally, we show that one consequence of inhibiting SPL is intracellular inhibition of histone deacetylases, thus linking our observations in sphingolipid metabolism to a well-characterized Huntington's disease pathway. Our approach is easily applied to any data with ordinal clinical measurements, and may deepen our understanding of disease processes.Identifying gene subsets affecting disease phenotypes from transcriptome data is challenge. Here, the authors develop a method that combines transcriptional data with disease ordinal clinical measurements to discover a sphingolipid metabolism regulator involving in Huntington's disease progression.