The team started with a wild-type Citrobacter strain originally isolated from rice rhizosphere soil, the thin zone of earth immediately surrounding living roots where chemistry, biology and plant physiology collide. From this native microorganism they constructed engineered variants by boosting the expression of the bacterium’s own functional genes rather than introducing foreign DNA. Two over-expression lines, designated Citrobacter sp. XT1-2-2-::APS and Citrobacter sp. XT1-2-2-::SiR, carry enhanced copies of key reductase genes, enzymes that drive reduction reactions capable of altering the chemical state of heavy metals. The choice of genes matters: reductive transformations are central to how certain bacteria immobilize cadmium, converting more mobile or bioavailable forms of the metal into species that plant roots struggle to absorb.

To test the concept, the researchers ran pot experiments comparing rice grown with the wild-type strain against rice inoculated with the two over-expression strains. At the maturity stage of the crop, the decisive measurement was the concentration of available cadmium in the soil, the fraction of the metal that plants can actually take up. Both engineered strains produced a significant reduction in available cadmium relative to the wild-type inoculant, indicating that amplifying the reductase genes translated directly into stronger immobilization of the metal in the rhizosphere. Because the engineered strains were built by enhancing endogenous genes, the biosafety profile of the construct is comparatively straightforward; the authors note that hemolysis assays were negative and that all work with the engineered microorganisms followed institutional containment guidelines, with animal toxicity studies and field trials planned before any environmental release.

The cadmium result alone would be noteworthy, but the study’s real ambition was to understand the mechanism, and that required looking at the soil as an ecosystem rather than a chemical reservoir. Using high-throughput sequencing of the 16S rRNA gene for bacteria and the internal transcribed spacer region for fungi, the team mapped the structure and diversity of the rhizosphere microbial communities under each treatment. The analysis, processed through standard pipelines such as QIIME with operational taxonomic units as the unit of counting, revealed that inoculation with the over-expressing strains significantly altered the community structure of both bacteria and fungi compared with the wild-type treatment, and that community abundance and diversity increased in the engineered-strain soils.

Those shifts were not random. Several bacterial phyla rose in relative abundance in soils treated with the over-expression strains, including Bacillota, Bacteroidota, Pseudomonadota, Thermodesulfobacteriota and Acidobacteriota. Each of these groups has known associations with soil health and metal chemistry: Bacillota include many spore-forming, plant-growth-promoting species; Pseudomonadota encompass numerous rhizosphere colonizers and metal-tolerant taxa; and sulfate- and sulfur-metabolizing lineages within the Thermodesulfobacteriota can influence metal speciation through the production of metal sulfides. On the fungal side, the community was dominated by Ascomycota, Chytridiomycota, Basidiomycota and unclassified fungi, with the engineered inoculants again shifting the balance. The picture that emerges is of a rhizosphere whose resident ecology is nudged into a configuration less favorable to cadmium mobility.

The second half of the investigation moved inside the plant. The researchers performed a proteomic analysis of rice roots using tandem mass tag labeling combined with liquid chromatography-tandem mass spectrometry, a technique that allows many protein samples to be compared quantitatively in a single experiment. TMT labels are attached to peptides from each sample, the samples are pooled and analyzed together, and the relative intensity of each tag reports how abundant the corresponding protein was in the original tissue. Applied to rice roots, this approach identified differentially expressed proteins between treatments, giving a systems-level view of how the plant responds to the engineered bacteria at the molecular level.

The differentially expressed proteins were enriched mainly in biosynthetic pathways, metabolic pathways, glutathione metabolism, cell wall organization and antioxidant defense pathways, as annotated against the Kyoto Encyclopedia of Genes and Genomes. Each of these categories tells a coherent story about cadmium stress. Glutathione is a tripeptide that binds heavy metals and feeds the synthesis of phytochelatins, peptides that sequester cadmium inside plant cells and limit its toxicity. Cell wall organization matters because the wall is the first barrier a metal ion encounters on its way into a root cell, and modifications to its composition can reduce metal entry. Antioxidant defenses counter the oxidative damage that cadmium provokes indirectly by disrupting electron transport and depleting cellular thiols. Together, the proteomic signature suggests that rice roots inoculated with the engineered strains mount a coordinated defensive and metabolic reprogramming, consistent with the observed reduction in cadmium uptake.

What makes the study conceptually interesting is the dual mechanism it documents. The engineered bacteria act externally, lowering the pool of available cadmium in the rhizosphere and recruiting or enriching a microbial community that supports that state, while acting internally, priming the plant’s own stress-response and sequestration machinery through changes in root protein expression. The authors frame this as illuminating the complex interactions between rice and its associated microorganisms, and specifically as showing how over-expression of key reductase genes in XT1-2-2 propagates outward to the endophytic and rhizosphere microbiome and inward to the host proteome. In other words, a targeted genetic edit in one soil bacterium produces cascading effects across a whole ecological network, and those cascades converge on the practical outcome of less cadmium in the crop.

The practical implications extend well beyond one pot experiment. Rice is the dietary staple for billions of people and is unusually efficient at accumulating cadmium from flooded soils, which is why cadmium-tainted rice remains a recurring public health scandal in several rice-growing regions. Conventional remediation, such as soil replacement, chemical immobilizers or growing non-food crops on contaminated land, is expensive and slow. A living inoculant that can be applied to paddies, derived from a native soil organism and engineered only to over-produce its own enzymes, represents a potentially cheaper and more scalable intervention. The authors position the work as theoretical support for applying functional microorganisms to promote rice growth while reducing cadmium accumulation, and the combined microbiome and proteome datasets provide the mechanistic grounding that regulators and agronomists would need before such an approach could move outdoors.

Caution is still warranted. The experiments were conducted under controlled laboratory conditions, and the authors themselves state that further animal toxicity studies and field trials are planned to evaluate the biosafety of the engineered strain prior to any potential environmental application. Rhizosphere communities vary enormously with soil type, water management, climate and cropping history, so the community shifts observed in pots will need to be reproduced in real paddies. Yet the study stands as a striking demonstration of synthetic biology meeting soil ecology: by tuning the expression of two reductase genes in a bacterium that rice roots already know, the researchers changed the chemistry of the soil, the composition of the microbial ecosystem and the protein inventory of the plant, and in doing so kept a toxic heavy metal out of the food chain. If the field results hold, the humble rhizosphere may become the most important battleground in the fight against cadmium in rice.

Subject of Research: Engineered Citrobacter bacteria reducing cadmium uptake in rice through rhizosphere microbiome modulation and root proteomic changes

Article Title: A genetically engineered Citrobacter sp. XT1-2-2 modulates rhizosphere soil microbiome and induces root proteomic changes to reduce cadmium uptake in rice

Article References: Liu, Z., Du, D., Cheng, W., Yan, Z., Wei, X., Zhang, M., Li, Y., Wu, S., & Shan, S. (2026). A genetically engineered Citrobacter sp. XT1-2-2 modulates rhizosphere soil microbiome and induces root proteomic changes to reduce cadmium uptake in rice. BMC Genomics. https://doi.org/10.1186/s12864-026-13465-5

Image Credits: AI Generated

DOI: 10.1186/s12864-026-13465-5

Keywords: Citrobacter sp. XT1-2-2, cadmium, rice, rhizosphere microbiome, soil microbiology, genetic engineering, proteomics, TMT labeling, glutathione metabolism, heavy metal immobilization, plant growth-promoting bacteria, BMC Genomics

