The study’s evolutionary analyses identified a small number of positive-selection sites in CsOSCA1 and weak codon usage bias. The authors also observed differing selection strategies among orthogroups. These findings characterize the composition and evolutionary patterns of the OSCA gene family in the sampled tea accessions; they do not establish whether the same patterns occur across the full range of cultivated tea.

Expression analysis found that CsOSCA1.1 was preferentially expressed in roots and increased in expression under polyethylene glycol stress. A proximal W-box fragment was reported to mediate this stress response, providing a candidate regulatory element associated with the observed expression change.

In a plant functional experiment, the researchers overexpressed tea CsOSCA1.1 in Arabidopsis. Plants carrying the construct showed enhanced tolerance to polyethylene glycol and sodium chloride compared with the experimental comparison, together with increased antioxidant enzyme activities and improved membrane integrity. These results support a role for CsOSCA1.1 in osmotic and salt-stress responses under the tested conditions, but the experiment did not demonstrate improved stress tolerance in tea crops themselves.

The available evidence is an abstract and an early-access version of a peer-reviewed BMC Plant Biology article that the publisher states remains subject to further edits. The reported conclusions are therefore bounded to the sampled tea accessions, the expression assays, and the Arabidopsis overexpression experiment described in the abstract.

Subject of Research: Genome-wide characterization and functional testing of the tea OSCA gene family

Article Title: Genome‑wide identification, evolutionary dynamics, and functional characterization of the CsOSCA gene family in tea plants

Article References: Wang, L., Liu, D., Sun, M., Jiao, X., Ruan, X., Wu, Q., & Wang, W. (2026). Genome‑wide identification, evolutionary dynamics, and functional characterization of the CsOSCA gene family in tea plants. BMC Plant Biology. https://doi.org/10.1186/s12870-026-09962-0

Image Credits: AI Generated

DOI: 10.1186/s12870-026-09962-0

Keywords: tea plant, CsOSCA1.1, OSCA gene family, osmotic stress, salt stress, gene family evolution, purifying selection, pangenome, gene expression, Arabidopsis, Genome, wide

The work addresses a central limitation of single-genome references: one assembled sequence does not represent the full range of genetic variation across a population. For the Emirati Reference Genome Project, researchers assembled multiple human genomes and combined their haplotypes in a population-specific pangenome graph. The central assembly was diploid and trio-based, while the broader collection included HiFi sequencing-based chromosome-scale assemblies. This design links individual chromosome assemblies to a graph intended to preserve both shared sequence and population-level variation.

Assembly quality was assessed through contiguity, base accuracy and the completeness of chromosome-scale scaffolds. Across the assembled haplotypes, the median contiguity N50 was 151.9 megabases and the median base-accuracy QV score was 59.31. Of all assembled chromosomes, 31.8% were telomere-to-telomere scaffolds, defined in the study as chromosome-spanning sequences flanked by telomeric repeats at both ends. Another 64.3% achieved chromosome-level continuity, comprising a single scaffold that spanned an entire chromosome but retained one or more sequence gaps or lacked a telomeric signature at one or both ends. Thus, the results show substantial chromosome-scale representation, while not every assembled chromosome met the stricter telomere-to-telomere criterion.

The researchers next applied the ERGP pangenome to short-read genotyping in an independent cohort of Emirati participants. This analysis identified structural-variant sites in linkage disequilibrium with phenotype-associated variants from genome-wide association studies. The observation shows that the population graph supported structural-variant inference in this cohort and recovered sites associated with previously reported phenotype-associated loci through statistical linkage. It does not establish that the structural variants cause the associated phenotypes, nor does it demonstrate clinical effectiveness.

The Emirati pangenome complements existing human pangenome resources by incorporating a population-specific collection of chromosome-scale assemblies and a graph that includes regionally enriched variation. The reported results provide a reference foundation for population-specific genomics in the United Arab Emirates. The authors identify potential relevance to population genomics and precision medicine, but the evidence presented here is limited to assembly quality, graph coverage and structural-variant inference in one independent Emirati cohort; no direct clinical benefit was tested or demonstrated.

Subject of Research: Construction and evaluation of a population-specific Emirati human pangenome and diploid chromosome-level reference

Article Title: The Emirati pangenome featuring a diploid chromosome-level reference

Article References: Olbrich, M., Mousa, M., Wohlers, I., Al-Aamri, A., Alnaqbi, H., Alsuwaidi, A. H., Manoharan, H. V., Marzouka, N. A.-D., Bartels, L., Ramachandran, S. E., Thomas, A. A., Gasso, M. P., Almehairbi, M., Alameri, M., Tay, G., Hamoudi, R., Ibrahim, S., Al Ghaithi, N., & Alsafar, H. (2026). The Emirati pangenome featuring a diploid chromosome-level reference. Nature Communications. https://doi.org/10.1038/s41467-026-78776-3

Image Credits: AI Generated

DOI: 10.1038/s41467-026-78776-3

Keywords: Emirati pangenome, human genomics, reference genome, diploid assembly, pangenome graph, structural variation, chromosome sequencing, HiFi sequencing, United Arab Emirates, Emirati, pangenome, featuring

The clinical story behind the study began with a 26-year-old man admitted to Changhai Hospital with fever, night sweats, a five-kilogram weight loss over two weeks, and leg pain that rest would not relieve. He was febrile at 39 degrees Celsius with a racing heart and laboratory markers of severe inflammation and anemia. Transthoracic echocardiography revealed vegetations on his mitral valve, and all four bottles from a blood culture set grew S. sinensis. Because of a documented penicillin allergy, clinicians initiated vancomycin, but the infection had already done serious damage. Surgeons found large vegetations on both mitral leaflets and a perforation of the anterior leaflet, requiring mechanical valve replacement. Notably, the organism was recovered again from the excised valvular tissue, and the patient later recovered after a 28-day hospitalization.

From that single patient, the team obtained three isolates: two from the same blood culture, designated 418A and 418B, and one from the valvular vegetation and necrotic tissue, designated 433. Whole-genome sequencing on an Illumina NovaSeq 6000 platform reached depths exceeding 100-fold for each sample. Single-nucleotide polymorphism analysis showed that the blood isolates differed from each other by zero SNPs and from the valve isolate by just one SNP. That near-clonal relationship matters epidemiologically: the three isolates represent repeated sampling of one infection rather than three independent observations, a point the authors emphasize throughout their discussion of the findings.

To place these isolates in evolutionary context, the researchers assembled a dataset of 202 Streptococcus genomes, combining their three clinical isolates with 199 public reference genomes spanning 20 named species. Orthology analysis with OrthoFinder identified 7,200 orthogroups, of which only 382 occurred in every genome, while 6,739 appeared in some but not all. The pan-genome accumulation curves kept rising as more genomes were added and never plateaued, confirming that Streptococcus species carry enormous accessory gene pools beyond a conserved backbone. This dynamic genome structure, shaped by gene gain, gene loss, and recombination, is exactly the kind of genomic backdrop in which lineage-specific virulence traits can arise and persist.

The phylogenetic reconstruction also helped settle a long-running taxonomic debate. Early analyses based on 16S rRNA and single genes had left S. sinensis precariously positioned between the Mitis, Anginosus, and Sanguinis groups of streptococci. Using alignments of core single-copy orthologous genes, the team showed that S. sinensis occupies a Mitis group-associated lineage and is most closely related to S. cristatus, followed by S. sanguinis and S. gordonii, both common oral streptococci. This placement is consistent with an oral reservoir for the species, an idea supported by clinical reports of odontogenic disease preceding endocarditis and by molecular detection of S. sinensis in subgingival plaque from patients with periodontitis.

