Whole-exome sequencing-based mutational profiling of hepatocellular adenoma malignant transformation to hepatocellular carcinoma.

Hepatocellular adenoma (HCA) represents a rare benign hepatic neoplasm with potential for malignant transformation into hepatocellular carcinoma (HCC), yet the underlying mechanism remains elusive. In this study, we investigated the genomic landscape of this process to identify therapeutic strategies for blocking malignant transformation. Using micro-detection techniques, we obtained specimens of adenoma, cancerous neoplasm and adjacent normal liver from three patients undergoing hepatic resection surgery. Whole-exome sequencing (WES) was performed, and genomic interactions between HCA and HCC components within the same tumour were evaluated using somatic variant calling, copy number variation (CNV) analysis, clonality evaluation and mutational signature analysis. Our results revealed genomic heterogeneity among patient cases, yet within each sample, HCA and HCC tissues exhibited a similar mutational landscape, suggesting a high degree of homology. Using nonnegative matrix factorization and phylogenetic trees, we identified shared and distinct mutational characteristics and uncovering necessary pathways associated with HCA-HCC malignant transformation. Remarkably, we found that HCA and HCC shared a common monoclonal origin while displaying significant genetic diversity within HCA-HCC tumours, indicating fundamental genetic connections or evolutionary pathways between the two. Moreover, elevated immune therapy-related markers in these patients suggested heightened sensitivity to immune therapy, providing novel avenues for the treatment of hepatic malignancies. This study sheds light on the genetic mechanisms underlying HCA-HCC progression, offering potential targets for therapeutic intervention and highlighting the promise of immune-based therapies in managing hepatic malignancies.

Patients with a new-onset cutaneous sebaceous neoplasm following immunosuppression should be evaluated for Muir-Torre syndrome with germline mismatch repair gene mutation analysis: case reports.

Patients with Muir-Torre syndrome may have a systemic malignancy and a sebaceous neoplasm such as an adenoma, epithelioma, and/or carcinoma. The syndrome usually results from a germline mutation in one or more mismatch repair genes. Iatrogenic or acquired immunosuppression can promote the appearance of sebaceous tumors, either as an isolated event or as a feature of Muir-Torre syndrome and may unmask individuals genetically predisposed to the syndrome. Two iatrogenically immunosuppressed men with Muir-Torre syndrome features are described. Similar to these immunocompromised men, Muir-Torre syndrome-associated sebaceous neoplasms have occurred in solid organ transplant recipients, human immunodeficiency virus-infected individuals, and patients with chronic diseases who are treated with immunosuppressive agents. Muir-Torre syndrome-associated sebaceous neoplasms occur more frequently and earlier in kidney recipients, who are receiving more post-transplant immunosuppressive agents, than in liver recipients. The development of sebaceous neoplasms is decreased by replacing cyclosporine or tacrolimus with sirolimus or everolimus. Specific anti-cancer vaccines or checkpoint blockade immunotherapy may merit exploration for immune-interception of Muir-Torre syndrome-associated sebaceous neoplasms and syndrome-related visceral cancers. We suggest germline testing for genomic aberrations of mismatch repair genes should routinely be performed in all patients-both immunocompetent and immunosuppressed-who develop a Muir-Torre syndrome-associated sebaceous neoplasm.

Combined DNA Analysis from Stool and Blood Samples Improves Tumor Tracking and Assessment of Clonal Heterogeneity in Localized Rectal Cancer Patients.

Objectives: In this study, stool samples were evaluated for tumor mutation analysis via a targeted next generation sequencing (NGS) approach in a small patient cohort suffering from localized rectal cancer. Introduction: Colorectal cancer (CRC) causes the second highest cancer-related death rate worldwide. Thus, improvements in disease assessment and monitoring that may facilitate treatment allocation and allow organ-sparing "watch-and-wait" treatment strategies are highly relevant for a significant number of CRC patients. Methods: Stool-based results were compared with mutation profiles derived from liquid biopsies and the gold standard procedure of tumor biopsy from the same patients. A workflow was established that enables the detection of de-novo tumor mutations in stool samples of CRC patients via ultra-sensitive cell-free tumor DNA target enrichment. Results: Notably, only a 19% overall concordance was found in mutational profiles across the compared sample specimens of stool, tumor, and liquid biopsies. Conclusion: Based on these results, the analysis of stool and liquid biopsy samples can provide important additional information on tumor heterogeneity and potentially on the assessment of minimal residual disease and clonal tumor evolution.

