The research team, led by Mamta Masand and Harmeet S. Chawla, assembled the genome of HEAR genotype 13R4502 using PacBio HiFi long-read sequencing, a technology that reads stretches of DNA tens of thousands of bases long with high accuracy. The resulting assembly, named UM_B. napus Ver1, spans roughly 1.03 gigabases—close to the expected size of the rapeseed genome—with about 97 percent of the sequence anchored to the crop’s 19 chromosomes. That near-complete anchoring matters enormously for breeders, because it means almost every gene can be placed in its correct physical location, allowing variation between varieties to be compared chromosome by chromosome rather than scattered across disconnected fragments.

Quality control metrics for the assembly are striking. The benchmarking tool BUSCO, which checks for the presence of universally conserved single-copy genes, found 99.6 percent complete, indicating that the assembly captures essentially the full gene repertoire of the organism. The researchers projected 107,468 protein-coding gene models, a high number that reflects rapeseed’s famously complex genome: B. napus is an allotetraploid, formed by hybridization between the ancestors of turnip and cabbage, and it carries overlapping sets of genes inherited from both progenitor species. Repetitive elements accounted for just over 60 percent of the assembly, a reminder that the majority of a plant genome consists of transposable elements and other repeated sequences that have historically made assembly difficult.

It is precisely this repetitive, duplicated complexity that makes long-read sequencing so transformative. Earlier reference genomes, often built from short-read data, collapsed duplicated genes together or broke through repetitive regions, leaving gaps exactly where biologically important variation could hide. Because long reads span these difficult regions, the Manitoba team’s assembly can distinguish gene copies that short-read approaches would have merged, and it can reveal large structural changes—deletions, insertions, inversions and duplications—that smaller fragments of DNA simply cannot span. Structural variants are increasingly recognized as major drivers of trait differences in crops, even though they are far harder to detect than the single-letter changes that older genotyping methods focused on.

To hunt for such variation, the researchers did something methodologically clever: they compared reads from two different sequencing platforms. Alignments of PacBio HiFi reads were assessed independently against Oxford Nanopore reads, and only structural variants supported by both platforms were retained. This cross-platform filtering is a powerful defense against technology-specific artifacts, because errors that plague one chemistry may not occur in another. The analysis yielded 47,182 structural variants confirmed by both datasets—a remarkable catalog of large-scale genomic differences in this HEAR genotype relative to existing references.

Within that catalog, the team focused on structural variants overlapping genes involved in fatty acid metabolism, the biochemical pathway that determines how much erucic acid the seed accumulates. Sixteen fatty-acid-related candidate variants emerged, affecting members of gene families central to lipid biosynthesis and elongation. Among them were FAD2, which encodes a desaturase that converts oleic acid into linoleic acid, and FAE1, also known as KCS18, a beta-ketoacyl-CoA synthase that catalyzes the chain-elongation step responsible for producing very-long-chain fatty acids like erucic acid. Variants were also found in other KCS family genes, in KAR2 and in ACX family genes, which participate in fatty acid elongation and peroxisomal beta-oxidation. Because these enzymes sit directly on the road from simple fatty acids to erucic acid, structural changes within or near them are prime candidates for shaping the HEAR phenotype.

The researchers strengthened their evidence through additional layers of validation. Comparisons against multiple existing B. napus reference genomes provided independent support for four of the candidate variants, showing that the differences were not simply artifacts of aligning against a single reference. More decisively, polymerase chain reaction genotyping confirmed two of the candidates directly in the lab: a 66-base-pair deletion within FAD2 and a 51-base-pair insertion within FAE1/KCS18 in genotype 13R4502 relative to the widely used reference genome ZS11. PCR validation is a gold-standard check in this field, since it physically amplifies the region in question and shows that the size difference is real, not a computational artifact.

The findings are particularly interesting given what is already known about erucic acid genetics. FAE1 has long been established as a major gene controlling erucic acid content, and FAD2 is central to the opposite trait—the low-erucic, high-oleic oil that dominates edible canola markets. Finding structural variants in both genes in a high-erucic-acid genotype suggests that large-scale DNA rearrangements may complement the well-studied point mutations and copy-number changes known to influence these traits. If confirmed functionally, such variants could become precise molecular markers, letting breeders screen seedlings for the DNA differences rather than waiting years to measure oil composition in field-grown plants.

The authors are careful to note that discovery is not the same as proof. As they state, the functional contributions of these candidate variants to erucic acid accumulation require further transcriptomic, biochemical and genetic validation. Showing that a structural variant changes gene expression, alters enzyme function, or tracks consistently with oil phenotype across breeding populations would be the next steps. Nevertheless, the new assembly fills a genuine resource gap: despite several B. napus reference genomes already available, structural variation linked to erucic acid biosynthesis in HEAR germplasm had remained insufficiently characterized until now.

The broader implications reach beyond one Canadian genotype. UM_B. napus Ver1 is openly accessible, and its combination of contiguity, high completeness and validated variants makes it a template for future pan-genome work in rapeseed, where multiple high-quality assemblies from diverse lines can be compared to capture the full spectrum of genomic diversity. As industry demand for erucic acid grows alongside markets for biodegradable lubricants and bio-based polymers, understanding the structural genomic architecture that drives oil composition could accelerate the development of HEAR varieties with higher yields and tailored fatty acid profiles. What began as a technical exercise in long-read assembly may ultimately help turn a Winnipeg breeding program’s elite lines into the industrial oil crops of the future.

Subject of Research: Structural genomic variation in fatty acid metabolism genes of high-erucic-acid rapeseed revealed by long-read genome assembly

Article Title: Long-read genome assembly of high-erucic-acid Brassica napus reveals candidate structural variants in fatty acid metabolism genes

Article References: Long-read genome assembly of high-erucic-acid Brassica napus reveals candidate structural variants in fatty acid metabolism genes. (n.d.). https://doi.org/10.1186/s12864-026-13469-1

Image Credits: AI Generated

DOI: 10.1186/s12864-026-13469-1

Keywords: Brassica napus, high-erucic-acid rapeseed, genome assembly, PacBio HiFi, long-read sequencing, structural variants, fatty acid biosynthesis, FAE1, FAD2, erucic acid, plant genetics, BMC Genomics

The study, published in Genome Biology, was led by Zhiliang Zhang and Jun Xu, who contributed equally, under the corresponding authorship of Huihui Li at the State Key Laboratory of Crop Gene Resources and Breeding in Beijing. The research team recognized a critical gap in the field. While structural variant detection has been extensively benchmarked on human data—where the tools were largely developed in the first place—the performance and cross-species generalizability of these pipelines in plants had been almost entirely overlooked. Plant genomes present a fundamentally different landscape: they can span many gigabases, harbor extraordinarily high repeat content, and carry structural complexity that far exceeds what is found in most mammalian genomes. A pipeline that performs beautifully on human chromosome 1 might falter badly when confronted with the repetitive expanses of a maize or wheat genome.

