The study’s foundation is a sweeping computational survey across multiple public cancer databases. Using platforms including OncoMX, KM Plotter, UALCAN, and GEPIA2, the researchers compared TTK expression between tumor and normal tissues and then asked whether expression levels correlated with patient survival. The numbers were striking. In lung cancer, TTK messenger RNA was elevated with a log2 fold change of 4.10, meaning transcript levels were roughly an order of magnitude or more above normal tissue. Across cancers generally, TTK overexpression was a recurring theme, consistent with its established role as a spindle assembly checkpoint kinase that proliferating cells lean on heavily.

Survival analysis gave that overexpression clinical weight. Patients with elevated TTK showed significantly worse overall survival, with a hazard ratio of 1.62 and a p-value of 3.7e-15, an association so statistically robust that it is unlikely to be a fluke of the dataset. The effect was particularly pronounced in female smokers, a subgroup that has attracted growing attention because lung adenocarcinoma in women who smoke shows distinct molecular features from the classical smoking-driven disease. The authors also examined post-progression survival and first-progression endpoints, reinforcing the picture of TTK as a marker of aggressive disease rather than a passive byproduct of rapid cell division.

A crucial question for any transcriptomic biomarker is whether the RNA signal translates into protein. Here the team turned to the Human Protein Atlas and OncoDB, confirming that TTK protein levels are higher in lung adenocarcinoma tissue and that messenger RNA and protein abundance correlate strongly, with a correlation coefficient of 0.845. That tight coupling matters because it suggests the overexpression is not an artifact of transcriptional noise but reflects genuine accumulation of the kinase in tumor cells. The researchers also probed the mechanisms behind the deregulation, finding promoter hypomethylation at the TTK gene and hotspot mutations in the broader genomic landscape, both of which are plausible drivers of constitutive activation in a cancer context.

Perhaps the most provocative part of the study concerns the tumor immune microenvironment. Using TIMER 2.0 and CellTracer to estimate the composition of immune cells infiltrating tumors, the researchers found that high TTK expression correlated with reduced populations of CD8-positive cytotoxic T cells, the immune cells best able to kill tumor cells, with a correlation of minus 0.208. CD4-positive helper T cells were similarly depleted, with a correlation of minus 0.333, and hematopoietic stem cells showed an even stronger negative association of minus 0.658. Meanwhile, suppressive immune cell populations expanded in lockstep with TTK, showing a positive correlation of 0.65. Taken together, the pattern describes an immune-evasive microenvironment, the kind of landscape in which checkpoint inhibitor immunotherapies often struggle.

To explain how TTK might be regulated in the first place, the authors mapped a competing endogenous RNA, or ceRNA, network around the kinase. The ceRNA hypothesis holds that long non-coding RNAs can act as molecular sponges, soaking up microRNAs that would otherwise repress target messenger RNAs. In this case, the team identified a putative axis involving the transcription factor E2F1, which drives TTK transcription; the microRNA hsa-let-7b-5p, which would normally suppress TTK; and the long non-coding RNA TMPO-AS1, which appears to sponge hsa-let-7b-5p away from its target. When TMPO-AS1 is abundant, the microRNA is sequestered, TTK messenger RNA escapes repression, and kinase levels climb. Quantitative real-time PCR experiments in lung cancer cell lines confirmed the predicted pattern: TMPO-AS1 and TTK were upregulated while hsa-let-7b-5p was downregulated, lending experimental support to the computationally derived axis.

The therapeutic angle of the paper comes from molecular docking, a computational technique that predicts how tightly small molecules fit into a protein’s binding pocket. The researchers focused on the ATP-binding pocket of TTK, the same site targeted by clinical inhibitors in development. Two established TTK inhibitors, BAY 1217389 and CFI-402257, served as benchmarks, achieving docking scores of minus 9.128 and minus 8.417 kilocalories per mole respectively. Hesperidin, a flavonoid abundant in citrus peel, posted a docking score of minus 9.407 kilocalories per mole, slightly better than both clinical compounds. While docking scores are far from proof of drug activity, the comparison suggests that hesperidin binds the TTK pocket with comparable predicted affinity, raising the possibility of a natural-product scaffold for TTK inhibition or a combination therapy candidate that could be paired with existing agents.

TTK is not an arbitrary target. As the monopolar spindle 1 kinase, it is a central component of the spindle assembly checkpoint, the surveillance system that delays chromosome segregation until every chromosome is properly attached to the mitotic spindle. Tumors with high TTK activity can tolerate the chromosomal instability that fuels their evolution, and inhibiting the kinase in cancer cells pushes them through catastrophic mitosis. Several TTK inhibitors have already reached clinical testing in solid tumors, which is why the new study’s expression, survival, and immune-correlation data are timely: they help define which patients are most likely to carry TTK-dependent disease and provide a rationale for biomarker-guided trials.

