The study began with a large-scale computational sweep through publicly available blood transcriptome data. Using the limma statistical framework, the team compared gene expression in severe COVID-19 patients against healthy controls, drawing on dataset GSE171110, which comprised 44 patient samples and 10 controls. A parallel analysis of the gout dataset GSE160170, with six patients and six controls, applied the same thresholds: an adjusted p-value below 0.05 and a log2 fold change of at least 1. The COVID-19 cohort yielded 3,520 differentially expressed genes, while the gout cohort produced 1,860. Intersecting the two lists revealed just 19 genes shared by both inflammatory disorders — a compact transcriptional module that the authors interpret as the core molecular signature driving the convergent inflammatory response.

Functional enrichment analysis of these 19 shared genes pointed strongly toward the immune system. Gene Ontology terms highlighted B cell proliferation, immunoglobulin biosynthesis, cytokine activity, and cytokine receptor binding. This B-cell signal is particularly intriguing. In severe COVID-19, aberrant B cell activation and autoantibody production have been recognized as key drivers of immunopathology and lung damage. Gout, by contrast, has traditionally been viewed as a disease of innate immunity, driven by neutrophils and macrophages responding to urate crystals. Yet emerging evidence implicates B cells and autoantibody responses in chronic tophaceous gout as well. The convergence on humoral immunity in both diseases suggests a shared mechanism that has been underappreciated until now. KEGG pathway analysis further linked the shared gene set to autoimmune conditions including rheumatoid arthritis, asthma, and inflammatory bowel disease, underscoring the systemic nature of the overlapping inflammatory program.

To distill the 19 shared genes down to actionable targets, the researchers constructed a protein-protein interaction network using the STRING database at a high-confidence threshold of 0.7. The resulting network contained 19 nodes and 41 edges, with CXCL8, JUN, and IL10 showing the highest connectivity. Applying the CytoNCA plugin in Cytoscape, the team evaluated five centrality metrics — degree, closeness, eigenvector, local average connectivity, and betweenness — and identified six genes that ranked highest across all measures: CXCL8, MPO, JUN, MMP1, IL4, and IL10. These hub genes were significantly upregulated in both disease cohorts compared with healthy controls, and correlation heatmaps revealed strong positive co-expression among them in each dataset, suggesting they operate as a synergistic regulatory unit rather than isolated players.

Enrichment analysis of the six hub genes reinforced their relevance to both pathologies. Gene Ontology terms included positive regulation of MHC class II biosynthetic process, myeloid leukocyte activation, and regulation of isotype switching — processes central to adaptive immunity and implicated in the dysregulated humoral responses of severe COVID-19. KEGG pathway analysis tied the hub genes to the IL-17 signaling pathway, T cell receptor signaling, and rheumatoid arthritis. The IL-17 axis, in particular, has been independently implicated in both gout pathogenesis and COVID-19-associated hyperinflammation, positioning these six genes as convergent nodes within a shared inflammatory network.

The pivotal moment came when the team intersected these six hub genes with the predicted protein targets of colchicine, the ancient anti-inflammatory drug long used to treat acute gout attacks and recently repurposed for COVID-19 with encouraging results in trials such as GRECCO-19. Using the SwissTargetPrediction platform, the researchers generated a list of 100 putative colchicine targets from its PubChem entry. Only one of the six hub genes appeared on that list: MMP1. This single overlap suggested that colchicine’s therapeutic benefit in both diseases might flow, at least in part, through direct engagement of this matrix metalloproteinase — a mechanism not previously described.

To test this hypothesis computationally, the team performed molecular docking between colchicine and the MMP1 protein structure obtained from the Protein Data Bank. The docking simulation, run with AutoDock Vina, revealed a stable binding conformation with a calculated binding energy of −6.4 kcal/mol, a value within the range typically considered favorable for drug-target interactions. Colchicine formed conventional hydrogen bonds with the ASN-180 and ALA-182 residues of MMP1, at distances of 2.3 and 2.6 angstroms respectively, and was further anchored by a network of hydrophobic contacts within the catalytic pocket. Because these hydrogen bonds sit within the enzyme’s catalytic domain, the authors propose that colchicine binding may directly interfere with MMP1’s proteolytic activity, thereby modulating extracellular matrix remodeling and inflammatory cell migration.

Docking alone cannot guarantee that a ligand remains bound under physiological conditions, so the team followed up with a 100-nanosecond all-atom molecular dynamics simulation using GROMACS 2023 and the AMBER99SB-ILDN force field. The root mean square deviation values equilibrated within the first 10 nanoseconds and maintained a stable plateau thereafter, indicating a steady conformational state. Root mean square fluctuation analysis showed consistently low local motion across most residues, and the radius of gyration declined gradually over the first 60 nanoseconds before stabilizing around 1.46 nanometers, reflecting progressive compaction of the complex. The number of intermolecular hydrogen bonds remained steady at one to two throughout the trajectory, and comparison of snapshots at 0 and 100 nanoseconds showed no substantial dissociation or deformation of colchicine within the binding pocket. Together, these results support a durable and specific interaction between the drug and its newly identified target.

