SEARCH · Search PubMed

Results for “somatic variants”

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

SourceResearch connectionFrom yearThrough year

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

At least 37 records · Page 2Linked to original sources

PubMed2025

Accurate somatic small variant discovery for multiple sequencing technologies with DeepSomatic.

Somatic variant detection is an integral part of cancer genomics analysis. While most methods have focused on short-read sequencing, long-read technologies offer potential advantages in repeat mapping and variant phasing. We present DeepSomatic, a deep-learning method for detecting somatic small nucleotide variations and insertions and deletions from both short-read and long-read data. The method has modes for whole-genome and whole-exome sequencing and can run on tumor-normal, tumor-only and formalin-fixed paraffin-embedded samples. To train DeepSomatic and help address the dearth of publicly available training and benchmarking data for somatic variant detection, we generated and make openly available the Cancer Standards Long-read Evaluation (CASTLE) dataset of six matched tumor-normal cell line pairs whole-genome sequenced with Illumina, PacBio HiFi and Oxford Nanopore Technologies, along with benchmark variant sets. Across samples, both cell line and patient-derived, and across short-read and long-read sequencing technologies, DeepSomatic consistently outperforms existing callers.

Humans↗

PubMed2025

Contributions of Common, Rare, and Somatic Genetic Variants to Incidence of Atrial Fibrillation.

IMPORTANCE: Atrial fibrillation (AF) has a complex genetic architecture involving common, rare, and somatic variants. The association between these components requires further investigation. OBJECTIVE: To examine the individual and combined contributions of polygenic, monogenic, and somatic genetic variants to AF incidence, and develop an integrated genomic model (IGM-AF) for improved risk prediction. DESIGN, SETTING, AND PARTICIPANTS: This cohort study used whole-genome sequence data from participants of the UK Biobank, with follow-up for AF events through hospital records, death registries, and self-report. The UK Biobank recruited participants aged 40 to 69 years in the UK between 2006 and 2010. Study data were analyzed from August 2022 to November 2024. EXPOSURES: IGM-AF comprising an AF polygenic risk score (PRS), a composite rare variant gene set (AFgeneset), and somatic variants associated with clonal hematopoiesis of indeterminate potential (CHIP). Clinical AF risk was estimated using the Cohorts for Heart and Aging Research in Genomic Epidemiology AF (CHARGE-AF) score. MAIN OUTCOMES AND MEASURES: The primary outcome was hazard ratios (HRs) for 5-year incident AF attributable to PRS, AFgeneset, CHIP, and their interactions. The predictive performance of IGM-AF and its components was quantified using HRs, C statistics, and reclassification indices. RESULTS: A total of 416 085 individuals (mean [SD] age, 56.6 [8.0] years; 224 642 female [54.0%]) with 30 797 AF cases were included. The PRS (HR per 1 SD, 1.65; 95% CI, 1.63-1.67; P < 1 × 10-8), AFgeneset (HR, 1.63; 95% CI, 1.52-1.75; P = 1.46 × 10-42), and CHIP (HR, 1.26; 95% CI, 1.15-1.38; P = 1.41 × 10-6) were associated with incident AF. The 5-year cumulative incidence of AF was at least 2-fold among individuals having all 3 genetic drivers (common, rare, and somatic drivers) compared with those with only 1 driver. Integration of IGM-AF with a clinical risk model (CHARGE-AF) showed higher predictive performance (C statistic, 0.80; 95% CI, 0.80-0.80) compared with IGM-AF and CHARGE-AF alone. The classification of the at-risk population for AF was improved when IGM-AF was added to CHARGE-AF (net reclassification index, 0.08; 95% CI, 0.07-0.09). CONCLUSIONS AND RELEVANCE: Results of this cohort study demonstrated the complementary value of common, rare, and somatic variants in shaping genomic AF risk. Leveraging comprehensive genetic information may enhance screening and preventive interventions for AF.

Humans↗

PubMed2026

Distinct mutational landscapes for germline and somatic cancer variants in forty tumor suppressor genes.