To address this problem, the researchers constructed something genuinely valuable: high-confidence structural variant truth sets at base-level resolution for ten representative angiosperm species. Rather than relying on a single variant-calling strategy, they used multiple assembly-based approaches to build these truth sets, triangulating across methods to identify variants that could be trusted with confidence. This multi-strategy approach is important because assembly-based calling itself carries method-specific biases, and by combining several strategies the team could filter out artifacts and retain only variants supported by independent evidence. The truth sets they produced are not merely an internal benchmarking resource—they represent a community asset that other researchers can use to evaluate their own pipelines and validate their own findings.

With truth sets in hand, the team then undertook a systematic evaluation of 91 different HiFi-based structural variant calling pipelines. Each pipeline is a combination of two components: a read aligner, which maps long sequencing reads onto a reference genome, and a variant caller, which interprets the alignment patterns to identify structural variants. The combinatorial space is large—aligners such as Minimap2, NGMLR, and pbmm2 can each be paired with callers such as Sniffles1, Sniffles2, cuteSV1, cuteSV2, and DeBreak, among others—and the choice of pairing matters enormously. The researchers evaluated every combination across all ten species, producing a comprehensive performance matrix that reveals how sensitive and precise each pipeline is for detecting different types of structural variants in different genomic contexts.

The results expose a striking heterogeneity in performance across both pipelines and species. No single pipeline emerged as universally superior, and the ranking of tools shifted depending on which plant genome was being analyzed. This cross-species variability is precisely why the team introduced a novel metric: a stability score that quantifies how well a pipeline generalizes across species. A pipeline that achieves excellent results on rice but mediocre results on barley might score high on individual species but poorly on stability. Conversely, a pipeline that performs consistently well—even if not spectacularly—across all ten species earns a high stability score and becomes a more reliable choice for researchers working with diverse or uncharacterized plant genomes.

Using this framework, the team identified clear winners for different variant classes. For deletion detection—arguably the most common and biologically significant class of structural variant—the combinations NGMLR/Sniffles2 and NGMLR/DeBreak emerged as the optimal cross-species pipelines. These pairings combine the NGMLR aligner, which was designed specifically to handle the characteristic error profiles of long reads, with callers that are particularly adept at identifying missing segments of DNA. For insertion detection, a different set of pipelines proved most robust: Minimap2/Sniffles1, Minimap2/cuteSV2, and pbmm2/cuteSV1. The fact that different aligner-caller combinations excel at different variant types underscores a key insight of the study—there is no single best pipeline, and researchers must choose their tools based on the specific biological questions they are asking.

Beyond the pipeline rankings, the study revealed several important patterns about where and when structural variant detection breaks down. Detection accuracy declined sharply for large structural variants and for variants located within repetitive regions, particularly simple sequence repeats. This finding has profound implications for plant genomics, because many of the most interesting structural variants in crops—those affecting agronomic traits—are precisely the large ones, and they often reside in the repetitive, gene-poor regions that are hardest to analyze. The researchers also discovered that detection precision follows a nonlinear relationship with sequencing depth for each individual haplotype, with F1 scores peaking at approximately 15-fold coverage. This is a critical practical finding: sequencing more deeply does not necessarily improve structural variant calling, and beyond the optimal depth, accuracy may actually decline as the sheer volume of data introduces more opportunities for alignment errors and false positive calls.

The implications of this work extend well beyond the immediate benchmarking results. For plant breeders and geneticists, the availability of validated pipeline recommendations means that structural variant discovery can now be conducted with much greater confidence, enabling more accurate association studies, more reliable pan-genome construction, and better-informed breeding decisions. Structural variants are increasingly recognized as major drivers of crop domestication and adaptation, and having the right analytical tools to detect them is essential for harnessing that variation in breeding programs. The truth sets generated in this study provide a standardized reference against which new tools can be evaluated as they emerge, ensuring that the field does not simply accumulate more data but also improves its ability to interpret it correctly.

For the broader community, this study also serves as a cautionary tale about the dangers of uncritically applying tools developed for one organism to another. The plant genomics community has sometimes borrowed human-focused analytical frameworks without fully validating them for plant-specific challenges, and this benchmark demonstrates that such borrowing can lead to substantially suboptimal results. The stability score framework introduced here could readily be adapted to other non-model organisms, providing a template for rigorous cross-species tool evaluation in fields ranging from fungal genomics to vertebrate biology. As long-read sequencing becomes increasingly affordable and accessible, the bottleneck in structural variant discovery is shifting from data generation to data interpretation—and studies like this one are essential for ensuring that the computational tools researchers rely upon are actually up to the task.

The research was supported by the Biological Breeding-National Science and Technology Major Project, the National Natural Science Foundation of China, the Project of Sanya Yazhou Bay Science and Technology City, the Research Program of Zhongyuan Research Center, and the Agricultural Science and Technology Innovation Program. The full article, including supplementary figures and tables, is available open access in Genome Biology, allowing researchers worldwide to consult the detailed performance data and select the pipelines most appropriate for their own species of interest.

Subject of Research: Benchmarking of long-read structural variant calling pipelines across ten angiosperm species using PacBio HiFi sequencing

Article Title: Benchmarking long-read aligner-caller pipelines for structural variant detection across ten angiosperm species with PacBio high-fidelity sequencing

Article References: Benchmarking long-read aligner-caller pipelines for structural variant detection across ten angiosperm species with PacBio high-fidelity sequencing. (n.d.). https://doi.org/10.1186/s13059-026-04300-x

Image Credits: AI Generated

DOI: 10.1186/s13059-026-04300-x

Keywords: structural variants, PacBio HiFi, long-read sequencing, plant genomics, benchmarking, angiosperms, variant calling, NGMLR, Sniffles2, Minimap2, Genome Biology, crop genetics

Telomere-to-telomere, or T2T, assembly represents the gold standard of modern genomics. Conventional genome assemblies leave behind hundreds or even thousands of unresolved stretches known as gaps, typically in regions crammed with repetitive DNA such as centromeres, ribosomal RNA gene clusters and other satellite-rich zones. These gaps matter. They hide genes, obscure structural variation and make it impossible to say exactly how many chromosomes a species truly carries. A T2T assembly eliminates those holes entirely, reading each chromosome as one continuous thread from the telomere at one end to the telomere at the other. Achieving this for a plant genome is notoriously difficult, since plant genomes are often loaded with repetitive transposable elements that confound sequencing algorithms.

The team overcame the challenge by combining three complementary sequencing strategies. Pacific Biosciences HiFi reads provided highly accurate long stretches of DNA, Oxford Nanopore ultra-long reads extended those stretches across the most stubborn repetitive regions, and Hi-C, a technique that captures the physical contacts between chromosomes inside the cell nucleus, was used to order and orient the sequences into full-length chromosomes. The final assembly spans roughly 890.50 megabase pairs and achieves a scaffold N50 of 44.26 megabases, meaning half of the assembled sequence resides in chromosomes longer than 44 million bases. The entire genome was organized into exactly 20 gap-free chromosomes, each capped by telomeres, with 40 telomeres and 20 putative centromere regions identified in total.