The central discovery concerned the cyl-like region. All three clinical isolates carried an identical 11,341-base-pair block of twelve genes on Scaffold10, comprising eleven cyl-like genes and acpC, arranged as cylX-cylD-cylG-acpC-cylZ-cylA-cylB-cylE-cylF-cylI-cylJ-cylK. In group B Streptococcus, the canonical cyl operon produces granadaene, a hemolytic pigment that damages tissue barriers and allows the bacterium to penetrate the human placenta. Gene-specific BLAST verification showed 100 percent nucleotide identity and full-length coverage for every one of the twelve coding sequences relative to a reference genome, and independent read mapping confirmed 100 percent mapped-base coverage with mean depths ranging from about 210-fold to over 700-fold across the isolates. Pfam domain analysis using HMMER supported enzyme- and transporter-associated assignments for eight of the encoded proteins.

Computational structural modeling added a layer of plausibility to these annotations. Using the ESMFold protein language model, the team predicted that the CylG-like and CylZ-like proteins adopt soluble enzyme-like folds, with high confidence scores of 98.3 and 84.3 respectively on the pLDDT scale, while the CylB-like protein showed multiple predicted transmembrane helices consistent with a membrane permease, echoing the hemolysin transporter of group B Streptococcus. The authors are careful to stress that these models visualize sequence annotations only and do not demonstrate protein expression, enzymatic activity, or virulence. Synteny comparisons showed that the local genomic organization in S. sinensis broadly resembles corresponding regions in S. gallolyticus and partially in S. gordonii, while differing from the canonical cyl cluster of S. agalactiae.

Intriguingly, the cyl-like block is not universal within the species. Across the seven available S. sinensis genomes, the complete twelve-gene profile appeared only in the three near-clonal clinical isolates and one public reference genome, forming a closely related clade on the core-gene phylogeny, while the type strain HKU4 and two other reference genomes lacked the block entirely. This lineage-associated distribution could reflect ancestral inheritance, lineage-specific retention, or unresolved gain-and-loss events, but it stops well short of proving a causal link to endocarditis. With only one patient represented, the finding identifies a candidate region for future study rather than a confirmed virulence determinant of the disease.

Phenotypic experiments provided suggestive but preliminary support. On Columbia sheep-blood agar, all three isolates produced visible hemolytic zones, and in a liquid sheep erythrocyte assay their relative hemolysis signals, normalized to a detergent control, ranged from roughly 20 to 45 percent, with a clinical S. agalactiae strain as a comparator. The authors candidly note the limitations: the assay was run once with three technical tubes, lacked a bacteria-only optical blank and a cyl-negative control strain, and therefore cannot establish reproducible erythrocyte lysis or tie the phenotype to the cyl-like genes. On the antimicrobial resistance front, screening against the CARD database painted a reassuring picture: all three isolates carried only a single marker, Spne_parC_FLO, and remained phenotypically susceptible to levofloxacin and the other tested agents, including penicillin G, consistent with the generally favorable beta-lactam responsiveness reported for this species.

The study also retained classic host-interaction genes, including eno, pavA, and plr/gapA, which in other streptococci mediate plasminogen binding, fibronectin attachment, and invasion, while screening found no qualifying matches for canonical oral adhesins such as gspB, hsa, or the antigen I/II family proteins. The authors caution that such database-dependent screens cannot prove gene absence in a poorly characterized species. What the work ultimately delivers is a roadmap: independent isolates from multiple patients, complete genome assemblies, and functional experiments involving transcriptional analysis, mutagenesis, and infection models are now needed to determine whether this quiet oral colonist’s cyl-like arsenal is what allows it, in a small number of people, to breach the bloodstream and destroy a heart valve.

Subject of Research: Comparative genomics of Streptococcus sinensis isolates from infective endocarditis identifying a cyl-like genomic region

Article Title: Comparative genomics identifies a cyl-like region in Streptococcus sinensis isolates from a patient with infective endocarditis

Article References: Ma, W., Zhang, Y., Huang, X., Wan, Y., Deng, C., Lin, J., Xu, L., & Qin, Q. (2026). Comparative genomics identifies a cyl-like region in Streptococcus sinensis isolates from a patient with infective endocarditis. Frontiers in Microbiology, 17, Article 1990134. https://doi.org/10.3389/fmicb.2026.1990134

Image Credits: AI Generated

DOI: 10.3389/fmicb.2026.1990134

Keywords: Streptococcus sinensis, infective endocarditis, comparative genomics, cyl-like region, viridans group streptococci, virulence genes, hemolysis, phylogenomics, pangenome, antimicrobial resistance, whole-genome sequencing, Comparative

To understand why GraTools matters, it helps to understand what a pangenome variation graph actually is. Traditional genomics has relied on a single reference genome, against which the genomes of individuals or varieties are compared. That approach works well enough when the population is closely related to the reference, but it systematically distorts our view of species that carry large amounts of variation. Sequences absent from the reference simply vanish from the analysis, and alignments involving insertions, inversions and other structural changes become cumbersome. A pangenome variation graph, or PVG, sidesteps this problem by representing the genomes of many individuals simultaneously as a graph, in which shared sequences form common paths and variable regions branch into alternative paths. The result is a single data structure that captures the full spectrum of genomic diversity within a species, from the core genome shared by everyone to the dispensable genome present in only some individuals.

The catch has always been usability. Building such graphs has become increasingly tractable thanks to graph construction pipelines, but once the graph exists, working with it is another matter. Existing tools for querying and manipulating PVGs tend to be scattered across different software packages, each with its own input formats, and researchers frequently find themselves converting files back and forth between formats before any analysis can proceed. These conversions are not merely annoying; they consume time, introduce opportunities for error, and often prove computationally expensive when graphs grow to represent dozens or hundreds of genomes. The authors of GraTools identified this lack of integration and efficiency as the central bottleneck hindering the wider adoption of graph-based pangenomics, and they designed their tool specifically to eliminate it.

GraTools takes a deliberately pragmatic approach. Its input is the GFA file, the standard Graphical Fragment Assembly format in which pangenome variation graphs are typically produced, and it works directly on that original file rather than requiring researchers to transform it into intermediate representations. The workflow begins with a one-time graph import, during which GraTools converts the graph data into standard formats, namely BAM and BED, that are already supported by a mature ecosystem of highly optimised bioinformatics tools. This design decision is central to the tool’s philosophy: rather than reinventing indexing and querying infrastructure, GraTools leverages what already works best, storing the graph-derived data in formats that allow efficient storage and rapid querying of even large and complex graph structures.

Once the import is complete, the tool opens up a range of downstream operations that can be executed rapidly and repeatedly. Users can extract subgraphs corresponding to regions of interest, retrieve the underlying DNA sequences in FASTA format, and run comprehensive analyses such as calculating the ratio between the core and dispensable genome, a fundamental statistic in pangenomics that describes how much of the species’ genetic repertoire is universal versus optional. GraTools can also identify group-specific segments, meaning stretches of sequence found exclusively in particular subsets of individuals, which is exactly the kind of signal breeders and population geneticists look for when hunting for markers associated with subspecies identity, adaptation or agronomic traits. Crucially, all of this is done through an intuitive command-line interface that provides real-time execution feedback and a detailed logging system, so users always know what the tool is doing and can audit their analyses afterwards.

Scalability was a clear engineering priority. Pangenome graphs are among the largest data structures in modern genomics, and a tool that grinds to a halt on realistic datasets would be of little use. GraTools therefore employs a modular architecture that supports parallel processing and asynchronous input-output operations, allowing it to take full advantage of modern multi-core computing environments. The authors acknowledge the ISO 9001 certified IRD i-Trop HPC facility in Montpellier and the Core Cluster of the Institut Français de Bioinformatique as the computing resources on which the work was performed, a hint that the tool was developed and tested in a genuine high-performance setting rather than only on toy examples.