The ceRNA axis adds a second layer of therapeutic logic. If TMPO-AS1 sustains TTK expression by neutralizing hsa-let-7b-5p, then the long non-coding RNA itself becomes a vulnerability. Antisense oligonucleotides or small molecules that strip away TMPO-AS1 would liberate the microRNA and allow endogenous repression of TTK, a strategy conceptually distinct from direct kinase inhibition. Conversely, restoring hsa-let-7b-5p, a member of a microRNA family with well-documented tumor-suppressive roles, could suppress TTK and potentially other oncogenic targets simultaneously. The study’s qRT-PCR confirmation of the axis in lung cancer cell lines is an early but necessary step toward validating these approaches in the laboratory.

As with any purely computational and cell-line-based study, the findings require caution in interpretation. Immune infiltration estimates are inferred from bulk tumor transcriptomes rather than directly measured, docking scores do not guarantee pharmacological activity, and the ceRNA hypothesis itself remains debated among RNA biologists. Nevertheless, the convergence of evidence in this paper, spanning pan-cancer expression, survival statistics, protein-level validation, epigenetic and mutational correlates, immune microenvironment mapping, regulatory network reconstruction, and inhibitor docking, makes a coherent case that TTK is a functional driver and actionable biomarker in lung adenocarcinoma rather than a bystander. If subsequent experimental and clinical work confirms the TMPO-AS1/hsa-let-7b-5p/TTK axis and the predicted activity of hesperidin, the study could open a genuinely new horizon for targeted therapy in a disease that urgently needs one, particularly for the female smoker subgroup in which the survival signal was strongest.

Subject of Research: TTK kinase as a prognostic biomarker and therapeutic target in lung adenocarcinoma

Article Title: Bioinformatic landscapes of TTK-mediated ceRNA networks: new horizons for targeted therapy in lung adenocarcinoma

Article References: Patidar, P., Baweja, B., Saini, C., Vats, P., Jangir, K., Kumar, N., & Nema, R. (2026). Bioinformatic landscapes of TTK-mediated ceRNA networks: new horizons for targeted therapy in lung adenocarcinoma. Journal of Cancer Research and Clinical Oncology. https://doi.org/10.1007/s00432-026-06641-3

Image Credits: AI Generated

DOI: 10.1007/s00432-026-06641-3

Keywords: TTK, lung adenocarcinoma, ceRNA network, TMPO-AS1, hsa-let-7b-5p, hesperidin, spindle assembly checkpoint, immune evasion, molecular docking, biomarker, non-small cell lung cancer, E2F1

The research team compared two carefully matched groups of nine calves each: one group of IVP calves with increased birth weight, and a control group of nine IVP calves with normal birth weight. Rather than relying on a single measurement, the investigators layered multiple analytical approaches on top of one another. They recorded physiological and biochemical blood parameters, profiled circulating hormones, sequenced the whole transcriptome to capture messenger RNAs, long non-coding RNAs, microRNAs and circular RNAs, and ran untargeted metabolomics to inventory the small molecules coursing through the calves’ systems. This multi-omics design matters because fetal overgrowth is unlikely to be governed by any single gene or metabolite; it is a systems-level phenomenon in which regulatory RNAs, protein-coding genes and metabolic pathways reinforce one another.

The first striking finding concerned gestation itself. Increased birth weight in the IVP calves was accompanied by prolonged gestation length in the recipient dams, meaning the oversized fetuses were not simply growing faster at a fixed endpoint but were also staying in utero longer. Alongside this, multiple physiological, biochemical and hormonal parameters shifted significantly between the two groups. The study tracked markers familiar to livestock physiologists, including albumin, alanine aminotransferase, creatinine, total cholesterol, triglycerides, high- and low-density lipoprotein cholesterol, non-esterified fatty acids and beta-hydroxybutyrate, as well as antioxidant measures such as catalase, glutathione peroxidase, superoxide dismutase and total antioxidant capacity. Hormonal profiling encompassed growth hormone and insulin-like growth factor 1, two central regulators of fetal growth, along with immune markers including immunoglobulins A and G and interleukins. The pattern that emerges is one of systemic metabolic and endocrine reprogramming, not merely a mechanical consequence of a large fetus occupying the uterus.

