SEARCH · PubMed Health

Results for “whole-transcriptome sequencing”

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.

17 recordsLinked to original sources

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

Whole-transcriptome RNA sequencing and ceRNA network analyses provide novel insights into the antibacterial immune response of Hippocampus abdominalis against Vibrio harveyi.

Long non-coding RNAs (lncRNAs) stand as newly-arisen molecular types that exert regulatory effects, able to operate as competitive endogenous RNAs (ceRNAs) to engage microRNAs (miRNAs) in interaction, resulting in the recovery of target mRNA expression and activity. Increasing evidences indicate that the ceRNA network affects various biological processes in mammals, including development, cellular differentiation, metabolism, immune response, and disease pathogenesis. In teleost fish, the lncRNA-miRNA-mRNA regulatory networks have been reported occasionally. However, up to now, the roles of lncRNAs in the big-belly seahorse (Hippocampus abdominalis) remains unclear. In this study, we reported for the first time, via whole-transcriptome RNA sequencing, the lncRNA mediated ceRNA regulatory network in Vibrio harveyi-infected H. abdominalis. A total of 4197 differentially expressed mRNAs (DE-mRNAs), 1317 DE-lncRNAs, and 183 DE-miRNAs were identified. Furthermore, the crosstalk between miRNAs and lncRNAs as well as between miRNAs and mRNAs was inferred based on the negative correlations between miRNAs and their target lncRNAs/mRNAs. A core immune associated lncRNA-miRNA-mRNA putative regulatory network was thus constructed, comprising 211 lncRNA-miRNA and 224 mRNA-miRNA pairs. In conclusion, our findings provide an integrative overview of the ceRNA regulatory networks on the underlying immune responses to V. harveyi infection in the big-belly seahorse, and offer a solid theoretical foundation for the comparative immunological research of teleost fish.

Animals↗

PubMed2026

Integrative Genomic Profiling of Newly Diagnosed Prostate Cancers Progressing on Surveillance.

OBJECTIVE: To identify molecular features associated with earlier progression to definitive therapy amongst patients with localized prostate cancer (PCa) managed on active surveillance (AS). METHODS: We performed a retrospective pilot study of 7 patients with low- to intermediate-risk PCa undergoing serial multiparametric MRI (mpMRI)-targeted biopsies of the same lesion while on AS, who all proceeded to definitive therapy. Time-to-treatment (TTT) was defined as years from first biopsy on AS to definitive therapy. Laser-capture microdissection was used to separate tumor epithelium, benign glands, high-grade prostatic intraepithelial neoplasia, and stroma in each biopsy specimen. DNA from the tumor and matched benign tissue underwent whole-exome sequencing, and RNA from all compartments underwent whole-transcriptome sequencing. Somatic mutations and copy-number alterations were compared across serial biopsies and used to reconstruct phylogenies and quantify clonal complexity. RESULTS: Tumors exhibited substantial intratumoral heterogeneity, and in 3 of 6 paired cases, serial mpMRI-targeted biopsies showed discordant somatic profiles consistent with sampling distinct major clones over time. By contrast, no single gene-level alteration, and few large-scale chromosomal events, were associated with TTT. High clonal complexity, defined as ≥3 subclones, was associated with significantly shorter TTT than low complexity (median 1.9 vs 7.2 years; P = .0082). Exploratory pathway analyses of individual tissue components suggested TTT-associated differences in inflammatory signaling and stromal-epithelial cross-talk. CONCLUSION: In this small, hypothesis-generating cohort, clonal complexity was more closely associated with earlier definitive therapy than individual genomic alterations. Larger prospective studies are needed to validate whether multiomic measures of clonal architecture can improve AS risk stratification.

Humans↗

PubMed2026

Comprehensive Analyses of SOX7 Provide Novel Insights on Its Tumor Suppressor Role and Its Target Genes with Therapeutic Implications in Multiple Myeloma.