Beyond colchicine, the researchers explored a dual screening strategy for additional therapeutic candidates targeting the six hub genes. Among the top 10 predicted small molecules were simvastatin, which has established anti-inflammatory properties and has been reported to attenuate NLRP3 inflammasome activation; ciprofloxacin, a fluoroquinolone with immunomodulatory effects; and nimesulide, a selective COX-2 inhibitor that could theoretically ease gout flares, though its hepatotoxicity profile warrants caution. The team also reviewed traditional Chinese medicine formulations reported to act on the hub genes, identifying six prescriptions including Ma Xing Shi Gan Decoction, Yi-Fei-Jian-Pi-Tang, and Taodan Granules. Yi-Fei-Jian-Pi-Tang showed the broadest targeting spectrum, modulating five of the six hub genes, while Taodan Granules targeted three. The authors caution that several computational hits, such as nickel sulfate and electrocorundum, are environmental agents rather than drugs, and that all predictions require rigorous experimental validation before clinical translation.

Perhaps the most striking finding came from single-cell RNA sequencing analysis of peripheral blood mononuclear cells, drawing on datasets GSE188172 and GSE211783 and encompassing more than 94,000 cells across severe COVID-19 patients, gout patients, and healthy controls. After quality control and batch correction using Harmony integration, the team annotated eight major immune cell subpopulations and found that MMP1 expression was significantly upregulated in COVID-19 patients, with the elevation localized predominantly to B cells and plasma cells rather than myeloid compartments. This B-cell-restricted pattern was further accentuated in ventilated COVID-19 patients compared with non-ventilated cases, implying a positive correlation between MMP1 expression in B cells and disease severity. In gout patients, by contrast, MMP1 showed a more distributed expression across monocytes and macrophages, consistent with the innate immune-driven pathology of crystal-induced inflammation. The authors suggest that while MMP1 serves as a shared molecular hub at the transcriptomic level, its cellular sources and functional roles may be disease-context-dependent, and that therapeutic targeting may ultimately require cell-type-specific strategies. Taken together, the findings position MMP1 as a promising candidate for patients facing the double burden of gout and severe COVID-19, while underscoring the substantial gap that remains between computational prediction and clinical application.

Subject of Research: Shared molecular pathways between gout and severe COVID-19, focusing on MMP1 as a common hub gene and drug target

Article Title: MMP1 serves as a common link and drug target between gout and severe COVID-19 infection

Article References: Ma, H., Geng, Y., Liu, M., Feng, T., Lv, Y., Wang, H., & Niu, Y. (2026). MMP1 serves as a common link and drug target between gout and severe COVID-19 infection. Results in Chemistry, 31, Article 103967. https://doi.org/10.1016/j.rechem.2026.103967

Image Credits: AI Generated

DOI: 10.1016/j.rechem.2026.103967

Keywords: MMP1, gout, severe COVID-19, colchicine, molecular docking, molecular dynamics simulation, single-cell RNA sequencing, NLRP3 inflammasome, hub genes, bioinformatics, inflammation, drug repurposing

The research team, led by investigators at the Hospital for Special Surgery in New York in collaboration with colleagues at Tianjin First Central Hospital and the Performance Medicine Institute, focused on Yes-associated protein, or YAP, a transcriptional regulator that sits at the heart of the Hippo signaling pathway. YAP is famous in mechanobiology because its location within the cell is dictated by physical forces: when a cell is stretched, YAP moves into the nucleus and switches on genes that help the cell adapt to its mechanical environment. Although this behavior has been documented in many cell types, no previous study had systematically deleted YAP in mature tenocytes and then interrogated the resulting transcriptomic landscape under controlled mechanical loading. That gap mattered, because tenocytes in a living tendon are never at rest; they are continuously exposed to tensile strain, and any gene whose regulation by YAP becomes critical under those conditions is a gene worth knowing about.

To answer the question, the researchers built a genetically engineered mouse strain carrying a tamoxifen-activated Cre recombinase at the Rosa26 locus crossed with mice whose YAP gene was flanked by loxP sites. Tenocytes were harvested from the tail tendons of three-month-old male mice and grown in culture. Treating the cells with 4-hydroxytamoxifen excised the critical exons of YAP, and Western blotting confirmed that the protein had been effectively eliminated. The cells were then embedded in collagen gels in Flexcell Tissue Train plates and subjected to a cyclic tensile loading regimen of 4 percent strain at 1 hertz for 24 hours, a protocol designed to mimic the physiological stretching that tenocytes experience in a working tendon. After loading, RNA was extracted and sequenced on an Illumina NovaSeq X Plus platform, with reads aligned to the mouse genome and analyzed using the edgeR package to identify differentially expressed genes.

