Single-cell transcriptomes are established by transcription factors (TFs), which determine a cell's gene-expression complement. Post-transcriptional aspects of single-cell transcriptomes, and the RNA binding proteins (RBPs) responsible, are more technically challenging to determine, and combinatorial TF-RBP coordination of single-cell transcriptomes remains unexplored. We used fluorescent reporters to visualize alternative splicing in single Caenorhabditis elegans neurons, identifying complex splicing patterns in the neuronal kinase
sad-1. Most neurons express both isoforms, but the ALM mechanosensory neuron expresses only the exon-included isoform, while its developmental sister cell the BDU neuron expresses only the exon-skipped isoform. A cascade of three cell-specific TFs (UNC-86, MEC-3, and ALR-1) and two RBPs (MEC-8 and MBL-1) are combinatorially required for
sad-1 exon inclusion. Mechanistically, TFs combinatorially ensure RBP expression, which interact with
sad-1 pre-mRNA. Thus a combinatorial TF-RBP code controls single-neuron
sad-1 splicing. Additionally, we find "phenotypic convergence," previously observed for TFs, also applies to RBPs: different RBP combinations generate similar splicing outcomes in different neurons.