Diagnostic value of FNAC combined with BRAFV600E mutation detection in Hashimoto's thyroiditis complicated with papillary thyroid carcinoma.

Background: This study aimed to analyze the effect of preoperative fine needle aspiration cytology (FNAC) combined with BRAFV600E mutation detection as compared to that of fine needle aspiration cytology alone on the diagnostic performance of papillary thyroid carcinoma (PTC) combined with Hashimoto's thyroiditis (HT). Method: Patients with thyroid nodules in Hashimoto's thyroiditis, who underwent fine-needle aspiration cytology examination and BRAFV600E mutation detection in the puncture eluate at the outpatient clinic, were selected. Finally, 122 patients received surgical treatment and were included in the study. We used postoperative pathological results as the gold standard. Accordingly, we compared the sensitivity, specificity and accuracy of preoperative FNAC alone and FNAC combined with BRAFV600E mutation detection in for the diagnosis of PTC combined with HT. Results: For PTC patients with HT, the sensitivity of FNAC diagnosis was 93.69%, the specificity was 90.90% and the accuracy was 93.44%. However, the sensitivity, specificity and accuracy of FNAC combined with BRAFV600E mutation detection were 97.30%, 90.90% and 96.72%, respectively. Therefore, combined detection can improve the sensitivity and accuracy of diagnosis (p<0.05). Conclusion: FNAC combined with eluent BRAFV600E mutation detection can improve the sensitivity and accuracy of diagnosis of PTC in the background of HT.

PIK3CA Hotspot Mutations in Saliva as a Diagnostic Marker in Oral Squamous Cell Carcinoma Patients.

BACKGROUND/AIM: This study aimed at the analogous detection of PIK3CA mutations, common in oral squamous cell carcinoma (OSCC), in matched tumor and saliva samples. PATIENTS AND METHODS: Tissue and saliva samples were obtained from 29 patients diagnosed with primary OSCC. Saliva samples were obtained preoperatively; tissue specimens were acquired during tumor resection. Tumor DNA was extracted from both tissue and saliva samples. All samples were controlled for DNA quantity and quality and genetic matching of sample pairs was confirmed using the iPlex Pro Exome QC Panel. Variant detection was performed using the MassARRAY® System, a mass-spectrometry based detection system. Mutational analysis in tissue tumor DNA was made using the multiplexed ClearSEEK™ PIK3CA v1.0 Panel covering 20 hotspot mutations in PIK3CA. In saliva samples, variants were analyzed using both the ClearSEEK™ and the UltraSEEK® Lung v1.1 Panel, with a higher limit of detection but covering less PIK3CA variants. RESULTS: Overall, a PIK3CA variant was found in seven of the 29 tumor tissue samples (24%) by ClearSEEK™; UltraSEEK® additionally confirmed the variant in four of these seven positive samples. Of the three variants not detected by UltraSEEK®, two were not included in the panel and one was included but not detected. Of the seven variants found in tissue, five could also be detected in the matching saliva samples (71%), either by utilizing ClearSEEK™ or UltraSEEK® Conclusion: The detection of PIK3CA hotspot mutations in OSCC and their simultaneous occurrence in saliva underline the potential benefit of liquid biopsies for non-invasive cancer detection and follow-up care of OSCC patients.

Mutational analysis of pulmonary large cell neuroendocrine carcinoma: APC gene mutations identify a good prognostic factor.

Pulmonary large cell neuroendocrine carcinoma (LCNEC) is a highly aggressive neoplasm with biological heterogeneity. Mutations in multiple genes have been identified in LCNEC. However, associations between gene alterations, histopathological characteristics, and prognosis remain ambiguous. Here, we investigated the clinicopathologic, immunohistochemical, and genomic characteristics of 19 patients with LCNEC and 9 patients with atypical carcinoid (AC). We revealed high mutation frequencies of TP53 (89.5 %), RB1 (42.1 %), APC (31.6 %), and MCL1 (31.6 %) in LCNEC, while genetic alterations were rarely found in AC. APC alterations mainly occurred to the exon 16 and were only identified in LCNEC with wild-type RB1. The 19 LCNEC were further subgrouped into APC wild-type (LCNEC-APCMT, 6/19) and APC-mutated (LCNEC-APCWT, 13/19) subgroups. In comparison with LCNEC-APCWT, LCNEC-APCMT displayed lower TMB (median: 12.64 vs 4.20, P = 0.045), and relatively mild cytologic atypia. In addition, LCNEC-APCMT distinguished itself from AC and LCNEC-APCWT by obviously downregulated expression of neuroendocrine markers (CD56 and Syn, P < 0.01) and significantly altered expression of genes downstream of APC (ß-catenin migrating into the cytoplasm and nucleus, P < 0.001; c-Myc upregulating, P = 0.005). The OS of LCNEC-APCMT was numerically intermediate between AC and LCNEC-APCWT. We first proposed that APC alterations were common in LCNEC with wild-type RB1 and that LCNEC-APCMT was associated with lower TMB and better OS in comparison with LCNEC-APCWT.

