SEARCH · PubMed Health

Results for “comparative transcriptomics”

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 217 records · Page 12Linked to original sources

PubMed2026

A High-Resolution Stereo-Seq Spatial Transcriptomic Resource for Adult Holstein Cattle Liver.

The bovine liver is a highly compartmentalized organ that plays essential roles in continuous gluconeogenesis and nitrogen recycling; however, its spatial molecular architecture has remained largely uncharacterized due to the limitations of traditional bulk and single-cell approaches. To address this gap, Spatial Enhanced Resolution Omics-sequencing (Stereo-seq) was utilized to generate a subcellular-resolution (500 nm) transcriptomic map of an adult Holstein cattle liver, and a refined reference-guided workflow was implemented to overcome standard annotation limitations in livestock. Raw sequencing data were processed using the Stereo-seq Analysis Workflow and analyzed with Stereopy, Seurat, SingleR, and reference-guided workflows. Spatial aggregation was evaluated at Bin20, Bin50, Bin100, Bin150, and Bin200. Increasing bin size increased molecular identifier counts and detected-gene complexity while progressively reducing spatial granularity. Bin50, corresponding to 50 × 50 DNA nanoballs and an approximate nominal footprint of 25 × 25 µm, was therefore selected as a practical intermediate aggregation level for the primary analyses. Quality-control assessment, Leiden clustering, UMAP visualization, reference-based cell-type annotation, cluster-marker analysis, and spatial mapping of canonical hepatic genes demonstrated preservation of biologically interpretable liver transcriptional organization. Raw sequencing data processed spatial matrices, annotated objects, and analysis code are publicly available to support reanalysis and computational benchmarking. In summary, we present a Stereo-seq spatial transcriptomic resource generated from liver tissue of an adult Holstein cow. This initial resource provides a valuable foundation for future studies of bovine liver biology, comparative genomics, and the spatial basis of livestock health and production traits.

Animals↗

PubMed2025

Integrated genomic, transcriptomic, and metabolomic analyses of Chrysanthemum aromaticum provide insights into the volatile terpene biosynthesis.

Chrysanthemum aromaticum is renowned for its uniformly emitted strong and attractive scent, primarily attributed to volatile terpenes. Despite its commercial and horticultural significance, the molecular mechanisms underlying volatile terpene production in C. aromaticum remain largely unexplored. Here, we present the haplotype-resolved genome assembly of C. aromaticum, with a total size of 3.10 Gb, comprising nine anchored chromosomes with a contig N50 of 30.66 Mb and a scaffold N50 of 350.58 Mb. Phylogenetic analyses revealed a distant relationship between C. aromaticum and C. indicum, suggesting that C. aromaticum likely represents a distinct species rather than a variety of C. indicum. Through integrated genomic, transcriptomic, metabolomic, and biochemical analyses, we identified seven TPS involved in monoterpene biosynthesis and six TPS for sesquiterpene biosynthesis. Notably, comparative genomic analysis revealed a gene cluster for α-bisabolol biosynthesis in C. aromaticum, which has specifically expanded in Chrysanthemum species through tandem gene duplications, contributing to the elevated accumulation of α-bisabolol in the leaves of C. aromaticum. Our study provides important insights into the biosynthesis of volatile terpenes, highlighting the genetic basis for C. aromaticum's unique aromatic profile.

Chrysanthemum↗

PubMed2026

Transcriptomic characterization of the intestine in Stichopus monotuberculatus under gradient temperature stress and HSP gene family-mediated molecular adaptation.