Another technically significant feature is GraTools’ handling of coordinate systems. In a pangenome graph, positions can be defined relative to the primary reference genome or relative to alternative genomes embedded within the graph, and switching between these perspectives has traditionally required additional processing. GraTools supports coordinates defined on both the primary reference and the alternative genomes within the graph without requiring the user to re-import the data, which means researchers can move fluidly between reference-based and graph-based views of their regions of interest. The tool’s outputs, meanwhile, are deliberately designed to be compatible with external software, so that results can be visualised or further manipulated using the visualisation and analysis tools researchers already have in their pipelines.

The demonstration case study will resonate with anyone working in crop genomics. The team applied GraTools to an Asian rice pangenome graph comprising 13 accessions, a dataset representative of the scale at which many research groups currently operate. Using this graph, they showed that GraTools could easily extract subgraphs, compute depth statistics across the graph structure, and identify segments specific to rice subspecies. Rice is one of the world’s most important staple crops and a model species for pangenome research, with indica and japonica subspecies carrying substantial differentiated sequence content, so the ability to pull out subspecies-specific segments directly from a graph, without elaborate format juggling, illustrates the practical value of the tool for real biological questions in breeding and diversity studies.

The potential applications extend well beyond rice. Because pangenome variation graphs are now being built for humans, livestock, pathogens and countless plant species, a general-purpose manipulation tool has relevance across population genetics, breeding and even genomic medicine, where graph-based representations are increasingly used to capture human diversity more faithfully than a single reference can. The authors position GraTools as suitable for both biologists and bioinformaticians, and its combination of an accessible command-line interface, transparent logging and standard-format outputs lowers the barrier for experimental groups that lack dedicated graph-genomics specialists. The fact that the software is open source under GNU GPLv3, with documentation hosted at gratools.readthedocs.io, should further encourage adoption and community contributions.

GraTools arrives at a moment when pangenomics is shifting from a specialist pursuit to a mainstream expectation. Major pangenome projects across multiple species have demonstrated that single references underrepresent true diversity, and funders and journals increasingly expect analyses that account for it. What has been missing is the connective tissue between graph construction and graph interpretation, and this is precisely the gap GraTools fills. By unifying subgraph extraction, sequence retrieval, core and dispensable genome statistics and group-specific segment detection behind a single, fast interface built on standard file formats, the tool offers what the authors describe as an interesting alternative to existing approaches for manipulating PVGs, facilitating rapid, efficient and flexible downstream analyses. For a field whose growth has been throttled less by ideas than by plumbing, that plumbing fix may prove to be the most consequential contribution of all.

Subject of Research: A software tool for exploring and manipulating pangenome variation graphs

Article Title: GraTools, a user-friendly tool for exploring and manipulating pangenome variation graphs

Article References: Ravel, S., Marthe, N., Carrette, C., Mohamed, M., Sabot, F., & Tranchant-Dubreuil, C. (2026). GraTools, a user-friendly tool for exploring and manipulating pangenome variation graphs. BMC Bioinformatics. https://doi.org/10.1186/s12859-026-06657-z

Image Credits: AI Generated

DOI: 10.1186/s12859-026-06657-z

Keywords: pangenome, variation graph, bioinformatics, GFA format, core genome, dispensable genome, subgraph extraction, rice pangenome, genomic diversity, open-source software, population genetics, plant breeding

The scale of the underlying analysis is considerable. The researchers assembled and examined 2,737 MAC genomes, drawn largely from publicly available data repositories such as the NCBI and the European Nucleotide Archive, and subjected them to exhaustive pairwise comparison using average nucleotide identity, or ANI. ANI is a standard metric in modern bacterial taxonomy: it measures the fraction of the two genomes that can be aligned with high sequence similarity, averaged across all shared genes. Closely related strains typically share more than 95 percent ANI, while organisms conventionally considered the same species usually exceed 98 percent. Rather than simply adopting these thresholds by convention, the team statistically modeled the distribution of pairwise ANI values across the entire complex, determining where natural gaps in genomic similarity occurred and using those gaps to define species and subspecies boundaries at the 95 and 98 percent levels respectively.

From this analysis, fifteen species-level clusters emerged, and the agreement with independent classification schemes was striking. When the researchers compared their clusters against assignments produced by GTDB-Tk, the Genome Taxonomy Database toolkit that has become a de facto standard for genome-based prokaryotic taxonomy, they found 99.9 percent concordance. This near-perfect alignment between a bottom-up, data-driven clustering approach and an established reference framework suggests that the MAC’s genomic structure is far more orderly than its tangled nomenclature would imply. The problem, in other words, is not that the bacteria defy classification; it is that the names attached to them have not kept pace with what the genomes actually show.

Some of the study’s most consequential findings concern the internal structure of the two best-known members of the complex. Within Mycobacterium avium, the subspecies avium, hominissuis, and silvaticum did not form three separate genomic groups but instead constituted a single genomic continuum, blurring the boundaries that clinical and veterinary microbiology have traditionally drawn between these lineages. Notably, the IS900-positive lineage responsible for paratuberculosis, a chronic intestinal disease of ruminants with debated links to human illness, remained genomically distinct from the rest of the continuum, yet its divergence did not reach the level that would justify subspecies status under the framework’s statistical criteria. This finding will resonate in veterinary medicine, where Mycobacterium avium subsp. paratuberculosis, commonly known as Map, is a pathogen of major economic importance, and it raises the question of whether the current subspecies designations within M. avium reflect biology or historical convention.

The picture within Mycobacterium intracellulare proved equally revealing. The subspecies chimaera, implicated in device-associated infections and notorious for the 2015 heater-cooler unit outbreaks that affected cardiac surgery patients worldwide, displayed deep internal structure, suggesting that it is not a monolithic entity. At the same time, the data supported formal recognition of M. intracellulare subsp. yongonense as a distinct lineage, while genomes assigned to the nominal species M. paraintracellulare clustered squarely within M. intracellulare subsp. intracellulare, effectively dissolving that species into its neighbor. The researchers also found that several recently described species names, including M. bouchedurhonense and M. timonense, lacked genomic coherence altogether, meaning that isolates carrying these labels do not form recognizable genetic groups and the names should be treated with skepticism in future work.

Beyond taxonomy, the study delivered a panoramic view of the complex’s collective gene pool. The MAC pangenome, the total set of gene clusters found across all 2,737 genomes, came to 191,904 gene clusters, a number that continues to grow as new genomes are added. Only 0.8 percent of these clusters constitute the universal core, the genes shared by every member of the complex. This open pangenome architecture is characteristic of bacterial groups with extensive horizontal gene transfer and diverse ecological niches, and it helps explain why members of the MAC can colonize everything from natural water sources and soil to bathroom drains, medical devices, and the lungs of patients with underlying respiratory disease. A tiny shared core combined with an enormous accessory genome means that individual isolates can differ dramatically in gene content even when their core sequences are nearly identical.

Defining the taxonomy was only half of the project. The team recognized that the biggest barrier to bringing whole-genome sequencing into routine clinical and epidemiological practice for the MAC is not the absence of data but the absence of standardized, accessible tools that translate raw genomic similarity into harmonized taxonomic assignments. Laboratories that sequence a MAC isolate today may run different pipelines, apply different thresholds, and arrive at different names for the same organism. To close that gap, the researchers developed a consensus genome-based framework that combines FastANI, a fast alignment-based ANI calculator, with dRep, a widely used genome clustering and dereplication tool, and then went a step further by training machine learning classifiers to reproduce the framework’s delineations automatically.

