I get in a bit of a huff at the vocal experimentalist crowd that seem to want to discount genomics studies as a blanket cultural norm without individual consideration.
10.10.2025 16:45 โ ๐ 1 ๐ 0 ๐ฌ 0 ๐ 0@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.
I get in a bit of a huff at the vocal experimentalist crowd that seem to want to discount genomics studies as a blanket cultural norm without individual consideration.
10.10.2025 16:45 โ ๐ 1 ๐ 0 ๐ฌ 0 ๐ 0Like, try to tell me a whole class of bacteria all maintain that amount of co-localized machinery without ever using it. Itโs possible, but incredibly unlikely. And in fact, SIP studies on DOM degradation were done after this that showed the same taxa we proposed to degrade DOM.
10.10.2025 16:42 โ ๐ 0 ๐ 0 ๐ฌ 1 ๐ 0But I agree, everything becomes much more interesting/defensible when you focus on genomes rather than individual genes. At least then I can say โthis group is enriched in whole classes of CAZymes, co-localized on the genome in PULs, with transporters for the degrading products, etc.โ
10.10.2025 16:27 โ ๐ 1 ๐ 0 ๐ฌ 1 ๐ 0Which is a baffling thought process to me. Gene presence and phylogeny are like 99.9% correlated. Itโs like, โis the organism selected because of its genes or are the genes enriched because the organism is?โ To which I think, whatโs the difference?
10.10.2025 16:23 โ ๐ 1 ๐ 0 ๐ฌ 1 ๐ 0Actually though, the genes vs organisms thing has been brought up by reviewers of multiple papers in the form of โhow do you know those genes arenโt enriched simply because those organisms are enriched? Is it really because the genes are important?โ
10.10.2025 16:20 โ ๐ 0 ๐ 0 ๐ฌ 1 ๐ 0But I agree. I think at its core, I just think a lot of genomics is treaty superficially and without a starting goal. We try to think of these things like in-silico experiments with expected outcomes based on current theory, rather than โI wonder what gene go up/down?โ
10.10.2025 16:17 โ ๐ 0 ๐ 0 ๐ฌ 1 ๐ 0Thanks man! Iโm actually proud of this one specifically. We got desk rejected 3 times from glam journals with a perfunctory โitโs only genomics, we donโt careโ from editors not in the field. Our ISME reviewers in the field had things nice enough to say I wanted to frame the comments.
10.10.2025 16:14 โ ๐ 1 ๐ 0 ๐ฌ 1 ๐ 0โFinishing a complete giraffe genome from telomere to telomere with Verkko-Filletโ walks you through the process of finishing a complete T2T genome using a new Python toolkit developed by Juhyun Kim that makes Verkko output more accessible... www.biorxiv.org/content/10.1...
10.10.2025 15:26 โ ๐ 7 ๐ 2 ๐ฌ 1 ๐ 0Wonโt get any complaints from me on that. Lots of people out there with sequence data making stretched claims based on cursory analysis of a few isolated genes.
10.10.2025 15:30 โ ๐ 0 ๐ 0 ๐ฌ 0 ๐ 0This was a great read. I feel similarly.
09.10.2025 16:52 โ ๐ 7 ๐ 3 ๐ฌ 1 ๐ 0Also throws away a lot of potential to identify new functions in organism with atypical pathways because weโre relying on whatโs been seen before in lab rat microbes, but I donโt think pure genomics should typically be the domain of proposing new pathways anyway. Leave that for wet lab experiments.
09.10.2025 16:01 โ ๐ 1 ๐ 0 ๐ฌ 1 ๐ 0Then if you go to the KEGG page for each KO, it lists experimental evidence for function.
Still room for error, proteins could still have promiscuous function, but at least if we require that most pathways/module genes are present, with key genes, narrows down to a much more trustworthy set.
Thatโs actually a great question. Itโs why we required the โkeyโ genes, and >70% of the pathway to be present, because many genes belong to multiple pathways. But for example, carbon fixation pathways tend to have key genes that arenโt part of other pathways,so we required those to be present.
09.10.2025 15:53 โ ๐ 0 ๐ 0 ๐ฌ 1 ๐ 0Accepted version of our article "Half of microbial eukaryote literature focuses on only twelve human parasites" now out in ISMEJ!
Awesome effort led by undergrad Joanna Lepper and with @hbrappap.bsky.social.
Perfect, thanks Florian! Should have just assumed Anvio would have something. Such a Swiss Army knife these days.
08.10.2025 15:05 โ ๐ 0 ๐ 0 ๐ฌ 0 ๐ 0hear*
07.10.2025 16:36 โ ๐ 0 ๐ 0 ๐ฌ 0 ๐ 0Unfortunately, the power of EnrichM was the consideration of KOs into modules/pathways and doing the stats on the different DBs with standard outputs. Really wish it got published but Joel got an industry job. Maybe one can hack the backend DBs, but my guess is the formats have changed.
07.10.2025 16:01 โ ๐ 1 ๐ 0 ๐ฌ 1 ๐ 0Thanks @jcamthrash.bsky.social and @acritschristoph.bsky.social! DRAM I know of, though it's slowwww from what I here--default to single threads for things that could be multi-threaded. Not sure if that's been fixed. There's also Metacerberus for faster annotation.
07.10.2025 15:59 โ ๐ 1 ๐ 0 ๐ฌ 2 ๐ 0I think that's why KEGG is so powerful, it ties its genes into machine-readable functional units so you're not looking at them in isolation. My opinion, anyway.
07.10.2025 15:56 โ ๐ 0 ๐ 0 ๐ฌ 1 ๐ 0Like here our criteria was that >70% of genes in a an enriched pathway had to be present, including the "key" genes required for proper functioning. One enriched gene says little, but we try to make it as unlikely that we're focusing on statistical accidents.
academic.oup.com/ismej/articl...
Yeah, EnrichM just had something like a Bonferoni correction applied to the p-values. And like any gene/based analysis like that, Iโd only consider it valid if, say, >50% of genes in a module or pathway were enriched. We donโt come up with big stories based on a single gene.
06.10.2025 20:24 โ ๐ 1 ๐ 0 ๐ฌ 2 ๐ 0In EnrichM, you could supply two sets of MAGs, it would annotate with KEGG/PFAM/CAZY, and give outputs showing which KOs and KEGG modules, PFAMS, and CAZY groups were statistically enriched in one set over another, giving a p-value for each KO/PFAM/Module and fold enrichment between the two groups.
06.10.2025 18:19 โ ๐ 0 ๐ 0 ๐ฌ 0 ๐ 0Curious, if someone wanted to annotate a large set of MAGs with KEGG, PFAM, CAZY, and then do some statistical comparison, how would you do that these days?
Back in the day I loved EnrichM, but the backend databases for KEGG, etc, are woefully out of date now.
To fill in a bit, the reason I think The Expanse is the GOAT is because itโs got an epic story arc, super interesting world building, but itโs not dense, has great character development. I get exhausted by sci-fi thatโs trying to look smart with the sci but donโt have characters I truly care about.
05.10.2025 18:34 โ ๐ 2 ๐ 0 ๐ฌ 0 ๐ 0โ ๏ธNew paper โก๏ธ isolating >600 freshwater microbes using dilution-to-extinction and low-nutrient, lake-mimicking media.
๐คฏExpands cultured diversity to ~70 % of detected genera and enables physiological studies of streamlined oligotrophs.
๐ www.nature.com/articles/s41...
#protistsonsky
Team, I need a new book/series to read. Iโm partial to Sci-Fi or fantasy, points for anything that โchanged your lifeโ or got you super reved.
Books/authors I love:
The Expanse series (personal GOAT)
The Bobiverse
Anything Kurt Vonnegut or Neil Gaiman
Already simul-reading 4 non-fiction.
Historically, viruses were thought to primarily use host cell's translational machinery. New work from @harvardcellbio.bsky.social faculty Amy Lee reveals that a giant DNA virus encodes its own IF4F initiation complex, suggesting an unexpected evolutionary innovation. www.biorxiv.org/content/10.1...
02.10.2025 18:10 โ ๐ 93 ๐ 39 ๐ฌ 0 ๐ 2Systematic comparative benchmarking of computational methods for the detection of transposable elements in long-read sequencing data www.biorxiv.org/content/10.1... ๐งฌ๐ฅ๏ธ๐งช
02.10.2025 15:55 โ ๐ 10 ๐ 2 ๐ฌ 0 ๐ 1Interested in aligning your long sequences or small genomes against huge reference databases containing millions of prokaryotic genomes ? A new tool has been released that can do this efficiently - LexicMap - www.nature.com/articles/s41...
01.10.2025 18:47 โ ๐ 13 ๐ 8 ๐ฌ 0 ๐ 0Theyโre coming yโallโฆthe megaplasmids are coming
(I say this in a slightly frustrated manner as someone who has been thoroughly unsuccessful in getting funding to study megaplasmids for the better part of 15 years)