Single nucleotide polymorphisms, or SNPs, are the workhorses of contemporary genomics. These single-letter variations scattered across the genome allow scientists to fingerprint individuals, reconstruct family trees, trace population ancestry and, most importantly for breeders, predict which juvenile fish carry the genetic potential for fast growth, efficient feed conversion or disease resistance. High-density SNP arrays containing tens of thousands of markers deliver all of this information, but at a price per sample that becomes prohibitive when a breeding program needs to genotype thousands or tens of thousands of fish each generation. The cost barrier has been a major constraint on the large-scale application of genomic approaches in large yellow croaker, and the same problem afflicts many other aquaculture species where profit margins per animal are thin.

The new array was not built from scratch. The team started with NingXin-III, a previously established 55,000-marker array for the species, and systematically distilled it down to roughly 10,000 of the most informative SNPs. The selection strategy was deliberate rather than arbitrary: representative SNPs were chosen from haplotype blocks, the stretches of DNA that tend to be inherited together, so that each block would be covered by at least one marker capturing its variation. On top of that skeleton, the researchers supplemented the panel with SNPs at loci previously associated with economically important traits, ensuring the compact array retained predictive power where it counts most for selection decisions.

The technical performance of the trimmed panel proved remarkably strong. Across validation testing, the detection rate and genotype concordance, a measure of how faithfully the array calls each genetic variant compared with a reference standard, both exceeded 99 percent. The retained SNPs were also distributed evenly across all 24 chromosomes of the croaker genome, a critical property for any array intended to support genome-wide analyses, because clustering of markers on particular chromosomes would leave large genomic regions invisible to the platform. Even coverage means the 10 K panel can serve as a genuinely representative sample of the species’ genetic variation rather than a patchy snapshot.

Where the array truly earns its keep is in the suite of applied tests the researchers ran. In population genetic analyses, the panel effectively distinguished different croaker populations, a capability central to germplasm identification, the process of verifying that a batch of fry or broodstock genuinely belongs to the improved strain a hatchery claims it is. Machine learning classifiers built on the 10 K dataset, including logistic regression and random forest models evaluated with metrics such as the Matthews correlation coefficient and area under the receiver operating characteristic curve, achieved an accuracy exceeding 0.99 for germplasm identification, with results highly consistent with those obtained from the full 55 K array. In other words, cutting the marker count by more than 80 percent cost almost nothing in classification power.

Pedigree tracking delivered equally striking results. The array achieved 100 percent accuracy in parentage assignment, pedigree reconstruction and sex prediction. For breeding companies, parentage verification is not a luxury: without reliable family records, selective breeding programs cannot estimate heritabilities, control inbreeding or make accurate selection decisions, and mixed-up family assignments can quietly erode years of genetic gain. A cheap molecular pedigree check that works perfectly on every tested individual removes one of the most tedious and error-prone bookkeeping tasks in fish hatcheries, where thousands of families may be reared in shared tanks and physical tagging is impractical at scale.

Genomic selection, the most ambitious application, required an extra computational step. Because a 10 K panel samples only a fraction of the genome, the researchers used genotype imputation, statistical inference that fills in the untyped markers based on patterns of correlation with the typed ones, to reconstruct a denser dataset before running genomic prediction with the genomic best linear unbiased prediction (GBLUP) method. Imputation accuracy ranged from 0.817 to 0.964 depending on the scenario. After imputation, the 10 K panel retained genomic prediction performance close to that of the 55 K panel, although the authors are candid that the absolute predictive ability of both datasets was moderate, a reminder that prediction accuracy for complex traits such as body weight, feed conversion ratio and critical swimming speed depends on many factors beyond marker density, including reference population size and trait architecture.

One especially practical metric for breeders is how well a low-density array identifies the top performers in a population, since selection programs act on the best animals rather than on average prediction accuracy across the board. Here the imputed 10 K data held up well: among individuals ranked in the top 10 percent by genomic estimated breeding values (GEBV), the compact panel overlapped with the 55 K panel by 86.21 to 87.93 percent. That means a hatchery using the cheap array would catch the overwhelming majority of the same elite fish it would have flagged with the expensive one, at a fraction of the genotyping cost, making whole-population genomic selection economically viable for the first time in this species.

Beyond its immediate value for large yellow croaker, the study offers a template for other aquaculture breeding programs wrestling with the same cost calculus. The design logic, starting from a validated high-density array, pruning markers by haplotype block, enriching for trait-associated loci and validating across identification, parentage and prediction tasks, is directly transferable to other farmed fish, shrimp and shellfish. The authors position NingXin-IV as an efficient and cost-effective tool for germplasm evaluation, parentage verification, breeding strain management and large-scale genomic selection, and explicitly frame it as a useful reference for developing low-density platforms elsewhere in aquaculture. The research was supported by funding including the National Science Fund for Distinguished Young Scholars and the China Postdoctoral Science Foundation, with fish maintained at Ningde Fufa Fisheries Company Limited under institutional animal care protocols, fin clips collected after anesthesia with MS-222 and every sampled fish returned alive to its culture tank.

For an industry built on a fish that has been farmed in China for decades, the arrival of a sub-cent-per-marker genotyping platform may prove to be a quiet revolution. Cheaper genotyping means more animals screened per generation, faster genetic progress for traits farmers care about, and stronger protection of certified breeding lines against fraudulent substitution. It also means genomic tools long confined to well-funded laboratories can move into the routine operations of coastal hatcheries, where the daily decisions about which fish to keep as broodstock ultimately shape the future of the species. If the NingXin-IV experience is any guide, the future of aquaculture genetics may belong not to ever-larger arrays, but to smartly designed small ones that deliver nearly all the answers at a price the industry can actually pay.

Subject of Research: Development of a low-density 10 K liquid SNP array for genetic improvement and genomic selection in large yellow croaker

Article Title: Development and application of a low-density 10 K liquid SNP array for genetic improvement in large yellow croaker (Larimichthys crocea)

Article References: Li, Y., Wang, J., Zhao, J., Ke, Q., Jiang, P., Zeng, J., Weng, H., Pu, F., Zhou, T., Li, N., & Xu, P. (2026). Development and application of a low-density 10 K liquid SNP array for genetic improvement in large yellow croaker (Larimichthys crocea). BMC Genomics. https://doi.org/10.1186/s12864-026-13387-2

Image Credits: AI Generated

DOI: 10.1186/s12864-026-13387-2

Keywords: large yellow croaker, SNP array, genomic selection, aquaculture, genotyping, Larimichthys crocea, parentage assignment, germplasm identification, genotype imputation, BMC Genomics, selective breeding, machine learning

The study, led by Desalegn Amsalu Gelayie of the Amhara Agricultural Research Institute in Bahir Dar, together with colleagues at Bahir Dar University and the University of Gondar, began in the most straightforward way possible for a modern protein investigation: by retrieving the amino acid sequence of zebu α-lactalbumin from the UniProt database in FASTA format. That sequence was then submitted to NCBI BLAST, the workhorse similarity-search tool of computational biology, to confirm its identity and find related proteins. From there, the researchers layered on a battery of analytical methods: physicochemical characterization with ExPASy ProtParam, secondary structure prediction, multiple sequence alignment, homology modeling, and rigorous quality evaluation of the resulting three-dimensional model, which they visualized and inspected using PyMOL and Discovery Studio.

The physicochemical profile that emerged paints a picture of a compact, stable protein. The molecular weight came out at 16,201.75 daltons, consistent with the small, single-chain nature of α-lactalbumin. The isoelectric point of 4.79 marks the protein as mildly acidic, a familiar feature of whey proteins that influences their solubility and behavior during milk processing. The extinction coefficient of 27,960 reflects the aromatic residue content that allows the protein to be quantified reliably by ultraviolet spectroscopy. Most telling, perhaps, is the instability index of 27.44. In the Gram-scale logic of the Guruprasad method, values below 40 predict a stable protein in solution, and 27.44 sits comfortably within that range, indicating that the zebu α-lactalbumin is unlikely to unravel under physiological conditions.

Secondary structure predictions, generated by multiple algorithms, showed the expected mixture of α-helices, β-strands, and random coil, though the exact proportions varied depending on which prediction tool was applied. This kind of tool-to-tool variation is a well-known feature of secondary structure prediction, particularly for small proteins like α-lactalbumin that adopt a molten-globule-like state under some conditions. The native fold of α-lactalbumin is a compact α-plus-β architecture, and the consensus of the predictions was consistent with that canonical arrangement, giving the researchers confidence that their downstream modeling rested on a sound structural foundation.

