Back

Single-cell profiling resolves gain- and loss-of-function mechanisms in CASR to advance mechanism-aware variant classification

Padigepati, S. R.; Stafford, D.; Kobayashi, Y.; Tan, C. A.; Kim, J.; Jiang, Y.; Ouyang, K.; Reuter, J.; Lacoste, A. M. B.

2026-09-29 genetics
10.64898/2026.09.26.754690 bioRxiv
Show abstract

Multiplexed assays of variant effect (MAVEs) and computational methods have advanced variant classification but have struggled to distinguish the molecular mechanism of variants with sufficient accuracy, limiting utility for multi-mechanism genes. Here, we coupled single-cell RNA-seq profiling with supervised machine learning to predict the mechanism of variants in CASR, a gene in which gain of function (GOF) variants cause hypocalcemia and loss of function (LOF) variants cause hypercalcemia. We engineered endogenous CASR-depleted HEK293 cells to express a single exogenous variant copy, measured their transcriptional profiles at single-cell resolution, trained a multi-class machine learning model on 96 expert-annotated variants (25 benign, 35 LOF, 36 GOF), and predicted mechanism for 157 variants of uncertain significance, of which 16 were predicted GOF and 54 LOF. The model showed high performance (macro F1 0.98, overall accuracy 0.98), and its predictions were corroborated by clinical phenotype and serum calcium measurements across more than 300,000 patients in a commercial genetic testing database. These results highlight the effectiveness of using molecular mechanism annotations with scRNA-seq to train models that distinguish molecular mechanisms of pathogenicity. The approach has immediate relevance for clinical variant classification and establishes a foundation for mechanism-guided clinical management in disorders with mechanism-dependent interventions, with potential application to other multi-mechanism genes.

Matching journals

The top 4 journals account for 50% of the predicted probability mass.

50% of probability mass above

"Similar papers" are the closest papers from that journal in the model's embedding space. They show what the match is built on, but the ranking comes mostly from a classifier over the whole training set, not from these examples alone.