SpatialSPM: statistical parametric mapping for the comparison of gene expression pattern images in multiple spatial transcriptomic datasets.

Spatial transcriptomic (ST) techniques help us understand the gene expression levels in specific parts of tissues and organs, providing insights into their biological functions. Even though ST dataset provides information on the gene expression and its location for each sample, it is challenging to compare spatial gene expression patterns across tissue samples with different shapes and coordinates. Here, we propose a method, SpatialSPM, that reconstructs ST data into multi-dimensional image matrices to ensure comparability across different samples through spatial registration process. We demonstrated the applicability of this method by kidney and mouse olfactory bulb datasets as well as mouse brain ST datasets to investigate and directly compare gene expression in a specific anatomical region of interest, pixel by pixel, across various biological statuses. Beyond traditional analyses, SpatialSPM is capable of generating statistical parametric maps, including T-scores and Pearson correlation coefficients. This feature enables the identification of specific regions exhibiting differentially expressed genes across tissue samples, enhancing the depth and specificity of ST studies. Our approach provides an efficient way to analyze ST datasets and may offer detailed insights into various biological conditions.

Immunological subtyping of salivary gland cancer identifies histological origin-specific tumor immune microenvironment.

Gene expression analysis enhances proper cancer subtyping, a better understanding of the molecular characteristics of cancer, and strategies for precision medicine. However, salivary gland cancer (SGC) subtyping remains largely unexplored because of its rarity and diverse histopathological and immunological characteristics. This study aimed to determine whether the histological origin and immunological characteristics of SGC subtypes are intrinsic tumor immunity factors. We performed immune profiling of 94 RNA-seq of SGC tissues and found that the SGCs that originated from the excretory duct (ED), such as the salivary duct and mucoepidermoid carcinomas, exhibit higher immunity than those from the intercalated duct (ID), such as the adenoid cystic and myoepithelial carcinomas, based on the computationally predicted immune score (p < 0.001), immune cell enrichment in the tumor immune microenvironment (TIME) (p < 0.001), T-cell receptor diversity (p < 0.001), and expression of signal I (major histocompatibility complex, MHC, p < 0.001) and signal II (co-stimulatory, p < 0.001 and co-inhibitory, p < 0.001) genes. Further analysis revealed that tolerogenic dendritic cell-induced dysfunctional T-cell populations and T-cell exclusion in the TIME are the major immune evasive mechanisms of the ED-and ID-derived SGCs, respectively.

Characterizing intrinsic molecular features of the immune subtypes of salivary mucoepidermoid carcinoma.

INTRODUCTION: Characterizing the tumor microenvironment (TME) and immune landscape of cancer has been a promising step towards discovering new therapeutic biomarkers and guiding precision medicine; however, its application in mucoepidermoid carcinoma (MEC) has been sparse. Here, we conducted a comprehensive study to understand the properties of the TME and immune profiles of MEC. METHOD: 20 patients with MEC were collected from Yonsei Head and Neck Cancer Centre, Yonsei University, South Korea. Total RNA sequencing was conducted to determine gene expression profiles. Bioinformatic and immunoinformatic analyses were applied to characterize the TME and identify immunophenotypic subgroups, and to investigate the molecular features that explain the distinct phenotypes. RESULTS: The MEC samples were subdivided into two groups, immune hot and immune cold, based on the heterogenous immune cell-infiltration and activation level. The immune-hot subgroup exhibited a higher level of immune activity, including T cell infiltration, cytolytic score, IFN-γ, antigen-presenting machinery, and immune modulator genes. Further characterizing molecular features of two subgroups, downregulation of lipid metabolic regulators, including MLXIPL and FASN, and the migration of chemokines and leukocytes were observed, respectively. And, Group-specific expression of immune checkpoint molecules, such as TIGIT, PD-L2, and CTLA-4, was observed in the immune-hot group, which can be exploited as a potential immunotherapeutic biomarker. CONCLUSIONS: Immunophenotypically heterogeneous MEC subgroups analysis has shown distinctive molecular characteristics and provided potential treatment options. These findings yield new insights into TME of MEC and may help next step to study this uncharted cancer.

HLA-I-restricted CD8+ T cell immunity may accelerate tumorigenesis in conjunction with VHL inactivation.

