The genes in question belong to the WRKY family, one of the largest families of transcription factors in plants. Transcription factors are proteins that bind to specific DNA sequences and switch other genes on or off, acting like the conductors of a genetic orchestra. WRKY factors are named for a highly conserved stretch of amino acids, the WRKY motif, and they are known to orchestrate a plant’s reactions to drought, cold, disease, and salinity. Until now, however, no systematic catalogue of WRKY genes existed for Codonopsis pilosula, leaving a gap in understanding how this economically significant edible and medicinal plant mounts its defense against soil stress.

To fill that gap, the researchers scanned the entire genome of the species and identified 47 WRKY genes, which they named CpWRKYs. The genes were not spread evenly across the plant’s genetic blueprint; instead, they were distributed unevenly across eight chromosomes, a pattern that often reflects the uneven pressures of evolution. Using phylogenetic analysis, which compares gene sequences to reconstruct family trees, together with examinations of conserved protein domains and motifs, the team sorted the 47 genes into seven subfamilies. Members of each subfamily shared similar conserved motifs but differed in their gene structures, the arrangement of coding regions and intervening sequences that shapes how each gene is read and regulated.

A second layer of analysis focused on how the family grew over evolutionary time. Collinearity analysis, which looks for stretches of DNA that appear in duplicated blocks across the genome, revealed 13 duplicated gene pairs among the CpWRKYs. This pointed to segmental duplication, the replication of large chromosomal segments, as the main engine driving the expansion of the WRKY family in this species. Such duplication events give evolution raw material to work with: one copy can keep performing the ancestral job while the other is free to mutate and take on new roles, including the fine-tuning of stress responses.

Identifying the genes was only the first step. The real question was which of them actually respond when the plant’s roots hit salty, alkaline soil. To answer it, the researchers turned to transcriptome analysis, a technique that measures which genes are actively being transcribed into RNA under given conditions. When they compared plants exposed to saline-alkali stress with unstressed controls, they found that the differentially expressed WRKY genes, the ones whose activity changed markedly, mostly belonged to three groups: I, IIc, and III. These groups are known in other species to be involved in stress signaling, suggesting that Codonopsis pilosula recruits a familiar cast of molecular players to deal with an unfamiliar challenge.

But gene expression alone does not tell the whole story. Plants respond to stress not only by changing which genes are active but also by altering their metabolism, the chemical reactions that keep cells alive. The team therefore performed metabolic profiling to see which small molecules accumulated or declined under saline-alkali conditions. The results were striking: lipid metabolites, the fat-derived molecules that make up cell membranes and serve as signaling compounds, were the main differentially accumulated metabolites under stress. This makes biological sense, because saline-alkali stress damages membranes and triggers a cascade of lipid-based signaling that helps the plant reorganize its cellular defenses.

The most intriguing finding came when the researchers connected the two data streams. Using gene co-expression network analysis, a computational approach that identifies groups of genes whose activity rises and falls in tandem, they discovered that the genes involved in lipid metabolism were linked to specific WRKY transcription factors. One gene in particular, Cpi05G005333-I, emerged as a hub in this network, suggesting that it may act as a master switch coordinating the lipid-based response to saline-alkali stress. In other words, the WRKY factors appear to sit upstream of the metabolic changes, turning on the machinery that reshapes the plant’s lipid profile when the soil turns hostile.

By combining phylogenetic classification, duplication analysis, transcriptomics, metabolomics, and network modeling, the study offers the first comprehensive portrait of the WRKY family in Codonopsis pilosula and identifies candidate genes that likely participate in the transcriptional response to saline-alkali stress. The work was carried out by researchers at Changzhi Medical College in China, with additional contributions from the Inner Mongolia Key Laboratory of Plants Adversity Adaptation and Genetic Improvement in Cold and Arid Regions at Inner Mongolia Agricultural University, and was supported by funding from Changzhi Medical College and the Fundamental Research Program of Shanxi Province.

