Search PubMed⌕ Search

SEARCH · Search PubMed

Results for “transcriptomic”

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.

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

At least 289 records · Page 16Linked to original sources
PubMed2024

Species-wide quantitative transcriptomes and proteomes reveal distinct genetic control of gene expression variation in yeast.

Gene expression varies between individuals and corresponds to a key step linking genotypes to phenotypes. However, our knowledge regarding the species-wide genetic control of protein abundance, including its dependency on transcript levels, is very limited. Here, we have determined quantitative proteomes of a large population of 942 diverse natural Saccharomyces cerevisiae yeast isolates. We found that mRNA and protein abundances are weakly correlated at the population gene level. While the protein coexpression network recapitulates major biological functions, differential expression patterns reveal proteomic signatures related to specific populations. Comprehensive genetic association analyses highlight that genetic variants associated with variation in protein (pQTL) and transcript (eQTL) levels poorly overlap (3%). Our results demonstrate that transcriptome and proteome are governed by distinct genetic bases, likely explained by protein turnover. It also highlights the importance of integrating these different levels of gene expression to better understand the genotype-phenotype relationship.

Saccharomyces cerevisiae↗
PubMed2026

Changes in the transcriptome and synthetic lethal dependencies following KRAS mutant expression reveal profound tissue specificity.

Oncogenic KRAS mutations exhibit a striking tissue-restricted tropism, occurring with high frequency in pancreatic, colorectal, and lung adenocarcinomas while remaining rare in other lineages. The molecular basis for why these specific tissues are uniquely permissive to KRAS transformation, and how this context shapes therapeutic vulnerabilities, remains poorly defined. Here, we utilized CRISPR-mediated genome engineering to generate endogenous, conditional KRAS-mutant isogenic cell line models across three primary permissive lineages (lung, colon, and pancreas) and the nonpermissive breast lineage. Integrated genome-wide CRISPR fitness screens and comparative transcriptome analyses revealed that KRAS-driven synthetic lethal (SL) dependencies are profoundly shaped by their tissue of origin. Strikingly, we observed minimal overlap in SL hits across lineages, with only three genes shared among the permissive lines, suggesting that the KRAS oncogene operates through divergent, context-specific genetic networks. Mechanistically, we show that KRAS activation induces a universal MYC-driven metabolic signature, but the specific machinery required to sustain this state is lineage-restricted. We identified a dependency on the diphthamide synthesis pathway to maintain translational fidelity amid a KRAS-induced hypertranslational state. These findings demonstrate that even when driven by the same oncogene, tumors exhibit distinct regulatory landscapes and unique genetic vulnerabilities. Our results provide a framework for developing lineage-aware therapeutic strategies, moving beyond universal KRAS inhibition toward targeted interventions tailored to a tumor's specific tissue context.

Proto-Oncogene Proteins p21(ras)↗
PubMed2026

Ionizing radiation induces bidirectional transcriptomic reprogramming and dynamic NOS2/TREM2 regulation in triple-negative breast cancer cells.

PURPOSE: To characterize irradiation-associated transcriptomic changes in murine triple-negative breast cancer cells and examine dose- and time-response patterns of selected radiation-responsive candidates. MATERIALS AND METHODS: RNA sequencing (RNA-seq) was performed in 4T1 cells collected 24 h after 4 Gy irradiation, followed by Reactome and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment and gene set enrichment analyses. Representative RNA-seq-derived genes were examined by reverse transcription quantitative PCR (RT-qPCR), and selected immune- and inflammation-related transcripts were further assessed across additional radiation doses and post-irradiation time points. Inducible nitric oxide synthase (NOS2) and triggering receptor expressed on myeloid cells 2 (TREM2) protein abundance was assessed by Western blotting, and nitrite accumulation in culture supernatants was measured using a Griess reagent-based assay as an indirect readout of nitric oxide production. RESULTS: RNA sequencing identified 757 differentially expressed genes, including 285 upregulated and 472 downregulated genes. Irradiation was associated with enrichment of inflammatory, interferon-related, immune-system, and cell-adhesion transcriptional signatures, whereas downregulated genes were enriched in cell-cycle-, chromosome-cohesion-, DNA-damage-response-, DNA-repair-, and SUMOylation-related pathways. Selected immune- and inflammation-related transcripts showed distinct temporal patterns. Nos2 mRNA increased across the examined 0-6 Gy dose range and at later post-irradiation time points, whereas NOS2 protein showed different kinetics, with an early peak after 4 Gy irradiation and no clear further increase above 6 Gy. Nitrite accumulation increased after irradiation. Trem2 showed the largest fold increase among strongly upregulated transcripts identified by RNA-seq, but RT-qPCR detected a significant increase only at 24 h, and TREM2 protein abundance remained unchanged across the examined doses and time points. CONCLUSIONS: Ionizing radiation was associated with broad bidirectional transcriptional remodeling in 4T1 cells, involving immune-, inflammatory-, and interferon-related signatures together with reduced representation of cell-cycle- and DNA-repair-related gene sets. The discordant mRNA and protein patterns of NOS2 and TREM2 indicate that transcript-level responses do not necessarily translate into corresponding protein-level changes. These findings define irradiation-associated molecular responses requiring further functional investigation.

