SCRAMseq用转录组读数取代定制试验
多重变异效应试验（MAVE）能一次检测大量变异，但多数版本需要为每个基因专门构建表型读数，例如细胞存活或某一种表面标志物的表达。正因如此，MAVE在癌症和常见病遗传学中的推广快于罕见病领域：罕见病确诊患者少，专门开发试验难以获得足够支持。论文的前提还包括，单一读数可能把蛋白质的完整活性压缩成一个数字。
SCRAMseq（Single Cell RNAseq Associated with MAVEs by sequencing）同时应对这两个问题。研究引入一个平衡的错义变异文库，使每个细胞最多携带一个基因拷贝；细胞扩增超过72小时后，对单细胞进行全长scRNA测序。每个细胞产生两条相互关联的记录：RNA基因分型确定它携带哪个变异，转录组则描述其表型。变异效应被汇总为由基因表达特征得出的耐受性评分，因此该试验不依赖任何预先选定的标志物。
2,265个TP63变异文库经角质形成细胞转化读出
作者聚焦TP63，这是一种驱动复层上皮发育的转录因子。ClinVar中TP63变异有35%意义未明，其中73%为错义变异；已知致病错义变异集中在DNA结合域（DBD）和 sterile alpha motif（SAM）结构域。团队对DBD第260–319位氨基酸和SAM结构域第513–568位氨基酸进行诱变，共产生2,265个变异，占这些区域所有可能替换的99.69%。
为让变异处于有生物学意义的背景中，他们优化了将人BJ-T成纤维细胞转化为诱导性角质形成细胞样细胞的过程，方法是用慢病毒载体共表达P63和KLF4。14天内转化效率超过70%，并通过流式分选ITGβ4（基底角质形成细胞的表面标志物，也是P63的靶标）测得75%。批量RNA测序发现超过1,400个差异表达基因，富集于角质形成细胞、皮肤和表皮相关条目。这一转化体系是基因特异性参照，用来评判可推广的SCRAMseq读数。
单变异试验确认中性与致病的划分
作者从每个预测类别中挑选单个变异，用独立方法检测。四个中性变异（DBD中的G293H和G310Q；SAM中的T527I和H556A）和十四个致病变异接受了KRT14定量免疫荧光、ITGβ4流式细胞术、RNA测序、非变性PAGE和荧光素酶置换报告试验分析。
区分结果一致。只有中性变异能激活KRT14并转化成纤维细胞，没有任何致病变异能做到。每个致病性SAM变异都出现聚集，聚集体丰度是单体的四倍，而DBD变异像野生型一样保持单体状态。中性变异达到85–90%的KRT14阳性细胞，与野生型P63相当，但G310Q约为60%。G310Q在转录组上也与众不同：它缺失了野生型P63诱导的整块基因，包括KRT5、KRT14和ITGβ4，并在置换试验中显示荧光素酶表达约为野生型的七倍。论文把G310Q视为活性可测降低的中性类变异，而非干净的阴性。
SCRAMseq与基因特异性试验一致率为84.2%
核心比较是SCRAMseq与MITE，后者是成熟流程：按ITGβ4分选转化细胞，并通过批量扩增子测序读取变异。作者把MITE作为功能参照而非临床金标准，报告总体一致率84.2%，其中致病-致病一致占53.8%，中性-中性一致占30.4%。在转化条件下分类一致性的AUC为0.91，在非转化条件下为0.92。
剩余的不一致被归因于两种方法测量的内容不同：MITE按ITGβ4表达分离细胞，而SCRAMseq使用广泛的转录组特征。论文还报告，在Spearman相关、基于AUC的相关和一致率上，实验评分优于包括大语言模型在内的计算机预测工具。这一比较支持功能数据的价值，但它是对预测工具的基准比较，并不证明SCRAMseq可以取代临床分类。
评分分布在全文库中区分变异类别
在所有已分类变异中，中性变异的耐受性评分分布与致病变异显著不同（Mann–Whitney U，P < 2.2e-16）。中性变异评分中位数高于零，致病变异中位数低于零。作者随后选取代表性最好的五个中性变异（L568A、W559F、Y535C、D524K、D543N）和五个致病变异（I562G、F513P、I537P、I558S、L545N），比较其评分分布，同样发现显著差异（P = 1.6e-15）。
模式并不对称。致病突变几乎完全无法激活该特征，而中性突变在较宽的评分范围内都能激活。作者把中性变异之间的这种分散与过表达时间有限以及体系的无偏设计联系起来。实际含义是，低评分是功能丧失的强信号，而高评分并不保证活性与野生型等同。
MYOD1显示该流程并非TP63专属
为检验SCRAMseq是否依赖TP63转化体系，作者构建了覆盖MYOD1（一种发育转录因子）前端区域的MITE文库，约1,300个变异。他们特意选择已知功能结构域之外的N端区域，并在与TP63非转化体系相同的非选择性过表达条件下运行。
在单细胞实验中，成功基因分型的细胞数为：P63转化组2,331个，P63非转化组1,216个，MYOD1非转化组1,836个。检测到的变异数为：P63转化组2,228个中的1,108个，P63非转化组1,056个中的719个，MYOD1非转化组1,306个中的790个。因此MYOD1实验证明了在第二个基因和没有已知功能试验的结构域中的可行性，但它展示的是适用性，而非对MYOD1变异解读的临床验证。
证据支持什么，又止步于哪里
数据支持的最强主张是方法学层面的：全长单细胞转录组可以作为MAVE的通用表型读数，并且在TP63中足够接近地重现基因特异性MITE分类，因而具有实用价值。两个基因、两种条件的设计，比单基因研究更好地支持了跨实验条件和基因的可移植性。
若干边界对解读很重要。MITE是功能参照而非临床金标准，因此84.2%的一致率衡量的是试验之间的一致性，而非诊断准确性。每次运行的变异回收都是部分的，每个实验成功基因分型的细胞数在1,216至2,331之间，因此覆盖度因条件而异。作者还指出，包括本研究在内的现有MAVE仍依赖相对简单的表型读数，中性类变异之间的分散反映了过表达窗口较短。论文自身的表述是，SCRAMseq拓宽了MAVE可应用的范围，而不是消除了对单个变异进行正交验证的需要。
论文确立的优先取舍是：放弃基因特异性选择试验，换取可指向任何致病驱动基因的转录组范围读数。对TP63而言，这一取舍的代价是与MITE约16%的分类不一致；对于没有成熟功能试验的基因，例如本文测试的MYOD1区域，则没有MITE基线可失去。读者真正应带走的实际判断阈值是这种不对称性，而不是标题式的一致率数字。
本文引用的文献
Genotype-phenotype single-cell transcriptomics for massive parallel assessment of genetic variants
SCRAMseq将全长单细胞RNA测序与MITE变异文库结合，把基因型映射到转录组表型，以AUC 0.91（转化条件）和0.92（非转化条件）区分TP63中性变异与致病变异，并与基因特异性MITE试验达到84.2%的一致率。