SEARCH · Search PubMed

Results for “transcriptomic analysis”

Search indexed PubMed citations on genomics, clinical trials, systematic reviews and public health. Explore titles, authors and supplied subject terms, then open the PubMed record.

SourceResearch connectionFrom yearThrough year

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

At least 433 records · Page 24Linked to original sources

PubMed2025

Pitavastatin, Procollagen Pathways, and Plaque Stabilization in Patients With HIV: A Secondary Analysis of the REPRIEVE Randomized Clinical Trial.

IMPORTANCE: In a mechanistic substudy of the Randomized Trial to Prevent Vascular Events in HIV (REPRIEVE) randomized clinical trial, pitavastatin reduced noncalcified plaque (NCP) volume, but specific protein and gene pathways contributing to changes in coronary plaque remain unknown. OBJECTIVE: To use targeted discovery proteomics and transcriptomics approaches to interrogate biological pathways beyond low-density lipoprotein cholesterol (LDL-C), relating statin outcomes to reduce NCP volume and promote plaque stabilization among people with HIV (PWH). DESIGN, SETTING, AND PARTICIPANTS: This was a post hoc analysis of the double-blind, placebo-controlled, REPRIEVE randomized clinical trial. Participants underwent coronary computed tomography angiography (CTA), plasma protein analysis, and transcriptomic analysis at baseline and 2-year follow-up. The trial enrolled PWH from April 2015 to February 2018 at 31 US research sites. PWH without known cardiovascular diseases taking antiretroviral therapy and with low to moderate 10-year cardiovascular risk were eligible. Data analyses were conducted from October 2023 to February 2024. INTERVENTION: Oral pitavastatin calcium, 4 mg per day. MAIN OUTCOMES AND MEASURES: Relative change in plasma proteomics, transcriptomics, and noncalcified plaque volume among those receiving treatment vs placebo. RESULTS: Among 558 individuals (mean [SD] age, 51 [6] years; 455 male [82%]) included in the proteomics assessment, 272 (48.7%) received pitavastatin and 286 (51.3%) received placebo. After adjusting for false discovery rates, pitavastatin increased abundance of procollagen C-endopeptidase enhancer 1 (PCOLCE), neuropilin 1 (NRP-1), major histocompatibility complex class I polypeptide-related sequence A (MIC-A) and B (MIC-B), and decreased abundance of tissue factor pathway inhibitor (TFPI), tumor necrosis factor ligand superfamily member 10 (TRAIL), angiopoietin-related protein 3 (ANGPTL3), and mannose-binding protein C (MBL2). Among these proteins, the association of pitavastatin with PCOLCE (a rate-limiting enzyme of collagen deposition) was greatest, with an effect size of 24.3% (95% CI, 18.0%-30.8%; P < .001). In a transcriptomic analysis, individual collagen genes and collagen gene sets showed increased expression. Among the 195 individuals with plaque at baseline (88 [45.1%] taking pitavastatin, 107 [54.9%] taking placebo), changes in NCP volume were most strongly associated with changes in PCOLCE (%change NCP volume/log2-fold change = -31.9%; 95% CI, -42.9% to -18.7%; P < .001), independent of changes in LDL-C level. Increases in PCOLCE related most strongly to change in the fibro-fatty (<130 Hounsfield units) component of NCP (%change fibro-fatty volume/log2-fold change = -38.5%; 95% CI, -58.1% to -9.7%; P = .01) with a directionally opposite, although nonsignificant, increase in calcified plaque (%change calcified volume/log2-fold change = 34.4%; 95% CI, -7.9% to 96.2%; P = .12). CONCLUSIONS AND RELEVANCE: Results of this secondary analysis of the REPRIEVE randomized clinical trial suggest that PCOLCE may be associated with the atherosclerotic plaque stabilization effects of statins by promoting collagen deposition in the extracellular matrix transforming vulnerable plaque phenotypes to more stable coronary lesions. TRIAL REGISTRATION: ClinicalTrials.gov Identifier: NCT02344290.

