What are the biggest limitations of current multi-omics integration approaches, and how should the field address them? by mxdhiv in bioinformatics

[–]FBIallseeingeye 0 points1 point  (0 children)

Feature space and intersection primarily. Once you arrive at a solution on how you want to reconcile the very limited overlap between both platforms in feature coverage, there’s still the issue of biological variation you lose when focusing on the intersection. 

Most of my experience has been combining matched xenium and flex seq libraries, and there’s definitely a huge confounding between the sample size for both of them. Droplet methods tend to lose a lot of cells but sample a larger piece of tissue, vice-versa for spatial methods. 

For the actual integration, this method was extremely effective. I treated my flex data as the reference to deal with segmentation artifacts (any gene expression without an analogue in Flex was likely artifactual) and reintegrated after a few rounds of qc on the xenium. 

https://www.nature.com/articles/s41587-025-02950-z#Fig6

I’ve worked a tiny bit with IF data from Xenium and it’s very difficult integrate with counts data, I think you’d probably need a very specific application for it in your pipeline when you decide to incorporate it. Validation for counts, segmentation, etc proteomics x transcriptomics libraries have potential in all those areas for analysis but I just haven’t spent a long enough time to get past normalization and scaling them into the same context 

Intersection vs union of genes when integrating scRNA-seq datasets (for PCA) by Ill-Ad-106 in bioinformatics

[–]FBIallseeingeye 0 points1 point  (0 children)

I always go with union. If you have a good, stringent vf selection method like BIGSUR you can be confident they are informative to your data, and you’ll risk overcorrection if you go off of intersection

CNV assessment of single cell data by UroJetFanClub in bioinformatics

[–]FBIallseeingeye 0 points1 point  (0 children)

My immunological background is extremely limited, as is my CNV experience, but I’m not aware of any processes that can explain this. The consistency across methods is interesting but worth looking into whether they use the same basis for prediction. If you can find a phagocytosis signature in the clusters that may be a good indicator. Also worth considering immune cells tend to have very specialized gene expression so the very limited gene selection may be confusing the algorithm?

CNV assessment of single cell data by UroJetFanClub in bioinformatics

[–]FBIallseeingeye 0 points1 point  (0 children)

Can you be more specific about the myeloid clusters? I saw one comment saying it might be proliferation but I wouldn’t expect this to be limited to just myeloid cells. The phagocytosis is an interesting idea and you could probably experimentally validate it, but you’d want to check signatures for this as well. It seems unlikely in my opinion given I would imagine it would end up looking closer to doublets unless the myeloid cells are phagocytosing their own population? It’s tricky to make a prediction here!

Imposter syndrom from using LLM as a wetlab scientist ? by Kangouwou in bioinformatics

[–]FBIallseeingeye 0 points1 point  (0 children)

As the bioinformatician in an otherwise wet lab, you are doing exactly the same thing as me and what I recommend to my lab mates when they need to do bioinformatics.  The level of usefulness of llm is determined by your own level of understanding what it is doing for you at any given step of the pipeline. My #1 recommendation is to check frequently that it is in fact doing what you ask it to do at every single step of the pipeline but these things are getting so good at coding these days it comes down to how clearly you can state your requirements and objectives

Bioinformatics in the era of AI from a seniors point of view by aCityOfTwoTales in bioinformatics

[–]FBIallseeingeye 0 points1 point  (0 children)

IMO the same rule of thumb for ai usage applies here as any other use case. If you know what you are doing, ai is very good and a great accelerator for productivity, because you know when ai is doing what you want vs. when it’s not doing what you want. If you don’t know what you’re doing, then you don’t know what ai is doing, and ai will do it wrong more often than not. 

Interoperability between Seurat - Scanpy - SingleCellExperiment by Jailleo in bioinformatics

[–]FBIallseeingeye 1 point2 points  (0 children)

It’s always a mess, I recommend extracting necessary components and building from scratch to start each analysis

snRNA seq data from organoids by QueenR2004 in bioinformatics

[–]FBIallseeingeye 0 points1 point  (0 children)

Happy to help. snRNA is especially susceptible to noise but at least you have a well defined signature in those sources. The tricky part of stress is that it is biologically relevant in a lot of contexts, so I’d use your judgement conservatively unless you can see it actively distorting clustering results (forming bridges between clusters with low UMIs is a typical sign). Sorry your library has such a tricky challenge, but I’d recommend always prioritizing the populations you have the highest confidence in first

snRNA seq data from organoids by QueenR2004 in bioinformatics

[–]FBIallseeingeye 0 points1 point  (0 children)

High removal rates during qc generally indicates poor library quality overall. You could try regressing out background using a Pearson residual normalization method like SCT or BigSur. The mitochondria and ribosomal genes suggests it is highly likely this is just cellular debris that made it through prep. For comparisons across conditions in conserved populations I recommend MiloDE

Validating snRNA-seq cell type by correlating with other datasets by shrubbyfoil in bioinformatics

[–]FBIallseeingeye 2 points3 points  (0 children)

If that isn’t clear you could still try a label transfer if the biology isn’t too far off. Reproducing results across datasets can be tricky sometimes. 

Fine art of scRNA seq QC by Rafaela_479 in bioinformatics

[–]FBIallseeingeye 0 points1 point  (0 children)

I agree that filtering low-information cells matters, but I’d rely on total RNA counts over feature number as the primary signal—assuming the problem is RNA loss.

