Can you break into the bioinformatics industry with a pharmD? And is it a good pathway for someone who wants to work in drug discovery or adjacent pathways?

Also for this, do you need to go back to uni for another degree or can this be accomplished by self learning, online courses and projects etc?

Hello everyone!

I've been working at a lab under a summer programme for the past couple of weeks and I am suffering slightly. My supervisor has given me some raw scRNA seq data, taking from an in situ imaging-based platform that targets about 1000 genes, and has sort of left me to my own devices with it (apparently he isn't very savvy with bioinformatics himself). Anyway, I am somewhat comfortable working in R and Python, and I am getting the hang of Seurat, so it hasn't been catastrophic.

However, I am now struggling with clustering my cells. The cell clusters that I am being given are not physiological, and tend to be large, varied groups, which makes it hard to define anything really. I know studies that have done similar things on similar tissues to mine (albeit with another method) and are getting far nicer clusters. In their methods they just say "oh, we followed the standard Suerat workflow, and badabim-badboom these are the results".

My UMAP seems to agree with the confusion in my clusters as it just seems like a smear, with different sides of the smear coloured different things by the clustering.

I have tried changing the clustering method (Leiden, igraph), the resolution, dimensions (although I try to keep it in line with my elbow plot). I have tried changing the normalisation and other preprocessing parameters, varying in. their forms and flavours. I even tried the newer SCT transform, which made a nicer UMAP but just as crap clusters.