The increasing frequency of extreme temperature events under climate change poses a growing threat to the stability of tropical sea cucumber aquaculture. To characterize the molecular responses of the tropical sea cucumber Stichopus monotuberculatus to acute temperature stress, juveniles were exposed for 96 h to 15 °C, 20 °C, 25 °C, 30 °C, and 35 °C, followed by transcriptomic profiling of the intestine. By transcriptomic analysis, 2258, 634, 1618, and 2980 differentially expressed genes (DEGs) were identified at 15, 20, 30, and 35 °C compared to control, respectively. More DEGs were generally detected at temperatures further from 25 °C, with the 35 °C group showing the largest transcriptional response. Although cold and heat stress both affected metabolism and protein homeostasis, their enrichment profiles differed. At 15 °C, DEGs were mainly enriched in the spliceosome and p53 signaling pathways, highlighting RNA processing and p53 signaling as prominent features of the cold-stress response. At 35 °C, DEGs were mainly enriched in the PI3K-Akt signaling pathway, ubiquitin-mediated proteolysis, and mitophagy, indicating enhanced regulation of cell survival, protein turnover, and mitochondrial quality control. HSP genes also responded differently to cold and heat stress. Most HSP70 and HSP90 family members were downregulated at low temperatures, whereas HSP70 genes and small heat shock proteins were markedly upregulated at high temperatures. Overall, the intestinal transcriptome showed distinct responses to cold and heat stress. These results identify pathways and HSP genes potentially involved in the temperature response of S. monotuberculatus and provide useful information for evaluating temperature tolerance and defining suitable temperatures for its aquaculture.

Heat shock protein↗

PubMed2026

Genome-wide characterization of NOD-like receptor genes links NLR repertoire evolution to spleen immune responses after Aeromonas hydrophila challenge in the Chinese spiny frog (Quasipaa spinosa).

NOD-like receptors (NLRs) are cytosolic pattern-recognition receptors that detect pathogen-associated and damage-associated molecular patterns and mediate innate immune signaling in vertebrates. However, the genomic repertoire, evolutionary diversification, and infection-associated expression of NLR genes remain poorly defined in non-model amphibians. In this study, 66 NLR genes were identified from the Chinese spiny frog (Quasipaa spinosa) genome and designated as QsNLR1-QsNLR66. These genes were unevenly distributed across chromosomes and were classified into three phylogenetic groups, with most members exhibiting conserved motif architectures. Gene duplication analysis indicated that dispersed duplication was the main contributor to QsNLR expansion. Synteny analysis detected five conserved orthologous gene pairs between Q. spinosa and Pelophylax nigromaculatus, suggesting partial conservation of NLR genomic organization between the two amphibians. Ka/Ks analysis showed that several duplicated gene pairs, including NLRC3-like/QsNLR36 and NLRC3-like/QsNLR50, exhibited Ka/Ks ratios greater than one, suggesting potential sequence divergence after duplication. Spleen RNA sequencing (RNA-seq) after Aeromonas hydrophila challenge revealed enrichment of immune-related Gene Ontology (GO) terms and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways. Weighted gene co-expression network analysis linked several QsNLRs to infection-associated modules, among which QsNLR57 was co-expressed with CYBB, ADAM17, SPI1, and HK2. RT-qPCR using time-matched phosphate-buffered saline (PBS) controls showed distinct temporal patterns, with stronger induction of QsNLR29, QsNLR57, and QsNLR66 and weaker or delayed responses of QsNLR50 and QsNLR56. These results characterize the NLR repertoire of Q. spinosa and identify infection-associated QsNLR candidates for future studies of antibacterial immunity in amphibians.

Animals↗

PubMed2026

CD36 Influences Leukemia Progression in MLL-AF9-Driven AML by Modulating the Leukemia Immune Microenvironment.