One of the most consequential results is the definitive resolution of the chromosome number of I. cornuta. Botanists have debated the basic chromosome number of this species for years, with conflicting counts circulating in the cytological literature. The new gap-free assembly shows unambiguously that the species carries 20 pairs of chromosomes, a diploid complement written as 2n = 2x = 40. Crucially, the team did not rely on sequencing data alone. An independent karyotype analysis of root-tip cells, in which chromosomes are stained and counted directly under the microscope, confirmed the genomic count. Agreement between the molecular assembly and classical cytogenetics settles a long-standing controversy and provides a firm foundation for evolutionary comparisons across the Aquifoliaceae family.

Genome quality was assessed with the metrics that the assembly community takes most seriously. The annotation, recovered using RNA sequencing data, predicted 35,437 protein-coding genes. Benchmarking with BUSCO, a test that searches for a set of universally conserved single-copy genes, returned a completeness score of 98.33 percent, indicating that almost the entire expected gene repertoire is present. The base-level accuracy of the assembly, expressed as a quality value or QV, reached 48.97, which corresponds to an error rate of roughly one mistake in every 79,000 bases. These figures place the Chinese holly assembly among the most accurate plant genomes published to date and demonstrate that the integration of HiFi, ultra-long and Hi-C data can deliver reference-grade sequences even for species with no prior genomic resources.

With the reference in hand, the researchers turned to evolutionary history. Comparative genomics against related plant species revealed two whole-genome duplication events in the ancestry of I. cornuta, ancient episodes in which the entire genetic complement was doubled and subsequently pared down through gene loss and divergence. Whole-genome duplications are a major engine of plant evolution, because duplicate genes are free to accumulate mutations and take on new functions. The analysis also uncovered lineage-specific expansions of gene families, with a pronounced enrichment in pathways governing secondary metabolism, the biochemical arsenal plants use to synthesize defensive compounds, pigments and medicinal molecules.

These expanded gene families are particularly intriguing given what the Chinese holly is known for. The species thrives under stresses that would defeat many ornamentals, and its tissues contain bioactive compounds of pharmaceutical interest. The authors are careful, however, to frame these expansions as candidate gene sets that could potentially contribute to environmental adaptability. Genomic evidence of expansion does not by itself prove function. Verifying the biological roles of these genes will require transcriptomic profiling across stress conditions and direct functional assays, such as gene knockout or overexpression experiments, work that the T2T reference now makes dramatically easier to design and interpret.

The significance of the assembly extends well beyond a single species. Ilex is a large and ecologically diverse genus of more than 400 evergreen and deciduous trees and shrubs distributed across temperate and tropical regions, and it had never before had a gap-free reference genome of any kind. The I. cornuta assembly therefore serves as a genomic anchor for the entire genus, enabling researchers to map other holly species against a complete, contiguous blueprint, reconstruct chromosome evolution within the family, and pinpoint the genetic basis of traits ranging from sex determination, a famously thorny question in hollies, to fruit color and cold hardiness.

For breeders, the practical implications are immediate. Ornamental holly breeding has traditionally depended on slow phenotypic selection and conventional crossing, approaches that are poorly suited to traits controlled by many genes or by structural variants buried in repetitive DNA. A telomere-to-telomere reference changes that calculus. It allows genome-wide association studies, marker-assisted selection and eventually genome editing to be pursued with confidence that the underlying sequence is complete and correctly ordered. Genes involved in stress tolerance, secondary metabolite production and ornamental quality can now be located precisely, tracked in breeding populations and, in time, targeted directly.

The work also stands as a showcase of how far genome science has traveled. Little more than a decade ago, completing a plant chromosome end to end was deemed unreachable; today, a coordinated pipeline of three sequencing technologies can deliver a gap-free reference for a species whose genome had never been assembled before. The Chinese holly now joins the small but growing roster of plants with fully finished genomes, and the resource, deposited and openly described by the Nanjing-led team, is freely available to researchers worldwide. For a genus that decorates gardens, supports traditional medicine and endures some of the toughest growing conditions in its native range, the holly has finally been given a complete instruction manual, and scientists are only beginning to read it.

Subject of Research: Telomere-to-telomere gap-free genome assembly of the Chinese holly, Ilex cornuta

Article Title: Telomere-to-telomere, gap-free genome assembly of the Chinese holly, Ilex cornuta

Article References: Zhou, P., Xu, K., Gong, W., Li, F., Sun, H., Zhang, Q., & Zhang, M. (2026). Telomere-to-telomere, gap-free genome assembly of the Chinese holly, Ilex cornuta. BMC Plant Biology. https://doi.org/10.1186/s12870-026-10103-w

Image Credits: AI Generated

DOI: 10.1186/s12870-026-10103-w

Keywords: Ilex cornuta, Chinese holly, telomere-to-telomere, genome assembly, centromere, chromosome number, whole-genome duplication, comparative genomics, PacBio HiFi, Oxford Nanopore, Hi-C, molecular breeding

The study, published in BMC Plant Biology, assembled the complete mitochondrial genomes of four species in the genus: Rehmannia glutinosa, the source of the famous herbal medicine dihuang, along with R. henryi, R. piasezkii, and R. chingii. All four belong to the broomrape family, Orobanchaceae, a group better known for its parasitic members but which also contains these therapeutic and horticulturally valuable plants. Despite decades of research into the bioactive compounds and metabolic pathways that make Rehmannia medicinally important, its mitochondrial genomes had remained essentially unexplored territory until now.

What the team found is a striking demonstration of structural plasticity. The mitochondrial genome of R. glutinosa spans 547,002 base pairs and exists as a typical single circular molecule, as does that of R. piasezkii, which is considerably larger at 742,657 base pairs. But the other two species broke the mold entirely. R. chingii, at 662,985 base pairs, is organized as one linear molecule accompanied by two circular structures, while R. henryi, at 663,966 base pairs, is split across four separate circular molecules. In other words, four species from the same genus carry mitochondrial genomes that differ not just in size but in their fundamental molecular architecture.

This kind of multipartite organization is not merely a curiosity. Plant mitochondrial genomes are known to exist as dynamic populations of molecules whose shapes depend on recombination between repeated sequences scattered throughout the DNA. When repeat-mediated recombination occurs at certain sites, a single master circle can be resolved into multiple smaller circles or linear fragments, and different configurations can coexist within the same cell. The assembly approach used in this study relied on long-read PacBio HiFi sequencing, which allowed the researchers to resolve the complex assembly graphs and confirm which alternative connections between contigs were genuinely supported by the read data, a critical step when the genome refuses to conform to a single circular model.