CRISPR/Cas13a-based supersensitive circulating tumor DNA assay for detecting EGFR mutations in plasma.

Despite recent technological advancements in cell tumor DNA (ctDNA) mutation detection, challenges persist in identifying low-frequency mutations due to inadequate sensitivity and coverage of current procedures. Herein, we introduce a super-sensitivity and specificity technique for detecting ctDNA mutations, named HiCASE. The method utilizes PCR-based CRISPR, coupled with the restriction enzyme. In this work, HiCASE focuses on testing a series of EGFR mutations to provide enhanced detection technology for non-small cell lung cancer (NSCLC), enabling a detection sensitivity of 0.01% with 40 ng cell free DNA standard. When applied to a panel of 140 plasma samples from 120 NSCLC patients, HiCASE exhibits 88.1% clinical sensitivity and 100% specificity with 40 µL of plasma, higher than ddPCR and Super-ARMS assay. In addition, HiCASE can also clearly distinguish T790M/C797S mutations in different positions at a 1% variant allele frequency, offering valuable guidance for drug utilization. Indeed, the established HiCASE assay shows potential for clinical applications.

Large-scale mutational analysis identifies UNC93B1 variants that drive TLR-mediated autoimmunity in mice and humans.

Nucleic acid-sensing Toll-like receptors (TLR) 3, 7/8, and 9 are key innate immune sensors whose activities must be tightly regulated to prevent systemic autoimmune or autoinflammatory disease or virus-associated immunopathology. Here, we report a systematic scanning-alanine mutagenesis screen of all cytosolic and luminal residues of the TLR chaperone protein UNC93B1, which identified both negative and positive regulatory regions affecting TLR3, TLR7, and TLR9 responses. We subsequently identified two families harboring heterozygous coding mutations in UNC93B1, UNC93B1+/T93I and UNC93B1+/R336C, both in key negative regulatory regions identified in our screen. These patients presented with cutaneous tumid lupus and juvenile idiopathic arthritis plus neuroinflammatory disease, respectively. Disruption of UNC93B1-mediated regulation by these mutations led to enhanced TLR7/8 responses, and both variants resulted in systemic autoimmune or inflammatory disease when introduced into mice via genome editing. Altogether, our results implicate the UNC93B1-TLR7/8 axis in human monogenic autoimmune diseases and provide a functional resource to assess the impact of yet-to-be-reported UNC93B1 mutations.

Retinal Dystrophies Associated With Peripherin-2: Genetic Spectrum and Novel Clinical Observations in 241 Patients.

Purpose: To describe the clinical, electrophysiological and genetic spectrum of inherited retinal diseases associated with variants in the PRPH2 gene. Methods: A total of 241 patients from 168 families across 15 sites in 9 countries with pathogenic or likely pathogenic variants in PRPH2 were included. Records were reviewed for age at symptom onset, visual acuity, full-field ERG, fundus colour photography, fundus autofluorescence (FAF), and SD-OCT. Images were graded into six phenotypes. Statistical analyses were performed to determine genotype-phenotype correlations. Results: The median age at symptom onset was 40 years (range, 4-78 years). FAF phenotypes included normal (5%), butterfly pattern dystrophy, or vitelliform macular dystrophy (11%), central areolar choroidal dystrophy (28%), pseudo-Stargardt pattern dystrophy (41%), and retinitis pigmentosa (25%). Symptom onset was earlier in retinitis pigmentosa as compared with pseudo-Stargardt pattern dystrophy (34 vs 44 years; P = 0.004). The median visual acuity was 0.18 logMAR (interquartile range, 0-0.54 logMAR) and 0.18 logMAR (interquartile range 0-0.42 logMAR) in the right and left eyes, respectively. ERG showed a significantly reduced amplitude across all components (P < 0.001) and a peak time delay in the light-adapted 30-Hz flicker and single-flash b-wave (P < 0.001). Twenty-two variants were novel. The central areolar choroidal dystrophy phenotype was associated with 13 missense variants. The remaining variants showed marked phenotypic variability. Conclusions: We described six distinct FAF phenotypes associated with variants in the PRPH2 gene. One FAF phenotype may have multiple ERG phenotypes, demonstrating a discordance between structure and function. Given the vast spectrum of PRPH2 disease our findings are useful for future clinical trials.