The headline result was stark. Comparing YAP-deficient tenocytes with controls under identical mechanical loading, the team identified 54 differentially expressed genes, of which only six were upregulated and 48 were downregulated. The overwhelming downward shift suggests that YAP acts predominantly as a transcriptional activator in mechanically stimulated tenocytes, and that its removal silences a large slice of the load-responsive program. Among the most strongly upregulated genes were Gzmc and Gm28151, while the most strongly downregulated included Car3, which encodes carbonic anhydrase 3, and Rgs2, a regulator of G-protein signaling. Both Car3 and Rgs2 have been linked to YAP signaling in other biological systems, which the authors note lends partial external consistency to their findings.

When the researchers ran the 54 genes through Gene Ontology and KEGG pathway enrichment analyses, a clear thematic signature emerged: lipid metabolism. The genes were significantly enriched in biological processes such as fat cell differentiation and fatty acid metabolism, in cellular components including lipid droplets and lipoprotein particles, and in molecular functions such as heparin and glycosaminoglycan binding. At the pathway level, the PPAR signaling pathway stood out, alongside regulation of lipolysis in adipocytes and the AMPK signaling pathway. Gene set enrichment analysis reinforced the metabolic theme, showing that YAP knockout cells were significantly enriched for hallmark gene sets covering adipogenesis, cholesterol homeostasis, and epithelial-mesenchymal transition, while control cells showed enrichment of E2F target genes associated with proliferation. The authors are careful to frame these observations as correlative: the experiments were performed in an in vitro culture system without injury or degeneration, so the findings generate hypotheses rather than establish mechanisms of tendon repair.

The lipid connection is biologically intriguing rather than incidental. Lipids serve as an energy source for tenocytes and help organize cell membranes, influencing the fluidity and lipid raft dynamics that modulate mechanosensitive channels and integrin-mediated force transmission. Lipid-derived mediators also participate in extracellular matrix remodeling by regulating collagen synthesis and the activity of matrix-degrading enzymes. PPARγ, the central regulator of adipogenic differentiation, has been shown to interact functionally with YAP in other cell types, and the new data raise the possibility that this crosstalk operates in tendon cells as well. Whether the metabolic transcriptional program observed here is a cause or a consequence of YAP deficiency, and whether it helps maintain the tenocyte phenotype under load, remains an open question that the authors say will require lipidomic analyses and targeted perturbation of PPAR signaling in future work.

To distill the 54-gene list down to the most consequential players, the team deployed two complementary machine learning approaches. LASSO regression identified seven candidate genes, and support vector machine recursive feature elimination identified 39; their intersection yielded seven hub genes: Car3, Cxcl13, Hvcn1, Ltc4s, Plin4, Rasd1, and Rgs2. Expression analysis confirmed significant differences between knockout and control groups for all seven, and correlation analysis revealed strong positive relationships among them, hinting that they operate as a coordinated module. The authors emphasize that these selections were derived from a limited number of biological replicates, three per group, which is the standard minimum for bulk RNA sequencing, and that validation in larger cohorts is essential before any general conclusions can be drawn.

One of the more provocative analyses involved CIBERSORTx, a computational method that estimates immune cell proportions from bulk gene expression data. Applied here as an exploratory check rather than a quantitative measurement, the analysis suggested significant differences between knockout and control samples in the signatures of regulatory T cells and CD4-positive naive T cells, and all seven hub genes showed positive correlations with the Treg signature. The authors are explicit that cultured tenocytes contain no immune cells and that these scores are indirect inferences about the immune landscape that might arise from tenocyte-immune crosstalk in a living tendon. Even so, the result dovetails with a growing appreciation that tendon healing depends on precisely orchestrated inflammatory responses, and it suggests that YAP-dependent paracrine signaling from tenocytes could help shape the recruitment and polarization of immune cells during repair.

The study also ventured into translational territory. Screening the Drug-Gene Interaction Database revealed that four of the hub genes, LTC4S, CA3, CXCL13, and RGS2, are associated with known drugs or small molecules, and twelve compounds emerged from the network analysis. Molecular docking simulations ranked candidate ligands against the LTC4S protein, with Compound 39 showing the strongest predicted binding energy at minus 7.78 kilocalories per mole, followed by TK05 at minus 6.24. The researchers stress that these are purely computational predictions; docking scores indicate geometric and energetic compatibility, not functional binding or therapeutic efficacy, and experimental validation would be required before any clinical relevance could be considered.