CD36, a fatty-acid translocase, is increasingly implicated in acute myeloid leukemia biology and treatment resistance, yet its contribution to leukemogenesis is still unclear. Using the MLL-AF9 model, we transduced hematopoietic stem/progenitor cells (HSPCs) from Cd36-knockout (KO) or wild-type (WT) mice and assessed leukemic potential with in vitro assays, transplants, and transcriptomic, metabolomic, and immune profiling. Both Cd36KO- and Cd36WT-HSPCs underwent efficient MA9-driven transformation, with comparable colony formation and Hox/Meis1 pathway activation, indicating Cd36 is dispensable for leukemic initiation. However, Cd36 deletion markedly attenuated disease progression, reducing leukemic burden and extending survival in irradiated mice (median 22 vs. 15 days, P = 0.001). Effects were strikingly amplified in immunocompetent, non-irradiated recipients (median 63 vs. 22 days, P = 0.002), revealing immune-dependent suppression. Immune profiling showed enhanced CD4⁺ and CD8⁺ T cell infiltration, reduced CD4⁺CD25⁺ regulatory-like cells, and lower Tim-3 expression in Cd36KO-MA9 spleens, consistent with a less exhausted, more effective anti-leukemic T cell response. Despite enhanced T cell infiltration, TCR repertoires remained conserved, indicating functional reprogramming rather than clonal selection. Consistent with a suppressive leukemia immune microenvironment, RNA-seq gene set enrichment analysis identified upregulation of inflammatory (TNFα/NF-κB) and hypoxic pathways in Cd36WT-MA9 cells. Untargeted metabolomics revealed metabolic shifts in Cd36KO cells, involving a reduction in three key metabolites, UDP-GlcNAc, UDP-Galactose/UDP-Glucose, and O-Phospho-L-Serine, that likely support an immune evasion mechanism. These findings demonstrate that while Cd36 is not essential for MLL-AF9-mediated transformation, its cell-intrinsic expression in leukemic cells suppresses anti-leukemic immunity and accelerates progression. This positions CD36 as a promising target to enhance immune surveillance and limit AML aggressiveness.

Acute Myeloid Leukemia (AML)↗

PubMed2026

Losartan shows limited benefit in preclinical models of Geleophysic dysplasia.

Geleophysic dysplasia (GD) is a rare genetic disorder characterized by short stature, joint contractures, and cardiopulmonary complications, with early mortality, and linked to mutations in ADAMTSL2 (GD1), FBN1 (GD2), or LTBP3 (GD3) genes. These mutations are hypothesized to disrupt extracellular matrix (ECM) organization and enhance transforming growth factor beta (TGF-β) signaling. Losartan, an angiotensin II receptor blocker, has been proposed to mitigate TGF-β-mediated pathologies. In this study we tested the efficacy of losartan as a therapeutic drug for GD. We evaluated losartan's therapeutic potential using Adamtsl2 p.A165T mutant mice and patient-derived fibroblasts. Survival, growth, TGF-β signaling, and ECM protein expression were assessed. Losartan did not improve survival or growth in our mutant mice. Compared with control fibroblasts, patient-derived fibroblasts showed reduced basal TGF-β1 secretion. Consistent with this finding, transcriptomic analyses did not reveal activation of the TGF-β signaling pathway, and no differences in SMAD phosphorylation were observed between patient and control cells. Losartan treatment failed to modulate TGF-β signaling or ECM protein incorporation. These results suggest limited benefits of losartan in GD and challenge the notion of TGF-β dysregulation in GD pathogenesis, indicating a need for alternative targeted therapies.

Losartan↗

PubMed2026

Perineuronal net degradation in aggressive glioblastomas with KANK1::NTRK2 fusions.

BACKGROUND: Approximately 10% of glioblastomas harbor targetable genomic fusions. NTRK2 participates in a variety of fusion events that drive tumorigenesis. Two previous reports have described KANK1::NTRK2 fusions in adult glioblastoma patients with poor survival. METHODS: We performed a retrospective analysis of glioblastoma patients treated at Dartmouth-Hitchcock Medical Center (DHMC) from 2020 to 2025 to identify cases harboring KANK1::NTRK2 fusions. Clinical presentation, treatment, histopathologic features, and outcomes were reviewed. In addition, we conducted GeoMx whole-transcriptome and high-plex proteomic digital spatial profiling of a KANK1::NTRK2-positive glioblastoma and a comparator tumor from a long-term survivor. Candidate biomarkers were orthogonally validated using immunohistochemistry and/or immunofluorescence. RESULTS: Two patients with KANK1::NTRK2 fusion glioblastoma were identified, both demonstrating rapid progression, therapeutic resistance, and survival of less than 7 months. Proteomic profiling showed increased expression and activation of canonical NTRK2 downstream signaling pathways, particularly MEK1/2 and ERK1/2. This was accompanied by upregulation of extracellular matrix remodeling enzymes, including MMP3, MMP14, and ADAM15, along with reduced expression of extracellular matrix-associated transcripts and perineuronal net components in particular compared to a non-fusion glioblastoma. CONCLUSIONS: These limited, hypothesis-generating findings suggest constitutive NTRK2 signaling may promote coordinated extracellular matrix degradation and remodeling, potentially facilitating rapid and aggressive tumor growth and invasion in a subset of glioblastomas.

