SEARCH · PubMed Health

Results for “co-expression”

Explore indexed PubMed citations for clinical trials, systematic reviews and public health research. Read source abstracts and follow each citation to its original PubMed record.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 37 records · Page 2Linked to original sources

PubMed2026

Uncovering molecular regulatory networks of low-temperature stress response in Trachinotus ovatus via integrated transcriptome and metabolome analyses.

Golden pompano (Trachinotus ovatus) is one of the most economically important marine fish species in China. It is susceptible to low-temperature stress, which significantly challenges its production and supply. Nevertheless, study on the regulatory mechanisms underlying low-temperature stress responses in golden pompano remains limited. Here, we firstly performed a time-series transcriptome analysis to reconstruct dynamic response patterns under low-temperature stress in golden pompano. Transcriptome profiling identified common differentially expressed genes (DEGs), including fos, hlf, and hmgb1, as well as condition-specific DEGs across distinct low-temperature stress groups. Based on cluster analysis, all DEGs were classified into five distinct expression patterns, reflecting diversified regulation of expression in golden pompano during low-temperature stress. Furthermore, condition-specific regulatory modules were explored via weighted gene co-expression network analysis (WGCNA), highlighting that the two module hub genes, serbf2 and lipc, might respond to low-temperature stress by regulating the lipid catabolic process. Subsequently, untargeted metabolomic analysis revealed that glycerophospholipid metabolism was a significantly enriched common pathway, highlighting its crucial role in mediating the response to low-temperature stress. Finally, by integrating transcriptomic and metabolomic analyses, a gene-metabolite interaction network associated with glycerophospholipid metabolism under low-temperature stress was established. These findings underscore the significance of multiple candidate genes and glycerophospholipid metabolism in golden pompano's response to low-temperature stress, thereby laying a solid molecular foundation for the development of low-temperature-tolerant fish strains.

Animals↗

PubMed2026

Sex- and development-specific transcriptomic profiling of venom and silk genes in the wolf spider Pardosa astrigera provides insights into ecological adaptation and predatory strategies.

Spider venom and silk glands are two major secretory systems that contribute to prey capture, defense, and reproduction, but their sex- and development-specific molecular regulation in wandering wolf spiders remains poorly understood. Here, the transcriptome of Pardosa astrigera, an important agricultural natural enemy in China, revealed significant sex- and development-associated molecular differentiation among adult females, adult males, and spiderlings. A total of 100,025 unigenes were obtained, of which 23,852 were functionally annotated, providing a comprehensive transcriptomic resource for this species. Differential expression patterns showed marked variation among groups, with 531, 1792, and 832 DEGs detected in PAF vs PAS, PAM vs PAS, and PAF vs PAM, respectively. These genes were mainly associated with metabolic, oxidation-reduction, cuticle development, MAPK signaling, and lysosome pathways. Fifteen co-expression modules revealed distinct expression patterns. The turquoise, pink, yellow, and red modules were development-related, whereas the blue module was male-biased. Venom- and spidroin-related genes were distributed across multiple modules, suggesting coordinated regulation. Overall, 42 venom peptides, 21 venom proteins, and 11 spidroins were identified. Representative genes showed strongly biased expression, including spiderling-biased U3_Pp1a and U5_Pp1e, female-biased U4_Pp1a, and male-biased SMase D_108750 and PaTuSp_108466. These findings reveal sex- and development-biased expression patterns of venom- and silk-related candidate genes in P. astrigera and may provide molecular insights into ecological adaptation and predatory strategies in wandering wolf spiders.

Animals↗

PubMed2026

Transcriptomics reveals species-specific adaptive strategies to calorie restriction in two Argopecten scallops with distinct lifespans.