The significance of the work lies in what it resolves about a plant that occupies a peculiar dual identity. Gardenia is simultaneously a beloved ornamental, prized for its glossy evergreen foliage and perfume-grade blooms, and a cornerstone of traditional Chinese medicine, where its dried fruits are used in formulations believed to clear heat and dampness. That combination of horticultural and pharmacological value has driven the species through millennia of domestication and geographic expansion, yet the genetic architecture underlying its population differentiation and trait variation had remained largely unresolved. Earlier studies relied on low-resolution molecular markers such as amplified fragment length polymorphisms and simple sequence repeats, which could sketch broad relationships but could not pinpoint the genes or pathways shaped by selection. Whole-genome resequencing changes that calculus entirely, allowing researchers to examine variation at single-base resolution across the entire genome rather than at a few hundred scattered landmarks.

To build the dataset, the team sampled accessions spanning nine populations that together cover the species’ cultivated and wild range in China. Genomic DNA was extracted and subjected to high-throughput sequencing, and reads were mapped to a reference genome to call variants. Quality control metrics reported in the study, including the transition-to-transversion ratio, a standard benchmark of variant-calling accuracy in population genomics, indicated that the SNP set was suitable for downstream evolutionary analyses. The resulting matrix of 7.19 million SNPs is an order of magnitude richer than anything previously available for the species, and the researchers have released it as a foundational resource for the community, with implications for marker-assisted breeding and germplasm conservation.

One of the study’s most consequential findings concerns the relationship between wild and cultivated gene pools. Despite millennia of domestication, the analyses showed limited divergence between wild-type and cultivated-type accessions, suggesting that gardenia’s journey alongside humans has involved gene flow and shared ancestry rather than a sharp genetic split. Domestication in many crops produces a pronounced bottleneck, with cultivated varieties carrying only a fraction of the diversity found in their wild relatives. Gardenia tells a softer story, one consistent with a plant whose cultivation history involved continuous exchange of planting material across regions rather than isolation on a single farm or in a single monastery garden. That continuity is good news for conservation, because it means the wild gene pool and the cultivated one remain genetically intertwined rather than drifting apart into separate evolutionary trajectories.

Where the species does show clear structure is along geographic lines. Population genomic analyses, including principal component analysis and neighbor-joining clustering, revealed that accessions group according to where they were collected, forming distinct regional clusters that mirror China’s varied landscapes and climates. Among these populations, one stands out with a profile unlike any other: Jiangxi, the province in southeastern China where the study’s lead institutions are based, displays a distinctive genetic signature characterized by the lowest nucleotide diversity, the shallowest linkage disequilibrium decay, and the highest genetic cohesion of any population examined. Each of those statistics tells a story. Low nucleotide diversity means the Jiangxi accessions carry fewer alternative variants at any given site, a hallmark of a population that has passed through a bottleneck or been founded by a small number of ancestors. Shallow linkage disequilibrium decay, by contrast, is a subtler signal: linkage disequilibrium describes the non-random association of variants along chromosomes, and its decay with distance reflects how quickly recombination scrambles those associations. A shallow decay pattern in Jiangxi suggests larger effective population sizes or stronger historical connectivity, allowing recombination to break down associations efficiently over generations.

The combination of low diversity and high cohesion in Jiangxi is intriguing precisely because it is unusual. A population can lose diversity through a founder event, but cohesion, measured as the degree to which individuals within a population are genetically similar to one another relative to outsiders, points to a lineage that has remained reproductively isolated from its neighbors while expanding from a limited ancestral pool. The researchers suggest that Jiangxi may represent a distinct evolutionary unit within the species, one whose history of isolation and local adaptation has carved out a genetic identity recognizable even in a dataset of nearly 200 individuals. For breeders, that distinctiveness matters: populations with unique allele combinations can carry traits absent elsewhere in the species, making Jiangxi germplasm a potentially valuable reservoir for future cultivar development, provided it is conserved before habitat loss erodes it.

Beyond describing structure, the team went hunting for the genomic regions shaped by divergent selection, the process by which different populations experience contrasting selective pressures that push their genomes in different directions. A genome-wide scan identified 44 such regions harboring 597 candidate genes. When the researchers subjected those candidates to pathway enrichment analysis using the Kyoto Encyclopedia of Genes and Genomes, three pathways rose to statistical significance: brassinosteroid biosynthesis, ascorbate and aldarate metabolism, and diterpenoid biosynthesis. The enrichment is not a random assortment of biochemistry. Brassinosteroids are plant steroid hormones that regulate growth, cell elongation, stress tolerance and reproductive development, making them plausible modulators of key agronomic traits such as plant architecture and flowering. Ascorbate, better known as vitamin C, sits at the center of antioxidant defenses, and its metabolism pathway connects directly to how plants cope with oxidative stress, a trait of obvious importance for a shrub cultivated across China’s climatic diversity.

The third enriched pathway, diterpenoid biosynthesis, may be the most pharmacologically tantalizing of the three. Diterpenoids are a large class of plant secondary metabolites built from four isoprene units, and in gardenia they connect directly to the species’ medicinal value. Gardenia fruits are rich in iridoid glycosides and crocins, pigmented and bioactive compounds whose biosynthesis flows through terpenoid pathways, and the identification of diterpenoid biosynthesis genes among regions under divergent selection suggests that different gardenia populations may have been shaped, in part, by selection on the very chemistry that makes the plant medicinally valuable. The researchers are careful to frame these as candidate genes whose roles plausibly modulate key agronomic and pharmacological traits, rather than proven targets, and follow-up functional validation will be needed to confirm the connections. But the enrichment analysis provides a concrete starting point for dissecting the molecular physiology of population differentiation in a species where chemistry is the product.

The study’s practical implications extend in two directions. For breeding, the dense SNP map enables marker-assisted selection, the technique of using DNA markers to track desirable traits without waiting years for plants to mature and express them. Breeders working on gardenia can now associate markers with traits such as flower quality, fruit chemistry and stress tolerance, accelerating the development of improved cultivars for both ornamental and medicinal markets. For conservation, the clear geographic structure revealed by the analyses provides a roadmap for prioritizing which populations to protect. Because regional lineages carry unique allele combinations shaped by local evolutionary history, conserving a single large population is not sufficient; the genetic resources of the species are distributed across its range, and each population, from the distinctive Jiangxi lineage to the others sampled, represents an irreplaceable slice of the species’ evolutionary legacy.

What makes the gardenia study resonate beyond its species is the methodological template it offers for other underexplored medicinal and ornamental plants. Whole-genome resequencing at this scale, once the preserve of major crop species with billion-dollar industries, is now accessible for specialty crops and wild relatives, and the analytical pipeline deployed here, from variant calling through population structure to selection scans and pathway enrichment, can be applied wherever a reference genome exists. As demand grows for plant-derived pharmaceuticals, natural pigments and resilient ornamentals, the genomic resources assembled for gardenia position the species as a model for understanding how millennia of cultivation and geographic expansion shape the genetic architecture of a plant that humans have valued, in different ways, for thousands of years. The white flowers will keep blooming, but behind them now stands a molecular map of everything that made them possible.

Subject of Research: Genomics

Article Title: Whole-genome resequencing reveals genetic diversity and population differentiation in Gardenia jasminoides J. Ellis

Article References: Yang, W., Li, K., Yu, L., Duan, Y., Zhu, S., & Deng, S. (2026). Whole-genome resequencing reveals genetic diversity and population differentiation in Gardenia jasminoides J. Ellis. BMC Genomics. https://doi.org/10.1186/s12864-026-13456-6

Image Credits: AI Generated

DOI: 10.1186/s12864-026-13456-6

Keywords: Gardenia jasminoides, whole-genome resequencing, genetic diversity, population differentiation, SNPs, InDels, brassinosteroid biosynthesis, diterpenoid biosynthesis, marker-assisted breeding, germplasm conservation, BMC Genomics, Whole-genome