Beneath the architectural drama, the researchers documented a detailed landscape of repetitive DNA. Simple sequence repeats, short tandem motifs that are useful as molecular markers, showed a broadly conserved composition across the four species. However, the numbers of dispersed repeats, longer duplicated sequences that can mediate recombination, and tandem repeats varied quantitatively among the species. This variation in repeat content is exactly what would be expected to drive the differences in genome structure, since each pair of sufficiently large repeats represents a potential recombination point that can reshape the molecule into alternative conformations. The findings thus connect the fine-scale repeat inventory directly to the large-scale structural diversity observed among the four genomes.

Perhaps the most surprising result is how little the functional content of these genomes varies despite their structural chaos. The team predicted 394 C-to-U RNA editing sites across 31 protein-coding genes in each of the four species, a profile that was highly conserved throughout the genus. RNA editing of this kind is a hallmark of plant mitochondria, where specific cytidines in the transcribed RNA are converted to uridines, often restoring codons that would otherwise be nonfunctional. The consistency of this editing landscape suggests that while the physical packaging of the genome is fluid, the information it encodes and the post-transcriptional processing it requires are under strong stabilizing selection.

The study also traced the ongoing conversation between two organelles. Mitochondrial genomes in plants frequently absorb DNA from the chloroplast, and these mitochondrial plastid DNA sequences, or MTPTs, were abundant in all four Rehmannia species, with 21 to 28 such transfers identified in each genome. Among them was a conserved region of approximately 4.5 kilobases containing complete sequences of ndhB, a photosynthesis-related gene, along with two plastid ribosomal protein genes, rps7 and rps12. These fragments are essentially molecular fossils of past DNA traffic from chloroplast to mitochondrion, and their conservation across the genus points to a shared history of inter-organelle gene transfer predating the divergence of these species.

To place these genomes in an evolutionary context, the researchers conducted phylogenetic analysis that strongly supported the monophyly of Rehmannia, confirming that the genus forms a single coherent lineage within Orobanchaceae. They then measured nucleotide diversity and calculated the ratio of nonsynonymous to synonymous substitution rates, known as Ka/Ks, across the shared protein-coding genes. The overall picture was one of conservation: most mitochondrial genes are evolving slowly and under purifying selection, as expected for genes encoding essential components of the respiratory chain. Yet a handful of genes stood out, with sdh3, atp6, and matR showing relatively higher sequence variation among the species, hinting at localized hotspots of divergence within an otherwise conservative genome.

The contrast between these two layers of evolution, structural and sequence-based, is what makes the study compelling. On one hand, the mitochondrial genomes of Rehmannia have diversified dramatically in size and organization over what is presumably a modest evolutionary timescale, with one species carrying its genome in four circles and another in a hybrid of linear and circular molecules. On the other hand, the gene content, RNA editing patterns, and rates of sequence evolution remain remarkably stable. This coexistence of pronounced structural plasticity with conserved genomic features appears to be a defining characteristic of plant mitochondrial genomes, and the Rehmannia data provide a particularly clean example because the four species are close relatives, minimizing confounding effects from deep evolutionary divergence.

For the broader research community, the work delivers a set of complete, well-annotated mitochondrial genomes for a medicinally significant genus, along with a framework for understanding how repeat-mediated recombination, chloroplast DNA influx, and selective constraint interact to shape mitochondrial genome evolution in flowering plants. It also underscores a practical lesson for genome assembly: assuming a single circular mitochondrial chromosome can lead researchers astray, and long-read sequencing with careful validation of alternative assembly paths is increasingly essential. As more species across Orobanchaceae and beyond receive the same treatment, comparisons like this one will help reveal whether the structural exuberance seen in Rehmannia is the rule or the exception in plant mitochondrial evolution, and whether the genes that vary most, such as sdh3 and atp6, hold clues to adaptation in these remarkable medicinal plants.

Subject of Research: Comparative assembly and evolutionary analysis of mitochondrial genomes in four Rehmannia species

Article Title: Assembly and comparative analysis of the four complete mitochondrial genomes in Rehmannia (Orobanchaceae): structural diversity and evolutionary insights

Article References: Jiang, X., Yan, C., Qiu, A., Hou, R., Wu, Q., Huang, S., Wang, J., Ren, M., Wang, S., Li, J., & Liu, Z.-L. (2026). Assembly and comparative analysis of the four complete mitochondrial genomes in Rehmannia (Orobanchaceae): structural diversity and evolutionary insights. BMC Plant Biology. https://doi.org/10.1186/s12870-026-10054-2

Image Credits: AI Generated

DOI: 10.1186/s12870-026-10054-2

Keywords: Rehmannia, mitochondrial genome, Orobanchaceae, genome assembly, RNA editing, MTPTs, chloroplast gene transfer, repeat-mediated recombination, plant evolution, PacBio HiFi, Ka/Ks selection, medicinal plants

The research team, led by Theresa Lüth and colleagues at the Institute of Neurogenetics at the University of Lübeck, together with collaborators at the University of Luxembourg and the Medical University of Innsbruck, set out to solve a persistent technical problem. Short-read sequencing, the current gold standard for mitochondrial variant detection, relies on reads of only about 50 to 300 base pairs. That brevity creates uneven coverage, mapping biases, and vulnerability to contamination from NUMTs, fragments of mitochondrial DNA that have been pasted into the nuclear genome over evolutionary time. Short reads also cannot directly reconstruct a complete mitochondrial haplotype, the full sequence of a single mitochondrial genome, because each read captures only a tiny slice. Previous long-read approaches using Oxford Nanopore sequencing had achieved reliable detection only for variants above roughly 5 percent, and many full-length mitochondrial sequencing protocols relied on multiple overlapping PCR amplicons, which introduce their own amplification biases.

The Lübeck team’s solution was elegantly simple: amplify the entire roughly 16,600-base-pair mitochondrial genome in one single long-range PCR reaction, and then read each complete molecule end to end on the PacBio Vega platform. PacBio HiFi sequencing generates circular consensus reads with Phred quality scores of at least 30, corresponding to 99.9 percent per-base accuracy, which is precisely the fidelity needed to distinguish genuine low-frequency variants from sequencing noise. In their runs, more than 95 percent of all sequenced reads met this HiFi standard, and after filtering, the majority of reads spanned the full-length mitochondrial amplicon at approximately 16.6 kilobases. Average sequencing depth across the 24 samples reached a staggering 83,180-fold coverage, with uniform coverage across every position of the mitochondrial genome, evidence that the single-amplicon approach avoided the amplification bias that plagues overlapping PCR strategies.

Validation was the crux of the study, and the researchers designed a rigorous benchmarking experiment. They prepared predefined mixtures of two mitochondrial haplotypes, one from haplogroup D4e1’3 as the major component and one from haplogroup J1c2 as the minor component, at ratios of 5, 2, 1, and 0.1 percent. Because both samples had previously been characterized with Illumina NextSeq short-read sequencing, the expected variant profile of each mixture was known in advance, allowing the team to calculate sensitivity, precision, and the F1 score, the harmonic mean of the two. The results were striking. For mixtures between 5 and 1 percent, the F1 score was a perfect 1.0, meaning every expected variant was detected and no spurious variants appeared. Detected minor-variant levels closely matched the expected values, averaging 0.057 for the 5 percent mixture, 0.022 for 2 percent, 0.011 for 1 percent, and 0.001 for the 0.1 percent mixture.