Humans↗

PubMed2001

Gene discovery using the serial analysis of gene expression technique: implications for cancer research.

Cancer is a genetic disease. As such, our understanding of the pathobiology of tumors derives from analyses of the genes whose mutations are responsible for those tumors. The cancer phenotype, however, likely reflects the changes in the expression patterns of hundreds or even thousands of genes that occur as a consequence of the primary mutation of an oncogene or a tumor suppressor gene. Recently developed functional genomic approaches, such as DNA microarrays and serial analysis of gene expression (SAGE), have enabled researchers to determine the expression level of every gene in a given cell population, which represents that cell population's entire transcriptome. The most attractive feature of SAGE is its ability to evaluate the expression pattern of thousands of genes in a quantitative manner without prior sequence information. This feature has been exploited in three extremely powerful applications of the technology: the definition of transcriptomes, the analysis of differences between the gene expression patterns of cancer cells and their normal counterparts, and the identification of downstream targets of oncogenes and tumor suppressor genes. Comprehensive analyses of gene expression not only will further understanding of growth regulatory pathways and the processes of tumorigenesis but also may identify new diagnostic and prognostic markers as well as potential targets for therapeutic intervention.

Gene Expression Profiling↗

PubMed2024

An isoform-resolution transcriptomic atlas of colorectal cancer from long-read single-cell sequencing.

Colorectal cancer (CRC) ranks as the second leading cause of cancer deaths globally. In recent years, short-read single-cell RNA sequencing (scRNA-seq) has been instrumental in deciphering tumor heterogeneities. However, these studies only enable gene-level quantification but neglect alterations in transcript structures arising from alternative end processing or splicing. In this study, we integrated short- and long-read scRNA-seq of CRC samples to build an isoform-resolution CRC transcriptomic atlas. We identified 394 dysregulated transcript structures in tumor epithelial cells, including 299 resulting from various combinations of splicing events. Second, we characterized genes and isoforms associated with epithelial lineages and subpopulations exhibiting distinct prognoses. Among 31,935 isoforms with novel junctions, 330 were supported by The Cancer Genome Atlas RNA-seq and mass spectrometry data. Finally, we built an algorithm that integrated novel peptides derived from open reading frames of recurrent tumor-specific transcripts with mass spectrometry data and identified recurring neoepitopes that may aid the development of cancer vaccines.

Colorectal Neoplasms↗

PubMed2025

Comprehensive Transcriptome Annotation of Thousands of HIV-1 Genomes.

