Altered Gene Expression Research Articles

Nucleolin lactylation contributes to intrahepatic cholangiocarcinoma pathogenesis via RNA splicing regulation of MADD

Background & aimsIntrahepatic cholangiocarcinoma (iCCA) is a fatal malignancy of the biliary system. The lack of a detailed understanding of oncogenic signaling or global gene expression alterations has impeded clinical iCCA diagnosis and therapy. The role of protein lactylation, a newly unraveled post-translational modification that orchestrates gene expression, remains largely elusive in the pathogenesis of iCCA. MethodsProteomics analysis of clinical iCCA specimens and adjacent tissues was performed to screen for proteins aberrantly lactylated in iCCA. Mass spectrometry, macromolecule interaction and cell behavioral studies were employed to identify the specific lactylation sites on the candidate protein(s) and to decipher the downstream mechanisms responsible for iCCA development, which were subsequently validated using a xenograft tumor model and clinical samples. ResultsNucleolin (NCL), the most abundant RNA-binding protein in the nucleolus, was identified as a functional lactylation target that correlates with iCCA occurrence and progression. NCL was lactylated predominantly at lysine 477 by the acyltransferase P300 in response to a hyperactivity of glycolysis, and promoted the proliferation and invasion of iCCA cells. Mechanistically, lactylated NCL bound to the primary transcript of MAP kinase-activating death domain protein (MADD) and warranted an efficient translation of MADD by circumventing alternative splicing that generates a premature termination codon. NCL lactylation, MADD and subsequent ERK activation promoted xenograft tumor growth, and were found to associate with the overall survival of iCCA patients. ConclusionNCL is lactylated to upregulate MADD through an RNA splicing-dependent mechanism, which potentiates iCCA pathogenesis via the MAPK pathway. Our findings reveal a novel link between metabolic reprogramming and canonical tumor-initiating events, and provide biomarkers that can be potentially used for prognostic evaluation or targeted treatment of iCCA. Impact and implicationIntrahepatic cholangiocarcinoma (iCCA) is a highly aggressive liver malignancy with largely uncharacterized pathogenetic mechanisms. Herein, we demonstrated that glycolysis promotes P300-catalyzed lactylation of NCL, which upregulates MAP kinase-activating death domain protein (MADD) through precise mRNA splicing, and activates ERK signaling to drive iCCA development. These findings unravel a novel link between metabolic rewiring and canonical oncogenic pathways, and provide new biomarkers for prognostic assessment and targeting of clinical iCCA.

Omics-based construction of regulatory variants can be applied to help decipher pig liver-related traits

Genetic variants can influence complex traits by altering gene expression through changes to regulatory elements. However, the genetic variants that affect the activity of regulatory elements in pigs are largely unknown, and the extent to which these variants influence gene expression and contribute to the understanding of complex phenotypes remains unclear. Here, we annotate 90,991 high-quality regulatory elements using acetylation of histone H3 on lysine 27 (H3K27ac) ChIP-seq of 292 pig livers. Combined with genome resequencing and RNA-seq data, we identify 28,425 H3K27ac quantitative trait loci (acQTLs) and 12,250 expression quantitative trait loci (eQTLs). Through the allelic imbalance analysis, we validate two causative acQTL variants in independent datasets. We observe substantial sharing of genetic controls between gene expression and H3K27ac, particularly within promoters. We infer that 46% of H3K27ac exhibit a concomitant rather than causative relationship with gene expression. By integrating GWAS, eQTLs, acQTLs, and transcription factor binding prediction, we further demonstrate their application, through metabolites dulcitol, phosphatidylcholine (PC) (16:0/16:0) and published phenotypes, in identifying likely causal variants and genes, and discovering sub-threshold GWAS loci. We provide insight into the relationship between regulatory elements and gene expression, and the genetic foundation for dissecting the molecular mechanism of phenotypes.

Mice lacking DIO3 exhibit sex-specific alterations in circadian patterns of corticosterone and gene expression in metabolic tissues