On the metabolomics side, the untargeted screen was remarkably deep. The platform detected 3,259 metabolites in total, from which the researchers distilled 132 differential metabolites that separated the heavy calves from the normal-weight controls. Of these, 82 were significantly up-regulated and 50 were significantly down-regulated. Among the metabolite classes represented were phosphatidylcholines, lysophosphatidylcholines, diacylglycerols and linoleic acid species, pointing toward perturbed lipid metabolism, a theme that resonates with the altered cholesterol, triglyceride and fatty acid measures in the blood biochemistry. Lipids are not only fuel; they are signaling molecules and structural components of every cell membrane, and their dysregulation during development can have lasting consequences for organ growth and function.

The transcriptomic analysis went beyond protein-coding genes to interrogate the non-coding RNA landscape, an area increasingly recognized as central to developmental programming. Using weighted gene co-expression network analysis, the team identified gene modules whose expression patterns correlated with birth weight, and they constructed a competing endogenous RNA, or ceRNA, network to map how circular RNAs might sponge microRNAs and thereby modulate messenger RNA targets. Two candidate axes stood out. The first links the circular RNA circ36279 to the microRNA bta-miR-431, with predicted target genes including SPON2, TNN, INSIG1 and AEN. The second connects circ35151 to bta-miR-154b, with predicted targets SHD, TBX2 and IL20RA. Each of these targets tells a plausible story: INSIG1 sits at the heart of cholesterol and lipid metabolism regulation, TBX2 is a transcription factor implicated in developmental patterning, IL20RA participates in inflammatory signaling, and SPON2 and TNN relate to extracellular matrix biology, the scaffolding that tissues remodel as they grow.

To bridge the gene-level and metabolite-level data, the researchers compared the KEGG pathways enriched independently among differentially expressed genes and among differential metabolites, and identified five overlapping pathways enriched on both sides of the omics divide. Within these shared pathways, Pearson correlation analysis of selected differentially expressed genes and differential metabolites revealed concrete cross-omics associations, pairs of genes and metabolites whose abundance moved together across the two groups of calves. The pathway that the authors single out as the strongest recommendation for further investigation is arginine and proline metabolism. This choice is biologically compelling: arginine is a precursor for nitric oxide, polyamines and creatine, all of which influence vascular tone, cell proliferation and energy metabolism in the developing fetus, while proline is central to collagen synthesis and thus to the growth of connective tissue and skeleton.

What makes this study notable in the broader context of assisted reproduction is its framing of fetal overgrowth as a molecular phenotype rather than a breeding-management nuisance. The phenomenon echoes large offspring syndrome, a well-documented complication of in vitro embryo technologies in ruminants, in which calves and lambs derived from cultured embryos show excessive birth weight, elongated gestation and metabolic abnormalities. By working exclusively within IVP calves, comparing heavy versus normal birth weight offspring produced under the same laboratory pipeline, the authors isolated the molecular correlates of overgrowth from the confounding effects of comparing IVP calves with artificially inseminated ones. That design choice strengthens the inference that the observed transcriptomic and metabolomic differences are genuinely tied to the overgrowth phenotype itself.

The clinical stakes for the dairy industry are tangible. Recipient dams carrying oversized IVP calves face higher rates of dystocia, and calving assistance carries costs in veterinary intervention, lost milk days, calf mortality and compromised welfare. If the ceRNA axes and metabolic pathways identified here can be validated as biomarkers, embryo production programs could potentially screen or select embryos, or manage recipient dams differently, to reduce the incidence of heavy births. More fundamentally, the identification of circ36279-bta-miR-431 and circ35151-bta-miR-154b as candidate regulatory circuits provides mechanistic entry points for experimental follow-up: manipulating these RNAs in embryo culture systems could test whether they causally influence growth trajectories rather than merely correlating with them.

The study also illustrates how far livestock genomics has moved in a decade. Whole-transcriptome sequencing that captures circular RNAs, once a niche pursuit in human biomedicine, is now being deployed on farm animals to solve practical production problems, and untargeted metabolomics capable of cataloguing more than three thousand compounds provides the chemical context in which genes act. The integration of these layers, anchored by physiological and hormonal measurements, reflects a maturing multi-omics paradigm in animal science. The authors, who also include Yao Chang, Lei Wang, Ruina Ma, Shanjiang Zhao, Qingshan Gao and Qing Xu, acknowledge support from the China National Modern Agricultural Industry Technology System, and the work was approved by the Animal Ethics Committee of China Agricultural University.