Triple-negative breast cancer↗
PubMed2026

Integrated transcriptomic and functional characterization of Claudin-1 reveals its oncogenic and immunomodulatory roles in pancreatic ductal adenocarcinoma.

Pancreatic ductal adenocarcinoma (PDAC) remains among the deadliest malignancies, driven by its invasive nature and lack of effective biomarkers. Disruption of the epithelial barrier, mediated by tight junction components, is a critical yet underexplored contributor to PDAC progression. Claudins, integral regulators of tight junction integrity, display altered expression across cancers, but their prognostic and immunomodulatory roles in PDAC remain unclear. We performed an integrative analysis of 177 RNA-Seq datasets from TCGA and GTEx to characterize Claudin family alterations in PDAC. Differential expression, copy number variation, methylation, and co-expression networks were analyzed alongside clinical and survival data. Prognostic significance was assessed using Kaplan - Meier and Cox regression analyses, while immune cell infiltration was examined using deconvolution algorithms. Functional validation of Claudin-1 was conducted in Capan-1 cells using CRISPR/Cas9 knockout, followed by proliferation, wound-healing, and Western blot assays. Ten Claudin genes were significantly dysregulated, with Claudin-1 and Claudin-4 frequently amplified and associated with advanced stage and poor survival. High Claudin-1 expression correlated with reduced immune infiltration, indicating an immune-excluded phenotype characterized by immune cells retained in the tumor stroma but largely absent from the tumor parenchyma. Claudin-1 knockout markedly inhibited proliferation, migration, and EMT, evidenced by downregulation of Snail and Slug and restoration of E-cadherin expression. This integrative transcriptomic and functional study identifies Claudin-1 as a key driver of PDAC aggressiveness and immune modulation. These findings establish Claudin-1 as a promising prognostic biomarker and therapeutic target for restoring epithelial integrity and counteracting immune evasion in pancreatic cancer.

Humans↗
PubMed2026

Integrating single-cell transcriptomics to construct an oncogene-driven prognostic model and elucidate metabolic-immune crosstalk in hepatocellular carcinoma.

Hepatocellular carcinoma (HCC) is a leading cause of cancer-related deaths, its progression and treatment heterogeneity are mainly influenced by driver gene and tumor micro-environment (TME) interactions. Nevertheless, the mechanisms of this process at the single-cell level remain unclear. This study integrated TCGA and multi-center single-cell transcriptome data to identify a 575 genes HCC-specific core set, developing a single-cell "oncogene scoring" system to quantify individual carcinogenic activity. This score is significantly elevated in malignant and proliferative T cells and is closely associated with metabolic reprogramming, aberrant cell‒cell communication, and immunosuppressive phenotypes. Based on these characteristics, we constructed a machine learning-based Random Survival Forest (RSF) prognostic model validated in multiple independent cohorts, which classifies patients into distinct risk subtypes. The high-risk group exhibits genomic instability, increased tumor stemness, and immune evasion, while the low-risk group was more sensitive to drugs such as sorafenib. This study highlights the potential pathways by which high oncogenic activity is associated with HCC progression, suggesting a profound link with single-cell metabolic‒immune crosstalk. The constructed RSF model offers a promising computational framework for risk stratification and provides hypothesis-generating insights that may inform future personalized treatment strategies for HCC patients.