Even at the extreme 0.1 percent level, where the researchers deliberately lowered the variant-calling threshold to 0.0005, all expected minor-haplotype variants were recovered, and sensitivity remained at 1.00. Ten additional variants discordant with the Illumina reference did pass the lowered threshold, dropping the F1 score to 0.91, but follow-up analysis suggested these were threshold-dependent noise rather than true signals. Two of them sat in the notoriously ambiguous homopolymeric D-loop region, and read-level assessment showed the discordant alleles were also present at low levels in the 5 percent mixture, simply below its higher calling threshold. Notably, all discordant variants occurred on molecules assigned to the major haplotype, not the minor component, and a graph-based tool called Himito found no evidence of NUMTs contamination anywhere in the dataset. The authors are careful to stress that these exploratory results should not be interpreted as establishing a general analytical detection limit of 0.1 percent.

The depth of sequencing also proved to be a resource that could be traded for throughput. By systematically downsampling the data from roughly 70,000-fold coverage down to 300-fold, the team showed that accuracy and sensitivity remained very high at depths above 1,500-fold for mixtures of 1 percent or more, and the F1 score stayed at 1.00 across all tested depths for the 5 percent mixture. The practical implication is considerable: with the 230 barcodes available on a single PacBio Vega SMRT Cell, each sample would still receive approximately 13,000-fold mitochondrial coverage, ample for accurate detection of heteroplasmy at the 1 percent level. In other words, large cohorts could be multiplexed without sacrificing the low-frequency sensitivity that makes this technique valuable.

Beyond counting individual variants, the long reads allowed something short reads fundamentally cannot: direct reconstruction of complete mitochondrial haplotypes at the single-molecule level. When the researchers classified individual full-length reads according to haplotype-specific marker variants, the observed proportions of minor-haplotype reads closely matched the expected mixture ratios, and visualization of representative reads from the 1 percent mixture showed that informative variants consistently co-occurred on the same sequencing molecule. This phasing capability matters because it confirms that variants arising from the same mitochondrial genome can be linked together, a prerequisite for understanding how mitochondrial populations evolve within cells and tissues.

With the workflow validated, the team applied it as a proof of principle to a biological question at the heart of Parkinson’s disease research. They studied fibroblast cell lines from five carriers of the LRRK2 p.Gly2019Ser variant, the most common genetic cause of familial Parkinson’s disease: three affected by the disease and two unaffected. The cells were cultured in parallel and sampled at passages 1, 4, 8, and 12 over 84 days, and mitochondrial DNA was extracted and sequenced at each time point. Restricting their analysis to variants above 1 percent heteroplasmy, based on the mixture validation, the researchers observed a decrease in the number of detected heteroplasmic variants across passages, confirmed by a linear mixed-effects model showing significantly fewer variants at passage 12 compared with passage 1. In exploratory analyses, the heteroplasmic variant load was higher in fibroblasts from affected carriers than unaffected carriers, and it was associated with older age at sample collection, though the authors caution that the small number of cell lines demands careful interpretation.

Perhaps most intriguing were the variant-specific trajectories the long-read data revealed. In one fibroblast line from an affected carrier, a synonymous variant in MT-CO2 fell from 20.5 to 7.7 percent heteroplasmy across passages, and a synonymous variant in MT-ND5 dropped from 23.3 to 8.3 percent, while a missense variant in MT-ND2 rose steadily from 8.4 to 57.5 percent. These opposing patterns, with some variants declining and others expanding, suggest that selective pressures or clonal expansion processes actively shape mitochondrial heteroplasmy during cell culture, consistent with recent single-cell work showing that selection, rather than random drift, governs heteroplasmy levels in dividing cells. No structural variants such as deletions were detected above the 1 percent threshold in the analyzed samples, and no consistent pattern distinguished affected from unaffected carriers, underscoring the substantial inter-individual variability in heteroplasmy dynamics.

The study is not without limitations, which the authors lay out candidly. Long-range PCR risks primer dropout when polymorphisms fall within primer-binding sites, variants in those regions cannot be reliably assessed, PCR amplification can introduce low-level polymerase errors that the benchmarking design does not fully capture, and the PCR-based enrichment erases native DNA methylation, precluding epigenetic analysis. The validation also focused on single-nucleotide variants rather than structural variants. Still, the achievement stands: this is the first demonstration that PacBio HiFi sequencing can detect mitochondrial heteroplasmy with high accuracy down to 1 percent using a haplotype mixture benchmarking framework, a substantial advance over the roughly 5 percent limit previously reported for Oxford Nanopore sequencing under comparable conditions. As the cost of long-read sequencing continues to fall, this workflow offers researchers a powerful new lens for studying mitochondrial variation in aging, cancer, mosaic somatic disease, and neurodegeneration, and future studies in larger cohorts will determine whether the elevated variant load seen in affected LRRK2 carriers holds up as a genuine molecular signature of Parkinson’s disease.

Subject of Research: High-accuracy detection of mitochondrial DNA heteroplasmy using PacBio HiFi long-read sequencing in Parkinson's disease

Article Title: Long-read sequencing enables high-accuracy mitochondrial heteroplasmy detection in Parkinson’s disease

Article References: Lüth, T., Schaake, S., Much, C., Belyea, M. M., Seibler, P., Grünewald, A., May, P., Klein, C., Weissensteiner, H., & Trinh, J. (2026). Long-read sequencing enables high-accuracy mitochondrial heteroplasmy detection in Parkinson’s disease. BMC Genomics, 27(1), Article 831. https://doi.org/10.1186/s12864-026-13426-y

Image Credits: AI Generated

DOI: 10.1186/s12864-026-13426-y

Keywords: mitochondrial DNA, heteroplasmy, PacBio HiFi, long-read sequencing, Parkinson's disease, LRRK2, variant calling, fibroblasts, haplotype phasing, BMC Genomics, NUMTs, benchmarking

The research team, led by Yusheng Wang and Cunyu Zhou of Yangtze University together with colleagues at Huazhong University of Science and Technology, the Wuhan Botanical Garden of the Chinese Academy of Sciences, and Southwest Minzu University, assembled complete mitochondrial genomes from nine accessions of D. fruticosa and one accession of D. glabra. Rather than relying on short-read sequencing, which often struggles with the repetitive architecture of plant mitogenomes, the team used PacBio HiFi long reads, a technology that produces highly accurate reads long enough to span the complex repeats that make plant mitochondrial DNA notoriously difficult to assemble. When combined with two previously published accessions, the dataset encompassed twelve complete mitogenomes, providing an unusually rich comparative resource for a single plant genus.