Disruption of circadian rhythms is associated with neurological, endocrine and metabolic pathologies. We have recently shown that mice lacking functional type 3 deiodinase (DIO3), the enzyme that clears thyroid hormones, exhibit a phase shift in locomotor activity, suggesting altered circadian rhythm. To better understand the physiological and molecular basis of this phenotype, we used Dio3+/+ and Dio3-/- mice of both sexes at different zeitgeber times (ZTs) and analyzed corticosterone and thyroxine (T4) levels, hypothalamic, hepatic, and adipose tissue expression of clock genes, as well as genes involved in the thyroid hormone action or physiology of liver and adipose tissues. Wild type mice exhibited sexually dimorphic circadian patterns of genes controlling thyroid hormone action, including Dio3. Dio3-/- mice exhibited altered hypothalamic expression of several clock genes at ZT12, but did not disrupt the overall circadian profile. Expression of clock genes in peripheral tissues was not disrupted by Dio3 deficiency. However, Dio3 loss in liver and adipose tissues disrupted circadian profiles of genes that determine tissue thyroid hormone action and physiology. We also observed circadian-specific changes in serum T4 and corticosterone as a result of DIO3 deficiency. The circadian alterations manifested sexual dimorphism. Most notable, the time curve of serum corticosterone was flattened in Dio3-/- females. We conclude that Dio3 exhibits circadian variations, influencing the circadian rhythmicity of thyroid hormone action and physiology in liver and adipose tissues in a sex-specific manner. Circadian disruptions in tissue physiology may then contribute to the metabolic phenotypes of DIO3-deficient mice.

Chronic sleep fragmentation reduces left ventricular contractile function and alters gene expression related to innate immune response and circadian rhythm in the mouse heart

Sleep disorders have emerged as a widespread public health concern, primarily due to their association with an increased risk of developing cardiovascular diseases. Our previous research indicated a potential direct impact of insufficient sleep duration on cardiac remodeling in children and adolescents. Nevertheless, the underlying mechanisms behind the link between sleep fragmentation (SF) and cardiac abnormalities remain unclear. In this study, we aimed to investigate the effects of SF interventions at various life stages on cardiac structure and function, as well as to identify genes associated with SF-induced cardiac dysfunction. To achieve this, we established mouse models of chronic SF and two-week sleep recovery (SR). Our results revealed that chronic SF significantly compromised left ventricular contractile function across different life stages, leading to alterations in cardiac structure and ventricular remodeling, particularly during early life stages. Moreover, microarray analysis of mouse heart tissue identified two significant modules and nine hub genes (Ddx60, Irf9, Oasl2, Rnf213, Cmpk2, Stat2, Parp14, Gbp3, and Herc6) through protein–protein interaction analysis. Notably, the interactome predominantly involved innate immune responses. Importantly, all hub genes lost significance following SR. The second module primarily consisted of circadian clock genes, and real-time PCR validation demonstrated significant upregulation of Arntl, Dbp, and Cry1 after SF, while subsequent SR restored normal Arntl expression. Furthermore, the expression levels of four hub genes (Ddx60, Irf9, Oasl2, and Cmpk2) and three circadian clock genes (Arntl, Dbp, and Cry1) exhibited correlations with structural and functional echocardiographic parameters. Overall, our findings suggest that SF impairs left ventricular contractile function and ventricular remodeling during early life stages, and this may be mediated by modulation of the innate immune response and circadian rhythm. Importantly, our findings suggest that a short period of SR can alleviate the detrimental effects of SF on the cardiac immune response, while the influence of SF on circadian rhythm appears to be more persistent. These findings underscore the importance of good sleep for maintaining cardiac health, particularly during early life stages.

Mendelian randomization analysis using GWAS and eQTL data to investigate the relationship between chronotype and neuropsychiatric disorders and their molecular basis.

Chronotype is a proxy sleep measure that has been associated with neuropsychiatric disorders. By investigating how chronotype influences risk for neuropsychiatric disorders and vice versa, we may identify modifiable risk factors for each phenotype. Here we used Mendelian randomization (MR), to explore causal effects by (1) studying the causal relationships between neuropsychiatric disorders and chronotype and (2) characterizing the genetic components of these phenotypes. Firstly, we investigated if a causal role exists between five neuropsychiatric disorders and chronotype using the largest genome-wide association studies (GWAS) available. Secondly, we integrated data from expression quantitative trait loci (eQTLs) to investigate the role of gene expression alterations on these phenotypes. Evening chronotype was causal for increased risk of schizophrenia and autism spectrum disorder and schizophrenia was causal for a tendency toward evening chronotype. We identified 12 eQTLs where gene expression changes in brain or blood were causal for one of the phenotypes, including two eQTLs for SNX19 in hippocampus and hypothalamus that were causal for schizophrenia. These findings provide important evidence for the complex, bidirectional relationship that exists between a sleep-based phenotype and neuropsychiatric disorders, and use gene expression data to identify causal roles for genes at associated loci.

In vitro aging alters the gene expression and secretome composition of canine adipose-derived mesenchymal stem cells.