Multiple myeloma (MM) is an incurable hematological malignancy. SOX7, located within the recurrently deleted 8p23.1 region in MM, is suggested to act as a tumor suppressor. We characterized SOX7 through genetic, epigenetic, and functional analyses in MM cell lines. SOX7 was frequently silenced due to deletion and/or promoter hypermethylation. Ectopic SOX7 expression in KMS-18 and MM.1S cell lines caused a progressive decline in SOX7-transduced cells and induced G1 cell cycle arrest and/or apoptosis. Although SOX7 re-expression did not enhance bortezomib efficacy, treatment with the pan-histone deacetylase inhibitor panobinostat induced G1 arrest, promoted apoptosis, and increased SOX7 expression in MM.1S cells. Whole-transcriptome sequencing identified G1/S progression-related Wnt/β-catenin pathway genes as major SOX7-regulated targets, while ChIP-Seq analysis revealed widespread genomic SOX7 occupancy in MM.1S. Flow cytometric analysis of permeabilized bone marrow tumor cells from newly diagnosed and relapsed MM patients demonstrated generally low SOX7 protein expression. Collectively, these results indicate that SOX7 functions as a tumor suppressor in MM, and its inactivation promotes cell cycle progression. The anti-myeloma effects of panobinostat in MM.1S cells may be partially mediated through SOX7 induction.

Multiple Myeloma↗

PubMed2026

Genomic Characterization of ETV6::RUNX1-Positive Childhood B-ALL in a Chinese Cohort: Novel Fusion Partners, Co-Occurring Mutations, and Risk-Stratifying Biomarkers.

BACKGROUND: ETV6::RUNX1 is the most common genetic abnormality in pediatric B-cell acute lymphoblastic leukemia (ALL; ∼25%), yet the comprehensive genetic architecture and molecular predictors of intermediate-risk (IR) stratification remain incompletely characterized. METHODS: We performed whole-transcriptome sequencing (Illumina NovaSeq 6000, rRNA depletion, 41.70 Gb/sample) on bone marrow samples from 93 pediatric ETV6::RUNX1-positive B-ALL patients. Bioinformatics analysis included STAR alignment, MuTect2 variant calling, FusionCatcher fusion detection, and VEP annotation. The Jaccard index with permutation testing assessed mutation co-occurrence; logistic regression identified independent predictors of IR classification. RESULTS: Beyond ETV6::RUNX1, we identified 51 distinct fusion genes across the cohort, including the reciprocal RUNX1-ETV6 (73.1%), chr8::KLF1210 (38.7%), and KLF12-chr8 (34.4%). Somatic mutations in 249 genes were detected; the most frequent were KIAA1715 (17.2%), KRAS (11.8%), and NSD2 (10.8%). Network analysis revealed significant chromatin modifier co-occurrence (KIAA1715-KMT2C: J = 0.136, p = 0.015) and KRAS-NRAS mutual exclusivity (J = 0.000, p = 0.042). PTCH1 (OR = 3.50, 95% CI 0.21-58.49, p = 0.41) and GNB1 (OR = 6.5, 95% CI 1.2-34.8, p = 0.029) mutations independently predicted IR classification. chr8::KLF1210 fusion correlated with higher Day-19 MRD levels (p = 0.038). CONCLUSIONS: GNB1 mutation represents a novel independent predictor of IR stratification in ETV6::RUNX1-positive B-ALL. The chromatin modifier co-occurrence module and extensive fusion architecture reveal biological heterogeneity within this favorable-risk subtype, with potential implications for risk-adapted therapeutic strategies.

B‐ALL↗

PubMed2026

Intraductal Papillary Squamous Neoplasm (IPSN) of the Pancreas: Histological and Molecular Characterization of a Novel and Distinct Intraductal Cancer Precursor.

We report 6 intraductal papillary squamous neoplasms (IPSNs) of the pancreas, a rare but distinctive tumor whose biological features remain largely unknown. Five cases were investigated using an integrated approach combining histomorphological evaluation, immunohistochemistry, and multiregional molecular profiling through whole-exome DNA sequencing and whole-transcriptome RNA sequencing. Only targeted DNA sequencing was available on a sixth recently diagnosed case. Histologically, the intraductal lesions were characterized by large, confluent papillae with fibrovascular cores lined by multilayered epithelial cells with diffuse squamous differentiation. All cases harbored a concomitant invasive carcinoma. The associated invasive carcinomas consistently included a pancreatic tubular/ductal adenocarcinoma; in 5 cases, a poorly differentiated squamous cell carcinoma was also present, the proportion/features of which met the diagnostic criteria of adenosquamous carcinoma in 2 of them. Genomic analyses revealed that IPSNs and their matched invasive carcinomas shared the majority of somatic alterations, supporting a shared clonal origin for the 2 components. Activating KRAS mutations and biallelic inactivation of CDKN2A were detected in all cases. Recurrent mutations involved members of the SWI/SNF chromatin-remodeling complex and KMT2D. Additionally, FGFR1 and MYC amplifications were identified in 2 distinct cases (1 case each). Molecular alterations restricted to the invasive component involved mediators of the transforming growth factor-β signaling pathway. Transcriptomic profiling demonstrated a basal-like expression pattern in all IPSNs and squamous cell carcinomas, although in 2 cases, the matched pancreatic tubular/ductal adenocarcinoma shifted toward a classical transcriptomic subtype. In conclusion, through integrated histological assessment and multiregional molecular sequencing, we demonstrate that IPSN represents a bona fide precursor of invasive pancreatic cancer, a new addition to the intraductal neoplasms category. This study challenges the current paradigm that pancreatic squamous epithelium plays no role in the initiation of pancreatic carcinogenesis, providing the first evidence of its involvement in early tumorigenic processes and yielding immediate implications for pancreatic tumor classification and biological understanding.

