Question on muscarinic receptor gene expression in WHB-10Xv3 dataset

Hi,

I retrieved expression data of muscarinic receptor gene CHRM1 - 5 in all neuron cells from the WHB-10Xv3 dataset. Attached figure 1 plotted gene expression in CPM and percent of detection, where CHRM4 expression was barely detected.

CHRM4 is known to be expressed at relatively abundant level comparing to other CHRM genes in selected brain regions. See attached figure 2 and 3 for gene expression in bulk brain tissues (GTEx data) and protein expression (Human Protein Atlas).

What may be the reason for the seeming lack of CHRM4 expression detection in the WHB-10Xv3 dataset. Does CHRM4 truly have lower expression than the other CHRM genes in neurons of the measured brain regions, or some factors contributed to its under estimation? There may be factors such as 10X Genomic platform, experiment QC, brain sample collection, CHRM4 sequence, … Can you please help assess it?

Thanks in advance,

Sean

Hi @xyaous ,

The short answer is that I don’t know, however you are not the first person to notice this happening in one of our data sets for a gene (e.g., this post: Known protein not detected by RNA explorer).

I think any of your suggestions are reasonable. In some cases the differences are technical, although it seems unlikely here since at least three human and one mouse data sets show similarly low expression of CHRM4 across most cells–so also not likely a cross-species difference. Another possibility (which I don’t think is likely given how transcripts are selected) is that the reads for CHRM4 might not be properly aligned with the location of the gene annotation in the genome. We can see in both MTG and M1 that reads pile up in both exonic and intronic locations.

What we can safely say is that CHRM4 does not have high expression in very many nuclei collected from human brains across a variety of studies. Finally, while it is possible that CHRM4 is expressed more highly in single cells rather than single nuclei, this does not appear to be the case in mouse (Single-nucleus and single-cell transcriptomes compared in matched cortical cell types).

I’m not sure what possibilities this leaves (maybe previous studies were mistaken?), but hopefully this post at least provides some directions to pursue and rules out some options.

Best,
Jeremy