The machine learning approach relies on k-mer frequency profiles, a technique that represents each genome by counting the occurrences of all possible nucleotide words of a fixed length. This representation captures the overall compositional signature of a genome without requiring computationally expensive whole-genome alignments, making it well suited for rapid classification of large numbers of isolates. Trained on the ANI-based delineations established in the study, the k-mer classifiers achieved near-perfect performance, with a Cohen’s kappa statistic of 0.994, a measure of agreement that corrects for chance and approaches the theoretical maximum of 1.0. In practical terms, the classifiers can reproduce the full genome-based taxonomic assignment of a MAC isolate from its sequence composition alone, with essentially no loss of accuracy relative to the slower reference pipeline.

All of this capability has been consolidated into MAC-Explorer, an open-access web tool available at mycobacteriaexplorer.fr. The tool integrates the k-mer-based classification with insertion sequence screening, allowing users to upload a MAC genome and receive a reproducible, harmonized assignment at both the species and lineage levels. The insertion sequence component is particularly relevant for surveillance, since elements such as IS900 and IS1245 have historically been used as markers for specific pathogenic lineages, and their systematic detection within a unified framework brings molecular epidemiology and taxonomy onto the same footing. Because the tool is web-based and free to use, it lowers the technical barrier for diagnostic laboratories and public health agencies that lack dedicated bioinformatics teams, potentially enabling genome-based identification and lineage tracking of MAC isolates as a routine rather than exceptional practice.

The study, which was supported by Santé publique France and by a Poste d’Accueil from the Assistance Publique-Hôpitaux de Paris in partnership with the Institut Pasteur, arrives at a moment when nontuberculous mycobacterial infections are increasingly recognized as a growing clinical challenge, particularly among patients with cystic fibrosis, bronchiectasis, and other structural lung diseases. By replacing an unstable, convention-driven nomenclature with a statistically modeled, genome-based framework, and by making that framework immediately usable through an open web tool, Poignon and colleagues have provided the kind of shared reference system that diagnostics, antimicrobial management, and outbreak investigation all require. The fifteen species clusters, the clarified status of the M. avium subspecies continuum, the dissolution of several questionable species names, and the near-perfect machine learning classifiers together form a foundation on which future MAC surveillance and research can be built, one genome at a time.

Subject of Research: Genome-based taxonomy and identification of the Mycobacterium avium complex

Article Title: MAC-explorer: bridging genome-based taxonomy and an identification tool for the Mycobacterium avium complex

Article References: Poignon, C., Matondo, M., Gianetto, Q. G., Douché, T., Saffarian, A., Veziris, N., Aubry, A., & Godmer, A. (2026). MAC-explorer: bridging genome-based taxonomy and an identification tool for the Mycobacterium avium complex. Genome Medicine. https://doi.org/10.1186/s13073-026-01783-y

Image Credits: AI Generated

DOI: 10.1186/s13073-026-01783-y

Keywords: Mycobacterium avium complex, nontuberculous mycobacteria, average nucleotide identity, genome-based taxonomy, pangenome, machine learning, k-mer classification, whole-genome sequencing, MAC-Explorer, subspecies delineation, bacterial taxonomy, genomic surveillance

The central problem the team set out to solve is a familiar one in crop genetics. Conventional reference genomes are linear, meaning they represent a single consensus sequence against which all other varieties are compared. Such references capture only a limited slice of the diversity that exists within a species, and variants that fall outside the reference sequence can be missed entirely. Graph-based pangenomes address this limitation by encoding the genomes of many individuals simultaneously, allowing single-nucleotide polymorphisms, insertions and deletions, and large structural variants to be represented side by side. But the very comprehensiveness that makes large graphs powerful also makes them computationally demanding, and in some cases they emphasize structural changes at the expense of smaller variants that nonetheless matter for breeding.

MiniTG1.0 took a deliberately economical approach. Rather than aggregating dozens or hundreds of genomes, the researchers built the graph from just two carefully chosen parents: Solanum lycopersicum TS-RS1, a cultivated tomato, and Solanum pimpinellifolium TS-RS2, its wild relative. The two lineages differ substantially in fruit weight, sugar and acid composition, and disease-related chemistry, which makes the contrasts between them informative. The graph integrated 5.76 million variants spanning SNPs, INDELs, and structural variants across the population the team would go on to analyze, providing a dense but tractable map of the differences separating the cultivated and wild genomes.

To connect that variation to traits, the researchers generated 230 seventh-generation recombinant inbred lines, a population in which chromosomes from the two parents have been shuffled and fixed through repeated selfing. Each line carries a unique mosaic of the parental genomes, so any trait difference among the lines can be traced back to the chromosomal segments inherited from each parent. By combining genome-wide association studies across this population with transcriptomic data, metabolite profiling, and targeted functional experiments, the team created a pipeline that moves from statistical association to experimentally verified gene function within a single framework.

The performance gains over a linear reference were substantial. Joint analysis of the three variant classes explained 18.18 percent more average trait heritability than the linear reference genome, and it improved the mapping of established loci controlling trichomes and growth habit. The graph also placed the H and SELF-PRUNING genes closer to their expected trait signals than the linear reference did, indicating that the compact graph preserved, and in some respects sharpened, the resolution needed to locate genes underlying complex phenotypes. For a breeding-oriented tool, that combination of efficiency and precision is the key selling point.

Fruit weight provided the most striking demonstration. The analysis recovered the well-known fw2.2 region, a classic locus in tomato genetics, and additionally uncovered a previously unknown fruit-weight region the team designated fw6.4. Within this new locus, the gene SlILL6 emerged as a strong candidate. When the researchers overexpressed SlILL6, the resulting plants produced fruits weighing roughly 62 to 63 grams, compared with approximately 101 grams in wild-type plants, a reduction of nearly forty percent. The evidence supports a role for SlILL6 in limiting fruit growth through altered cell expansion, adding a new entry to the short list of genes known to govern tomato fruit size.

The study also extended its reach into fruit chemistry. Metabolite profiling detected 1,258 compounds across the population and identified 128 that accumulated at different levels in the two parental backgrounds. A metabolome-based genome-wide association study then connected variation in chlorogenic acid, an antioxidant phenolic compound, to the gene SlCGT. This connection proved functionally meaningful in the context of disease. Using CRISPR-Cas9 gene editing, the team created knockout lines that accumulated more chlorogenic acid and, when challenged with Botrytis cinerea, the fungus responsible for gray mold, developed smaller lesions. Overexpression lines showed the opposite pattern, accumulating less chlorogenic acid and suffering larger lesions. The researchers further showed that chlorogenic acid directly suppresses fungal growth, closing the loop between gene, metabolite, and phenotype.

The authors emphasized that the study demonstrates how a smaller, strategically designed genome graph can deliver both computational efficiency and biological precision. By pairing a wild parent with a cultivated parent and a segregating population, the framework allowed the team to follow inherited variation directly into fruit traits rather than searching through broad diversity panels alone. The discovery of SlILL6 and SlCGT, they noted, shows how the same pipeline can move from statistical association to experimentally tested function, creating a clearer route from pangenomic data to concrete breeding decisions involving fruit size, quality-related metabolites, and resistance to postharvest disease.

The practical implications for molecular design breeding are considerable. SlILL6 could be evaluated as a target for adjusting fruit size, giving breeders a lever on yield that complements the classical fw2.2 locus. SlCGT offers a route to raising chlorogenic acid levels and strengthening resistance to gray mold, and notably the researchers found this could be done without reducing soluble solids in the tested lines, addressing a common trade-off between defense chemistry and fruit quality. MiniTG1.0 itself supplies a compact analytical resource for tracking small and large variants together, which may help breeders recover useful alleles from wild germplasm more efficiently than current approaches allow.