The assembled genomes ranged in size from 263,889 to 294,600 base pairs, a span of roughly thirty thousand base pairs that underscores how variable mitochondrial genomes can be even among closely related individuals of the same species. Despite this size variation, the overall architecture proved remarkably stable. Guanine-cytosine content was conserved across all twelve accessions, and the core gene set remained essentially unchanged: thirty protein-coding genes, three ribosomal RNA genes, and between seventeen and eighteen transfer RNA genes. This combination of a stable functional core with a fluctuating genome size is a hallmark of plant mitochondrial evolution, where most of the size differences arise not from changes in gene content but from the gain and loss of non-coding sequence, repeats, and foreign DNA.

Codon usage analysis added another layer to the picture. Across all twelve mitogenomes, the researchers detected a clear preference for codons ending in A or U, a pattern consistent with the nucleotide composition biases that characterize many plant mitochondrial genomes. Such preferences matter beyond mere bookkeeping: they influence translational efficiency, the evolution of synonymous sites, and the choice of markers for phylogenetic reconstruction. For researchers hoping to use mitochondrial genes to untangle relationships within Dasiphora and among its relatives in the Rosaceae, knowing which codons and which genes carry the strongest compositional signals is essential groundwork.

Perhaps the most visually striking findings came from the repeat analysis. The team catalogued abundant simple sequence repeats and dispersed repeats scattered throughout the genomes, including a set of reverse repeats that appear to be specific to the D. glabra lineage. In plant mitochondria, large repeats are not passive passengers; pairs of repeats can recombine with one another, flipping genome segments, generating alternative genomic conformations, and occasionally triggering structural rearrangements. The presence of lineage-specific reverse repeats in D. glabra suggests that this species may harbor distinct recombination dynamics, potentially shaping its mitochondrial genome architecture in ways that differ from D. fruticosa even though the two species share a recent common ancestor.

The study also documented extensive traffic of DNA between the cell’s three genome-containing compartments, a phenomenon that has fascinated plant genomicists for decades. The researchers identified multiple chloroplast-derived fragments embedded within the mitochondrial genomes, evidence of DNA flowing from the plastid into the mitochondrion over evolutionary time. They also detected large nuclear mitochondrial DNA segments, known as NUMTs, in which mitochondrial sequence has been transferred into the nucleus. Strikingly, D. glabra harbored especially large mitochondrial-derived segments, indicating that the scale of inter-compartmental DNA transfer can differ substantially even between sister species. These transfers are more than curiosities; they can disrupt genes at insertion sites, inflate genome size, and complicate phylogenetic analyses by mixing signals from different compartments.

To gauge how much raw evolutionary variation the mitogenomes contain, the team measured nucleotide diversity across the genomes. The analysis pinpointed several highly variable mitochondrial genes as well as a set of hypervariable intergenic regions in the plastid genomes. Such variable loci are exactly what population geneticists and conservation biologists need: they serve as fine-grained markers for tracking gene flow, identifying population structure, and reconstructing the demographic history of species that colonized high-altitude habitats. For a genus like Dasiphora, whose members occupy fragmented alpine landscapes sensitive to climate change, a menu of variable markers is a practical tool for future monitoring and conservation work.

Phylogenetic trees built from the mitochondrial genomes and, separately, from the plastid genomes told a consistent and intriguing story: D. glabra was consistently placed within the D. fruticosa clade rather than as a clearly diverged sister lineage. This result echoes the reticulate, or network-like, evolutionary history that nuclear data had previously suggested for the pair, and it raises questions about how species boundaries are maintained, or blurred, in this group. Whether the pattern reflects recent divergence, incomplete lineage sorting, or historical hybridization cannot be settled by organellar genomes alone, but the congruence between mitochondrial and plastid signals strengthens the case that the two species are separated by a very shallow evolutionary divide.

Collinearity analysis, which aligns the genomes to one another to detect shared blocks of sequence, revealed high structural conservation across the twelve mitogenomes, punctuated by localized rearrangements. In other words, the overall gene order and large-scale architecture have been preserved, but specific regions have been shuffled, inverted, or repositioned, likely through the repeat-mediated recombination described above. This mosaic of conservation and local disruption is characteristic of plant mitochondrial genomes, which evolve slowly at the sequence level yet can reorganize their structure abruptly. Documenting where those breakpoints occur in Dasiphora provides a baseline for future comparisons with other Rosaceae genera and for studying the mechanisms that drive mitochondrial genome rearrangement in plants.

The significance of this work extends well beyond two shrubs. Dasiphora fruticosa is an ecologically important species across vast tracts of the Northern Hemisphere’s mountainous terrain, and it serves as a model for studying adaptation to cold, high-radiation environments. By delivering the first complete mitogenomic dataset for the genus, the study fills a conspicuous gap in the genomic resources available to plant scientists and lays a foundation for research into how organellar genomes contribute, or fail to contribute, to high-altitude adaptation. The work was supported by China’s National Key Research and Development Program, the Chinese Academy of Sciences’ Light of West China Program, the Second Tibetan Plateau Scientific Expedition and Research program, and the China Postdoctoral Science Foundation, reflecting the strategic importance of alpine ecosystem research in the region. As long-read sequencing becomes routine, studies of this kind are transforming plant mitochondria from genomic backwaters into front-line systems for exploring how genomes move, merge, and evolve inside the cell, and the shrubby cinquefoils of the world’s mountains have now joined that conversation.

Subject of Research: Comparative mitochondrial genomics of the alpine shrubs Dasiphora fruticosa and D. glabra

Article Title: Comparative mitogenomics of Dasiphora fruticosa and D. glabra: insights into genome structure, inter-organellar DNA transfer, and evolutionary dynamics

Article References: Wang, Y., Zhou, H., Jiang, Z., Wang, M., Wang, H., Zhang, H., & Zhou, C. (2026). Comparative mitogenomics of Dasiphora fruticosa and D. glabra: insights into genome structure, inter-organellar DNA transfer, and evolutionary dynamics. BMC Plant Biology. https://doi.org/10.1186/s12870-026-09452-3

Image Credits: AI Generated

DOI: 10.1186/s12870-026-09452-3

Keywords: Dasiphora, mitochondrial genome, mitogenomics, DNA transfer, chloroplast, NUMTs, repeats, phylogeny, nucleotide diversity, alpine plants, Rosaceae, PacBio HiFi

Glucose-6-phosphate isomerase, or GPI, sits near the top of the glycolytic pathway, catalyzing the interconversion of glucose-6-phosphate and fructose-6-phosphate. Every cell in the body depends on glycolysis to some degree, but mature red blood cells are uniquely exposed. Once they extrude their nuclei during development, erythrocytes lose the ability to manufacture replacement enzymes, so their entire energy supply hangs on whatever GPI molecules they started with. When the GPI gene is knocked out on both copies of chromosome 19, the result is hereditary nonspherocytic hemolytic anemia, a condition whose severity ranges from mild, compensated hemolysis to life-threatening hydrops fetalis and, in some families, neurological impairment. More than 40 pathogenic variants have been described, yet clinicians still cannot reliably predict how sick a given patient will become from their genotype alone.