Germline and somatic cancer variants in tumor suppressor genes (TSGs) share loss-of-function mechanisms, but studies of a few genes (DICER1 and CEBPA) have demonstrated differences in variant consequence and location. To systematically assess whether TSGs display distinct mutational patterns, we leveraged large public genetic databases and compared 32,941 high-quality pathogenic/likely pathogenic (P/LP) germline variants in ClinVar, with 12,907 oncogenic/likely oncogenic (O/LO) somatic tumor variants from cBioPortal across 40 TSGs. Only 3,863 (9.2%) variants were shared. Eighteen TSGs showed significantly different distributions of variant occurrences by molecular consequence, replicated with non-overlapping somatic data from the COSMIC database (chi-squared tests, false discovery rate = 5%). DICER1, TP53, and SMAD4 displayed excess somatic missense events, while nine TSGs (e.g., RB1 and APC) contained excess somatic stop-gain events throughout the coding sequence. Analysis by tumor type revealed excess stop-gain events in tissues exposed to environmental mutagens with corresponding mutation signatures. For several TSGs (WT1), germline variants predispose to tumors (Wilms' tumor) distinct from the majority source of somatic data (myeloid leukemia). Germline and somatic events are also distributed unevenly across cDNA locations, with 103 regions of preferential clustering in 39 TSGs (78 somatic and 25 germline). Twenty somatic clusters contained recurring frameshifts in homopolymer runs, many in tumors with microsatellite instability. Germline clusters contain more germline-exclusive variants, some driving non-cancer phenotypes reflecting genetic pleiotropy. Altogether, germline and somatic variants of TSGs represent unique sets with substantially different patterns shaped by selection pressures from gene-specific and somatic mutational mechanisms. Characterizing these distinctions enables more accurate clinical interpretation of TSG variants.

Humans↗

PubMed1974

Cytogenetic studies on the mechanism of formation of isoantigenic variants in somatic cell hybrids. I. Banding analyses of isoantigenic variant sublines derived from the fusion of TA3Ha carcinoma with MSWBS sarcoma cells.

Explore the source record for details and available documents.

Absorption↗

PubMed2026

Molecular residual disease assessment in colorectal and bladder cancer by somatic structural variant analysis of cell-free DNA whole-genome sequencing data.

BACKGROUND: Whole-genome sequencing (WGS)-based methods for circulating tumor DNA (ctDNA) detection typically rely on tumor-informed identification of somatic single nucleotide variants (SNVs). Somatic structural variants (SVs) are another type of cancer-specific genomic alteration, which owing to their larger genomic footprint and unique breakpoint junctions, are easier to distinguish from sequencing noise than SNVs. They are, however, rarely used for ctDNA detection because of (1) artifacts from WGS procedures that SV callers may falsely interpret as genuine SVs. This makes it difficult to establish high-confidence SV catalogos from short-read tumor WGS and can cause false-positive ctDNA detections. (2) Lack of robust strategies to quantify SV-supporting reads in plasma WGS. To address these barriers and enable integration of SV biomarkers into WGS-based ctDNA detection, we present a bioinformatic framework for algorithmic curation of somatic SV calls from fresh-frozen and formalin-fixed paraffin-embedded (FFPE) tumors, coupled with a novel approach for sensitive, accurate mapping and quantification of SV breakpoint-supporting reads in plasma WGS. METHODS: Tumor, normal and plasma WGS data from 144 patients with stage III colorectal cancer was used to establish the bioinformatic framework. This included ~30x WGS data from 1564 serially collected plasma samples. The framework was validated using tumor/normal/plasma WGS data from 32 patients with muscle-invasive bladder cancer. SV-based ctDNA detection was benchmarked against previously published SNV-based ctDNA results for the same samples. RESULTS: After curation of SV calls and quantification in plasma WGS, our SV-based approach enabled robust ctDNA detection with overall specificity exceeding 99% in plasma samples. Furthermore, we observed strong concordance (Pearson’s r > 0.93, p < 2.2 × 10− 16) between ctDNA-positive samples identified by our SV-based method and previous SNV-based analyses, validating the reliability of our approach. Finally, we demonstrated application of the method in an independent bladder cancer cohort, highlighting its generalizability and potential clinical use. CONCLUSIONS: We provide a bioinformatic framework that establishes somatic SVs as ultra-specific biomarkers for WGS-based, tumor-informed ctDNA detection. The approach delivers specific detection even when the SV catalogos are established from FFPE samples. The SV framework can stand alone or enhance SNV-based analysis pipelines.

Humans↗

PubMed1994

Histone-induced condensation of rat testis chromatin: testis-specific H1t versus somatic H1 variants.

Due the likely role of H1 histone variants in inducing the formation of folded DNA filaments with different stabilities, the condensing capacity of the testis-specific H1t versus the somatic variants was tested. Circular dichroism analyses of rat testis H1-depleted oligonucleosomes (5-2kbp) revealed that H1t, which appears in germ cells during the meiotic prophase of mammalian spermatogenesis, exerts the lowest condensing effect as compared to the other variants. The distribution of H1 subtypes among different chromatin fractions was also investigated and gave evidence that H1t is more abundant in chromatin regions which are more sensitive to DNAase I digestion.

Animals↗

PubMed1980

Identification of a trans-acting function regulation HLA-DR expression in a DR-negative B cell variant.

Somatic cell hybridizations were performed between an HLA-DR negative variant of a human B lymphoid cell line (B-LCL) and normal unrelated B-LCLs. The HLA-DR codes for polymorphic determinants on a heterodimeric cell surface lymphocyte differentiation glycoprotein. A variant subline which was selected in a single step from a diploid heterozygous DR-1 DR-3 B-LCL had lost expression of both DR-1 and DR-3 and the heretodimer; it has been described earlier. In a fusion with a DR-2 B-LCL, the hybrids expressed DR-2 and reexpressed the DR-1 and DR-3 alleles. Similar results were seen in a fusion with a different normal B-LCL. Hybrid clones from both fusions were tested with a large number of alloantisera and essentially all informative sera showed reexpression. The results show that (1) the variant did not arise by mutations in the structural genes for DR-1 and DR-3; (2) the normal cells are supplying a missing gene product needed for expression of DR; (3) this gene product is capable of acting in trans. Chromosome counts showed that the apparent recessiveness of the variant in the hybrids was not due to chromosomal segregation.

B-Lymphocytes↗

PubMed1988

Properties of anti-idiotypic T cell lines propagated with syngeneic B lymphocytes. I. T cells bind intact idiotypes and discriminate between the somatic idiotypic variants in a manner similar to the anti-idiotopic antibodies.

We describe the development of T cell lines possessing a binding specificity for syngeneic T15 idiotopes (Id) expressed on phosphorylcholine (PC)-reactive Ig. The lines were obtained by cultivation of BALB/c splenic T cells with T15 Id+ stimulator cells BCg3R-1d, a BALB/c lymphoma transfected with genomic sequences mu and kappa with S107 (T15) variable regions. Resulting Thyl-2+, L3T4+ cell lines depend on the T15 Id+ BCg3R-1d cells for growth and demonstrate the ability to bind TEPC15, a S107 germline-encoded, PC-specific Ig alpha. The specificity of the 125I-TEPC15 binding was studied in a competitive RIA with various unlabeled Ig. The isolated H and L chains of TEPC15 failed to inhibit the 125I-TEPC15 binding, and the T15-, PC-binding proteins M603 (alpha) and M511 (alpha) inhibited the binding either poorly or not at all. Moreover, the T cell lines had a discriminatory binding specificity for various T15+ Ig that are somatic variants of TEPC15 and that differ from each other in discrete, conformational Id (epitopes) detectable with specific monoclonal anti-Id. The T cell lines could be grouped according to their binding patterns, which were comparable to the recognition patterns of certain monoclonal anti-Id. These data suggest the existence of T cells with specificity for serologically-defined determinants of syngeneic idiotypes.

Animals↗

PubMed2025

Comprehensive benchmarking of somatic structural variant detection at ultra-low allele fractions.

Postzygotic mosaicism gives rise to somatic structural variants (SVs) at ultra-low variant allele fractions (VAFs), which pose challenges for detection due to the high-coverage sequencing required and noise introduced by sequencing artifacts. Although somatic SV detection has been extensively studied in cancer, these studies are not directly applicable to the study of tissue mosaicism, as they rely on matched normals, target higher VAF ranges, and are enriched for different types of SVs. We present comprehensive benchmark data and best practices for non-cancer somatic SV detection. We created a synthetic mosaic sample by combining six HapMap individuals at varying proportions, generating allele fractions as low as 0.25%. This sample was sequenced to ~2,300x total coverage using Illumina, PacBio, and Nanopore technologies across multiple sequencing centers. A high-confidence benchmark SV set containing over 21,000 pseudo-somatic insertions and deletions ≥50bp was derived from haplotype-resolved assemblies. We evaluated 12 SV discovery pipelines and identified caller-specific strengths and sequencing platform-specific shortcomings. We find that short read-based approaches show reduced recall for insertions and repeat-associated SVs, whereas long-read sequencing achieves high accuracy throughout the genome, increasing linearly with coverage. The best algorithm's sensitivity exceeded 80% for VAFs ≥4% and 15% for VAFs of 0.5-1% with 60x coverage. The publicly available benchmarking data and comparative analysis of current methods provide a foundation for robust discovery of SV mosaicism in non-cancer tissues..

Journal Article↗

PubMed2026

Somatic-only SDHD variant with tumor-specific loss of heterozygosity in metastatic carotid body tumor: a case report with review of literature.

Carotid body tumors (CBTs) are rare paragangliomas in which genetic predisposition, particularly pathogenic variants in succinate dehydrogenase (SDHx) genes, plays an important role in tumorigenesis. Previous genetic studies of CBTs have primarily focused on germline SDHx variants, whereas somatic alterations remain poorly characterized. Among SDHx genes, SDHB variants are known to be associated with a higher metastatic risk, while SDHD variants are generally linked to a lower metastatic rate. We report the case of a 41-year-old man who presented with a painless right-sided neck mass. Imaging studies demonstrated a hypervascular tumor located at the carotid bifurcation with circumferential encasement of the carotid artery. During surgery, the tumor was classified as a Shamblin type III CBT, and complete surgical resection with vascular reconstruction was performed. Histopathological examination confirmed paraganglioma with metastasis to a single cervical lymph node. Germline genetic testing did not reveal any pathogenic variants. However, comprehensive tumor genomic profiling identified a somatic SDHD c.304C>G (p.His102Asp) variant accompanied by tumor-specific loss of heterozygosity (LOH) and copy-number loss at the SDHD locus, findings compatible with biallelic SDHD inactivation. The patient remained free of recurrence during follow-up. To our knowledge, this represents the first reported case of metastatic CBT harboring a somatic-only SDHD variant with tumor-specific LOH. This case suggests that reliance on germline testing alone may underestimate the molecular drivers of CBT and highlights the potential clinical value of tumor-based genomic profiling for risk stratification and prognostic assessment.

Carotid body tumor↗

PubMed2026

Deep tissue sequencing improves genetic diagnostic yield in focal cortical dysplasia.

Focal cortical dysplasias (FCDs) are malformations of cortical development associated with drug-resistant focal epilepsy. We analyzed surgical tissue from 25 consecutive cases recruited from adult and pediatric epilepsy surgery programs. We performed high-depth sequencing of lesional tissue, validated somatic variants using droplet digital PCR or amplicon sequencing, and investigated genotype-phenotype correlations. A pathogenic or likely pathogenic variant was detected in 64% (n = 16/25) of cases. Of these, five cases with FCDIIa or FCDIIb had germline variants in NPRL3 (n = 3) or DEPDC5 (n = 2). Somatic variants were identified in 44% (n = 11/25) of cases. The genetic yield for FCDIIb was 77% of cases having a pathogenic mTOR pathway variant detected (n = 10/13), and for FCDIIa 66% (n = 6/9). High depth sequencing approaches allowed detection of somatic variants with very low (down to 0.4%) variant allele fractions (VAFs). No pathogenic variants were detected in 3 cases with FCDI. 62% (n = 15/24) of the cases with ≥12 months follow up experienced a favourable seizure outcome (Engel 1-2) following surgery. Of note, n = 9 patients required repeat surgery to resect residual dysplasia. Determining a genetic diagnosis reveals aetiology and paves the way to precision therapies that may benefit those with FCD who do not respond to current treatments.

Humans↗

PubMed2026

Improving long-read somatic structural variant calling with pangenome and de novo personal genome assembly.

Accurate detection of mosaic and somatic structural variants (SVs) provides early diagnostic and therapeutic evidence for cancers. While long-read whole-genome sequencing leads to more accurate SV detection than short read sequencing, existing long-read SV callers only look at alignment against a single reference genome and are susceptible to systematic false discovery caused by germline differences between the individual genome and the reference genome. Here we develop a new SV filtering method that jointly considers the alignment against a pangenome and the de novo assembly of the germline genome. It dramatically reduces false positive mosaic and somatic SVs in cancer cell lines with little loss in sensitivity for existing long read SV callers. Our study highlights the essential need for pangenome or personal genome assembly to integrate SV calls for both SV discoveries and clinical diagnostics.

Journal Article↗

PubMed2026

Somatic mosaicism in the brain: linking development, ageing and neurodegeneration.

Somatic mosaicism is increasingly recognized as a pervasive feature of the human brain and a potential contributor to neurological disease across the lifespan. Unlike germline variants, somatic variants arise post-zygotically and are unevenly distributed across regions, cell types and even individual neurons, enabling focal biological effects that can scale to network-level dysfunction. In this Review, we synthesize current evidence that developmental timing, clonal architecture and cell-type-specific selective pressures shape how somatic variants influence brain structure and function. Early embryonic variants can produce broad regional clones and severe phenotypes, whereas later events are usually more restricted; with ageing, ongoing DNA damage and imperfect repair generate private variants that might cumulatively reduce cellular resilience. We also summarize advances in detection approaches, including bulk, error-corrected and single-cell sequencing, and discuss their strengths and current limitations for clinical translation. Emerging data link brain somatic variants to neurodevelopmental and neurodegenerative phenotypes, supporting a unified framework in which mosaic genetics bridges focal lesions and distributed neurological syndromes. Integrating genomic, cellular and physiological analyses in longitudinal human studies will be essential to define causality, identify biomarkers and guide future targeted interventions.

Journal Article↗

PubMed2001

Mutation analysis of NTRK2 and NTRK3, encoding 2 tyrosine kinase receptors, in sporadic human medullary thyroid carcinoma reveals novel sequence variants.

Somatic mutations in the proto-oncogene RET are found in 25% to 80% of sporadic medullary thyroid carcinomas (MTCs). The significance of somatic RET mutation in MTC initiation and progression, however, remains unknown. Like RET, TRK is a neurotrophic receptor tyrosine kinase. Immunostaining has shown that only a subset of normal C cells expresses Trk family receptors, but in C-cell hyperplasia, they consistently express NTRK2, with variable expression of NTRK1 and NTRK3. In later stages of MTC, NTRK2 expression was reduced while NTRK3 expression was increased. In the context of these data, we sought to determine whether sequence variants in NTRK2 and NTRK3 are responsible for these differences in protein expression. We determined the genomic structure of NTRK2 and found that it consists of at least 17 exons varying in size from 36 to 306 bp. Mutation analysis of sporadic MTC did not reveal any sequence variants in NTRK2 but did reveal 3 variants in NTRK3, c.573C >T (N191N, exon 5), c.678T > C (N226N, exon 6) and c.1488C > G (A496A, exon 12) occurring among 19 chromosomes (31%), 1 chromosome (2%) and 24 chromosomes (39%), respectively. Corresponding germline also harbored these variants. There was a trend toward excess association of the NTRK3 variant c.1488C > G (A496A) in cases (24/62 chromosomes, 39%) compared to controls (18/62, 29%), but this difference did not reach significance (p > 0.05). The remaining 2 NTRK3 variants occurred with similar frequencies between MTC cases and population-matched controls (19 vs. 17 and 1 vs. 0, p > 0.05). We conclude that sequence variants in NTRK2 and NTRK3 are not likely to be responsible for large differences in expression at the protein level, but we cannot exclude very low penetrance effects.

Carcinoma, Medullary↗

PubMed1992

A human systemic lupus erythematosus-related anti-cardiolipin/single-stranded DNA autoantibody is encoded by a somatically mutated variant of the developmentally restricted 51P1 VH gene.

We report the Ig H and L chain V region sequences from the cDNAs encoding a monoclonal human IgG anti-cardiolipin/ssDNA autoantibody (R149) derived from a patient with active SLE. Comparison with the germ-line V-gene repertoire of this patient revealed that R149 likely arose as a consequence of an Ag-driven selection process. The Ag-binding portions of the V regions were characterized by a high number of arginine residues, a property that has been associated with anti-dsDNA autoantibodies from lupus-prone mice and patients with SLE. The VH gene encoding autoantibody R149 was a somatically mutated variant of the 51P1 gene segment, which is frequently associated with the restricted fetal B cell repertoire, malignant CD5 B cells, and natural autoantibodies. These data suggest that in SLE patients a common antigenic stimulus may evoke anti-DNA and anti-cardiolipin autoantibodies and provide further evidence that a small set of developmentally restricted VH genes can give rise to disease-associated autoantibodies through Ag-selected somatic mutations.

Amino Acid Sequence↗

PubMed1996

Co-operative interactions of oligonucleosomal DNA with the H1e histone variant and its poly(ADP-ribosyl)ated isoform.

H1 histone somatic variants from L929 mouse fibroblasts were purified by reverse-phase HPLC. We analysed the ability of each H1 histone variant to allow the H1-H1 interactions that are essential for the formation of the higher levels of chromatin structure, and we investigated the role played by the poly(ADP-ribosyl)ation process. Cross-linking analysis showed that H1e is the only somatic variant which, when bound to DNA, is able to produce H1-H1 polymers; the size of polymers was decreased when H1e was enriched in its poly(ADP-ribosyl)ated isoform. Measurement of the methyl-accepting ability in native nuclei compared with nuclei in which poly(ADP-ribosyl)ation was induced showed that the poly(ADP-ribosyl)ated H1 histone had not been removed from linker regions, in spite of its different interaction with DNA.

Animals↗

PubMed1980

Linked marker analysis of spontaneous HLA variants of somatic cells.

Fourteen independent, spontaneously arising HLA-A1 and -B27 variant clones were isolated from the pseudodiploid B lymphoid line T5-1 by selection using A1 and B27 alloantiserum, respectively, and complement. T5-1 is heterozygous for the tightly linked loci HLA-DR, -B, and -A and probably for -C as well. Following recloning, each of the variants was tested for the HLA specificities of T5-1. None of these spontaneous variants had a genetic lesion which extended to the nearest flanking HLA marker, less than 1 cM distant. On the other hand, a variant isolated from mutagenized cells had a lesion which extended completely through one HLA region in cis, suggesting that haploid expression in the HLA region is compatible with viability, that there are no recessive lethals in the opposite HLA region, and that spontaneously arising variants with extended lesions could have been recovered had they occurred. From these results, we conclude that extended genetic lesions of 0.8 cM or longer contribute less than 10% to the overall rate of spontaneous variation for HLA alleles, which we previously estimated at about 5 x 10(-7) per cell per generation.

B-Lymphocytes↗

PubMed2026

A comprehensive survey of genetic variants in neuroblastoma.

BACKGROUND: Neuroblastoma (NB) is the most common extracranial solid tumor in children and is characterized by marked clinical and molecular heterogeneity. Genomic alterations play a critical role in NB pathogenesis; however, population-specific mutational features remain insufficiently characterized, particularly among Chinese patients. METHODS: Whole-exome sequencing (WES) was performed on tumor, para-tumor, and matched peripheral blood samples from nine pathologically confirmed Chinese patients with NB. Somatic variant profiles were compared with four publicly available NB datasets from cBioPortal, published in 2012, 2013, 2015, and 2023. Mutational patterns, recurrently altered genes, and Gene Ontology (GO) enrichment were analyzed using R version 4.3.2 and clusterProfiler version 4.10.0. RESULTS: A total of 77 missense variants were identified in our cohort. Single-nucleotide polymorphisms (SNPs) represented the predominant variant type, and C > T substitutions were the most frequent nucleotide change. MAP1A variants, comprising two missense variants in one patient, and RBM33 variants, comprising two distinct variants in two patients, were detected in our cohort and, to the best of our knowledge, have not been previously reported in NB, although their frequencies were low. No MYCN amplification or variants in ALK, ATRX, or DAXX were detected. Comparative analysis with the cBioPortal datasets revealed no somatic variants universally shared across all cohorts. In addition, high-risk patients exhibited distinct mutational patterns, with enrichment of the Gene Ontology term "collagen-containing extracellular matrix." CONCLUSIONS: These findings highlight the molecular diversity of NB and suggest the presence of potential population-specific genetic features in Chinese patients. The low-frequency MAP1A and RBM33 variants identified in this cohort warrant further validation in larger, independent cohorts. Moreover, the enrichment of extracellular matrix-related pathways in high-risk NB supports further investigation of tumor-microenvironment interactions as potential therapeutic targets.

Extracellular matrix↗