Humans↗

PubMed2026

Comprehensive analysis of differentially expressed mRNAs, lncRNAs, and miRNAs involved in ovarian differentiation and development in Qihe gibel carp (Carassius gibelio var. Qihe).

Qihe gibel carp (Carassius gibelio var. Qihe) exhibits diverse reproductive modes including gynogenesis and sexual reproduction, yet the molecular mechanisms of ovarian differentiation remain poorly understood. Ovarian tissues at 20, 30, and 60 days after hatching (dah), representing key stages covering early ovarian differentiation and primary oocyte growth, were subjected to whole-transcriptome sequencing. A total of 27,259 mRNAs, 2622 lncRNAs, and 2467 miRNAs were differentially expressed. Cell cycle, transcription, translation, and DNA replication pathways were significantly upregulated from 20 to 60 dah. Oocyte meiosis was enriched from 20 and 30 dah, whereas metabolic pathways (lipid, carbohydrate, and nucleotide metabolism) were enriched from 30 to 60 dah, indicating sequential progression from meiosis initiation to primary oocyte growth with nutrient synthesis. Hub lncRNAs and key ceRNA networks (e.g., MSTRG.28669.5-miR-221-ccnb2) were identified. This study provides the first comprehensive characterization of ncRNA-mediated regulation and ceRNA networks during ovarian development in Qihe gibel carp, establishing a foundation for understanding ovarian differentiation in this species.

Animals↗

PubMed2026

BCL11B enhancer hijacking by t(14;16)(q32;q24) translocation defines a novel high-risk subtype of T-ALL.

The molecular classification of T-cell acute lymphoblastic leukemia (T-ALL) remains incomplete, limiting risk stratification and the development of targeted therapies. Enhancer hijacking is a critical oncogenic mechanism that deregulates proto-oncogenes by repositioning cisregulatory regions via structural variants. Here, we performed an integrated analysis of pediatric and adult T-ALL and mixed-phenotype acute leukemias (MPALs), using whole-genome and whole-transcriptome sequencing. This analysis identified a group of 14 patients with predominantly T-lineage neoplasms driven by a t(14;16)(q32;q24) translocation, harboring universal GATA3 mutations and CDKN2A/B deletions. Mechanistically, this translocation repositions the ThymoD locus downstream of BCL11B, causing monoallelic, ectopic overexpression of FENDRR and mesenchymal transcription factor genes FOXF1 and FOXC2 and activating epithelial-mesenchymal transition transcription signatures. Immunophenotypic and single-cell RNA sequencing analyses revealed marked lineage ambiguity with myeloid and B-cell differentiation potentials specific to this subtype. Furthermore, functional analyses in CD34+ cord blood cells demonstrated that FOXF1 overexpression promotes myeloid differentiation while suppressing T-cell differentiation, serving as a key factor for lineage specification. Clinically, this subtype was detected in 0.15% to 4.0% of T-ALL/MPAL cases depending on the cohort, showing a median age of 15 years and enrichment in adolescents and young adults. Importantly, patients with t(14;16)(q32;q24) have an extremely poor prognosis, showing a trend toward worse outcomes than high-risk groups such as KMT2A-rearranged early T-cell progenitor-like, SPI1-rearranged, and LMO2 γδ-like T-ALLs. The unique molecular landscape and poor prognosis of patients with the t(14;16)(q32;q24) translocation underscore the need for the development of novel subtype-specific therapeutic approaches.

Humans↗

PubMed2026

Screening and identification of the ncRNA-mRNA regulatory network associated with DNA methylation in goose embryonic myoblasts.