The authors are refreshingly candid about the limitations of their work. There is no in vivo validation, the differentially expressed genes have not yet been confirmed by qRT-PCR or Western blotting, the sample size is small, and the cultured tenocytes lack the three-dimensional extracellular matrix architecture, vasculature, and immune environment of native tendon tissue. The immune findings in particular remain indirect predictions rather than measurements. Still, the study delivers something the field has lacked: a systematic, mechanically contextualized map of what YAP does in mature tenocytes at the transcriptomic level. The data, deposited in the Gene Expression Omnibus under accession GSE346125, provide a hypothesis-generating resource for tenocyte mechanobiology. Future experiments using conditional YAP knockout mice in tendon injury models, combined with functional validation of genes such as Car3 and Rgs2 and targeted manipulation of PPAR and lipid signaling, will determine whether the molecular circuitry uncovered here matters for real tendon healing, and whether the immune and pharmacological angles it opens can be turned into therapies for the millions of people sidelined by tendon injuries each year.

Subject of Research: YAP-dependent mechanotransduction and transcriptional responses in tenocytes under cyclic tensile loading

Article Title: YAP-dependent transcriptional responses in tenocytes under cyclic tensile loading: a conditional knockout RNA-sequencing study

Article References: Xing, X., Piacentini, A., Suzuki, Y., Sanchez, L. A., Mendias, C., & Rodeo, S. A. (2026). YAP-dependent transcriptional responses in tenocytes under cyclic tensile loading: a conditional knockout RNA-sequencing study. Frontiers in Bioengineering and Biotechnology, 14, Article 1901935. https://doi.org/10.3389/fbioe.2026.1901935

Image Credits: AI Generated

DOI: 10.3389/fbioe.2026.1901935

Keywords: YAP, tenocytes, mechanotransduction, tendon injury, conditional knockout, RNA sequencing, PPAR signaling, lipid metabolism, CIBERSORTx, hub genes, molecular docking, tendon repair

Until this study, the genetic machinery of the porcine carpal gland was essentially terra incognita. The researchers note that only a single genome-wide association study had previously touched on the organ, and no transcriptomic investigation had ever been carried out. That gap is striking given the gland’s ecological and practical importance. Beyond its role in pig behavior, the gland’s secretions have even found value in wild boar population management, where chemical signaling can be exploited to monitor or influence free-ranging populations. Understanding which genes build and operate the gland could therefore matter not just for fundamental biology, but for livestock breeding and wildlife management alike.

To lay the groundwork, the team began with classical histology. Carpal glands were dissected from the forelegs of pigs and immediately fixed in 10 percent paraformaldehyde for two to three days, a step that preserves cellular architecture in lifelike detail. The fixed tissues were then embedded in paraffin, sliced into sections just five micrometers thick, and stained with hematoxylin and eosin, the workhorse dye pair that has served histologists for more than a century. Under the light microscope, the stained sections revealed the gland’s secretory architecture in fine detail, allowing the researchers to count secretory portions and secretory cells and to compare structures across animals.

The histological comparisons produced a clear and intriguing pattern. Both the visible phenotypes and the underlying tissue structures of the carpal glands differed significantly among the four populations studied: wild boars, Min pigs, Duroc pigs, and HS pigs. Wild boars and Min pigs, the more extensively characterized of the groups, possessed more carpal gland secretory portions and more secretory cells than their Duroc and HS counterparts. This gradient, running from wild and indigenous Chinese breeds toward commercially selected Western breeds, hints that intensive domestication and selection for production traits may have reshaped an organ whose primary business is chemical communication rather than meat yield.

With the anatomical baseline established, the researchers turned to the transcriptome. They sequenced the RNA of 32 carpal glands using RNA-seq, a technique that captures a snapshot of every gene being actively transcribed in the tissue at the moment of sampling. The sequencing data revealed a total of 5,205 differentially expressed genes, or DEGs, when the team compared populations with one another and males with females. That such a large fraction of the transcriptome varies across breeds and sexes underscores how dynamic the gland’s molecular program is, and it provided a rich dataset for the next stage of analysis.

Raw lists of differentially expressed genes can be overwhelming, so the team applied weighted gene co-expression network analysis, or WGCNA, a statistical framework that groups thousands of genes into modules based on the correlated patterns of their expression across samples. Genes that rise and fall together often participate in the same biological processes, so modules can act as functional units rather than isolated data points. The analysis identified two modules whose expression profiles showed the highest correlation coefficients with carpal gland-related traits, effectively narrowing thousands of genes down to a handful of coordinated gene networks most likely to be shaping the organ’s structure and output.

To interrogate those modules further, the researchers layered on additional computational tools. Gene set enrichment against Gene Ontology and KEGG databases flagged seven pathways associated with the carpal gland, including ‘gland development’ (GO:0048732) and ‘tube development’ (GO:0035295), terms that fit neatly with the organ’s role as a secretory skin structure whose ducts and acini must be built during development. The team also constructed protein-protein interaction networks, which map the physical and functional relationships among the proteins encoded by the module genes, and predicted transcription factors, the regulatory proteins that switch genes on and off, acting on the hub genes within the WGCNA modules.