One of the more technically interesting components of the analysis was the prediction of intrinsic disorder, the tendency of certain sequence segments to remain flexible and unstructured even in an otherwise folded protein. The analysis identified four disordered segments, with residues at positions 1 to 3, 82 to 89, and 139 to 142 predicted to be disordered. The remainder of the chain was predicted to be ordered, meaning those residues hold stable positions within the folded structure and contribute directly to its stability. Short disordered stretches at the termini and in a single internal loop are not unusual for small globular proteins and can even be functionally important, providing the flexibility needed for the conformational changes that α-lactalbumin undergoes when it partners with β-1,4-galactosyltransferase to form the lactose synthase complex.

The centerpiece of the study was the homology model. Using known structures as templates, the researchers built a three-dimensional model of the zebu protein and validated it against standard quality metrics. The model showed strong structural similarity to the bovine taurine α-lactalbumin structure deposited in the Protein Data Bank under the identifier 1f6s.6.A, with a sequence identity of 99.19 percent and a query coverage of 87 percent. In plain terms, the zebu and taurine proteins are nearly identical at the amino acid level, and the modeled zebu structure reproduces the taurine fold with high fidelity. The valid metrics and qualified prediction characteristics reported for the model indicate that it meets the accepted standards for a reliable comparative structure.

Perhaps the most consequential structural feature confirmed in the model is the calcium-binding site. α-Lactalbumin is a calcium metalloprotein, and its bound Ca2+ ion is essential for the protein’s correct folding and biological activity. The predicted zebu structure contained a calcium ion coordinated by five residues within a four-angstrom radius: lysine 98, aspartate 101, aspartate 103, aspartate 106, and aspartate 107 on chain A. The cluster of acidic aspartate residues is exactly what one expects for a high-affinity calcium site, and the model also displayed the intra-protein interactions and metal complexes that stabilize the native fold. The conservation of this site across the zebu sequence is a strong signal that the protein’s functional chemistry is intact.

Why does this matter beyond the structural biology community? Milk yield in cattle is fundamentally limited by lactose synthesis, because lactose is the principal osmole that draws water into mammary secretions. If α-lactalbumin were structurally compromised in zebu cattle, it could contribute to the lower milk volumes typically observed in these breeds compared with specialized Bos taurus dairy cattle. The new results suggest that this is not the case: the zebu protein is stable, its disordered regions are limited, and its conserved regions essential for lactose synthase binding are preserved, indicating functional similarity with the taurine protein. The differences in milk production between the two cattle types are therefore more plausibly attributed to the broader genetic and environmental factors that the authors highlight, rather than to any defect in this key milk protein.

The practical implications flow from that conclusion. A structurally sound, functionally equivalent α-lactalbumin in zebu cattle means that breeding programs aimed at improving milk production in tropical regions can work with the protein’s existing architecture rather than needing to compensate for structural shortcomings. It also opens the door to biotechnological applications, from optimizing milk processing characteristics to potentially engineering or supplementing the protein in dairy products. The authors are careful, however, to frame their work as a foundation rather than a final word. They recommend follow-up studies using molecular dynamics simulations to watch the protein move over time, functional pathway analysis to place it in its metabolic context, and protein network analysis to map its interaction partners within the mammary cell.

Published as an open-access paper in BMC Genomics on 2 October 2026, the study is a reminder of how much can be learned from sequence databases, careful modeling, and disciplined validation, even without a laboratory full of crystallography equipment. For the millions of farmers who depend on zebu cattle across Africa and South Asia, the message embedded in this small protein is quietly encouraging: the molecular machinery of lactose synthesis in their herds is fundamentally sound. What remains is to translate that structural assurance into the breeding strategies and biotechnological tools that could help close the persistent gap between the milk yields of the tropics and those of the temperate-zone dairy industry.

Subject of Research: Computational structural modeling of the α-lactalbumin milk protein in zebu cattle

Article Title: In silico analysis and structural modeling of zebu cattle α-lactalbumin milk protein

Article References: Gelayie, D. A., Kerisew, B., Gessese, A. T., & Bergushe, Y. (2026). In silico analysis and structural modeling of zebu cattle α-lactalbumin milk protein. BMC Genomics. https://doi.org/10.1186/s12864-026-13388-1

Image Credits: AI Generated

DOI: 10.1186/s12864-026-13388-1

Keywords: α-lactalbumin, zebu cattle, Bos indicus, milk protein, in silico analysis, homology modeling, protein structure, lactose synthesis, calcium binding, dairy breeding, bioinformatics, BMC Genomics

The study, published on 2 October 2026, tackles a problem that has limited previous work in the field. Although H. manillensis has recently emerged as a key genomic reference among medicinal leeches, researchers had not systematically combined whole-genome read mapping with RNA sequencing data at the level of individual samples. The new workflow does exactly that, layering together genome alignment, transcriptome alignment and counting, DNA barcoding to confirm species identity, phylogenetic placement, reference-guided cataloging of known antithrombotic genes, transcript-supported discovery of new candidates, and targeted sequence comparison of selected loci. By integrating these strands of evidence, the authors could distinguish genes that are genuinely present and expressed from those that merely appear in a reference annotation, a distinction that matters enormously when the goal is to prioritize candidates for drug development.

A crucial early step was making sure the samples really were what the researchers thought they were. Species misidentification is a surprisingly common hazard in leech research, where closely related species can look nearly identical. The team used the mitochondrial cytochrome c oxidase subunit I gene, the standard DNA barcode for animals, to verify all four samples. BLAST searches against reference databases returned a top hit of MN882684 with 99.842 percent identity across 1,266 base pairs for every sample, and phylogenetic analysis placed each one squarely within the H. manillensis reference clade. That confirmation matters because any downstream claim about antithrombotic gene content would be undermined if the animals had been mislabeled, and it gives future users of the dataset confidence that the gene catalog truly describes this species.

The sequencing data itself was of high quality. Whole-genome mapping rates ranged from 72.17 percent to 92.75 percent across the four samples, with mean sequencing depths between 108.0-fold and 134.4-fold, deep enough to support confident read mapping and variant calling. The RNA-seq alignment rates ranged from 67.62 percent to 70.49 percent, providing a solid transcriptomic backbone for expression analysis. With these datasets in hand, the researchers turned to the central question: how many antithrombotic genes does the Asian buffalo leech carry, and which of them are actually switched on?

The reference-guided analysis recovered 72 known antithrombotic genes, matching the published annotation framework for the species and confirming that the pipeline was reproducing established results. The more exciting outcome came from the transcript-supported discovery stage, which added 38 new candidate transcripts. Of these, 30 represent novel isoforms of already-known gene loci, flagged by the annotation tool gffcompare with class code “j”, meaning they are splice variants or alternative transcripts of characterized genes. The remaining 8 are overlap-type candidates with class codes “o” and “x”, indicating transcripts that overlap known loci in ways that suggest either alternative structures or potentially distinct genes sharing genomic space. Notably, after filtering, no high-confidence class code “u” candidates, which would represent entirely unknown intergenic transcripts, survived, suggesting that the antithrombotic repertoire of this leech is largely built from elaborations of known gene families rather than from wholly novel loci hiding between annotated genes.

Among the gene families, one stood out for its sheer size. The H. manillensis elastase inhibitor family, abbreviated HMEI, was the largest supported family in the catalog, followed by hirustasin-related genes. Elastase inhibitors are of particular interest because elastase plays roles in inflammation and tissue remodeling as well as in clot-related processes, so a leech family that has expanded this class of inhibitors may be modulating more than just coagulation. Hirustasin, meanwhile, is a known serine protease inhibitor first characterized from the medicinal leech, and its relatives in H. manillensis reinforce the picture of a leech that deploys multiple independent inhibitory strategies against its host’s hemostatic system. The expansion of these families points to strong evolutionary pressure on blood-feeding animals to maintain a diverse and redundant toolkit for keeping blood liquid.

Expression data revealed a strikingly uneven picture across the four samples. The most abundant antithrombotic transcript overall was lefaxin_Hman2, a member of a family of anticoagulant proteins, yet expression of other loci varied dramatically from animal to animal, with clear sample-specific patterns among the HMEI genes in particular. This heterogeneity suggests that individual leeches may deploy somewhat different pharmacological cocktails when they feed, whether due to feeding status, physiological condition, or individual genetic variation. For drug discovery, that variability is a double-edged sword: it complicates any attempt to harvest a uniform product from leeches, but it also means that different individuals may naturally emphasize different bioactive molecules, broadening the pool of candidates accessible from this species.