For now, the findings are associative, and the authors are appropriately measured in their conclusions, presenting the two ceRNA axes and the arginine and proline metabolism pathway as candidates for further examination rather than proven drivers of overgrowth. But the roadmap they provide is clear. Follow-up studies can validate the expression of SPON2, TNN, INSIG1, AEN, TBX2 and IL20RA, quantify the metabolites in the overlapping KEGG pathways, and test whether interventions in arginine availability during embryo culture or gestation moderate fetal growth. If those experiments succeed, the oversized IVP calf, long a stubborn side effect of a technology the dairy world increasingly depends upon, may finally become a problem that can be predicted, and perhaps prevented, at the molecular level.

Subject of Research: Molecular mechanisms of fetal overgrowth in Holstein calves produced by in vitro embryo production

Article Title: Exploring the transcriptomic and metabolomic associations with fetal overgrowth in Holstein cattle derived from in vitro embryo production

Article References: Qu, L., Hang, Z., Chang, Y., Wang, L., Ma, R., Zhao, S., Gao, Q., Xu, Q., & Wang, Y. (2026). Exploring the transcriptomic and metabolomic associations with fetal overgrowth in Holstein cattle derived from in vitro embryo production. BMC Genomics. https://doi.org/10.1186/s12864-026-13408-0

Image Credits: AI Generated

DOI: 10.1186/s12864-026-13408-0

Keywords: in vitro embryo production, Holstein cattle, fetal overgrowth, large offspring syndrome, transcriptomics, metabolomics, ceRNA network, circular RNA, microRNA, arginine and proline metabolism, lipid metabolism, multi-omics

T(6;11) translocation renal cell carcinoma is one of the more unusual members of the kidney cancer family. It arises when a piece of chromosome 6 carrying the TFEB gene is swapped with a piece of chromosome 11 carrying the Alpha gene, producing an abnormal fusion transcript called Alpha-TFEB. The fusion gene drives abnormally high levels of TFEB activity inside tumor cells, and clinicians have long observed that this subtype, first described in 2001 as a distinctive pediatric renal neoplasm, can behave aggressively in some patients while remaining indolent in others. Because the tumor is rare, accounting for only a small fraction of renal cell carcinomas, large clinical trials have been difficult to organize, and the molecular details of how the fusion gene actually promotes malignancy have remained only partially mapped. The new study was designed to fill in one crucial stretch of that map.

The research team, led by He-Qin Zhan of the Department of Pathology at Anhui Medical University, with Hai-Yan Ni, Xiao-Xun Zhang and Xin-Ran Wang as co-first authors, built on cell lines their group had established in earlier work. These lines, designated HK-2-TFEB and Caki-2-TFEB, carry stable expression of the Alpha-TFEB fusion gene and serve as laboratory models of the translocation cancer. Western blotting confirmed that TFEB protein was abundant in both engineered lines compared with matched control cells. When the investigators measured RNA levels, they found that LINC00612, a long intergenic non-protein coding RNA, was significantly elevated in the fusion-positive cells, while miR-31-5p was reduced and CTNNB1, the gene encoding beta-catenin, was increased. That reciprocal pattern, high lncRNA, low microRNA, high target gene, is the classic signature of a competing endogenous RNA, or ceRNA, network.

The ceRNA hypothesis, first articulated by Salmena and colleagues in 2011, proposes that long non-coding RNAs can communicate with protein-coding messenger RNAs by competing for the same microRNA molecules. MicroRNAs are short RNA strands that normally bind to target messenger RNAs and silence them, so any RNA that binds and sequesters a microRNA effectively lifts the brake on that microRNA’s real targets. LINC00612 has previously been implicated in this kind of crosstalk in osteosarcoma, bladder cancer and colon adenocarcinoma, and a separate study linked a LINC00612/miR-31-5p axis to endothelial cell injury in cigarette smoke exposure. What the new work adds is the upstream driver: in this kidney cancer, the entire circuit appears to be switched on by the transcription factor TFEB itself.

Establishing that upstream connection required two complementary molecular techniques. Dual-luciferase reporter assays, in which the LINC00612 promoter is wired to a light-producing enzyme, showed that TFEB binding specifically activated the promoter, and chromatin immunoprecipitation experiments physically pulled down the promoter DNA along with TFEB protein, confirming that the transcription factor docks directly onto the LINC00612 regulatory region inside cells. In other words, TFEB is not merely correlated with LINC00612 expression; it is transcriptionally commanding it. Because the Alpha-TFEB fusion gene is the defining genetic lesion of this tumor subtype, the result provides a satisfying causal chain: the chromosome translocation creates the fusion, the fusion elevates TFEB activity, and TFEB turns up the volume on LINC00612.