The practical implications could be significant. Codonopsis pilosula is not only a staple of traditional herbal medicine but also a crop whose cultivation is expanding into marginal, salt-affected soils. Understanding which WRKY genes drive salt tolerance opens the door to breeding or engineering varieties that can withstand saline-alkali conditions more effectively, whether through conventional selection informed by genetic markers or through targeted genome editing. The lipid connection is particularly interesting because membrane remodeling is a shared feature of stress responses across many plant species, meaning the findings could inform broader efforts to develop stress-resilient crops.

There is still much to confirm. The study identifies candidate genes and networks, but functional validation, the laboratory work of silencing or overexpressing individual genes to see how the plant responds, will be needed to prove which CpWRKYs are truly indispensable for salt tolerance. Nevertheless, the research lays a solid foundation. It demonstrates how a genome-wide, multi-omics approach can turn a medicinal herb into a model for understanding stress adaptation, and it brings scientists one step closer to the goal of growing valuable crops on land that would otherwise remain barren.

Subject of Research: Genome-wide identification of WRKY transcription factors and their role in saline-alkali stress response in Codonopsis pilosula roots

Article Title: Genome-wide identification of CpWRKYs and elucidation of the saline-alkali response mechanism in Codonopsis pilosula roots

Article References: Liu, J., Zhang, X., Jing, X., Zhao, L., Yan, Y., Song, X., Li, Y., Su, W., & Wan, Y. (2026). Genome-wide identification of CpWRKYs and elucidation of the saline-alkali response mechanism in Codonopsis pilosula roots. BMC Plant Biology. https://doi.org/10.1186/s12870-026-09770-6

Image Credits: AI Generated

DOI: 10.1186/s12870-026-09770-6

Keywords: Codonopsis pilosula, WRKY transcription factors, saline-alkali stress, genome-wide analysis, transcriptomics, lipid metabolism, gene co-expression network, segmental duplication, plant stress response, medicinal plants, BMC Plant Biology, transcriptional regulation

The study, published in the journal Veterinary Oncology, was led by Ingrid Marie Moberg and colleagues at the Norwegian University of Life Sciences, Oslo University Hospital, and the University of Oslo. The team analyzed RNA sequencing data from 179 canine mammary tumors, drawn from a cohort of naturally occurring tumors removed from companion dogs during routine surgery. Their goal was twofold: to define molecular subgroups of canine mammary tumors without reference to histology, and then to ask whether tumors from the same dog share those molecular identities. The answer to the second question, in most cases, was no.

To classify the tumors, the researchers first reduced their high-dimensional gene expression data using principal component analysis, then applied unsupervised consensus clustering, an algorithm that groups samples based purely on expression similarity and stability across repeated resampling. This revealed five robust transcriptomic clusters. Crucially, the clusters did not map cleanly onto histological diagnosis. Each cluster contained a mixture of benign and malignant tumors, confirming that the microscopic appearance of these tumors tells only part of the story of their underlying biology.

Gene set enrichment analysis then revealed that the five clusters bear a remarkable resemblance to molecular subtypes long recognized in human breast cancer. Three of the clusters showed enrichment of hormone-related gene programs, including estrogen response, and were found to harbor the transcription factor GATA3, a classic marker of luminal breast cancer in humans. One cluster combined cell-cycle proliferation with immune and interferon signaling, evoking the aggressive basal-like or triple-negative phenotype, while a fifth cluster was defined by epithelial-mesenchymal transition, angiogenesis, and inflammatory modules, reminiscent of the claudin-low subtype. The parallel with human breast cancer taxonomy is not merely cosmetic; it strengthens the case for dogs as a comparative model in which the biology of mammary cancer can be studied in a spontaneously arising disease.

Beyond clustering, the team constructed a gene co-expression network using the hCoCena framework, identifying ten modules of genes that are expressed together and that encode distinct biological processes. One module captured hormone signaling, two captured immune and interferon responses, three reflected metabolism and proliferation through glycolysis, PI3K-AKT-mTOR signaling, oxidative phosphorylation, and MYC targets, and others traced epithelial differentiation through WNT signaling, epithelial-mesenchymal transition, cell-cycle regulation through E2F targets and the G2M checkpoint, and tumor microenvironment features such as angiogenesis. Transcription factor enrichment within the modules pointed to GATA3 as a regulator of the hormonal program and E2F1 as a driver of the proliferative module, providing candidate master switches behind the observed phenotypes.