Calorie restriction (CR) is a well-established non-genetic intervention for lifespan extension in multiple model organisms. Seasonal food shortage in cold and temperate seas may mimic CR, inducing in bivalves a response similar to that in vertebrates and thereby prolonging life expectancy. However, the relationship and the mechanism underlying the food availability and lifespan in bivalves remain largely unexplored. Two closely related scallop species the short-lived warm-water Argopecten irradians (lifespan <2 years) and the longer-lived cold-water Argopecten purpuratus (7-10 years) provide an ideal comparative system to investigate species-specific adaptive strategies. In this study, we subjected both species to CR for 30 and 56 days and performed comparative transcriptomic profiling, weighted gene co-expression network analysis (WGCNA), and physiological assays to elucidate their distinct molecular responses. Transcriptomic analysis revealed that A. purpuratus exhibited substantially more DEGs than A. irradians at both time points under CR, with both species showing downregulation of metabolic pathways but to different extents. A. irradians mounted an early nutrient-sensing response at 30 days (IGF1R, PIK3R3, INSR suppression), indicating acute sensitivity to limitation; by contrast, A. purpuratus displayed delayed FoxO activation at 56 days, along with its downstream effectors NFKBIA, CREB3L4, and SMAD4, suggesting a gradual adaptive program may link to its extended lifespan. WGCNA identified three negatively correlated modules in each species, with coral2 being the most prominent in A. irradians and darkolivegreen in A. purpuratus. The former was dominated by ciliary motility genes, whereas the latter featured coordinated repression of oxidative phosphorylation. Additionally, both species exhibited conserved suppression of mTOR/S6K growth signaling and activation of cellular maintenance programs. Collectively, these findings expand the understanding of CR-mediated longevity regulation in bivalves and provide candidate gene resources for future functional studies and breeding programs.

Pectinidae↗

PubMed2026

Multi-omics analysis reveals coordinated epigenetic dysregulation in atrazine-induced dopaminergic neurotoxicity.

Atrazine (ATR), a widely used triazine herbicide, has been linked to neurotoxicity, yet the epigenetic mechanisms underlying its dopaminergic effects remain unclear. This study investigated whether coordinated miRNA dysregulation and DNA methylation alterations contribute to ATR-induced Parkinson's disease (PD)-like neurotoxicity. Male Sprague-Dawley rats were administered ATR (50 mg/kg/day) for 90 days, resulting in motor and cognitive deficits with dopaminergic dysfunction, including increased α-synuclein and reduced tyrosine hydroxylase expression. Small RNA sequencing identified 72 differentially expressed miRNAs in the substantia nigra, enriched in PI3K-Akt, MAPK, and Ras signaling pathways. In a cohort of six PD patients and six matched controls, genome-wide DNA methylation profiling revealed 4694 differentially methylated positions, predominantly hypomethylated, with overlapping enrichment in neuronal signaling pathways. Weighted gene co-expression network analysis identified a PD-associated module strongly correlated with disease status (r = -0.95, P < 0.001). Multi-omics integration identified CASP3 as a central hub gene. External validation supported CASP3 relevance in PD (AUC = 0.833), and molecular docking suggested potential ATR-CASP3 interaction. Further analysis predicted upregulated miR-3552 as a potential upstream regulator of CASP3. These findings indicate that ATR-induced neurotoxicity may be mediated through the miR-3552/CASP3 signaling axis, ultimately regulating apoptosis and contributing to neurodegeneration.

Animals↗

PubMed2026

An oxidative stress - and immunotherapy-related six-gene signature defines immune subtypes and predicts prognosis and immunotherapy response in hepatocellular carcinoma.

BACKGROUND: Oxidative stress and the tumor immune microenvironment jointly shape hepatocellular carcinoma (HCC) progression and response to immunotherapy, yet integrated biomarkers linking these processes are lacking. METHODS: Transcriptomic and clinical data from The Cancer Genome Atlas (TCGA) and Gene Expression Omnibus (GEO) datasets were used to identify oxidative stress- and immunotherapyrelated differentially expressed genes (OSIRDEGs). Functional enrichment, weighted gene co-expression network analysis (WGCNA) and LASSO-Cox regression were used to construct a prognostic signature. Consensus clustering, TIDE, CIBERSORT and ssGSEA characterized immune phenotypes. Somatic mutation, copy-number and drug-response data were integrated to assess genomic alterations and drug sensitivity. Expression of model genes was validated by qRT-PCR and western blotting in HCC cell lines. RESULTS: We identified 24 OSIRDEGs enriched in cell-cycle and mitotic pathways. WGCNA intersection yielded 18 module genes, from which a six-gene signature (BUB1B, CDKN2A, CENPE, HMMR, PTTG1, SPP1) was derived. The signature robustly stratified patients into high- and low-risk groups with significantly different progression-free and disease-free survival in both TCGA-LIHC and GSE14520. Based on signature expression, two molecular subtypes were defined, exhibiting distinct survival, immune landscapes and predicted immunotherapy responsiveness. Model genes harbored recurrent alterations and showed significant correlations with anticancer agents. All six genes were upregulated at mRNA and protein levels in metastatic HCC cell lines versus normal hepatocytes. CONCLUSIONS: We systematically explored the landscape of OSIRDEGs in HCC, and proposed a validated six-gene signature that refines prognostic stratification, delineates immunerelevant HCC subtypes and highlights candidate biomarkers for therapeutic selection and mechanistic investigation.