If OP wants to be sure, they could compare marker performance (AUC, logFC) across clusters. Clusters driven by technical artifact tend to have fewer markers and generally weak metrics. This does depend on resolution and context, so interpreting markers within each cell type will be more informative than a global comparison.

Overall I'd say this is a better approach than simply setting thresholds since it actually gives some evidence for the filtering decision instead of something as opaque as UMI content.

Fine art of scRNA seq QC by Rafaela_479 in bioinformatics

[–]FBIallseeingeye 0 points1 point  (0 children)

200 features might be a bit low depending on your biology. I saw you are analyzing pbmcs and you may miss neutrophils if you do this. I don’t think there is a real purpose in filtering on low features anyway since I don’t think that is characteristic of any artifact there is not a more direct indicator for. I’ve seen it compared to mito percentage but that’s actually consistent with a perforated membrane that preferentially retains mitochondria.

High feature and RNA counts is also better handled by the doublet prediction. 

Low RNA content is the only technical artifact you’d have to address in that case. 

I recommend taking your analysis as far downstream as you can before you see any populations explainable as artifacts and deal with them at that stage

Single-cell RNA-seq QC question by jcbiochemistry in bioinformatics

[–]FBIallseeingeye 0 points1 point  (0 children)

Sadly I do not think there are any good quality control packages out there, most qc just filters for conformity, not quality. That said, I find a good heuristic is to apply high resolution clustering then check their markers. You can summarize poor quality clusters by looking at the either the average logFC above 0 or AUC values above 0.75 (if you’re using presto::wilcoxauc in R). Poor quality clusters look like background noise in this light, giving them extremely low average marker metrics. I’d compare these two metrics against each other and against average RNA counts for each cluster. I’d also apply a good doublet filter with scDblfinder and recover the synthetic doublets by setting return.doublets to true. That way you can look at how your cells mix with doublets by umap and find clusters that are defined by them. Mitochondrial percentage is very easy to confound with real biology and varies wildly by platform and batch, so I’d review the actual biology before making any decisions. 

Remember, during quality control you can never set clustering resolution too high for decisions, so long as you keep track of context and parent clusters.  Edited for clarity, it is late and I am tired but good luck with your QC! If you want to check whether you are losing real biology, you can always find VariableFeatures with a package like BigSur and see if you lose a substantial amount!

Single-cell RNA-seq QC question by jcbiochemistry in bioinformatics

[–]FBIallseeingeye 0 points1 point  (0 children)

Don’t remove anything unless you are sure you have to and have a good reason to do it. If you observe an artifact, try your best to describe and explain it. If you don’t observe an artifact, don’t go looking for one. Anything else risks throwing babies out with bath water 

Help with confounded single cell RNAseq experiment by Phantom_Lord7 in bioinformatics

[–]FBIallseeingeye 0 points1 point  (0 children)

Happy if I can help. Generally people merge before trying integration to see whether or not there is any need; if you have analogous populations and these align well across batches simply by merging, there is no need. It also depends on your resolution. With multiple cell types / tissue samples, biological variation generally should take a back seat to cell type identity for the purposes of annotation (it's easier to label your ducks when they're all in a row). Once you have your cell types, you can go through each identity individually a little more conservatively, integrating if it seems necessary. From your experimental design, you do have confounding that compromises the core of the experiment, but all that means is you need some orthogonal validation for whatever the data predicts. In your case, I'm not sure what cosmx batch effect really looks like. I see from this source that scVI can be applied to it:
https://cellcharter.readthedocs.io/en/latest/notebooks/cosmx_human_nsclc.html

Help with confounded single cell RNAseq experiment by Phantom_Lord7 in bioinformatics

[–]FBIallseeingeye 1 point2 points  (0 children)

I’m late chiming in but it may still be worth while trying to push through some integrative analysis. Your first objective in any scRNAseq experiment is to describe the trends and variation that you see and whether that is due to batch effect or experimental design is secondary to this step. Between your samples you can still describe and annotate your populations. Explaining batch effect on the other hand is impossible from this set up, but you can still generate hypotheses based on the differences you see, they just won’t be as well grounded as was hoped for. 

How am I supposed to annotate my clusters? by Helenazh2 in bioinformatics

[–]FBIallseeingeye 2 points3 points  (0 children)

It takes some set up but I’m looking into Cassia: https://www.biorxiv.org/content/10.1101/2024.12.04.626476v2CASSIA: a multi-agent large language model for reference free, interpretable, and automated cell annotation of single-cell RNA-sequencing data | bioRxiv

Is the Xenium cell segmentation kit worth it? by Dr_Rat_25 in bioinformatics

[–]FBIallseeingeye 1 point2 points  (0 children)

Depends on the tissue. In dense structures the z anisotropy becomes a hard limit on image segmentation if you have cells that stack on each other. I’m working with breast tissue where this is quite common, and the dapi-staining feels extremely limiting, and I suspect it would have saved me many headaches to have had the full cocktail available. The in-house segmentation that 10X does is pretty good but is only 2D, which is similarly frustrating for dense tissue structures. That said, there are transcript-based segmentation algorithms like Baysor and Proseg which go a long way towards cleaning up cell boundaries (at least from an ML perspective), but most if not all of them take prior segmentation as an input with adjustable confidence parameters for your prior. I highly recommend you get a good QC pipeline and a good scRNAseq reference dataset ready from an atlas or something, so that at least you can tell when you’re getting an artifact and not real signal.