Genomic islands are discrete regions of a bacterial chromosome that often originate from other microbes. They typically differ in nucleotide composition from the rest of the genome, are flanked by repeats or insertion sequences, and can encode anything from secretion systems to metabolic enzymes. Because they move between lineages, they act as hotspots of genomic plasticity, allowing bacteria to acquire new functions in single evolutionary leaps rather than through slow accumulation of point mutations. For S. thermophilus, whose transition from a commensal ancestor living in mammalian guts to an industrial workhorse in fermenters involved extensive genomic remodeling, these islands are a natural place to look for the traits that set high-performing strains apart from the rest.

The research team assembled a dataset of 109 S. thermophilus genomes and subjected them to a battery of comparative and predictive bioinformatics analyses. Phylogenomic trees built from single-copy orthologs and pangenome reconstructions revealed that the species divides into two major lineages, which the authors labeled STH1 and STH2. This clustering mattered, because the two branches turned out to differ in their evolutionary dynamics. Strains belonging to STH2 showed markedly greater genomic plasticity, with twenty equivalent genomic islands shared among members of this lineage, suggesting that this branch has been more receptive to horizontal gene acquisition or has exchanged DNA more frequently among its members.

To capture the full functional landscape of the species, the researchers analyzed the pangenome, the sum of all genes found across all strains. They found that genes associated with probiotic-related functions, including cell wall and membrane biosynthesis, amino acid and carbohydrate metabolism, and secondary metabolite production, were concentrated predominantly in the accessory and unique fractions of the genome rather than the core. In other words, the traits most relevant to probiotic potential are not universal properties of the species but strain-specific endowments, and genomic islands appear to be a major vehicle delivering them. Among island categories, metabolic islands were the most abundant and accounted for the majority of probiotic-associated functions, followed by virulence islands, whose limited cargo is itself informative about the species’ biosafety.

That biosafety picture is one of the study’s most reassuring findings. Although S. thermophilus is a close relative of streptococci that cause disease, the team found that both the panvirulome, the complete set of virulence-related genes across all 109 strains, and the panresistome, the corresponding repertoire of antimicrobial resistance genes, were small. Features such as capsule-associated genes, rather than classical virulence factors, may contribute to probiotic performance by helping the bacterium persist and interact with host tissues. The limited virulence and resistance gene repertoire, the authors conclude, further supports the biosafety of the species and reinforces its suitability for industrial and probiotic applications, a conclusion consistent with its long GRAS and QPS regulatory standing.

The researchers went beyond gene prediction to test actual metabolic behavior. Using Biolog phenotyping, which measures a bacterium’s ability to metabolize a panel of substrates under different oxygen conditions, they examined three candidate strains: CIRM-BIA 16, CIRM-BIA 336 and CIRM-BIA 961. Overall, metabolic activity was limited to a few substrates under both aerobic and anaerobic conditions, reflecting the well-known nutritional specialization of S. thermophilus to the lactose-rich environment of milk. However, one strain stood out. The metabolic island identified in CIRM-BIA 336 carries genes that may support enhanced galactose utilization under anaerobic conditions through the tagatose-6-phosphate pathway, a route for fermenting this milk sugar that many dairy strains lack.

This detail has practical resonance for the dairy industry. Galactose accumulates when S. thermophilus ferments lactose in milk, and its incomplete metabolism can cause problems in products such as yogurt and reduced-lactose dairy foods, where residual galactose affects texture, sweetness and shelf life. A strain equipped with a genomic island enabling tagatose-6-phosphate pathway-dependent galactose consumption could therefore offer both a technological advantage in fermentation and a route to gentler products for galactose-sensitive consumers. The finding illustrates how a single mobile genetic element can translate directly into a phenotype of commercial interest, and it supports the authors’ broader conclusion that genomic islands are associated with strain-specific probiotic potential and metabolic versatility.

The mobilome analysis extended beyond islands proper. The team also mapped integrative conjugative elements and integrative mobilizable elements, two classes of mobile DNA that can carry entire functional modules between cells, sometimes across species boundaries. A striking pattern emerged: most genomic islands, as well as these integrative elements, occupied conserved positions in the chromosome across strains. This suggests that integration is not random but prefers specific hotspots, likely reflecting insertion-site preferences of the phage and plasmid machinery that mediates transfer, as well as selection against disrupting essential genes. Conserved landing pads for mobile DNA mean that parallel acquisitions can occur in different lineages, generating equivalent islands with related cargo that shape the species’ evolution in reproducible ways.

The study also probed the potential for bacteriocin production and secondary metabolite biosynthesis, traits prized in probiotics because bacteriocins can inhibit competing pathogens. Using antiSMASH to detect biosynthetic gene clusters and the DRAMP database to identify putative bacteriocins, the researchers documented these features across the dataset and linked them to the accessory genome and its mobile elements. Combined with the capsule-associated genes and the enriched cargo of metabolic islands, the picture is of a species whose probiotic-relevant toolkit is distributed unevenly among strains and is built, to a substantial degree, from horizontally acquired parts. This aligns with the emerging discipline of probiogenomics, which the authors list among the study’s keywords alongside horizontal gene transfer, lactic acid bacteria, comparative genomics and bioinformatics.

For strain designers and fermentation technologists, the practical message is that selecting S. thermophilus strains for next-generation fermented foods or probiotic formulations should involve screening their genomic islands, not just their core genomes. Two strains with nearly identical housekeeping genes may differ sharply in carbohydrate metabolism, defense molecules and surface features because of what their islands carry. The work also offers a caution: genomic islands act as major hotspots of genomic plasticity, particularly within the STH2 lineage, so genetic stability over industrial propagation and storage deserves attention when a desirable trait is island-encoded. With this study, the humble yogurt bacterium reveals itself as a genomic patchwork, stitched together over millennia of horizontal gene transfer, and the seams of that patchwork turn out to hold precisely the traits that make it valuable in our food and potentially in our guts.

Subject of Research: Genomic islands underlying host adaptation and probiotic potential in Streptococcus thermophilus

Article Title: Genomic islands associated with host adaptation and probiotic-related functions in Streptococcus thermophilus

Article References: Ariute, J. C., Roland, N., Azevedo, V., Le Loir, Y., Aburjaile, F. F., & Guédon, E. (2026). Genomic islands associated with host adaptation and probiotic-related functions in Streptococcus thermophilus. BMC Genomics. https://doi.org/10.1186/s12864-026-13419-x

Image Credits: AI Generated

DOI: 10.1186/s12864-026-13419-x

Keywords: Streptococcus thermophilus, genomic islands, probiotics, horizontal gene transfer, pangenomics, lactic acid bacteria, panvirulome, panresistome, galactose metabolism, dairy fermentation, BMC Genomics, biosafety

The experiment was deliberately simple in design but technically ambitious in execution. Whole peripheral blood from three healthy donors was exposed to 4 Gy of X-rays, a dose comparable to what might be encountered in a serious radiation accident or in the upper range of clinical radiotherapy fractions, while sham-irradiated samples from the same donors served as controls. Six hours after exposure, the researchers isolated the leukocytes and profiled their RNA using Oxford Nanopore Technologies’ long-read platform. Unlike short-read sequencers that chop RNA molecules into fragments of a hundred or so bases and then computationally stitch them back together, nanopore sequencers read entire RNA transcripts from end to end, preserving the full-length structure of each isoform.

At the gene level, the analysis confirmed what decades of radiation biology have established. Compared with the sham-irradiated controls, 117 genes were significantly upregulated and 66 were downregulated, and the enriched pathways were the canonical ones: DNA damage repair and inflammatory responses. Genes involved in detecting double-strand breaks, signaling their presence, and orchestrating their repair were switched on, while the cell’s routine metabolic housekeeping was dialed back. This validation matters because it demonstrates that the long-read approach recovers the same radiation-response signatures that biodosimetry researchers have relied on for years, meaning the new method is not a departure from established science but an extension of it.