That genotype-phenotype discordance is precisely what motivated the new study. Traditional diagnostic workflows lean heavily on identifying coding-sequence mutations and assessing their likely structural damage to the encoded protein. But occasionally, patients with two defective alleles fare better than those with one, or genetically similar patients diverge sharply in clinical course. Several research groups have proposed that layers of regulation beyond the protein sequence, including allele-specific expression, cis-regulatory elements, and epigenetic modifications, might modulate disease severity in red cell enzymopathies. Direct evidence, however, has been scarce, largely because standard short-read sequencing struggles to resolve these features. Short reads often cannot phase variants across long stretches of DNA, and they cannot simultaneously capture the native methylation marks that decorate the genome.

The Polish-led team behind the new report turned to PacBio HiFi long-read sequencing to break through that barrier. Their patient was a 47-year-old woman with transfusion-dependent hemolytic anemia. Short-read Illumina whole-genome sequencing had already identified two rare heterozygous missense variants in the GPI gene: c.572A>G, which substitutes arginine for histidine at position 191, and c.1414C>T, which substitutes cysteine for arginine at position 472. Both variants are vanishingly rare in the European population, with frequencies below 0.0001 in the gnomAD database. Critically, short reads alone could not determine whether the two mutations sat on the same chromosome or on opposite ones, a distinction that matters enormously for interpreting residual enzyme function.

The long-read data settled the question directly. By phasing the variants across their shared haplotypes, the researchers showed that the mutations sit in trans, one on each parental copy of the gene. Haplotype 1 carries p.His191Arg, while Haplotype 2 carries p.Arg472Cys, and no wild-type allele remains. The team also catalogued 15 non-coding variants within the GPI locus, though deep-learning splice predictions from SpliceAI found no disruption of canonical splicing motifs above the standard reporting threshold. In parallel, full-length isoform sequencing of blood RNA, performed in two independent technical replicates, assembled the GPI transcriptome in unprecedented detail and revealed two previously undescribed minor isoforms alongside the canonical transcript, a descriptive bonus that the authors deliberately kept out of their mechanistic model.

The most provocative observation came from layering transcript counts and methylation calls onto the phased haplotypes. Quantifying allele-resolved reads at the discriminating c.1414 position, the researchers found that Haplotype 2, the one bearing p.Arg472Cys, was modestly over-represented in the transcript pool in both replicates: roughly 62 percent in the first and 55 percent in the second, pooling to about 60 percent overall. The direction was consistent, but the effect was statistically non-significant, and the authors are explicit that this constitutes a trend toward allele-specific expression, not an established one. Isoform-level allelic quantification was underpowered, with too few full-length reads to draw any conclusion at that resolution.

Methylation added a second intriguing thread. Using the polymerase kinetics of PacBio HiFi sequencing, the team obtained native 5-methylcytosine calls across the GPI locus without the bisulfite conversion that conventional methylation mapping requires. The wild-type cytosine at c.1414, preserved on Haplotype 1, sits within a CpG dinucleotide, the genomic context most prone to methylation in human cells. In the small number of informative reads spanning that position, the site was called as methylated on every read carrying the C allele. On Haplotype 2, the C-to-T transition physically destroys the CpG, so the methylated state there is abolished by definition. The co-occurrence of CpG loss on the same allele that shows a directional excess of transcripts is exactly the kind of correlation that fuels a hypothesis, and the authors are careful to label it as precisely that, noting that only a handful of reads support the methylation call.

What makes the model worth testing is what is already known about the protein on Haplotype 2. Prior biochemical work has established p.Arg472Cys as a thermally unstable variant that degrades rapidly, a property consistent with its comparatively low AlphaMissense pathogenicity score of 0.257, compared with 0.914 for p.His191Arg, since structure-based predictors do not necessarily capture instability effects. If the same allele that produces an unstable protein also carries a CpG whose loss slightly lifts transcriptional repression, then the cell might partially compensate by simply making more of a fragile enzyme from that allele. That compensation would depend on a nucleus and would be lost when red cells enucleate and stop transcribing, potentially explaining the age-dependent decline in GPI activity that earlier kinetic studies of hemolytic anemia have proposed. In nucleated tissues, by contrast, continuous transcription from the trans allele might explain why this patient, carrying a p.His191Arg variant previously linked to neurological features, presented with an exclusively hematological phenotype.

The authors are refreshingly transparent about the limits of a single-patient, correlative dataset. They did not measure GPI protein abundance, enzymatic activity, or stability in this patient, and the allelic imbalance never reached statistical significance. To convert their candidate dual-effect variant into an established mechanism, they lay out a demanding experimental agenda: allele-specific transcript quantification at far greater depth with formal correction for mapping bias, statistically powered haplotype-resolved methylation comparison across cell fractions, a direct test of whether the disrupted CpG carries regulatory activity through reporter assays or targeted demethylation, and direct measurement of GPI protein levels and enzymatic activity comparing reticulocyte-enriched and mature erythrocyte fractions. Only such evidence could show that transcriptional output, epigenetic state, and protein stability genuinely interact to shape disease.

Even as a hypothesis, the case marks a turning point in how rare enzymopathies can be examined. One sequencing workflow simultaneously phased compound heterozygous variants, mapped native methylation without chemical conversion, resolved novel transcript isoforms, and generated a specific, testable account of why one patient’s disease unfolded as it did. Short-read exomes and genomes, the workhorses of modern diagnostics, simply cannot deliver that integrated view of a locus. As long-read multi-omics becomes cheaper and more routine, the field will be able to ask systematically whether dual-effect variants, single changes that sabotage both a protein and its own regulatory context, recur across other patients and other genes. If they do, the humble one-letter mutation may turn out to be a far more layered act of molecular sabotage than anyone suspected.

Subject of Research: Integrative long-read multi-omics characterization of a compound heterozygous GPI deficiency patient revealing a candidate dual-effect coding and cis-regulatory variant

Article Title: Integrative Long‐Read Multi‐Omics of a Patient With GPI Deficiency: A Molecular Case Study of a Candidate Dual‐Effect GPI Variant

Article References: Stolarek, I., Delimata‐Raczek, J., Koralewska, N., Sikora, K., Rakoczy, M., Marcinkowska‐Swojak, M., Handschuh, L., Czyż, J., & Figlerowicz, M. (2026). Integrative Long‐Read Multi‐Omics of a Patient With GPI Deficiency: A Molecular Case Study of a Candidate Dual‐Effect GPI Variant. Journal of Cellular and Molecular Medicine, 30(17), Article e71338. https://doi.org/10.1111/jcmm.71338

Image Credits: AI Generated

DOI: 10.1111/jcmm.71338

Keywords: GPI deficiency, long-read sequencing, PacBio HiFi, hemolytic anemia, allele-specific expression, DNA methylation, glycolysis, red blood cells, compound heterozygosity, multi-omics, Iso-Seq, rare disease