The researchers are careful to note the boundaries of the work. The findings were generated from two parental genomes, greenhouse-grown recombinant lines, and targeted functional tests, so field trials across diverse environments and genetic backgrounds will be needed to confirm agronomic performance, disease protection, and any trade-offs before these targets enter commercial breeding programs. Even so, the study sketches a template that other crop teams may follow: instead of building ever-larger pangenomes, choose informative parents, build a compact graph, and let a purpose-designed population carry the variation into traits that breeders can actually use. If that template proves portable, the quiet efficiency of MiniTG1.0 may matter as much as the individual genes it revealed.

Subject of Research: Graph-based pangenomics of tomato fruit size, metabolite composition, and fungal disease resistance

Article Title: Tomato pangenome maps genes for fruit size and defense

Article References: Tomato pangenome maps genes for fruit size and defense. (n.d.). Original publication

Image Credits: AI Generated

DOI: Not provided

Keywords: tomato, pangenome, graph genome, fruit weight, SlILL6, SlCGT, chlorogenic acid, gray mold, Botrytis cinerea, recombinant inbred lines, CRISPR-Cas9, molecular breeding

Quinoa’s credentials as a nutritional powerhouse are well established, but the review emphasizes that its agronomic versatility is equally striking. The seeds contain 12 to 23 percent protein, comparable to milk and superior to wheat, maize, and rice, and they deliver all essential amino acids, including lysine, methionine, and threonine, which are typically scarce in true cereals. The crop is naturally gluten-free, making it a lifeline for people with celiac disease, and its leaves, sprouts, and microgreens are emerging as novel plant foods in their own right, with dried quinoa leaves reaching protein levels of 37.7 percent, exceeding even spinach. Beyond human nutrition, the whole plant serves as green fodder or silage, and harvest residues feed livestock, making quinoa a genuinely multi-purpose crop.

The agronomic headline, however, is salt. Quinoa is a facultative halophyte, a rare classification among crops that means it can complete its entire life cycle at salinity levels up to 400 millimolar sodium chloride, equivalent to around 40 decisiemens per meter, conditions that would devastate wheat or rice. Remarkably, its growth is often optimal at moderate salinity of 100 to 200 millimolar NaCl, suggesting an inherent appetite for salt. Field trials in the United Arab Emirates demonstrated that quinoa produced seed yields between 0.7 and 1.05 tonnes per hectare using irrigation water at 14 to 15 decisiemens per meter, comparable to yields in non-saline traditional growing regions. Yet the tolerance has limits: one study identified a soil salinity threshold of roughly 12 decisiemens per meter above which seed yield declines rapidly, and germination remains a particularly sensitive stage.

At the heart of quinoa’s salt strategy lies one of plant biology’s most debated structures: the epidermal bladder cell. These balloon-like cells, which cloak the aerial surfaces of leaves, stems, and inflorescences, can reach volumes up to 1,000 times larger than ordinary epidermal cells. For decades they were regarded as salt glands, dumping toxic sodium and chloride ions away from metabolically active tissue. Supporting evidence came from experiments in which mechanically brushing off the bladder cells produced a salt-sensitive phenotype with increased sodium accumulation in the leaf lamina. Transcriptomic work identified key transporters, including the sodium-selective channel CqHKT1.2 and the anion transporter CqClC-c, proposed to load these ions into bladder vacuoles, with stalk cells acting as traffic controllers directing polar ion transport.

But the story has recently grown more complicated. Mutant quinoa lines that completely lack epidermal bladder cells show no loss of salt tolerance compared with wild-type plants, challenging the sequestration hypothesis. Instead, researchers have proposed that oxalic acid accumulated in the bladders deters arthropod herbivores. Intriguingly, multiple studies indicate that bladder cells accumulate potassium in quantities exceeding sodium even under saline conditions, and a recent study found epidermal bladder cells in some spinach varieties with potassium as the dominant ion. Knocking down a WD40 gene reduced bladder cell density and increased water loss, while removing stem bladder cells, but not leaf ones, reduced transpiration and growth. These findings point toward an evolutionarily conserved role in water management rather than a simple salt-dumping function, and the review makes clear that the debate is far from settled.

Whatever the bladder cells’ exact role, quinoa’s osmotic strategy is well characterized. Unlike glycophytes that burn energy synthesizing organic osmolytes, quinoa preferentially uses inorganic ions as cheap osmotica, with an estimated 80 to 95 percent of leaf osmotic adjustment achieved through accumulation of sodium, potassium, and chloride. Toxic ions are sequestered into vacuoles by tonoplast antiporters driven by proton pumps, while tolerant genotypes excel at retaining potassium in the cytosol, maintaining the high cytosolic potassium-to-sodium ratio essential for enzymatic function. Organic solutes play supporting roles: proline acts as a reactive oxygen species scavenger and stabilizer of subcellular structures, betalain pigments in colored varieties protect the photosynthetic machinery from oxidative damage, and soluble sugars contribute to osmotic balance and stress signaling.

Drought tolerance follows a complementary logic. Quinoa behaves as an isohydric species, rapidly closing its stomata in response to water deficits to preserve leaf water potential, a response mediated by abscisic acid transported from roots to shoots. Although this closure restricts carbon dioxide uptake, quinoa maintains photosynthetic capacity even under severe water deficits, hinting at robust non-stomatal resilience. Root architecture adds another layer: genotypes from arid environments develop faster taproot elongation and longer, coarser root segments than their humid-environment counterparts, a phenotypic flexibility that lets the plant forage for deep soil moisture. The crop also tolerates frost down to roughly minus 8 degrees Celsius depending on developmental stage, using soluble sugars and proline as cryoprotectants and accumulating dehydrin proteins to shield seeds and embryos from desiccation and cold.

The review is candid about quinoa’s vulnerabilities, which matter enormously for its expansion beyond the Andes. Heat is the most serious threat: temperatures exceeding 32 to 35 degrees Celsius can slash yields by 31 to 85 percent depending on genotype and timing, and high temperatures reduce pollen viability by 30 to 80 percent, with the most sensitive window falling 8 to 10 days before flowering. Heat also alters seed composition, changing mineral concentrations of calcium, iron, and zinc in ways that persist after the stress ends. Waterlogging is equally damaging, reducing yields by more than 70 percent when it strikes during anthesis, and comparative studies suggest it is more detrimental than drought or salinity. Preharvest sprouting, driven by weak seed dormancy selected against during domestication, threatens harvests in increasingly erratic rainy seasons. Wild relatives such as Chenopodium berlandieri, which show superior yield stability and pollen viability under extreme heat, are now viewed as valuable resources for introgression breeding.

On the genomics front, the past decade has transformed quinoa from an orphan crop into a genetically well-resourced species. High-quality chromosome-scale reference genomes for the coastal Chilean accession QQ74 and the Real cultivar enabled the identification of TSARL1, a transcription factor whose mutation underlies saponin-free sweet quinoa, and revealed expansions of gene families involved in ion transport and abscisic acid homeostasis. An updated QQ74-V2 assembly placed 90.5 percent of the sequence on 18 chromosome-scale scaffolds and uncovered a 52-megabase pericentromeric inversion on chromosome 3B distinguishing coastal from highland ecotypes. A newly developed Chenopodium pangenome spanning 12 species and all eight known genome types shows that Gypsy retrotransposon expansion drives genome size differences and that roughly 65 percent of gene families in the genus are variable shell genes enriched for environmental adaptation. Functional tools have kept pace: the Apple latent spherical virus system for virus-induced gene silencing has confirmed the roles of sodium transporters CqHKT1 and CqSOS1, and a stable Agrobacterium-mediated floral culture transformation system now permits non-chimeric transgenic plants.