The real surprises emerged when the team descended from genes to transcripts. A single gene can produce multiple RNA isoforms through alternative splicing, and these isoforms can be regulated independently of one another. The long-read data revealed 102 transcripts that were significantly upregulated and 17 that were downregulated, a pattern of isoform-specific regulation that was completely masked at the gene level. The most striking example was IL32, an inflammatory cytokine gene that showed no significant change when measured as a whole gene, yet two of its transcript isoforms were strongly upregulated, with log2 fold-changes of 1.19 and 1.99. In other words, radiation was selectively boosting particular molecular versions of IL32 while the gene-level average diluted the signal into statistical invisibility.

IL32 is not a random bystander. It is deeply involved in inflammatory and immune signaling, and its isoform-specific induction after irradiation suggests that the immune reaction to radiation damage is more finely tuned than gene-level studies imply. The researchers also identified changes in transcript usage for WDR74, ITM2B, AK2, and RPS19, meaning that irradiation altered which isoform of these genes the cells preferentially produced rather than how much total gene product was made. WDR74 participates in DNA damage response signaling, AK2 is a kinase with roles in apoptosis and energy metabolism, and RPS19 is a ribosomal protein essential for red blood cell production, so these usage shifts point to radiation rewiring not just expression levels but the functional architecture of the transcriptome itself.

Beyond confirming known signatures and revealing hidden regulation, the long-read approach enabled genuine discovery. The team identified novel candidate isoforms that had never been annotated in reference databases, including protein-coding isoforms of PHPT1 and FBXO22 as well as several non-coding candidates. PHPT1 encodes a phosphohistidine phosphatase, and FBXO22 is part of the ubiquitin ligase machinery that tags proteins for degradation, a process with known links to DNA damage response. The newly found FBXO22 isoform showed evidence of altered protein domain architecture relative to the reference isoform, meaning that if this isoform is translated, it could produce a protein with a different functional structure than the canonical version. To assess coding potential, the researchers ran the candidates through tools including TransDecoder, the Coding Potential Calculator 2, and CPAT, and used SQANTI3 for structural quality annotation of the novel isoforms, a standard workflow for separating real transcripts from sequencing artifacts.

The clinical and biodosimetric implications of this work are substantial. Biodosimetry, the estimation of absorbed radiation dose from biological markers, is critical after radiological accidents when physical dosimeters are absent, and gene expression profiling of blood is one of the most promising approaches because the transcriptomic response scales with dose. The Mainz study shows that transcript-level markers can carry information that gene-level markers miss, potentially improving the sensitivity and specificity of dose estimates. The IL32 finding is a case in point: a biodosimetric panel built only on gene-level data would overlook a radiation-responsive signal that isoform-level analysis captures clearly.

Portability adds another dimension. Oxford Nanopore sequencers are compact enough to be carried to the field, which opens the possibility of on-site gene expression studies after radiation emergencies, in disaster zones, or in military settings where samples cannot be shipped to a central laboratory. The authors note that the diffusion of long-read RNA sequencing in radiation biology has been limited despite these advantages, and their study provides a systematic demonstration of what the technology can extract from irradiated blood. The workflow, from whole blood exposure through leukocyte isolation to full-length transcript profiling, could serve as a template for larger biodosimetric studies.

Radiotherapy is another arena where the findings could prove useful. Evaluating the transcriptomic radiation response of blood provides insight into absorbed dose, hematotoxicities, and immune reactions, and transcript-level resolution could help clinicians understand why some patients develop severe immune reactions or blood toxicity while others tolerate the same dose. Isoform-specific markers might eventually distinguish inflammatory signatures that predict toxicity from those that reflect routine tissue response, refining the personalization of radiation treatment.

The study has clear limitations that the authors acknowledge through their design choices. With only three donors, the analysis captures individual variability but cannot establish population-level biomarkers, and the single 4 Gy dose at a single 6-hour time point leaves the dose-response and time-course dimensions unexplored. The researchers addressed donor effects with leave-one-out sensitivity analyses and compared their long-read results against short-read pipelines, finding that 34 differentially expressed genes were shared between approaches while others were platform-specific. Still, the conceptual advance stands: radiation-induced transcriptomic alterations include isoform-level regulatory events and novel candidate isoforms that are masked in gene-level analyses and conventional short-read pipelines. As long-read sequencing matures, the molecular record of radiation exposure may prove far richer, and far more clinically actionable, than the gene-level view has allowed.

Subject of Research: Radiation-induced transcriptomic signatures in human leukocytes revealed by long-read RNA sequencing

Article Title: Decoding radiation-induced transcriptomic signatures of human leukocytes using long-read RNA-seq: clinical and biodosimetric implications

Article References: Salah, A., Schmidberger, H., Marini, F., & Zahnreich, S. (2026). Decoding radiation-induced transcriptomic signatures of human leukocytes using long-read RNA-seq: clinical and biodosimetric implications. BMC Genomics. https://doi.org/10.1186/s12864-026-13475-3

Image Credits: AI Generated

DOI: 10.1186/s12864-026-13475-3

Keywords: long-read RNA sequencing, biodosimetry, ionizing radiation, transcriptomics, Oxford Nanopore, human leukocytes, isoform discovery, BMC Genomics, Decoding, radiation-induced, transcriptomic, signatures

Propolis is a composite substance that worker bees manufacture from plant resins they collect in the field and blend with their own secretions. It seals cracks, reinforces the hive structure, and modulates the microbial environment inside the colony, making it one of the most functionally important products of bee labor. Yet propolis yield varies enormously between colonies, and the biological basis of that variation has remained poorly understood. Because the antennae are the primary sensory organs that detect signals from nectar, pollen, and propolis sources during foraging, the researchers reasoned that resin collection behavior might be tightly linked to olfaction at the colony level, a relationship that had never been directly tested between colony types.

To investigate, the team compared foragers from high propolis-producing colonies with foragers from ordinary colonies of A. mellifera. The first striking result came from morphology: the total length of both the left and right antennae was significantly shorter in high-yield colonies than in ordinary colonies. When the researchers broke the antenna down segment by segment, they found that the lengths of two flagellomeres, FL7 and FL10, were also significantly shorter in the high-yield bees. The flagellum carries the densest array of olfactory sensilla on the bee antenna, so differences in the dimensions of specific flagellomeres suggest that the sensory hardware itself differs between bees that collect resin aggressively and those that do not.

The morphological differences were only the beginning. The researchers sequenced the antennal transcriptomes of the two colony types and analyzed FPKM heatmaps to visualize gene expression patterns across samples. They then comprehensively characterized three gene families central to insect olfaction: the odorant binding proteins, or OBPs; the chemosensory proteins, or CSPs; and the odorant receptors, or ORs. In total, the analysis covered six CSPs, 19 OBPs, and 34 ORs, amounting to 58 genes whose expression could be compared between the high-yield and ordinary colonies.

Characterizing these gene families required more than simple counting. The team mapped the chromosomal locations of each gene and examined their conserved protein domains to confirm family membership, then performed phylogenetic analyses that placed the honeybee genes in evolutionary context with four other species: the eastern honeybee Apis cerana, the parasitoid wasp Chouioia cunea, the jewel wasp Nasonia vitripennis, and bumblebees of the genus Bombus. These comparative trees allow researchers to distinguish genes that have been conserved across social insects from those that have expanded or contracted specifically in the western honeybee lineage, providing clues about which olfactory components might be specialized for resin detection and other colony-level tasks.

