01  Cell Genomics · 17 August 2026_

Defining and cataloging variants in pangenome graphs

Pouria Salehi Nowbandegani, Shenghan Zhang, Haoyang Hu, Heng Li, Luke J. O'Connor

  • Pangenome graph
  • Structural variation
  • Human genetics
  • Methods

Defines what a genetic variant is inside a pangenome graph, then catalogues 29.6 million of them in the draft human pangenome, 3.5 million of which have a reference allele that is not in GRCh38.

Summary

Structural variation makes some human haplotypes align badly to the linear reference, which is reference bias. A graph reference fixes that but breaks the definition of a variant, because there is no single coordinate to call against. The authors define variants against a reference tree: it keeps every node of the pangenome graph but only a subset of its edges, and any non-reference edge is a variant edge. Applied to the Minigraph-Cactus draft human pangenome they find 29.6 million variants. 3.5 million of them, 11.7 percent, have a reference allele absent from GRCh38, so they are effectively invisible without a graph. Those sit in tangled multiallelic regions. HLA-A and RHD are worked through in detail, showing thousands of small variants entangled with structural variants. They release pantree and a VCF catalogue.

What to take from it

  • The reference tree is the useful idea: a spanning subset of graph edges that restores a single coordinate system without flattening the graph.
  • 11.7 percent of variants have no GRCh38 reference allele. That is the size of the blind spot a linear reference leaves.
  • Small variants and structural variants are not separable in HLA-A and RHD; they are entangled and have to be called together.
  • Output is a VCF catalogue, so the graph result plugs into existing tooling rather than requiring a new stack.

Review

This is the paper to read this week. The variant-definition problem has been the quiet obstacle to using pangenome graphs in practice: everyone agrees the graph is a better reference, and nobody could say what a variant call against it means. A reference tree is a clean answer, and shipping a VCF makes it adoptable. The 3.5 million figure is the number to remember when arguing that a graph earns its cost.

doi:10.1016/j.xgen.2026.101327 · PubMed 42607688

02  Trends in Genetics · 14 September 2026_

African genomes are not a subset

Wubalem D. Seifu, Zewdu E. Bedada, Changqing Zeng

  • Population genetics
  • Human genetics

A forum piece arguing that African genomes making up roughly 1 percent of major database entries, while Africa holds more human genetic variation than the rest of the world combined, is a scientific accuracy problem and not only an equity one.

Summary

Africa harbours more genetic variation than the rest of the world combined, yet about 1 percent of genomes in major databases come from individuals of African ancestry. The authors argue this is a measurement error rather than a fairness gap: reference panels and models built on that skewed sampling distort drug dosing, degrade polygenic risk scores, and weaken precision medicine everywhere, not only for African populations.

What to take from it

  • Roughly 1 percent representation against the largest share of human diversity is the whole argument in one number.
  • Recasting under-representation as measurement error changes who has a stake in fixing it.
  • Direct consequence for pangenome construction: a graph built from the same skewed sample inherits the same blind spot, at higher cost.

Review

Short and worth the five minutes. The framing move is the content. For anyone choosing assemblies to put into a pangenome graph, it is a reminder that sample choice, not graph construction, sets the ceiling on what the graph can represent.

doi:10.1016/j.tig.2026.08.008 · PubMed 42736083

03  Gut Microbes · 15 September 2026_

Pangenome analysis reveals both niche-specific specialists and microbial side hustlers in colorectal cancer microbiomes

Justin Lee, Samuel Minot, Neelendu Dey

  • Pangenome
  • Microbiome
  • Cancer

Scores 3,711 co-associated gene bins across 13 gut species for how well they track phylogeny, separating traits confined to one lineage from traits scattered across distant lineages.

Summary

Species-level associations between the gut microbiome and colorectal cancer replicate poorly between studies. The authors propose that what is conserved is function, not taxonomy, and test it by measuring monophyly of 3,711 co-associated gene bins across 13 species in three cohorts. Both cancer-associated and health-associated bins span the full range of monophyly. In one species with a mixed relationship to cancer, health-associated bins were significantly more lineage-conserved than cancer-associated ones, consistent with cancer-linked traits being mobile. Across species, convergence shows up in Type IV secretion systems, TonB-dependent receptors and RagB/SusD nutrient uptake proteins. Some pairs of gene bins co-occur at a cancer probability above 90 percent.