The path forward, the authors argue, lies in integrating these resources with modern breeding technologies. High-throughput phenomics using drones equipped with thermal and hyperspectral sensors can capture dynamic stress responses in the field, while artificial intelligence-driven genomic selection predicts breeding values for complex polygenic traits such as yield and stress tolerance. CRISPR/Cas9 editing, potentially delivered through transgene-free viral systems, could introduce targeted mutations for saponin reduction and heat tolerance, and speed breeding techniques that shorten generation times promise to compress cultivar development timelines dramatically. The review also flags institutional hurdles: quinoa’s absence from the Multilateral System of the International Treaty on Plant Genetic Resources and the Nagoya Protocol’s restrictions hamper cross-border germplasm exchange, and researchers have proposed a Global Collaborative Network on Quinoa to standardize descriptors and phenotyping protocols. If these biological and institutional challenges can be met, the authors conclude, quinoa will serve not merely as a promising pseudo-cereal but as a vital model system for translating halophytic mechanisms into resilient cropping systems, turning the world’s degraded, salt-blighted soils into productive farmland.

Subject of Research: Physiological mechanisms and genomic resources underlying abiotic stress resistance in quinoa

Article Title: Quinoa as a naturally stress-resistant crop: current status and future promises

Article References: Zhang, H., Feng, G., & Feng, Y. (2026). Quinoa as a naturally stress-resistant crop: current status and future promises. Stress Biology, 6(1), Article 12. https://doi.org/10.1007/s44154-025-00283-0

Image Credits: AI Generated

DOI: 10.1007/s44154-025-00283-0

Keywords: quinoa, halophyte, salinity tolerance, drought resistance, epidermal bladder cells, pangenome, genomic selection, CRISPR, heat stress, waterlogging, preharvest sprouting, food security

The research team, working under the ethical approval of the Bangladesh Livestock Research Institute, collected 400 pooled fecal samples from commercial poultry farms across four districts: Dhaka, Tangail, Mymensingh, and Sylhet. The sampling strategy deliberately spanned the three major poultry production systems in the country, with 160 samples each from broiler and layer farms and 80 from Sonali farms, a locally important dual-purpose breed. Each composite sample was gathered from five separate locations within a single farm, then pre-enriched in buffered peptone water and cultured on selective media. Molecular confirmation by polymerase chain reaction targeting the 16S rRNA gene identified E. coli in 280 of the 400 samples, a striking 70 percent detection rate that underscores how pervasive the organism is in these production environments.

Antimicrobial susceptibility testing against twelve antibiotics from nine classes revealed alarming phenotypic patterns. Resistance to ampicillin reached 97.86 percent, tetracycline resistance stood at 97.50 percent, and trimethoprim-sulfamethoxazole resistance was recorded at 91.79 percent, while nearly 69 percent of isolates resisted ciprofloxacin, a fluoroquinolone classified as critically important for human medicine. Only imipenem, a last-resort carbapenem, retained high effectiveness, with 96.43 percent of isolates remaining susceptible. These figures reflect decades of unrestricted over-the-counter antibiotic availability in Bangladesh, where antimicrobials have long been used as growth promoters and prophylactic agents, creating sustained selective pressure that favors the emergence and persistence of multidrug-resistant strains.

From this collection, the researchers selected three multidrug-resistant isolates, one from each production system, for whole-genome sequencing on an Illumina NextSeq 2000 platform. The sequencing generated roughly 8.2 to 9.8 million paired-end reads per isolate, achieving approximately 300-fold coverage. Assembled genome sizes ranged from 4.4 to 5.4 megabases with GC contents between 50 and 51 percent. These three genomes were then compared against 83 publicly available genomes from NCBI GenBank, drawn from poultry, human, and environmental sources across the world, producing a comparative dataset of 86 genomes designed to place the Bangladeshi isolates within a genuinely global context.

The resistance gene analysis revealed dense and varied arsenals. The layer isolate proved the most heavily burdened, harboring blaCTX-M-15, blaOXA-1, blaTEM-1B, qnrS1, sul2, sul3, tet(A), tet(M), and aac(6′)-Ib-cr, a combination conferring resistance across beta-lactams, fluoroquinolones, sulfonamides, tetracyclines, and aminoglycosides through mechanisms including antibiotic inactivation, target protection, target alteration, and efflux. The broiler isolate carried blaCTX-M-27 alongside floR, cmlA7, aph genes, dfrA14, sul2, and tet(A), while the Sonali isolate combined blaCTX-M-15 and blaTEM-1B with mph(A), qnrS1, sul3, and dfrA14. The CTX-M family of extended-spectrum beta-lactamases, particularly CTX-M-15 and CTX-M-27, is among the most clinically consequential resistance determinants worldwide, and its consistent presence in poultry-associated isolates signals a reservoir of genes with direct human health implications.

Country-level comparison using the Jaccard similarity index showed that the Bangladeshi isolates shared the highest average AMR gene profile similarity with isolates from Japan at approximately 0.32, followed by South Korea and Norway at around 0.30, with lower similarity to isolates from Argentina, Cambodia, Vietnam, and India. At the individual strain level, the closest matches were strain 45E from Switzerland with a similarity of 0.55, strain 1500 from Norway at 0.50, and EC10 from Italy at 0.49. While these values are moderate rather than close, they indicate that the resistance gene content of the Bangladeshi poultry isolates overlaps meaningfully with globally distributed strains rather than representing an isolated local phenomenon.

Virulence gene profiling painted an equally detailed picture. All three isolates carried the conserved fimbrial adhesin cluster fimA through fimI and curli-associated genes csgA through csgC, indicating strong colonization capacity. The broiler isolate carried the largest virulence repertoire, including extensive type II and type VI secretion system genes, the invasion-associated ibeB and ibeC genes, a full kps capsule gene cluster, and ompA. The layer isolate showed a reduced set dominated by adhesion and iron acquisition genes, including the yersiniabactin system, while the Sonali isolate displayed an intermediate profile with additional immune modulation genes gndA and gtrA. Shared virulence genes with global isolates were highest for Italy at roughly 110 genes, China and Japan at around 109 and 108, and the United States at 108, while Jaccard similarity of virulence profiles peaked with Norway at 0.63, Vietnam at 0.60, and Pakistan at 0.59. A conserved core of adhesion, motility, iron acquisition, and regulatory genes, including fur, phoP, pmrA, rcsB, and rpoS, was present across all three isolates, reflecting the stable machinery required for host adaptation.

Plasmid analysis added another layer of complexity. The broiler isolate harbored the most elaborate plasmid architecture, carrying Col(pHAD28), Col440I, IncFIB, and IncFIB(pLF82) replicon types, a family of plasmids well recognized as major vehicles for the co-dissemination of resistance and virulence genes. The layer isolate carried Col(MG828), IncFIB, and IncFIC(FII), while the Sonali isolate carried only IncX1, the simplest profile. Plasmid profile similarity with global strains was highest for Denmark at approximately 0.40 and the Philippines at 0.39, and lowest for India, Thailand, the United Kingdom, and Germany, underscoring the heterogeneous and regionally variable nature of plasmid content. Because IncF-type plasmids are notorious for linking multiple resistance genes onto single mobile elements, their prominence raises concerns about co-selection, in which the use of any single antibiotic can maintain an entire multidrug resistance payload.

Pangenome analysis across the three poultry isolates identified 7550 gene clusters, of which 2238, or 29.6 percent, were core genes shared by all three, while the remaining 70.4 percent formed the accessory genome, with the layer isolate alone contributing 5078 unique genes. Scaling up to the full 86-genome dataset revealed a pangenome of 29,366 gene clusters with only 4.6 percent core genes and nearly 80 percent cloud genes, a signature of the famously open pangenome of E. coli, which expands continuously as new genomes are added. Multilocus sequence typing assigned all three Bangladeshi isolates to sequence type ST457, a type also reported in Italy. In the minimum spanning tree, ST457 sat within a central clonal complex connected to globally dominant lineages including ST10, ST167, ST648, and the notorious ST131, which was documented across Argentina, Denmark, Italy, Malaysia, the United Kingdom, and the United States in the comparison dataset.

