Advertisement · 728 × 90

Posts by Vikram Shivakumar

Thank you! And thanks to all the wonderful people at JHU I’ve had the pleasure of working with during my PhD. And of course, to my advisor @benlangmead.bsky.social and my super supportive thesis committee members @aphillippy.bsky.social and @mikeschatz.bsky.social.

Excited to 🇺🇸 -> 🇬🇧 soon!

1 month ago 11 0 0 0

1/ Excited to share my first first-author preprint from my PhD!

We introduce Perseus, a lineage-aware confidence estimation framework for taxonomic classification in long-read metagenomics.

Preprint: www.biorxiv.org/content/10.6...
Code: github.com/matnguyen/Pe...

1 month ago 14 8 1 0

I wish I had this a few months back while searching for a postdoc 😅

1 month ago 0 0 0 0

Awesome! Thanks for sharing the code too!

1 month ago 0 0 0 0

Does this take into account how many pubs the first author had at publication time or what they have currently, when applying the filter? This is a really nifty measure and tool!

1 month ago 0 0 1 0
Post image

As an alternative to the h-index, I made the Mentorship Index (M-index) to proxy a scientist's contribution to mentoring junior scientists. Ex. M10-index = # last-author publications where the first author had < 10 pubs.

Calculate yours: jef.works/Mentorship-I...
Blog: jef.works/blog/2026/03...

🧵👇

1 month ago 10 3 4 2
Post image

New tool from @alexsweeten.bsky.social to find and classify all your satellites: "AniAnn's: alignment-free annotation of tandem repeat arrays using fast average nucleotide identity estimates"
📄 www.biorxiv.org/content/10.6...
📦 github.com/marbl/anianns

2 months ago 44 19 4 1
Jobs | University of Utah Founded in 1850, The University of Utah is the flagship institution of higher learning in Utah, and offers over 100 undergraduate and more than 90 graduate degree programs to over 30,000 students. Uni...

I am hiring a staff bioinformatician for my new lab at the University of Utah! Please consider applying if you are on the hunt:
employment.utah.edu/salt-lake-ci...

3 months ago 31 32 1 1
HLi Lab - Vacancies Openings

I am looking for a postdoc to develop high-performance algorithms in computational genomics. Email or DM me if interested. For more information, see hlilab.github.io/vacancies. RTs appreciated!

3 months ago 44 64 1 0

(t)rust the process?

4 months ago 1 0 0 0
Advertisement

Really excited to see our new work in scaling Mumemto to any size pangenome published in Genome Research this morning. And right on cue with the great opportunity to present this work at #GI2025 this week.

5 months ago 16 5 1 0
Post image Post image Post image

Nicole Brown gave a fantastic talk on Identifying introgressions across pangenomes with Panagram

It uses k-mer conservation to annotate genomic variation across hundreds of genomes, followed by normalization of k-mer profiles to identify introgression events
github.com/kjenike/pana... #GI2025

5 months ago 9 5 1 0
Figure 1: (A) Anchor-based merging requires a common sequence (red) present in each partition. Multi-MUMs are merged by identifying overlaps between partition-specific matches in the anchor coordinate space, and a uniqueness threshold determines if a MUM is still unique in each partition after truncation. (B) String-based merging enables compu- tation of multi-MUMs between partitions without a common sequence. An example tree (left) is shown, highlighting the use case where partial multi-MUMs specific to internal nodes (starred) can be computed by merging subclade-based partitions up a tree. (right) MUM overlaps are computed by running Mumemto on the MUM sequences, and the uniqueness threshold array ensures overlaps remain unique across the merged dataset. (C) An example Burrows-Wheeler Transform (BWT), matrix (BWM), and Longest Com- mon Prefix (LCP) array, with sequence IDs for each suffix shown (ID). A non-maximal unique match (UM) is shown, and the uniqueness threshold for this match is found us- ing the flanking LCP values. (D) A partial multi-MUM (in blue) is found in all-but-one sequence (excluded in red). Using two anchor sequences (red and orange), all-but-one partial MUMs can be computed using an augmented anchor-based merging method (sec- tion 2.6).