Hepatocellular carcinoma↗
PubMed2026

SpatialRNA: a Python package for easy application of Graph Neural Network models on single-molecule spatial transcriptomics dataset.

SUMMARY: Image-based spatial transcriptomics (iST) deliver gene expression measurements of RNA transcripts in tissue slices with single-molecule resolution and spatial context preserved. Modern Graph Neural Network (GNN) models are promising methods for capturing the complex molecular and cellular phenotypes in tissues at single-transcript and single-cell levels. A key application of GNNs is the detection of spatial domains or niches, that is, groups of molecules and/or cells that collaboratively work together to produce complex phenotypes. Due to the vast number of detected transcripts in (iST) dataset, applying GNNs on RNA molecule graphs is not trivial. We present a Python package, SpatialRNA, for easy (sub)graph generation from tissue samples and provide comprehensive tutorials for convenient and efficient application of Graph Neural Network models under the PyG framework. This highly scalable tool comprehensively segments tissue into spatial domains, aiding in biological interpretation of iST data and its underlying molecular microenvironments. AVAILABILITY AND IMPLEMENTATION: The SpatialRNA package is freely accessible from online repository https://github.com/ruqianl/spatialrna and can be installed via pip. Comprehensive tutorials, guidance on parameter selection, and complete workflows of case studies are available from the documentation website https://ruqianl.github.io/spatialrna_docs/, and uploaded on Zenodo with a DOI 10.5281/zenodo.17339575.

Neural Networks, Computer↗
PubMed2026

stDyer-image improves clustering analysis of spatially resolved transcriptomics and proteomics with morphological images.

MOTIVATION: Spatially resolved transcriptomics (SRT) and spatially resolved proteomics (SRP) data enable the study of gene expression and protein abundances within their precise spatial and cellular contexts in tissues. Certain SRT and SRP technologies also capture corresponding morphology images, adding another layer of valuable information. However, few existing methods developed for SRT data effectively leverage these supplementary images to enhance clustering performance. RESULTS: Here, we introduce stDyer-image, an end-to-end deep learning framework designed for clustering for SRT and SRP datasets with images. Unlike existing methods that utilize images to complement gene expression data, stDyer-image directly links image features to cluster labels. This approach draws inspiration from pathologists, who can visually identify specific cell types or tumor regions from morphological images without relying on gene expression or protein abundances. Benchmarks against state-of-the-art tools demonstrate that stDyer-image achieves superior performance in clustering. Moreover, it is capable of handling large-scale datasets across diverse technologies, making it a versatile and powerful tool for spatial omics analysis. AVAILABILITY AND IMPLEMENTATION: The source code of stDyer-image and detailed tutorials are available at https://github.com/ericcombiolab/stDyer-image.

Proteomics↗
PubMed2026

Identification of immune cell type-specific susceptibility genes in multiple cancers using transcriptome-wide association studies.

BACKGROUND: Transcriptome-wide association studies (TWAS) integrate gene expression and genome-wide association studies (GWAS) to identify disease susceptibility genes. Because gene expression varies substantially across cell types within tissues, cell type-specific prediction models may enhance the power of TWAS. METHODS: We conducted cell type-specific TWAS leveraging single-cell RNA sequencing data from the OneK1K cohort (14 immune cell types, 1.27 million cells) and GWAS summary statistics for 7 cancers (>290 000 cases in total). To improve prediction accuracy, we developed a modeling framework that incorporates shared gene expression effects across cell types. RESULTS: At a false discovery rate of 5%, we identified 106 (Bonferroni 5%: 13) previously unreported loci for breast cancer, 51 (4) loci for prostate cancer, 11 (4) loci for lung cancer, 39 (5) loci for melanoma, 9 (1) loci for ovarian cancer, and 2 (1) loci for diffuse large B-cell lymphoma, with most genes exhibiting cell type specificity. Gene set analyses confirmed joint associations of unreported genes with breast and prostate cancer risk in UK Biobank data. Additional lung tissue single-cell RNA sequencing data with 113 individuals validated 18 of 32 (56.3%) statistically significant genes for lung cancer. Across cancers, 139 statistically significant genes were shared by at least 2 cancer types and were primarily enriched in specific immune cell types. CONCLUSION: Cell type-specific TWAS improve the identification of novel cancer susceptibility loci and provide insights into the immune landscape of cancer etiology.