The authors caution that the small number of sequenced Bangladeshi isolates, the heavily fragmented assembly of the layer isolate, and the absence of simultaneously sampled human and environmental isolates limit direct inference of transmission routes. The genetic similarity with international strains indicates shared ancestry and gene exchange patterns rather than proven transmission. Nevertheless, the study delivers a clear message: poultry production systems in Bangladesh harbor E. coli carrying clinically relevant resistance and virulence determinants embedded in globally connected lineages. Strengthening antimicrobial stewardship in veterinary sectors, improving farm-level biosecurity, and building integrated One Health surveillance that links human, animal, and environmental sampling are highlighted as urgent priorities to interrupt the circulation of these genes before their passage into the food chain and human populations becomes still more efficient.

Subject of Research: Comparative genomics of multidrug-resistant poultry-associated Escherichia coli in Bangladesh

Article Title: Comparative Genomic Analysis of Multidrug‐Resistant Escherichia coli Across Poultry–Human–Environmental Interfaces

Article References: Shanto, M. R. H., Ashab Uddin, A. S. M., Supto, M. S. M., Mahim, N. J., Howlader, M. M. R., Ahmed, S. S. U., & Uddin, M. B. (2026). Comparative Genomic Analysis of Multidrug‐Resistant Escherichia coli Across Poultry–Human–Environmental Interfaces. MicrobiologyOpen, 15(5), Article e70406. https://doi.org/10.1002/mbo3.70406

Image Credits: AI Generated

DOI: 10.1002/mbo3.70406

Keywords: antimicrobial resistance, Escherichia coli, whole-genome sequencing, poultry, Bangladesh, pangenome, blaCTX-M, virulence genes, plasmids, ST457, One Health, genomic surveillance

The core idea behind SVPG is straightforward in concept but demanding in execution. Instead of aligning long sequencing reads against a single reference genome, SVPG uses a pangenome graph — a network of genomic sequences that together represent many individuals — as the backdrop against which structural variants are identified. This design choice matters because variants that are common in some populations are often absent from linear reference genomes. When such variants are forced into alignment against a reference that lacks them, reads spanning those regions can map poorly or not at all, causing true variants to be missed and spurious ones to be called. A pangenome graph, by contrast, provides a natural home for alternative alleles, allowing reads to align to whichever path in the graph best explains them.

Long-read sequencing technologies, which produce reads spanning thousands to tens of thousands of DNA bases, have transformed structural variant discovery because most structural variants are large enough that a single read can bridge the entire event. Long reads, however, bring their own computational complications: base-level errors, uneven coverage and repetitive regions of the genome can all confound alignment and variant interpretation. SVPG is built to work directly with these long-read data, combining the sensitivity of read-based analysis with the representational completeness of a pangenome. The result, according to the study, is a method that performs accurate pangenome-facilitated structural variant detection while simultaneously accelerating a second, related task that has become a growing bottleneck for the field.

That second task is pangenome graph augmentation. As sequencing projects continue to enroll participants from increasingly diverse populations, curators of pangenome resources face a persistent dilemma: when a new sample arrives, should the entire pangenome be rebuilt from scratch, an enormously expensive computational undertaking, or should the new sample simply be analyzed independently, forgoing the benefit of the collective reference? Neither option is ideal. Rebuilding everything from the ground up can take days or weeks of computation on large cohorts, while analyzing new samples in isolation discards precisely the advantage that pangenomes were designed to provide. SVPG tackles this problem head-on by providing a rapid procedure for augmenting an existing pangenome graph with new samples, updating the graph to incorporate novel variants without requiring a complete reconstruction.

The practical implications of this dual capability extend well beyond methodological elegance. Large-scale genomics initiatives, such as the Human Pangenome Reference Consortium and national genome sequencing programs around the world, are generating hundreds of thousands of genomes and continually adding to their reference structures. Every improvement in the speed and accuracy of graph augmentation translates directly into faster turnaround for researchers who depend on these resources. A laboratory that can incorporate a newly sequenced individual into the pangenome within hours rather than weeks can provide variant interpretations to clinical collaborators, population geneticists and evolutionary biologists far more quickly. In a clinical context, where structural variants can underlie developmental disorders, cancer and other conditions, timely and accurate detection has direct consequences for patients.

Underlying the method’s performance is a careful treatment of how reads interact with graph structures. Graph-based alignment is computationally more demanding than alignment to a linear reference because each read may map to multiple possible paths through the network. SVPG navigates this complexity by identifying regions of the graph where a new sample’s reads diverge from existing paths and by distinguishing genuinely novel variation from artifacts of sequencing and alignment. Variants confirmed in this way can then be added as new paths in the graph, enriching the reference with each additional sample while preserving the integrity of previously incorporated data. This incremental philosophy stands in contrast to the rebuild-from-scratch workflows that have historically dominated pangenome construction.

Accuracy in structural variant calling is typically measured against sets of variants that have been independently validated, often through a combination of multiple sequencing technologies, optical mapping and painstaking manual curation. Structural variant discovery has historically suffered from high false discovery rates in part because the signatures of these events — split reads, discordant read pairs, changes in read depth — can be mimicked by other genomic features such as segmental duplications and tandem repeats. By anchoring detection in a pangenome that already accounts for much of the common variation, SVPG reduces the ambiguity that produces such false positives, while its long-read foundation preserves the sensitivity needed to detect rare and sample-specific events that no reference graph could anticipate.

The study’s emphasis on pangenome-facilitated detection also reflects a broader conceptual shift in genomics. For two decades, essentially all genome analysis has been conducted against the same reference framework, originally derived from a small number of donors. That framework has been enormously productive, but it embeds biases: populations and individuals whose genomes differ most from the reference tend to be analyzed least accurately. Pangenome graphs are the field’s principal remedy, and tools like SVPG represent the software layer needed to make that remedy practical at scale. A pangenome is only as useful as the analytical ecosystem built around it, and each new tool that treats the graph as a first-class object rather than a curiosity expands the range of questions that can be asked of diverse genomic data.

For researchers working on nonhuman species, the significance is equally clear. Pangenomes are being constructed for crop plants, livestock, model organisms and pathogens, in many cases with even more dramatic structural diversity than found in humans. Polyploid crops, in particular, present alignment challenges that make structural variant detection exceptionally difficult. A method that can rapidly augment a pangenome graph with new samples and call structural variants accurately against long-read data offers a template that can be adapted across these domains, supporting applications from crop breeding to the surveillance of structural variation in evolving pathogen populations.

The publication of SVPG in Nature Methods signals that graph-based structural variant detection has moved from an aspirational goal to a maturing set of practical techniques. As long-read sequencing costs continue to fall and pangenome resources grow, the bottleneck in genomic analysis is increasingly shifting from data generation to interpretation and integration. Methods that combine sensitivity, specificity and computational efficiency — and that keep reference structures current without prohibitive recomputation — address precisely that bottleneck. If the trajectory of the field holds, approaches like SVPG are likely to become standard components of the genomic toolkit, quietly reshaping how structural variation is discovered, catalogued and understood across the tree of life.

Subject of Research: Pangenome-based structural variant detection and rapid pangenome graph augmentation using long-read sequencing data.

Article Title: SVPG: a pangenome-based structural variant detection approach and rapid augmentation of pangenome graphs with new samples

Article References: Jiang, T., Hu, H., Gao, R., Cao, S., Jiang, Z., Zhou, M., Gao, W., Zhou, S., & Wang, G. (2026). SVPG: a pangenome-based structural variant detection approach and rapid augmentation of pangenome graphs with new samples. Nature Methods. https://doi.org/10.1038/s41592-026-03219-2