Figure 1: (A) Anchor-based merging requires a common sequence (red) present in each partition. Multi-MUMs are merged by identifying overlaps between partition-specific matches in the anchor coordinate space, and a uniqueness threshold determines if a MUM is still unique in each partition after truncation. (B) String-based merging enables compu- tation of multi-MUMs between partitions without a common sequence. An example tree (left) is shown, highlighting the use case where partial multi-MUMs specific to internal nodes (starred) can be computed by merging subclade-based partitions up a tree. (right) MUM overlaps are computed by running Mumemto on the MUM sequences, and the uniqueness threshold array ensures overlaps remain unique across the merged dataset. (C) An example Burrows-Wheeler Transform (BWT), matrix (BWM), and Longest Com- mon Prefix (LCP) array, with sequence IDs for each suffix shown (ID). A non-maximal unique match (UM) is shown, and the uniqueness threshold for this match is found us- ing the flanking LCP values. (D) A partial multi-MUM (in blue) is found in all-but-one sequence (excluded in red). Using two anchor sequences (red and orange), all-but-one partial MUMs can be computed using an augmented anchor-based merging method (sec- tion 2.6).

Post image Post image

Fantastic talk by @vikramshivakumar.bsky.social Mumemto—Scalable multi-MUM finding for pangenomes
Papers biorxiv.org/content/10.1101/2025.05.20.654611 & doi.org/10.1186/s13059-025-03644-0
Code: github.com/vikshiv/mume...
Very efficient pangenome visualization tool, revealing synteny and variations!

5 months ago 23 12 1 1

Looking forward to lots of great talks from JHU folks at CSHL this week!

5 months ago 4 0 0 0

Now this, undergrads, is how you cold email a professor.

5 months ago 2 0 0 0
Post image

If that’s not enough, we threw in a complete, T2T giraffe genome! Giraffe genomes are pretty cool. Almost all of their chromosomes are Robertsonian fusions of the typically telocentric ruminant chromosomes. 🐄 vs. 🦒...

6 months ago 9 4 2 0

Last week we were in the Washington Post for our characterization of Robertsonian chromosomes. This week we are entering our 10th day of being shut down and all of our research is on hold. To help me feel not-so-bad, here is a thread of some studies we released right before the shutdown 🧵 [1/n]...

6 months ago 48 17 1 0
Figure 1(A) ANI quantifies the similarity between two genomes. ANI can be defined as the number of aligned positions where the two aligned bases are identical, divided by the total number of aligned bases. Historically, ANI was calculated using a single gene family for multiple sequence alignment. Another approach finds orthologous genes between two genomes and reports the average similarity between their CDSs. This method was later extended to whole-genome alignment by identifying local alignments and excluding supplementary alignments with lower similarity. (B) Different ANI tools employ various approaches in calculating ANI values. ANIm, OrthoANI, and FastANI use aligners to identify homologous regions, whereas Mash uses k-mer hashing to estimate similarities. Only alignments with higher similarity represented by green arrows are included in ANI calculations, while red arrows, corresponding to paralogs, are excluded. (C) The proposed benchmarking method evaluates the performance of different tools using both real and simulated data. It assumes that more distantly related species on the phylogenetic tree should have lower ANI similarities. This is measured by calculating the statistics of Spearman rank correlation. We expect a negative correlation between ANI and the tree distance (scatter plot on the right).
https://academic.oup.com/bib/article/doi/10.1093/bib/bbaf267/8160681