Rosace: a robust deep mutational scanning analysis framework employing position and mean-variance shrinkage.

Deep mutational scanning (DMS) measures the effects of thousands of genetic variants in a protein simultaneously. The small sample size renders classical statistical methods ineffective. For example, p-values cannot be correctly calibrated when treating variants independently. We propose Rosace, a Bayesian framework for analyzing growth-based DMS data. Rosace leverages amino acid position information to increase power and control the false discovery rate by sharing information across parameters via shrinkage. We also developed Rosette for simulating the distributional properties of DMS. We show that Rosace is robust to the violation of model assumptions and is more powerful than existing tools.

Resolving Discrepancies in Idylla BRAF Mutational Assay Results Using Targeted Next-Generation Sequencing.

BRAF mutation identification is important for the diagnosis and treatment of several tumor types, both solid and hematologic. Rapid identification of BRAF mutations is required to determine eligibility for targeted BRAF inhibitor therapy. The Idylla BRAF mutation assay is a rapid, multiplex allele-specific PCR test designed to detect the most common oncogenic BRAF V600 mutations in formalin-fixed paraffin-embedded (FFPE) tissue samples. Here, we describe the validation of the Idylla BRAF mutation assay in our laboratory. During routine clinical practice, we noticed cases in which BRAF V600 mutations were identified with unusual amplification curves, with three cases displaying a delayed amplification within a double amplification pattern and two false-positive calls. We therefore initiated a quality improvement effort to systematically and retrospectively evaluate next-generation sequencing (NGS)-tested cases with BRAF mutations identified within five amino acids of BRAF codon V600 and did not identify additional false-positive cases. We hypothesize that late amplification in a double amplification pattern may represent non-specific amplification, whereas cases displaying single delayed amplification curves may stem from the presence of either non-V600 variants, very low-level V600 variants, cytosine deamination artifacts, and/or non-specific amplification by an allele-specific PCR primer. Regardless, we recommend that Idylla BRAF cases with non-classical amplification curves undergo reflex NGS testing. These findings are likely relevant for other Idylla assays interrogating hotspot mutations in genes such as EGFR, IDH1/2, KRAS, and NRAS.

Liquid Biopsies as Non-Invasive Tools for Mutation Profiling in Multiple Myeloma: Application Potential, Challenges, and Opportunities.

Over the last decades, the survival of multiple myeloma (MM) patients has considerably improved. However, despite the availability of new treatments, most patients still relapse and become therapy-resistant at some point in the disease evolution. The mutation profile has an impact on MM patients' outcome, while typically evolving over time. Because of the patchy bone marrow (BM) infiltration pattern, the analysis of a single bone marrow sample can lead to an underestimation of the known genetic heterogeneity in MM. As a result, interest is shifting towards blood-derived liquid biopsies, which allow for a more comprehensive and non-invasive genetic interrogation without the discomfort of repeated BM aspirations. In this review, we compare the application potential for mutation profiling in MM of circulating-tumor-cell-derived DNA, cell-free DNA and extracellular-vesicle-derived DNA, while also addressing the challenges associated with their use.

Trio-based whole-exome sequencing reveals mutations in early-onset high myopia.

PURPOSE: Myopia, especially high myopia (HM), represents a widespread visual impairment with a globally escalating prevalence. This study aimed to elucidate the genetic foundations associated with early-onset HM (eoHM) while delineating the genetic landscape specific to Shaanxi province, China. METHODS: A comprehensive analysis of whole-exome sequencing was conducted involving 26 familial trios displaying eoHM. An exacting filtration protocol identified potential candidate mutations within acknowledged myopia-related genes and susceptibility loci. Subsequently, computational methodologies were employed for functional annotations and pathogenicity assessments. RESULTS: Our investigation identified 7 genes and 10 variants associated with HM across 7 families, including a novel mutation in the ARR3 gene (c.139C>T, p.Arg47*) and two mutations in the P3H2 gene (c.1865T>C, p.Phe622Ser and c.212T>C, p.Leu71Pro). Pathogenic mutations were found in syndromic myopia genes, notably encompassing VPS13B, TRPM1, RPGR, NYX and RP2. Additionally, a thorough comparison of previously reported causative genes of syndromic myopia and myopia risk genes with the negative sequencing results pinpointed various types of mutations within risk genes. CONCLUSIONS: This investigation into eoHM within Shaanxi province adds to the current understanding of myopic genetic factors. Our results warrant further functional validation and ocular examinations, yet they provide foundational insights for future genetic research and therapeutic innovations in HM.