Genetic variation added another layer of insight. The researchers assessed the variant burden across the antithrombotic genes and found it was highest for hyaluronidase_Hman2 and progranulin_Hman, two loci involved in tissue penetration and inflammation-related processes respectively. Intriguingly, the most highly expressed genes were not always the most variable, decoupling expression level from sequence diversity in a way that hints at different evolutionary constraints acting on different parts of the arsenal. Genes that interact directly with host immune and clotting proteins may be under pressure to diversify, while core anticoagulants may be conserved because their function cannot be tinkered with without cost. Selected-gene alignments confirmed localized substitutions and indel-rich regions in informative loci, providing concrete sequence differences that future functional studies can test.

The practical significance of the work lies in its framework as much as its findings. By combining species confirmation, deep genome mapping, transcriptome evidence, and rigorous annotation-based filtering, the study provides a reproducible pipeline for prioritizing antithrombotic candidates not just in H. manillensis but potentially in other blood-feeding species. The 38 transcript-supported candidates, particularly the 30 novel isoforms of known loci, represent a concrete shortlist for laboratory follow-up: expressing these variants recombinantly and testing their activity against thrombin, factor Xa, elastase, and platelet pathways would be the logical next step. The authors are affiliated with the Yunnan International Joint Laboratory with South and Southeast Asia for the Integrated Development of Animal-Derived Anti-Thrombosis Chinese Medicine at Qujing Normal University, and the work was funded under grant number 202503AP140025, reflecting a broader national effort to translate leech-derived molecules into modern antithrombotic therapies.

For a field that began with the discovery of hirudin in medicinal leech saliva more than a century ago, this study marks a shift from serendipity to systematic genomics. Cardiovascular disease remains the leading cause of death worldwide, and existing anticoagulants carry well-known risks of bleeding complications, so the demand for molecules with refined, targeted activity continues to grow. The Asian buffalo leech, now genomically characterized at the sample level with a verified catalog of 72 known genes and 38 new transcript-supported candidates, offers researchers a richer starting point than ever before. What remains is the functional work: confirming what these newly identified isoforms actually do, and whether any of them can be developed into the next generation of clot-preventing drugs. The leech, it turns out, still has plenty to teach modern medicine.

Subject of Research: Whole-genome and transcriptome profiling of antithrombotic genes in the Asian buffalo leech Hirudinaria manillensis

Article Title: Integrated whole-genome and transcriptome profiling of antithrombotic genes in Asian buffalo leech (Hirudinaria manillensis)

Article References: Farooq, T., Naeem, M., Bashir, N. H., Liu, Z., Chen, Y., & Chen, H. (2026). Integrated whole-genome and transcriptome profiling of antithrombotic genes in Asian buffalo leech (Hirudinaria manillensis). BMC Genomics. https://doi.org/10.1186/s12864-026-13354-x

Image Credits: AI Generated

DOI: 10.1186/s12864-026-13354-x

Keywords: Hirudinaria manillensis, antithrombotic genes, leech genomics, whole-genome sequencing, RNA-seq, COI barcoding, hirudin, medicinal leech, transcriptomics, drug discovery, anticoagulants, BMC Genomics

The research team, led by Abdulwasid Abubakari and Charity Ahiabor of Accra Technical University in Ghana, together with George Osei-Adjei and corresponding author Hizbullah Khan of Guangdong Medical University in China, assembled and analyzed 467 high-quality P. aeruginosa genome assemblies. Every genome included in the analysis met a strict quality threshold of greater than 95 percent completeness, and the collection spanned 17 African countries over a remarkable 26-year window, from 1998 to 2024. All of the data were drawn from publicly available assemblies in the NCBI database, meaning the study required no new sampling or ethics approvals, but its systematic reanalysis with standardized tools allowed comparisons that had previously been impossible across such a geographically and temporally scattered dataset.

The methodological backbone of the study was a pangenome analysis, an approach that partitions the collective gene repertoire of a bacterial species into the core genome shared by all isolates and the accessory genome carried by only some. The researchers used iterative pangenome clustering across seven amino acid identity thresholds ranging from 50 to 98 percent, a strategy that captures gene families at different levels of evolutionary relatedness. They complemented this with exploratory core-genome single nucleotide polymorphism phylogenetics to reconstruct the evolutionary relationships among isolates, and with pangenome-wide association studies, or Pan-GWAS, to detect accessory genes whose presence correlates with specific geographic regions.

The results revealed an exceptionally “open” pangenome architecture, quantified by a scaling exponent of gamma equal to 0.23. In pangenome mathematics, a low gamma value signals that each newly sequenced genome is likely to bring a substantial number of previously unseen genes into the catalogue, indicating enormous genetic diversity and extensive gene acquisition. Across the African population, the analysis catalogued 25,501 distinct gene families, a figure that underscores how genetically versatile this pathogen is on the continent. An open pangenome also has practical implications for surveillance: it suggests that local African isolates may carry functional repertoires that global reference-based analyses routinely miss, and that continued sampling will keep yielding new genetic material rather than approaching a saturation point.

Phylogenomic reconstruction resolved 466 unique core genome SNP profiles, and the resulting tree was characterized by strong regional clustering. In other words, isolates from the same part of the continent tended to be more closely related to each other than to isolates from distant regions, a pattern consistent with largely local transmission and evolution rather than a single continental epidemic lineage sweeping across borders. This regional structure matters for public health planning, because it implies that resistance control strategies may need to be tailored to regional population dynamics rather than applied as one uniform continental blueprint.

Perhaps the most striking finding is what the authors describe as an “African Clonal Inversion.” Globally, the sequence type known as ST235 is regarded as the archetypal P. aeruginosa super-clone, a multidrug-resistant lineage that has spread through hospitals worldwide and dominates high-risk clone surveys on other continents. In the African dataset, however, the picture is reversed. Among the 71 isolates identified as belonging to high-risk clones, ST111 accounted for 46.5 percent, or 33 of 71 isolates, outnumbering ST235, which represented 14.1 percent or 10 of 71, by more than three-fold. This inversion suggests that the forces shaping P. aeruginosa success in African healthcare settings differ from those driving the global spread of ST235, and it raises the possibility that ST111 possesses ecological or resistance advantages that are particularly effective in African hospital environments.

The study also mapped the resistome, the complete set of antimicrobial resistance genes carried by the isolates, and found it to be geographically stratified in a way that has direct clinical consequences. The gene blaNDM-1, which encodes the New Delhi metallo-beta-lactamase and confers resistance to some of the most powerful last-line carbapenem antibiotics, dominated in Northern and Eastern African country subsets. Meanwhile, blaVIM-2, another metallo-beta-lactamase gene but from a distinct enzymatic family, was concentrated in Southern Africa. Because both genes threaten the carbapenem class that clinicians rely on when treating severe P. aeruginosa infections, their uneven continental distribution means that empirical treatment guidelines and diagnostic panels may need to account for which resistance determinants are actually circulating in a given region rather than assuming a homogeneous African resistance landscape.

Beyond resistance genes, the Pan-GWAS analysis uncovered regional associations in the accessory genome that hint at how local P. aeruginosa populations differ in their broader functional biology. The researchers found enrichment of a Type VI Secretion System component, the tla3 gene, in North African isolates. The Type VI Secretion System is a molecular spear-like apparatus that bacteria use to inject effector proteins into competing microbes and host cells, playing a major role in interbacterial competition and virulence. In West Africa, the analysis identified enrichment of phenazine biosynthesis clusters, including the phzA2 gene. Phenazines are redox-active secondary metabolites that contribute to P. aeruginosa’s survival, biofilm formation and pathogenicity. These regional functional signatures suggest that different African populations may have evolved distinct ecological strategies, potentially shaped by local hospital conditions, antibiotic prescribing practices, microbial competition and environmental reservoirs.

To make these findings actionable rather than merely archival, the team integrated the entire dataset into a live interactive dashboard hosted on Microreact, a widely used platform for visualizing genomic epidemiology. This means that researchers, public health officials and clinicians across Africa and beyond can explore the phylogenetic trees, geographic distributions and resistance gene patterns themselves, updating the picture as new genomes are deposited. Such open infrastructure is particularly valuable on a continent where surveillance capacity varies widely between countries, because it lowers the technical barrier for local laboratories to place their own isolates in continental and global context.