With the top of the cascade established, the team tested what each link actually does to cancer cell behavior. When they overexpressed LINC00612 in ordinary HK-2 kidney cells, the cells became more proliferative, more migratory in wound healing assays and more invasive through membrane barriers, as measured by CCK-8 proliferation and transwell invasion assays. Conversely, when they knocked LINC00612 down in the fusion-positive HK-2-TFEB cells, those malignant capabilities were suppressed. The researchers also verified the middle and bottom of the relay: LINC00612 was shown to bind miR-31-5p directly, and miR-31-5p was shown to target CTNNB1. Reducing CTNNB1 expression likewise dampened the malignant functions of the cells, consistent with beta-catenin, a central player in the Wnt signaling pathway, being the downstream effector that ultimately drives the cancer phenotype.

Cell culture, however convincing, is only half the story, so the team turned to animal models to see whether the same axis governs tumor growth in living tissue. In mice, overexpression of either TFEB or LINC00612 promoted tumor growth, while overexpression of miR-31-5p inhibited it. The animal experiments, approved by the Animal Ethics Committee of Anhui Medical University, thus recapitulated the in vitro findings and demonstrated that manipulating any single node of the pathway, the transcription factor, the long non-coding RNA or the microRNA, is sufficient to change the trajectory of tumor growth in vivo. The consistency across molecular, cellular and animal levels gives the proposed TFEB/LINC00612/miR-31-5p/CTNNB1 axis unusual robustness for a study of such a rare malignancy.

The therapeutic implications are potentially significant. Because each component of the axis is a distinct molecular species, the pathway offers several points of attack. Antisense oligonucleotides or small interfering RNAs could theoretically silence LINC00612, microRNA mimics could restore miR-31-5p levels, and agents already in development against Wnt/beta-catenin signaling could blunt the downstream effect of CTNNB1. Notably, miR-31-5p has previously been described as a tumor suppressor in renal cell carcinoma through its targeting of cyclin-dependent kinase 1, which fits neatly with the new finding that its depletion by LINC00612 removes a brake on tumor progression. Prior work from the same group had also shown that the Alpha gene upregulates TFEB and promotes malignant transformation, making the present study a natural extension of a decade-long investigation into this fusion-driven cancer.

The study also contributes to a broader shift in cancer biology. Long non-coding RNAs, once dismissed as transcriptional noise, are increasingly recognized as master regulators of tumor behavior, and ceRNA networks have now been documented in hepatocellular carcinoma, gastric cancer, cervical cancer, prostate cancer and clear cell renal cell carcinoma, among others. By tying a specific lncRNA to a defining chromosomal translocation, the Anhui team has shown how a structural genetic lesion can be translated, through non-coding RNA circuitry, into the everyday business of cancer cell proliferation and spread. That perspective may also help explain why some t(6;11) tumors behave aggressively while others do not: variation in the activity of the LINC00612/miR-31-5p module could modulate the malignant potential that the fusion gene confers.

Caveats remain, as they always do at this stage of translation. The work rests on engineered cell lines and mouse models rather than large panels of patient tumors, and the rarity of t(6;11) renal cell carcinoma means that validating the axis in clinical specimens will require multi-institution collaboration. Whether blocking LINC00612 in patients is feasible and safe is unknown, and microRNA-based therapies have historically faced delivery challenges. Nevertheless, the identification of a complete, experimentally verified signaling relay, from the Alpha-TFEB fusion gene through LINC00612 and miR-31-5p to CTNNB1, gives researchers a concrete molecular target where previously there was only a correlation. For patients with this rare kidney cancer, whose treatment options have been limited by the tumor’s resistance to conventional therapies, the study offers something precious: a map of the enemy’s command structure, drawn molecule by molecule, and several new places to strike.

Subject of Research: The TFEB/LINC00612/miR-31-5p/CTNNB1 signaling axis in t(6;11) translocation renal cell carcinoma progression

Article Title: TFEB promotes the progression of t(6; 11) translocation renal cell carcinoma by upregulating LINC00612 expression

Article References: Ni, H.-Y., Zhang, X.-X., Wang, X.-R., Ma, J., Yu, X.-J., & Zhan, H.-Q. (2026). TFEB promotes the progression of t(6; 11) translocation renal cell carcinoma by upregulating LINC00612 expression. Medical Oncology, 43(11), Article 304. https://doi.org/10.1007/s12032-026-03411-6

Image Credits: AI Generated

DOI: 10.1007/s12032-026-03411-6