The convergence of these analyses distilled the dataset to a remarkably compact candidate list. Ten hub genes emerged as plausible candidates for the porcine carpal gland: CCL2, GJA1, S100A4, PLAUR, LUM, ANXA1, BGN, BMP1, GPNMB, and SERPINE1. Several of these are familiar faces in secretory and structural biology. LUM and BGN encode small leucine-rich proteoglycans that scaffold the extracellular matrix, while BMP1 is a procollagen-processing enzyme central to tissue remodeling. GJA1 encodes connexin 43, a gap junction protein that allows secretory cells to communicate directly, and CCL2 and ANXA1 participate in signaling and inflammatory pathways that often intersect with glandular secretion. Alongside the ten hub genes, the team pinpointed two transcription factors, SNAI1 and SNAI2, as important candidate regulators of the gland.

Critical to the credibility of any genomics study is validation, and the researchers did not stop at computational prediction. They confirmed their key findings using quantitative real-time PCR, or qRT-PCR, an independent laboratory technique that measures the expression of specific genes with high sensitivity. They also performed quantitative trait-transcript analysis, known as QTT analysis, which links transcript abundance directly to measured quantitative traits across individuals, providing a statistical bridge between gene expression and the physical characteristics of the gland. This two-pronged validation strategy strengthens the case that the identified hub genes and transcription factors are genuinely tied to carpal gland biology rather than artifacts of a single analytical pipeline.

The study’s authors, led by Dengshuai Cui and corresponding author Yuanmei Guo at Jiangxi Agricultural University’s National Key Laboratory for Pig Genetic Improvement and Germplasm Innovation, describe their work as outlining the transcriptome profiling of the carpal gland across four populations and providing new insights into the genetic mechanism of the organ. The identification of two highly correlated co-expression modules, ten hub genes, and the SNAI1 and SNAI2 transcription factors gives researchers concrete molecular entry points for future experiments, whether those involve functional studies of individual genes, breeding programs that account for chemical signaling traits, or comparative work on scent glands in other mammals. For an organ that has spent centuries quietly shaping pig societies while escaping molecular scrutiny, the carpal gland has finally begun to reveal its genetic secrets, and the implications may ripple from the pigsty to wildlife management well beyond it.

Subject of Research: Transcriptomic analysis of the porcine carpal gland to identify genes and regulatory modules underlying its structure and pheromone secretion

Article Title: Identification of crucial modules and genes associated with carpal gland using RNA-seq and WGCNA in four pig (Sus scrofa) populations

Article References: Cui, D., Xiong, S., Xie, L., Yu, N., Wei, J., He, R., Xiao, S., Li, L., & Guo, Y. (2026). Identification of crucial modules and genes associated with carpal gland using RNA-seq and WGCNA in four pig (Sus scrofa) populations. BMC Genomics. https://doi.org/10.1186/s12864-026-13445-9

Image Credits: AI Generated

DOI: 10.1186/s12864-026-13445-9

Keywords: pig, carpal gland, RNA-seq, WGCNA, transcriptomics, hub genes, transcription factors, pheromones, Sus scrofa, differentially expressed genes, histology, BMC Genomics

The study, published in the open-access journal Stress Biology, focused on a commercially available fungicide formulation containing two beneficial fungal strains, Trichoderma asperellum ICC012 and Trichoderma gamsii ICC080. Trichoderma species are filamentous fungi famous among plant scientists for their mycoparasitic and antagonistic activities against a wide range of plant pathogens, and they are already registered and widely used as microbial fungicides. What has remained poorly understood, however, is precisely what happens at the level of gene expression inside host roots when these beneficial fungi are deployed ahead of a pathogen attack. The new research set out to close that gap by combining classical disease screening with high-throughput RNA sequencing and an advanced statistical framework known as Weighted Gene Co-expression Network Analysis, or WGCNA.

The phenotypic results were striking. In greenhouse trials using the Greek melon variety ‘Xrisi-kefali’, seedlings inoculated with FOM alone developed the full suite of Fusarium wilt symptoms, and by the end of the observation period the disease incidence had reached one hundred percent, with roughly 79 percent of leaves lost to defoliation. Seedlings pre-treated with the Trichoderma formulation before pathogen inoculation told a very different story. Disease incidence in this group was limited to about 40 percent, defoliation dropped to approximately 28 percent, and weekly disease index assessments over four consecutive weeks showed significantly lower severity compared with pathogen-only plants. Beyond disease suppression, the treated seedlings also grew better, showing greater shoot height, thicker stems, and more internodes, consistent with the well-documented growth-promoting activity of Trichoderma in the rhizosphere, which is thought to involve improved nutrient availability, mineral solubilization, iron chelation, and the modulation of hormonal signaling pathways.

