Steven Robbins's Avatar

Steven Robbins

@stevenjrobbins.bsky.social

Do my science @ace_uq studying coral reef microbiomes. Data wrangler, meta-omics and long-read wonk, clean energy enthusiast, Saganist zealot, collector of weird zoology facts, other nonsense.

3,476 Followers  |  514 Following  |  818 Posts  |  Joined: 20.11.2023  |  1.9794

Latest posts by stevenjrobbins.bsky.social on Bluesky

Preview
Theย Microflora Danica atlas of Danish environmental microbiomes - Nature Microflora Danicaโ€”an atlas of Danish environmental microbiomesโ€”reveals that although human-disturbed habitats have high alpha diversity, species reoccur, revealing hidden homogeneity.

Microflora Danica: What can you learn from collecting and sequencing 10,000+ samples from a single country? Check out our new paper in @nature.com to find out. Incredible work led by Caitlin Singleton, Thomas B. N. Jensen, and Mads Albertsen from @aau.dk. ๐Ÿฆ ๐Ÿงซ๐Ÿงฌ
www.nature.com/articles/s41...

03.12.2025 20:50 โ€” ๐Ÿ‘ 64    ๐Ÿ” 29    ๐Ÿ’ฌ 0    ๐Ÿ“Œ 2
Preview
Meta-virus resource (MetaVR): expanding the frontiers of viral diversity with 24 million uncultivated virus genomes Abstract. Viruses are ubiquitous in all environments and impact host metabolism, evolution, and ecology, although our knowledge of their biodiversity is st

๐Ÿฆ ๐Ÿงช๐Ÿงฌ๐Ÿšจ New paper and database alert: the new IMG/VR release is now MetaVR ! We have a new website - meta-virome.org - with quick search capabilities for the >24M viruses, >12M vOTUs, and >42M protein clusters (including >790k with predicted structures !). academic.oup.com/nar/advance-...

03.12.2025 02:34 โ€” ๐Ÿ‘ 62    ๐Ÿ” 41    ๐Ÿ’ฌ 1    ๐Ÿ“Œ 1
Preview
Hologenomic insights into the molecular adaptation of deep-sea coral Bathypathes pseudoalternata Wei et al. reveal how deep-sea corals partner with symbionts in the extreme environment using the hologenome of Bathypathes pseudoalternata. A division of labor shows that the host provides a protecte...

Symbiont-coral partnerships

Hologenome of B. pseudoalternata reveals a division of labor between deep-sea coral & symbionts. Host provides habitat & manages immunity, while symbionts recycle waste, supply nutrients, resist viruses & produce antioxidants
www.cell.com/cell-host-mi...

20.11.2025 20:03 โ€” ๐Ÿ‘ 5    ๐Ÿ” 1    ๐Ÿ’ฌ 0    ๐Ÿ“Œ 1

Enormous milestone. All thanks to vaccinations!!!

27.11.2025 03:02 โ€” ๐Ÿ‘ 15    ๐Ÿ” 4    ๐Ÿ’ฌ 1    ๐Ÿ“Œ 0
Preview
Long-read metagenomics reveals phage dynamics in the human gut microbiome - Nature Complex prophage integration dynamics, including low-level induction, cross-family host range and transposase-mediated mobilization, challenge existing paradigms and deepen our understanding of phageโ€“...

Long read Metagenomics, #phage and #prophage in the gut by Ami Bhatt's group. Beautiful data showing changes in phages over two years

#phagesky

www.nature.com/articles/s41...

26.11.2025 21:55 โ€” ๐Ÿ‘ 64    ๐Ÿ” 34    ๐Ÿ’ฌ 0    ๐Ÿ“Œ 1

Thanks Greg. Quarto is cool, but I don't really need to render the code or make it overly fancy. I really just want a notebook to markdown-style application to document what i'm doing, with the code that i'm using to execute things on the server.

26.11.2025 23:00 โ€” ๐Ÿ‘ 0    ๐Ÿ” 0    ๐Ÿ’ฌ 1    ๐Ÿ“Œ 0

Hey Team! What do you think are the best free digital notebooks these days.

I've used:
Evernote--great, but now it's paid

TiddlyWiki wasn't bad, but got rid of scrollable code blocks

Obsidian--basic, and I don't like not being able to scroll my code blocks like OG Tiddlywiki.

Suggestions?

