Sign in

iiinova is part of the Atmosphere — an open, public network where you own your identity.

Li Song

@mourisl.bsky.social
264 followers 370 following 10 posts

毛利元光. Assistant Prof at the Department of Biomedical Data Science at Dartmouth College. Research on bioinformatics, algorithms. Lab page: mourisl.github.io

PostsRepliesMedia
Reposted by Li Song
Alessio Campanelli @ale-campa.bsky.social · 20/04/2026
[Construction] With an experimental pipeline building directly to m-Fulgor, we indexed 98% of the AllTheBacteria v0.2 dataset (all species excluding unknown and dustbin) in 4 days. The resulting index takes just 130GB, for 1.8M genomes. Figure: basic stats for some indexes. (4/6)
A table showing the base statistics for four datasets: ATB/no dustbin, no unknown; ATB/Dustbin; ATB/Salmonella Enterica; and ATB/Escherichia Coli. For each dataset, the table provides the size in GB, the number of references, and the number of kmers. The ATB/no dustbin, no unknown dataset is the largest at 130.51 GB with over 1.8 million references, while the ATB/Dustbin dataset has the highest kmer count at over 55 billion.
143 4
Reposted by Li Song
Rob Patro @robp.bsky.social · 17/04/2026
QCatch is now published in Bioinformatics (academic.oup.com/bioinformati...)! Great work from Yuan and Dongze for quality control and analysis downstream of simpleaf/alevin-fry (taking advantage of its structured AnnData output). Give it a try: github.com/COMBINE-lab/...
academic.oup.com
QCatch: A framework for quality control assessment and analysis of single-cell sequencing data
AbstractMotivation. Single-cell sequencing data analysis requires robust quality control (QC) to mitigate technical artifacts and ensure reliable downstrea
0184 18
Reposted by Li Song
ace-gtdb.bsky.social @ace-gtdb.bsky.social · 15/04/2026
GTDB release 11 based on RefSeq 232 (R11-RS232) is live at gtdb.ecogenomic.org. This release covers 901,341 genomes (23% increase) and has 199,923 species clusters (39% increase). Release notes at: forum.gtdb.ecogenomic.org/t/announcing.... Release statistics at: gtdb.ecogenomic.org/stats/r232.
gtdb.ecogenomic.org
GTDB - Genome Taxonomy Database
The Genome Taxonomy Database (GTDB) is an initiative to establish a standardised microbial taxonomy based on genome phylogeny.
15331 53
Reposted by Li Song
Steven Salzberg @stevensalzberg.bsky.social · 04/06/2025
Introns have to come from somewhere, right? @celineh2ooo.bsky.social and I looked at multiple genome alignments with 1000s of genomes and found 342 cases where humans (and our relatives) had gained a new intron. Still not sure where these come from, but it's a fascinating question
24212 42
Reposted by Li Song
Heng Li @lh3lh3.bsky.social · 30/05/2025
Neng Huang developed longcallR for joint SNP calling and phasing from long RNA-seq reads, AND for identifying allele-specific splicing/junctions (ASJ). Although ASJs of statistical significance are rare, a large fraction involve unannotated junctions. In Rust!
0167 16
Reposted by Li Song
Ben Langmead @benlangmead.bsky.social · 29/05/2025
Industry friends, now is the time for MUCH more speaking out on behalf of academic colleagues under duress. Here are core open source methods that many of your products doubtlessly depend on either directly or indirectly (see en.wikipedia.org/wiki/HMMER) being abruptly defunded. Make noise.
en.wikipedia.org
HMMER - Wikipedia
17448 74
Reposted by Li Song
Jim Shaw @jimshaw.bsky.social · 28/05/2025
Announcing myloasm, a new long-read (ONT R10/PacBio) metagenome assembler that I've been working on during my postdoc in the Heng Li lab (@lh3lh3.bsky.social). myloasm-docs.github.io
myloasm-docs.github.io
myloasm - metagenomic assembly with (noisy) long reads
513378 133
Reposted by Li Song
Vikram Shivakumar @vikramshivakumar.bsky.social · 27/05/2025
Excited to share a new update to Mumemto, scaling MUM and conserved element finding to any size pangenome! Preprint out now w/ @benlangmead.bsky.social. Mumemto scales to the new HPRC v2 release and beyond, and can merge in future assemblies without any recomputation! 1/n
biorxiv.org
Partitioned Multi-MUM finding for scalable pangenomics
Pangenome collections are growing to hundreds of high-quality genomes. This necessitates scalable methods for constructing pangenome alignments that can incorporate newly-sequenced assemblies. We prev...
12715 27
Li Song @mourisl.bsky.social · 27/05/2025
Centrifuger has updated the pre-built index list to include this exciting GTDB new release r226 for taxonomic classification of sequencing data: github.com/mourisl/cent.... There is also a gtdb+refseq human/virus/fungi/contaminants index, hopefully will be useful for human microbiome studies.
030 3
Reposted by Li Song
Rob Patro @robp.bsky.social · 27/05/2025
Great 🧵 by Pierre on the Kaminari paper! In short, Kaminari is a simple and elegant, but highly effective index for approximate colored k-mer queries. The simplicity leads to very fast query, but with accuracy consistent with (or exceeding) best-in-class solutions; a very fun collaboration indeed!
0102 10
Reposted by Li Song
Steven Salzberg @stevensalzberg.bsky.social · 13/05/2025
Bioinformatics folks: check out our @biorxivpreprint on a new, very efficient and accurate system for automated genome annotation, EviAnn, led by my colleague Aleksey Zimin: www.biorxiv.org/content/10.1...
biorxiv.org
Efficient evidence-based genome annotation with EviAnn
For many years, machine learning-based ab initio gene finding approaches have been the central components of eukaryotic genome annotation pipelines, and they remain so today. The reliance on these app...
15525 55
Reposted by Li Song
iedb.bsky.social @iedb.bsky.social · 09/05/2025
Check out our latest collaboration with UniProt, who has integrated over 700,000 experimentally validated epitopes to enhance its protein entries with detailed immune response information. This data is accessible via the UniProt Feature Viewer and API! 💻🔬🧪 #collaboration #immunology #proteins
insideuniprot.blogspot.com
Inside UniProt
Rich Epitope Information Comes to UniProt Mammalian immune responses are mediated by interactions between antigens and immune system compo...
021 2
Reposted by Li Song
Rob Patro @robp.bsky.social · 07/05/2025
The deadline for WABI 2025 has been extended (but is still rapidly approaching) wabiconf.github.io/2025/ * abstract deadline: May 12 (AoE) * paper deadline: May 15 (AoE) Consider submitting your exciting algorithmic bioinformatics work to the WABI conference!
wabiconf.github.io
WABI 2025
WABI Conference on Algorithms in Bioinformatics
01011 11
Li Song @mourisl.bsky.social · 04/05/2025
Forgot to dustmasker the genomes before creating a Centrifuger index and indeed saw some misclassifications. Took a while to figure out and lessons learned... Need to implement a built-in masking step like Kraken2 in case forget doing it in the future..
100 1
Reposted by Li Song
Stephen Turner @stephenturner.us · 29/04/2025
Parsing GTF and FASTA files using the eccLib Library www.biorxiv.org/content/10.1... 🧬🖥️🧪 gitlab.platinum.edu.pl/eccdna/eccLib
biorxiv.org
Parsing GTF and FASTA files using the eccLib Library
Summary: Leveraging the Python/C API, eccLib was developed as a high-performance library designed for parsing genomic files and analysing genomic contexts. To the best of the authors' knowledge, it…
062 6
Reposted by Li Song
Noam Teyssier @noamteyssier.bsky.social · 29/04/2025
Extracting @NCBI SRA files with fasterq-dump can require 17x the size of the accession while decompressing. Our new tool xsra extracts sequences at 5x throughput with significantly less disk usage, built-in compression, and optional BINSEQ outputs github.com/arcInstitute...
github.com
GitHub - ArcInstitute/xsra: An efficient CLI to extract sequences from the SRA
An efficient CLI to extract sequences from the SRA - ArcInstitute/xsra
24015 40
Reposted by Li Song
Zamin Iqbal @zaminiqbal.bsky.social · 29/04/2025
Small update from AllTheBacteria (allthebacteria.org). Assemblies can be bulk downloaded from OSF as before, or you can now get individual assemblies from AWS. We now also have a LexicMap index on AWS, so you can align your favourite gene against 2.4million bacteria (next post for price estimates)
allthebacteria.org
AllTheBacteria
14723 47
Reposted by Li Song
Aaron Quinlan (he/him) @aaronquinlan.bsky.social · 28/04/2025
The Department of Human Genetics at the University of Utah is sponsoring the Rising Stars in Genetics and Genomics symposium! - We are seeking nominations bu June 1. - September 18-19, 2025 - Please share with the star postdocs that you know. docs.google.com/forms/d/e/1F...
15543 55
Reposted by Li Song
Igor Martayan @imartayan.bsky.social · 27/04/2025
The sequence analysis session of #RECOMB2025 is off to a great start with @jimshaw.bsky.social presenting devider, a new algorithm for haplotyping small sequences from long-read sequencing. www.biorxiv.org/content/10.1...
1266 26
Reposted by Li Song
Heng Li @lh3lh3.bsky.social · 26/04/2025
If you want to check if a human gene has copy-number changes or lands in a complex region, try pangene.bioinweb.org. Recently updated with more and better assemblies.
24414 44
Li Song @mourisl.bsky.social · 24/04/2025
Time to build a new index!!
071 7
Reposted by Li Song
Heng Li @lh3lh3.bsky.social · 18/04/2025
Minimap2-2.29 released with the support of short RNA-seq read alignment. More explanation and results here: lh3.github.io/2025/04/18/s...
lh3.github.io
Short RNA-seq read alignment with minimap2
0307 30
Reposted by Li Song
Heng Li @lh3lh3.bsky.social · 18/04/2025
Preprint on hifiasm Nanopore-only assembly. Led by Haoyu Cheng: www.biorxiv.org/content/10.1...
biorxiv.org
Efficient near telomere-to-telomere assembly of Nanopore Simplex reads
Telomere-to-telomere (T2T) assembly is the ultimate goal for de novo genome assembly. Existing algorithms capable of near T2T assembly all require Oxford Nanopore Technologies (ONT) ultra-long reads w...
513977 139
Reposted by Li Song
Rob Patro @robp.bsky.social · 18/04/2025
minimap2 adds support for short read spliced RNA-seq alignment! lh3.github.io/2025/04/18/s...
lh3.github.io
Short RNA-seq read alignment with minimap2
1348 34
Reposted by Li Song
Ragnar {Groot Koerkamp} @curiouscoding.nl · 27/03/2025
New set of thesis figures on pairwise alignment just dropped! - schematic and worked example for many algorithms - alignment modes
Schematic figures showing global pairwise alignment algorithmsA worked example for each algorithmSchemetic figures showing various modes, such as semi-global, local, and extension alignment
2339 33
Reposted by Li Song
Nils Homer @nilshomer.com · 14/03/2025
fqgrep release 1.1.0 now speeds up searching FASTQ files! Thank-you to both Markus Schlegel from @activegroupgmbh.bsky.social for updating seq_io and Nicholas D. Crosbie of grepq for some competition and inspiration. See more: github.com/fulcrumgenom...
github.com
GitHub - fulcrumgenomics/fqgrep: Grep for FASTQ files
Grep for FASTQ files. Contribute to fulcrumgenomics/fqgrep development by creating an account on GitHub.
1125 12
Reposted by Li Song
Ragnar {Groot Koerkamp} @curiouscoding.nl · 14/03/2025
Just published simd-sketch, a crate for fast bucket sketches. It's 7x to 30x faster than BinDash, by using the simd-minimizers crate for fast hashing, and a nearly branch-free implementation. Here's a blogpost with a survey of minhash history & methods, and evals: curiouscoding.nl/posts/simd-s...
Schematic of the mod-bucket algorithm: all k-mer hashes are partitioned into s buckets via their remainder mod s. Then, in each bucket the smallest hash is selected.
1129 12
Reposted by Li Song
Anshul Kundaje @anshulkundaje.bsky.social · 21/02/2025
Disease diagnostics using machine learning of B cell and T cell receptor sequences www.science.org/doi/10.1126/... TL;DR: BCRs ARE ALL YOU NEED! (Well actually .... keep reading) 1/
science.org
Disease diagnostics using machine learning of B cell and T cell receptor sequences
Clinical diagnosis typically incorporates physical examination, patient history, various laboratory tests, and imaging studies but makes limited use of the human immune system’s own record of antigen ...
69725 97
Reposted by Li Song
Rob Patro @robp.bsky.social · 24/01/2025
We just released a new version of piscem (should hit bioconda soon & a new version of simpleaf is imminent). The biggest improvement in this new version is increased mapping speed for scATAC-seq data, with the new version being 30-50% faster than the last release! github.com/COMBINE-lab/...
github.com
GitHub - COMBINE-lab/piscem: Rust wrapper for the next generation (still currently in C++)
Rust wrapper for the next generation (still currently in C++) - COMBINE-lab/piscem
0245 24
Li Song @mourisl.bsky.social · 15/01/2025
My PhD adviser Liliana Florea has developed a Coursera course "Bioinformatics Methods for Transcrptomics". A great resource to learn cutting-edge short- and long-read RNA-seq data analysis techniques: www.coursera.org/learn/bioinf...
coursera.org
Bioinformatics Methods for Transcriptomics
Offered by Johns Hopkins University. This course will cover bioinformatics methods for analyzing transcriptomic RNA sequencing data ... Enroll for free.
12713 27
Reposted by Li Song
Adam Phillippy @aphillippy.bsky.social · 02/01/2025
We may have cut our last Canu release in 2024, but we are starting off 2025 with Verkko2! 🎉 Not only is it 4x faster than Verkko1, this version integrates Hi-C data for both phasing and scaffolding, enabling the automated assembly of acrocentric chromosomes! Preprint: www.biorxiv.org/content/10.1...
58631 86
Li Song @mourisl.bsky.social · 31/12/2024
Wrapping up this year's work with a new release of Centrifuger: github.com/mourisl/cent... ! Finally implemented the abundance estimation module. I also found issues in Centrifuge's quantification procedure, so I add a script that uses "centrifuger-quant" to requantify Centrifuge's output.
github.com
Release Centrifuger v1.0.6 · mourisl/centrifuger
Add the module "centrifuger-quant" for abundance estimation. Fix an issue in creating gtdb's taxonomy tree. (#20 ) Fix a bug that causes hit length more than read length (query length)
054 5
Reposted by Li Song
Nick Borcherding @thehumanborch.bsky.social · 31/12/2024
scRepertoire 2: Enhanced and Efficient Toolkit for Single-Cell Immune Profiling www.biorxiv.org/content/10.1...
24713 47
Reposted by Li Song
Anshul Kundaje @anshulkundaje.bsky.social · 25/12/2024
Our ChromBPNet preprint out! www.biorxiv.org/content/10.1... Huge congrats to Anusri! This was quite a slog (for both of us) but we r very proud of this one! It is a long read but worth it IMHO. Methods r in the supp. materials. Bluetorial coming soon below 1/
biorxiv.org
ChromBPNet: bias factorized, base-resolution deep learning models of chromatin accessibility reveal cis-regulatory sequence syntax, transcription factor footprints and regulatory variants
Despite extensive mapping of cis-regulatory elements (cREs) across cellular contexts with chromatin accessibility assays, the sequence syntax and genetic variants that regulate transcription factor (T...
723189 231
Reposted by Li Song
Camille Marchet ⚡ @camillemrcht.bsky.social · 24/12/2024
Merry lexicographic minimiz... Christmas arxiv.org/abs/2412.17492
arxiv.org
On the number of $k$-mers admitting a given lexicographical minimizer
The minimizer of a word of size $k$ (a $k$-mer) is defined as its smallest substring of size $m$ (with $m\leq k$), according to some ordering on $m$-mers. minimizers have been used in bioinformatics -...
02514 25
Reposted by Li Song
Martin Steinegger 🇺🇦 @martinsteinegger.bsky.social · 23/12/2024
Unicore identifies single-copy protein structures across genomes using Foldseek, bypassing slow structure predictions by utilizing 3Di predictions from ProstT5, enabling rapid phylogenetic inference at the tree-of-life scale. 1/n 📄 www.biorxiv.org/content/10.1... 💾 github.com/steineggerla...
212256 122
Reposted by Li Song
Corey Watson @cdubsig.bsky.social · 18/12/2024
Check out long-read sequencing results in human IGH! HUGE duplications in constant gene regions...and 7 copies of IGHV3-23 on one haplotype! @lossiusandreas.bsky.social @williamlees.bsky.social www.biorxiv.org/content/10.1...
biorxiv.org
Ultra-long sequencing for contiguous haplotype resolution of the human immunoglobulin heavy chain locus
Genetic diversity within the human immunoglobulin heavy chain (IGH) locus influences the expressed antibody repertoire and susceptibility to infectious and autoimmune diseases. However, repetitive seq...
12813 28
Reposted by Li Song
Jim Shaw @jimshaw.bsky.social · 17/12/2024
We built a new tool for disentangling local haplotypes from long-read sequencing: check out devider! github.com/bluenote-157... www.biorxiv.org/content/10.1... 1/5
biorxiv.org
devider: long-read reconstruction of many diverse haplotypes
Reconstructing haplotypes is important when sequencing a mixture of similar sequences. Long-read sequencing can connect distant alleles to disentangle similar haplotypes, but handling se-quencing erro...
15630 56
Reposted by Li Song
Ragnar {Groot Koerkamp} @curiouscoding.nl · 17/12/2024
I extended my fast-minimizers blogpost with the lastest results, providing proper deduplicated minimizers and super-k-mers in 1.65ns/kmer, and canonical kmers in 1.90ns/kmer. Or 5.4s and 6.4s respectively on a human genome. Let me know if you want to use this :) curiouscoding.nl/posts/fast-m...
curiouscoding.nl
Computing random minimizers, fast
Table of Contents 1 Introduction 1.1 Results 2 Random minimizers 3 Algorithms 3.1 Problem statement Problem A: Only the set of minimizers Problem B: The minimizer of each window Problem C: Super-k-mer...
1155 15
Reposted by Li Song
Ben Langmead @benlangmead.bsky.social · 11/12/2024
Very excited to see Movi (by @mohsenzakeri.bsky.social) now out in iScience: www.cell.com/iscience/ful.... Movi builds on the "move structure" pangenome index, a compressed full-text index and close cousin to r-index. Compared to r-index, the move structure is simpler and more cache-efficient.
cell.com
Movi: A fast and cache-efficient full-text pangenome index
Biocomputational method; Classification of bioinformatical subject; Genomic analysis
14417 44
Reposted by Li Song
Jim Shaw @jimshaw.bsky.social · 05/12/2024
Thrilled to receive the award! I highly recommend the Microbiome Virtual International Forum @microbiomevif.bsky.social to folks interested in bioinfo and microbiome, free online seminars with flexible times + very high quality talks.
0113 11
Reposted by Li Song
Rob Patro @robp.bsky.social · 02/12/2024
Noor (@noorpratap.bsky.social) has been working on this for quite some time! It took a bit of work to figure out, but now you can map your scATAC-seq reads directly to the genome and sort and dedup them with alevin-fry-atac. You get a BED file with counts ready for downstream analysis! 1/
43914 39
Reposted by Li Song
Jim Shaw @jimshaw.bsky.social · 19/11/2024
I wrote about the backstory of our recently published metagenomic profiler called sylph (www.nature.com/articles/s41...), partially to celebrate the new migration to bluesky Check out the blog here: jim-shaw-bluenote.github.io/blog/2024/de... -- and I apologize in advance for the lack of brevity :)
jim-shaw-bluenote.github.io
Developing sylph - a look into bioinformatics tool development | Jim Shaw
Jim Shaw's academic website.
110436 104