To understand the molecular machinery behind this protection, the researchers sequenced the transcriptomes of melon roots collected at one, two, and three days after pathogen inoculation, comparing three treatments: plants treated with Trichoderma and then challenged with FOM, plants challenged with the pathogen alone, and untreated healthy controls. The sequencing effort generated approximately 944 million clean reads, averaging around 34 million per sample, with more than 88 percent mapping successfully to the melon reference genome. After stringent filtering, the team identified a non-redundant set of 9254 differentially expressed genes across all comparisons. Notably, the highest number of upregulated genes appeared in Trichoderma-pretreated, pathogen-challenged roots at the first and third days after inoculation, suggesting that the beneficial fungi trigger a time-dependent surge of defensive gene activity that pathogen infection alone does not produce.

The centerpiece of the study was the network analysis. Using WGCNA, the researchers grouped the 9254 genes into co-expression modules and eventually refined them into fourteen modules, four of which emerged as decisively correlated with the treatments. Two modules, dubbed midnightblue and cyan, were most strongly associated with pathogen infection alone at two and three days after inoculation, respectively. Two others, the blue and green modules, were tightly linked to the Trichoderma-pretreated plants at the first and third days. This separation proved conceptually important: the modules tied to pathogen infection alone reflected a stressed, pathogen-driven transcriptional state, while the modules tied to the combined treatment captured the active, Trichoderma-primed defense machinery that appears to establish resistance before substantial pathogen ingress.

Within the blue module, which represents the early Trichoderma-induced response, the researchers found a remarkable concentration of defense-related hub genes, the highly connected nodes that act as central regulators of the network. Four transcription factor families stood out: NAC domain-containing protein 2, the dehydration-responsive element-binding protein DREB1A, trihelix transcription factor GT-3b, and a heat stress transcription factor. These families are known orchestrators of plant immunity in other crops; a trihelix GT-3b gene has been linked to defense against Fusarium graminearum in maize, a pathogen-induced NAC factor modulates rust resistance in barley, and DREB1A overexpression enhances biotic stress tolerance in transgenic potato. Alongside these regulators, the module contained a gene encoding caffeoylshikimate esterase, or CSE, which participates in lignin biosynthesis, a branch of the phenylpropanoid pathway that plants use to reinforce cell walls and produce antimicrobial phytoalexins against fungal invaders. KEGG pathway enrichment confirmed the module’s defensive character, highlighting phenylpropanoid biosynthesis, the MAPK signaling cascade, and plant-pathogen interaction pathways.

The blue module also harbored genes encoding a chitin receptor kinase called LYK5, which works with CERK1 to detect fungal cell wall fragments and launch immune signaling, a phloem protein 2 involved in phloem-based defense, a nerolidol synthase-like gene producing volatile defense signals, and a subtilisin-like protease implicated in immune responses. The green module, correlated with the later Trichoderma-induced response, contributed its own arsenal: a two-component response regulator ARR14-like gene involved in hormone crosstalk, Kunitz trypsin inhibitor and miraculin-like proteins with documented antifungal activity, a MYB27 transcription factor, translation-related factors, and even an aquaporin that may influence defensive signaling. Together, these findings sketch a coherent picture of Trichoderma priming the melon root to mount a faster, stronger, and more coordinated defense than it could achieve on its own.

The dual nature of the protection became even clearer when the researchers examined the pathogen’s own transcripts within the same sequencing data. Across all three time points, FOM gene expression was consistently and dramatically suppressed in the presence of the Trichoderma formulation, often approaching near-zero levels. The repressed genes included multiple ribosomal proteins, pointing to a collapse in the pathogen’s protein synthesis capacity, as well as key metabolic enzymes such as ATP synthase, malate dehydrogenase, and cytochrome c peroxidase, indicating disruption of respiration and energy production. Most tellingly, the virulence arsenal itself was silenced: genes encoding xyloglucanases, pectate lyases, endoglucanases, and exoglucanases, the carbohydrate-degrading enzymes FOM needs to pierce plant cell walls, were markedly reduced, particularly at later time points. While the authors note that formal differential expression analysis would be needed to statistically confirm these patterns, the profiles strongly suggest that Trichoderma attacks the pathogen on two fronts simultaneously, starving it at a fundamental physiological level while arming the plant against it.

Several intriguing regulatory players emerged from the pathogen-associated modules as well. The midnightblue module’s hub genes were predominantly ribosomal RNA genes, upregulated after two days of FOM infection and inverted under Trichoderma pretreatment, hinting at a pathogen-driven reprogramming of the host’s translational machinery that the beneficial fungi can counteract. The cyan module contained upregulated long non-coding RNAs, increasingly recognized as key modulators of plant immunity through their effects on reactive oxygen species, calcium signaling, and transcription factor regulation, plus two uncharacterized genes homologous to retrotransposon proteins, suggesting that transposable elements may participate in the stress response. A pyridoxal kinase gene involved in vitamin B6 biosynthesis, an essential antioxidant in stress responses, was also highly expressed in this module.