Keywords: TFEB, LINC00612, miR-31-5p, CTNNB1, t(6;11) translocation renal cell carcinoma, long non-coding RNA, ceRNA network, Alpha-TFEB fusion gene, renal cancer, Wnt signaling, tumor progression, Medical Oncology

The research team, led by scientists at Guangxi University in collaboration with the Guangxi Buffalo Research Institute and international partners, performed an integrated transcriptomic analysis on milk-derived RNA from Murrah buffaloes at mid-lactation. The animals were divided into two clearly separated groups based on lactation performance. High-yield buffaloes produced an average of roughly 3,804 kilograms of milk per lactation, while low-yield animals averaged about 1,928 kilograms, a difference of nearly a factor of two. By sequencing the complete RNA content of milk from both groups, the researchers set out to identify the transcriptional features that accompany, and possibly underpin, this dramatic divergence in productivity.

The choice of milk as the sampling medium is one of the study’s most practical innovations. Mammary epithelial cells and other milk-associated cells are shed continuously into milk, meaning a simple sample of the fluid carries a molecular record of the secretory activity inside the udder. Unlike a mammary gland biopsy, milk collection is non-invasive, repeatable, and stress-free for the animal, which matters both for animal welfare and for the feasibility of scaling such measurements to breeding programs. The study’s authors note that no animals were euthanised, anaesthetised, or subjected to invasive procedures, and that the work was approved by the Animal Experiment Ethics Review Committee of Guangxi University.

Sequencing the milk RNA revealed a remarkably complex transcriptomic landscape. In total, the team identified 17,735 messenger RNAs, the protein-coding workhorses of the cell, alongside 2,641 long non-coding RNAs, 6,668 circular RNAs, and 578 microRNAs. This catalogue of four distinct RNA classes is what makes the study integrated rather than a conventional single-layer gene expression analysis. Long non-coding RNAs and circular RNAs can act as molecular sponges that bind microRNAs, thereby modulating how much protein gets produced from messenger RNA transcripts, a regulatory architecture known as the competing endogenous RNA, or ceRNA, network. Capturing all of these layers simultaneously from the same samples allowed the researchers to begin mapping how they interact.

The comparative analysis between high-yield and low-yield animals produced a coherent biological story. Buffaloes with higher milk output showed coordinated upregulation of genes involved in oxidative phosphorylation, the mitochondrial process that converts nutrients into cellular energy in the form of ATP. Genes tied to glucose metabolism, cellular metabolism in general, and protein synthesis were also elevated in the high-yield group. Taken together, these patterns point to an enhanced bioenergetic and biosynthetic state in the milk-associated cells of productive animals. In essence, the secretory cells of high-yielding buffaloes appear to run a hotter metabolic engine, burning more glucose through their mitochondria and building more protein, which is exactly what one would expect from tissue tasked with secreting roughly twice the volume of milk.

Equally telling was what was turned down. Transcripts associated with signal transduction, transporter activity, and apoptosis-related processes showed reduced expression in the high-yield group. Lower expression of apoptosis-related genes suggests that the milk cells of productive animals may resist programmed cell death, potentially sustaining secretory activity for longer. Reduced signalling and transporter transcripts hint at a shift in cellular priorities: rather than communicating with their environment or shuttling molecules across membranes, the cells of high-yield animals appear to channel resources into synthesis and secretion. The researchers are careful to frame these as associations rather than proven causes, describing the work as an exploratory framework, but the consistency of the metabolic signature across functional categories lends it weight.

By integrating the messenger RNA data with the non-coding RNA layers, the team constructed putative ceRNA networks that connect microRNAs, their messenger RNA targets, and the long non-coding and circular RNAs that sequester them. Within these networks, two candidate genes emerged as potentially associated with milk yield-related cellular pathways: SNX25, which encodes sorting nexin 25, a protein involved in intracellular trafficking, and DHPR, which encodes dihydropyridine reductase, an enzyme in the metabolic recycling pathway. These genes now stand as priority targets for future functional validation, the experimental step needed to determine whether they merely correlate with high yield or actively contribute to it.

The implications for buffalo breeding are potentially significant. Traditional genetic selection for milk yield in buffalo has been slowed by long generation intervals, limited pedigree and genomic records in many producing countries, and the sheer difficulty of measuring lactation traits early in an animal’s life. A transcriptomic marker panel derived from milk, measurable without harming the animal and potentially early in a lactation, could complement genomic selection by providing a functional readout of the actual secretory machinery. If the metabolic and ceRNA signatures identified here hold up in larger cohorts and across breeds, they could help identify high-potential animals sooner and guide the development of genetic improvement strategies tailored to buffalo biology rather than borrowed wholesale from dairy cattle.