Canine adipose-derived mesenchymal stem cells (cAD-MSCs) hold therapeutic promise due to their regenerative potential, particularly within their secretome. However, concerns arise regarding the impact of in vitro cultivation necessitated for storing therapeutic doses, prompting this study to comprehensively explore the impact of in vitro aging on gene expression and secretome composition. The study involved collecting abdominal adipose tissue samples from nine healthy female dogs, from which cAD-MSCs were extracted and cultured. Stem cells were validated through trilineage differentiation assays and flow cytometry immunophenotyping. Gene expression profiling using RT-qPCR array, and cAD-MSCs secretome LC-MS/MS analysis, were conducted at passages 3 and 6 to reveal gene expression and protein composition alterations during in vitro culture. The results demonstrate that the gene expression and secretome composition of cAD-MSCs were impacted by in vitro aging. Among many alterations in gene expression between two passages, two significant downregulations were noted in the MSC-associated PTPRC and IL10 genes. While the majority of proteins and their functional characteristics were shared between passages, the influence of cell aging on secretome composition is highlighted by 10% of proteins being distinctively expressed in each passage, along with 21 significant up- and downregulations. The functional attributes of proteins detected in passage 3 demonstrated a greater inclination towards supporting the regenerative capacity of cAD-MSCs. Moreover, proteins in passage 6 exhibited a noteworthy correlation with the blood coagulation pathway, suggesting an elevated likelihood of coagulation events. To the best of our knowledge, this study presents the first original perspective on the changes in secretome composition that occur when cAD-MSCs age in vitro. Furthermore, it contributes to broadening the currently restricted knowledge base concerning the secretome of cAD-MSCs. In conclusion, our findings show that the regenerative potential of cAD-MSCs, as well as their secretome, may be compromised by in vitro aging. Therefore, our study suggests a preference for earlier passages when considering these cells for therapeutic applications.

Glutamic Acid Modified Gold Nanorod Sensor for the Detection of Calcium ions in Neuronal Cells.

Calcium (Ca2+) ions play a crucial role in the functioning of neurons, governing various aspects of neuronal activity such as rapid modulation and alterations in gene expression. Ca2+ signaling has a significant impact on the development of diseases and the impairment of neuronal functions. Herein, the study reports a Ca2+ ion sensor in neuronal cells using a gold nanorod. The gold nanorod (GA-GNR) conjugated glutamic acid developed in the study was used as a nano-bio probe for the experimental and in vitro detection of calcium. The nanosensor is colloidally stable, preserves plasmonic properties, and shows good viability in neuronal cells, as well as promoting neuron cell line growth. The cytotoxicity and cell penetration of the nanosensor are studied using Raman spectroscopy, brightfield and darkfield microscopy imaging, and MTT assays. The quantification of Ca2+ ions in neuronal cells is determined by monitoring the surface plasmon resonance (SPR) of the GA-GNR. The change in the intensity profile in the presence of Ca2+ incubated neurons was effectively used to develop a portable prototype of an optical Ca2+ sensor, proposing it as a tool for neurodegenerative disease diagnosis and neuromodulation evaluation.

Wemics: A Single-Base Resolution Methylation Quantification Method for Enhanced Prediction of Epigenetic Regulation.

DNA methylation, an epigenetic mechanism that alters gene expression without changing DNA sequence, is essential for organism development and key biological processes like genomic imprinting and X-chromosome inactivation. Despite tremendous efforts in DNA methylation research, accurate quantification of cytosine methylation remains a challenge. Here, a single-base methylation quantification approach is introduced by weighting methylation of consecutive CpG sites (Wemics) in genomic regions. Wemics quantification of DNA methylation better predicts its regulatory impact on gene transcription and identifies differentially methylated regions (DMRs) with more biological relevance. Most Wemics-quantified DMRs in lung cancer are epigenetically conserved and recurrently occurred in other primary cancers from The Cancer Genome Atlas (TCGA), and their aberrant alterations can serve as promising pan-cancer diagnostic markers. It is further revealed that these detected DMRs are enriched in transcription factor (TF) binding motifs, and methylation of these TF binding motifs and TF expression synergistically regulate target gene expression. Using Wemics on epigenomic-transcriptomic data from the large lung cancer cohort, a dozen novel genes with oncogenic potential are discovered that are upregulated by hypomethylation but overlooked by other quantification methods. These findings increase the understanding of the epigenetic mechanism by which DNA methylation regulates gene expression.

Licensing effects of inflammatory factors and TLR ligands on the regenerative capacity of adipose-derived mesenchymal stem cells.