NTRK gene fusion↗

PubMed2025

Ancestry and somatic profile predict acral melanoma origin and prognosis.

Acral melanoma, which is not ultraviolet (UV)-associated, is the most common type of melanoma in several low- and middle-income countries including Mexico. Latin American samples are significantly underrepresented in global cancer genomics studies, which directly affects patients in these regions as it is known that cancer risk and incidence may be influenced by ancestry and environmental exposures. To address this, we characterise the genome and transcriptome of 123 acral melanoma tumours from 92 Mexican patients, a population notable because of its genetic admixture. Compared with other studies of melanoma, we found fewer frequent mutations in classical driver genes such as BRAF, NRAS or NF1. While most patients had predominantly Amerindian genetic ancestry, those with higher European ancestry had increased frequency of BRAF mutations and a lower median number of structural variants. The tumours with activating BRAF mutations have a transcriptional profile more similar to cutaneous non-volar melanocytes, suggesting that acral melanomas in these patients may arise from a distinct cell of origin compared to other tumours arising in these locations. KIT mutations were found in a subset of these tumours, and quadruple wild-type samples (non BRAF/NRAS/NF1/KIT) differed from mutated samples in their structural genomic profile and overall and recurrence-free survival patterns. Transcriptional profiling defined three expression clusters; these characteristics were associated with recurrence-free and overall survival. We highlight potential novel low-frequency drivers, such as PTPRJ, NF2 and RDH5. Our study enhances knowledge of this understudied disease and underscores the importance of including samples from diverse ancestries in cancer genomics studies.

Journal Article↗

PubMed2026

A chromosome-scale genome of Capsicum pubescens provides insights into candidate terpene-associated gene clusters and pan variation of terpene synthases.

A chromosome-scale genome of Capsicum pubescens and comparative pan-TPS analysis support structural characterization and gene-level prioritization of a chromosome-9 terpene-associated candidate locus in this accession. Capsicum pubescens is one of the five domesticated Capsicum species, mainly cultivated in mid- to high-elevation regions of the Americas. Despite its distinctive morphology and fruit traits, genomic resources for C. pubescens remain less developed than those for the widely cultivated C. annuum. Here, we assembled a chromosome-scale reference genome for accession HNUCP0001, spanning 3.70 Gb with a scaffold N50 of 278.01 Mb. Comparative genomics revealed 679 significantly expanded gene families enriched in sesquiterpenoid and triterpenoid biosynthesis. Genome-wide biosynthetic gene-cluster mining identified multiple terpene-associated candidate loci, which were subsequently prioritized using genome-derived structural criteria and Capsicum pubescens-specific expression evidence. Subsequently, we curated the terpene synthase (TPS) repertoire and, across 16 Capsicum genomes, resolved 36 TPS orthogroups with pronounced presence/absence variation, highlighting dynamic lineage-specific diversification. Together, these analyses establish HNUCP0001 as an accession-specific genomic resource and provide a comparative framework for prioritizing terpene-associated TPS genes and candidate BGCs in Capsicum. These candidate loci, together with accession-level transcriptomic and metabolomic evidence, offer testable hypotheses for future functional studies of specialized terpenoid metabolism in C. pubescens.

Alkyl and Aryl Transferases↗

PubMed2026

Multi-omics integrative analysis provides insight into potential molecular responses to sustained high water flow in common carp (Cyprinus carpio) cultured in recirculating aquaculture.