Humans↗
PubMed2026

Precision-Based Filtering Facilitates Cross-Referencing of Conventional and Single-Nucleus Transcriptomes to Identify Time- and Temperature-Sensitive Cell Populations.

Transcriptome analysis via RNA sequencing (RNAseq) has become a ubiquitous method of molecular characterization from whole organisms, dissected tissues, and single cells. These experiments continue to provide an extraordinary volume of data describing molecular states and responses to many conditions. However, standard approaches to RNAseq analysis commonly use expression level filters that eliminate potentially useful data in the service of decreasing noise. Here we describe the implementation of a coefficient of variation-based filter for RNAseq gene expression data. This filter prioritizes consistent data across replicates, allowing lowly-expressed genes with low-variation measurements to be retained for downstream analysis. We show, using two independent Arabidopsis RNAseq datasets, that this filter allows for the inclusion of many more transcription factors than even a low-stringency expression level filter. This effect is independent of sequencing depth. We find that these lowly-expressed genes mark specific cell clusters in our single-nucleus (sn)RNAseq dataset and may facilitate future characterization of currently unknown cell types or states. We further characterize communities of co-expressed genes, sampled across the day at two growth temperatures, in relation to snRNAseq cell clusters, finding evidence for a highly photosynthetic cell population, and a cell state marked by high cell division and translation. These methods can be expanded to RNAseq analysis in many systems, facilitating the construction of more detailed models of tissue-specific gene regulatory networks.

Transcriptome analysis↗
PubMed2026

Identifying JAK2 and ANXA5 as Key Genes Linking Obstructive Sleep Apnea and Oxidative Stress via Machine Learning and Multilayer Transcriptomic Integration With Functional Validation.

Obstructive sleep apnea (OSA) is a common and severe sleep disorder closely associated with oxidative stress (OS). This study aims to identify and validate potential OS-related genes associated with OSA through bioinformatics methods. We successfully identified OS-related differentially expressed genes (OS-DEGs) by combining the limma test, weighted correlation network analysis (WGCNA), and OS-related genes from the GeneCards database. Key genes and potential biological roles were further identified using Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG), enrichment analysis, protein-protein interaction (PPI) network analysis, Lasso regression analysis, random forest algorithm, and support vector machine recursive feature elimination (SVM-RFE) method. Evaluate and validate the accuracy of key genes through receiver operating characteristic (ROC) curve analysis. The human single-cell RNA sequencing (scRNA-seq) dataset is used for cell classification annotation, analysis of key gene single-cell expression profiles, and virtual gene knockout experiments based on the scTenifoldKnk algorithm. Integrating scRNA-seq sequencing, pseudotime trajectory inference, cell-cell communication analysis, and bulk immune infiltration deconvolution reveals monocyte subtype remodeling in OSA. Finally, the expression levels of key genes in clinical samples were validated using real-time quantitative PCR (RT-qPCR) and Western blotting. A total of 57 common DEGs, indicating significant enrichment in OS, inflammation, and tumor pathways, particularly prominent in the immunometabolism pathway. By integrating DEGs, WGCNA, PPI results, and machine learning methods, key genes Janus kinase 2 (JAK2) and ANXA5 were screened out. JAK2 was significantly upregulated under disease conditions, while ANXA5 was significantly downregulated. ROC curve exhibited high accuracy (area under the curve [AUC] > 0.85). Human scRNA-seq analysis revealed that key genes were predominantly highly expressed in monocytes. Virtual knockout experiments demonstrated that these key genes play a crucial role in regulating immune responses and inflammatory reactions. PPI networks and enrichment analysis verified that downstream genes S100P, ALOX5AP, PROK2, and PADI4 may collaboratively participate in immune response and inflammation regulation. Finally, clinical sample experiment further validated the results of bioinformatics analysis. This study provides new research insights for the diagnosis, mechanism research, and treatment development of OSA in the future by integrating multilayer transcriptomic and machine learning techniques.