Introduction: Adipose tissue-derived mesenchymal stem cells are promising contributors to regenerative medicine, exhibiting the ability to regenerate tissues and modulate the immune system, which is particularly beneficial for addressing chronic inflammatory ulcers and wounds. Despite their inherent capabilities, research suggests that pretreatment amplifies therapeutic effectiveness. Methods: Our experimental design exposed adipose-derived mesenchymal stem cells to six inflammatory factors for 24h. We subsequently evaluated gene expression and proteome profile alterations and observed the wound closure rate post-treatment. Results: Specific pretreatments, such as IL-1β, notably demonstrated an accelerated wound-healing process. Analysis of gene and protein expression profiles revealed alterations in pathways associated with tissue regeneration. Discussion: This suggests that licensed cells exhibit potentially higher therapeutic efficiency than untreated cells, shedding light on optimizing regenerative strategies using adipose tissue-derived stem cells.

Chronic Opioid Treatment Arrests Neurodevelopment and Alters Synaptic Activity in Human Midbrain Organoids.

Understanding the impact of long-term opioid exposure on the embryonic brain is critical due to the surging number of pregnant mothers with opioid dependency. However, this has been limited by human brain inaccessibility and cross-species differences in animal models. Here, a human midbrain model is established that uses hiPSC-derived midbrain organoids to assess cell-type-specific responses to acute and chronic fentanyl treatment and fentanyl withdrawal. Single-cell mRNA sequencing of 25,510 cells from organoids in different treatment groups reveals that chronic fentanyl treatment arrests neuronal subtype specification during early midbrain development and alters synaptic activity and neuron projection. In contrast, acute fentanyl treatment increases dopamine release but does not significantly alter gene expression related to cell lineage development. These results provide the first examination of the effects of opioid exposure on human midbrain development at the single-cell level.

VoyAGEr, a free web interface for the analysis of age-related gene expression alterations in human tissues

We herein introduce voyAGEr, an online graphical interface to explore age-related gene expression alterations in 49 human tissues. voyAGEr offers a visualisation and statistical toolkit for the finding and functional exploration of sex- and tissue-specific transcriptomic changes with age. In its conception, we developed a novel bioinformatics pipeline leveraging RNA sequencing data, from the GTEx project, encompassing more than 900 individuals. voyAGEr reveals transcriptomic signatures of the known asynchronous ageing between tissues, allowing the observation of tissue-specific age periods of major transcriptional changes, associated with alterations in different biological pathways, cellular composition, and disease conditions. Notably, voyAGEr was created to assist researchers with no expertise in bioinformatics, providing a supportive framework for elaborating, testing and refining their hypotheses on the molecular nature of human ageing and its association with pathologies, thereby also aiding in the discovery of novel therapeutic targets. voyAGEr is freely available at https://compbio.imm.medicina.ulisboa.pt/app/voyAGEr.

VoyAGEr, a free web interface for the analysis of age-related gene expression alterations in human tissues.

We herein introduce voyAGEr, an online graphical interface to explore age-related gene expression alterations in 49 human tissues. voyAGEr offers a visualisation and statistical toolkit for the finding and functional exploration of sex- and tissue-specific transcriptomic changes with age. In its conception, we developed a novel bioinformatics pipeline leveraging RNA sequencing data, from the GTEx project, encompassing more than 900 individuals. voyAGEr reveals transcriptomic signatures of the known asynchronous ageing between tissues, allowing the observation of tissue-specific age periods of major transcriptional changes, associated with alterations in different biological pathways, cellular composition, and disease conditions. Notably, voyAGEr was created to assist researchers with no expertise in bioinformatics, providing a supportive framework for elaborating, testing and refining their hypotheses on the molecular nature of human ageing and its association with pathologies, thereby also aiding in the discovery of novel therapeutic targets. voyAGEr is freely available at https://compbio.imm.medicina.ulisboa.pt/app/voyAGEr.

Effects of microgravity on neural crest stem cells.

Exposure to microgravity (μg) results in a range of systemic changes in the organism, but may also have beneficial cellular effects. In a previous study we detected increased proliferation capacity and upregulation of genes related to proliferation and survival in boundary cap neural crest stem cells (BC) after MASER14 sounding rocket flight compared to ground-based controls. However, whether these changes were due to μg or hypergravity was not clarified. In the current MASER15 experiment BCs were exposed simultaneously to μg and 1 g conditions provided by an onboard centrifuge. BCs exposed to μg displayed a markedly increased proliferation capacity compared to 1 g on board controls, and genetic analysis of BCs harvested 5 h after flight revealed an upregulation, specifically in μg-exposed BCs, of Zfp462 transcription factor, a key regulator of cell pluripotency and neuronal fate. This was associated with alterations in exosome microRNA content between μg and 1 g exposed MASER15 specimens. Since the specimens from MASER14 were obtained for analysis with 1 week's delay, we examined whether gene expression and exosome content were different compared to the current MASER15 experiments, in which specimens were harvested 5 h after flight. The overall pattern of gene expression was different and Zfp462 expression was down-regulated in MASER14 BC μg compared to directly harvested specimens (MASER15). MicroRNA exosome content was markedly altered in medium harvested with delay compared to directly collected samples. In conclusion, our analysis indicates that even short exposure to μg alters gene expression, leading to increased BC capacity for proliferation and survival, lasting for a long time after μg exposure. With delayed harvest of specimens, a situation which may occur due to special post-flight circumstances, the exosome microRNA content is modified compared to fast specimen harvest, and the direct effects from μg exposure may be partially attenuated, whereas other effects can last for a long time after return to ground conditions.