To investigate the potential molecular responses by which water flow intensity affects the growth of common carp (Cyprinus carpio) in a recirculating aquaculture system (RAS), a control group (CG, actual water velocity 0.3 cm/s) and three sustained flow treatment groups were established, including a low-flow group (LF, 1 body length per second, bl/s), a medium-flow group (MF, 2 bl/s), and a high-flow group (HF, 3 bl/s). After 12 weeks of culture in the RAS, growth performance was compared among groups under different flow intensities. The best-performing group and the control group were then selected for the determination of intestinal digestive enzyme activities, as well as transcriptomic and whole-genome bisulfite sequencing analyses of muscle tissue. The results showed that the specific growth rate and feed intake of the HF group were significantly higher than those of the other groups (P < 0.05), whereas no significant difference in feed conversion ratio was observed among groups. Compared with the CG group, lipase activity was significantly higher in the HF group (P < 0.05), while α-amylase and trypsin activities showed increasing trends without significant differences. RNA-seq identified a total of 273 differentially expressed genes, including 72 upregulated genes and 201 downregulated genes in the HF group relative to the CG group. These genes were mainly enriched in glycolysis, pyruvate metabolism, ATP metabolism, the pentose phosphate pathway, the insulin signaling pathway, the PPAR signaling pathway, and the adipocytokine signaling pathway, indicating that sustained high water flow induced a muscle transcriptional response characterized by remodeling of energy metabolism and substrate utilization. Whole-genome bisulfite sequencing analysis showed that DNA methylation in common carp muscle occurred predominantly in the CpG context. Differentially methylated regions between the HF and CG groups were mainly distributed in transcription-related regulatory regions, including promoters, CpG islands, and CpG island shores. In promoter regions, the number of hypermethylated regions in the HF group relative to the CG group was markedly higher than that of hypomethylated regions. Integrated analysis further identified two candidate genes showing both promoter differential methylation and differential expression, namely LOC109094644 and bcorl1, suggesting that adaptation to high water flow may involve IGF-related growth regulation and remodeling of upstream transcriptional programs. The qPCR results were consistent with the transcriptomic data. Taken together, within the tested range, a sustained water flow of 3 bl/s was more conducive to the growth of common carp in the RAS, which may be associated with enhanced lipid digestion and utilization, remodeling of the muscle energy metabolic network, changes in promoter methylation, and the coordinated regulation of key candidate genes. This study provides a theoretical basis for clarifying the exercise adaptation mechanism of common carp in recirculating aquaculture and for optimizing flow velocity parameters.

Animals↗

PubMed2025

Genomic and Immune Landscape of Pancreatic Ductal Adenocarcinoma Associated with Germline Pathogenic Variants in ATM.

PURPOSE: Germline pathogenic variants (PV) in ATM increase the risk of pancreatic ductal adenocarcinoma (PDAC), but the underlying tumor biology of PDAC associated with germline PV in ATM has not been adequately explored. EXPERIMENTAL DESIGN: Whole-genome, whole-exome, and RNA sequencing were performed on PDAC tumors from 25 germline ATM PV carriers diagnosed at Mayo Clinic between 2007 and 2017. Somatic and copy-number alterations, mutational signatures, transcriptomic subtypes, and the immune landscape were evaluated. RESULTS: High-quality whole-exome and whole-genome sequencing were obtained from 21 and 15 tumors, respectively. Biallelic inactivation of ATM was observed in 87%, KRAS PV in 90%, CDKN2A homozygous loss in 60%, and TP53 alterations in <10% of these tumors. A predominant clock-like mutational signature was present in all samples. Whole-transcriptome analysis identified that the aberrantly differentiated endocrine exocrine subtype accounted for 18% of PDAC and was consistently associated with >5-year overall survival. In addition, a 28-gene expression-based signature associated with overall survival was identified and further validated in The Cancer Genome Atlas cohort. Immune landscape analysis through CODEX identified enriched CD4 T-helper cell/tumor interactions and reduced B7H3-high cell/tumor interactions in ATM PV carriers compared with noncarriers. CONCLUSIONS: The observed absence of TP53 PV and enrichment for CDKN2A alterations in ATM tumors, along with differences in the mutational signatures, transcriptomic subtypes and immune landscape, improve our understanding of the mechanistic pathways involved in PDAC development in germline ATM PV carriers and help identify potential targeted therapeutic strategies.

Humans↗

PubMed2025

HighDimMixedModels.jl: Robust high-dimensional mixed-effects models across omics data.