The study’s authors are careful to frame their work as a foundation rather than a finished map. The dataset, while the largest of its kind for the continent, still reflects the uneven distribution of sequencing capacity and data deposition across Africa, and the researchers themselves categorized some countries with limited sampling as yielding only exploratory Pan-GWAS results. Nevertheless, by demonstrating that African P. aeruginosa populations have their own clonal hierarchy, their own resistance geography and their own accessory gene signatures, the analysis makes a compelling case that the continent can no longer be treated as a footnote in global pathogen genomics. As antimicrobial resistance continues to escalate worldwide, understanding how Priority 1 pathogens evolve in undermapped regions is not just an African concern but a global one, and this study provides both the evidence and the tools to begin that work in earnest.

Subject of Research: Pan-African genomic population structure and antimicrobial resistance distribution of Pseudomonas aeruginosa

Article Title: Pan-African genomics of Pseudomonas aeruginosa highlights regional population structure and AMR stratification

Article References: Abubakari, A., Ahiabor, C., Osei-Adjei, G., & Khan, H. (2026). Pan-African genomics of Pseudomonas aeruginosa highlights regional population structure and AMR stratification. BMC Genomics. https://doi.org/10.1186/s12864-026-13365-8

Image Credits: AI Generated

DOI: 10.1186/s12864-026-13365-8

Keywords: Pseudomonas aeruginosa, pangenomics, antimicrobial resistance, Africa, genomic epidemiology, high-risk clones, resistome, ST111, ST235, Pan-GWAS, MLST, BMC Genomics

The research, led by Mohammad Taghi Zarrinnia together with Karim Hasanpur, Sadegh Alijani, Arash Javanmard, and Majid Olyayee, focused on 153 Marandi hens. Performance and reproductive records were collected across a suite of traits that matter enormously to layer farms: hatch weight, estimated hatch weight, weight and age at first laying, maturity weight, average egg weight, cumulative weight of the first fifty eggs, average weight of those first fifty eggs, egg-laying rate, egg mass, total egg count, feed conversion ratio, and egg weight loss, among others. The team also modeled the trajectory of egg weight over time, a repeated measurement that captures how a hen’s eggs change as her laying cycle progresses.

At the heart of the study lies RAPGEF6, which encodes Rap Guanine Nucleotide Exchange Factor 6, a protein involved in intracellular signaling pathways. The researchers genotyped single nucleotide polymorphisms, or SNPs, within this gene using a well-established laboratory technique called PCR–RFLP, or polymerase chain reaction followed by restriction fragment length polymorphism analysis. This method uses restriction enzymes that cut DNA at specific sequence motifs; when a SNP alters the recognition site, the resulting fragment lengths differ, allowing researchers to infer each bird’s genotype. Genotype frequencies within the population were then determined, providing a picture of genetic variation at the studied markers in the Marandi breed.

To analyze the data, the team deployed two complementary statistical strategies in SAS software. For the remaining growth and performance traits, they used Generalized Linear Models, a flexible framework that relates trait values to genotype categories while accounting for other sources of variation. The repeated egg weight measurements, however, demanded a more sophisticated treatment, because successive records from the same hen are not statistically independent — a hen that lays heavier-than-average eggs in one period tends to do so in the next. Ignoring this correlation inflates false confidence, so the researchers fitted eight different covariance structures to model how measurements within each bird relate to one another over time.

Those eight candidate structures spanned the standard toolkit of repeated-measures modeling: variance components, compound symmetry and its heterogeneous variant, first-order autoregressive and heterogeneous autoregressive structures, Toeplitz forms including banded and heterogeneous versions, a non-diagonal factor analytic structure, an unstructured model, and the autoregressive moving average model known as ARMA(1,1). Using SAS’s PROC MIXED procedure, the team compared how well each structure accounted for the repeated egg weight records. The ARMA(1,1) model emerged as the best performer, combining an autoregressive component, in which each measurement depends on the previous one, with a moving average component that captures short-lived shocks. Choosing the optimal covariance structure matters because it sharpens the precision of every downstream estimate, allowing a more consistent and reliable assessment of egg production patterns in the flock.

The headline association came from the Generalized Linear Model analysis: one SNP in RAPGEF6 showed a significant relationship with maturity weight, the body weight a hen reaches at maturity, with a p-value of 0.014. For the other sixteen traits examined, no significant associations emerged, with all p-values exceeding 0.05. On its face, a significant hit on a growth trait in a signaling gene might seem like a promising lead for breeders hoping to use DNA markers to accelerate genetic improvement in indigenous poultry.

But the researchers did not stop there, and this is where the study distinguishes itself. They conducted a post-hoc statistical power analysis using the GLMPOWER procedure in SAS, asking a blunt question: given the sample size of 153 hens, how likely was the study to detect an effect of the size actually observed? The answer was sobering. The analysis revealed only 5 percent power to detect the maturity weight effect — meaning that if the true association were exactly as strong as the one observed, the study would have missed it 95 percent of the time. A finding detected under such conditions is statistically unstable, and the authors explicitly label it as such.

This transparency carries lessons far beyond one Iranian chicken breed. In genetic association studies, underpowered designs are notorious for producing findings that fail to replicate, because a significant result in a low-power study is disproportionately likely to reflect an inflated effect estimate or simple chance rather than a genuine biological signal. The Tabriz team confronted this reality directly in their conclusions, stating that given the single-SNP design, the small sample size, and the severely limited statistical power, all findings should be interpreted as exploratory and preliminary. They further emphasized that no recommendations for marker-assisted selection — the breeding strategy of using DNA markers to guide selection decisions — can be made on the basis of the current data, and that validation in independent, well-powered populations is required before any practical application.

Why does this matter for the future of the Marandi chicken and indigenous livestock genetics more broadly? Indigenous breeds are reservoirs of locally adapted genetic variation, often carrying alleles suited to harsh environments, modest feed inputs, and disease pressure that commercial lines lack. Yet their genetic potential remains underexplored, and studies like this one represent early reconnaissance missions into that variation. Egg production sustainability is economically critical for layer farms, and identifying genes that influence growth and laying traits could eventually inform breeding programs that improve productivity while preserving the breed’s distinctive adaptive identity. RAPGEF6 polymorphisms, the authors suggest, may influence growth traits in Marandi chickens, offering preliminary insights into the genetic underpinnings of performance — insights that future, larger studies can now test rigorously.

The study also demonstrates the quiet power of careful statistical modeling in agricultural genomics. The finding that the ARMA(1,1) covariance structure best captured the repeated egg weight measurements is a practical contribution in its own right: it shows that for laying hens measured repeatedly across a production cycle, accounting for both the carryover between consecutive records and short-term fluctuations yields more precise estimates than simpler structures like compound symmetry. For researchers analyzing longitudinal production data in poultry and other livestock, such methodological guidance can meaningfully improve the reliability of genetic evaluations. Taken together, the study is a model of scientific candor — a report that celebrates a modest lead, quantifies exactly how fragile that lead is, and lays out a disciplined statistical framework for the larger, better-powered investigations that indigenous breed genomics will need next.

Subject of Research: Association between RAPGEF6 gene polymorphisms and egg production and growth traits in Marandi chickens

Article Title: Genetic marker analysis of RAPGEF6 gene and covariance structure modelling for egg

Article References: Genetic marker analysis of RAPGEF6 gene and covariance structure modelling for egg. (n.d.). https://doi.org/10.1186/s12864-026-13373-8

Image Credits: AI Generated

DOI: 10.1186/s12864-026-13373-8

Keywords: RAPGEF6, Marandi chicken, SNP, genetic association study, egg production, covariance structure, ARMA(1,1), PCR-RFLP, statistical power, marker-assisted selection, indigenous poultry, BMC Genomics

Coelomocytes are the cellular workhorses of the echinoderm immune system. They circulate in the coelomic fluids that fill the body cavity, performing functions that in vertebrates are distributed among blood cells and immune cells of many kinds. For decades, biologists have classified these cells by looking at their shapes under the microscope, and sea cucumbers, or holothuroids, stand out among echinoderms for displaying the greatest variety of coelomocyte morphotypes of any class in the phylum. Yet morphology alone has proved a poor guide to function. While the overall immune role of these cells is broadly accepted, the specific jobs performed by each cell type have remained murky, and molecular data specific to the different morphotypes have been scarce in the literature.