Image Credits: AI Generated

DOI: 10.1038/s41592-026-03219-2

Keywords: pangenome, structural variants, long-read sequencing, genome graphs, variant calling, graph augmentation, pangenome reference, genomics, reference bias, Nature Methods, SVPG, pangenome-based

The research team, led by Claudia Pérez-Calles and Jingtao Lilue, with senior authors Ashlee H. Rowe and Thomas M. Keane, built the new genomes using modern long-read sequencing and assembly approaches, producing reference-quality resources for Onychomys and its non-resistant relative. This comparative framework matters because earlier work on grasshopper mice had focused largely on a single gene family: the voltage-gated sodium channels, and in particular Nav1.8, a channel expressed in pain-sensing neurons. Previous studies had shown that structural and functional modifications of Nav1.8 in grasshopper mice block the effect of the painful toxins delivered by their scorpion prey, effectively converting a weapon of deterrence into dinner. Yet important questions remained about whether other molecular systems contribute to toxin resistance, and about how the animals respond physiologically when they are actually exposed to toxins in the wild.

By scanning the newly assembled genomes for signatures of positive selection, the researchers identified Onychomys-specific mutations in another sodium channel gene, Scn3a, which encodes the Nav1.3 channel. Unlike Nav1.8, Nav1.3 is expressed in the central nervous system and in the peripheral sensory system after nerve damage, where it is associated with heightened neuronal excitability. The discovery that grasshopper mice carry their own distinctive variants of this channel suggests that toxin resistance is not the work of a single molecular tweak but of a broader rewiring of the excitability machinery of sensory and central neurons. If Nav1.8 keeps painful toxins from firing up peripheral nociceptors, the newly identified Scn3a mutations may help ensure that any signals that do get through, or that arise from tissue damage during a scorpion battle, are blunted as well. The authors are careful to frame these as candidate adaptations whose precise functional consequences will require further experimental testing, but their recurrence in a gene so central to pain signalling makes them a compelling addition to the story.

The pangenome approach also revealed something entirely unexpected outside the nervous system: an Onychomys-specific tandem duplication of a gene called Cblif, which encodes a glycoprotein crucial for the absorption of vitamin B12 in the gut. At first glance, a vitamin transport gene seems an odd companion to toxin resistance, but the researchers connect it to the grasshopper mouse’s unusual dietary and anatomical specialisation. These rodents are carnivorous in a family dominated by seed eaters, and their stomach morphology is distinctive, characterized by an abundance of parietal cells, the acid-secreting cells of the gastric lining. The study found that Cblif is expressed in these parietal cells, and the researchers propose that the gene duplication may represent an adaptation supporting the species’ meat-rich diet and its modified stomach architecture. In an animal that consumes large quantities of arthropod prey, efficient handling of nutrients and gastric chemistry could be as important to survival as venom resistance itself, and the duplication illustrates how a pangenome can surface copy-number changes that single-reference comparisons routinely miss.

To move beyond static genome comparisons and observe the molecular response to toxin exposure as it happens, the team generated RNA sequencing data from two key sensory structures: the dorsal root ganglion, which houses the cell bodies of sensory neurons that relay pain signals from the body to the spinal cord, and the trigeminal ganglion, which performs the equivalent role for the face and head. Grasshopper mice were exposed to toxic sprays from pinacate beetles, another chemically defended prey item that these rodents subdue with gusto. The transcriptomic readout revealed that in Onychomys torridus, the southern grasshopper mouse, genes associated with transportation processes were upregulated following toxin exposure. This pattern suggests that when the animals encounter defensive chemicals, they mobilize cellular transport machinery, potentially to move toxins or their metabolites away from sensitive tissues, to shuttle protective molecules to where they are needed, or to manage the cellular stress caused by chemical assault.

Taken together, the findings sketch a multi-layered model of toxin resistance in grasshopper mice. At the first layer sit the sodium channel adaptations, with Nav1.8 modifications established by earlier work and Scn3a mutations newly identified here, that alter how pain and danger signals are generated and propagated. At the second layer sits an inducible transcriptional response, in which transport-related genes are dialled up in sensory ganglia after chemical exposure, providing a dynamic defence that complements the hard-wired channel changes. And at a third layer, the Cblif duplication hints at the digestive and metabolic specialisations that accompany a predatory lifestyle, reminding readers that adaptation to toxic prey reshapes far more than the nervous system. The authors conclude that Onychomys-specific Scn3a mutations, together with the upregulation of transport-related genes following toxin exposure, may contribute to the evolution of toxin resistance in these rodents, a formulation that is deliberately measured but nonetheless expands the known genetic substrate of the trait.

The significance of the work extends beyond one charismatic genus of desert mice. Pangenomics, the practice of building reference resources that capture the full spectrum of variation within a group of related organisms, is rapidly becoming the standard for species where a single reference genome conceals important structural differences, such as duplications, inversions and gene losses. The Onychomys study demonstrates how such resources can expose candidate adaptations that would be invisible in a conventional single-genome alignment, from tandem duplications to lineage-specific amino acid changes in physiologically critical genes. It also provides a template for studying evolutionary convergence, since toxin resistance has evolved independently in several lineages that interact with venomous prey, including other rodents, hedgehogs and opossums. Comparative work across these systems, enabled by reference-quality genomes, can help distinguish shared molecular solutions from lineage-specific innovations.

There are also potential biomedical implications. Sodium channelopathies are central to chronic pain, epilepsy and cardiac arrhythmia, and the variants that grasshopper mice have evolved under natural selection represent natural experiments in how channel function can be modified without catastrophic side effects. Understanding precisely how the Onychomys Nav1.8 and Nav1.3 variants alter channel gating, drug binding and toxin sensitivity could inform the design of new analgesics that mimic the mouse’s painlessness. Similarly, the transport genes upregulated after toxin exposure may point to endogenous detoxification and cytoprotective pathways that could be pharmacologically enhanced. The researchers caution that their results are a foundation rather than a finished edifice: functional validation, including electrophysiological characterization of the mutant channels and mechanistic dissection of the transport response, remains to be done.

For now, the grasshopper mouse has yielded another instalment of its secrets, and the picture that emerges is one of layered, genome-wide specialization. High-quality reference genomes for Onychomys species and Peromyscus eremicus now stand as public resources for the community, and the specific discoveries, from Scn3a mutations to the Cblif duplication to the toxin-induced transport response, provide a rich set of hypotheses for the next generation of experiments. What was once framed as a single-gene story about a pain-proof sodium channel has become a systems-level account of how a mammalian predator co-evolves with dangerous prey, rewiring its neurons, its digestive physiology and its stress responses in the process. As the authors and their collaborators continue to probe these adaptations, the howling little hunter of the desert may yet teach medicine a thing or two about the molecular biology of pain.

Subject of Research: Genomic adaptations underlying toxin resistance in predatory grasshopper mice (Onychomys)

Article Title: The Onychomys pangenome reveals the unique molecular adaptations that confer toxin resistance

Article References: Pérez-Calles, C., Lilue, J., Anderson, E., Haggerty, L., Martin, F. J., McCarthy, S. A., Adams, D. J., Thybert, D., Rowe, A. H., & Keane, T. M. (2026). The Onychomys pangenome reveals the unique molecular adaptations that confer toxin resistance. Genome Biology. https://doi.org/10.1186/s13059-026-04288-4

Image Credits: AI Generated

DOI: 10.1186/s13059-026-04288-4

Keywords: grasshopper mouse, Onychomys, pangenome, toxin resistance, Nav1.8, Scn3a, sodium channels, Cblif, gene duplication, scorpion venom, dorsal root ganglion, positive selection