Humans↗
PubMed2025

Inflammatory pathways and immune dysregulation in pediatric postoperative septic shock: A study integrating transcriptomics, machine learning and molecular docking.

This study elucidates the molecular and immune regulatory mechanisms of pediatric postoperative septic shock. Transcriptomic data were obtained from the Gene Expression Omnibus database. Differentially expressed genes were identified using the limma package, and gene co-expression modules were constructed using Weighted Gene Co-expression Network Analysis. Functional enrichment was performed via gene set enrichment analysis, Gene Ontology, and Kyoto Encyclopedia of Genes and Genomes analyses. Immune cell infiltration was assessed using ESTIMATE and CIBERSORT. Mendelian randomization was applied to explore causal relationships between gene expression and septic shock. Feature genes were selected using machine learning algorithms, and a diagnostic nomogram model was constructed. Finally, molecular docking analysis was performed to screen and evaluate the binding affinity of traditional Chinese medicine monomers to core target proteins. A total of 1331 differentially expressed genes were identified, and the turquoise module was strongly correlated with septic shock. Enrichment analysis revealed significant activation of IL-6/JAK/STAT3, TNF-α/NF-κB, and PI3K/Akt/mTOR pathways. Immune infiltration analysis indicated suppressed immune scores and imbalances in neutrophils, macrophages, T cells, and B cells. Mendelian randomization confirmed causal associations for 6 genes, including PIM3. The predictive model based on feature genes demonstrated high diagnostic performance. Molecular docking suggested that quercetin and astramembrannin I could stably bind PIM3. This study systematically identified core genes, dysregulated immune pathways, and candidate small-molecule interventions in pediatric septic shock, providing novel insights for early diagnosis and targeted therapy.

Humans↗
PubMed2024

PoweREST: Statistical Power Estimation for Spatial Transcriptomics Experiments to Detect Differentially Expressed Genes Between Two Conditions.