The study also contributes a substantial public resource to the emerging field of livestock transcriptomics. Cataloguing nearly 18,000 messenger RNAs and thousands of non-coding transcripts from buffalo milk provides a reference for other researchers studying mammary biology, lactation physiology, and even milk quality traits, since many of the same metabolic pathways that drive yield also influence fat and protein composition. The raw sequencing data have been deposited in the Sequence Read Archive, allowing independent verification and reuse. The work was supported by the Guangxi Key Research and Development Program and the National Natural Science Foundation of China, reflecting the strategic importance of buffalo dairy production in southern China and across South Asia.

As with any exploratory transcriptomic study, caveats remain. The findings describe molecular features associated with divergent yield phenotypes at a single lactation stage in one breed, and the ceRNA networks are computational predictions that require experimental confirmation in mammary cells. The authors themselves position the work as a framework and a source of candidate targets rather than a definitive mechanism. Even so, the study marks a step toward a future in which a routine milk sample, already collected on every dairy farm every day, doubles as a molecular diagnostic window into the animal’s productive capacity. For an industry built on the biology of secretion, reading the transcriptome in the milk itself may prove to be one of the most elegant shortcuts breeding science has found in decades.

Subject of Research: Transcriptomic analysis of milk-derived RNA linked to divergent milk yield in buffalo

Article Title: Integrated transcriptomic analysis of milk-derived RNA associated with divergent milk yield phenotypes in buffalo

Article References: Song, J., Wang, S., Jiang, H., Li, L., Jiang, J., Teng, C., Xiong, J., Wang, L., Lu, Y., Pauciullo, A., Liu, Q., & Li, Z. (2026). Integrated transcriptomic analysis of milk-derived RNA associated with divergent milk yield phenotypes in buffalo. BMC Genomics. https://doi.org/10.1186/s12864-026-13363-w

Image Credits: AI Generated

DOI: 10.1186/s12864-026-13363-w

Keywords: buffalo, milk yield, transcriptomics, RNA sequencing, non-coding RNAs, ceRNA network, oxidative phosphorylation, lactation, molecular breeding, gene expression, Murrah buffalo, BMC Genomics

Circular RNAs are an unusual class of molecules. Unlike ordinary linear messenger RNAs, which are transcribed, translated, and degraded in a straightforward sequence, circRNAs are formed when the splicing machinery of the cell joins the downstream end of an RNA molecule back to its upstream end, a process known as back-splicing. The result is a covalently closed loop with no free ends, which makes the molecule remarkably resistant to degradation by the exonucleases that normally chew up RNA. First noticed decades ago as rare splicing accidents, circRNAs are now recognized as abundant, conserved, and often functional components of the transcriptomes of animals and plants alike, with roles that include sponging microRNAs, modulating transcription, and in some cases even serving as templates for translation.

In the new work, Behzad Hajieghrari of Jahrom University and Mousa Torabi Giglou of the University of Mohaghegh Ardabili in Iran systematically reanalyzed strand-specific RNA sequencing data from rice exposed to simultaneous heat and drought stress. Because circRNAs lack the poly-A tails and defined ends of linear transcripts, detecting them requires specialized computational approaches. The researchers subjected the sequencing reads to rigorous quality control, mapped them to the rice genome, and then ran two independent circRNA prediction algorithms, CIRI2 and CIRCexplorer2, in parallel. Only candidates supported by both methods were retained, a dual-algorithm strategy designed to filter out false positives arising from repetitive sequence or misaligned reads. The screen yielded 208 high-confidence circRNAs distributed across all twelve rice chromosomes.

The comparative profiling revealed a striking pattern of stress-dependent circularization. Eighty-three circRNAs were detected exclusively in unstressed control samples, fifty-one appeared only in stressed samples, and seventy-four were shared between the two conditions. In other words, the circular transcriptome is not a static backdrop; it is remodeled when the plant senses combined stress. Junction-read analysis, which counts the sequencing reads that span the diagnostic back-splice junction, exposed a spectrum of circularization strength, from highly abundant circRNAs whose junction reads dominate their genomic loci to low-confidence candidates whose signals are partially masked by the background of linear transcripts from the same genes.

Genomic annotation showed that most of the identified circRNAs originated from exonic regions, with a substantial contribution from intergenic regions, and displayed a pronounced bias toward the negative DNA strand. Several host genes produced multiple distinct circRNA isoforms through alternative back-splicing, meaning that a single gene can generate a small family of circular molecules with potentially different regulatory partners. This isoform diversity adds a layer of complexity to the rice transcriptome that linear RNA analysis alone cannot capture, and it hints that alternative circularization may itself be a regulated process that the plant tunes under stress.