Humans↗

PubMed2026

Single-allele chromatin tracing reveals cytokine-dependent super-enhancer repositioning in CD4+ T cells.

Naive CD4+ T cells interpret cytokine cues to commit to T helper lineages. Here, we examined the impact of cytokines on the Ets1-Fli1 locus, which encodes paralogous transcription factors essential for T cell effector responses. Epigenomic and chromatin interaction profiling in double-positive (DP) thymocytes defined a T cell super-enhancer and a CTCF-bound boundary. Chromatin tracing at single-allele resolution revealed concurrent multi-way interactions among Ets1, Fli1, and the super-enhancer coupled to transcription. Deleting the CTCF boundary decompacted the locus without altering Ets1 expression or T cell development, whereas deleting the super-enhancer drew Ets1 and Fli1 closer; promoter proximity increased Ets1-Fli1 co-expression. In CD4+ Th1 cells, cytokines increased multi-way interactions and repositioned the super-enhancer toward the geometric center, activating both paralogs. Super-enhancer deletion rendered Th1 cells structurally and transcriptionally DP like. Thus, cytokines can drive lineage-specific gene activation by repositioning a super-enhancer, providing a mechanistic framework for how noncoding variants at the Ets1-Fli1 locus contribute to CD4+ T cell-mediated immune disorders.

Animals↗

PubMed2026

Synovial short-lived plasma cells mediate adalimumab resistance in rheumatoid arthritis via MIF-CD74 axis-driven, partially TNF-α-independent inflammation.

OBJECTIVE: Synovial plasma cell infiltration predicts inadequate response to adalimumab in patients with rheumatoid arthritis (RA), yet the cellular and molecular mechanisms underlying this association remain unclear. This study aimed to dissect the functional heterogeneity of synovial plasma cells between adalimumab responders and non-responders at single-cell resolution, and to identify the molecular pathways driving treatment resistance. METHODS: This study was based on a prospective clinical cohort of 101 RA patients receiving adalimumab, from which synovial tissues of 8 patients (4 ACR20 responders and 4 non-responders) were profiled by 10x Genomics single-cell RNA sequencing (66,539 high-quality cells). A systematic ligand-receptor screening was performed to identify candidate signaling axes. Core findings were validated at four levels: an independent single-cell validation cohort (n = 4), external bulk RNA-seq cohorts (GSE15602, GSE47726), multiplex immunofluorescence on synovial tissues (n = 9 per group), and in vitro functional experiments using patient-derived peripheral blood monocyte-derived macrophages stimulated with recombinant human MIF under pharmacological intervention with adalimumab, the MIF inhibitor ISO-1, and an anti-CD74 neutralizing antibody. RESULTS: Plasma cells were significantly enriched in non-responder synovium, with a heterogeneous pattern characterized by quantitative accumulation of long-lived plasma cells (LLPCs) and functional dominance of short-lived plasma cells (SLPCs): SLPCs contributed 58.15% of total ribosomal module activity and preferentially overexpressed MIF. Systematic screening of 145 candidate ligand-receptor pairs identified MIF-CD74 as the only axis satisfying all four independent evidence layers. Tissue-level immunofluorescence confirmed that approximately 95% of synovial CD138+ plasma cells in non-responders co-expressed MIF, compared with approximately 45% in responders. In vitro, rh-MIF upregulated macrophage activation markers (CD74, CD80, CD86, HLA-DR) and induced IL-6 and TNF-α secretion. Adalimumab neutralized supernatant TNF-α but failed to suppress MIF-driven IL-6 and IL-1β activation, whereas ISO-1 and anti-CD74 effectively blocked MIF-induced effects at all levels examined. These findings were replicated in patient-derived PBMC macrophages. CONCLUSION: In adalimumab-resistant RA, a functionally active SLPC subset drives partially TNF-α-independent macrophage inflammation through the MIF-CD74 axis, representing a resistance pathway not fully addressed by anti-TNF therapy. Targeting MIF or CD74 blocked this axis in vitro, supporting MIF-CD74-directed precision intervention.