Single-cell RNA sequencing, or scRNA-seq, has transformed how biologists untangle cellular heterogeneity, allowing researchers to profile the gene expression of thousands of individual cells and group them into transcriptionally distinct populations. The technique, however, has been largely underutilised in studies of non-model organisms, particularly among invertebrates. The new study set out to close that gap for sea cucumbers by applying the method to coelomocytes harvested from the perivisceral fluid of Holothuria forskali, a species collected with the help of the collection service of the Roscoff Biological Station in France. Sequencing was carried out at the Genomic Platform of the GIGA Institute in Liège.

The analysis identified ten distinct clusters of cells, each assumed to correspond to a separate transcriptional coelomocyte population. The clustering revealed a striking architecture: one cluster, designated cluster 0, occupied a central position relative to all the others on the two-dimensional map of gene expression, suggesting that it may represent undifferentiated cells, a reservoir of progenitors from which more specialised populations arise. At the opposite extreme, cluster 6 sat markedly apart from every other group, hinting at a fundamentally different biological identity and function.

To work out what each cluster might actually do, the team performed functional enrichment analyses and hunted for immune marker genes among the transcripts enriched in each population. The results showed that several clusters carry out key immune functions, including the recognition of pathogens, phagocytosis, the engulfing and destruction of foreign particles, the activation of complement-like pathways, and the regulation of redox balance, the delicate chemistry of oxidation and reduction that cells must manage to survive an immune response. These findings provide the first tentative clues about which transcriptional populations correspond to which defensive duties, offering a molecular handle on a system that has until now been described mostly by cell shape.

One of the most intriguing results concerns a recently discovered cell type called the carotenocyte. These cells, rich in carotenoids, the pigments responsible for many of the reds, oranges and yellows in nature, were only recently described in sea cucumbers, and their presence in the perivisceral fluid was confirmed when the researchers examined their processed samples under the microscope. By drawing on transcriptomic data previously generated for this cell type using bulk RNA sequencing, the team was able to confidently assign cluster 6, the most divergent group on the map, to the carotenocyte lineage, and to deepen the picture of what these pigment-laden cells express.

The convergence of two independent datasets proved decisive. Marker genes supported by both the bulk RNA sequencing of carotenocyte-enriched samples and the new single-cell data gave the researchers confidence in the identity of cluster 6, an approach that illustrates how single-cell and bulk methods can be combined to anchor cell identities in non-model organisms where no reference atlas exists. The study’s supplementary tables catalogue the marker genes for each cluster, their functional annotations against multiple databases, and the results of enrichment analyses against KEGG pathways and gene ontology categories, providing a resource for other researchers working on echinoderm immunity.

Rigorous quality control underpinned the analysis. The authors filtered cells on the basis of the number of unique molecular identifiers and detected genes per cell, examined the expression of mitochondrial genes to flag stressed or damaged cells, and used the DoubletFinder tool to detect and remove doublets, artefacts in which two cells are captured and sequenced as one. Sensitivity analyses testing different filtering parameters showed that while individual clusters could merge or split depending on the settings, the overall configuration of the map remained stable, with cluster 0 staying central and cluster 6 remaining the most divergent population, a reassuring sign that the biological signal is robust.

The implications reach well beyond sea cucumbers. Echinoderms are deuterostomes, the same major lineage as chordates, which makes their immune cells valuable comparators for understanding how immune cell lineages evolved on the branch of the tree that eventually produced our own adaptive immune system. A better grasp of holothuroid coelomocyte diversity could also help interpret stress responses in these animals, which are ecologically important grazers and the target of a rapidly expanding aquaculture industry in Asia. Knowing which cells respond to pathogens or environmental stressors, and through which genes, could inform disease management in farmed populations and shed light on how wild sea cucumbers cope with a changing ocean.

The authors are careful to frame the study as a pioneering first step. Linking the remaining transcriptional clusters to the morphotypes described in the historical literature will require further analyses, and the functional assignments drawn from enrichment analyses remain hypotheses to be tested experimentally. Still, the study demonstrates that single-cell transcriptomics can be successfully deployed on a non-model invertebrate with no reference-grade tools, and it delivers a preliminary but detailed map of the functional diversity of holothuroid coelomocytes. For a phylum whose immune system has been studied mostly through a microscope lens, the jump to single-cell genomics marks the beginning of a new era, one in which the hidden division of labour among sea cucumber immune cells can finally be read directly from their genes.

Subject of Research: Single-cell transcriptomic profiling of coelomocyte immune cell diversity in the sea cucumber Holothuria forskali

Article Title: Single-cell transcriptomics reveals transcriptional diversity of sea cucumber perivisceral fluid coelomocytes

Article References: Wambreuse, N., Lavergne, A., Fievez, L., Bureau, F., Zhang, L., Deng, B., Caulier, G., Eeckhaut, I., & Delroisse, J. (2026). Single-cell transcriptomics reveals transcriptional diversity of sea cucumber perivisceral fluid coelomocytes. BMC Genomics. https://doi.org/10.1186/s12864-026-13356-9

Image Credits: AI Generated

DOI: 10.1186/s12864-026-13356-9

Keywords: sea cucumber, coelomocytes, single-cell RNA sequencing, echinoderm immunity, carotenocytes, Holothuria forskali, transcriptomics, immune cells, deuterostomes, phagocytosis, BMC Genomics, non-model organisms

Chemosensation is the invisible language of insect life. From locating a host plant across a forest clearing to recognizing a potential mate by pheromone, nearly every critical decision an insect makes depends on its ability to detect, bind, and transduce chemical signals. These functions are carried out by families of proteins expressed in sensory organs, chiefly the antennae. Odorant-binding proteins, or OBPs, and chemosensory proteins, or CSPs, capture volatile molecules in the fluid-filled pores of sensory hairs and shuttle them to receptors embedded in nerve cell membranes. Odorant receptors, or ORs, ionotropic receptors, or IRs, and gustatory receptors, or GRs, then convert chemical binding events into electrical signals the brain can interpret. Sensory neuron membrane proteins, or SNMPs, assist in the detection of fatty acid-derived compounds such as pheromones, while Niemann-Pick type C2 proteins, or NPC2s, contribute to the transport of hydrophobic molecules. Until this study, not a single one of these gene families had been reported for R. cyatheae, leaving researchers with no molecular handle on how the wasp tracks its unusual fern host.

To build that handle, the team assembled a transcriptome database from male and female adult wasps, sequencing the complete set of messenger RNA molecules being expressed in their bodies. The effort yielded 30,296 unigenes, the distinct assembled sequences representing the insect’s active genes, with an N50 length of 3,286 base pairs, a measure indicating that half of the assembled sequences were longer than that threshold and reflecting a high-quality assembly. The researchers then compared these sequences against six major public annotation databases: NR, Swiss-Prot, Pfam, eggNOG, GO, and KEGG. In total, 11,109 unigenes, or 36.67 percent of the full set, could be functionally annotated. The NR database proved the most productive, matching 10,774 unigenes, while the KEGG pathway database annotated the fewest, 6,300. This annotation step is essential because it allows researchers to distinguish genuine chemosensory genes from the thousands of other sequences involved in basic cellular housekeeping.

Screening the annotated transcriptome against known insect chemosensory gene families revealed 90 candidate chemosensory genes in R. cyatheae. The inventory includes 11 OBPs, 10 CSPs, 6 NPC2s, 24 ORs, 20 IRs, 15 GRs, and 4 SNMPs. Within the OR family, the researchers distinguished 23 typical odorant receptor genes and a single Orco gene, the highly conserved co-receptor that partners with typical ORs across virtually all insect species and is essential for odor detection. The numbers themselves tell an evolutionary story. Compared with insects that rely heavily on chemical communication, such as moths with dozens or hundreds of ORs, the wasp’s relatively compact receptor repertoire may reflect its narrow ecological niche, though the authors present the counts as a baseline for future comparison rather than drawing firm ecological conclusions.

To place these genes in evolutionary context, the team constructed phylogenetic trees relating the chemosensory genes of R. cyatheae to their counterparts in other insect species. Phylogenetic analysis groups genes by sequence similarity, allowing researchers to infer which genes are orthologs, descended from a common ancestor, and which may have expanded or contracted in particular lineages. Such trees also help flag genes that cluster with well-characterized relatives in other insects, providing hypotheses about function that can later be tested experimentally. For a species whose chemosensory biology had never been examined at the molecular level, this comparative framework is the first map connecting the wasp’s sensory genes to the broader landscape of insect olfaction research.