CD8+ T cells recognize and kill tumor cells with HLA-I tumor antigens in early tumorigenesis, the efficiency of which differs according to antigen-recognition coverage, as shown in earlier tumor onset in HLA-I homozygosity. However, the universality of these associations remains unknown. Here, we assessed the tumor type and driver mutation specificity in the association between tumor onset age and HLA-I zygosity. Statistical analyses identified an unexpected negative relationship in tumors with VHL biallelic loss, wherein HLA-I heterozygosity was associated with earlier tumor onset, while all others showed either no or a positive association. Testing on an independent dataset reproduced the VHL-dependent acceleration of tumor onset in the HLA-I heterozygous group, confirming the association. Further speculation proposed VEGF-A-mediated T cell exhaustion under VHL inactivation as a potential mechanism. Our findings suggest that CD8+ T cell immunity in early tumor suppression can be conditional to the genetic status of tumors and may even lead to adverse consequences.

Tumor microenvironment-aware, single-transcriptome prediction of microsatellite instability in colorectal cancer using meta-analysis.

Detecting microsatellite instability (MSI) in colorectal cancers (CRCs) is essential because it is the determinant of treatment strategies, including immunotherapy and chemotherapy. Yet, no attempt has been made to exploit transcriptomic profile and tumor microenvironment (TME) of it to unveil MSI status in CRC. Hence, we developed a novel TME-aware, single-transcriptome predictor of MSI for CRC, called MAP (Microsatellite instability Absolute single sample Predictor). MAP was developed utilizing recursive feature elimination-random forest with 466 CRC samples from The Cancer Genome Atlas, and its performance was validated in independent cohorts, including 1118 samples. MAP showed robustness and predictive power in predicting MSI status in CRC. Additional advantages for MAP were demonstrated through comparative analysis with existing MSI classifier and other cancer types. Our novel approach will provide access to untouched vast amounts of publicly available transcriptomic data and widen the door for MSI CRC research and be useful for gaining insights to help with translational medicine.

LRRK2 Inhibition Mitigates the Neuroinflammation Caused by TLR2-Specific α-Synuclein and Alleviates Neuroinflammation-Derived Dopaminergic Neuronal Loss.

Evidence suggests that crosstalk occurs between microglial leucine-rich repeat kinase 2 (LRRK2)-a regulator of neuroinflammation-and neuron-released α-synuclein (αSyn)-a promoter of microglial activation and neuroinflammatory responses-in neuroinflammation-mediated Parkinson's disease (PD) progression. Therefore, we examined whether LRRK2 inhibition reduces the responses of microglia to neuroinflammation caused by neuron-released αSyn. We examined the neuroinflammatory responses provoked by Toll-like receptor 2 (TLR2)-positive αSyn of neuronal cells using an LRRK2 inhibitor in the mouse glioma cells, rat primary microglia, and human microglia cell line; and the effects of LRRK2 inhibitor in the co-culture of ectopic αSyn-expressing human neuroblastoma cells and human microglia cells and in mouse models by injecting αSyn. We analyzed the association between LRRK2 activity and αSyn oligomer and TLR2 levels in the substantia nigra tissues of human patients with idiopathic PD (iPD). The TLR2-specific αSyn elevated LRRK2 activity and neuroinflammation, and the LRRK2 inhibitor ameliorated neuroinflammatory responses in various microglia cells, alleviated neuronal degeneration along with neuroinflammation in the co-culture, and blocked the further progression of locomotor failure and dopaminergic neuronal degeneration caused by TLR2-specific αSyn in mice. Furthermore, LRRK2 phosphorylation was increased in patients with iPD showing αSyn-specific high TLR2 level. These results suggest the application of LRRK2 inhibitors as a novel therapeutic approach against αSyn-mediated PD progression.

Comprehensive clinicopathologic, molecular, and immunologic characterization of colorectal carcinomas with loss of three intestinal markers, CDX2, SATB2, and KRT20.