The technical achievement at the heart of the study is considerable. Building a truly complete genome has long been one of genomics’ most stubborn challenges, because repetitive DNA sequences near the centers and ends of chromosomes resist standard sequencing methods. The research team overcame this by combining PacBio HiFi long-read sequencing, which produces highly accurate reads spanning tens of thousands of DNA letters, with Hi-C technology, which captures the physical contacts between distant stretches of chromosomes and allows the sequences to be anchored into their correct chromosomal positions. The final assembly spans 2.42 gigabases, with a scaffold N50 of 72.35 megabases and a contig N50 of 66.56 megabases, meaning that half of the assembled genome resides in stretches longer than those figures. The guanine-cytosine content of the assembly sits at 42.36 percent. In practical terms, the contiguity of this assembly approaches the ideal of a single gap-free sequence per chromosome, giving researchers an unprecedented view of the bear’s genome architecture.

With that high-resolution reference in hand, the team turned to comparative genomics, aligning the Tibetan brown bear’s chromosomes against those of its closest living relative, the polar bear. The analysis revealed high chromosomal synteny between the two species, meaning that large blocks of genes retain the same order and orientation across both genomes despite millions of years of separate evolution. This structural conservation is scientifically valuable for two reasons. It confirms that the assembly is accurate, since genuine chromosome-scale structure should be preserved between closely related ursids, and it provides a stable framework for pinpointing the regions that have instead diverged, which are precisely the locations where adaptation to very different environments, Arctic sea ice versus high-altitude steppe, is most likely to leave its mark.

The genomic record also preserves a vivid portrait of the species’ demographic turbulence. By reconstructing changes in effective population size through time, the researchers found that the Tibetan brown bear expanded during the early Pleistocene, a period of dynamic climate and habitat change. That expansion, however, was followed by two severe population bottlenecks. The first occurred approximately two million years ago, coinciding with the pre-Poyang Glaciation, and the second struck around one hundred thousand years ago during the Last Glacial Period. Together, these crashes reduced the bear’s cumulative effective population size by nearly eighty percent. For a species that today persists only on the Qinghai-Xizang Plateau, this history of repeated near-collapses carries a sobering message about the genetic fragility that may still lurk beneath its wild, unbroken landscape.

Yet survival through those bottlenecks suggests the species carries more than vulnerability. Enduring glacial cycles on a plateau whose average elevation exceeds 4,000 meters demands a physiology unlike that of lowland relatives, and the most striking discovery of the study concerns exactly that. When the researchers examined candidate hibernation-related pathways, they found that the Tibetan brown bear displays a distinctly carbohydrate-dominated metabolic signature. This stands in sharp contrast to the hibernation strategy documented in North American brown bears, which rely primarily on lipid, or fat-based, metabolism to fuel their long winter dormancy. In other words, two populations of the same widespread species appear to have evolved fundamentally different biochemical approaches to the same survival problem.

The implications of that finding extend well beyond bears. Fat is the standard fuel of hibernation across most studied mammals, because it stores more energy per gram and spares blood sugar during months of fasting. A sugar-first strategy on the Tibetan Plateau hints at selective pressures unique to that environment, where the short growing season, the composition of available foods, and the metabolic demands of life at low oxygen levels may have favored a different balance of carbohydrate and lipid pathways. The authors emphasize that this distinct metabolic adaptation highlights the bear’s unique mechanisms for surviving the extreme plateau environment, and it gives physiologists a natural experiment in alternative hibernation biochemistry encoded within a single species.

For conservation biologists, the new genome arrives at a critical moment. The Tibetan brown bear sits at the top of the plateau’s food web, playing a crucial ecological role, yet high-quality genomic resources for the species have been scarce until now. A chromosome-level reference genome transforms what conservation managers can do: it enables accurate estimates of genetic diversity and inbreeding, allows the tracking of gene flow between populations, and provides the resolution needed to identify locally adapted lineages that deserve special protection. Given the roughly eighty percent historical reduction in effective population size documented in the study, such tools are not a luxury. They are the baseline data upon which any serious plan for the species’ long-term management on the plateau must rest.

The study also strengthens the broader evolutionary picture of the bear family. Ursids have long fascinated evolutionary biologists because the family contains species with radically different ecologies, from the omnivorous brown bear to the hypercarnivorous, ice-bound polar bear. A near telomere-to-telomere assembly for the Tibetan brown bear adds a critical high-quality data point for reconstructing the phylogenomic relationships within Ursus and for dating the divergences that produced today’s brown bear lineages. Because structural variants, gene duplications, and regulatory regions can now be examined in their full chromosomal context rather than through fragmented assemblies, questions about how bears colonized and adapted to some of Earth’s harshest habitats can be addressed with far greater precision than before.

The research was carried out by Muran Zhao, Anmin Wang, Hai Liu, Chenxing Yu, Yanlin Liu, Nan Sun, and Guogang Li of Qinghai Normal University, with fieldwork, sample collection, and laboratory analysis conducted under permits granted by the university. The work was supported by the National Natural Science Foundation of China. The team notes that the new assembly is intended as a resource for the wider community, facilitating future studies of ursid evolutionary history and supporting conservation and management of the plateau’s wildlife. Because the genome is derived from a female bear, it also enables improved analysis of the sex chromosomes, an area where earlier fragmented assemblies often fell short.

What began as a technical sequencing project has ended with a discovery that may echo well beyond the Qinghai-Xizang Plateau. If a close relative of the brown bears studied across North America and Eurasia can hibernate on a carbohydrate-dominated metabolic program, then hibernation is not a single fixed solution that evolution produced once and reused everywhere. It is a flexible toolkit, and different populations have assembled it from different parts. As climate change alters the length and severity of winters across the world’s mountains, understanding that flexibility, and the genes that underpin it, could prove essential not only for protecting the Tibetan brown bear but for predicting which hibernating species can adjust and which cannot. The near-complete genome of this high-altitude survivor is now available as the roadmap for answering those questions.

Subject of Research: Chromosome-level genome assembly and hibernation-related evolution of the Tibetan brown bear (Ursus arctos pruinosus)

Article Title: Near telomere-to-telomere genome reveals the phylogenomics and hibernation-related evolution in Tibetan brown bear (Ursus arctos pruinosus)

Article References: Zhao, M., Wang, A., Liu, H., Yu, C., Liu, Y., Sun, N., & Li, G. (2026). Near telomere-to-telomere genome reveals the phylogenomics and hibernation-related evolution in Tibetan brown bear (Ursus arctos pruinosus). BMC Genomics. https://doi.org/10.1186/s12864-026-13324-3

Image Credits: AI Generated

DOI: 10.1186/s12864-026-13324-3

Keywords: Tibetan brown bear, telomere-to-telomere genome, genomics, hibernation, Qinghai-Xizang Plateau, population bottlenecks, PacBio HiFi, Hi-C sequencing, phylogenomics, conservation genomics, carbohydrate metabolism, Ursidae evolution