Recent advancements in Spatial Transcriptomics (ST) have significantly enhanced biological research in various domains. However, the high cost of current ST data generation techniques restricts its application in large-scale population studies. Consequently, there is a pressing need to maximize the use of available resources to achieve robust statistical power. One fundamental question in ST analysis is to detect differentially expressed genes (DEGs) among different conditions using ST data. Such DEG analysis is often performed but the associated power calculation is rarely discussed in the literature. To address this gap, we introduce, PoweREST (https://github.com/lanshui98/PoweREST), a power estimation tool designed to support power calculation of DEG detection with 10X Genomics Visium data. PoweREST enables power estimation both before any ST experiments or after preliminary data are collected, making it suitable for a wide variety of power analyses in ST studies. We also provide a user-friendly, program-free web application (https://lanshui.shinyapps.io/PoweREST/), allowing users to interactively calculate and visualize the study power along with relevant the parameters.

Differentially expressed genes↗
PubMed2026

Spatial Transcriptomics Identifies Characteristic Immunological Niches in Atopic Dermatitis.

BACKGROUND: Atopic dermatitis (AD) is primarily driven by a Type 2 immune response, with T helper (TH2) cells producing IL-4 and IL-13, thereby promoting inflammation, itch, and a compromised skin barrier. Yet, the spatial organization of pathogenic immune cells and their interactions with stromal and epithelial compartments in human AD skin remain incompletely understood. METHODS: We performed 10× Genomics Visium spatial transcriptomics on FFPE skin biopsies from patients with AD (n = 6), psoriasis (n = 2), and healthy controls (n = 5). Data were integrated with AD single-cell RNA sequencing (scRNA-seq) datasets and complemented by imaging mass cytometry (IMC) and multiplex immunofluorescence (IF) to validate the spatial localization of immune cells. Cell-cell communication analysis revealed putative signaling interactions within immune niches. RESULTS: Spatial clustering resolved tissue compartments and demonstrated transcriptional dysregulation in keratinocytes in AD and psoriasis. AD lesions showed a conserved spatial organization of immune aggregates within the superficial dermis. Integration of scRNA-seq signatures revealed spatially organized co-localization of T cells and mature migratory dendritic cells (mmDCs). We developed a ring-based neighborhood analysis to characterize the cellular organization of the immune-stromal niches, revealing T cell-enriched regions surrounded by inflammatory fibroblasts and activated keratinocytes. Intercellular communication analysis further identified putative signaling within mmDC-T cell niches that may promote pathogenic T cell recruitment and activation. Application of tertiary lymphoid structure (TLS) signatures indicated the presence of TLS-like regions. IMC and IF validated the close spatial proximity between activated TH2 cells and mmDCs. CONCLUSION: AD lesions contain spatially organized TLS-like immune niches at the dermal-epidermal junction, characterized by the close association of T cells and mmDCs and coordinated interactions with surrounding stromal and epithelial compartments. These mmDC-T cell niches may represent potential targets for future therapeutic strategies aimed at disrupting persistent local inflammatory pathways and improving long-term disease control.

atopic dermatitis↗
PubMed2026

Dysregulated Sheddase Signalling as a Molecular Driver of Plaque Instability Revealed by Integrative Transcriptomics.

Atherosclerosis is a major cause of mortality due to chronic and progressive low-grade inflammation and fibroproliferative remodelling of the intima of arteries. Comprehensive understanding of the interplay between plaque biology and the mechanisms underlying plaque vulnerability and rupture is essential. Here, we aimed to investigate the transcriptomic profiles of stable and unstable atherosclerotic plaques using RNA sequencing data from human carotid atherosclerotic plaque samples based on next-generation knowledge discovery (NGKD) methods. High-throughput RNA-seq data from plaques dissected in stable and unstable regions of four patients were obtained from the Gene Expression Omnibus (GEO) database. GEO RNA-seq Experiments Interactive Navigator (GREIN) software was used to obtain raw gene-level counts and filtered metadata for this dataset. The data were further filtered and normalized using Express analyst to derive differentially expressed genes (DEGs) in unstable plaques compared to stable plaques. The DEGs were further analysed using WebGestalt, STRING DB, preranked gene set enrichment analysis (GSEA), and Ingenuity Pathway Analysis (IPA) software. We identified 4792 DEGs in unstable plaques based on a p-value cutoff of <&#x2009;0.05. NGKD analysis revealed that the sheddase pathway, collagen degradation, activation of matrix metalloproteinases (MMPs), and extracellular matrix (ECM) degradation ranked among the top five upregulated pathways, whereas the inhibition of MMPs and smooth muscle contraction pathways were identified as the most prominent downregulated pathways in unstable plaques. We found that the sheddase pathway was one of the most significantly upregulated canonical pathways in unstable plaques and this finding opens new avenues for potential therapeutic interventions in patients with atherosclerosis.

Humans↗
PubMed2026

Comparative Transcriptomic Analyses Identify Candidate Genes for Convergent Reproductive Shifts in a Bimodal Viviparous Amphibian.

Shifts in reproductive mode represent key evolutionary innovations that shape species' life histories and evolutionary trajectories. Species showing bimodal reproductive strategies with multiple independent origins offer a rare opportunity to gain insights into the adaptive processes and mechanisms underlying convergent traits. The fire salamander, Salamandra salamandra, is the only amphibian exhibiting intraspecific variation in reproductive mode across multiple independent reproductive shifts, enabling investigation of the transition between larviparity (females give birth to aquatic larvae) and pueriparity (females give birth to fully developed terrestrial juveniles) within a single species and across different timescales. Pueriparity is an adaptive innovation that skips the aquatic larval stage, allowing individuals to exploit habitats with no available water bodies. The fire salamander is larviparous across most of its range, but pueriparity has evolved independently at least three times: once in the early Pleistocene within S. s. bernardezi in the mountains of northern Spain, and more recently on two land-bridge islands (NW Spain) inhabited by S. s. gallaica. To identify candidate genes associated with these distinct reproductive modes, we compared gene expression profiles of the uterus and oviduct of pregnant females across two independent evolutionary transitions using RNA-sequencing. We detected shared changes in maternal gene expression among pueriparous S. s. bernardezi and S. s. gallaica relative to their larviparous counterparts, in addition to differences unique to each independent evolutionary transition. Functional enrichment analyses indicated that differentially expressed genes were associated with reproductive timing, angiogenesis, and maternal signalling, consistent with the phenotypic differences observed in the uterine environment and embryonic development between the two reproductive modes. This study represents an important first step towards understanding the genomic basis of the evolution of pueriparity in a remarkable bimodal reproductive system, and provides transcriptomic resources and candidate genes for future research into the genomic architecture underlying this poorly understood adaptive trait.

Animals↗
PubMed2026

Comparative Transcriptomic Analysis of Human Macrophages During Mycobacterium avium Versus Mycobacterium tuberculosis Infection.

The treatment of Mycobacterium avium (Mav) infection, responsible for over 80% of nontuberculous mycobacterial pulmonary disease, remains challenging due to rising antibiotic resistance and unsatisfactory success rates. Hence, there is a need for a deeper understanding of host-pathogen interactions to inform the development of alternative therapeutic approaches, like host-directed therapy (HDT), aimed at improving host antimycobacterial defenses. However, compared to Mycobacterium tuberculosis (Mtb) infections, knowledge of host-pathogen interactions for Mav infection is still limited. To address this knowledge gap, we performed a genome-wide host transcriptomic analysis of Mav-infected primary human macrophages-the primary host cell-alongside Mtb-infected macrophages to leverage insights from Mtb research. Our findings show substantial overlap in the gene expression patterns between Mav-infected and Mtb-infected macrophages, including induction of cytokine responses and modulation of various G-protein coupled receptors (GPCRs) involved in (lipid-mediated) macrophage immune functions. Notable differences were observed in the expression of immediate early genes (IEGs), phospholipases, and genes of the GTPase of immunity-associated protein (GIMAP) family. This study laid a foundation for identifying both shared and Mav-specific host response pathways, providing direction for future investigations into host-pathogen interactions during Mav infection and the identification of novel targets for HDT.

Humans↗
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&#x2009;Gb, comprising nine anchored chromosomes with a contig N50 of 30.66&#x2009;Mb and a scaffold N50 of 350.58&#x2009;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 &#x3b1;-bisabolol biosynthesis in C. aromaticum, which has specifically expanded in Chrysanthemum species through tandem gene duplications, contributing to the elevated accumulation of &#x3b1;-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

Single-cell transcriptomic atlas of Alzheimer's disease middle temporal gyrus reveals region, cell type, and sex specificity of gene expression with novel genetic risk for MERTK in female.

BackgroundAlzheimer's disease (AD), the most common age-related neurodegenerative disease, is closely associated with both amyloid-&#x3b2; plaque and neuroinflammation. Two thirds of AD patients are female, and they have a higher disease risk; women with AD have more extensive brain histological changes than men along with more severe cognitive symptoms and neurodegeneration.ObjectiveThis study aimed to determine how sex difference induces structural brain changes and molecular cell vulnerabilities in AD, with a focus on identifying sex-specific transcriptional alterations and genetic risk factors.MethodsWe performed single nucleus RNA sequencing on postmortem brains from individuals with AD and age- and sex-matched controls, focusing on the middle temporal gyrus, a cortical brain region strongly affected by the disease, and integrated single nucleus RNA sequencing results with genome-wide association study (GWAS) data using cell type-specific enrichment and generalized gene-set analysis approaches. The analysis pipeline is provided with threshold information.ResultsWe identified a selectively vulnerable subpopulation of layer 2/3 excitatory neurons that were RORB-negative and CDH9-expressing in both males and females. Disease-associated, but sex-independent, reactive astrocyte signatures were also present. In clear contrast, the microglia signatures of AD brains differed between males and females. Integrating single cell transcriptomic data with results from GWAS, we identified MERTK genetic variation as a candidate novel risk factor for AD selectively in females.ConclusionsTaken together, our single cell atlas of middle temporal gyrus revealed a unique cellular-level view of sex-specific transcriptional changes in AD, illuminating GWAS identification of sex-specific AD genes. These data serve as a rich resource for interrogation of the molecular and cellular basis of AD.

Alzheimer's disease↗