BACKGROUND: Local goose breeds Shitou and Wuzong exhibit distinct growth rates, implying divergent embryonic muscle development. This study used embryonic myoblasts from the Magang goose, an established model with superior growth traits, to explore the underlying common regulatory mechanisms. Extending our previous findings that 5-AZA (DNA methylation inhibitor) and BC339 (DNA hydroxylation inhibitor) oppositely affect myoblast proliferation and differentiation, we performed whole-transcriptome sequencing on inhibitor-treated goose embryonic myoblasts. This aimed to identify DNA methylation-mediated ncRNA-mRNA networks governing myoblast fate, with key interactions being functionally validated. RESULT: 5-AZA significantly promotes cell proliferation and differentiation by inhibiting DNA methyltransferase activity and reducing DNA methylation levels, whereas BC339 significantly suppresses cell proliferation and differentiation by inhibiting demethylation and increasing DNA methylation levels. Specifically, we identified 6,309 mRNAs, 579 lncRNAs, 194 miRNAs, and 825 circRNAs that were differentially expressed in response to 5-AZA and BC339 treatment. Based on GO and KEGG enrichment analyses, differentially expressed genes related to muscle development were selected to construct a ceRNA network. This network comprises 5 differentially expressed lncRNAs (DELs: MSTRG.17572.1, XR_001211738.1, MSTRG.1886.1, XR_001212555.1, MSTRG.8995.2), 2 differentially expressed circRNAs (DECs: novel_circ_029953, novel_circ_017636), 11 differentially expressed miRNAs (DEMs: miR-383-x, miR-10174-y, miR-191-x, miR-24-x, miR-9619-y, novel-m0303-5p, novel-m0105-3p, miR-204-x, miR-211-z, novel-m0075, miR-26-y), 5 differentially expressed genes (DEGs: KIF3A, CCND1, PPM1A, Table 2, TGFBR1), forming a total of 24 interactions. This study identified miR-9619-y as a critical negative regulator of goose embryonic myoblast development through targeted inhibition of CCND1. Dual-luciferase reporter assays confirmed the direct binding of miR-9619-y to the 3'-untranslated region of CCND1. Functional experiments demonstrated that overexpression of miR-9619-y significantly reduced the EdU-positive cell ratio and myotube area percentage, accompanied by cell cycle arrest at the G0/G1 phase. Conversely, inhibition of miR-9619-y promoted myoblast proliferation and differentiation while decreasing the proportion of cells in G0/G1 phase. During the proliferation stage, miR-9619-y overexpression significantly suppressed CCND1 expression at both mRNA and protein levels, down-regulated MyoD expression, and reduced Myf5 mRNA abundance; whereas miR-9619-y inhibition up-regulated these genes and their corresponding proteins. During the differentiation stage, overexpression of miR-9619-y similarly decreased the mRNA levels of CCND1, Myh1, and MyoG, as well as the protein levels of MyHC and CCND1, with inhibition producing the opposite effects. CONCLUSION: In this study, we predicted a ceRNA network based on bioinformatics analysis governing goose embryonic myoblast development, identifying key molecular components including mRNAs, miRNAs, lncRNAs, and circRNAs, along with 24 regulatory axes. Functional experiments further demonstrated that miR-9619-y arrests cell cycle progression and negatively regulates the proliferation and differentiation of goose embryonic myoblasts, as evidenced by its impact on both the mRNA and protein expression of key myogenic factors through targeted inhibition of CCND1. These findings, together with the bioinformatically predicted ceRNA network, suggest potential complex post-transcriptional regulatory mechanisms underlying myogenesis in geese and offer candidate molecular targets for genetic improvement of meat production performance in waterfowl breeding programs.

Animals↗

PubMed2026

The maternal-to-zygotic transition is a critical window for PFOA-induced disruption of developmental programming.