The expression results revealed clear differences between the two colony types. Genes from all three families showed differential expression between high propolis-producing and ordinary colonies, indicating that the olfactory gene repertoire is deployed differently in bees whose colonies yield drastically different amounts of resin. This is the first time that differences in antennal morphology and antennal RNA-seq results have been reported between colonies with differential propolis yield, as well as the first report of differences in the expression of these three gene families between the two colony types. The finding implies that the propensity to collect propolis is not merely a matter of opportunity or environment but is associated with measurable molecular differences in the sensory organs of the foragers themselves.

To make sure the sequencing results were robust, the researchers validated the expression patterns of all 58 genes by quantitative real-time PCR, or qRT-PCR. The qRT-PCR expression levels were consistent with the RNA-seq expression patterns across the board, confirming that the transcriptomic differences were not artifacts of library preparation or normalization. This dual-method validation is technically important because antennal transcriptomes can be challenging: the tissue is small, expression levels of individual chemosensory genes are often low, and subtle differences between biological groups can be swamped by noise. Agreement between FPKM-based sequencing quantification and independent qRT-PCR measurement substantially strengthens the conclusion that the two colony types differ in their antennal gene expression.

The study also carries implications for beekeeping practice. Propolis is prized both as a hive product and as a raw material for biomedicine, with documented antimicrobial and anti-inflammatory properties, and colonies that produce more of it are valuable for breeding programs. If the olfactory gene families identified here turn out to be reliable markers of resin collection propensity, they could inform selective breeding of high-yield lines, complementing traditional selection based on colony performance alone. The morphological markers, such as antennal length and the dimensions of specific flagellomeres, might likewise serve as rapid phenotypic indicators that breeders can measure without molecular testing.

From a broader scientific perspective, the results provide a molecular basis for further research into how worker bee antennae function within honeybee colonies. The antennae regulate foraging behavior by detecting signals from nectar, pollen, and propolis sources, and the new data suggest that this regulatory role varies between colonies in ways that correlate with behavioral differentiation. The authors argue that their findings offer new perspectives on the molecular mechanisms underlying the interaction between worker bee behavioral differentiation and the olfactory system, a theme that resonates with long-standing questions in social insect biology about how genetically similar individuals come to perform such different tasks for the same colony.

The research was funded by the China Agriculture Research System under grant CARS-44-SYZ07 and by a Henan Institute of Science and Technology High-level Talent Research Project. The work was conducted by the Bee Products and Biomedicine Research and Development Team at the Henan Institute of Science and Technology, with corresponding authors Dongbiao Lü and Yunchao Kan leading the collaboration. Published open access in BMC Genomics, the study makes its transcriptomic datasets, phylogenetic analyses, and validation data available to the research community, laying the groundwork for functional studies that could pin down exactly which olfactory genes drive the dramatic differences in propolis yield between honeybee colonies.

Subject of Research: Genomics

Article Title: RNA sequencing of antennae and comparative analysis of key gene families in Apis mellifera workers exhibiting drastic differences in propolis yield

Article References: Wang, G., Zhang, Z., Lü, D., Ge, X., Wang, L., Sun, S., Guo, C., Tian, J., Zhang, P., Yang, M., Du, K., Li, D., Zhang, Z., Yu, H., & Kan, Y. (2026). RNA sequencing of antennae and comparative analysis of key gene families in Apis mellifera workers exhibiting drastic differences in propolis yield. BMC Genomics. https://doi.org/10.1186/s12864-026-13462-8

Image Credits: AI Generated

DOI: 10.1186/s12864-026-13462-8

Keywords: Apis mellifera, propolis, antennal transcriptome, odorant binding proteins, chemosensory proteins, odorant receptors, RNA sequencing, BMC Genomics, sequencing, antennae, comparative, analysis

ADAMTS15 belongs to the ADAMTS family of secreted enzymes, which carry a distinctive thrombospondin type 1 motif and cut through components of the extracellular matrix, the scaffold that surrounds cells. The enzyme participates in tissue remodeling through the hydrolysis of hyaluronic acid, a major structural molecule in that matrix. Earlier work in mice showed that ADAMTS15 is expressed during myogenesis, the process by which muscle fibers form, and that it helps regulate myoblast fusion, the step in which immature muscle cells merge into multinucleated fibers. But in goats, the mechanisms governing when and where the gene is switched on remained largely unknown. The new study set out to fill that gap by characterizing the gene’s expression profile, mapping its regulatory elements, and testing whether variants in those elements associate with growth traits in a local Chinese breed.

The team worked with Anhui White goats, a breed raised in Anhui Province, and examined ADAMTS15 expression across a range of tissues as well as in the longissimus dorsi, the large back muscle that runs along the spine and is a primary cut in meat production. Their measurements confirmed that the gene is expressed in goat muscle and revealed a progressive decline during development, with embryonic muscle tissues showing significantly higher expression than postnatal stages. This pattern mirrors what has been observed in mice and reinforces the idea that ADAMTS15 plays a role in the early construction of muscle tissue, when myoblasts are actively fusing and fibers are being assembled, rather than in the maintenance of mature muscle.

With the expression pattern established, the researchers turned to the central technical challenge of the study: locating the sequences that drive the gene’s transcription. They used dual-luciferase reporter assays, a technique in which candidate DNA fragments are cloned upstream of a luciferase gene and introduced into cells; if a fragment contains promoter or enhancer activity, the luciferase enzyme lights up, providing a quantitative readout of regulatory strength. By testing a series of overlapping fragments, they mapped the promoter activity region to a stretch between 678 and 385 base pairs upstream of the transcription start site. The enhancer screen identified five active fragments, labeled E2, E4, E6, E7, and E8, each capable of boosting transcription when tested in the reporter system.

A closer look at the mapped promoter region revealed that it overlaps with a CpG island, a cluster of cytosine-guanine dinucleotides that often serves as a landing zone for methylation, a chemical modification that can silence genes. DNA methylation analysis of this island produced an unexpected result: all 14 CpG sites within it remained unmethylated in longissimus dorsi muscle from 75-day-old fetuses, 1-day-old newborns, and 10-month-old adult goats. The complete absence of methylation across three developmental stages suggests that this CpG island functions as a stable, open regulatory element rather than a dynamic methylation-controlled switch, and it hints that other mechanisms, such as transcription factor binding, may be the primary drivers of the gene’s developmental expression pattern.

The researchers then screened for genetic variants across the promoter, the enhancer regions, and the gene itself using Sanger sequencing, the gold-standard method for reading DNA sequence directly. They identified single nucleotide polymorphisms, or SNPs, positions where a single DNA letter differs between individuals, and tested whether these variants associate with six growth traits measured in the goat population, including body weight, body length, chest circumference, and cannon circumference, the girth of the lower leg bone that reflects skeletal development. The association analysis produced a set of significant hits that span the gene’s entire regulatory architecture.

In the promoter region, two variants stood out. Promoter-SNP2 was significantly associated with body weight, while Promoter-SNP5 showed a significant association with cannon circumference. In the enhancer regions, Enhancer-SNP1 was significantly associated with both chest circumference and cannon circumference. Within the gene itself, SNP2 was significantly associated with body weight, and SNP3 showed significant associations with all six growth traits, making it the broadest hit in the study. These results indicate that variation in the sequences that control ADAMTS15 expression, as well as variation in the gene’s coding and genic regions, contributes to measurable differences in how goats grow.