High-dimensional mixed-effects models are an increasingly important form of regression in which the number of covariates rivals or exceeds the number of samples, which are collected in groups or clusters. The penalized likelihood approach to fitting these models relies on a coordinate descent algorithm that lacks guarantees of convergence to a global optimum. Here, we empirically study the behavior of this algorithm on simulated and real examples of three types of data that are common in modern biology: transcriptome, genome-wide association, and microbiome data. Our simulations provide new insights into the algorithm's behavior in these settings, and, comparing the performance of two popular penalties, we demonstrate that the smoothly clipped absolute deviation (SCAD) penalty consistently outperforms the least absolute shrinkage and selection operator (LASSO) penalty in terms of both variable selection and estimation accuracy across omics data. To empower researchers in biology and other fields to fit models with the SCAD penalty, we implement the algorithm in a Julia package, HighDimMixedModels.jl.

Algorithms↗

PubMed2026

Integrated single-cell and spatial transcriptomic analyses reveal malignant epithelial glycolytic heterogeneity and spatial niche remodeling during colorectal cancer progression.

Colorectal cancer (CRC) progression is shaped by metabolic reprogramming and complex interactions within the tumor microenvironment. However, the cellular heterogeneity, spatial organization, and clinical relevance of glycolytic activity in CRC remain incompletely understood. In this study, we integrated single-cell RNA sequencing, bulk transcriptomics, and spatial transcriptomics data to systematically characterize glycolytic heterogeneity in CRC. Glycolytic activity was quantified using five independent scoring methods, consistently showing that epithelial cells exhibited the highest glycolytic activity across the two single-cell cohorts. Stratification of CopyKAT-verified aneuploid malignant epithelial cells into high-glycolysis (HG) and low-glycolysis (LG) subgroups by glycolysis scores revealed that HG cells exhibited higher stemness scores and chromosomal copy number variations. Cell-cell communication analysis revealed that, compared with LG cells, HG cells exhibited increased interaction frequency and strength with immune and stromal populations, indicating enhanced malignant epithelial-microenvironment crosstalk. Spatial transcriptomics analyses further revealed that glycolytic activity varied across normal colorectal tissue, primary CRC, and colorectal liver metastases, accompanied by progressive remodeling of epithelial-associated spatial niches and MIF-mediated intercellular communication. Bulk transcriptomic analysis identified a glycolysis-related prognostic signature with robust predictive performance, which served as an independent prognostic factor for overall survival in CRC cohorts. Collectively, these findings indicate that glycolytic heterogeneity is a key feature of CRC malignant epithelial cells and is closely associated with tumor progression, microenvironmental remodeling, and clinical outcomes.

Humans↗

PubMed2026

First-in-Class Small Molecule Inhibitor of Oncogene AVIL in Glioblastoma.

Glioblastoma multiforme (GBM) is the most prevalent and aggressive malignant primary brain tumor, marked by rapid growth, extensive invasiveness, and a median survival of only ∼15 months despite current multimodal therapy. To identify new therapeutic vulnerabilities, we investigated the actin-regulatory protein AVIL, previously implicated through a MARS-AVIL gene fusion in rhabdomyosarcoma. Comprehensive genomic and transcriptomic analyses across REMBRANDT, TCGA, and CGGA datasets revealed recurrent AVIL amplification and consistently elevated AVIL expression in GBM compared with normal brain tissue. AVIL was overexpressed across all GBM molecular subtypes and glioma stem cell (GSC) states but was nearly undetectable in normal astrocytes, neural stem cells, and brain tissues. Functional studies demonstrated that AVIL is both necessary and sufficient for glioma genesis: AVIL silencing eradicated GBM cells in vitro and suppressed xenograft growth in vivo, while AVIL overexpression enhanced proliferation, migration, and transformation. Mechanistically, AVIL drives tumor progression through actin cytoskeleton remodeling and activation of the FOXM1-LIN28B oncogenic pathway. Using a small molecule microarray screen, we identified a selective AVIL-binding compound (compound A) that potently inhibited GBM cell growth with minimal toxicity to normal astrocytes. Gene expression changes induced by compound A mirrored those following AVIL knockdown, indicating on-target activity. Compound A demonstrated robust antitumor efficacy in multiple preclinical GBM models, including orthotopic xenografts, GSC-derived tumors, patient-derived xenografts, and temozolomide-resistant GBM with favorable pharmacokinetics and blood-brain barrier penetration. The minimal AVIL expression in normal tissues and lack of phenotype in AVIL-deficient mice underscore its potential as a low-toxicity therapeutic target. Together, these findings establish AVIL as a critical oncogenic driver in GBM and introduce a first-in-class AVIL inhibitor with strong translational promise for precision neuro-oncology.

