Nuclei were included for clustering analysis if they passed all of the following QC thresholds:
>30% cDNA longer than 400 base pairs
>500,000 reads aligned to exonic or intronic sequence
>40% of total reads aligned
>50% unique reads
TA nucleotide ratio > 0.7
After clustering (see below), clusters were identified as outliers if more than half of nuclei co-expressed markers of inhibitory (GAD1, GAD2) and excitatory (SLC17A7) neurons or were NeuN+ but did not express the pan-neuronal marker SNAP25. Median values of QC metrics listed above were calculated for each cluster and used to compute the median and inter-quartile range (IQR) of all cluster medians. Clusters were also identified as outliers if the cluster median QC metrics deviated by more than three times the IQRs from the median of all clusters. In total, 15,928 nuclei passed QC criteria and were split into three broad classes of cells (10,708 excitatory neurons, 4,297 inhibitory neurons, and 923 non-neuronal cells) based on NeuN staining and cell class marker gene expression
Clusters were identified as donor-specific if they included fewer nuclei sampled from donors than expected by chance. For each cluster, the expected proportion of nuclei from each donor was calculated based on the laminar composition of the cluster and laminar sampling of the donor. For example, if 30% of layer 3 nuclei were sampled from a donor, then a layer 3-enriched cluster should contain approximately 30% of nuclei from this donor. In contrast, if only layer 5 were sampled from a donor, then the expected sampling from this donor for a layer 1-enriched cluster was zero. If the difference between the observed and expected sampling was greater than 50% of the number of nuclei in the cluster, then the cluster was flagged as donor-specific and excluded. In total, 325 nuclei were assigned to donor-specific or outlier clusters that contained marginal quality nuclei and were excluded from further analysis. Three donor-specific clusters came from neurosurgical donors (n=95 nuclei) and were similar to other layer 5 types reported in our analysis, but had higher expression of activity-dependent genes.
To confirm exclusion, clusters automatically flagged as outliers or donor-specific were manually inspected for expression of broad cell class marker genes, mitochondrial genes related to quality, and known activity-dependent genes.
>30% cDNA longer than 400 base pairs
>500,000 reads aligned to exonic or intronic sequence
>40% of total reads aligned
>50% unique reads
TA nucleotide ratio > 0.7
After clustering (see below), clusters were identified as outliers if more than half of nuclei co-expressed markers of inhibitory (GAD1, GAD2) and excitatory (SLC17A7) neurons or were NeuN+ but did not express the pan-neuronal marker SNAP25. Median values of QC metrics listed above were calculated for each cluster and used to compute the median and inter-quartile range (IQR) of all cluster medians. Clusters were also identified as outliers if the cluster median QC metrics deviated by more than three times the IQRs from the median of all clusters. In total, 15,928 nuclei passed QC criteria and were split into three broad classes of cells (10,708 excitatory neurons, 4,297 inhibitory neurons, and 923 non-neuronal cells) based on NeuN staining and cell class marker gene expression
Clusters were identified as donor-specific if they included fewer nuclei sampled from donors than expected by chance. For each cluster, the expected proportion of nuclei from each donor was calculated based on the laminar composition of the cluster and laminar sampling of the donor. For example, if 30% of layer 3 nuclei were sampled from a donor, then a layer 3-enriched cluster should contain approximately 30% of nuclei from this donor. In contrast, if only layer 5 were sampled from a donor, then the expected sampling from this donor for a layer 1-enriched cluster was zero. If the difference between the observed and expected sampling was greater than 50% of the number of nuclei in the cluster, then the cluster was flagged as donor-specific and excluded. In total, 325 nuclei were assigned to donor-specific or outlier clusters that contained marginal quality nuclei and were excluded from further analysis. Three donor-specific clusters came from neurosurgical donors (n=95 nuclei) and were similar to other layer 5 types reported in our analysis, but had higher expression of activity-dependent genes.
To confirm exclusion, clusters automatically flagged as outliers or donor-specific were manually inspected for expression of broad cell class marker genes, mitochondrial genes related to quality, and known activity-dependent genes.