text: user7113: I am analysing a human single cell RNA seq experiment, where we have 4 groups, four samples each. Data has been analysed using Seurat, with the canonical workflow. I have tried DE using various methods (Wilcox, MAST, LR), and for each of them the top genes (biggest fold change, lowest adjusted p) are genes over expressed by specific samples, e.g.:
  
  Looking at the average expression for these genes per each subject, and plotting on a Boxplot, it does not seem to be a genuine difference (e.g. for the gene above):
  
  Besides the subject, there is no other parameter (e.g. age, sex) that is biased across groups and that shows up as differentially represented in the PCA plot.
  Is there a DE method that takes into account this behaviour and corrects for it?
  How can these results be presented? 
  Besides increasing number of subjects (not really feasible), is there a solution to analyse these data?
  Thanks a lot!
  
  user8777: One potential problem here is that, as far as I can see, you're not taking the hierarchical structure of your scRNA-seq data into account. I.e. the fact that cells from the same patient / sample will be more similar than cells from different samples, even within the same group. This is also the problem addressed in the paper mentioned by @alan. Also, it usually makes more sense to do this kind of group vs. group comparison for each subpopulation (cluster/cell type) separately (assuming you have multiple subpopulations in your data and did some kind of clustering).
  One way of doing this is by using a pseudobulk-level approach, where you basically aggregate (by either averaging or summing) your data over each sample, for each cluster separately. So you end up with one counts matrix per subpopulation with one row for each gene and one column for each sample. The sample columns are then the sums (or averages) of all cells belonging to that cluster-sample combination.
  You can then use bulk-RNAseq methods such as edgeR or DESeq2 on the aggregated data to test differences between groups of samples for each cluster. The muscat package implements all these steps and has quite extensive documentation. So you could check that and test it out on your data.
  The paper describing the method of muscat can be found here: https://www.biorxiv.org/content/10.1101/713412v2
  
  user4507: Leave-one-out cross validation is an iterative method that leaves out one sample until each sample has been left out once. You can then identify the intersection of the differentially expressed genes from the sample space of repeated analyses to produce a set of genes that are consistently differentially expressed independent of individual sample biases.
  more info here:
  http://efavdb.com/leave-one-out-cross-validation/
  
original_text: 7113: <p>I am analysing a human single cell RNA seq experiment, where we have 4 groups, four samples each. Data has been analysed using Seurat, with the canonical workflow. I have tried DE using various methods (Wilcox, MAST, LR), and for each of them the top genes (biggest fold change, lowest adjusted p) are genes over expressed by specific samples, e.g.:
  <a href="https://i.stack.imgur.com/PCyTH.png" rel="nofollow noreferrer"><img src="https://i.stack.imgur.com/PCyTH.png" alt="enter image description here"></a></p>
  
  <p>Looking at the average expression for these genes per each subject, and plotting on a Boxplot, it does not seem to be a genuine difference (e.g. for the gene above):
  <a href="https://i.stack.imgur.com/l0041.png" rel="nofollow noreferrer"><img src="https://i.stack.imgur.com/l0041.png" alt="enter image description here"></a></p>
  
  <p>Besides the subject, there is no other parameter (e.g. age, sex) that is biased across groups and that shows up as differentially represented in the PCA plot.</p>
  
  <p>Is there a DE method that takes into account this behaviour and corrects for it?
  How can these results be presented? 
  Besides increasing number of subjects (not really feasible), is there a solution to analyse these data?</p>
  
  <p>Thanks a lot!</p>
   8777: <p>One potential problem here is that, as far as I can see, you're not taking the hierarchical structure of your scRNA-seq data into account. I.e. the fact that cells from the same patient / sample will be more similar than cells from different samples, even within the same group. This is also the problem addressed in the paper mentioned by @alan. Also, it usually makes more sense to do this kind of group vs. group comparison <strong>for each subpopulation (cluster/cell type) separately</strong> (assuming you have multiple subpopulations in your data and did some kind of clustering).</p>
  
  <p>One way of doing this is by using a <strong>pseudobulk-level approach</strong>, where you basically aggregate (by either averaging or summing) your data over each sample, for each cluster separately. So you end up with one counts matrix per subpopulation with one row for each gene and one column for each <strong>sample</strong>. The sample columns are then the sums (or averages) of all cells belonging to that cluster-sample combination.</p>
  
  <p>You can then use bulk-RNAseq methods such as <strong>edgeR</strong> or <strong>DESeq2</strong> on the aggregated data to test differences between groups of samples for each cluster. The <a href="https://bioconductor.org/packages/release/bioc/html/muscat.html" rel="nofollow noreferrer">muscat</a> package implements all these steps and has quite extensive documentation. So you could check that and test it out on your data.</p>
  
  <p>The paper describing the method of <em>muscat</em> can be found here: <a href="https://www.biorxiv.org/content/10.1101/713412v2" rel="nofollow noreferrer">https://www.biorxiv.org/content/10.1101/713412v2</a></p>
   4507: <p>Leave-one-out cross validation is an iterative method that leaves out one sample until each sample has been left out once. You can then identify the intersection of the differentially expressed genes from the sample space of repeated analyses to produce a set of genes that are consistently differentially expressed independent of individual sample biases.
  more info here:
  <a href="http://efavdb.com/leave-one-out-cross-validation/" rel="nofollow noreferrer">http://efavdb.com/leave-one-out-cross-validation/</a></p>
  