Assessment of the mode of action of perchloroethylene-induced mouse liver tumors.

Perchloroethylene (PCE) is used as a solvent and chemical intermediate. Following chronic inhalation exposure, PCE selectively induced liver tumors in mice. Understanding the mode of action (MOA) for PCE carcinogenesis in mice is important in defining its possible human cancer risk. The proposed MOA is based on the extensive examination of the peer-reviewed studies that have assessed the mouse liver effects of PCE and its major oxidative metabolite trichloroacetic acid (TCA). Similar to PCE, TCA has also been demonstrated to liver tumors selectively in mice following chronic exposure. The Key Events (KE) of the proposed PCE MOA involve oxidative metabolism of PCE to TCA [KE 1]; activation of the peroxisome proliferator-activated receptor alpha (PPARα) [KE 2]; alteration in hepatic gene expression including cell growth pathways [KE 3]; increase in cell proliferation [KE 4]; selective clonal expansion of hepatic preneoplastic foci [KE 5]; and formation of hepatic neoplasms [KE 6]. The scientific evidence supporting the PPARα MOA for PCE is strong and satisfies the requirements for a MOA analysis. The PPARα liver tumor MOA in rodents has been demonstrated not to occur in humans; thus, human liver cancer risk to PCE is not likely.

Gene Expression Patterns Associated with Survival in Glioblastoma.

The aim of this study was to investigate gene expression alterations associated with overall survival (OS) in glioblastoma (GBM). Using the Nanostring nCounter platform, we identified four genes (COL1A2, IGFBP3, NGFR, and WIF1) that achieved statistical significance when comparing GBM with non-neoplastic brain tissue. The four genes were included in a multivariate Cox Proportional Hazard model, along with age, extent of resection, and O6-methylguanine-DNA methyltransferase (MGMT) promotor methylation, to create a unique glioblastoma prognostic index (GPI). The GPI score inversely correlated with survival: patient with a high GPI had a median OS of 7.5 months (18-month OS = 9.7%) whereas patients with a low GPI had a median OS of 20.1 months (18-month OS = 54.5%; log rank p-value = 0.004). The GPI score was then validated in 188 GBM patients from The Cancer Genome Atlas (TCGA) from a national data base; similarly, patients with a high GPI had a median OS of 10.5 months (18-month OS = 12.4%) versus 16.9 months (18-month OS = 41.5%) for low GPI (log rank p-value = 0.0003). We conclude that this novel mRNA-based prognostic index could be useful in classifying GBM patients into risk groups and refine prognosis estimates to better inform treatment decisions or stratification into clinical trials.

Abstract 2523: Predictive biomarkers for neoadjuvant FLOT4 chemotherapy response in gastroesophageal cancer: Results of multi-omics approach