Adalimumab↗

PubMed2026

Temporal DIA-MS proteomics reveals coordinated metabolic reprogramming associated with oil accumulation in oil palm mesocarp.

Oil palm (Elaeis guineensis Jacq.) is the most productive oil-bearing crop globally, yet the molecular basis of mesocarp development and lipid accumulation remains poorly understood. Ultra-deep data-independent acquisition mass spectrometry (DIA-MS) was applied to characterize proteome dynamics in two contrasting genotypes, seedless (KS) and thin-shelled (TS), across five developmental stages (P1-P5) spanning fruit development to mature oil accumulation. Phenotypic analysis revealed higher mesocarp proportion and oil content in KS during late maturation. A total of 137,615 peptides corresponding to 12,163 protein groups were identified, providing a temporal proteomic landscape of mesocarp development. Multivariate analysis indicated that developmental progression was the primary contributor to proteomic variation, whereas genotype-associated differences increased during lipid accumulation. Differentially abundant proteins were mainly associated with carbohydrate metabolism, photosynthesis, proteolysis, antioxidant responses, and lipid biosynthesis. Gene Ontology (GO), Kyoto Encyclopedia of Genes and Genomes (KEGG), and KOG analyses suggested extensive remodeling of metabolic networks, including developmental changes in photosynthesis-associated proteins and increased representation of lipid-associated pathways during maturation. Weighted protein co-expression network analysis identified 17 modules associated with developmental progression and lipid accumulation, highlighting candidate proteins involved in carbon metabolism, energy production, and cellular protection. Genes encoding selected hub protein candidates were further examined by RT-qPCR. Biochemical analyses supported these proteomic patterns, showing increased acetyl-CoA availability, enhanced antioxidant enzyme activities (SOD, CAT, APX, and GR), improved GSH/GSSG balance, and reduced oxidative damage in KS. Together, these findings provide a temporal proteomic and biochemical framework for understanding genotype-associated differences in oil accumulation and identify candidate metabolic networks for functional studies.

Carbon metabolism↗

PubMed2026

Dynamic evolution of chaperone-mediated autophagy is associated with tumor microenvironment remodeling and prognostic stratification in lung adenocarcinoma: insights from single-cell transcriptomics, ensemble machine learning, and experimental validation.

BACKGROUND: Lung adenocarcinoma (LUAD) shows prognostic heterogeneity, and tumor-node-metastasis (TNM) staging is limited for individualized management. Chaperone-mediated autophagy (CMA) maintains proteostasis, but its role during adenocarcinoma in situ (AIS)-minimally invasive adenocarcinoma (MIA)-invasive adenocarcinoma (IAC) progression remains unclear. METHODS: Single-cell RNA sequencing (scRNA-seq) data from GSE189357 and bulk transcriptomes from The Cancer Genome Atlas (TCGA)-LUAD and Gene Expression Omnibus (GEO) cohorts were integrated. CMA activity, cell-cell communication, weighted gene co-expression network analysis (WGCNA), tumor-normal differential expression, machine-learning survival modeling, tumor microenvironment (TME) features, drug sensitivity, and EPC1 function were analyzed. RESULTS: CMA-high tumor epithelial cells increased from AIS (58.1%) to MIA (65.7%) but declined in IAC (44.4%; p < 0.001). CMA-low cells preferentially received fibroblast-derived extracellular matrix cues. A CMA-negatively correlated module identified 69 core genes. Random survival forest (RSF) performed best among 117 machine-learning combinations (mean concordance index > 0.873). High-risk patients had worse survival across cohorts, and the risk score was independently associated with overall survival (hazard ratio = 16.013, 95% confidence interval: 9.579-26.768, p < 0.001). High-risk tumors showed proliferative activation and M0 macrophage enrichment, whereas low-risk tumors showed stronger immune-related signaling. EPC1 overexpression suppressed malignant phenotypes in A549 cells. CONCLUSION: CMA dynamics are associated with stromal and immune remodeling during LUAD progression. A CMA-based model provides robust prognostic stratification and may offer a basis for future TME-guided studies.