The authors are careful to frame their findings within both the promise and the caveats of biological control. The experimental design mirrored real-world practice, applying Trichoderma preventively so it can establish in the rhizosphere before pathogen challenge, and the priming-based mechanism of induced resistance depends on that prior colonization. They also acknowledge that microbial biocontrol agents can reshape soil fungal communities, sometimes reducing fungal diversity or shifting competitive equilibria in ways that depend on soil type and inoculant persistence, a knowledge gap that matters for durable, biodiversity-conscious disease management. Still, the study delivers a molecular roadmap that could inform targeted breeding strategies and refined industrial formulations. By pinpointing specific hub genes and pathways, from CSE-driven lignin reinforcement to NAC and DREB transcriptional cascades, the work transforms Trichoderma biocontrol from an empirical black box into a mechanistically understood, engineerable system for protecting one of the world’s favorite fruits.

Subject of Research: Transcriptomic mechanisms of Trichoderma-mediated biological control of Fusarium wilt in melon roots

Article Title: Transcriptomic insights into the biocontrol mechanism of Trichoderma spp. against Fusarium wilt in melon

Article References: Aci, M. M., Tsalgatidou, P. C., Krommydas, K., Boutsika, A., Delis, C., Pavli, O. I., Schena, L., & Zambounis, A. (2026). Transcriptomic insights into the biocontrol mechanism of Trichoderma spp. against Fusarium wilt in melon. Stress Biology, 6(1), Article 44. https://doi.org/10.1007/s44154-026-00315-3

Image Credits: AI Generated

DOI: 10.1007/s44154-026-00315-3

Keywords: Trichoderma, Fusarium wilt, melon, biocontrol, transcriptomics, RNA-seq, WGCNA, plant immunity, hub genes, phenylpropanoid biosynthesis, MAPK signaling, sustainable agriculture

The clinical backdrop to the work is well established. Epidemiological studies have repeatedly shown that patients with Parkinson’s disease face a strikingly elevated risk of osteoporosis and osteoporotic fractures, a burden that exceeds what can be explained by poor mobility, falls, or age alone. Meta-analyses cited by the researchers indicate high rates of osteoporotic fracture in Parkinson’s disease, and cross-sectional clinical work has associated biomarkers such as serum uric acid with reduced bone mineral density in affected patients. Yet the molecular mechanisms that bind neurodegeneration to bone loss have remained murky. The Henan University team set out to interrogate that relationship experimentally, asking whether a single genetic factor—absence of ACE2—could simultaneously perturb dopaminergic neuron survival in the brain and bone homeostasis in the skeleton.

To do so, the researchers used a well-characterized mouse model of parkinsonism in which the neurotoxin MPTP, or 1-methyl-4-phenyl-1,2,3,6-tetrahydropyridine, selectively destroys dopamine-producing neurons in the substantia nigra. Crucially, they crossed this challenge with mice lacking a functional ACE2 gene—specifically Ace2-null males, designated Ace2−/y—because ACE2 sits at the center of the protective arm of the renin-angiotensin system, converting angiotensin II into angiotensin-(1–7), a peptide with documented anti-inflammatory and neuroprotective actions. Previous studies from other groups and from the same team had shown that ACE2 activation mitigates behavioral deficits and neuroinflammation in chemically induced Parkinson’s models, and that the ACE2/Ang-(1–7)/Mas cascade strengthens bone structure and metabolism. The new study asked what happens when this protective factor is removed entirely under parkinsonian stress.

Behavioral testing revealed a clear aggravating effect. MPTP exposure significantly worsened motor dysfunction and depression-like behaviors in the mice, and the combination of MPTP toxicity with ACE2 deficiency produced a particularly severe pathological picture in the brain. Immunohistochemistry, Western blotting, and histopathological staining showed reduced activity of dopaminergic neurons and heightened microglial activation—the inflammatory response of the brain’s resident immune cells. At the molecular level, the researchers measured elevated levels of total α-synuclein, the misfolding-prone protein that defines Parkinson’s pathology, alongside increased abundance of Caspase-3 and Bax, two canonical executioners of programmed cell death. Together, these markers indicate that ACE2 loss intensifies both the protein aggregation burden and the apoptotic pressure on vulnerable neurons.

The bone findings were equally striking, and notably they emerged in parallel rather than secondarily. In the skeletal tissue of the Ace2-deficient mice, the team documented reduced abundance of total Wnt ligands, β-catenin, bone morphogenetic proteins (BMP), and insulin-like growth factor 1 (IGF-1), along with diminished phosphorylation ratios of the downstream kinases that transmit these signals. This matters because each of these cascades is a cornerstone of bone formation: Wnt/β-catenin signaling drives osteoblast differentiation and bone accrual, BMPs are potent inducers of bone formation used clinically in spine fusion and fracture repair, and IGF-1 couples muscle and bone metabolism through mTOR-dependent pathways. The researchers are careful to note that the parallel reduction of these signaling proteins suggests potential perturbation of the cascades rather than definitive proof of pathway failure, a distinction that reflects appropriate scientific caution.