Caudal-type homeobox 2 (CDX2), special AT-rich sequence-binding protein 2 (SATB2), and keratin 20 (KRT20) are frequently used as intestinal epithelium-specific markers in immunohistochemical studies. However, subsets of colorectal carcinomas (CRCs) show loss of these markers. We analyzed The Cancer Genome Atlas data to explore molecular correlates of CDX2, SATB2, and KRT20 genes in 390 CRCs. The decreased mRNA expression of each of the three genes commonly correlated with microsatellite instability-high (MSI-H), CpG island methylator phenotype-high (CIMP-H), BRAF/RNF43 mutations, consensus molecular subtype 1, and high tumor mutational burden. The downregulation of CDX2 or SATB2 was dependent on both MSI-H and CIMP-H, whereas that of KRT20 was more dependent on MSI-H than on CIMP-H. Next, we evaluated the immunohistochemical expression of CDX2, SATB2, and KRT20 in 436 primary CRCs. In contrast to RNA-level expression, decreased expression of CDX2 and SATB2 was more dependent on CIMP-H than on MSI-H. However, consistent with RNA-level expression, decreased expression of KRT20 was more dependent on MSI-H than on CIMP-H. CIMP-H and lymphatic invasion were consistently associated with both CDX2 loss and SATB2 loss in CRCs, regardless of MSI status. In microsatellite stable CRCs, CDX2 loss correlated with BRAF mutation, whereas SATB2 loss was associated with KRAS mutations and decreased T-cell infiltration. Cases with concurrent loss of all three markers were found exclusively in MLH1-methylated MSI-H/CIMP-H CRCs. In conclusion, MSI-H and/or CIMP-H are major common correlates of decreased CDX2/SATB2/KRT20 expression in CRCs, but the specific features associated with the loss of each marker are different in CRCs.

Genomic and transcriptomic characterization of heterogeneous immune subgroups of microsatellite instability-high colorectal cancers.

BACKGROUND: Colorectal cancers (CRCs) with microsatellite instability-high (MSI-H) are hypermutated tumors and are generally regarded as immunogenic. However, their heterogeneous immune responses and underlying molecular characteristics remain largely unexplained. METHODS: We conducted a retrospective analysis of 73 primary MSI-H CRC tissues to characterize heterogeneous immune subgroups. Based on combined tumor-infiltrating lymphocyte (TIL) immunoscore and tertiary lymphoid structure (TLS) activity, MSI-H CRCs were classified into immune-high, immune-intermediate, and immune-low subgroups. Of these, the immune-high and immune-low subgroups were further analyzed using whole-exome and transcriptome sequencing. RESULTS: We found considerable variations in immune parameters between MSI-H CRCs, and immune subgrouping of MSI-H CRCs was performed accordingly. The TIL densities and TLS activities of immune-low MSI-H CRCs were comparable to those of an immune-low or immune-intermediate subgroup of microsatellite-stable CRCs. There were remarkable differences between immune-high and immune-low MSI-H CRCs, including their pathological features (medullary vs mucinous), genomic alterations (tyrosine kinase fusions vs KRAS mutations), and activated signaling pathways (immune-related vs Wnt and Notch signaling), whereas no significant differences were found in tumor mutational burden (TMB) and neoantigen load. The immune-low MSI-H CRCs were subdivided by the consensus molecular subtype (CMS1 vs CMS3) with different gene expression signatures (mesenchymal/stem-like vs epithelial/goblet-like), suggesting distinct immune evasion mechanisms. Angiogenesis and CD200 were identified as potential therapeutic targets in immune-low CMS1 and CMS3 MSI-H CRCs, respectively. CONCLUSIONS: MSI-H CRCs are immunologically heterogeneous, regardless of TMB. The unusual immune-low MSI-H CRCs are characterized by mucinous histology, KRAS mutations, and Wnt/Notch activation, and can be further divided into distinct gene expression subtypes, including CMS4-like CMS1 and CMS3. Our data provide novel insights into precise immunotherapeutic strategies for subtypes of MSI-H tumors.

NTRK oncogenic fusions are exclusively associated with the serrated neoplasia pathway in the colorectum and begin to occur in sessile serrated lesions.