Chaperone-mediated autophagy↗

PubMed2026

Single-cell spatial transcriptomic atlas of the mouse adrenal gland reveals sexual dimorphism in steroidogenic enzyme and hormone receptor expression.

The adrenal cortex shows sexual dimorphism in structure and function. We analysed adrenal glands from 7-week-old male and female BALB/c mice using Visium HD with Cellpose 3 segmentation, comprising 236,077 cells across eleven populations, including four cortical zones. Using curated marker-gene-based zonal annotation, we focused on steroidogenic enzymes and hormone receptors, complementing our companion study based on the same primary dataset. The X-zone was nearly absent in males but prominent in females. Females showed higher Hsd3b1 expression across cortical zones and higher Cyp11b1 expression in outer cortical compartments. The strongest sex difference involved Srd5a2, with markedly higher expression in male zona fasciculata (inner: 77.1% vs. 28.9%), independently supported by RNAscope and immunohistochemistry. Mc2r and Mrap showed discordant spatial distributions, with limited co-expression, suggesting potential MC2R-independent MRAP roles. Agtr1a dominated angiotensin II receptor expression in zona glomerulosa without major sex differences, providing a zone-resolved reference for adrenal sexual dimorphism.

Adrenal cortex↗

PubMed2026

A dual-dimensional CRISPR toolkit enables one-step high-efficiency multiplex genome editing in Komagataella phaffii.

Against the backdrop of green biomanufacturing, engineering methanol-utilizing Komagataella phaffii (K. phaffii) represents an effective strategy to expand the one carbon (C1) product profile and speed up the industrialization of C1-based bioeconomy. To address the technical challenges of low efficiency and cumbersome experimental procedures for multiplex gene editing and precise large-fragment integration during the reconstruction of complex metabolic pathways in K. phaffii, this study established a CRISPR toolkit - Efficient Multi-Gene Editing System 3.0 (EMGES 3.0) - which enabled one-step large-fragment integration coupled with multiplex gene knockout. EMGES 3.0 was constructed through the synergistic optimization of a repair-engineered chassis and an episomal CRISPR vector. For chassis engineering, five DNA repair modules: Δlig4 (DNA Ligase IV, non-homologous end joining end ligation), ppMRE11(The endogenous MRE11 gene from Pichia pastoris) overexpression (The Meiotic Recombination 11, DNA double-strand break end resection), Δrad9 (Radiation-Sensitive 9, DNA damage checkpoint regulation), Δmph1 (Mutator Phenotype Helicase 1, improvement of homologous recombinant strand extension), and PapRecT-PaSSB co-expression (stabilization of recombination intermediates) were integrated to generate the highly recombinogenic strain Y09. For vector engineering, cenARS was replaced by panARS and the endogenous promoter PGAP was employed to drive the double hammerhead ribozyme-single guide RNA-hepatitis delta virus ribozyme (double HH-sgRNA-HDV: dHgH)-mediated sgRNA expression, yielding the optimized vector Nov_pGAP_panARS_pLAT1_Cas9. These two features on K. phaffii together enhanced the EMGES 3.0 to a higher standard of transformation rate and editing efficiency. According to our results, EMGES 3.0 achieved dual-functional gene knockout efficiencies between 76.6% and 100%. For insertion of medium-long fragments (>4.5 kb), the efficiency achieved 93.3%. In addition, the one-step integration of ultra-long fragments (>16 kb) achieved 14.8%, which was reported for the first time. Furthermore, the efficiency of simultaneous long-fragment integration at three neutral loci reached 38.4% (>15 kb). We applied the system for one-step production of free fatty acids (FFAs, yield: 5.82 ∼ 7.30 mg/L/OD600) and resveratrol (yield: 1.14 ∼ 1.28 mg/L) using methanol as the sole carbon source. EMGES 3.0 provides a robust technical foundation for complex compounds biosynthesis and high-yield industrial strains, while also advancing K. phaffii as an industrial synthetic biology chassis for efficient C1 utilization.