Figure 1(A) ANI quantifies the similarity between two genomes. ANI can be defined as the number of aligned positions where the two aligned bases are identical, divided by the total number of aligned bases. Historically, ANI was calculated using a single gene family for multiple sequence alignment. Another approach finds orthologous genes between two genomes and reports the average similarity between their CDSs. This method was later extended to whole-genome alignment by identifying local alignments and excluding supplementary alignments with lower similarity. (B) Different ANI tools employ various approaches in calculating ANI values. ANIm, OrthoANI, and FastANI use aligners to identify homologous regions, whereas Mash uses k-mer hashing to estimate similarities. Only alignments with higher similarity represented by green arrows are included in ANI calculations, while red arrows, corresponding to paralogs, are excluded. (C) The proposed benchmarking method evaluates the performance of different tools using both real and simulated data. It assumes that more distantly related species on the phylogenetic tree should have lower ANI similarities. This is measured by calculating the statistics of Spearman rank correlation. We expect a negative correlation between ANI and the tree distance (scatter plot on the right). https://academic.oup.com/bib/article/doi/10.1093/bib/bbaf267/8160681

Excited to share our EvANI benchmarking workflow, published in Briefings in Bioinformatics doi.org/10.1093/bib/...
Computing average nucleotide identity (ANI) is neither conceptually nor computationally trivial. Its definition has evolved over years, with different meanings and assumptions (1/5)

7 months ago 30 12 1 1
Advertisement

10/10 tool name 👌

7 months ago 3 0 1 0
Anchor-based merging requires a common sequence (red) present in each partition. Multi-MUMs are merged by identifying overlaps between partition-specific matches in the anchor coordinate space, and a uniqueness threshold determines if a MUM is still unique in each partition after truncation. (B) String-based merging enables computation of multi-MUMs between partitions without a common sequence. An example tree (left) is shown, highlighting the use case where partial multi-MUMs specific to internal nodes (starred) can be computed by merging subclade- based partitions up a tree. (right) MUM overlaps are computed by running Mumemto on the MUM sequences, and the uniqueness threshold array ensures overlaps remain unique across the merged dataset. (C) An example Burrows-Wheeler Transform (BWT), matrix (BWM), and Longest Common Prefix (LCP) array, with sequence IDs for each suffix shown (ID). A non-maximal unique match (UM) is shown, and the uniqueness threshold for this match is found using the flanking LCP values. (D) A partial multi-MUM (in blue) is found in all-but-one sequence (excluded in red). Using two anchor sequences (red and orange), all-but-one partial MUMs can be computed using an augmented anchor-based merging method.

Anchor-based merging requires a common sequence (red) present in each partition. Multi-MUMs are merged by identifying overlaps between partition-specific matches in the anchor coordinate space, and a uniqueness threshold determines if a MUM is still unique in each partition after truncation. (B) String-based merging enables computation of multi-MUMs between partitions without a common sequence. An example tree (left) is shown, highlighting the use case where partial multi-MUMs specific to internal nodes (starred) can be computed by merging subclade- based partitions up a tree. (right) MUM overlaps are computed by running Mumemto on the MUM sequences, and the uniqueness threshold array ensures overlaps remain unique across the merged dataset. (C) An example Burrows-Wheeler Transform (BWT), matrix (BWM), and Longest Common Prefix (LCP) array, with sequence IDs for each suffix shown (ID). A non-maximal unique match (UM) is shown, and the uniqueness threshold for this match is found using the flanking LCP values. (D) A partial multi-MUM (in blue) is found in all-but-one sequence (excluded in red). Using two anchor sequences (red and orange), all-but-one partial MUMs can be computed using an augmented anchor-based merging method.

Post image (A) Phylogeny of geographically diverse A. thaliana accessions (Lian et al. 2024), with broad geographical regions colored. Internal nodes are labeled with the coverage of partial multi-MUMs across the leaves of each node. Internal node partial MUMs are computed by merging subtree-based partitions progressively up the phylogeny. (B) Global multi-MUM synteny across the full dataset shown in blue (with inversions in green). Global MUMs are computed by merging all partitions together (representing the root node). Additionally, three geographically distinct subgroups are highlighted and partition-specific multi-MUMs (in purple, with inversions in pink) reveal local structural variation in centromeric regions.