26.11.2025 21:49 โ€” ๐Ÿ‘ 5    ๐Ÿ” 2    ๐Ÿ’ฌ 2    ๐Ÿ“Œ 0

๐Ÿ“Œ

26.11.2025 15:32 โ€” ๐Ÿ‘ 0    ๐Ÿ” 0    ๐Ÿ’ฌ 0    ๐Ÿ“Œ 0

Thank you! Can I ask, do you find it better than ChatGPT5 at the moment?

22.11.2025 15:22 โ€” ๐Ÿ‘ 0    ๐Ÿ” 0    ๐Ÿ’ฌ 0    ๐Ÿ“Œ 0

Super helpful. Thanks Sebastian!

21.11.2025 16:54 โ€” ๐Ÿ‘ 0    ๐Ÿ” 0    ๐Ÿ’ฌ 0    ๐Ÿ“Œ 0

Thanks man! Granted, what I saw with my eyes was a much fainter version of these, and these were points of the night that were strongest. Most of the night you couldn't see much. So these are definitely the enhanced highlight reel.

21.11.2025 16:54 โ€” ๐Ÿ‘ 1    ๐Ÿ” 0    ๐Ÿ’ฌ 0    ๐Ÿ“Œ 0
Preview
proGenomes4: providing 2 million accurately and consistently annotated high-quality prokaryotic genomes Abstract. The pervasive availability of publicly available microbial genomes has opened many new avenues for microbiology research, yet it also demands rob

proGenomes4: providing 2 million accurately and consistently annotated high-quality prokaryotic genomes academic.oup.com/nar/advance-... #jcampubs

20.11.2025 14:02 โ€” ๐Ÿ‘ 18    ๐Ÿ” 8    ๐Ÿ’ฌ 0    ๐Ÿ“Œ 0
Preview
A Grammar of Graphics for Comparative Genomics An extension of ggplot2 for creating complex genomic maps. It builds on the power of ggplot2 and tidyverse adding new ggplot2-style geoms & positions and dplyr-style verbs to manipulate theโ€ฆ

gggenomes: A Grammar of Graphics for Comparative Genomics thackl.github.io/gggenomes/ #Rstats

18.11.2025 19:02 โ€” ๐Ÿ‘ 68    ๐Ÿ” 31    ๐Ÿ’ฌ 0    ๐Ÿ“Œ 0

Them too, for sure!

18.11.2025 21:06 โ€” ๐Ÿ‘ 1    ๐Ÿ” 0    ๐Ÿ’ฌ 0    ๐Ÿ“Œ 0

For mmseqs, I would think assembled contigs make more sense, but I don't know who i'd trust more than Titus to know where and when to trust an assembly.

18.11.2025 20:02 โ€” ๐Ÿ‘ 2    ๐Ÿ” 0    ๐Ÿ’ฌ 1    ๐Ÿ“Œ 0

That's the nice thing about long reads. The reads themselves are the primary data, and long enough to be their own evidence. It's not underrated. The issue now is it's kind of easy at the moment to find skeletons hiding in the assembly/binning closets, we just couldn't see them before.

18.11.2025 19:57 โ€” ๐Ÿ‘ 2    ๐Ÿ” 0    ๐Ÿ’ฌ 1    ๐Ÿ“Œ 0

Just from experience, it's the way i'd go with, say, long read data or contigs with at least a few genes on them. If there's half a gene on each "read," LCA classification won't go that far...maybe at least to domain, though, which is the goal here. Curious to see if it works.

18.11.2025 18:36 โ€” ๐Ÿ‘ 1    ๐Ÿ” 0    ๐Ÿ’ฌ 0    ๐Ÿ“Œ 0

If they're like 150bp reads, you're not going to get good classification though. It'll get better the more genes are on those contigs. There's really just no good way to classify short reads that's reliable, as far as i'm aware.

18.11.2025 18:14 โ€” ๐Ÿ‘ 2    ๐Ÿ” 0    ๐Ÿ’ฌ 0    ๐Ÿ“Œ 0

Ah ok, so by this I think you mean you have a subset of reads in a MAG that you want to classify as euk or prok. If you have something like a full gene or multiple genes, LCA might be what you want. Works better on assembled contigs since it's looking at what multiple genes on the sequence blast to.

18.11.2025 18:13 โ€” ๐Ÿ‘ 0    ๐Ÿ” 0    ๐Ÿ’ฌ 1    ๐Ÿ“Œ 0