CRISPR-Cas Systems↗

PubMed2026

THE1 repeats: Ancient endogenous retroviruses rampaging behind sarcoid myopathy.

Transposon-like human element (THE1) repeats are classified as part of the ERVL (endogenous retrovirus-like) -MaLR family and were previously thought to be remnants derived from retroviral infections in the germ cells of ancient primates. At present, at least 60,561 THE1 repeats have been mapped in the human genome, and they are divided into four subfamilies: THE1A, THE1B, THE1C, and THE1D. Although these repeats are largely silenced in somatic cells through multiple mechanisms, some of them are known to influence gene functions and genomic structures, particularly under disease conditions. Recently, we happened to identify the CIR1 driven from THE1B in the individuals with sarcoid myopathy (SM), which is characterized by the non-caseating granulomatous inflammation recognized in skeletal muscles biopsy. Retrotrans-genomics based on the COFFEE method showed that various THE1B fusion transcripts (THE1B/FTs) and retrotransposon-associated genomic dynamics in muscle of SM, including (1) splicing of THE1B to neighborhood exons, (2) expression of THE1B associated long noncoding RNAs in the intergenic regions, (3) ectopic expression of tissue specific transcripts driven by THE1B in SM, and (4) identification of a read through transcript, SIRPB1-SIRPD. Seven of eight THE1B/FTs downregulated by tofacitinib in cutaneous sarcoidosis were confined to the granuloma associated macrophage in sarcoidosis. The remaining transcript, SMFT4 exhibited stronger expression in mesenchymal stem cells. In addition, ectopic co-expression of canonical and LINE1 and THE1A-exonizing IL23R transcripts were identified in SM. This directory of THE1B/FTs in SM may serve as a fundamental molecular basis for understanding the pathogenesis of sarcoidosis.

Journal Article↗

PubMed2026

Integrated bioinformatics analysis reveals cross-talking hub genes and therapeutic agents between sepsis and acute myocardial infarction.

BACKGROUND: Sepsis and acute myocardial infarction (AMI) are two significant diseases that may share overlapping etiological mechanisms. This study aims to systematically identify core genes common to both conditions and to explore their potential as therapeutic targets and drug candidates through an integrative analysis of clinical data and bioinformatics. METHODS: The AMI dataset was obtained from the GEO database, and RNA sequencing data were collected from blood samples of patients with sepsis at our hospital. Common genes were identified using differential expression gene analysis (DEG) and weighted gene co-expression network analysis (WGCNA). Functional enrichment analyses, including Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis, were performed. A protein-protein interaction (PPI) network was constructed, and hub genes were identified using the MCC/Degree algorithm. Diagnostic value was assessed via receiver operating characteristic curve analysis. Immune infiltration patterns, single-cell sequencing data, and molecular docking simulations were employed to evaluate immune relevance and identify potential therapeutic compounds. RESULTS: A total of 417 genes were identified between sepsis and AMI, with enrichment analysis revealing significant involvement in inflammatory responses. Three hub genes-JAK2, MYD88, and TIMP1-were selected for further investigation. ROC curves confirmed their strong diagnostic performance for both diseases. Immune infiltration analysis showed that these core genes were significantly correlated with the infiltration levels of various immune cell types. Molecular docking indicated that quercetin exhibited stable binding affinity with the proteins encoded by these genes. qPCR validation further confirmed the upregulation of these three genes, supporting the anti-inflammatory effects of quercetin as a potential targeted therapy. CONCLUSION: JAK2, MYD88, and TIMP1 were identified as shared core genes in sepsis and AMI. These genes not only serve as potential diagnostic biomarkers but also offer novel targets for developing common therapeutic strategies for both conditions. Furthermore, quercetin emerges as a promising candidate for targeted treatment.

Humans↗

PubMed2026

Integrated Multi-Omics Analysis Reveals the Genetic Basis of Phenotypic Variation in Tibetan Sheep.