To move beyond correlation, the team performed functional analysis of the two most compelling regulatory variants. They found that Promoter-SNP2 genotypes significantly affected promoter activity in the reporter assay, meaning that the different DNA letters at this position change how strongly the promoter drives transcription. Similarly, Enhancer-SNP1 genotypes significantly affected enhancer activity, altering the strength of the enhancer fragment in which it sits. This combination of statistical association with growth traits and direct functional effects on regulatory activity provides strong evidence that these variants are not merely linked to growth differences but actively participate in producing them, by tuning the expression level of ADAMTS15 during muscle development.

The study offers a detailed foundation for understanding the transcriptional regulatory mechanisms of ADAMTS15 in goat muscle development, and it demonstrates a workflow that other researchers can apply to any gene of interest in livestock: characterize expression, map regulatory elements with reporter assays, examine methylation status, screen for variants with Sanger sequencing, and test both associations with traits and functional effects on regulatory activity. The identified variants may serve as potential molecular markers for goat molecular breeding, allowing breeders to select animals carrying favorable alleles at these positions to improve growth performance. For a breed like the Anhui White goat, which is valued as a local genetic resource, such markers could accelerate genetic improvement while preserving the breed’s distinctive characteristics.

Beyond the practical implications for breeding, the work adds to a growing body of evidence that the ADAMTS family of enzymes, long studied for their roles in tissue remodeling and matrix biology, also participates in the genetic architecture of growth in livestock. The finding that a gene involved in myoblast fusion in mice shows a developmentally regulated expression pattern and trait-associated variants in goats suggests that the underlying biology of muscle formation is conserved across species, even as the specific regulatory sequences evolve. As genome editing and precision breeding technologies mature, understanding the regulatory elements and variants that govern genes like ADAMTS15 will become increasingly important for translating basic biology into improved animal health and productivity.

Subject of Research: Genomics

Article Title: Identification of regulatory elements and association analysis of genetic variants with growth traits of goat ADAMTS15 gene

Article References: Zhang, S., Yan, N., Peng, J., Chen, J., Zhu, L., Li, S., Ji, K., Xu, H., Huang, D., & Ling, Y. (2026). Identification of regulatory elements and association analysis of genetic variants with growth traits of goat ADAMTS15 gene. BMC Genomics. https://doi.org/10.1186/s12864-026-13467-3

Image Credits: AI Generated

DOI: 10.1186/s12864-026-13467-3

Keywords: goat, ADAMTS15, gene regulation, promoter, enhancer, growth traits, genetic variation, molecular breeding, BMC Genomics, Identification, regulatory, elements

The answer matters more than it might first appear. Quantitative real-time PCR, or RT-qPCR, is the workhorse technique for measuring gene expression. It works by amplifying target DNA sequences and measuring how many cycles of amplification are needed before the signal crosses a detection threshold, a value known as the Ct. Lower Ct values mean more of the original messenger RNA was present, and therefore the gene was more active. But the raw numbers from an RT-qPCR experiment are meaningless on their own, because they are influenced by everything from the amount of starting material to the efficiency of the enzymatic reaction. To correct for these technical variables, researchers normalize their results against internal reference genes, sometimes called housekeeping genes, which are assumed to be expressed at roughly constant levels regardless of experimental conditions.

The problem, and the reason this study exists, is that the assumption of constant expression is frequently wrong. Reference genes are real genes with real biological functions, and under stress their activity can shift dramatically. If a chosen reference gene changes its expression in response to salt stress or drought, every downstream measurement becomes distorted, sometimes reversing the apparent direction of a result. A comprehensive literature search conducted by the team revealed that no suitable reference genes had ever been reported for Rubia cordifolia, either in China or internationally. For a species routinely used in clinical practice as a traditional Chinese medicinal herb, and increasingly studied at the molecular level for the pharmacologically active compounds it produces, this was a conspicuous gap.

To fill it, the researchers assembled a panel of twelve candidate reference genes and profiled their expression across a demanding set of experimental conditions. Samples were taken from various tissues of the plant, and whole plants were subjected to five distinct abiotic stress treatments: salt stress induced by sodium chloride, dehydration stress, and treatments with three chemically distinct signaling molecules that plants use to orchestrate stress responses, namely salicylic acid, methyl jasmonate, and abscisic acid. This design is important because reference gene stability is condition-specific. A gene that holds steady in leaves may wobble in roots, and a gene unperturbed by salt may fluctuate during dehydration. By spanning tissues and stress types, the study aimed to identify genes whose stability is robust rather than accidental.

Measuring the candidates, however, is only half the task. The raw stability data must be interpreted, and different analytical algorithms can weight the evidence differently, sometimes reaching different conclusions from the same numbers. The team therefore ran their dataset through four widely used algorithms. geNorm evaluates pairwise variation among candidate genes and calculates an expression stability measure, on the logic that truly stable genes should vary in concert with one another. NormFinder builds a statistical model that estimates expression variation while accounting for subgroup structure, such as different tissues or treatments. BestKeeper takes a correlation-based approach, comparing each candidate directly against an index constructed from the whole panel. The comparative delta Ct method, as its name suggests, compares relative expression values between pairs of candidates to rank their stability.

Each algorithm contributed to a consensus picture. Across the tissues and stress conditions, hnRNP, UBQ, HIS3.3, PP2A, and PTBP2 emerged as the most stable of the twelve candidates, while ACT7 and EF-1α proved to be the least stable. That last finding carries a cautionary lesson for the wider field: ACT and EF-1α are among the most frequently used reference genes in plant studies precisely because of their housekeeping reputations, yet reputations do not transfer automatically across species or conditions. To synthesize the four individual analyses, the team also employed RefFinder, an integrative tool that aggregates the rankings from geNorm, NormFinder, BestKeeper, and the delta Ct method into a single comprehensive assessment. The integrated analysis confirmed hnRNP, HIS3.3, PTBP2, PP2A, and RPL5 as the top five most stable reference genes overall.

Identifying candidates is one thing; proving they actually improve experimental accuracy is another. For validation, the researchers turned to a target gene called MVK, which encodes mevalonate kinase, an enzyme embedded in the mevalonate pathway that underpins the biosynthesis of isoprenoid compounds. This pathway is directly relevant to the pharmacological interest in Rubia cordifolia, since many of the plant’s bioactive secondary metabolites are built from isoprenoid precursors. By measuring MVK expression and normalizing the results against each candidate reference gene in turn, the team could see in practice how the choice of reference gene changed the apparent answer. Normalization against the most stable references produced consistent, reliable MVK profiles, confirming that the screened genes are not merely statistically stable but functionally dependable as baselines for real experiments.

The technical significance of the work extends beyond a single species. Reference gene validation studies follow a well-established template, but each new species requires its own screen because orthologous genes can behave very differently across evolutionary lineages, and even within a species, stability varies with tissue, developmental stage, and stressor. By establishing a validated reference set for Rubia cordifolia, the study removes a fundamental source of error from all future gene expression work in this plant, whether the goal is understanding how the herb responds to salinity in cultivation, deciphering the enzymatic routes to its medicinal compounds, or engineering enhanced production of valuable metabolites. The authors describe the work as a foundational step toward reliable functional genomics research in the species, paving the way for accurate elucidation of gene function.

The study also arrived with appropriate scientific and ethical diligence. The work was funded by the High-Level Key Discipline Construction Project of Traditional Chinese Medicine of the National Administration of Traditional Chinese Medicine under grant number SATCM-RH[2023]85. In collecting plant and seed samples, the researchers state that they strictly abided by the policy statements of the International Union for Conservation of Nature concerning research involving endangered species, as well as the guiding principles of the Convention on International Trade in Endangered Species of Wild Fauna and Flora. The authors declared no competing interests, and the article is published open access under a Creative Commons license, making the data and conclusions freely available to researchers worldwide.