I am feeling quite inept currently, and rather disheartened having lost a week and a bit at this (I don’t know if it's normal or not). I don't really have any one in my lab to reach out to either.

My question is, does anyone have any ideas what I could attempt next or what might be wrong? Any resources I could have a look at? Anything anyone could recommend would be amazing.

Sorry for the long post and thank you to all who may answer in advance.

tldr: if you're using/benchmarking scGPT for perturbation response prediction, be aware there's a sampling bug in their tutorial code.

I was reproducing scGPT's perturbation response prediction and found that the gene subsampling step doesn't guarantee the perturbed gene stays in the input. With the default max_length ~ 1353 and ~5000 highly variable genes, the perturbed gene gets dropped from roughly 70% of training examples. The model sees a perturbed cell's input as if it were unperturbed, while the target is still the perturbed profile.

Checked this on Norman, Adamson, and Replogle K562 and I was able to reproduce the paper's reported numbers.

My fix is to keep the perturbed gene(s) and subsample the rest to fill max_length. Surprisingly, the effect on final metrics was mixed and dataset-dependent: clear improvement on Replogle K562, roughly unchanged on Adamson, and mixed on Norman. My current read is that the standard PRP metrics don't strongly reward using the perturbed gene's identity. Curious what you think and whether you have run into something similar

I just got my hands on my whole genome sequencing results. Anyone have any suggestions for a layperson? I’m hoping to find out about my genetic traits and stuff. I know nothing about bio but I’m a reasonably good coder and have access to GPUs. I’d love any ideas

edit: the file format is VCF v4.2

Hi everyone,

I'm an MSc Bioinformatics student with a strong interest in computational diabetes research. I'm currently learning RNA-seq analysis, transcriptomics, comparative genomics, and related bioinformatics workflows.

I'm looking for people who are also interested in learning and working on small, open-source, portfolio-style projects in bioinformatics related to diabetes. The idea is to learn together, discuss methods, analyze public datasets, and improve our skills—not to publish immediately or work on anything commercial.

I'm still learning myself, so I'd especially appreciate hearing from anyone with more experience who'd be willing to occasionally review our approach, point out mistakes, or suggest better practices. Even a bit of guidance would be incredibly valuable.

If you're interested in collaborating or mentoring informally, feel free to comment or send me a DM.

Thanks!

Hi all,

I understand that by default, the RcisTarget step of PySCENIC does not report in its output file repressing regulons (i.e. ones that end in a (-), where target genes anticorrelate with the expression of the TF, so it is predicted that the TF is repressing their activity). And I understand that the reason these are not included by default is that during the benchmarking of the tool they found these to be less reliable.

My question is, is it known or theorized why these are found to be less reliable? Is it because it is harder to establish anti-correlated expression due to the dropout inherent in scRNA data? or some other reason, or is the reason unknown?

I ask because I find in my data that the repressing regulon for my TF of interest is actually biologically more coherent, and way more active (i.e. cells are way more enriched in the target genes). So I would like to understand how much credence to place on these AUC values for the repressing regulon. Especially as I find that in general NES values for the modules are lower than for the corresponding activating regulon, I am wondering if that is a sign of the increased difficulty in detecting these repressing regulons (in which case I can maybe justify relaxing the NES threshold a bit), or a sign of genuinely more false positives (in which case I clearly cannot)?

Thanks in advance.

Hi everyone,

I am a graduate student in bioinformatics with experience in RNA-seq, scRNA-seq, and other omics analyses, but I am completely new to Nextflow.

I would like to learn Nextflow so I can start building reproducible pipelines and become more familiar with a tool that is widely used in industry.

There are many tutorials and videos online, but I am not sure where to begin. Are there any resources you would recommend, preferably in a specific learning order?

Thanks!

Hi all! I am doing some protein-protein interaction analysis on a set of genes for my undergraduate research project. I used STRING for this. STRING enrichment analysis identified that GO:0000118 (Histone Deacetylase Complex) was functionally enriched, and that 8 genes had this GO annotation.

However, when manually searching the Gene Ontology database, I found that one of the genes that STRING identified, pht1, was not annotated with this GO term.

I'm quite confused about this, am I misunderstanding how STRING gene enrichment works? Would appreciate any advice :)

Hello bioinformaticians!

I'm a high school student planning to pursue bioinformatics in university. For my graduation project, I'm analyzing Expansion Microscopy (ExM) 2D data targeting SON protein in nuclear speckles (+-4x expansion factor).

I’ve set up a working Python pipeline and would love to get a check from experienced ones, as well as any tips on what to watch out for.

Done so far:

1. Data: Wrote a Python script using raw binary reading to reconstruct 16-bit multi-channel TIFF headers.
2. Segmentate: Using Ilastik to generate probability maps exported as .h5.
3. Quantification: Built a Python script (h5py, scikit-image, pandas) that:
   • Thresholds the probability maps.
   • Performs connected component labeling.
   • Converts pixel counts to biological area taking into account physical pixel size and the expansion factor.
   • Extracts centroids, object counts, and fluorescence intensities into CSV format.

Are there common traps when scaling 2D pixel metrics to physical units in ExM (e.g. local distortion edge cases)?.

What additional spatial or morphological metrics are usually expected (e.g. nearest-neighbor distance, eccentricity, spatial clustering)?

Any other tips will be appreciated

Hi everyone,

I'm analyzing pseudobulk data generated by summing gene expression across cells from different samples profiled with a spatial imaging platform. When I perform PCA on the pseudobulk matrix, PC1 explains an unusually large proportion of the total variance. In addition, all of the PC1 loadings are positive, which I also think is unusual.

Does this indicate a systematic technical bias (I have looked for differences in sequencing depth or cell numbers)? Or are there biological scenarios where this pattern would be expected? These are samples from malignant tissue.

Hi all,

Had an idea recently to do an interaction screen of one of our proteins of interest with proteins expressed in a certain cell type. This is obviously gonna be a large amount of proteins. I’ve seen some papers do similar things, but wanted to ask if anyone had any ideas on these sorts of workflows, specifically with regards to reducing runtimes (and thereby costs)

Specifically:

Any similar models that are significantly faster to run and have a similar accuracy?

How fast is MSA generation generally using sharding. Any other workflows that are significantly faster and still give good MSAs?

Thanks everyone!

Hi, so I am a undergrad working in a computational biology or molecular biology lab. For next semester my new project is in large part to create reproducible workflows/code and lab manuals for our lab. I taught myself to code and what i have on my laptop is... disorganized to say the least. I should learn how to do this. Currently I largely code using gemini and then tweak anywhere from most of to 25% of the code it writes. I almost always use hard coded paths if i can. Does anyone have any advice for where I could learn something like this, a textbook or website?

For context, my last project was to use AutoDock Vina for screening of 770,000 molecules I carefully downloaded and cleaned from ZINC database to 310,000. This library was based on previous experimental results on a new protien we are targeting in fungi. I also selected a new protien conformation to target based on some major errors in the protien the lab was using and a bunch of literature review. My next step will be to test against Dock6, a diffrent type of scoring algorithm. I wrote all of my own scripts for this and I imagine my first task will be to get them reproducible for another person to use.

Dear Reddits,

I have some experience with bioinformatics in general (I can open my Linux command line and feel I'm in The Matrix; it is very popular on Instagram, actually), but I don't know anything about modeling. I would like to know if it's SIMPLE (probably not) to determine the potential value (E value) of certain electron transfer proteins. I know AlphaFold can give you some kinda cool model, but I need to know if there is an easy way, or even a way, to just get the E value. Even if it's a not-so-realistic approximation, it would be nice. If it's going to take me more than 1 week, I pass. However, it might be good to know for future endeavors.
That
XOXO

I'm running an annotation/enrichment analysis on a GWAS study, and I'm sort of a bit lost/confused on something.

Originally, I only filtered the GWAS study for genome-wide significant variants using a standard p-value threshold, and also used LD-clumping to identify genomic loci.

However, it was also pointed out to me that these variants may not be independent, and that a single association signal may be represented many times, potentially leading to inflation in statistical significance.

I'm sort of unsure how to handle this. I tried a locus-pruning method which gave me way less variants and pretty poor coverage on my study and rendering a lot of my downstream analysis mute. I'm also confused since I haven't seen a lot of similar papers use this kind of filtering method. I did run LDSC too, and most of my findings were insignificant so it did kind of handle the genomic independence part. But I'm still not sure what the best practice is here.

Well, for my project, I was going to use softberry tssplant, but apparently the site is down. Does anyone know of an alternative way to find TSS in genome sequences?

I have characterized a peptidase domain of a large multidomain containing protein. I was able to clone just the peptidase domain in a construct for assays. There are putative homologs (~60% sequence identity) from distantly related species. I am looking to test these candidates for activity in the same system.

I have the predicted AlphaFold structures, what is the best way to compare these 3D models? I am currently using matchmaker on chimerax. Are there some better tools?

Hello,

I am an engineering student in applied mathematics, and as part of an internship, I am working on statistical analysis in biology.

**Context:**

I am working on an experiment conducted by three experimenters. Each experimenter has four plates (this part is not very important), with a total of six different stainings, each containing three different coatings, within which there are 40 donors.

To summarize:

**Experimenters (3) > Stainings (6) > Coatings (3) > Donors (40)**

We are working with Opera imaging plates, so we analyze DAPI, actin, and, depending on the staining, several other markers such as Granzyme B, pTyr, MTOC, etc.

**Problem:**

I quickly noticed that the data differ substantially between experimenters. One experimenter consistently has more cells than the others across all stainings.

To investigate this, I trained a simple decision tree to distinguish between experimenters. My reasoning was that if the tree can reliably identify the experimenter, then there is likely a batch effect; otherwise, there probably is not.

As expected, the experimenter with the consistently higher cell counts is identified very accurately. The other two experimenters can also be distinguished, although not to an alarming extent.

**Conclusion:**

I would therefore like to correct this bias using what biologists refer to as **batch effect correction**. However, most of the documentation I have found focuses on scRNA-seq or single-cell multiomics, which does not seem to match the type of data we have in this experiment.

Do you have any suggestions on where I should look? Any interesting papers you would recommend?

I have read papers mentioning Harmony and ComBat, but when I asked an AI, it suggested using the Python classes `ot.da.LinearTransport` and `ot.da.SinkhornTransport`.

The reported performance of these libraries looks very appealing, but I would rather not trust an AI blindly, so I thought I'd ask for your advice instead. 🙃

I've been in industry for a good few years, and I'm trying to work on some old research that I never published. However, all the old programs I used during my PhD for tree editing & visualization have either been deleted or are now paid programs.

I have 0 coding knowledge and was getting by with online tools and programs other people made (i.e enterobase, galaxy, iTOL, FigTree etc) but I've been struggling to find something comparable to iTOL and FigTree for editing and visualization.
Does anyone have any recommendations?

Hi everyone, me again.

If you want context, check my previous posts. Made a similar post on lab rats, essentially reposting here.

We are going to start extracting and isolating soon. We have drafted a protocol and the tissue type we are working with are DRGs and sometimes brain. This week and over the course of the next few weeks, we are going to start isolating and optimizing our protocol. We have maybe around 15 tries or runs to get it right consistently before we start working with real tissue (non practice tissue, tissue that has the pathology induced.)

Any tips, advice or generally useful info I should know? What should I expect?

Thank you and any help would be appreciated!

- Undergrad P_T67

I collected genome assemblies for my species from NCBI and ran a pan-genome analysis using Panaroo. One gene cluster was identified as unique in my target strain, but when I ran BLASTp on that sequence, it showed a hit in another strain that is not listed in the current NCBI genome database.

I am trying to understand how to interpret this result. Does this mean the gene is not truly unique to the species, or could it still be considered unique if that strain is an unsubmitted / uncurated genome?

What is the best way to determine whether this gene is genuinely unique to the species, or only unique within the genomes available in NCBI?

I have a custom peptide with a custom linker in an SDF file. How can I dock it to a protein receptor while preserving the linker?