Body size is a key economic trait influencing the profitability of farmed animals. This study used genome-wide association studies (GWAS) to identify five single nucleotide polymorphisms (SNPs) significantly associated with body size in the Tibetan sheep population, advancing molecular breeding and providing a basis for genomic selection. These SNPs are located within five candidate genes. SNaPshot validated GWAS results, demonstrating significant correlations between candidate SNPs and body size traits in Tibetan sheep. Concurrently, hematoxylin and eosin staining, alongside muscle fiber analysis, confirmed pronounced morphological differences in muscle tissue between sheep of varying conformation. Therefore, transcriptome and proteomics were performed on the longest dorsi muscle from large and small Tibetan sheep of both sexes. The transcriptome, together with weighted gene co-expression network analysis (WGCNA), identified VEPH1 and PRKG1 as core genes regulating body characteristics in Tibetan sheep through their involvement in the PI3K-Akt signaling pathway and pathways related to fat deposition. The integrative analyses demonstrated significantly different expression of CARNS1 and CRYAB at both transcriptional and protein levels between the muscles of large- and small-sized Tibetan sheep of both sexes, suggesting their importance in body size traits by influencing muscle morphology. This study provides valuable genomic resources that advance sheep genetics research.

GWAS↗

PubMed2026

Microglial GRB2 is essential for brain ventriculogenesis and CSF homeostasis.

Microglia play essential yet poorly understood roles in brain development, including axon guidance, regulation of neurogenesis, and pruning of neuronal projections. Congenital hydrocephalus (CH), characterized by enlarged cerebrospinal fluid (CSF)-filled ventricles, is a leading cause of pediatric brain surgery, but its molecular mechanisms remain unclear. We have identified what we believe to be novel, recurrent, damaging missense variants in the SH3-binding domain of the adaptor protein Growth Factor Receptor-Bound Protein 2 (GRB2) in unrelated patients with CH. GRB2 is significantly co-expressed with one of its known upstream receptor tyrosine kinase partners, CSF1R, in the developing human brain, particularly in a microglial subtype associated with regulation of neural stem cells. Immunoprecipitation validated GRB2-CSF1R binding in mouse microglial cells and human monocyte cell line. Cx3cr1-Grb2fl/fl mice engineered with conditional deletion of Grb2 in microglia exhibit congenital absence of microglia and early postnatal severe communicating (non-obstructive) hydrocephalus, mimicking GRB2-mutant patients. The severe ventriculomegaly of Cx3cr1-Grb2fl/fl mice is associated with both depletion of cerebral cortical neurons and impairment of glia-lymphatic-mediated CSF flow. Together, these findings implicate a role of GRB2 in microglia that could be essential for brain development and CSF homeostasis.

Genetics↗

PubMed2026

Hepatic metabolic adaptation to endurance exercise: temporal and sex differences by multiomics integration and validation.

BACKGROUND: Although endurance exercise benefits liver health, sex-specific adaptive trajectories remain unclear. This study mapped dynamic liver adaptation in males and females during prolonged training and identified underlying molecular programs. METHODS: Using publicly available time-resolved liver multi-omics data generated by the Molecular Transducers of Physical Activity Consortium (MoTrPAC), we established a computational pipeline for differential analysis of transcriptomic, proteomic, phosphoproteomic, and metabolomic data with FDR correction, followed by FGSEA pathway enrichment. Kinase activities were inferred through ortholog mapping and PhosphoSitePlus. Cross-omics co-expression networks were constructed using WGCNA and topological overlap to link omics features with physiological phenotypes. For experimental validation, liver tissues were collected from endurance-trained Sprague-Dawley rats, and key nodes were confirmed by Western blotting, qRT-PCR, and immunofluorescence/immunohistochemical staining. Public scRNA-seq data were further integrated to map multi-omics signals to single-cell resolution and assess functional changes in specific cell types. RESULTS: The hepatic response to exercise stress was stage-specific, shifting from early transcriptional activation to later proteomic and metabolic remodeling. Multi-omics integration revealed distinct sex-associated adaptive trajectories: males were more strongly associated with energy metabolism, redox-related programs, and amino acid/organic acid catabolism, whereas females showed prominent membrane lipid remodeling, proteostasis -related programs, and mitochondrial/ribosomal translational features. Single-cell analysis showed that tissue remodeling occurred without major lineage turnover, instead involving altered communication among pre-existing cell communities. Validation of PPP1R3G identified a protein-dominant exercise-responsive marker, supporting the contribution of post-transcriptional or protein-level regulation. CONCLUSIONS: Hepatic adaptation to endurance stress follows a cross-omics evolutionary pattern with sex-specific reprogramming of energy supply and homeostatic maintenance. This time-resolved framework clarifies how exercise improves liver function and supports sex-oriented metabolic interventions and therapeutic target discovery.