(A) Phylogeny of geographically diverse A. thaliana accessions (Lian et al. 2024), with broad geographical regions colored. Internal nodes are labeled with the coverage of partial multi-MUMs across the leaves of each node. Internal node partial MUMs are computed by merging subtree-based partitions progressively up the phylogeny. (B) Global multi-MUM synteny across the full dataset shown in blue (with inversions in green). Global MUMs are computed by merging all partitions together (representing the root node). Additionally, three geographically distinct subgroups are highlighted and partition-specific multi-MUMs (in purple, with inversions in pink) reveal local structural variation in centromeric regions.

Great talk by Vikram @vikramshivakumar.bsky.social on studying pangenomes and synteny visualization in #WABI25
Github: github.com/vikshiv/mume...
First paper: genomebiology.biomedcentral.com/articles/10....
Second: www.biorxiv.org/content/10.1... #WABI2025

8 months ago 22 8 0 0
Post image

Vikram Shivakumar telling us about "Partitioned Multi-MUM finding for scalable pangenomics" #WABI25! So many kinds of matches!

8 months ago 8 3 0 0
Preview
Protein language models reveal evolutionary constraints on synonymous codon choice Evolution has shaped the genetic code, with subtle pressures leading to preferences for some synonymous codons over others. Codons are translated at different speeds by the ribosome, imposing constrai...

This preprint from Helen Sakharova is one of the coolest things to come out of my lab: “Protein language models reveal evolutionary constraints on synonymous codon choice.” Codon choice is a big puzzle in how information is encoded in genomes, and we have a new angle. www.biorxiv.org/content/10.1...

8 months ago 216 84 6 4

Not saying I agree either way, but one pro for text-based file formats are less dependencies needed for viewing files

8 months ago 0 0 1 0

This is so amazing, thank you!

8 months ago 1 0 0 0
comic doodle of Vikram Shivakumar in a sweater and checkered shirt on a pink gradient background, with various elements of the talk to the left: two old moms pointing at MUMs, below an explanation of what those are (large chunks of the same DNA sequence through the genome), at the bottom a few of the organisms worked on: a tomato, a potato, an arabidopsis weed.

comic doodle of Vikram Shivakumar in a sweater and checkered shirt on a pink gradient background, with various elements of the talk to the left: two old moms pointing at MUMs, below an explanation of what those are (large chunks of the same DNA sequence through the genome), at the bottom a few of the organisms worked on: a tomato, a potato, an arabidopsis weed.

#SciArt doodle of @vikramshivakumar.bsky.social's talk yesterday at the @sangerinstitute.bsky.social on MUMs*

*maximal unique matches in pangenomes, now if you did that on sequenced moms you could do mummoms

8 months ago 18 2 2 0

Excited to share our new preprint on detecting foldback artifacts in long reads with my advisors Matthew Meyerson and @lh3lh3.bsky.social ! Stop by poster C-180 on Wednesday at ISMB/ECCB2025 to learn more and chat!

9 months ago 4 2 0 0
Post image

And of course, the poster itself:

9 months ago 3 0 0 0

If you’re in Liverpool, stop by my poster A217 at ISMB/EECB 2025, and chat about all things pangenomes, MUMs, and alignment (and the Beatles or Oasis-mania)

9 months ago 16 6 2 0
Advertisement

Really excited to see this published! To more mum-finding 🍻

10 months ago 21 7 2 0
We are what we index; a primer for the Wheeler Graph era Talk by Ben Langmead - WABI 2025

🖥️🧬We're thrilled to announce that one of our keynote speakers at #WABI2025 will be the inimitable @benlangmead.bsky.social! wabiconf.github.io/2025/talks/t... Ben's keynote is titled "We are what we index; a primer for the Wheeler Graph era", & it's sure to be a whirlwind tour of full-text indexing!

10 months ago 20 5 1 0