Glioblastoma↗

PubMed2025

Tensor decomposition of multi-dimensional splicing events across multiple tissues to identify splicing-mediated risk genes associated with complex traits.

Identifying risk genes associated with complex traits remains challenging. Integrating gene expression data with Genome-Wide Association Study (GWAS) through Transcriptome-Wide Association Study (TWAS) methods has discovered candidate risk genes for various complex traits. Splicing, which explains a comparable heritability of complex traits as gene expression, is under-explored due to its multidimensionality. To leverage multiple splicing events in a gene and shared splicing across tissues, we develop Multi-tissue Splicing Gene (MTSG), which employs tensor decomposition and sparse Canonical Correlation Analysis (sCCA) to extract meaningful information from high-dimensional multiple splicing events across multiple tissues. We build MTSG models using GTEx data and apply them to GWAS summary statistics of Alzheimer's disease (AD) (111,326 cases and 677,663 controls) and schizophrenia (SCZ) (36,989 cases and 113,075 controls). We identify 174 and 497 significant splicing-mediated risk genes for AD and SCZ, respectively, at Bonferroni correction. For AD, our results demonstrate significant enrichment of AD related pathways and identify additional AD risk genes not detected in the single-tissue analysis, while preserving most top genes identified in the brain frontal cortex. Consistently, for SCZ, genes identified by our brain-wide MTSG model, built from a cluster of 13 brain tissues, exhibit stronger enrichment in SCZ-relevant genes and MTSG identifies unique SCZ risk genes compared to single-tissue models. These results showcase that our MTSG models capture distinctive splicing events across tissues, which might be overlooked when using single tissue alone. Our MTSG models can be applied to other complex traits to help identify splicing-mediated disease risk genes.

Humans↗

PubMed2026

Convergent IGF2 overexpression in pheochromocytoma/paraganglioma: insights from Beckwith-Wiedemann syndrome.

Beckwith-Wiedemann syndrome (BWS) is an imprinting disorder characterized by overgrowth and tumor predisposition, caused by dysregulated expression of genes on chromosome 11p15.5. An association between BWS and pheochromocytoma/paraganglioma (PPGL) has been suggested in isolated case reports over the past fifty years, but the molecular basis for this link remains unclear. We identified four patients with BWS who developed metastatic PPGL and investigated IGF2 pathway activation in these tumors and in PPGL across various genotypes. Pan-cancer transcriptomic analysis of The Cancer Genome Atlas (TCGA) demonstrated that PPGL overexpresses IGF2, with pseudohypoxic tumors exhibiting higher expression compared to other molecular clusters. Loss of heterozygosity and loss of imprinting at 11p15.5 partially explain this overexpression, with PPGL additionally demonstrating globally elevated expression of imprinted genes compared to most other tumor types, suggesting a broader relaxation of genomic imprinting. Cognate receptor profiling revealed that PPGLs are equipped to respond to IGF2 signaling, with high expression of IGF1R and insulin receptor isoform A (IR-A). Immunohistochemistry confirmed IGF2 protein overexpression in both BWS-associated and genotypically diverse sporadic PPGLs. Our results indicate that IGF2 overexpression is a convergent molecular feature of PPGL across genotypes and suggest the IGF2 pathway as a potential diagnostic and therapeutic target.

Humans↗

PubMed2026

Transcriptomic insights into thermal stress reveal physiological trade-off between thermal stress adaptation and reproductive investment in Spodoptera litura.