To infer what these circRNAs might be doing, the team performed functional enrichment analysis on their host genes. The results pointed to involvement in protein folding, nutrient reservoir activity, RNA degradation, and branched-chain amino acid catabolism. Each of these categories makes biological sense in the context of combined heat and drought. Protein folding machinery, including heat shock proteins, is central to surviving thermal damage; nutrient reservoir proteins reflect the metabolic reallocation that stress demands; RNA degradation pathways govern how quickly stress transcripts turn over; and amino acid catabolism connects to nitrogen mobilization and osmotic adjustment. The enrichment pattern suggests that circRNAs are not random byproducts but are embedded in the metabolic and proteostatic circuits that determine whether a rice plant tolerates or succumbs to compound stress.

Differential expression analysis between control and stressed libraries identified seven circRNAs specifically induced under combined heat and drought conditions. These stress-responsive candidates represent the most direct leads for future experimental work, since their induction implies that the plant actively upregulates their production as part of its adaptive program. Whether they act by sequestering microRNAs, interacting with RNA-binding proteins, or influencing the splicing of their own host genes remains to be determined, but their stress-specific behavior marks them as priority targets for functional validation.

The most intriguing part of the study concerns the predicted regulatory network built around these molecules. Drawing on the competing endogenous RNA, or ceRNA, hypothesis, the researchers predicted which microRNAs could bind each circRNA and which messenger RNAs those microRNAs could in turn regulate. In the ceRNA framework, a circRNA with binding sites for a particular microRNA can act as a molecular sponge, soaking up that microRNA and thereby relieving repression of the microRNA’s genuine mRNA targets. Network topology analysis of the resulting three-layer circRNA-microRNA-mRNA circuit pinpointed several microRNAs, including osa-miR414, osa-miR1439, and osa-miR2919, as candidate topological hubs, meaning they occupy central positions with many connections and could exert outsized influence over the network’s behavior. The predicted targets of these hub microRNAs encode stress-responsive transcription factors and signaling proteins, suggesting a plausible route by which circRNA abundance changes could ripple outward to reshape the expression of entire stress-response gene programs.

The authors are careful to frame these network findings as predictive. The ceRNA relationships were inferred computationally rather than demonstrated experimentally, and microRNA target prediction in plants, while reasonably accurate due to near-perfect base pairing requirements, still generates false positives. Nevertheless, the study delivers the first comprehensive map of circRNAs in rice under combined heat and drought stress, and it does so with a methodological transparency that should make follow-up work straightforward: the full lists of predicted circRNAs, their genomic coordinates, sequences, host gene annotations, differential expression results, and predicted interaction networks are all provided in supplementary data files. The researchers also note that the work received no external funding and was carried out with resources covered by the authors themselves.

The broader significance lies in what this means for crop improvement. Extreme weather events increasingly combine heat and water deficit during the rice growing season, and breeding for tolerance to each stress individually has not reliably produced varieties that withstand the combination. If circRNA-mediated regulation proves to be a genuine coordinating mechanism, it opens a new class of molecular markers and, eventually, engineering targets: circRNAs or the splicing elements that control their production could be tuned to bolster the plant’s post-transcriptional defenses. For now, the seven stress-induced circRNAs and the hub microRNAs osa-miR414, osa-miR1439, and osa-miR2919 constitute a concrete experimental agenda. Validating their interactions, confirming their sponging activity, and testing their effects on stress tolerance in living rice plants will determine whether these molecular rings are merely correlates of stress or true architects of the crop’s resilience.

Subject of Research: Circular RNA-mediated post-transcriptional regulation in rice under combined heat and drought stress

Article Title: Circular RNAs orchestrate integrated post-transcriptional responses to combined heat and drought stress in rice

Article References: Hajieghrari, B., & Giglou, M. T. (2026). Circular RNAs orchestrate integrated post-transcriptional responses to combined heat and drought stress in rice. Molecular Genetics and Genomics, 301(1), Article 183. https://doi.org/10.1007/s00438-026-02516-x

Image Credits: AI Generated

DOI: 10.1007/s00438-026-02516-x

Keywords: circular RNAs, rice, heat stress, drought stress, ceRNA network, microRNAs, back-splicing, RNA-seq, Oryza sativa, post-transcriptional regulation, plant stress responses, Molecular Genetics and Genomics