Microfluidics-based EGFR mutation detection and its implication in the resource-limited clinical setting.

Management of lung cancer today obligates a mutational analysis of the epidermal growth factor receptor (EGFR) gene particularly when Tyrosine Kinase Inhibitor (TKI) therapy is being considered as part of prognostic stratification. This study evaluates the performance of automated microfluidics-based EGFR mutation detection and its significance in clinical diagnostic settings. Formalin-fixed, paraffin-embedded (FFPE) samples from NSCLC patients (n = 174) were included in a two-phase study. Phase I: Validation of the platform by comparing the results with conventional real-time PCR and next-generation sequencing (NGS) platform. Phase II: EGFR mutation detection on microfluidics-based platform as part of routine diagnostics workup. The microfluidics-based platform demonstrates 96.5% and 89.2% concordance with conventional real-time PCR and NGS, respectively. The system efficiently detects mutations across the EGFR gene with 88.23% sensitivity and 100% specificity. Out of 144 samples analysed in phase II, the platform generated valid results in 94% with mutation detected in 41% of samples. This microfluidics-based platform can detect as low as 5% mutant allele fractions from the FFPE samples. Therefore the microfluidics-based platform is a rapid, complete walkaway, with minimum tissue requirement (two sections of 5 µ thickness) and technical skill requirement. The method can detect clinically actionable EGFR mutations efficiently and can be considered a reliable diagnostic platform in resource-limited settings. From receiving samples to reporting the results this platform provides accurate data without much manual intervention. The study helped to devise an algorithm that emphasizes effective screening of the NSCLC cases for EGFR mutations with varying tumour content. Thus it helps in triaging the cases judiciously before proceeding with multigene testing.

A Novel COL2A1 Gene Pathogenic Variant in a Turkish Family With Ocular Stickler Syndrome.

A 6-month-old female infant with megalophthalmos was referred with the suspicion of congenital glaucoma. Refractive measurements obtained with handheld autorefractometry were -7.00 -2.00 × 90° in the right eye and -6.00 -2.00 × 100° in the left eye and ultrasonic axial lengths were 22.50 mm in both eyes. Intraocular pressures and vertical and horizontal corneal diameters of the proband were 11 mm Hg, 11 mm, and 11.50 mm in both eyes, respectively. She was diagnosed as having early-onset high myopia. Her father also had degenerative high myopia (-12.00 diopters) in the right eye, bilateral congenital lens opacities, and retinal detachment in the left eye. Her mother was emmetropic with normal eye examination results. Clinical exome sequencing analysis revealed a novel ENST00000380518.3 c.3528_3530 delins GACCATTAGCA (Chr12:48369813: GCA > TGCTAATGGTC) variant in the collagen type II alpha 1 chain (COL2A1) on chromosome 12q13 (OMIM 108300), consistent with the Stickler syndrome type 1. Subsequent segregation analysis revealed paternal inheritance. Although many pathogenic null variants have been described within the COL2A1 gene, there is currently no documented literature pertaining to this specific variant, making this the inaugural report of its manifestation in scientific discourse. [J Pediatr Ophthalmol Strabismus. 2024;61(3):e23-e27.].

Genotypic and Phenotypic Characterization of a Cohort of Patients Affected by Rod Cyclic Nucleotide Channel-Associated Retinitis Pigmentosa.