Yup...I would just say there's no good answer to this. Best I can think of is LCA-assignment in mmseqs. At least you might get a "most" likely it's in this domain of life. I tend to stay away from pure classifiers for things like this except in special cases.

18.11.2025 18:08 โ€” ๐Ÿ‘ 4    ๐Ÿ” 0    ๐Ÿ’ฌ 1    ๐Ÿ“Œ 1

Those suggesting below to look for ribosomal genes, etc, work for the first case, not for the second. If they're random contigs you want to classify, the only good way to really do it is something like Lowest common ancestor-assignment as implemented by mmseqs.

18.11.2025 18:04 โ€” ๐Ÿ‘ 0    ๐Ÿ” 0    ๐Ÿ’ฌ 0    ๐Ÿ“Œ 0

Need a bit more info on what you mean. By "binned into a bacterial genome," do you mean that you want to make sure the MAG is prokaryote and not euk, or a subset of reads that got binned into a MAG are not euk. The first is easier, the second probably doesn't have a good answer.

18.11.2025 18:01 โ€” ๐Ÿ‘ 0    ๐Ÿ” 0    ๐Ÿ’ฌ 2    ๐Ÿ“Œ 0

I think heโ€™s saying heโ€™s got some random genes in a bin, right? For graftM or singlem to work theyโ€™d have to be a core marker gene or have a specifically curated graftM package to classify with. Neither seem likely?

18.11.2025 17:57 โ€” ๐Ÿ‘ 0    ๐Ÿ” 0    ๐Ÿ’ฌ 1    ๐Ÿ“Œ 0
Preview
ECR connection: Meet Meren when you need to A means for ECRs to get advice from a senior scientist outside of their support network

Starting this week, I set aside one hour each week to meet ECRs outside my group who want to discuss career development, mentorship, or any non-technical professional questions.

Here is a blog that explains my motivation for this and how to schedule a meeting:

merenlab.org/2025/11/16/E...

16.11.2025 16:56 โ€” ๐Ÿ‘ 51    ๐Ÿ” 29    ๐Ÿ’ฌ 1    ๐Ÿ“Œ 0

โ€œMy magic beans can work miracles, theyโ€™re different from all of the fake magic beans out thereโ€ has been the sales pitch of charlatans since time immemorial.

17.11.2025 01:46 โ€” ๐Ÿ‘ 146    ๐Ÿ” 21    ๐Ÿ’ฌ 4    ๐Ÿ“Œ 2

My colleague saw no reads mapping to some regions and other reads all truncating all at the same bp position like their ends didn't map anywhere.

But really, first i'd just see if your 5S is really there by blasting some closely related species' 5S to it. Then go down a rabbit hole if it's there.

16.11.2025 04:50 โ€” ๐Ÿ‘ 0    ๐Ÿ” 0    ๐Ÿ’ฌ 0    ๐Ÿ“Œ 0

So, my (probs incomplete) understanding of minimap is it's producing HSPs, best alignments of parts of each read, not trying to align the whole read. What we wanted to see with our Hifiasm misassemblies was whether reads spanned a whole region, and they didn't--i.e. they didn't support the assembly.

16.11.2025 04:46 โ€” ๐Ÿ‘ 0    ๐Ÿ” 0    ๐Ÿ’ฌ 2    ๐Ÿ“Œ 0

@jung-gt.bsky.social I know itโ€™s academic, but it doesnโ€™t *have* to be made up. These companies make crazy profits margins. United healthcare has a 22-26% profit margin. Itโ€™s not like theyโ€™d go under without the subsidies. Theyโ€™re choosing to extort their customers because healthcare is necessary.

16.11.2025 03:27 โ€” ๐Ÿ‘ 11    ๐Ÿ” 0    ๐Ÿ’ฌ 0    ๐Ÿ“Œ 0

Hence the mauve alignment to see if things look fishy.

But could also just be barnap not picking it up.

16.11.2025 03:21 โ€” ๐Ÿ‘ 1    ๐Ÿ” 0    ๐Ÿ’ฌ 1    ๐Ÿ“Œ 0

Iโ€™ve seen some recent HiFiasm isolate data from a collaborator that produced some confusing misassembly, like duplicated sequences not present in a separate Flye assembly, stretches of sequence with no reads mapping to them.

16.11.2025 03:21 โ€” ๐Ÿ‘ 0    ๐Ÿ” 0    ๐Ÿ’ฌ 1    ๐Ÿ“Œ 0

@stevenjrobbins is following 20 prominent accounts