In the opposite direction, ACE2 deficiency upregulated mediators of the RANKL/RANK/OPG axis, a triad that governs osteoclast formation and bone resorption. RANKL binding to RANK on osteoclast precursors drives the differentiation of bone-resorbing cells, while OPG acts as a soluble decoy receptor that restrains the process. Dysregulation of this axis tilts bone turnover toward net loss. Intriguingly, the same axis operates in the brain, where it has been identified as a critical inflammatory signaling system in ischemic injury, and Rho GTPases downstream of these pathways modulate osteoclast differentiation directly. The coordinated shift of this inflammatory skeletal axis in both brain and bone tissue under ACE2 deficiency is one of the study’s most suggestive observations, hinting at a common pathological language spoken by the two organs.

To move from candidate pathways to gene-level targets, the team turned to transcriptomics. They mined public GEO datasets and applied weighted gene co-expression network analysis, or WGCNA, a computational method that groups genes into modules based on correlated expression patterns and identifies the hub genes most central to disease-associated modules. This analysis pinpointed ten hub genes, including DNM1, which encodes dynamin 1, a protein essential for synaptic vesicle recycling; OCRL, a phosphatidylinositol phosphate phosphatase linked to the oculocerebrorenal syndrome; and OPA1, a mitochondrial fusion protein whose mutations cause dominant optic atrophy and which has been implicated in mitochondrial parkinsonism through stem cell modeling. The dysregulation of these genes was linked to synaptic dysfunction and inflammation—two processes squarely at the heart of Parkinson’s pathophysiology.

The team then evaluated whether these hub genes could serve as diagnostic biomarkers. Using receiver operating characteristic, or ROC, analysis on public single-disease transcriptome datasets for Parkinson’s disease and osteoporosis separately, they found that the core gene signatures achieved areas under the curve ranging from 0.683 to 0.981, indicating diagnostic accuracy that spans moderate to near-perfect discrimination. Functional enrichment of the core genes pointed to involvement in synaptic signaling, MAPK signaling, and the Rap1 and Ras pathways—small GTPase cascades that regulate cell proliferation, differentiation, and cytoskeletal dynamics in both neurons and bone cells. Such dual-diagnostic performance, if replicated in human cohorts, would suggest that a shared molecular signature underlies both conditions and could be exploited clinically to identify patients at risk of combined neurodegenerative and skeletal decline.

The authors are appropriately measured in their claims. They emphasize that these are preclinical findings obtained under short-term MPTP treatment in growing young male mice, meaning that the observed bone metabolic disturbance was transient and that the results may not translate directly to elderly human patients, in whom Parkinson’s disease typically manifests and in whom bone loss is chronic and sex-dependent. They also stress that the coordinated dysregulation observed in brain and bone is consistent with a bone-brain axis pathological phenotype, but that the current experimental design cannot confirm causal bidirectional cross-talk between the tissues. Distinguishing whether ACE2 deficiency independently damages both organs, or whether pathology in one propagates to the other—perhaps through circulating osteocalcin, sympathetic nervous system output, or inflammatory mediators—will require interventional studies that manipulate one tissue and measure the other.

Even with those caveats, the study offers a compelling framework and a set of concrete targets for follow-up. Restoring ACE2 activity or mimicking its product, angiotensin-(1–7), has already shown neuroprotective effects in experimental Parkinson’s models, including reduced α-synuclein expression through the NEAT1/miR-153-3p axis, and ACE2 activation has been reported to promote hippocampal neurogenesis via Wnt/β-catenin signaling. The present results raise the possibility that such therapies could carry a skeletal benefit as well, protecting against the osteoporosis that so often compounds the disability of Parkinson’s disease. The ten hub genes, meanwhile, provide a molecular shortlist for mechanistic validation, and their diagnostic AUC values justify testing in human blood or tissue datasets. As the bone-brain axis matures from a descriptive concept into a mechanistic research program, work like this demonstrates how a single molecule, studied across two organs at once, can illuminate disease connections that medicine has long observed clinically but struggled to explain at the level of genes and signaling pathways.

Subject of Research: ACE2-dependent molecular mechanisms linking Parkinson's disease neurodegeneration and bone metabolic alterations via the bone-brain axis

Article Title: ACE2 and Parkinsonism‑related bone metabolic alterations: signaling pathways and hub gene analysis

Article References: Liu, T., Ren, Y., Tian, X., & Wei, J. (2026). ACE2 and Parkinsonism‑related bone metabolic alterations: signaling pathways and hub gene analysis. Molecular Genetics and Genomics, 301(1), Article 185. https://doi.org/10.1007/s00438-026-02511-2

Image Credits: AI Generated

DOI: 10.1007/s00438-026-02511-2

Keywords: Parkinson's disease, ACE2, bone metabolism, bone-brain axis, osteoporosis, Wnt/beta-catenin signaling, RANKL/RANK/OPG, IGF-1, hub genes, WGCNA, alpha-synuclein, MPTP mouse model