Animals↗

PubMed2026

New insights into diagnostic values and mechanisms of ferroptosis associated with immune infiltration in diabetic kidney disease.

The pathogenesis of diabetic kidney disease (DKD) is complex and closely related to ferroptosis and immune dysregulation, but the relevance is unclear. The present study investigates the potential mechanisms of ferroptosis-related genes (FRGs) in DKD and their relationship with the immune-inflammatory response. It searches for new diagnostic biomarkers to help diagnose and treat DKD. Four Gene Expression Omnibus (GEO) datasets, GSE30528, GSE30529 and GSE30122 as the test set, and GSE96804 for validation, were analyzed. FRGs were obtained from GeneCards, and 47 ferroptosis-related differentially expressed genes (FRDEGs) were identified by intersecting with DKD-related differentially expressed genes. Functional enrichment analyses, including Gene Ontology, Kyoto Encyclopedia of Genes and Genomes, Gene Set Enrichment Analysis and Gene Set Variation Analysis, revealed that these FRDEGs are primarily associated with ferroptosis, hypoxia response and immune inflammation. Subsequently, the weighted gene co-expression network analysis (WGCNA) was employed to expand the ferroptosis-related gene network, and intersection of the 47 FRDEGs with key WGCNA module genes yielded 10 key genes. Based on the 10 key genes, the least absolute shrinkage and selection operator and support vector machine algorithms identified three hub genes [chemokine ligand 5 (CCL5), forkhead box C1 (FOXC1) and lactotransferrin (LTF)] for DKD diagnosis. Receiver operating characteristic curves confirmed their diagnostic value, with FOXC1 and LTF validated in the independent dataset. Immune infiltration analysis via CIBERSORT revealed eight immune cell types with significantly different infiltration levels between the DKD and control group in the integrated GEO datasets. Notably, both LTF and CCL5 showed a significant positive correlation with gamma delta T cells (γδT). Quantitative PCR results confirmed differential expression of the three hub genes in the DKD group, with elevated expression observed in DKD mice following intervention with rosiglitazone and hyperoside.

bioinformatics analysis↗

PubMed2026

Gastric amphicrine carcinoma in the stomach: an unexpected presentation of MUTYH-associated polyposis.

Amphicrine carcinomas of the stomach, defined by dual exocrine and neuroendocrine differentiation within the same neoplastic cell, are exceedingly rare. MUTYH-associated polyposis (MAP) is an autosomal recessive polyposis syndrome characterized by multiple colorectal adenomas and variable upper gastrointestinal involvement; however, amphicrine carcinomas have not been previously documented in this setting. We report a gastric amphicrine carcinoma arising in the background of extensive fundic gland polyposis in a patient with MAP. Endoscopy revealed a 3.5-cm flat elevated lesion in the gastric fundus amid extensive fundic gland polyposis. Histologically, the tumor consisted of a single population of cells exhibiting combined glandular and neuroendocrine differentiation without zonal or biphasic architecture, and many of these cells demonstrated true amphicrine morphology. Immunohistochemistry confirmed co-expression of cytokeratin and the neuroendocrine markers chromogranin A and synaptophysin in the same cell population. Germline targeted next-generation sequencing identified biallelic MUTYH variants in trans (c.733C>T, p.Arg245Cys [likely pathogenic]; c.842C>T, p.Ala281Val [variant of uncertain significance]), supporting a diagnosis of MAP. To our knowledge, this is the first reported case of a gastric amphicrine carcinoma in a MAP patient, expanding the spectrum of MAP-associated upper gastrointestinal neoplasia and underscoring the importance of vigilant endoscopic surveillance in hereditary polyposis syndromes.

Carcinoma↗