The heart of the study, however, lies in its analysis of synchronous tumors. The researchers focused on 45 dogs that each carried exactly two tumors, yielding 90 paired samples classified as benign-benign, malignant-benign, or malignant-malignant. When they compared cluster assignments within each pair, concordance was low. Only about 45 percent of tumor pairs landed in the same transcriptomic cluster overall, with malignant-malignant pairs showing the highest agreement at 56 percent, benign-benign pairs at 44 percent, and mixed malignant-benign pairs at just 36 percent. In other words, the majority of dogs carried two tumors with fundamentally different molecular identities.

To quantify this divergence at the level of individual genes, the team calculated intraclass correlation coefficients for every gene across the paired samples. Genes were scored as low, moderate, or high in correlation between a dog’s two tumors. The results were unambiguous: between roughly 76 and 90 percent of genes showed low correlation across all diagnostic categories, and fewer than one percent of genes were highly correlated within any category. This pattern of widespread discordance held regardless of whether both tumors were benign, both malignant, or one of each, indicating that molecular independence between synchronous tumors is the norm rather than the exception.

A small set of exceptions proved informative. The analysis identified a handful of genes, including OMD and EN1, whose expression is strongly correlated within certain categories of synchronous tumor pairs and which have been reported as prognostic markers in human cancers. The authors suggest that these genes may point to shared disease processes or protective mechanisms against malignant transformation, and that they deserve further investigation, potentially at the DNA level, to uncover any genetic factors underlying their coordinated behavior. Meanwhile, examination of module-level variation showed that programs linked to cell-cycle activity, hormone signaling, and immune responses fluctuated most between paired tumors, while modules tied to epithelial differentiation and metabolism remained comparatively stable within individuals.

The clinical implications of this work reach in two directions at once. For veterinary medicine, the findings suggest that each tumor in a multi-tumor patient should be evaluated as a biologically independent lesion rather than assumed to be representative of its neighbors. Because molecular subtypes in human breast cancer drive dramatically different treatment decisions, from endocrine therapy for hormone receptor-positive disease to chemotherapy and PARP inhibitors for triple-negative tumors, a biology-driven approach could eventually refine the limited therapeutic options currently available for canine patients. The authors note, for instance, that identifying hormone-positive subtypes could inform whether ovariohysterectomy offers real benefit at the time of tumor removal, a decision veterinarians currently make without molecular guidance.

For comparative oncology, the study reinforces the value of the canine model in a way that human cohorts cannot easily replicate. Synchronous bilateral breast cancer occurs in only around one percent of human patients, whereas synchronous mammary tumors affect the majority of affected dogs. Studying these paired tumors within the same genetic background eliminates many host-specific confounders, offering a natural experiment in tumor evolution and inter-individual heterogeneity. The authors acknowledge limitations inherent to bulk RNA sequencing, which averages expression across all cells in a tissue and may reflect tissue composition as much as tumor biology, and they caution that a single RNA sample may not represent an entire tumor. Even so, their conclusion stands firm: histopathology alone does not capture the molecular reality of canine mammary tumors, and expression-based profiling offers a more faithful map of the biological terrain that clinicians and researchers alike will need to navigate.

Subject of Research: Transcriptomic heterogeneity of synchronous canine mammary tumors and their molecular resemblance to human breast cancer subtypes

Article Title: High molecular heterogeneity in synchronous canine mammary tumors detected by transcriptomic analysis

Article References: Moberg, I. M., Murphy, S. L., Hansen, N., Borge, K. S., Gunnes, G., Sørlie, T., Lingaas, F., Bergholtz, H., & Solbakken, M. H. (2026). High molecular heterogeneity in synchronous canine mammary tumors detected by transcriptomic analysis. Veterinary Oncology, 3(1), Article 11. https://doi.org/10.1186/s44356-026-00066-3

Image Credits: AI Generated

DOI: 10.1186/s44356-026-00066-3

Keywords: canine mammary tumors, transcriptomics, RNA sequencing, tumor heterogeneity, synchronous tumors, breast cancer subtypes, gene co-expression network, consensus clustering, comparative oncology, veterinary pathology, High, molecular