Spodoptera litura, a highly polyphagous lepidopteran pest, poses a major threat to agricultural productivity due to its remarkable adaptability to diverse environmental conditions. Although heat stress is known to trigger transcriptional reprogramming in insects, the molecular mechanisms underlying thermal stress responses in S. litura remain poorly understood. In the present study, fourth-instar larvae were exposed to acute heat stress (44 °C) and compared with control conditions (27 ± 1 °C) to investigate heat-induced transcriptional alterations affecting physiology and reproduction. High-quality RNA-Seq data achieved more than 80% mapping efficiency, with a total of 15,782 transcripts were identified. Transcriptome analysis of S. litura larvae showed 323 differentially expressed genes (DEGs), of which 262 genes were significantly upregulated and 61 were downregulated in heat-stressed larvae compared to the control group. The DEGs were associated with stress response, reproduction, signalling, proteostasis, detoxification, oxidative stress, metabolism, development, and chromatin regulation. Heat shock proteins genes, including HSP70, HSP90, and HSP27, together with co-chaperones such as TRET-1, STIP1, and Starvin, were strongly upregulated, indicating enhanced cellular protection against protein damage and oxidative stress under heat stress. Conversely, key reproductive and cell cycle-related genes, including BARR, CAPD2, FEO, CDK2 and MORULA, were significantly downregulated, suggesting reproductive impairment and developmental arrest. RT-qPCR validation corroborated the RNA-Seq findings, demonstrating a heat-induced physiological trade-off that prioritizes survival over reproduction. Consistent with these molecular responses, heat-stressed insects exhibited marked reproductive impairment, including significant reductions in gonadosomatic index, eupyrene sperm bundle count, mating frequency, mating success, female calling behaviour, copulation duration, fecundity, and egg fertility. Collectively, these findings provide comprehensive insights into the molecular basis of thermal adaptation in S. litura and demonstrate that acute heat stress compromises reproductive fitness while activating conserved stress-response pathways that promote short-term survival.

Animals↗

PubMed2026

Antibacterial mechanisms and pathogen-dependent protective effects of the golden pompano LEAP2-derived peptide TroLEAP2-21.

Antimicrobial peptides (AMPs) are essential components of the innate immune system, with liver-expressed antimicrobial peptide 2 (LEAP2) playing a pivotal role in fish immunity. This study investigated the antimicrobial activity and mechanisms of TroLEAP2-21, a 21-amino-acid short peptide from golden pompano (Trachinotus ovatus), against Gram-positive (Lactococcus garvieae, Staphylococcus epidermidis) and Gram-negative (Vibrio alginolyticus, Vibrio harveyi) bacteria. The predicted three-dimensional structure and helical wheel projection of TroLEAP2-21 suggested typical AMP-like physicochemical features. troleap2 expression in the liver and intestine of T. ovatus was significantly upregulated post L. garvieae or V. harveyi infection, suggesting its potential involvement in antibacterial defense. In vitro, TroLEAP2-21 exhibited antibacterial activity against the tested bacterial strains, with membrane disruption, increased membrane permeability, cytoplasmic leakage, and membrane depolarization observed after peptide treatment. Gel retardation assays further indicated species-dependent association of TroLEAP2-21 with bacterial genomic DNA. In vivo, under the tested intraperitoneal injection conditions, TroLEAP2-21 was associated with improved survival and reduced tissue damage in V. harveyi-infected T. ovatus, whereas no significant survival benefit was observed against L. garvieae. Transcriptomic analysis at 48 h post-infection showed transcriptional changes in immune-related DEGs (rsad2, mx1/mx2, il-8) and enrichment of TLR and Jak-STAT signaling pathways at the transcriptional level in peptide-treated fish. FISH showed the tissue localization of tnf-α and nf-κb transcripts and revealed treatment-associated changes in fluorescence signals, and qRT-PCR of eight immune genes supported transcriptomic results in tissues at 48 h post-infection. Collectively, these findings characterize TroLEAP2-21 as a short LEAP2-derived peptide with antibacterial and immunomodulatory activities. Its comparative advantages over other LEAP2-related peptides and its practical application potential remain to be further investigated.

AMPs↗