One of the most intriguing findings emerged when the researchers compared gene expression between male and female adults. Thirteen of the candidate chemosensory genes were differentially expressed between the sexes, meaning their activity levels differed significantly. Sex-biased expression in chemosensory genes often points to roles in behaviors that differ between males and females, such as mate location through pheromone detection in males or host plant selection for oviposition in females. Because R. cyatheae infestations coincide with the sprouting period of its host fern, understanding which sex drives host finding, and which molecular sensors are involved, could reveal weak points in the infestation cycle that might be targeted by behavioral interventions or attractant-based monitoring tools.

To confirm that the transcriptome data accurately reflected biology rather than sequencing artifacts, the team verified the tissue expression profiles of the differentially expressed genes using reverse transcription quantitative polymerase chain reaction, or RT-qPCR. This laboratory technique measures the abundance of specific messenger RNA molecules with high sensitivity and is the standard method for validating transcriptome-wide expression patterns. The RT-qPCR results corroborated the expression differences observed in the sequencing data, strengthening confidence that the 13 sex-biased genes are genuinely regulated differently in males and females. This kind of validation is a critical quality checkpoint in genomics studies, where computational identification of candidate genes must be backed by independent experimental evidence before functional work can proceed.

The significance of the study extends beyond a single pest species. Tree ferns of the genus Alsophila are living fossils, and A. spinulosa is protected within the Chishui Alsophila National Nature Reserve, whose administration supported the fieldwork. Conservation of an endangered fern depends on managing the insects that attack it, and effective management in turn depends on understanding the chemical ecology of the pest. The gene catalogue now makes it possible to pursue questions that were previously unanswerable: which odorant receptors respond to fern volatiles, whether males detect female-produced pheromones, and whether synthetic lures could be designed to monitor or disrupt mating. Similar transcriptome-based approaches have enabled the development of semiochemical-based control tools for other forest pests, and the new data provide the raw material for analogous efforts against R. cyatheae.

The work also fills a taxonomic gap. Hymenoptera, the order that includes wasps, bees, and ants, is enormously diverse, yet chemosensory gene repertoires have been characterized for only a fraction of its members, and herbivorous hymenopterans that attack ferns are especially poorly studied. By documenting OBPs, CSPs, NPC2s, ORs, IRs, GRs, and SNMPs in a fern-feeding sawfly relative, the study offers comparative data that will help researchers trace how chemosensory gene families evolve as insects shift onto unusual host plants. The authors note that their findings lay a solid molecular foundation for future investigations of gene function and the mechanisms of olfactory perception in R. cyatheae.

From a technical standpoint, the study demonstrates a now-standard but powerful pipeline for working with non-model insects: assemble a transcriptome, annotate it against multiple databases, mine it for gene families of interest, place the candidates in phylogenetic context, and validate expression patterns experimentally. The resulting catalogue of 90 chemosensory genes, anchored by a high-quality assembly of 30,296 unigenes and confirmed sex-specific expression patterns, transforms R. cyatheae from a molecularly unknown pest into a tractable research system. For the endangered tree fern it threatens, that transformation may ultimately matter most: every gene identified is a potential target for understanding, predicting, and eventually managing the outbreaks that strike when the fern’s tender new fronds emerge. The research was funded by a project on pest resistance mechanisms of different Cyathea species based on multi-generation transcriptomes and protein metabolomes, and by a Science and Technology Innovation Talent Team project of Guizhou Province, reflecting sustained regional investment in protecting this unique fragment of prehistoric flora.

Subject of Research: Identification of chemosensory genes in the adult transcriptome of the tree fern pest wasp Rhoptroceros cyatheae

Article Title: Analysis and identification of chemosensory genes in the transcriptome of adult Rhoptroceros cyatheae

Article References: Analysis and identification of chemosensory genes in the transcriptome of adult Rhoptroceros cyatheae. (n.d.). https://doi.org/10.1186/s12864-026-13393-4

Image Credits: AI Generated

DOI: 10.1186/s12864-026-13393-4

Keywords: Rhoptroceros cyatheae, chemosensory genes, transcriptome, odorant-binding proteins, odorant receptors, Alsophila spinulosa, tree fern, olfaction, RT-qPCR, phylogenetic analysis, BMC Genomics, entomology

A new study published in BMC Genomics by Wei Fan, Hongli Liu, Jian Su and colleagues, working across the Neijiang Academy of Agricultural Sciences, Sichuan Agricultural University, Neijiang Normal University and an industry partner, set out to explain that growth disparity at the molecular level. Rather than focusing on a single organ, the team took a multi-tissue approach, sequencing the transcriptomes of the brain, liver and skeletal muscle of leucistic northern snakeheads whose growth rates differed significantly. The logic behind this design is straightforward but powerful: growth is not the product of one tissue acting alone. It emerges from a conversation between the brain, which coordinates appetite and endocrine signals; the liver, which processes nutrients and secretes growth-regulating factors; and the muscle, where the actual deposition of new protein takes place.

Transcriptome sequencing, the core technology of the study, captures a snapshot of which genes are actively being transcribed into messenger RNA in each tissue at the moment of sampling. By comparing these expression profiles between fast-growing and slow-growing leucistic individuals, researchers can identify genes and signaling pathways whose activity correlates with rapid growth. High-throughput sequencing platforms now make it possible to quantify the activity of virtually every expressed gene simultaneously, and sophisticated bioinformatics pipelines then sort the resulting data, flagging genes whose expression differs significantly between groups and mapping them onto known biological pathways.

The analysis revealed a series of differential pathways closely associated with growth rate. Among the most notable was the TGF-beta signaling pathway, which the authors found to be activated in the larger, faster-growing group. This pathway family is a central regulator of tissue development across the animal kingdom, and in fish it has well-documented roles in myogenesis, the formation and maturation of muscle fibers. Members of the TGF-beta superfamily, including myostatin, act as molecular brakes on muscle growth, while other related signals promote proliferation and differentiation of muscle precursor cells. An imbalance in this finely tuned system can therefore translate directly into differences in how much muscle an animal builds from the same amount of feed.

According to the study’s functional analysis, the differential pathways identified in the fast-growing fish appear to influence key biological processes such as muscle growth and metabolism. In practical terms, this suggests that the growth advantage of certain leucistic snakeheads is not simply a matter of eating more, but of running an internally different program: one in which the signals governing muscle fiber recruitment, protein turnover and energy allocation are tuned toward more efficient somatic development. The liver’s role in this program is particularly important, because it integrates dietary nutrients, synthesizes the insulin-like growth factors that drive muscle anabolism, and manages lipid and glucose metabolism that fuels tissue construction.

The northern snakehead itself is a biologically remarkable animal, which makes these findings all the more interesting. Native to East Asian freshwater systems, Channa argus is an obligate air-breather equipped with a suprabranchial organ that allows it to survive in oxygen-poor waters and even tolerate brief periods out of water. It is famous for its hardiness, its voracious appetite and its rapid growth, traits that have made it a valuable aquaculture species in China even as its ecological vigor has raised alarm in regions where it has been introduced elsewhere. A leucistic morph of such a robust species, combining unusual appearance with the species’ inherent growth potential, represents exactly the kind of variant that fish farmers dream of, provided the genetics of its performance can be understood and managed.

Leucism itself is worth distinguishing from albinism, since the two conditions are often confused. Albino animals lack melanin entirely because of defects in melanin synthesis, and they typically have red or pink eyes. Leucistic animals, by contrast, suffer a partial loss of pigmentation across all pigment cell types, producing pale or white coloration while often leaving the eyes normally pigmented. In fish, the genetic and cellular bases of leucism are diverse, and pigmentation anomalies can sometimes be linked to broader developmental differences. Understanding whether the leucistic morph’s growth traits are independent of, or connected to, its pigmentation genetics remains an open question, but the current study’s focus on growth pathways provides a crucial piece of the puzzle.

The multi-tissue design of the research deserves particular emphasis, because single-organ studies can miss the systemic nature of growth regulation. The brain, as the command center, integrates environmental cues and hormonal feedback to modulate feeding behavior and the secretion of pituitary growth hormone. The liver responds to that growth hormone by producing insulin-like growth factors, and it also channels dietary energy toward either storage or construction of body tissue. The muscle, finally, is where those signals are executed, with satellite cells and myofibrillar protein synthesis determining how much flesh is actually added. By sampling all three tissues from the same individuals, the researchers could look for coordinated shifts in gene expression across this axis, a far more informative picture than any single tissue could provide.