Early embryogenesis is governed by precisely timed gene regulatory programs that coordinate cell fate specification, tissue patterning, and morphogenesis. The maternal-to-zygotic transition (MZT) represents a pivotal developmental milestone during which regulatory control shifts from maternally deposited transcripts to activation of the zygotic genome. Disruption of this transition has the potential to alter developmental trajectories with lasting consequences. Per- and polyfluoroalkyl substances (PFAS), environmentally persistent contaminants, have been linked to developmental abnormalities, yet their impact on core embryonic gene regulatory networks especially with exposure during MZT is not well understood. Using zebrafish (Danio rerio), a tractable vertebrate model and New Approach Methodology (NAM), we investigated how PFAS exposure during the MZT alters early developmental programming. Embryos were exposed starting at different times before and within the MZT time window and collected at 24 h post-fertilization (hpf) for transcriptomic analysis. Targeted qRT-PCR revealed dysregulation of genes controlling transcriptional activation, lineage specification, proliferation, and differentiation. Whole-transcriptome RNA sequencing (RNA-seq) further identified widespread perturbations in gene networks governing transcriptional regulation, cell signaling, and embryonic morphogenesis. Temporal analysis revealed that exposure beginning at 3.5 hpf, followed by 8 hpf, corresponding to early zygotic genome activation and near completion of zygotic activation, respectively, resulted in the greatest differential gene expression changes at 24 hpf. Consistent with these early gene regulatory perturbations, larvae exposed starting at 8 hpf also exhibited altered behavior at 5 days post-fertilization. Together, these findings demonstrate that PFAS exposure during MZT disrupts the establishment of embryonic gene regulatory networks, linking environmental toxicant exposure to altered developmental patterning and organismal outcomes. This work underscores the vulnerability of early developmental transitions to environmental perturbation and positions MZT as a critical window of susceptibility during development.

NAMs (new approach methodologies)↗

PubMed2025

Expression regulation network in papillae of sea cucumbers: Whole-transcriptome and DNA methylation datasets.

To elucidate the expression regulation network of papilla size of sea cucumbers (Apostichopus japonicus), the whole-transcriptome and DNA methylome datasets of different sizes of papillae in sea cucumbers were generated. Average clean bases of whole-transcriptome (16.35 G) and DNA methylome (28.92 G) were obtained using RNA sequencing and whole-genome bisulfite sequencing techniques. A total of 3,188 ceRNA networks were also identified including 3,081 long non-coding RNAs (lncRNA)/microRNAs (miRNA)/mRNA networks and 107 circular RNA (circRNA)/miRNA/mRNA networks. Methylome data indicate that there were 3,307 and 3,776 differentially methylated regions (DMRs) with high-level methylation as well as 3,125 and 3,016 DMRs with low-level methylation in big papillae compared to small papillae. The identified DMRs were mainly distributed in introns, promotors, or exons. The whole-transcriptome and DNA methylome datasets generated from this study not only established a robust theoretical foundation (especially from the epigenetic aspect) for elucidating expression regulation network determining papilla size in sea cucumbers but also can be a valuable resource of biomarker mining for papilla appearance-based selective breeding in sea cucumbers.

DNA Methylation↗

PubMed2026

CoxFormer enables spatial omics inference with multimodal generative modeling.

Gene co-expression maps transcriptome-wide gene-gene relationships, yet high-quality estimates cover less than half the genome. Meanwhile, spatial omics either profiles restricted in situ panels or lacks cellular resolution. Extending co-expression transcriptome-wide could overcome these limitations by inferring unassayed gene expression at subcellular resolution. Here we show that CoxFormer integrates literature-derived gene knowledge with co-expression networks from bulk tissues and large-scale single-cell atlases to learn 512-dimensional representations for 32,016 human genes. These embeddings capture functional gene relationships and serve as a generative prior for spatial inference across platforms and modalities. Without requiring a matched single-cell RNA-sequencing reference, CoxFormer supports four applications beyond measured genes: histology-based expression imputation, gene activity prediction from chromatin accessibility, subcellular super-resolution inference, and pathological region detection. Together, CoxFormer extends gene embedding from gene- and cell-level tasks to whole-transcriptome spatial inference, providing a unified framework for biological analysis beyond the limited gene coverage of current spatial omics technologies.

Humans↗

PubMed2026

Parent-of-origin effects on allelic expression bias in interspecific poplar hybrids.

In hybrid plants, phenotypic outcomes are governed by interactions between the two parental genomes. However, the mechanisms underlying the interplay of divergent regulatory networks from these genomes remain poorly understood. In this study, we compared gene-level and allele-specific expression patterns, as well as differentially enriched pathways between F₁ and complex backcross (CBC) lines derived from a natural interspecific hybrid population of Populus fremontii (Pf) and P. angustifolia (Pa). Metabolic differences between Pf and Pa which exhibit low and high levels respectively of phenylpropanoid-derived condensed tannins were leveraged. Using individualized transcriptome references, differential expression and clustering analyses revealed CBC-biased and F₁-biased expression for genes involved in phenylpropanoid metabolism and photosynthesis, respectively. Biased expression of these genes at the allele level was also observed in F1. At the whole-transcriptome level, Pa-biased genes predominated in F₁ hybrids, and Pa alleles displayed more conserved expression patterns than Pf alleles across examined samples. Further analyses indicated that allelic expression bias was significantly associated with parental origin, which could be driven by sequence variations in cis-regulatory elements and differences in CpG island length. Our findings demonstrate strong parent-of-origin effects on divergent regulatory networks governing gene expression in poplar hybrids and provide clues for strategic parental selection tailored to specific metabolic pathways of interest.