What to take from it

  • Monophyly score is the instrument: it turns lineage-conserved versus mobile into a number you can rank gene bins by.
  • Lineage-conserved traits point to strain-targeted biomarkers; mobile traits point to therapeutic targets robust to who the patient is.
  • Interactions between gene bins carry signal that species-level analysis cannot see.

Review

A good example of pangenome thinking used as an ecological instrument rather than a catalogue exercise. The gene-bin-pair result is the part that would not have been found any other way. Proof of principle, three cohorts, so treat the effect sizes as provisional.

doi:10.1080/19490976.2026.2728331 · PubMed 42740607

04  Vaccine · 31 August 2026_

An immunoinformatics-based multi-epitope vaccine candidate confers cross-protection against two Actinobacillus pleuropneumoniae serovars

Libin Tan, Yilin Fang, Peifan Liu, Mengmeng Su, Xinkun Zhao, Buyun Xu, Chenxi Li, Xiaoyue Li, Yuan Wang, Wei Zhang

  • Pangenome
  • Vaccine
  • Veterinary

Uses a pangenome core to pick vaccine antigens: 11 candidates screened from the core genome, 7 immunoreactive, fused into one multi-epitope protein that partly protects mice against two serovars.

Summary

Existing inactivated vaccines against Actinobacillus pleuropneumoniae give poor cross-serovar protection. The authors combined pangenome analysis with subtractive proteomics to screen the core genome and identified 11 candidate antigens; 7 were immunoreactive by ELISA and Western blot. Those plus the ApxI-III toxins were used for T and B cell epitope prediction and assembled into a fusion protein, MVAPP. Docking and immune simulation predicted a response; in mice MVAPP raised specific antibodies, shifted splenic T cells toward CD4, and gave partial protection against challenge from two serovars.

What to take from it

  • The core genome is the natural antigen pool when you want cross-serovar coverage - this is what a pangenome is for.
  • Partial protection in mice, two serovars. A candidate, not a vaccine.
  • Subtractive proteomics after the core call is what removes host-homologous candidates.

Review

Standard reverse vaccinology, but a clean illustration of a pangenome collection doing real work: you need the core to know which antigens every serovar carries. No graph needed anywhere in this pipeline.

doi:10.1016/j.vaccine.2026.129104 · PubMed 42673711

05  Computational Biology and Chemistry · 26 August 2026_

Comparative genomics of ESKAPE pathogen species: integrating pan-genome architecture, antimicrobial resistance, and virulence factor repertoires

Srishti Singh, Pankaj Kumar

  • Pangenome
  • Amr
  • Comparative genomics

Builds six species pangenomes from 20 genomes each and asks whether resistance genes sit in the accessory genome. The answer is species-specific, not general.

Summary

120 genomes, 20 per ESKAPE species, quality-checked with CheckM2, pangenomes with Roary, resistance genes with AMRFinderPlus, virulence against VFDB via DIAMOND. Pangenome sizes ran from 4,720 to 17,272 genes, with Enterobacter and Pseudomonas holding the largest accessory genomes. Multidrug resistance was present in 93.3 percent of strains. After FDR correction, resistance genes were enriched in the accessory genome of Enterobacter, Enterococcus, Klebsiella and Staphylococcus, but not Acinetobacter or Pseudomonas. Accessory genome size correlated with resistance class burden in Staphylococcus, Enterococcus and Enterobacter only. Virulence repertoires varied widely, highest in Pseudomonas and lowest in Enterococcus.

What to take from it

  • The accessory-genome-carries-resistance story holds in four of six species and fails in two. Worth knowing before assuming it.
  • 20 genomes per species is small for a pangenome; the accessory curve will not be saturated.
  • FDR correction within each test family is the reason several apparent correlations do not survive.

Review

Useful as a negative result. The interesting line is Pseudomonas and Acinetobacter showing no enrichment in either compartment, which argues their resistance is chromosomal and regulatory rather than acquired. The sample size caps how far that can be pushed.

doi:10.1016/j.compbiolchem.2026.109359 · PubMed 42659948

06  Virulence · 10 August 2026_

