Seurat’s AddModuleScore function was applied to quantify differentiation status of each thyrocyte in our scRNA-seq dataset. We calculated the mean abundance levels of the abovementioned 13 thyroid function-related mRNA genes against the aggregated abundance of random control gene sets as the TDS score of each thyrocyte. Similarly, the BRAF and RAS score of each thyrocyte were quantified with AddModuleScore() function based on the upregulated signature genes in the BRAFV600E-mutated and RAS-mutated samples from the bulk TCGA profiles (Supplementary Data 6).
Do you have any questions about this protocol?
Post your question to gather feedback from the community. We will also invite the authors of this article to respond.