The implications for molecular breeding are significant. Traditional selective breeding of farmed fish relies on visible performance traits, measuring which individuals grow fastest and breeding from them, a slow process that captures only part of the underlying genetic variation. Transcriptomic markers, by contrast, can reveal which molecular pathways distinguish high-performing animals, potentially allowing breeders to select for favorable expression profiles or underlying genetic variants at a much earlier stage. For the leucistic northern snakehead, whose market value depends on both its unusual coloration and its meat quality, being able to identify fast-growing individuals early, or to breed lines in which the growth-promoting pathways are consistently active, could substantially improve the economics of farming this niche product.

The study also contributes more broadly to aquaculture science, a field under growing pressure to produce protein efficiently as wild fisheries stagnate and demand for farmed fish rises. Growth traits are among the most economically important characteristics in virtually every farmed species, and the signaling pathways identified here, particularly TGF-beta signaling, are conserved across fish lineages. Insights gained from the snakehead may therefore inform research on other cultured species, from carp to tilapia, where manipulating the same regulatory networks could enhance growth performance. The authors note that their work lays a theoretical foundation for the molecular breeding of the leucistic northern snakehead and for the sustainable development of the aquaculture industry, framing the research not merely as a description of one unusual fish but as a step toward a more precise, biology-driven approach to fish farming.

There remain, of course, important caveats and next steps. Transcriptomics measures gene expression, not the final functional output of proteins and metabolism, and the causal role of any single pathway must be confirmed through further experiments, whether functional validation of candidate genes, proteomic follow-up or breeding trials that test whether the expression signatures predict real-world growth. The study was conducted under institutional animal care approval at Sichuan Agricultural University and reported in accordance with the ARRIVE guidelines, reflecting the growing rigor expected of animal research. Supported by Sichuan provincial science and technology programs, the work exemplifies the collaboration between academic institutions and industry that increasingly characterizes modern aquaculture genomics. For now, the pale ghost of the snakehead world has offered up its secrets a little further: beneath its unusual skin, the leucistic northern snakehead carries a molecular growth program whose dynamics researchers are only beginning to read, and whose full exploitation could reshape how this remarkable species is farmed.

Subject of Research: Molecular mechanisms underlying growth rate differences in the leucistic morph of the northern snakehead, Channa argus

Article Title: Multi-tissue transcriptomics elucidate growth disparity mechanisms in the leucistic morph of Channa argus

Article References: Fan, W., Liu, H., He, Y., Zhuo, T., Wang, J., Geng, Y., Huang, X., Chen, H., Jiao, X., Wang, Q., Yang, H., Deng, Y., & Su, J. (2026). Multi-tissue transcriptomics elucidate growth disparity mechanisms in the leucistic morph of Channa argus. BMC Genomics. https://doi.org/10.1186/s12864-026-13392-5

Image Credits: AI Generated

DOI: 10.1186/s12864-026-13392-5

Keywords: Channa argus, leucism, transcriptomics, growth rate, TGF-beta signaling, aquaculture, molecular breeding, muscle growth, gene expression, northern snakehead, BMC Genomics, fish genetics

The research, led by first authors Leila Namvar, Kamilla Sedov and Madison James Yang under the corresponding authorship of Birgitt Schüle, examined 73 low-passage induced pluripotent stem cell clones generated from 25 parental lines. The parental cells came from two distinct sources: peripheral blood mononuclear cells, the mixed population of white blood cells that can be drawn from a simple blood sample, and fibroblasts, the connective-tissue cells typically harvested from a skin biopsy. Both cell types are widely used starting materials for making induced pluripotent stem cells, and until now there has been little systematic evidence about whether the choice between them leaves a genomic fingerprint on the resulting stem cell lines.

Optical genome mapping works in a fundamentally different way from the sequencing-based methods that dominate genomics today. Instead of breaking DNA into millions of short fragments and reassembling them computationally, the technique preserves ultra-high-molecular-weight DNA molecules, each hundreds of thousands of base pairs long. These molecules are fluorescently labeled at specific sequence motifs, creating a barcode pattern along each DNA strand. The labeled molecules are then run through a nanochannel array that stretches them out in a linear fashion, allowing cameras to image the barcode patterns directly. By aligning these observed barcodes against a reference genome, the software can detect where the actual genome deviates from the reference: deletions, duplications, inversions, insertions and translocations all leave characteristic signatures in the barcode pattern.

The practical advantage of this approach lies in its resolution and sensitivity for large, complex events. Conventional karyotyping, the century-old technique of staining and photographing chromosomes under a microscope, can only reliably detect rearrangements on the order of several million base pairs. Copy-number microarrays detect dosage changes but miss balanced events such as inversions and reciprocal translocations that preserve the total amount of DNA. Short-read whole-genome sequencing struggles with repetitive regions and with rearrangements that span large distances. Optical genome mapping, by contrast, can detect structural variants down to approximately two kilobase pairs at sufficient coverage, and because each variant call is supported by evidence at the level of individual DNA molecules, the method can also identify mosaic events present in only a fraction of the cells in a culture, a common situation in stem cell lines where a single aberrant clone can expand during passaging.

When the Stanford team applied this technology to their panel of clones, the differences between the two starting cell types were striking. Fibroblast-derived induced pluripotent stem cells showed a higher overall burden of structural variants, and this burden was enriched for duplications of at least 100 kilobase pairs. These duplications were also more likely to overlap protein-coding genes, meaning they had the potential to alter gene dosage in ways that could affect cellular function. In addition, the structural variants found in fibroblast-derived clones overlapped more frequently with fragile sites, specific regions of the genome that are prone to breakage during replication, and with recurrent chromosomal hotspot regions that have been repeatedly associated with culture-induced rearrangements in pluripotent stem cells.

Peripheral blood mononuclear cell-derived clones told a different story. These lines showed fewer structural variants overall, and a higher proportion of the clones carried no detectable clone-specific structural variants at all. In other words, when the starting material was blood, the reprogramming process more often produced genomically clean stem cell lines, at least by the standards of what optical genome mapping can see. The authors are careful to frame this as an association rather than a proven causal mechanism, but the pattern across 73 clones from 25 parental lines is difficult to dismiss as chance.

Why might fibroblasts be more prone to acquiring structural variants during reprogramming? The study does not definitively answer this question, but several biological considerations are consistent with the findings. Fibroblasts are adherent cells that typically require more extensive expansion in culture before and during reprogramming, and prolonged cell division is a well-known driver of copy-number changes, particularly duplications that arise through replication errors and DNA breakage at fragile sites. Blood cells, by contrast, are naturally short-lived and are isolated from a sample that can be processed relatively quickly. It is also possible that the epigenetic and replicative state of the two cell types at the moment of reprogramming influences how faithfully the genome is copied and reorganized as the cells pass through the pluripotent state. The enrichment of variants at fragile sites and known culture-associated hotspots in the fibroblast-derived lines points toward replication stress and cell culture itself as contributing factors.

The clinical and research implications of this work are substantial. Induced pluripotent stem cells are used to model diseases ranging from Parkinson’s disease to Alzheimer’s disease, and the Stanford study drew its donor samples from the Stanford Alzheimer’s Disease Research Center, with funding from the National Institutes of Health, the Michael J. Fox Foundation for Parkinson’s Research and the California Institute for Regenerative Medicine. In disease modeling, a structural variant that disrupts a gene could be mistaken for a disease-relevant finding or could mask the true cellular phenotype of interest. In cell therapy, where reprogrammed cells may eventually be transplanted into patients, a duplication affecting an oncogene or a deletion affecting a tumor suppressor could raise safety concerns. The finding that structural variants frequently overlapped protein-coding genes, fragile sites and recurrent chromosomal hotspots underscores why clone-level genomic assessment matters before any downstream application.

What makes the study particularly actionable is its argument for optical genome mapping as a routine quality-control step during stem cell line generation and selection. The authors position the technology as a post-reprogramming screening tool capable of catching clone-specific events below the resolution of conventional cytogenetic and SNP-array-based assays. Because the method detects variants at approximately two kilobase pairs and provides molecule-level evidence for mosaic events, it occupies a sweet spot between the coarse resolution of karyotyping and the computational complexity of long-read sequencing. For laboratories producing banks of induced pluripotent stem cell lines, the practical message is that screening each clone with a genome-wide structural variant assay before committing it to experiments or therapies could prevent costly surprises downstream.