Characterizing the pangenome-based virulome and resistome of uropathogenic Escherichia coli from Ecuador

Isaias Mejia-Limones, Ricardo David Avellan-Llaguno, Gabriel Morey-Leon, Derly Andrade-Molina

  • Pangenome
  • Amr
  • Epidemiology

142 uropathogenic E. coli genomes from Ecuador: 3,385 core genes against a 13,347-gene accessory genome, with ST131 at 42 percent.

Summary

Whole-genome sequencing of 142 uropathogenic E. coli isolates from Ecuador gave a 16,732-gene pangenome, split into a 3,385-gene core and a 13,347-gene accessory set. ST131 accounted for 42.3 percent of isolates and serotype O25:H4 for 42.96 percent. The resistome mapped onto the accessory repertoire with 3,712 significant associations across 29 resistance genes. The virulome sorted into 23 profiles, dominated by adhesion and siderophore genes.

What to take from it

  • Core is 20 percent of the pangenome. That ratio is the standard picture for E. coli and a reason the species needs a pangenome at all.
  • ST131 at 42 percent means one high-risk global lineage dominates local disease.
  • 3,712 gene-to-resistance associations is a screening result, not a causal one.

Review

A regional surveillance paper, valuable mostly because Ecuador has not been sampled much. The pangenome here is the analysis frame rather than the finding. Read the core/accessory ratio and the ST131 share; the association counts need replication.

doi:10.1080/21505594.2026.2707776 · PubMed 42576353

07  Food Microbiology · 13 July 2026_

Genomic insights into Salmonella enterica serovar Corvallis from China 2014-2023

Shaofei Yan, Ke Wu, Meng Zhang, Dongmei Tan, Li Zhan, Bosheng Li, Xiaorong Yang, Zushun Yang, Weiwei Li, Yunchang Guo, Changwei Lei

  • Pangenome
  • Amr
  • Food safety

260 Chinese isolates plus 598 global genomes trace a multidrug-resistant Salmonella Corvallis clade moving between the UK and China, with resistance genes riding specific plasmid replicons.

Summary

Salmonella Corvallis has emerged in China with an elevated resistance profile. 260 isolates from the national foodborne surveillance programme 2014-2023 were analysed with 598 public genomes. Isolates from the chicken production chain carried significantly more resistance genes than human isolates. Correlation networks tied replicons to gene sets: IncHI2 with aadA16, dfrA27 and sul1; IncQ1 with floR; IncA/C2 with dfrA12 and aadA2, with IS26 likely mediating co-mobilization. SNP distances ran 0 to 350, a substantial share at 10 or below, indicating transmission. Strains clustered around China, the UK and the Americas, with ST1541 dominant.

What to take from it

  • Poultry chain isolates carry more resistance genes than human ones, which points at where selection is happening.
  • The replicon-to-gene-set pairings are specific enough to be usable as markers.
  • IS26 as the co-mobilization mechanism explains why these genes travel as blocks.

Review

Solid One Health surveillance. Pangenome dynamics is one component among several rather than the point. The UK-China link inferred from SNP distance is the claim that would most benefit from isolation dates and travel data.

doi:10.1016/j.fm.2026.105236 · PubMed 42668177

08  Food Research International · 13 July 2026_

Lineage-aware comparative genomics shows no consistent reservoir-specific plasmidome in low-moisture food associated Salmonella Enteritidis

Kinley Tenzin, Suhan Bheemaiah Balyatanda, Shivaprasad Doddabematti Prakash, Kaliramesh Siliveru, Thomas G. Platt

  • Pangenome
  • Food safety
  • Negative result

Tests whether dry-food Salmonella Enteritidis carries a distinctive plasmid set and finds it does not: plasmid architecture tracks lineage, not reservoir.

Summary

Outbreaks of Salmonella Enteritidis from low-moisture foods are increasing, and plasmids carry stress tolerance traits, so the hypothesis was a reservoir-specific plasmidome. 2,824 non-clinical records from NCBI Pathogen Detection were organised into 51 SNP clusters. Comparisons were restricted to mixed clusters, using 64 assemblies. IncF replicons were near-universal, and no replicon family was enriched in low-moisture isolates after FDR-adjusted, lineage-stratified testing. MOB-suite reconstructed 88 plasmids into 11 backbone clusters, one dominating. Plasmid counts, mobility and function tracked SNP cluster rather than source. Pangenome analysis of 4,899 gene clusters likewise found no accessory genes enriched by reservoir.