For practitioners of plant molecular biology, the practical takeaway is straightforward but consequential: never assume a reference gene is stable, and never borrow one on reputation alone. The Henan team’s finding that ACT7 and EF-1α, two of the plant kingdom’s most popular normalizers, rank as the least stable candidates in Rubia cordifolia is a vivid reminder that every species tells its own molecular story. With hnRNP, UBQ, HIS3.3, PP2A, PTBP2, and RPL5 now validated as dependable anchors, researchers studying this ancient medicinal herb can finally measure its genes with confidence, and the quiet, meticulous science of normalization has delivered one of its characteristically unglamorous but genuinely enabling victories.

Subject of Research: Validation of stable internal reference genes for RT-qPCR gene expression analysis in Rubia cordifolia under abiotic stress

Article Title: Screening and validation of suitable internal reference genes for gene expression analysis in Rubia cordifolia L. under different abiotic stresses

Article References: Wenwen, L., Ke, W., Mingzhu, X., Conglong, L., & Suiqing, C. (2026). Screening and validation of suitable internal reference genes for gene expression analysis in Rubia cordifolia L. under different abiotic stresses. BMC Genomics. https://doi.org/10.1186/s12864-026-13415-1

Image Credits: AI Generated

DOI: 10.1186/s12864-026-13415-1

Keywords: Rubia cordifolia, RT-qPCR, internal reference genes, gene expression, abiotic stress, geNorm, NormFinder, BestKeeper, RefFinder, BMC Genomics, traditional Chinese medicine, functional genomics

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: Masand, M., Kaur, R., Mukthar, M. M., Youssef, M. S., Walkowiak, S., Duncan, R. W., Zelmer, C., McVetty, P. B. E., & Chawla, H. S. (2026). Long-read genome assembly of high-erucic-acid Brassica napus reveals candidate structural variants in fatty acid metabolism genes. BMC Genomics. 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 newly assembled mitochondrial genome of P. tuberculosus measures 115,079 base pairs with a GC content of 26.1 percent, a size and composition that fit comfortably within the range reported for closely related Phellinus species. For fungi, whose mitochondrial genomes can vary enormously in size and architectural complexity, this figure immediately signals both conservation and room for the kind of expansion that makes fungal mitogenomics so intriguing. Unlike the compact animal mitochondria most people learn about in school, fungal mitochondrial genomes routinely swell with introns, repeats and intergenic spacers, and the Phellinus genomes are no exception.

Annotation of the P. tuberculosus mitogenome revealed a core genetic toolkit that appears to be a hallmark of the group: fourteen core protein-coding genes, twenty-seven transfer RNA genes and two ribosomal RNA genes. Those fourteen protein-coding genes correspond to the standard set required for oxidative phosphorylation, the electron transport chain machinery that mitochondria use to generate cellular energy. Their consistent presence across the compared species underscores how tightly evolution preserves the essential respiratory apparatus, even as the surrounding genomic landscape churns with rearrangements, duplications and mobile elements.

Beyond the genes themselves, the researchers catalogued 79 simple sequence repeats along with multiple larger repeat sequences scattered across the genome. Simple sequence repeats, sometimes called microsatellites, are short tandemly repeated motifs that mutate rapidly and can serve as informative markers for population-level studies. Larger repeat elements, by contrast, are the raw material of genome rearrangement: because recombination between repeats can invert, shuffle or duplicate genomic segments, their abundance helps explain why fungal mitochondrial genomes can differ dramatically in gene order even between closely related species. The repeat landscape of P. tuberculosus thus provides both a record of past genomic turbulence and a potential engine of future change.

One of the most technically interesting findings concerns RNA editing, the post-transcriptional process that alters the sequence of RNA transcripts relative to the DNA template. The team predicted a total of 42 C-to-U RNA editing sites across the mitochondrial protein-coding transcripts, with the nad2 gene carrying the highest number of editing sites. C-to-U editing in fungal mitochondria changes cytidine residues to uridine in the RNA, sometimes correcting codons so that the final proteins match functional expectations even when the gene sequence alone would suggest otherwise. The concentration of editing sites in nad2, which encodes a subunit of the respiratory chain’s NADH dehydrogenase complex, hints that this gene in particular may rely on post-transcriptional correction to produce a properly functioning protein.

Comparative genomic analysis with the two previously characterized Phellinus mitochondrial genomes showed that the P. tuberculosus lineage shares conserved mitochondrial genome features with its relatives, from gene content to overall organization. At the same time, the study’s title points to dynamic genome evolution within Hymenochaetaceae more broadly, reflecting the tension that defines fungal mitogenomics: the core is stable, but the periphery is restless. Introns, in particular, are notorious agents of mitochondrial genome expansion in basidiomycetes, and their gain and loss across lineages can inflate genome size without adding protein-coding capacity. The Phellinus comparison captures this dynamic in miniature, showing how genomes that carry the same essential genes can nonetheless diverge substantially in size and structure.

The phylogenomic component of the study used mitochondrial sequence data to test how P. tuberculosus relates to other members of its genus and family. The analyses showed that P. tuberculosus clusters with the two previously characterized Phellinus species in a strongly supported clade, a result that both validates the mitochondrial genome as a phylogenetic marker and reinforces the monophyly of the genus as currently sampled. Strong statistical support for this grouping matters because Phellinus, like many wood-inhabiting polypore lineages, has a taxonomic history complicated by morphological similarity, and molecular data have repeatedly reshaped how these fungi are classified. A well-supported mitochondrial phylogeny offers an independent line of evidence for species boundaries and evolutionary relationships.

Why does a mitochondrial genome from a single mushroom species matter beyond taxonomy? Phellinus species are white-rot fungi, meaning they possess the enzymatic machinery to degrade lignin, the tough aromatic polymer that gives wood its rigidity. That capability has enormous ecological significance in forests, where white-rot fungi drive carbon cycling by decomposing dead wood, and equally significant biotechnological potential, since lignin-degrading enzymes are of interest for biofuels, bioremediation and green chemistry. Many Phellinus species also produce medicinally functional compounds and nutrients, and several have been used in traditional medicine systems for centuries. Understanding the genetic architecture of these organisms, including their mitochondrial genomes, contributes to a fuller picture of their biology, their evolution and their potential applications.

The authors frame their work as filling a genuine gap: despite the medicinal importance of Phellinus, a comprehensive genus-level analysis of mitochondrial genomes had been unavailable until now. By providing the first complete mitogenome for P. tuberculosus and comparing it with its congeners, the study establishes a baseline for future phylogenetic, ecological and biotechnological investigations of medicinally important species. The annotated genome, with its catalogued repeats and predicted RNA editing sites, becomes a resource that other researchers can mine for markers, for insights into mitochondrial gene expression, and for reconstructing the evolutionary history of a lineage that shapes forests and pharmacies alike.

The study also illustrates a broader trend in genomics: mitochondrial genomes, once treated as an afterthought compared to nuclear genomes, have become workhorses for evolutionary biology in fungi. Their relatively small size makes them cheap to assemble completely, their uniparental inheritance simplifies comparisons between lineages, and their combination of conserved cores and variable architecture makes them ideal for resolving relationships at shallow and intermediate evolutionary depths. For the Hymenochaetaceae, a family that includes numerous pathogens and decomposers of standing trees as well as medicinal species, each new complete mitogenome adds a data point that sharpens both phylogeny and our understanding of how fungal genomes evolve. The P. tuberculosus genome, published open access and citable with a permanent DOI, now joins that growing library, and the strongly supported Phellinus clade it anchors suggests that more genus-level mitochondrial comparisons across Hymenochaetaceae will be well worth watching.

Subject of Research: Comparative mitochondrial genomics of the medicinal fungus Phellinus tuberculosus within the family Hymenochaetaceae