Alternative splicing in HIV-1 has been a central focus of decades of research, uncovering key mechanisms of viral gene regulation, immune evasion, and therapeutic response - yet, no reference resource has existed to support transcriptome-wide analysis, limiting adoption of modern computational methods. We present HIV Atlas (https://ccb.jhu.edu/HIV_Atlas), the first reference-quality annotation of HIV-1 and SIV transcriptional diversity. We manually curated transcriptomes for HIV-1HXB2 and SIVmac239 and developed Vira, an automated annotation-transfer method specifically designed to address unique challenges of viral genome biology, to generate high-quality annotations for 2,077 complete HIV-1 genomes. Using the resources presented in our work, we evaluated conservation of splice sites, revealing near-perfect preservation of major donors and acceptors. Furthermore, using several public datasets, we demonstrate how HIV Atlas enhances methodology, improves the quality and novelty of results, and opens novel avenues for research, supporting more accurate and comprehensive analyses of bulk, single-cell, and spatial RNA-seq in HIV-1 studies.

Journal Article↗

PubMed2024

Genetic association between epilepsy and gliomas: Insights from Mendelian randomization and single-cell transcriptomic analyses.

BACKGROUND: Seizures are prevalent in glioma patients, especially in those with low-grade gliomas. The interaction between gliomas and epilepsy involves complex biological mechanisms that are not fully understood. METHODS: We collected Genome-Wide Association Study data for epilepsy and gliomas, performed differential expression analysis, and conducted Gene Ontology (GO) enrichment analysis on the identified genes. Single-cell RNA sequencing data (scRNA-seq) from GSE221534 dataset in Gene Expression Omnibus (GEO) were used to analyze cell-cell interactions within glioma samples from patients with and without epilepsy. RESULTS: Mendelian Randomization (MR) analysis revealed significant associations between genetic variants related to epilepsy and glioma risk, suggesting a potential causal relationship, especially in astrocytomas. Differential expression analysis identified epilepsy-related genes that were significantly upregulated in astrocytoma tissues compared to normal brain tissues. GO enrichment analysis indicated that these genes are involved in critical biological processes such as neurogenesis and cellular signaling. The scRNA-seq analysis showed, compared to non-epileptic samples, glioma stem cells, microglia, and NK cells are increased in the core regions of astrocytomas in epileptic patients. Additionally, intercellular communication between tumor cells and other non-tumor cells is markedly enhanced in astrocytoma samples from epileptic patients. CONCLUSION: This study provides evidence of a genetic association between epilepsy and gliomas and elucidates the biological mechanisms through which epilepsy may influence glioma progression.

Humans↗

PubMed2025

Unsupervised multiscale clustering of single-cell transcriptomes to identify hierarchical structures of cell subtypes.

BACKGROUND: Cell clustering is an essential step in uncovering cellular architectures in single-cell RNA sequencing (scRNA-seq) data. However, the existing cell clustering approaches are not well designed to dissect complex structures of cellular landscapes at a finer resolution. RESULTS: Here, we develop a multiscale clustering (MSC) approach to construct a sparse cell-cell correlation network for unsupervised identification of de novo cell types and subtypes across multiple resolutions. Based upon simulated silver- and gold-standard data as well as real scRNA-seq data in diseases, MSC demonstrates significantly improved performance compared to established benchmark methods and reveals a biologically meaningful cell hierarchy to facilitate the discovery of novel disease-associated cell subtypes and mechanisms. CONCLUSIONS: We present MSC as a new single-cell multiscale clustering framework as a powerful tool for advancing discoveries in disease-associated cell populations using single-cell sequencing data.

Single-Cell Analysis↗

PubMed2026

RNF43 Mutations Are Associated With the Classical Molecular Subtype, Vigorous Antitumor Immune Responses, and Prolonged Survival in Pancreatic Adenocarcinoma.

RNF43 mutations were correlated with microsatellite status in colorectal cancer and with fewer and later recurrences in pancreatic ductal adenocarcinoma (PDAC). Here, we undertake a detailed assessment of RNF43 mutations in PDAC. A total of 313 PDACs (308 microsatellite stable [MSS] and 5 microsatellite-instable [MSI] cases) underwent next-generation sequencing (Oncomine Tumor Mutation Load assay; Thermo Fisher). Spatial analyses (NanoString) classified PDACs according to their transcriptomic and proteomic immune signaling. Fluorescent imaging was used to define spatial compartments (tumor: pancytokeratin+/CD45- and leukocytes: pancytokeratin-/CD45+). Each of 20 PDACs with RNF43 mutations (RNF43mut) and without RNF43 mutations (RNF43wt) underwent multiplex immunofluorescence analysis to determine immune status. A total of 153 PDACs (22 RNF43mut and 131 RNF43wt cases) underwent bulk RNA sequencing to assign into molecular subtypes. Overall, 24 RNF43 mutations were identified (22 MSS PDACs and 2 MSI PDACs). The incidence of RNF43 mutations in MSS PDACs (7.1%) was consistent with The Cancer Genome Atlas (6.7%). However, RNF43 mutations were more frequent among MSI PDACs (40%). Additionally, RNF43mut had differential frequencies of other mutations (including Wnt pathway genes), higher tumor mutational burden values (5.5 mut/mb vs 1.67 mut/mb; P < .01), and significantly longer overall survival (47 vs 18 months; P < .0001) than RNF43wt. Moreover, RNF43mut exhibited significantly higher densities of CD8+ T lymphocytes, dendritic cells, and B lymphocytes (P < .001) and an upregulation of ITGAX, CD11c, CD8, and HLA-DR compared with RNF43wt. Patients with RNF43mut PDACs were more often of the classical molecular subtype (20/22, 90.9%). RNF43mut PDACs showed high tumor mutational burden values, suggesting increased neoantigen load coupled with an abundance of antigen-presenting immune cells and an upregulation of immune determinants promoting antigen presentation. All this contributes to stronger antitumor immune responses and improved clinical outcomes.

Humans↗

PubMed2026

Predicting cellular responses to perturbation across diverse contexts with State.

While machine learning models offer potential for predicting transcriptomic effects of perturbation, they currently struggle to generalize across cellular contexts. Here, we introduce State, a machine learning model that predicts perturbation effects while accounting for cellular heterogeneity within and across experiments. State is trained using single-cell gene expression data to predict perturbation effects across sets of cells. State improved discrimination of effects on large datasets by more than 30% and identified differentially expressed genes across genetic, signaling, and chemical perturbations with significantly improved accuracy compared with baselines. Its cell embeddings trained on observational data from 167 million cells enable the identification of strong perturbations in cellular contexts where no perturbations were observed during training. We further introduce Cell-Eval, a comprehensive evaluation framework that can be used to evaluate future models. Overall, the performance and flexibility of State set the stage for scaling the development of AI models of cell state.

Machine Learning↗

PubMed2026

Genome-wide identification and cold-stress-responsive expression analysis of the NOX gene family in Cucumis melo.

NADPH oxidases (NOXs) are crucial enzymes for reactive oxygen species (ROS) generation in plants and play vital roles in growth, development, and stress responses. To elucidate the sequence characteristics of the NOX gene family and its low-temperature response patterns in melon (Cucumis melo L.), this study conducted genome-wide identification and expression profiling of NOX family members using bioinformatics analysis, RNA-seq transcriptome sequencing, and real-time quantitative PCR (RT-qPCR). The results revealed that eight NOX members were identified in the melon genome, distributed across six chromosomes. All members harbored conserved domains including Ferric_reductase, FAD_binding_8, NAD_binding_6, and NADPH_Ox, and the encoded proteins were generally basic and hydrophilic. Phylogenetic analysis classified the NOX proteins into five subgroups. Synteny analysis indicated the presence of only one pair of intraspecific duplicated genes in melon, which was under purifying selection. The promoter regions contained multiple hormone- and stress-responsive cis-acting elements, with CmNOX2 and CmNOX4 harboring low-temperature responsive elements. Following treatment at 4℃ for 24 h and 48 h, leaf relative electrolyte leakage (REL) increased from 28.33% to 42.67% and 52.67%, respectively; transcriptome analysis identified 5,633 and 6,882 differentially expressed genes (DEGs), respectively. Cold-responsive genes exhibited significant differential expression, with SLAC1 and CPK19 showing sustained upregulation. RT-qPCR results demonstrated that the expression of CmNOX2, CmNOX5, CmNOX6, and CmNOX7 was significantly downregulated after low-temperature treatment, whereas CmNOX4 expression was significantly upregulated at 48 h. Integrating promoter elements and expression characteristics, CmNOX4 may represent an important candidate gene involved in melon low-temperature response. This study systematically characterized the structure, evolution, and expression patterns of the melon NOX gene family, identified candidate genes responsive to low temperature, and provides a reference for further investigation into the mechanisms underlying melon cold adaptation.

Cucumis melo↗

PubMed2025

Comparative analysis of DREB gene family in buckwheat: the role of FtDREB02 in the delphinidin biosynthesis and drought stress response.

Dehydration response element binding (DREB) transcription factors play a pivotal role in plant abiotic stress responses, but its evolutionary and functional characterization in buckwheat remains unexplored. Here, we conducted a comprehensive analysis of the DREB gene family across three buckwheat species, revealing segmental duplication as the primary driver of family expansion and potential purifying selection during evolution. A FtDREB02 gene, classified as group A2, was identified through genome-wide association analysis (GWAS) on drought tolerance and delphinidin content. Functional validation in Arabidopsis thaliana and the hairy root of Tartary buckwheat (Fagopyrum tataricum) demonstrated that overexpression of this gene promotes delphinidin biosynthesis and enhances plant resistance to water scarcity. Through the integration of DAP-seq and PEG transcriptome cluster analysis, a FtANS candidate was screened. Functional studies showed that FtDREB02 regulates delphinidin content by binding directly to DRE elements of the FtANS promoter. This research identifies and comprehensively analyzes the DREB family within buckwheat species, elucidating the regulatory mechanisms of FtDREB02 in controlling flavonoid biosynthesis and drought resistance, providing potential genetic resources for breeding buckwheat varieties with excellent agronomic traits.

Anthocyanins↗

PubMed2003

Microarray analysis of E2Fa-DPa-overexpressing plants uncovers a cross-talking genetic network between DNA replication and nitrogen assimilation.

Previously we have shown that overexpression of the heterodimeric E2Fa-DPa transcription factor in Arabidopsis thaliana results in ectopic cell division, increased endoreduplication, and an early arrest in development. To gain a better insight into the phenotypic behavior of E2Fa-DPa transgenic plants and to identify E2Fa-DPa target genes, a transcriptomic microarray analysis was performed. Out of 4,390 unique genes, a total of 188 had a twofold or more up- (84) or down-regulated (104) expression level in E2Fa-DPa transgenic plants compared to wild-type lines. Detailed promoter analysis allowed the identification of novel E2Fa-DPa target genes, mainly involved in DNA replication. Secondarily induced genes encoded proteins involved in cell wall biosynthesis, transcription and signal transduction or had an unknown function. A large number of metabolic genes were modified as well, among which, surprisingly, many genes were involved in nitrate assimilation. Our data suggest that the growth arrest observed upon E2Fa-DPa overexpression results at least partly from a nitrogen drain to the nucleotide synthesis pathway, causing decreased synthesis of other nitrogen compounds, such as amino acids and storage proteins.

Arabidopsis↗

PubMed2026

Single-cell transcriptomics reveals heterogeneous stress responses and Mg2+-mediated survival mechanisms in Lactobacillus delbrueckii subsp. bulgaricus during freeze-drying and storage.

Maintaining the viability of lactic acid bacteria during dehydration and subsequent storage remains a significant challenge. Here, we employed single-cell RNA sequencing to reveal the heterogeneous stress responses of Lactobacillus delbrueckii subsp. bulgaricus, identifying seven distinct transcriptional clusters across the liquid culture, freeze-drying, and storage phases. The dominant clusters in the freeze-drying and storage were not completely consistent, showing significant functional differentiation. Genomic stability may be important for survival during freeze-drying and storage, while intracellular energy homeostasis appears important for viability during storage. The magnesium transporter mgtB was highly expressed in clusters tolerant to freeze-drying and storage, suggesting a critical role for Mg2+ homeostasis. Further experimental validation confirmed that Mg2+ treatment significantly bolstered stress resistance, increasing immediate post-freeze-drying survival by over 2-fold (up to 92.90%) and post-storage survival by over 5-fold (up to 5.98%). Proteomic data indicated that Mg2+ supplementation correlated with the maintenance of several biological functions potentially relevant to bacterial survival during freeze-drying and storage, including DNA repair, translation, and central carbon metabolism. These findings provide a map of microbial stress resistance through population heterogeneity and offer a potential strategy that may be adapted for enhancing the stability of other industrial lactic acid bacteria products.

Freeze Drying↗

PubMed2026

Multilevel genomic, transcriptomic, and epidemiologic evidence linking diabetic retinopathy to Alzheimer disease.

BACKGROUND: Diabetic retinopathy (DR) and Alzheimer disease (AD) share metabolic and vascular dysfunctions, but the extent to which they reflect overlapping genetic susceptibility and neurovascular-metabolic regulatory pathways remains unclear. We combined multi-omics analyses with population-based data to examine the genetic convergence, cellular pathways, and longitudinal association between DR and AD. METHODS: We performed a two-sample Mendelian randomisation (MR) to estimate the association between genetically predicted DR liability and AD risk. We used Bayesian colocalisation analysis to identify shared genomic loci, and summary-data-based MR (SMR) to detect expression-mediated genes jointly associated with DR and AD. We analysed single-cell RNA sequencing data to characterise shared cellular features and related biological pathways. We also conducted an MR-based mediation analysis to explore whether lipid-related, metabolic, or inflammatory traits mediated the observed DR-AD association, and a longitudinal analysis of the UK Biobank cohort to assess the association between DR and incident AD. RESULTS: With the MR analysis, we found that genetically predicted liability to DR was associated with a modest increase in AD risk. Colocalisation analysis supported a shared genetic signal. We identified three genes with shared expression-mediated associations across DR and AD through SMR. Functional enrichment analyses revealed partially overlapping neurovascular and metabolic pathways. Using MR-based mediation analysis, we found no significant intermediary traits linking DR and AD. Findings from the UK Biobank cohort were directionally consistent with the genetic analyses. CONCLUSIONS: Genetic liability to DR is associated with an increased risk of AD and is accompanied by shared expression-mediated effects and convergent neurovascular-metabolic pathways. These findings support the possibility that DR may serve as a clinically accessible indicator of increased neurodegenerative vulnerability.

Humans↗

PubMed2026

Repeated emergence and fitness heterogeneity of KPC-33 in ST11 Klebsiella pneumoniae under ceftazidime-avibactam pressure.

Ceftazidime-avibactam (CZA) is an important therapeutic option for infections caused by Klebsiella pneumoniae carbapenemase (KPC)-producing Klebsiella pneumoniae. However, CZA exposure also selects for emergent KPC variants. Their in vivo evolutionary patterns, fitness consequences, and underlying molecular mechanisms remain unclear. We performed a longitudinal multiomics analysis of 35 clonally related ST11 KPC-producing K. pneumoniae isolates collected from eight hospitalized patients during clinical follow-up, most of whom had received CZA therapy. Whole-genome sequencing, antimicrobial susceptibility testing, in vitro competition assays, enzyme kinetic analysis, and transcriptomic sequencing were used to systematically characterize the within-host evolutionary dynamics of KPC variants and the fitness heterogeneity of KPC-33. Multiple KPC variants were identified during longitudinal follow-up, among which KPC-33 was the most frequently detected. Among the seven patients who received CZA treatment, KPC-33 was detected in longitudinal isolates from four patients. It was also identified in patient P3, who had not received CZA, whereas other variants were only sporadically identified. Biochemical analysis showed that KPC-33 exhibited an altered kinetic profile relative to KPC-2, characterized by reduced catalytic turnover and altered substrate affinity. KPC-33 did not exhibit a uniform and pronounced fitness defect but instead showed marked strain-dependent heterogeneity. Strains with higher competitive fitness generally showed only limited transcriptional changes, whereas those with lower fitness were accompanied by broader transcriptional remodeling. In this longitudinal cohort, KPC-33 was repeatedly detected, predominantly under CZA-associated selective conditions. Its fitness consequences were clearly strain background dependent and may be associated with the extent of transcriptional remodeling. These findings provide new evidence for understanding the in vivo evolution of CZA resistance.

KPC-33↗

PubMed2001

Parallel and comparative analysis of the proteome and transcriptome of sorbic acid-stressed Saccharomyces cerevisiae.

Exposure of Saccharomyces cerevisiae to 0.9 mM sorbic acid at pH 4.5 resulted in the upregulation of 10 proteins; Hsp42, Atp2, Hsp26, Ssa1 or Ssa2, Ssb1 or Ssb2, Ssc1, Ssa4, Ach1, Zwf1 and Tdh1; and the downregulation of three proteins; Ade16, Adh3 and Eno2. In parallel, of 6144 ORFs, 94 (1.53%) showed greater than a 1.4-fold increase in transcript level after exposure to sorbic acid and five of these were increased greater than two-fold; MFA1, AGA2, HSP26, SIP18 and YDR533C. Similarly, of 6144 ORFs, 72 (1.17%) showed greater than a 1.4-fold decrease in transcript level and only one of these, PCK1, was decreased greater than two-fold Functional categories of genes that were induced by sorbic acid stress included cell stress (particularly oxidative stress), transposon function, mating response and energy generation. We found that proteomic analysis yielded distinct information from transcript analysis. Only the upregulation of Hsp26 was detected by both methods. Subsequently, we demonstrated that a deletion mutant of Hsp26 was sensitive to sorbic acid. Thus, the induction of Hsp26, which occurs during adaptation to sorbic acid, confers resistance to the inhibitory effects of this compound.

DNA, Complementary↗

PubMed2025

PLAID: ultrafast single-sample gene set enrichment scoring.

SUMMARY: In recent years, computational methods have emerged that calculate enrichment of gene signatures within individual samples. These signatures offer critical insights into the coordinated activity of functionally related genes, proteins or metabolites, enabling the identification of unique molecular profiles in individual cells and patients. This strategy is pivotal for patient stratification and advancement of personalized medicine. However, the rise of large-scale datasets, including single-cell profiles and population biobanks, has exposed significant computational inefficiencies in existing methods. Current methods often demand excessive runtime and memory resources, becoming impractical for large datasets. Overcoming these limitations is a focus of current efforts by bioinformatics teams in academia and the pharmaceutical industry, as essential to support basic and clinical biomedical research. To address this critical need, we developed PLAID (Pathway Level Average Intensity Detection), an ultrafast and memory optimized single sample gene set enrichment algorithm that utilizes sparse matrix computation. PLAID delivers highly accurate gene set scoring and surpasses the performance of current methods in single-cell and bulk transcriptomics, and proteomics data. PLAID uniquely integrates the most widely used gene set scoring algorithms, enabling researchers to apply multiple methods for cross-validation with outstanding runtime efficiency and minimal memory requirement. AVAILABILITY AND IMPLEMENTATION: PLAID is implemented in the R language for statistical computing. PLAID source code and installation instructions are available with no restrictions at https://github.com/bigomics/plaid.

Algorithms↗

PubMed2025

Blood phenylalanine lowering partially reverses white matter changes in a mouse model of phenylketonuria.

Phenylketonuria (PKU) is a genetic defect caused by lack of the liver enzyme phenylalanine hydroxylase (PAH). This deficiency results in elevated blood phenylalanine (Phe) levels and neurotoxicity, which is manifested by reduced brain size, lower neurotransmitter levels, and reduced myelination. The goal of this study was to investigate brain myelination defects and their reversibility upon blood Phe lowering by analyzing the corpus callosum (CC) of adult Pahenu2 (PAH-deficient) mice. MRI and immunostaining demonstrated a significant reduction in CC volume in Pahenu2 mice. Treatment with an adeno-associated vector (AAV) encoding mouse PAH for 3.5 months improved but did not completely normalize CC volume. Total cholesterol, a major component of myelin, was unchanged in the CC of Pahenu2 mouse, while some sterol intermediates were significantly reduced by treatment. Single-nuclei transcriptomics showed an upregulation of oxidative stress-related pathways and increased expression of transthyretin, ApoE, Cst3, and Cd81 in CC in Pahenu2 mice. Normalization of blood Phe restored gene expression to levels comparable to those of heterozygous mice and was associated with the generation of differentiated myelin-producing oligodendrocyte subtypes and neuroprotective astrocytes. In summary, Pahenu2 mice showed white matter abnormalities and changes in transcriptome and sterol profiles, which were partially corrected by the normalization of blood Phe.

Animals↗

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↗