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 361 records · Page 20Linked to original sources

PubMed2002

Large-scale analysis of the human and mouse transcriptomes.

High-throughput gene expression profiling has become an important tool for investigating transcriptional activity in a variety of biological samples. To date, the vast majority of these experiments have focused on specific biological processes and perturbations. Here, we have generated and analyzed gene expression from a set of samples spanning a broad range of biological conditions. Specifically, we profiled gene expression from 91 human and mouse samples across a diverse array of tissues, organs, and cell lines. Because these samples predominantly come from the normal physiological state in the human and mouse, this dataset represents a preliminary, but substantial, description of the normal mammalian transcriptome. We have used this dataset to illustrate methods of mining these data, and to reveal insights into molecular and physiological gene function, mechanisms of transcriptional regulation, disease etiology, and comparative genomics. Finally, to allow the scientific community to use this resource, we have built a free and publicly accessible website (http://expression.gnf.org) that integrates data visualization and curation of current gene annotations.

Animals↗

PubMed2026

Molecular and immune profiling of HER2-low, HER2 ultra-low, and HER2-null male breast cancer.

BACKGROUND: HER2 expression is described along a biological continuum from null to positive and serves as a critical biomarker for therapeutic guidance in breast cancer (BC). While HER2-low and ultra-low categories have emerged as actionable targets for antibody-drug conjugates (ADCs) in female BC, their molecular and immune characteristics remain largely unexplored in male breast cancer. METHODS: We profiled 214 male breast tumors using next-generation sequencing and whole-transcriptome sequencing to assess mutational, transcriptomic, and immune landscapes. Tumor mutational burden (TMB) was defined as high if > 10 mutations/Mb. Immune cell fractions were inferred using Quantiseq deconvolution. RESULTS: Among 214 samples, 66 (30.8%) were HER2-null, 53 (24.8%) HER2 ultra-low, 80 (37.4%) HER2-low, and 15 (7.0%) HER2-positive. HER2 ultra-low tumors exhibited a higher prevalence of PIK3CA mutations (39.2% vs 22.6%, p ≤ 0.05) compared to HER2-null. No significant differences were observed in TMB-high frequency or PD-L1 expression across subgroups. Immune composition differed primarily between HER2-null and HER2-expressing subgroups: HER2-ultra-low tumors showed higher B-cell infiltration, whereas HER2-null tumors were enriched in neutrophils. Transcriptomic analysis revealed upregulation of selected stemness-associated genes (NANOG, KLF4, POU5F1) and CEACAM1 in HER2-null tumors, while HER2-low and HER2-ultra-low tumors were largely similar across most molecular and immune readouts in this cohort. CONCLUSIONS: HER2-null male breast cancer appears to represent the most biologically divergent subgroup within the HER2-negative spectrum, whereas HER2-low and HER2-ultra-low tumors were largely similar in this cohort. These findings support further investigation of HER2-null disease as a distinct biological state and provide hypothesis-generating data for biomarker development in this rare population.

Male↗

PubMed2026

Distinct cell morphotypes of Aureobasidium melanogenum ZN exhibit differential functional profiles in promoting maize growth.

Black yeast-like fungi of the genus Aureobasidium exhibit morphological plasticity, but whether distinct cellular states within the same genetic background are associated with different plant growth-promoting functions remains unclear. Here, yeast-like cells (YL), swollen cells (SC), and chlamydospores (CH) of Aureobasidium melanogenum ZN were characterized. YL was associated mainly with siderophore production and laccase activity, SC with extracellular polysaccharide accumulation, and CH with phosphate mobilization and higher ammonia and IAA production. Whole-genome and comparative genomic analyses revealed a shared repertoire related to nutrient acquisition, auxin-associated metabolism, extracellular oxidation, and carbohydrate remodeling, with expansions in nutrient- and cell-surface-related gene families. Transcriptomic and metabolomic analyses showed distinct deployment of these capacities, with CH exhibiting broad reprogramming of tryptophan-associated, nitrogen, phosphate, central-carbon, and amino-acid metabolism. In maize, CH at the optimal inoculation concentration of 105 CFU·mL-1 produced the strongest growth promotion, increasing plant height, dry biomass, root length, root surface area, and root volume by 58.6%, 365.1%, 191.0%, 194.3%, and 222.4%, respectively. Consistent with this pronounced growth phenotype, maize root transcriptomics showed coordinated CH-induced responses involving root development, nutrient transport, redox regulation, and root-interface remodeling. Root-zone tracking showed greater short-term stability and persistence of CH. These findings identify cellular state as an important functional dimension of Aureobasidium-plant interactions and provide a basis for developing fungal inoculants with defined beneficial cellular states.

Zea mays↗

PubMed2026

Blood Bile Acids for Inflammatory Bowel Disease Diagnosis and Disease Activity Assessment: A Metabolomics Meta-Analysis.

Alterations in circulating bile acids (BAs) have been reported in inflammatory bowel disease (IBD), but the consistency of these changes across clinically relevant comparisons remains unclear. Our goal was to investigate systemic BA alterations in IBD using a metabolomics meta-analysis with an exploratory analysis of BA-related gene expression as a supporting context. A systematic review and meta-analysis of 28 metabolomics studies examined blood BA profiles associated with IBD, IBD diagnosis, and disease activity assessment. Univariate analysis and logistic regression modeling of two independent IBD cohorts explored the blood BA-related genes and IBD. Across 28 studies that comprised 5056 IBD patients, 1721 healthy controls, and 314 non-IBD patients, 131 BAs were reported. Eight predefined clinical comparisons were eligible for the meta-analysis. Lower secondary BA levels were consistently observed in IBD patients compared with controls, between UC and CD, and in active versus remission patients. Deoxycholic acid, glycodeoxycholic acid, and taurodeoxycholic acid were frequently decreased, whereas glycocholic acid was increased in certain comparisons. Transcriptomics analyses revealed differential expression of several BA-related genes in blood, including SLC51A, ABCB4, and ACOT8, across the comparisons. Our findings identify consistent circulating BA alterations in IBD and highlight the relevance of blood BA for future biomarker research in the diagnosis and disease activity assessment.

Humans↗

PubMed2001

Molecular neuropathology of transgenic mouse models of Down syndrome.

Down syndrome (DS) is a complex, clinically heterogeneous disorder which shows both impairment of neurodevelopement and the neurodegenerative changes of Alzheimer's disease (AD). The phenotype of DS is caused by triplication of chromosome 21 and transgenic mouse models have been developed, and are being created, that carry single genes and chromosomal segments to excess. For example, transgenic mice containing additional copies of the amyloid precursor protein (APP) gene, have been useful in producing the Abeta deposition characteristic of AD and DS, but not the cytoskeletal changes that are the hallmarks of these human disorders. Such models are useful in replicating aspects of pathogenesis and allow for the testing of therapeutic agents to restore impaired function. Segmental trisomic mouse models, which survive to adulthood and possess three copies of multiple genes responsible for the DS phenotype, such as Ts1Cje and Ts65Dn, have been used to explore aspects of neurodevelopment and neurodegeneration. These animal models show some but not all the pathological, biochemical, and transcriptional changes seen in DS. They also have the advantage of allowing for the testing of therapeutic agents to restore impaired function. Analysis of the transcriptome and proteome of fetal and adult DS indicates that there is a complex relationship between gene dosage, gene and protein expression, and that data from animal models will need to be compared and evaluated in the light of data obtained from DS tissue.

Animals↗

PubMed2026

A critical appraisal of base-resolution m6A profiling techniques.

N6-methyladenosine (m6A) is the most prevalent internal modification in eukaryotic mRNA, influencing RNA fate and gene regulation. Early antibody-based approaches enabled transcriptome-wide profiling but lacked resolution and quantitative accuracy. Newer approaches now achieve base-resolution m6A detection using improved crosslinking, chemical or enzymatic conversion, and single-molecule sequencing. Antibody-free methods provide quantitative stoichiometry from minimal input, while nanopore direct RNA sequencing offers real-time, single-molecule readouts across entire transcriptomes. Collectively, these methods form a versatile toolkit that integrates global mapping with precise site-level analysis, advancing knowledge of context-dependent m6A regulation in physiology and disease. This review compares their principles, strengths, and limitations to guide method selection and highlight how next-generation epitranscriptomic tools are paving the way for clinical and therapeutic applications.

Humans↗

PubMed2026

GRNContext: an interactive web platform for contextualized gene regulatory networks visualization across human cancers.

SUMMARY: While current Gene Regulatory Network (GRN) databases provide comprehensive reference maps of potential interactions between transcription factors and target genes, they do not specify which regulatory interactions are active within specific biological contexts. This limitation is particularly critical in cancer, where transcriptional programs are inherently tissue-specific. To address this gap, we developed GRNContext, an interactive web platform designed for the visualization, exploration, and comparative analysis of gene regulatory networks contextualized across 33 cancer types from The Cancer Genome Atlas (TCGA). Our approach uses the TFLink human reference GRN as a starting point and integrates TCGA transcriptomic profiles to infer cancer-specific regulatory activity. Regulatory relevance was assessed using complementary machine learning and statistical methods, which were unified into a consensus score to prioritize and filter the most relevant candidate regulators for each target gene. By providing both curated context-specific GRNs and a user-friendly platform, GRNContext constitutes a comprehensive and accessible resource that supports mechanistic investigations, hypothesis generation, and translational research focused on transcriptional regulation in cancer. AVAILABILITY AND IMPLEMENTATION: GRNContext is supported by all major browsers and freely available on the web at https://apps.cienciavida.org/grncontext. It is implemented as a client-server web application featuring a FastAPI backend and a React frontend utilizing Cytoscape.js for interactive network visualization, all containerized via Docker for cross-platform compatibility.

Humans↗

PubMed2026

Inheritance of the epigenetic signature and reduced intermuscular bone phenotype acquired via DNA methylation editing of the runx2 b promoter in zebrafish.

The presence of intermuscular bones (IBs) can directly affect the economic value of aquaculture fish. Although genome editing can create IB-free fish by knocking out key IB-related genes, such as runx2b, the associated DNA sequence alterations raise food safety and health concerns, limiting its breeding applications. In this study, we used CRISPR/dCas9-mediated epigenome-editing technology targeting the runx2 b promoter in zebrafish to alter DNA methylation patterns without changing the DNA sequence. Our results showed that higher runx2 b promoter methylation patterns significantly inhibited eGFP mRNA expression levels in the recombinant plasmid. Using the CRISPR/dCas9-Dnmt7 system to enhance methylation of the zebrafish runx2b promoter, we observed a significant decrease in runx2 b mRNA expression levels in the F0 generation. The IBs in the 11 th-16 th muscle segments of the adult F0 fish were significantly shorter compared with the controls. Inbreeding of fish was used to produce F1 and F2 offspring that retained these high promoter methylation levels, along with persistent runx2b expression suppression and IB development inhibition. Transcriptome sequencing analysis suggested that increasing runx2 b promoter methylation levels may synergistically induce additional epigenetic modifications, potentially affecting the PPAR signaling pathway and FoxO transcription factor regulation, which appears to inhibit osteoblast proliferation and differentiation. Overall, this study demonstrates an innovative application of epigenetic editing technology for aquaculture breeding. By precisely regulating the expression patterns of key genes for economically important traits while preserving genomic DNA integrity, this approach provides a theoretical foundation and technical support for improving fish economic traits.

Animals↗

PubMed2026

Genomic, transcriptomic, and molecular predictors of response to neoadjuvant therapy in locally advanced rectal cancer: a narrative review.

Total neoadjuvant therapy (TNT) has emerged as a key treatment paradigm for locally advanced rectal cancer, reducing distant metastasis rates and facilitating organ preservation in selected patients. However, treatment response remains heterogeneous, highlighting the need for biomarkers that can guide treatment selection and optimise outcomes. This narrative review synthesises the current evidence regarding tumour-intrinsic genomic biomarkers associated with response to neoadjuvant therapy, encompassing somatic mutations, germline polymorphisms, gene expression profiles, mismatch repair (MMR) status, protein expression, epigenetic markers, and circulating tumour-derived biomarkers across conventional chemoradiotherapy (CRT) and TNT paradigms. Across the reviewed literature, individual somatic mutations, including KRAS, TP53, and BRAF, demonstrated limited reproducibility as predictive biomarkers, although KRAS mutations were recurrently associated with lower pathological complete response (pCR) rates in CRT-era cohorts. Germline polymorphisms in DNA repair (XRCC1) and folate metabolism (MTHFR) genes showed inconsistent associations with treatment response. In contrast, transcriptomic biomarkers demonstrated greater biological coherence, with proliferative, epithelial-mesenchymal transition, and metabolic signatures frequently associated with treatment resistance, while multi-gene classifiers generally outperformed single-gene markers. Among currently available tumour-intrinsic biomarkers, MMR deficiency was the most consistently reported biomarker associated with reduced response to fluoropyrimidine-based regimens, including TNT, although TNT-specific evidence remains comparatively limited. Dynamic circulating tumour DNA (ctDNA) monitoring, particularly ctDNA clearance during or after therapy, was consistently associated with pathological response and long-term oncologic outcomes across reviewed studies, whereas baseline ctDNA levels showed limited predictive value. Overall, the reviewed literature suggests that biomarker research in rectal cancer has evolved from single-gene analyses towards pathway-level and dynamic biomarkers. The integration of transcriptomic signatures, MMR status, and dynamic ctDNA monitoring may represent a promising strategy for personalising neoadjuvant therapy, improving patient selection for organ-preserving approaches, and enhancing oncologic outcomes in locally advanced rectal cancer. Nevertheless, the evidence base remains heterogeneous, and further prospective validation, assay standardisation, and evaluation within contemporary TNT cohorts are required before these biomarkers can be routinely incorporated into clinical decision-making.

Humans↗

PubMed2026

HIF-1 signaling contributes to lenvatinib resistance in patient-derived HCC organoids.

Resistance to lenvatinib remains an important limitation in hepatocellular carcinoma treatment. Six patient-derived organoid lines were established and classified as sensitive or resistant according to ex vivo drug responses, retaining histological and immunophenotypic features of matched parental tumors. Resistant organoids showed unchanged ATP activity, whereas sensitive ones exhibited pronounced morphological changes and reduced ATP activity at higher concentrations. Transcriptome sequencing identified 408 upregulated and 269 downregulated genes in resistant versus sensitive organoids, with HIF-1 signaling among altered pathways. In resistant organoids, lenvatinib increased HIF-1α, ANGPT2, and HK3 mRNA, whereas comparable changes were not detected in sensitive organoids. KC7F2 reduced these transcripts and further decreased ATP activity when combined with lenvatinib. In organoid-derived xenografts, this combination suppressed tumor growth and HIF-1 target expression more than lenvatinib alone, indicating HIF-1 signaling contributes to the resistant phenotype and its inhibition may enhance response.

Drug resistance↗

PubMed2024

Progressive cardiomyopathy with intercalated disc disorganization in a rat model of Becker dystrophy.

Becker muscular dystrophy (BMD) is an X-linked disorder due to in-frame mutations in the DMD gene, leading to a less abundant and truncated dystrophin. BMD is less common and severe than Duchenne muscular dystrophy (DMD) as well as less investigated. To accelerate the search for innovative treatments, we developed a rat model of BMD by deleting the exons 45-47 of the Dmd gene. Here, we report a functional and histopathological evaluation of these rats during their first year of life, compared to DMD and control littermates. BMD rats exhibit moderate damage to locomotor and diaphragmatic muscles but suffer from a progressive cardiomyopathy. Single nuclei RNA-seq analysis of cardiac samples revealed shared transcriptomic abnormalities in BMD and DMD rats and highlighted an altered end-addressing of TMEM65 and Connexin-43 at the intercalated disc, along with electrocardiographic abnormalities. Our study documents the natural history of a translational preclinical model of BMD and reports a cellular mechanism for the cardiac dysfunction in BMD and DMD offering opportunities to further investigate the organization role of dystrophin in intercellular communication.

Animals↗

PubMed2025

Spatial transcriptomics-aided localization for single-cell transcriptomics with STALocator.

Single-cell RNA-sequencing (scRNA-seq) techniques can measure gene expression at single-cell resolution but lack spatial information. Spatial transcriptomics (ST) techniques simultaneously provide gene expression data and spatial information. However, the data quality of the spatial resolution or gene coverage is still much lower than the quality of the single-cell transcriptomics data. To this end, we develop a ST-Aided Locator for single-cell transcriptomics (STALocator) to localize single cells to corresponding ST data. Applications on simulated data showed that STALocator performed better than other localization methods. When applied to the human brain and squamous cell carcinoma data, STALocator could robustly reconstruct the relative spatial organization of critical cell populations. Moreover, STALocator could enhance gene expression patterns for Slide-seqV2 data and predict genome-wide gene expression data for fluorescence in situ hybridization (FISH) and Xenium data, leading to the identification of more spatially variable genes and more biologically relevant Gene Ontology (GO) terms compared with the raw data. A record of this paper's transparent peer review process is included in the supplemental information.

Single-Cell Analysis↗

PubMed2026

Contrasting regulation of protein-coding genes and lncRNA homeologs in allotetraploid Coffea arabica.

A chromosome-level Bourbon assembly revealed that protein-coding homeologs are predominantly co-regulated between subgenomes. In contrast, intergenic lncRNAs display a modest, but statistically consistent bias toward subgenome E across diverse developmental and stress contexts. Coffea arabica is an allotetraploid species derived from natural hybridization between C. canephora and C. eugenioides, which contributed the C and E subgenomes, respectively. This genomic origin poses major challenges for genome assembly, annotation, and the interpretation of gene regulation. In this study, a high-quality genome assembly of C. arabica was generated and annotated, with particular emphasis on identifying protein-coding genes and intergenic long non-coding RNAs (lincRNAs). Homeologous relationships between genes from the C and E subgenomes were established, providing a robust framework to investigate subgenomic conservation and regulatory divergence. Using an extensive collection of publicly available RNA-seq libraries spanning multiple developmental stages, tissues, and environmental conditions, the relative transcriptional contribution of each subgenome was evaluated. On a global scale, gene expression was largely balanced between subgenomes, with no consistent evidence of subgenome dominance. While protein-coding genes showed comparable regulatory behavior across subgenomes, lincRNAs exhibited a more asymmetric expression pattern, suggesting higher subgenome-specific expression that is interpreted here as a consistent directional tendency rather than as evidence of subgenome dominance. Together, these results provide new insights into the regulatory architecture of the C. arabica genome and establish a foundational genomic and transcriptomic resource for future functional studies and crop improvement efforts.

Coffea↗

PubMed2003

Comparative analyses of potato expressed sequence tag libraries.

The cultivated potato (Solanum tuberosum) shares similar biology with other members of the Solanaceae, yet has features unique within the family, such as modified stems (stolons) that develop into edible tubers. To better understand potato biology, we have undertaken a survey of the potato transcriptome using expressed sequence tags (ESTs) from diverse tissues. A total of 61,940 ESTs were generated from aerial tissues, below-ground tissues, and tissues challenged with the late-blight pathogen (Phytophthora infestans). Clustering and assembly of these ESTs resulted in a total of 19,892 unique sequences with 8,741 tentative consensus sequences and 11,151 singleton ESTs. We were able to identify a putative function for 43.7% of these sequences. A number of sequences (48) were expressed throughout the libraries sampled, representing constitutively expressed sequences. Other sequences (13,068, 21%) were uniquely expressed and were detected only in a single library. Using hierarchal and k means clustering of the EST sequences, we were able to correlate changes in gene expression with major physiological events in potato biology. Using pair-wise comparisons of tuber-related tissues, we were able to associate genes with tuber initiation, dormancy, and sprouting. We also were able to identify a number of characterized as well as novel sequences that were unique to the incompatible interaction of late-blight pathogen, thereby providing a foundation for further understanding the mechanism of resistance.

Cluster Analysis↗

PubMed2026

Phenotypic and transcriptomic characterization of biallelic RNU2-2 developmental and epileptic encephalopathy.

OBJECTIVE: A significant proportion of individuals with suspected genetic developmental and epileptic encephalopathies (DEEs) remain unsolved following whole genome sequencing (WGS). Here we describe biallelic RNU2-2 variants causing a recently reported, severe, recessive DEE. METHODS: We screened individuals who have received WGS analyses at the Genomic Medicine Centre Karolinska for Rare Diseases for biallelic RNU2-2 variants. Deep phenotyping was performed through reviewing entire medical histories and phenotypic traits were transcribed to their corresponding Human Phenotype Ontology (HPO) term. HPO terms were used to generate pairwise phenotypic similarity scores and assess for significantly shared phenotype enrichment in the RNU2-2 sub-cohort. RNA sequencing analyses were performed in fibroblast and blood tissues to compare splicing events between RNU2-2 individuals and two independent control groups. RESULTS: We identified 14 individuals from nine families with 12 ultra-rare biallelic RNU2-2 variants clustering in the conserved 5' domains. Genotype data from 13 of 14 individuals has been reported previously as part of a larger cohort. All individuals presented with a highly concordant, severe DEE, characterized by severe to profound intellectual disability, inability to walk or communicate, hyperkinesia, and refractory seizures. Infantile spasms and tonic seizures were the predominant seizure types and a Lennox-Gastaut syndrome-like phenotype was common. These individuals had a significantly similar phenotypic signature when compared with 703 individuals with complex pediatric epilepsies (two-sided Monte Carlo permutation test, p = .005). RNA sequencing analyses showed aberrant splicing, with the most pronounced effects in fibroblast tissues in mutually exclusive exon and alternate 3' splice-site events, which were not detectable in blood. SIGNIFICANCE: We present deep phenotyping data and transcriptomic analyses that provide support for rare, 5' clustering biallelic RNU2-2 variants causing this novel, severe DEE. We propose an RNA sequencing methodology on fibroblast tissue for future validation of RNU2-2 variants.

autosomal recessive disease↗

PubMed2026

Tree Killer, Qu'est-ce Que C'est? Insights From Forest Pathogen Genomes.

Forests are central to planetary health but are increasingly challenged by emerging diseases driven by climate change, global trade, and anthropogenic disturbance. Despite the apparent resilience of long-lived, genetically diverse tree hosts, forest ecosystems have repeatedly experienced landscape-level pathogen-driven transformations. Advances in genomics, transcriptomics, and functional biology have transformed our understanding of how fungal and oomycete pathogens interact with their hosts across a continuum of lifestyles, from saprotrophy and necrotrophy to biotrophy. Here, we synthesize insights from comparative and population genomics and functional studies across diverse forest pathosystems to examine the traits that characterize successful tree pathogens. We highlight how lifestyle plasticity, adaptations to woody tissues, vector-mediated transmission, and biotrophic stealth enable pathogens to colonize perennial hosts and persist over long temporal scales. We further examine how genome plasticity, hybridization, and horizontal gene transfer generate adaptive potential that often outpaces host evolutionary responses under current environmental change. Finally, we discuss emerging genomic tools, including biosurveillance, machine learning-based classification, and genome editing, that are beginning to link genotype to phenotype and inform assessments of disease risk. By integrating genomic, ecological, and evolutionary perspectives, this review outlines general principles governing forest pathogen success and identifies priorities for future research aimed at improving understanding, early detection, and management of forest diseases in a changing world.

Trees↗

PubMed2003

Discovery of imprinted transcripts in the mouse transcriptome using large-scale expression profiling.

Candidate imprinted transcriptional units in the mouse genome were identified systematically from 27,663 FANTOM2 full-length mouse cDNA clones by expression profiling. Large-scale cDNA microarrays were used to detect differential expression dependent upon chromosomal parent of origin by comparing the mRNA levels in the total tissue of 9.5 dpc parthenogenote and androgenote mouse embryos. Of the FANTOM2 transcripts, 2114 were identified as candidates on the basis of the array data. Of these, 39 mapped to known imprinted regions of the mouse genome, 56 were considered as nonprotein-coding RNAs, and 159 were natural antisense transcripts. The imprinted expression of two transcripts located in the mouse chromosomal region syntenic to the human Prader-Willi syndrome region was confirmed experimentally. We further mapped all candidate imprinted transcripts to the mouse and human genome and were shown in correlation with the imprinting disease loci. These data provide a major resource for understanding the role of imprinting in mammalian inherited traits.

Animals↗

PubMed2026

Multistage Genetic, Transcriptomic, and Single-Cell Evidence Prioritizes MAP1LC3A among Ferroptosis-Related Genes in Glioblastoma.

Glioblastoma (GBM) remains a highly aggressive malignancy, and the contribution of ferroptosis-related genes to disease susceptibility remains incompletely understood. A genetically anchored, multistage framework was applied to prioritize ferroptosis-related genes associated with GBM. Among 483 genes curated from FerrDb V2, 315 had candidate cis-expression quantitative trait loci (cis-eQTLs) in eQTLGen, 250 retained at least three independent instruments after linkage disequilibrium clumping, and 226 yielded valid inverse-variance weighted (IVW) Mendelian randomization estimates using a GBM genome-wide association study comprising 6,183 cases and 18,169 controls. Thirty-four genes met the exploratory discovery criteria of P < 0.05 and a Benjamini-Hochberg false discovery rate (BH-FDR) < 0.20, with directionally concordant Bayesian weighted Mendelian randomization (BWMR) estimates. Replication-stage Mendelian randomization using GTEx V10 whole-blood cis-eQTLs supported four genes: ATG7, RPTOR, MAP1LC3A, and CHMP6. Evaluation across three independent tumor-control transcriptomic cohorts demonstrated that MAP1LC3A was consistently downregulated in tumor tissue and showed a significant random-effects pooled estimate (log₂ fold change, -1.273; 95% confidence interval, -1.625 to -0.920; false discovery rate = 0.016), whereas the other three genes lacked comparable cross-cohort statistical support. Single-cell virtual knockout analysis was subsequently performed in a patient-balanced subset of 2,400 malignant cells selected from 4,916 eligible cells across 20 adult IDH-wild-type GBM tumors. Across five independently seeded runs, 3, 15, 4, and 7 robust downstream genes were identified for ATG7, RPTOR, MAP1LC3A, and CHMP6, respectively. The resulting consensus sets comprised 17 unique genes, with RND3 shared across all four targets. Gene Ontology analysis indicated enrichment of cell-adhesion and cell-surface processes, whereas no KEGG or Reactome pathways remained significant after multiple-testing correction. Collectively, these findings prioritize MAP1LC3A for future experimental investigation while distinguishing genetic association, tumor-expression concordance, and computational perturbation from definitive evidence of causality or mechanism.

Humans↗