Neurotrophic tropomyosin receptor kinase (NTRK) gene fusions are emerging tissue-agnostic drug targets in malignancies including colorectal carcinomas (CRCs), but their detailed landscape in the context of various colorectal carcinogenesis pathways remains to be investigated. In this study, pan-tropomyosin receptor kinase (TRK) protein expression was assessed by immunohistochemistry (IHC) in retrospectively collected colorectal epithelial tumor tissues, including 441 CRCs [133 microsatellite instability-high (MSI-high) and 308 microsatellite stable (MSS)] and 595 premalignant colorectal lesions (330 serrated lesions and 265 conventional adenomas). TRK-positive cases were then subjected to next-generation sequencing and/or fluorescence in situ hybridization to confirm NTRK rearrangements. TRK IHC positivity was not observed in any of the MSS CRCs, conventional adenomas, traditional serrated adenomas, or hyperplastic polyps, whereas TRK positivity was observed in 11 of 58 (19%) MLH1-methylated MSI-high CRCs, 4 of 23 (17%) sessile serrated lesions with dysplasia (SSLDs), and 5 of 132 (4%) sessile serrated lesions (SSLs). The 11 TRK-positive MSI-high CRCs commonly harbored CpG island methylator phenotype-high (CIMP-high), MLH1 methylation, BRAF/KRAS wild-type, and NTRK1 or NTRK3 fusion (TPM3-NTRK1, TPR-NTRK1, LMNA-NTRK1, SFPQ-NTRK1, ETV6-NTRK3, or EML4-NTRK3). Both NTRK1 or NTRK3 rearrangement and BRAF/KRAS wild-type were detected in all nine TRK-positive SSL(D)s, seven of which demonstrated MSS and/or CIMP-low. TRK expression was selectively observed in distorted serrated crypts within SSLs and was occasionally localized at the base of serrated crypts. NTRK fusions were detected only in SSLs of patients aged ≥50 years, whereas BRAF mutation was found in younger age-onset SSLs. In conclusion, NTRK-rearranged colorectal tumors develop exclusively through the serrated neoplasia pathway and can be initiated from non-dysplastic SSLs without BRAF/KRAS mutations prior to full occurrence of MSI-high/CIMP-high. © 2021 The Pathological Society of Great Britain and Ireland. Published by John Wiley & Sons, Ltd.

Author Correction: NOVA1 induction by inflammation and NOVA1 suppression by epigenetic regulation in head and neck squamous cell carcinoma.

An amendment to this paper has been published and can be accessed via a link at the top of the paper.

An Improved, Assay Platform Agnostic, Absolute Single Sample Breast Cancer Subtype Classifier.

While intrinsic molecular subtypes provide important biological classification of breast cancer, the subtype assignment of individuals is influenced by assay technology and study cohort composition. We sought to develop a platform-independent absolute single-sample subtype classifier based on a minimal number of genes. Pairwise ratios for subtype-specific differentially expressed genes from un-normalized expression data from 432 breast cancer (BC) samples of The Cancer Genome Atlas (TCGA) were used as inputs for machine learning. The subtype classifier with the fewest number of genes and maximal classification power was selected during cross-validation. The final model was evaluated on 5816 samples from 10 independent studies profiled with four different assay platforms. Upon cross-validation within the TCGA cohort, a random forest classifier (MiniABS) with 11 genes achieved the best accuracy of 88.2%. Applying MiniABS to five validation sets of RNA-seq and microarray data showed an average accuracy of 85.15% (vs. 77.72% for Absolute Intrinsic Molecular Subtype (AIMS)). Only MiniABS could be applied to five low-throughput datasets, showing an average accuracy of 87.93%. The MiniABS can absolutely subtype BC using the raw expression levels of only 11 genes, regardless of assay platform, with higher accuracy than existing methods.

Cross-species oncogenic signatures of breast cancer in canine mammary tumors.

Genomic and precision medicine research has afforded notable advances in human cancer treatment, yet applicability to other species remains uncertain. Through whole-exome and transcriptome analyses of 191 spontaneous canine mammary tumors (CMTs) that exhibit the archetypal features of human breast cancers, we found a striking resemblance of genomic characteristics including frequent PIK3CA mutations (43.1%), aberrations of the PI3K-Akt pathway (61.7%), and key genes involved in cancer initiation and progression. We also identified three gene expression-based CMT subtypes, one of which segregated with basal-like human breast cancer subtypes with activated epithelial-to-mesenchymal transition, low claudin expression, and unfavorable disease prognosis. A relative lack of ERBB2 amplification and Her2-enrichment subtype in CMT denoted species-specific molecular mechanisms. Taken together, our results elucidate cross-species oncogenic signatures for a better understanding of universal and context-dependent mechanisms in breast cancer development and provide a basis for precision diagnostics and therapeutics for domestic dogs.