cis-regulation↗

PubMed2025

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

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

Humans↗

PubMed2025

Combined somatic mutation and transcriptome analysis reveals region-specific differences in clonal architecture in human cortex.

The human cerebral cortex is specialized into regions, but little is known about how human cellular lineages shape cortical regional variation and neuronal cell-type distribution during development. Here, we map single-cell lineages of human cortical regions and neuronal subtypes using >1,000 somatic single-nucleotide variants (sSNVs) identified from deep bulk whole-genome sequencing and analyzed over 25 regions and >72,000 single cells. In the fronto-parietal cortex, sSNVs are rarely restricted, marking neuron-generating clones that disperse into neighboring regions. In contrast, the primary visual cortex harbors 30%-70% more sSNVs than the neighboring secondary visual cortex. Clones at this border exhibit more restricted dispersion, suggesting late developmental lineage segregation. Single-nucleus sSNV and whole-transcriptome analysis reveal glutamatergic neuron clones with modest regional restrictions that share low-mosaic sSNVs with some GABAergic neurons, suggesting a recent dorsal cortical progenitor. Our analysis reveals human-specific cortical lineage patterns, regional differences in clonal patterns, and late divergence of some glutamatergic/GABAergic lineages.

Humans↗

PubMed2026

Single-cell and spatial transcriptomic technologies for lung cancer tumor microenvironment analysis.

Lung cancer remains one of the leading causes of cancer-related mortality worldwide; beyond its rising incidence, its marked molecular heterogeneity and complex tumor microenvironment (TME) hinder treatment response and drive resistance, contributing directly to its high mortality rate. Single-cell RNA sequencing (scRNA-seq) and spatial transcriptomics (ST) provide complementary approaches for dissecting these features. scRNA-seq enables high-resolution analysis of cellular diversity and transcriptional states but requires tissue dissociation and therefore loses spatial context. In contrast, ST preserves tissue architecture and provides insights into how gene-expression programs within the TME are organized, although no currently available spatial platform combines whole-transcriptome coverage with true single-cell resolution over large tissue areas. Together, these technologies have enabled detailed mapping of tumor, immune and stromal populations and of their spatial interactions, revealing functionally distinct cellular niches that contribute to immune evasion, metastasis and response to therapy. In this narrative review we organize the primary literature around a single question, how spatially structured cellular ecosystems, rather than individual cell types, determine therapeutic response and resistance in lung cancer - and we explicitly separate observations that are reproducible across independent cohorts and platforms from those that remain confined to single studies. We further summarize the technical, analytical and logistic barriers that currently prevent spatially resolved signatures from entering routine diagnostic pathology. Understanding dysregulated pathways and spatially constrained intercellular communication within the TME helps identify candidate biomarkers and may support the identification of therapeutic approaches directed at tumor-intrinsic programs as well as at microenvironment-driven resistance mechanisms.

Cell-cell communication↗

PubMed2026

Whole-transcriptome-scale isoform-resolved spatial imaging of single cells in tissues.

Cell and tissue functions arise from complex interactions among numerous genes, and a systematic understanding of these functions requires isoform-resolved transcriptomic analysis of single cells with high spatial resolution. Here, we introduce an in situ RNA amplification method and its integration with multiplexed error-robust fluorescence in situ hybridization (MERFISH) to detect short RNA sequences and enable whole-transcriptome-scale, isoform-resolved spatial transcriptomics of individual cells in intact tissues. Using this approach, we imaged ∼33,000 distinct RNAs-including ∼23,000 genes and ∼10,000 isoforms-in the mouse brain. Our data enabled systematic analyses of region- and cell-type-specific gene programs and ligand-receptor-based cell-cell communications. These data further revealed rich spatial diversity and cell-type specificity in isoform usage across numerous genes, as well as brain structures particularly rich in isoform specificity. We anticipate broad application of this method for characterizing the molecular and cellular basis of tissue functions, unlocking previously inaccessible discoveries in cell and organismal biology.

Animals↗