What to take from it

  • Restricting comparisons to mixed clusters is the methodological move: it stops population structure from manufacturing a reservoir signal.
  • A clean negative on both plasmid and accessory-gene levels.
  • 24 low-moisture assemblies is small, so absence of evidence here is limited evidence of absence.

Review

Worth reading for the design more than the result. Lineage stratification before testing is exactly what most reservoir-association studies skip, and it is probably why this one found nothing. Negative results at this quality are rare in food genomics.

doi:10.1016/j.foodres.2026.120003 · PubMed 42637401

09  Plant Molecular Biology · 10 September 2026_

Pan-genome characterization of the maize 4CL gene family and its dynamic responses to abiotic stress

Jingbing Zhao, Lingshuang Ren, Lili Li, Dongnan Shao

  • Pangenome
  • Plant genomics
  • Structural variation

Across 26 maize inbred lines, 13 4CL genes split into nine core and four near-core, with structural variants tied to expression and to changes in conserved protein domains.

Summary

Using a maize pangenome of 26 inbred lines, the authors identified 13 Zm4CL genes, nine core and four near-core, in three clades. Ka/Ks indicated mostly purifying selection with some members more relaxed. Structural variation was significantly associated with expression of Zm4CL2 and Zm4CL3, and sequence comparison suggested SVs also altered conserved protein domains in some genotypes. Cold stress raised 4CL enzyme activity at 12 hours; heat, salt and alkali caused an initial drop then recovery; drought had no significant effect. RT-qPCR confirmed the cold and drought expression changes.

What to take from it

  • Core versus near-core across 26 lines is a presence/absence result a single reference genome cannot produce.
  • SVs linked both to expression level and to domain structure, so the variation is not only regulatory.
  • Enzyme activity and transcript level do not move together under drought, which is a caution against reading expression as function.

Review

A conventional gene-family paper upgraded by using a pangenome instead of one reference. That upgrade is becoming standard in crop genomics and this is a representative example. The SV-to-domain link is asserted from sequence comparison, not tested, so treat it as a hypothesis.

doi:10.1007/s11103-026-01760-4 · PubMed 42720865

10  Plants · 26 August 2026_

Genus-wide pan-genome analysis of Populus bZIP transcription factors with reanalysis of public salt-stress transcriptomes

Qi Lin, Yayu Guo, Hou-Ling Wang

  • Pangenome
  • Plant genomics
  • Transcription factors

19 Populus genomes, 1,764 bZIP proteins sorted into 79 orthologous groups, and a reanalysis of public salt-stress RNA-seq that finds no consistent response across taxa.

Summary

The authors analysed 19 Populus genomes, keeping 1,764 bZIP proteins including 21 new loci and four annotation corrections. 1,762 were assigned to 79 orthologous gene groups: 43 core, 20 soft-core, 15 shell, one cloud, with 59 showing copy-number variation. Whole-genome and segmental duplication accounted for 81 percent of assigned proteins and mostly built the conserved component. 72.2 percent of bZIP proteins overlapped a transposable element within the gene body or 2 kb flanks, slightly below matched non-bZIP genes, and copy-number-variable groups showed no greater TE coverage. 97.6 percent of comparisons had Ka/Ks at or below 1. Across the public salt-stress RNA-seq datasets, no orthologous group showed a significant, directionally concordant response in two or more taxa.

What to take from it

  • 21 new loci and four annotation corrections found by looking across 19 genomes instead of one - the practical payoff of a pangenome.
  • The TE result is a negative: transposable elements do not explain which groups vary in copy number.
  • No concordant salt response across taxa. Either the public datasets are too heterogeneous or the response is genuinely lineage-specific.

Review

The annotation corrections are the quietly useful part. The salt-stress reanalysis is honest about finding nothing, and the authors name dataset heterogeneity as the likely cause rather than claiming lineage specificity. Read the core/soft-core/shell breakdown; skip the motif figures.

doi:10.3390/plants15172600 · PubMed 42739367

← All reading notes