Article Title: Comparative mitochondrial genomics of the medicinal mushroom Phellinus tuberculosus reveals conserved core genes and dynamic genome evolution within Hymenochaetaceae

Article References: Jin, C., Yuan, M.-X., Ma, J.-X., Ji, Y.-Y., Tang, L.-X., Li, X., Chen, M.-Y., Si, J., & Cui, B.-K. (2026). Comparative mitochondrial genomics of the medicinal mushroom Phellinus tuberculosus reveals conserved core genes and dynamic genome evolution within Hymenochaetaceae. BMC Genomics. https://doi.org/10.1186/s12864-026-13454-8

Image Credits: AI Generated

DOI: 10.1186/s12864-026-13454-8

Keywords: Phellinus tuberculosus, mitochondrial genome, Hymenochaetaceae, comparative genomics, RNA editing, white-rot fungi, medicinal mushrooms, phylogenomics, simple sequence repeats, fungal evolution, BMC Genomics, genome annotation

The research focused on four fundamental measurements used in livestock breeding to assess body conformation: body weight, body height, body oblique length, and chest circumference. Each of these traits captures a different aspect of an animal’s physical development. Body weight reflects overall mass and is often the single most economically important figure for meat production. Body height and body oblique length describe skeletal frame size, while chest circumference provides a proxy for body depth and volume. Together, they form the standard toolkit of quantitative measures that breeders use to judge whether an animal is maturing well. Understanding which genes influence these numbers is crucial for designing breeding programs that can improve productivity while preserving the breed’s adaptation to its extreme environment.

To gather the raw material for the analysis, the researchers sampled 536 Tianzhu White Yaks representing three distinct growth stages: animals at one year of age, animals at two and a half years, and animals aged four years or older. This staged design matters because growth is not a single event but a trajectory. Genes that govern early skeletal development may not be the same ones that control later muscle deposition and fat accumulation, and including animals from different life stages gives the study a window onto that changing biology. Blood samples were collected with ethical approval from the Lanzhou Institute, in strict compliance with China’s Guidelines for Ethical Treatment of Experimental Animals, and the four body measurements were recorded for each animal.

The genomic strategy the team chose was low-coverage resequencing, a technique that reads the genome at relatively shallow depth across hundreds of individuals rather than deeply sequencing a small number. The advantage is scale: with enough animals, statistical methods can reconstruct each individual’s full genotype computationally in a process called genotype imputation. After imputation and rigorous quality control, the study retained a staggering total of 24,383,269 single nucleotide polymorphisms, or SNPs, the single-letter differences in DNA sequence that pepper the yak genome. That depth of marker information, nearly two dozen million variants, is what gives a GWAS the resolving power to home in on genomic regions that genuinely influence a trait rather than hitchhiking along with them.

Armed with this dense variant map and the body measurements, the researchers applied a mixed linear model for the association testing. Mixed linear models are the workhorse of modern GWAS because they simultaneously test the link between each SNP and a trait while statistically accounting for the complex family structure and relatedness that are almost unavoidable in livestock populations. Without such correction, related animals sharing large stretches of genome would generate false associations. The team also used standard quality checks such as principal component analysis to characterize population structure and quantile-quantile plots to verify that their model was behaving as expected, with the observed test statistics aligning closely with what would be predicted under their assumptions.

The scan paid off. Across the four traits, the study identified a total of 366 significant SNP loci, which the researchers annotated to 32 candidate genes. For body weight, the candidates included PLCB1, SAR1B, and SRSF6. Body height pointed toward BRPF3, MAPK13, and ZNF76. Body oblique length was associated with variants near ABCD3 and NTM, while chest circumference flagged LMNTD1 and ACVR1. The latter gene is particularly notable to anyone familiar with developmental biology: ACVR1 encodes a receptor in the bone morphogenetic protein signaling pathway, a family of pathways long known to regulate skeletal growth, bone formation, and patterning during development. Its appearance in a chest circumference analysis, a measure closely tied to skeletal dimensions, fits neatly with what biologists already understand about how the body frame is built.

To move from a list of gene names toward biological meaning, the team performed functional enrichment analysis, asking whether the candidate genes cluster in particular cellular processes more than chance would predict. Several themes emerged. Multiple genes were implicated in the negative regulation of the G1/S transition of the mitotic cell cycle, the checkpoint that decides whether a cell commits to dividing. That process is fundamental to how organisms grow, since the number and timing of cell divisions in developing tissues directly shape final organ and body size. Another enriched process was fatty acid beta-oxidation, the metabolic pathway by which cells break down fats to generate energy, a fitting connection for traits like body weight and chest circumference that are influenced by both lean tissue and fat reserves, especially in an animal that must survive long, cold winters on sparse forage.

One of the most intriguing findings involved three genes, TEAD3, ZNF76, and RBL1, that were significantly enriched in the transcription regulator complex, a cellular component consisting of assemblies of proteins that control which genes are switched on or off. Transcription regulators sit high in the hierarchy of cellular control, and their involvement suggests that variation in these loci could influence cell proliferation, differentiation, and tissue formation, the underlying processes that ultimately determine how large and well-formed an animal becomes. RBL1, for instance, is part of the retinoblastoma family of proteins that act as gatekeepers of the cell cycle, while TEAD proteins are effectors of the Hippo signaling pathway, an evolutionarily conserved organ-size control system. Such connections, if confirmed by functional studies, would tie the breed’s growth biology to some of the deepest mechanisms known in developmental biology.

The practical implications reach well beyond the laboratory. Tianzhu White Yaks are a valuable genetic resource, prized on the plateau for their white coat and their role in the livelihoods of herding communities, and the study’s authors emphasize that identifying markers associated with growth traits is crucial for breeding. Marker-assisted selection and genomic selection both depend on precisely this kind of catalog: a validated set of DNA variants that predict how an animal will grow before it reaches maturity. By genotyping young calves for the significant SNPs identified here, breeders could, in principle, make earlier and more accurate selection decisions, accelerating genetic improvement while maintaining the genetic diversity that makes the breed resilient. The authors note that the identified SNPs and candidate genes will facilitate the development and utilization of Tianzhu White Yak resources.

The study also carries broader scientific weight. Yaks occupy a unique position among livestock as the large mammals best adapted to life above 3,000 meters, and every well-powered GWAS in the species adds to a comparative framework that spans cattle, yak, and their hybrids. Because the analysis spans animals of different ages, the same dataset could support future work on how gene effects change across development, and the 24 million SNP resource itself is a legacy asset for the yak genetics community. For now, the 366 loci and 32 candidate genes stand as a hypothesis list, the starting point for the functional experiments the researchers say should follow. From a blood draw on a windswept plateau to nearly two dozen million genetic letters, the study shows how modern genomics can translate the centuries-old practice of herding into a molecular map of what makes an animal grow.

Subject of Research: Genome-wide association analysis of growth traits in the Tianzhu White Yak

Article Title: Genome-wide association analysis of growth traits in Tianzhu White Yak

Article References: Yang, G., Qi, X., La, Y., Ma, X., Liu, Y., Wang, Z., Luo, W., Ren, W., Zhang, Z., Chu, M., Wu, X., Guo, X., Qi, W., & Liang, C. (2026). Genome-wide association analysis of growth traits in Tianzhu White Yak. BMC Genomics. https://doi.org/10.1186/s12864-026-13444-w

Image Credits: AI Generated

DOI: 10.1186/s12864-026-13444-w

Keywords: Tianzhu White Yak, GWAS, growth traits, SNPs, low-coverage sequencing, genotype imputation, livestock breeding, BMC Genomics, body conformation, ACVR1, transcription regulation, Tibetan Plateau