NOVA1 induction by inflammation and NOVA1 suppression by epigenetic regulation in head and neck squamous cell carcinoma.

Neuro-oncological ventral antigen 1 (NOVA1) is known as a neuron-specific pre-mRNA binding splicing factor. Previously, it was shown to be highly upregulated in T lymphocytes, as well as fibroblasts/stromal spindle cells, in tertiary lymphoid tissues formed by the benign immune-inflammatory process, while it was frequently downregulated in tumor cells and other cells within the tumor microenvironment. Here, we sought to identify the mechanisms of NOVA1 modulation in head and neck squamous cell carcinoma (HNSCC). NOVA1 was induced by inflammatory-immune signals within the tumor microenvironment and was suppressed by epigenetic dysregulation, such as that with miR-146. We found attenuated expression of NOVA1 to be associated with non-oropharynx sites such as oral cavity, hypopharynx, and larynx, human papilloma virus (HPV)-negative SCC defined by immunohistochemistry for p16INK4a expression, fewer tumor infiltrating lymphocytes, and poor patient outcomes. Moreover, changes were discovered in epithelial mesenchymal transition-associated markers according to NOVA1 status. This study provides some insights to the underlying mechanism of NOVA1 regulation and suggests that NOVA1 may serve as a prognostic biomarker and potential therapeutic target for HNSCC in the future.

A Therapeutic Strategy for Chemotherapy-Resistant Gastric Cancer via Destabilization of Both ß-Catenin and RAS.

: Treatment of advanced gastric cancer patients with current standard chemotherapeutic agents frequently results in resistance, leading to poor overall survival. However, there has been no success in developing strategies to overcome it. We showed the expression levels of both ß-catenin and RAS were significantly increased and correlated in tissues of 756 gastric cancer (GC) patients and tissues of primary- and acquired-resistance patient-derived xenograft tumors treated with 5-fluorouracil and oxaliplatin modulated with leucovorin (FOLFOX). On the basis of our previous studies, where small molecules to suppress colorectal cancer (CRC) via degrading both ß-catenin and RAS were developed, we tested the effectiveness of KYA1797K, a representative compound functioning by binding axin, in the growth of GC cells. The efficacy test of the drugs using gastric tumor organoids of Apc1638N mice showed that the CD44 and ALDH1A3 cancer stem cell markers were induced by FOLFOX, but not by KYA1797K. KYA1797K also efficiently suppressed tumors generated by re-engrafting the FOLFOX-resistant patient-derived xenograft (PDX) tumors, which also showed resistance to paclitaxel. Overall, the small-molecule approach degrading both ß-catenin and RAS has potential as a therapeutic strategy for treating GC patients resistant to current standard chemotherapies.

Role of cysteinyl leukotriene signaling in a mouse model of noise-induced cochlear injury.

Noise-induced hearing loss is one of the most common types of sensorineural hearing loss. In this study, we examined the expression and localization of leukotriene receptors and their respective changes in the cochlea after hazardous noise exposure. We found that the expression of cysteinyl leukotriene type 1 receptor (CysLTR1) was increased until 3 d after noise exposure and enhanced CysLTR1 expression was mainly observed in the spiral ligament and the organ of Corti. Expression of 5-lipoxygenase was increased similar to that of CysLTR1, and there was an accompanying elevation of CysLT concentration. Posttreatment with leukotriene receptor antagonist (LTRA), montelukast, for 4 consecutive days after noise exposure significantly decreased the permanent threshold shift and also reduced the hair cell death in the cochlea. Using RNA-sequencing, we found that the expression of matrix metalloproteinase-3 (MMP-3) was up-regulated after noise exposure, and it was significantly inhibited by montelukast. Posttreatment with a MMP-3 inhibitor also protected the hair cells and reduced the permanent threshold shift. These findings suggest that acoustic injury up-regulated CysLT signaling in the cochlea and cochlear injury could be attenuated by LTRA through regulation of MMP-3 expression. This study provides mechanistic insights into the role of CysLTs signaling in noise-induced hearing loss and the therapeutic benefit of LTRA.