Abstract Background: The FLOT4 chemotherapy regimen is a standard preoperative treatment for gastroesophageal junction (GEJ) and gastric cancer (GC), yet the response to this treatment is variable. Identification of predictive biomarkers is crucial for optimizing therapeutic strategies. This study aims to identify predictive biomarkers for FLOT4 treatment efficacy using a multi-omics approach. Methods: We conducted a preliminary differential gene expression analysis in diagnostic tissue pre-treatment and surgical resection post-treatment in six patients (two GEJ and four GC) who all exhibited a partial pathological response (G2) to neoadjuvant FLOT4 chemotherapy. Utilizing DESeq2 for differential gene expression and David pathway analysis for functional insights, we have begun to unravel the complex molecular response to treatment. Upcoming assays include blood gDNA and tumor tissue whole exome sequencing (WES), comprehensive methylation profiling, miRNA characterization, and proteomics via LC-MSMS. Results: Our preliminary data suggest alterations in gene expression associated with cell adhesion, extracellular matrix interaction, and immune response pathways. Three common genes (CXCL2, CXCL3, CXCL1) were identified among upregulated pathways in GEJ patients, while one common gene (IL1A) was observed in GC patients. These findings suggest a complex network of biological processes influenced by FLOT4. In-depth results, including whole-exome sequencing, methylation patterns, and miRNA profiles, will be completed and presented at the conference. Conclusion: These initial findings propose that changes in gene expression related to cellular architecture and immune response are indicative of a partial response to FLOT4 in both GEJ and GC patients. The identified genes and pathways may serve as potential biomarkers for predicting treatment outcomes. Future Directions: By the time of the conference, we anticipate including a more extensive patient cohort with complete response or no response to FLOT4, and comparative analyses with a FOLFOX treatment group. The integration of forthcoming WES, methylation, miRNA, and proteomics data will provide a more comprehensive molecular predictive model. Keywords: FLOT4 chemotherapy, gastroesophageal junction cancer, gastric cancer, predictive biomarkers, differential gene expression. Citation Format: Amira Al-Fahdi, Shoaib Al Zadjali, Ibrahim Al Haddabi, Anoopa Pullanhi, Muna Al Jabri, Aida Al Yayahee, Nada Al Shuaili, Nahad Al Mahrouqi, Mansour Al Moundhri. Predictive biomarkers for neoadjuvant FLOT4 chemotherapy response in gastroesophageal cancer: Results of multi-omics approach [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2024; Part 1 (Regular Abstracts); 2024 Apr 5-10; San Diego, CA. Philadelphia (PA): AACR; Cancer Res 2024;84(6_Suppl):Abstract nr 2523.

Abstract 1670: Mechanistic insight into mTOR signaling identifies RNA splicing machinery as a potential therapeutic target for cancer

Abstract Aberrant signal transduction regulating cell growth and proliferation is a hallmark of many cancers. Activation of mTOR signaling is common and is an attractive therapeutic target with mechanistically different mTOR inhibitors in clinical use. One consequence of exposure to the mTOR inhibitor rapamycin is an alteration of gene expression through a global increase in alternative exon inclusion (Cheng et al. 2022, doi.org/10.3390/ijms232012416). To examine the mechanisms underlying this effect, we utilized the TempO-Seq targeted sequencing platform to profile both the whole transcriptome and a selection of alternatively spliced RNA targets in the rapamycin sensitive and resistant breast cancer cell lines MCF-7 and MDA MB 231 following exposure to 13 mTOR inhibitors, some closely related and some from distinct chemical classes. Recent publications have demonstrated that high quality dose response data can be obtained using this sequencing platform (Harrill et al. 2021, doi: 10.1093/toxsci/kfab009; Ramaiahgari et al. 2019, doi: 10.1093/toxsci/kfz065) and TempO-Seq has been widely adopted in toxicology for deriving compound-specific expression signatures to identify mechanisms of action and long-term toxicity. Treated cells were also profiled using TempO-LINC, a novel high-throughput single-cell assay based on the same platform. We observed that certain mTOR inhibitors induced both common and cell-specific differential expression as well as widespread alternative exon inclusion, as expected. Intriguingly, transcripts with the most significant changes in splicing upon treatment were those encoding members of the SR family of splicing factors, often overexpressed or amplified in breast cancers (Park, et al. 2019, doi: 10.1016/j.celrep.2019.10.110). This is significant because SR protein expression is autoregulated by alternative inclusion of an exon containing a premature stop codon in a feedback loop leading to nonsense-mediated decay (NMD) and lowered protein levels. Indeed, NMD is known to be inhibited by rapamycin (Martinez-Nunez, et al. 2017, doi: 10.1093/nar/gkw1109), which would lead to the observed accumulation of these transcripts. Thus, monitoring a targeted selection of alternative splicing events and regulators revealed mTOR regulation of mRNA isoform abundance consistent with regulation of the NMD RNA turnover process and with overexpression of SR proteins in breast cancers. Importantly, targeted bulk and single-cell TempO-Seq assays enable efficient sequencing and simplified data analysis without the expense and bioinformatic burden of RNAseq, making them particularly well-suited for drug discovery. These results support the use of this approach for further characterization of cell-specific effects of mTOR inhibitors and their mechanisms of action. Citation Format: Gaelle Rondeau, Megan Opichka, Kevin White, Joanne M. Yeakley, Dennis J. Eastburn, Bruce Seligmann. Mechanistic insight into mTOR signaling identifies RNA splicing machinery as a potential therapeutic target for cancer [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2024; Part 1 (Regular Abstracts); 2024 Apr 5-10; San Diego, CA. Philadelphia (PA): AACR; Cancer Res 2024;84(6_Suppl):Abstract nr 1670.

Abstract 6092: Three-dimensional multi-omic analysis of early invasion of human pancreatic ductal adenocarcinoma from IPMN

Abstract Background: Pancreatic ductal adenocarcinoma (PDAC) is a deadly cancer, with only 11% of patients surviving beyond 5 years. Intraductal papillary mucinous neoplasm (IPMN), a precancerous lesion and the most common type of neoplastic pancreatic cyst, presents a critical opportunity for cancer interception, but the drivers of malignant progression in IPMN are still largely unknown. Although the immunosuppressive microenvironment of PDAC has been well characterized, the timing of the development of these immune alterations in IPMNs is not well studied and could provide a rational foundation for cancer immuno-prevention. Moreover, the inception of invasive PDAC and its spread into the microenvironment from IPMN is complex, multifactorial, and, due to its unique anatomy, challenging to analyze with traditional 2-dimensional visualization. We sought to capture the precise point of transition from IPMN to invasive carcinoma utilizing CODA, a supervised deep learning tool for three-dimensional (3D) reconstruction of serially sectioned human tissue, in order to quantitatively assess the molecular and cellular alterations associated with malignant progression. Method: Formalin-fixed paraffin-embedded (FFPE) tissue blocks with PDAC arising from IPMN were serially sectioned, every third slide was stained with H&E, and digitized at 20x magnification. Following pathologist-guided annotations on a subset of H&E slides, CODA generated 3D models of each tissue block, including automated annotation of 9 pancreatic tissue components. Using these annotations, we identified the transition from IPMN to invasive carcinoma in each sample and selected regions for multi-omic profiling based on quantitative features of the cellular microenvironment. Multi-omic profiling of IPMN, transition zone, and PDAC included laser capture microdissection followed by whole exome sequencing to identify somatic DNA alterations, as well as spatial transcriptomics to identify alterations in gene expression in neoplastic cells and spatial proteomics to identify cellular alterations in the surrounding microenvironment. Results: CODA generated accurate annotation on all the H&E slides from each case and robustly identified regions of interest for multi-omic profiling. Spatial proteomic profiling using imaging mass cytometry revealed a decrease in T cell density as IPMNs transition into PDAC. Sub- clustering of the T cell compartment identified decreases in multiple T cell subsets, including activated cytotoxic T cells and helper T cells, in the transition zone compared to IPMN. Conclusion: In this study, we integrated multi-omic profiling with CODA-generated high-resolution 3D tissue maps to identify molecular and cellular drivers of malignant progression in human PDAC. Citation Format: Shalini Datta, Sarah M. Shin, Jessie Kanacharoen, Michael Johannes Pflüger, Katsuya Hirose, André Forjaz, Sarah Graham, Pei-Hsun Wu, Ralph H. Hruban, Denis Wirtz, Won Jin Ho, Ashley L. Kiemen, Laura D. Wood. Three-dimensional multi-omic analysis of early invasion of human pancreatic ductal adenocarcinoma from IPMN [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2024; Part 1 (Regular Abstracts); 2024 Apr 5-10; San Diego, CA. Philadelphia (PA): AACR; Cancer Res 2024;84(6_Suppl):Abstract nr 6092.

Abstract 6930: Spatial analysis of multi-ancestral tumors to understand how intercellular signaling impacts tumor invasion

Abstract Colorectal cancer (CRC) persists as a leading cause of cancer-related deaths in the United States. Intratumoral heterogeneity stems from a multi-ancestral origin or clonal divergence and contributes to treatment resistance. Unfortunately, tumors in existing mouse models are largely homogenous, rendering predictions of drug efficacy inaccurate. To investigate whether signaling between distinct neoplastic clones in a tumor and cells in the surrounding microenvironment impacts CRC biology, we performed spatial transcriptomic analysis on samples collected from a transgenic mouse model previously established in the Halberg laboratory. This model (referred to as B6.FCTG3KM) develops three types of tumors: homotypic Adenomatous Polyposis Coli (APC)-deficient tumors that express red fluorescent protein (RFP), homotypic APC-deficient tumors that express constitutively active phosphoinositide 3-kinase (PI3K) and green fluorescent protein (GFP), and heterotypic tumors. Previous characterization of the model clearly demonstrated that heterotypic tumors are significantly more aggressive than either type of homotypic tumor. Analysis was conducted on four serially sectioned slides containing three tumors and adjacent normal tissue obtained from the colon of a male B6.FCTG3KM mouse using the 10x Genomics Visium CytAssist platform with the mouse whole transcriptome probe set. Collected tumors were staged using traditional histopathology techniques by a board-certified pathologist (Matkowskyj). Immunohistochemistry and immunofluorescence were used to assess β-catenin localization and to re-probe for RFP and GFP expression. Sequencing was performed on the Illumina NovaSeqX platform to an average depth of 255,528,294 reads per sample. Data analysis was performed using the 10x Genomics SpaceRanger and Kendziorski laboratory R/SpatialView packages. A fold-change threshold of 1.5 and an adjusted p-value threshold of 0.05 were used when establishing marker genes for each of the 8 clusters identified. One distinct cluster was identified by marker genes including Notum, Wif1, and Dkk2, factors previously implicated in CRC tumorigenesis, and corresponded to the center of each region of surface-level dysplastic epithelia. The remainder of the surface-level dysplastic epithelia corresponded to a second cluster, with marker genes including those coding for TRIM family proteins (Trim5, Trim10, Trim40, Trim12c), and matrix metalloproteinases (Mmp12, Mmp9, Mmp10, Mmp13). A third cluster corresponded to the desmoplastic stroma underlying each neoplasm. Ongoing analyses are underway to identify clonal interactions that alter gene expression and tumor aggressiveness. Citation Format: Mary M. Stangis, Lingxin Cheng, Dawn Albrecht, Santina Snow, Tessa Wirtz, Kristina A. Matkowskyj, Huy Q. Dinh, Christina Kendziorski, Richard B. Halberg. Spatial analysis of multi-ancestral tumors to understand how intercellular signaling impacts tumor invasion [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2024; Part 1 (Regular Abstracts); 2024 Apr 5-10; San Diego, CA. Philadelphia (PA): AACR; Cancer Res 2024;84(6_Suppl):Abstract nr 6930.

Abstract 6139: Genetic modifiers of KRAS-mutant colorectal cancer

Abstract Somatic driver mutations in KRAS are found in ~40% of colorectal cancers (CRCs) and are associated with worse outcomes. In the US, individuals of African ancestry with CRC are more likely to have tumors with KRAS mutations than individuals of European ancestry, even after considering differences in socioeconomic status and other risk factors. Germline variants can influence somatic mutation events and may provide a selective advantage such that a mutated cell is more likely to progress to a cancer. Germline variants have been found to associate with somatic mutations in a variety of tumors. Based on ancestry differences in KRAS tumor mutation frequency, we hypothesized that germline variants altering molecular pathways that support KRAS tumorigenesis will associate with KRAS-mutant status. To test this, we performed a two-step association study using KRAS mutation status as our phenotype in individuals with microsatellite stable CRC. First, we conducted a discovery analysis using genome-wide germline genotypes and somatic mutation data from 6,386 non-Hispanic White CRC patients from the Cancer Genome Atlas and the GECCO consortium, considering sex and tumor location as covariates. To select variants for additional study, we used in-silico tools, such as expression quantitative trait loci (eQTL), to predict how associated variants may alter gene expression and chromatin. After reviewing relevant genomic features, excluding SNVs with minor allele frequency <10% and pruning SNVs in linkage disequilibrium (r2>0.75), we identified candidate loci for a multi-ancestry independent validation analysis. In the discovery analysis, no variants met genome-wide significance thresholds (p-value <5x10-8), though 3 SNVs associated with KRAS mutation status with p-values <1x10-6 and 50 with p-values <1x10-5. Of these 50 SNVs, 21 met criteria for validation. We chose 80 additional SNVs (p-value <1x10-4) with predicted functional effect, for a total of 101 SNVs for analyses. In our preliminary validation analysis of 1,487 individuals grouped by self-reported race, no variants had a significant association in all populations. The rs726800 variant showed nominal association with KRAS mutation status (p-value=0.028) in self-reported non-Hispanic White individuals. This variant is within potential cis-regulatory distance of HS3ST3A1 and COX10-DT, both of which have been reported as cancer biomarkers. Nominal associations were identified in self-reported Black (n=7 SNVs) and Hispanic or Latino (n=3 SNVs) individuals, though both groups had limited sample sizes. Additional validation samples are being analyzed, as well as analysis with adjustment for genetic ancestry. Future functional studies of top candidate variants will add mechanistic insight. In summary, we identified germline SNVs that associate with KRAS somatic mutations in CRC and may inform the cellular context that supports development of colorectal tumors with KRAS mutations. Citation Format: Nijole P. Tjader, Johnny Ramroop, Tanish Gandhi, Peter T. Campbell, Andrew T. Chan, Steven Gallinger, Graham G. Giles, Marc J. Gunter, Tabitha A. Harrison, Michael Hoffmeister, Polly A. Newcomb, Shuji Ogino, Amanda I. Phipps, Conghui Qu, Robert E. Schoen, Andrew Pellatt, Michael O. Woods, Bethany Van Guelpen, Patrick Stevens, Heather Hampel, Ulrike Peters, Joseph P. McElroy, Amanda E. Toland. Genetic modifiers of KRAS-mutant colorectal cancer [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2024; Part 1 (Regular Abstracts); 2024 Apr 5-10; San Diego, CA. Philadelphia (PA): AACR; Cancer Res 2024;84(6_Suppl):Abstract nr 6139.