We generated subcluster-level profiles by summing the integer counts of each cell in each subcluster together to form 565 distinct profiles (Figure S7A,B). To compare pairwise correlations across nAChR subunits, all 565 profiles were selected and normalized to 100K UMIs. Pearson correlations were calculated across all genes using linear expression data and grouped via hierarchical clustering using the heatmap.2 function of the gplots R package. Note, Chrna3/Chrnb4/Chrna5 are adjacently located on chromosome 9; Chrna4/Chrnb2 genes are on chromosome 2 and 3, respectively. For voltage gated Na and K alpha subunit analysis, only the 323 profiles corresponding to neurons were used and the gene set was thresholded to include only those genes with > 1 out of 100K transcripts. Scatterplots of profile expression were displayed in log10 space.