INTRODUCTION: Retinitis pigmentosa (RP), a heterogeneous inherited retinal disorder causing gradual vision loss, affects over 1 million people worldwide. Pathogenic variants in CNGA1 and CNGB1 genes, respectively, accounting for 1% and 4% of cases, impact the cyclic nucleotide-gated channel in rod photoreceptor cells. The aim of this study was to describe and compare genotypic and clinical characteristics of a cohort of patients with CNGA1- or CNGB1-related RP and to explore potential genotype-phenotype correlations. METHODS: The following data from patients with CNGA1- or CNGB1-related RP, followed in five Italian inherited retinal degenerations services, were retrospectively collected: genetic variants in CNGA1 and CNGB1, best-corrected visual acuity (BCVA), ellipsoid zone (EZ) width, fundus photographs, and short-wavelength fundus autofluorescence (SW-AF) images. Comparisons and correlation analyses were performed by first dividing the cohort in two groups according to the gene responsible for the disease (CNGA1 and CNGB1 groups). In parallel, the whole cohort of RP patients was divided into two other groups, according to the expected impact of the variants at protein level (low and high group). RESULTS: In total, 29 patients were recruited, 11 with CNGA1- and 18 with CNGB1-related RP. In both CNGA1 and CNGB1, 5 novel variants in CNGA1 and 5 in CNGB1 were found. BCVA was comparable between CNGA1 and CNGB1 groups, as well as between low and high groups. CNGA1 group had a larger mean EZ width compared to CNGB1 group, albeit not statistically significant, while EZ width did not differ between low and high groups A statistically significant correlation between EZ width and BCVA as well as between EZ width and age were observed in the whole cohort of RP patients. Fundus photographs of all patients in the cohort showed classic RP pattern, and in SW-AF images an hyperautofluorescent ring was observed in 14/21 patients. CONCLUSION: Rod CNG channel-associated RP was demonstrated to be a slowly progressive disease in both CNGA1- and CNGB1-related forms, making it an ideal candidate for gene augmentation therapies.

Mutational Analysis of RIP Type I Dianthin-30 Suggests a Role for Arg24 in Endocytosis.

Saponin-mediated endosomal escape is a mechanism that increases the cytotoxicity of type I ribosome-inactivating proteins (type I RIPs). In order to actualize their cytotoxicity, type I RIPs must be released into the cytosol after endocytosis. Without release from the endosomes, type I RIPs are largely degraded and cannot exert their cytotoxic effects. Certain triterpene saponins are able to induce the endosomal escape of these type I RIPs, thus increasing their cytotoxicity. However, the molecular mechanism underlying the endosomal escape enhancement of type I RIPs by triterpene saponins has not been fully elucidated. In this report, we investigate the involvement of the basic amino acid residues of dianthin-30, a type I RIP isolated from the plant Dianthus caryophyllus L., in endosomal escape enhancement using alanine scanning. Therefore, we designed 19 alanine mutants of dianthin-30. Each mutant was combined with SO1861, a triterpene saponin isolated from the roots of Saponaria officinalis L., and subjected to a cytotoxicity screening in Neuro-2A cells. Cytotoxic screening revealed that dianthin-30 mutants with lysine substitutions did not impair the endosomal escape enhancement. There was one particular mutant dianthin, Arg24Ala, that exhibited significantly reduced synergistic cytotoxicity in three mammalian cell lines. However, this reduction was not based on an altered interaction with SO1861. It was, rather, due to the impaired endocytosis of dianthin Arg24Ala into the cells.

Genomic progression for local invasion of cutaneous squamous cell carcinoma from the superficial to the deep portion.

Cutaneous squamous cell carcinoma (cSCC) is the second most common skin cancer. While many treatments exist, our understanding of its genomic progression, especially from the epidermis to the deep dermis, remains limited. This study aims to identify genetic mutations associated with the progression of cSCC into the deep dermis, providing insights into its aggressive behavior and high-risk features. We performed high-depth whole-exome sequencing on 12 cSCC tissues, along with paired normal tissues from six patients, using microdissection techniques. The mutational analysis focused on identifying alterations enriched during cSCC progression. Gene Ontology enrichment analysis, immunohistochemical assays, and external single-cell RNA data were utilized for validation. A total of 8863 non-synonymous somatic mutations were identified in 4092 genes across the superficial and deep portions of cSCCs. Analysis of deep portion mutations revealed a significant correlation with gene ontology biological processes, particularly cell junction organization, and cell-cell adhesion. Clonal mutations in these processes were more prevalent in the deep portions, indicating their impact on the cSCC mutation landscape. Genetic evolution analysis identified 29 causal genes associated with dermal invasion in cSCC. We highlight somatic mutations in cSCC, revealing heterogeneity between superficial and deep regions. Altered genes in cell junction organization and cell-cell adhesion emerged as pivotal in dermal invasion. We identified 29 causal genes primarily in deep tumor regions. Our findings emphasize analyzing multiple tumor regions to capture varied mutational landscapes. These insights advance our understanding of cSCC progression, emphasizing genetic and cellular changes during tumor evolution.