The study also carries a subtler message about experimental design. When researchers plan a new panel of induced pluripotent stem cell lines, the choice between blood and skin as the source tissue is often made for practical reasons, such as patient accessibility or biopsy logistics. This work suggests that the choice may have genomic consequences that persist into the pluripotent state even at early passages, before extended culture has had time to accumulate additional changes. Fibroblast-derived lines may require closer scrutiny, particularly for large duplications affecting genes, while blood-derived lines may more often yield genomically clean clones. As induced pluripotent stem cell technology moves from the research bench toward clinical applications, understanding how the origin of a stem cell line shapes its genome is becoming an essential part of ensuring that the cells scientists study and eventually transplant are as faithful to the original genome as possible. Optical genome mapping, this study suggests, offers a practical and high-resolution way to make that assessment.

Subject of Research: Structural variant differences between blood- and fibroblast-derived human induced pluripotent stem cells detected by optical genome mapping

Article Title: Optical genome mapping identifies source-associated structural variant differences across early-passage human iPSCs

Article References: Namvar, L., Sedov, K., Yang, M. J., Hermosillo, R., Zafar, F., & Schüle, B. (2026). Optical genome mapping identifies source-associated structural variant differences across early-passage human iPSCs. BMC Genomics. https://doi.org/10.1186/s12864-026-13395-2

Image Credits: AI Generated

DOI: 10.1186/s12864-026-13395-2

Keywords: induced pluripotent stem cells, optical genome mapping, structural variants, fibroblasts, peripheral blood mononuclear cells, nuclear reprogramming, genomic quality control, copy number variation, fragile sites, chromosomal hotspots, disease modeling, BMC Genomics

The team, led by Jun Gao and Lingwei Sun of the Institute of Animal Science and Veterinary Medicine at the Shanghai Academy of Agricultural Sciences, with corresponding authors Caifeng Wu and Jianjun Dai, built the assembly from scratch using PacBio third-generation long-read sequencing combined with Hi-C technology. Long reads capture stretches of DNA far longer than older short-read methods, which is essential for resolving repetitive regions, structural variants, and other difficult stretches of the genome. Hi-C then provides the crucial scaffolding step: by chemically fixing chromosomes inside the cell nucleus and sequencing the physically adjacent DNA fragments, Hi-C data reveal which sequences belong to the same chromosome and how they are ordered, allowing the contigs to be anchored into complete chromosome-scale assemblies.

After a de-redundancy step to remove duplicated sequences, the assembled Pudong White Pig genome totaled approximately 2.68 gigabases, in line with the expected size of the pig genome. The quality metrics that matter most to genome builders told a favorable story: the new assembly showed higher sequence integrity and fewer genomic gaps than Sscrofa11.1, the Duroc-derived reference genome that is currently the most widely used standard in pig genomics. That comparison is significant because reference genomes shape nearly every downstream analysis, from mapping reads in association studies to calling variants; gaps and missing segments in a reference can hide biologically important variation, particularly for breeds that diverged from the reference animal’s lineage long ago.

A reference genome for a single breed, however, only tells part of the story. To capture the fuller spectrum of pig genetic diversity, the researchers constructed a pan-genome from 30 domestic and international pig breeds. A pan-genome represents the union of all genomic content across a set of individuals, distinguishing the core sequences shared by all breeds from the accessory sequences present in only some. In this analysis, the team identified a total of 135,859 non-reference sequences, or NRSs, stretches of DNA that simply do not exist in the standard reference genome. These absent sequences can include functional genes, regulatory elements, and structural variants that influence traits ranging from disease resistance to reproduction, which is precisely why they are invisible to studies that rely on a single reference.

Stringent quality control filtering validated 65,193 of those non-reference sequences, with a combined length of 108.33 megabases, roughly 4 percent of the pig genome’s total size. That is a substantial amount of genetic material that the standard reference fails to capture. Within this validated set, the researchers identified 594 non-reference sequences specific to Shanghai’s indigenous pig breeds, found across the four Shanghai native pig populations included in the pan-genome analysis. These breed-specific insertions represent a molecular signature of local genetic heritage, shaped by centuries of adaptation to the climate, feeding practices, and selection pressures of the Yangtze River delta region.

Some of these sequences land in scientifically intriguing places. One insertion site, specific to Shanghai indigenous pigs, sits at the BMPR1B gene, which is linked to reproductive performance in Taihu lake pigs, a famous group of Chinese native breeds renowned for their exceptional fecundity. The presence of a breed-specific insertion at a locus associated with such an economically vital trait suggests that local breeds may harbor unique regulatory or coding variation affecting litter size, variation that commercial reference genomes cannot reveal. For breeders working to conserve and improve native stocks, pinpointing such loci is a first step toward understanding and potentially exploiting the genetic basis of traits that make indigenous pigs valuable.

The analysis also revealed a striking pattern in the distribution of non-reference sequences: a hotspot region concentrated in the central portion of chromosome X, spanning roughly the 59 to 89 megabase interval. Hotspots of this kind can indicate regions of unusual structural dynamism, perhaps shaped by selection, recombination patterns, or the history of the X chromosome in domestication. Because the X chromosome plays a distinctive role in inheritance and can carry loci affecting fertility and other traits, a concentrated reservoir of reference-missing sequence there is a finding that will likely prompt closer examination in future studies of pig genetics and breeding.

Beyond the pan-genome work, the assembly allowed the team to confirm a previously reported splice mutation at position 41,486,012 on chromosome 8 of the Sscrofa11.1 reference genome. This locus is central to the genetics of coat color: at this site, both the Pudong White Pig and the Shanghai White breed carry the same dominant white genotype characteristics seen in the Landrace and Yorkshire breeds, the two dominant commercial white pig lines worldwide. The confirmation ties the Pudong White Pig’s defining visual trait to a known molecular mechanism and demonstrates that the new assembly can be used to validate and contextualize variants originally discovered on the older reference, an important cross-check for the genomics community.

The broader significance of the work lies in what it provides for conservation and breeding. The Pudong White Pig is endangered, and endangered livestock breeds face a double risk: the loss of the animals themselves and the loss of the genetic adaptations they carry, adaptations that may prove valuable as agriculture confronts disease pressure, climate change, and shifting market demands. A chromosome-level reference genome is a permanent, digital record of that heritage. It enables accurate genotyping of remaining herds, supports marker-assisted selection and genomic selection programs aimed at rebuilding population numbers without losing genetic diversity, and allows researchers to detect inbreeding and manage breeding pairs to maximize retained variation.

The study also adds to a growing recognition that global livestock genomics needs more than one reference per species. As this work shows, a single Duroc-based reference leaves out more than 100 megabases of real, validated sequence found across just 30 breeds, including hundreds of sequences unique to Shanghai’s native pigs. The Pudong White Pig assembly, funded by China’s National Key R&D Program and the Shanghai Agricultural Science and Technology Innovation Project, and conducted with ethics approval from the Shanghai Academy of Agricultural Sciences, offers a template for similar efforts: combine long reads with Hi-C, anchor every chromosome, then embed the result in a multi-breed pan-genome to reveal what single references miss. For a white-coated survivor of Shanghai’s agricultural landscape, the new genome is both a scientific milestone and a lifeline.

Subject of Research: Chromosome-level de novo genome assembly and pan-genome analysis of the endangered Shanghai Pudong White Pig

Article Title: Chromosome-level genome assembly of the Shanghai Pudong White Pig (Sus scrofa domesticus‌)

Article References: Gao, J., Sun, L., Cao, M., Tu, W., Zhang, H., Zhang, S., Xu, J., He, M., Zhang, D., Wu, C., & Dai, J. (2026). Chromosome-level genome assembly of the Shanghai Pudong White Pig (Sus scrofa domesticus‌). BMC Genomics. https://doi.org/10.1186/s12864-026-13407-1

Image Credits: AI Generated

DOI: 10.1186/s12864-026-13407-1

Keywords: Pudong White Pig, genome assembly, PacBio sequencing, Hi-C, pan-genome, non-reference sequences, Sus scrofa domesticus, endangered livestock breed, BMPR1B, dominant white coat color, chromosome X, BMC Genomics