Pancreatic Cancer Cells Require the Transcription Factor MYRF to Maintain ER Homeostasis

Splicing variants in MYRF cause partial loss of function in the retinal pigment epithelium

Myelin Regulatory Factor (MYRF) regulates retinal pigment epithelial (RPE) development and variants in the C-terminus are linked to isolated nanophthalmos, while loss-of-function variants cause syndromic disease. To define the molecular mechanism of this discrepancy, in vitro and animal studies were performed on a pathogenic C-terminal variant (p.Gly1126fs30* or dG-MYRF). ARPE-19 cells transduced with dG-MYRF revealed reduced target gene expression compared to WT-MYRF, with reduced steady state levels of C-terminal MYRF cleavage product, but intact cleavage and localization. A homozygous humanized MYRF C-terminal ( Myrf humdG/humdG ) mouse model was embryonic lethal by embryonic day (E) 18.5, while humanized wildtype ( Myrf humWT/humWT ) showed normal expression and survival. Bioinformatic analysis on integrated single cell RNA-seq from humanized E17.5 and knockout Rx-Cre;Myrf fl/fl (E15.5 and P0) mice supported shared differentially expressed genes with decreased effect size in Myrf humdG/humdG eyes. These findings, and the viability differences, support that dG-MYRF is a hypomorphic allele. Further, two novel MYRF splicing variants were identified in families with isolated nanophthalmos, with one confirmed to alter 40% of spliced transcripts, creating a nonfunctional isoform. These cases corroborate that isolated nanophthalmos results from hypomorphic alleles of MYRF, supporting a tissue-specific threshold effect and suggests that the C-terminus has unique roles in the RPE.

A mucus production programme promotes classical pancreatic ductal adenocarcinoma

OBJECTIVE

The optimal therapeutic response in cancer patients is highly dependent upon the differentiation state of their tumours. Pancreatic ductal adenocarcinoma (PDA) is a lethal cancer that harbours distinct phenotypic subtypes with preferential sensitivities to standard therapies. This study aimed to investigate intratumour heterogeneity and plasticity of cancer cell states in PDA in order to reveal cell state-specific regulators.

DESIGN

We analysed single-cell expression profiling of mouse PDAs, revealing intratumour heterogeneity and cell plasticity and identified pathways activated in the different cell states. We performed comparative analysis of murine and human expression states and confirmed their phenotypic diversity in specimens by immunolabeling. We assessed the function of phenotypic regulators using mouse models of PDA, organoids, cell lines and orthotopically grafted tumour models.

RESULTS

Our expression analysis and immunolabeling analysis show that a mucus production programme regulated by the transcription factor SPDEF is highly active in precancerous lesions and the classical subtype of PDA - the most common differentiation state. SPDEF maintains the classical differentiation and supports PDA transformation in vivo. The SPDEF tumour-promoting function is mediated by its target genes AGR2 and ERN2/IRE1β that regulate mucus production, and inactivation of the SPDEF programme impairs tumour growth and facilitates subtype interconversion from classical towards basal-like differentiation.

CONCLUSIONS

Our findings expand our understanding of the transcriptional programmes active in precancerous lesions and PDAs of classical differentiation, determine the regulators of mucus production as specific vulnerabilities in these cell states and reveal phenotype switching as a response mechanism to inactivation of differentiation states determinants.

MYRF: A unique transmembrane transcription factor- from proteolytic self-processing to its multifaceted roles in animal development

The Myelin Regulator Factor (MYRF) is a master regulator governing myelin formation and maintenance in the central nervous system. The conservation of MYRF across metazoans and its broad tissue expression suggest it has functions extending beyond the well-established role in myelination. Loss of MYRF results in developmental lethality in both invertebrates and vertebrates, and MYRF haploinsufficiency in humans causes MYRF-related Cardiac Urogenital Syndrome, underscoring its importance in animal development; however, these mechanisms are largely unexplored. MYRF, an unconventional transcription factor, begins embedded in the membrane and undergoes intramolecular chaperone mediated trimerization, which triggers self-cleavage, allowing its N-terminal segment with an Ig-fold DNA-binding domain to enter the nucleus for transcriptional regulation. Recent research suggests developmental regulation of cleavage, yet the mechanisms remain enigmatic. While some parts of MYRF's structure have been elucidated, others remain obscure, leaving questions about how these motifs are linked to its intricate processing and function.

Activation of endogenous retroviruses and induction of viral mimicry by MEK1/2 inhibition in pancreatic cancer

While pancreatic ductal adenocarcinomas (PDACs) are addicted to KRAS-activating mutations, inhibitors of downstream KRAS effectors, such as the MEK1/2 kinase inhibitor trametinib, are devoid of therapeutic effects. However, the extensive rewiring of regulatory circuits driven by the attenuation of the KRAS pathway may induce vulnerabilities of therapeutic relevance. An in-depth molecular analysis of the transcriptional and epigenomic alterations occurring in PDAC cells in the initial hours after MEK1/2 inhibition by trametinib unveiled the induction of endogenous retroviruses (ERVs) escaping epigenetic silencing, leading to the production of double-stranded RNAs and the increased expression of interferon (IFN) genes. We tracked ERV activation to the early induction of the transcription factor ELF3, which extensively bound and activated nonsilenced retroelements and synergized with IRF1 (interferon regulatory factor 1) in the activation of IFNs and IFN-stimulated genes. Trametinib-induced viral mimicry in PDAC may be exploited in the rational design of combination therapies in immuno-oncology.

Systemic deuteration of SCID mice using the water-isotopologue deuterium oxide (D2 O) inhibits tumor growth in an orthotopic bioluminescent model of human pancreatic ductal adenocarcinoma

Since its initial discovery as a natural isotopologue of dihydrogen oxide (1 H2 O), extensive research has focused on the biophysical, biochemical, and pharmacological effects of deuterated water (2 H2 O [D2 O, also referred to as "heavy water"]). Using a panel of cultured human pancreatic ductal adenocarcinoma (PDAC) cells we have profiled (i) D2 O-induced phenotypic antiproliferative and apoptogenic effects, (ii) redox- and proteotoxicity-directed stress response gene expression, and (iii) phosphoprotein-signaling related to endoplasmic reticulum (ER) and MAP-kinase stress response pathways. Differential array analysis revealed early modulation of stress response gene expression in both BxPC-3 and PANC-1 PDAC cells elicited by D2 O (90%; ≤6 h; upregulated: HMOX1, NOS2, CYP2E1, CRYAB, DDIT3, NFKBIA, PTGS1, SOD2, PTGS2; downregulated: RUNX1, MYC, HSPA8, HSPA1A) confirmed by independent RT-qPCR analysis. Immunoblot-analysis revealed rapid (≤6 h) onset of D2 O-induced MAP-kinase signaling (p-JNK, p-p38) together with ER stress response upregulation (p-eIF2α, ATF4, XBP1s, DDIT3/CHOP). Next, we tested the chemotherapeutic efficacy of D2 O-based drinking water supplementation in an orthotopic PDAC model employing firefly luciferase-expressing BxPC-3-FLuc cells in SCID mice. First, feasibility and time course of systemic deuteration (30% D2 O in drinking water; 21 days) were established using time-resolved whole-body proton magnetic resonance imaging and isotope-ratio mass spectrometry-based plasma (D/H)-analysis. D2 O-supplementation suppressed tumor growth by almost 80% with downregulated expression of PCNA, MYC, RUNX1, and HSP70 while increasing tumor levels of DDIT3/CHOP, HO-1, and p-eIF2α. Taken together, these data demonstrate for the first time that pharmacological induction of systemic deuteration significantly reduces orthotopic tumor burden in a murine PDAC xenograft model.

Interferon regulatory factor 1 (IRF1) controls the metabolic programmes of low-grade pancreatic cancer cells

OBJECTIVE

Pancreatic ductal adenocarcinomas (PDACs) include heterogeneous mixtures of low-grade cells forming pseudoglandular structures and compact nests of high-grade cells organised in non-glandular patterns. We previously reported that low-grade PDAC cells display high expression of interferon regulatory factor 1 (IRF1), a pivotal transcription factor of the interferon (IFN) system, suggesting grade-specific, cell-intrinsic activation of IFN responses. Here, we set out to determine the molecular bases and the functional impact of the activation of IFN-regulated responses in human PDACs.

DESIGN

We first confirmed the correlation between glandular differentiation and molecular subtypes of PDAC on the one hand, and the expression of IRF1 and IFN-stimulated genes on the other. We next used unbiased omics approaches to systematically analyse basal and IFN-regulated responses in low-grade and high-grade PDAC cells, as well as the impact of IRF1 on gene expression programmes and metabolic profiles of PDAC cells.

RESULTS

High-level expression of IRF1 in low-grade PDAC cells was controlled by endodermal lineage-determining transcription factors. IRF1-regulated gene expression equipped low-grade PDAC cells with distinctive properties related to antigen presentation and processing as well as responsiveness to IFN stimulation. Notably, IRF1 also controlled the characteristic metabolic profile of low-grade PDAC cells, suppressing both mitochondrial respiration and fatty acid synthesis, which may in part explain its growth-inhibiting activity.

CONCLUSION

IRF1 links endodermal differentiation to the expression of genes controlling antigen presentation and processing as well as to the specification of the metabolic profile characteristic of classical PDAC cells.

Aberrant transcription factors in the cancers of the pancreas

Transcription factors (TFs) are essential for proper activation of gene set during the process of organogenesis, differentiation, lineage specificity. Reactivation or dysregulation of TFs regulatory networks could lead to deformation of organs, diseases including various malignancies. Currently, understanding the mechanism of oncogenesis became necessity for the development of targeted therapeutic strategy for different cancer types. It is evident that many TFs go awry in cancers of the pancreas such as pancreatic ductal adenocarcinoma (PDAC) and pancreatic neuroendocrine neoplasms (PanNENs). These mutated or dysregulated TFs abnormally controls various signaling pathways in PDAC and PanNENs including RTK, PI3K-PTEN-AKT-mTOR, JNK, TGF-β/SMAD, WNT/β-catenin, SHH, NOTCH and VEGF which in turn regulate different hallmarks of cancer. Aberrant regulation of such pathways have been linked to the initiation, progression, metastasis, and resistance in pancreatic cancer. As of today, a number of TFs has been identified as crucial regulators of pancreatic cancer and a handful of them shown to have potential as therapeutic targets in pre-clinical and clinical settings. In this review, we have summarized the current knowledge on the role and therapeutic usefulness of TFs in PDAC and PanNENs.

Proteolytic Activation of Plant Membrane-Bound Transcription Factors

Due to the presence of a transmembrane domain, the subcellular mobility plan of membrane-bound or membrane-tethered transcription factors (MB-TFs) differs from that of their cytosolic counterparts. The MB-TFs are mostly locked in (sub)cellular membranes, until they are released by a proteolytic cleavage event or when the transmembrane domain (TMD) is omitted from the transcript due to alternative splicing. Here, we review the current knowledge on the proteolytic activation mechanisms of MB-TFs in plants, with a particular focus on regulated intramembrane proteolysis (RIP), and discuss the analogy with the proteolytic cleavage of MB-TFs in animal systems. We present a comprehensive inventory of all known and predicted MB-TFs in the model plant Arabidopsis thaliana and examine their experimentally determined or anticipated subcellular localizations and membrane topologies. We predict proteolytically activated MB-TFs by the mapping of protease recognition sequences and structural features that facilitate RIP in and around the TMD, based on data from metazoan intramembrane proteases. Finally, the MB-TF functions in plant responses to environmental stresses and in plant development are considered and novel functions for still uncharacterized MB-TFs are forecasted by means of a regulatory network-based approach.

Maintaining the thyroid gland in mutant thyroglobulin-induced hypothyroidism requires thyroid cell proliferation that must continue in adulthood

Congenital hypothyroidism with biallelic thyroglobulin (Tg protein, encoded by the TG gene) mutation is an endoplasmic reticulum (ER) storage disease. Many patients (and animal models) grow an enlarged thyroid (goiter), yet some do not. In adulthood, hypothyroid TG^cog/cog mice (bearing a Tg-L2263P mutation) exhibit a large goiter, whereas adult WIC rats bearing the TG^rdw/rdw mutation (Tg-G2298R) exhibit a hypoplastic thyroid. Homozygous TG mutation has been linked to thyroid cell death, and cytotoxicity of the Tg-G2298R protein was previously thought to explain the lack of goiter in WIC-TG^rdw/rdw rats. However, recent studies revealed that TG^cog/cog mice also exhibit widespread ER stress–mediated thyrocyte death, yet under continuous feedback stimulation, thyroid cells proliferate in excess of their demise. Here, to examine the relative proteotoxicity of the Tg-G2298R protein, we have used CRISPR–CRISPR-associated protein 9 technology to generate homozygous TG^rdw/rdw knock-in mice in a strain background identical to that of TG^cog/cog mice. TG^rdw/rdw mice exhibit similar phenotypes of defective Tg protein folding, thyroid histological abnormalities, hypothyroidism, and growth retardation. TG^rdw/rdw mice do not show evidence of greater ER stress response or stress-mediated cell death than TG^cog/cog mice, and both mouse models exhibit sustained thyrocyte proliferation, with comparable goiter growth. In contrast, in WIC-TG^rdw/rdw rats, as a function of aging, the thyrocyte proliferation rate declines precipitously. We conclude that the mutant Tg-G2298R protein is not intrinsically more proteotoxic than Tg-L2263P; rather, aging-dependent difference in maintenance of cell proliferation is the limiting factor, which accounts for the absence of goiter in adult WIC-TG^rdw/rdw rats.

Yin Yang 1-induced activation of LINC01133 facilitates the progression of pancreatic cancer by sponging miR-199b-5p to upregulate myelin regulatory factor expression

Increasing evidence has reported that long non-coding RNA (lncRNA) plays a vital role in the development of pancreatic cancer (PC). However, the function and mechanism of LINC01133 in PC tumorigenesis are still unknown. Herein, we found that LINC01133 was highly expressed in PC tissues and cell lines, and LINC01133 knockdown could block the growth and metastasis of PC cells. Besides, upregulated LINC01133 in PC cells was induced by Yin Yang 1 (YY1). Furthermore, LINC01133 directly targeted miR-199b-5p and promoted cancer malignancy by suppressing miR-199b-5p. It was also discovered that myelin regulatory factor (MYRF) was targeted by miR-199b-5p and positively correlated with LINC01133 expression in PC, and LINC01133 modulated PC progression through miR-199b-5p/MYRF pathway. In conclusion, we demonstrated that YY1-mediated the upregulation of LINC0113 increased MYRF expression by sponging miR-199b-5p, resulting in the accelerated development of PC. These findings might offer a novel insight into the development of efficient therapeutics for PC patients.

Potential Immune Biomarker Candidates and Immune Subtypes of Lung Adenocarcinoma for Developing mRNA Vaccines

mRNA vaccines against cancer have advantages in safety, improved therapeutic efficacy, and large-scale production. Therefore, our purpose is to identify immune biomarkers and to analyze immune status for developing mRNA vaccines and selecting appropriate patients for vaccination. We downloaded clinical information and RNA-seq data of 494 LUAD patients from TCGA. LUAD mutational information was hierarchically clustered by NMF package (Version 0.23.0). DeconstructSigs package (Version 1.8.0) and NMF consistency clustering were used to identify mutation signatures. Maftools package (Version 2.6.05) was used to select LUAD-related immune biomarkers. TIMER was used to discuss the correlation between genetic mutations and cellular components. Unsupervised clustering Pam method was used to identify LUAD immune subtypes. Log-rank test and univariate/multivariate cox regression were used to predict the prognosis of immune subtypes. Dimensionality reduction analysis was dedicated to the description of LUAD immune landscape. LUAD patients are classified into four signatures: T >C, APOBEC mutation, age, and tobacco. Then, GPRIN1, MYRF, PLXNB2, SLC9A4, TRIM29, UBA6, and XDH are potential LUAD-related immune biomarker candidates to activate the immune response. Next, we clustered five LUAD-related immune subtypes (IS1–IS5) by prognostic prediction. IS3 showed prolonged survival. The reliability of our five immune subtypes was validated by Thorsson’s results. IS2 and IS4 patients had high tumor mutation burden and large number of somatic mutations. Besides, we identified that immune subtypes of cold immunity (patients with IS2 and IS4) are ideal mRNA vaccination recipients. Finally, LUAD immune landscape revealed immune cells and prognostic conditions, which provides important information to select patients for vaccination. GPRIN1, MYRF, PLXNB2, SLC9A4, TRIM29, UBA6, and XDH are potential LUAD-related immune biomarker candidates to activate the immune response. Patients with IS2 and IS4 might potentially be immunization-sensitive patients for vaccination.

Construction of a novel methylation-related prognostic model for colorectal cancer based on microsatellite status

The present study aimed to construct a novel methylation-related prognostic model based on microsatellite status that may enhance the prognosis of colorectal cancer (CRC) from methylation and microsatellite status perspective. DNA methylation and mRNA expression data with clinical information were downloaded from The Cancer Genome Atlas (TCGA) data set. The samples were divided into microsatellite stability and microsatellite instability group, and CIBERSORT was used to assess the immune cell infiltration characteristics. After identifying the differentially methylated genes and differentially expression genes using R packages, the methylation-driven genes were further identified. Prognostic genes that were used to establish the methylation-related risk score model were generated by the univariate and multivariate Cox regression model. Finally, we established and evaluated the methylation-related prognostic model for CRC patients. A total of 69 MDGs were obtained and three of these genes (MIOX, TH, DKFZP434K028) were selected to construct the prognostic model. Patients in the low-risk score group had a conspicuously better overall survival than those in the high-risk score group (p < .0001). The area under the receiver operating characteristic curve for this model was 0.689 at 3 years, 0.674 at 4 years, and 0.658 at 5 years. The Wilcoxon test showed that higher risk score was associated with higher T stage (p = .01), N stages (p = .0028), metastasis (p = .013), and advanced pathological stage (p = .0013). However, the more instability of microsatellite status, the lower risk score of CRC patients (p = .0048). Our constructed methylation-related prognostic model based on microsatellite status presents potential significance in assessing recurrence risk stratification, tumor staging, and immunotherapy for CRC patients.

Tumor cell heterogeneity and its transcriptional bases in pancreatic cancer: a tale of two cell types and their many variants

Pancreatic ductal adenocarcinoma (PDAC), one of the most highly lethal tumors, is characterized by complex histology, with a massive fibrotic stroma in which both pseudo-glandular structures and compact nests of abnormally differentiated tumor cells are embedded, in different proportions and with different mutual relationships in space. This complexity and the heterogeneity of the tumor component have hindered the development of a broadly accepted, clinically actionable classification of PDACs, either on a morphological or a molecular basis. Here, we discuss evidence suggesting that such heterogeneity can to a large extent, albeit not exclusively, be traced back to two main classes of PDAC cells that commonly coexist in the same tumor: cells that maintained their ability to differentiate toward endodermal, mucin-producing epithelia and epithelial cells unable to form glandular structures and instead characterized by various levels of squamous differentiation and the expression of mesenchymal lineage genes. The underlying gene regulatory networks and how they are controlled by distinct transcription factors, as well as the practical implications of these two different populations of tumor cells, are discussed.

STIM1 Mediates Calcium-Dependent Epigenetic Reprogramming in Pancreatic Cancer

Pancreatic ductal adenocarcinoma (PDAC) displays a dismal prognosis due to late diagnosis and high chemoresistance incidence. For advanced disease stages or patients with comorbidities, treatment options are limited to gemcitabine alone or in combination with other drugs. While gemcitabine resistance has been widely attributed to the levels of one of its targets, RRM1, the molecular consequences of gemcitabine resistance in PDAC remain largely elusive. Here we sought to identify genomic, epigenomic, and transcriptomic events associated with gemcitabine resistance in PDAC and their potential clinical relevance. We found that gemcitabine-resistant cells displayed a coamplification of the adjacent RRM1 and STIM1 genes. Interestingly, RRM1, but not STIM1, was required for gemcitabine resistance, while high STIM1 levels caused an increase in cytosolic calcium concentration. Higher STIM1-dependent calcium influx led to an impaired endoplasmic reticulum stress response and a heightened nuclear factor of activated T-cell activity. Importantly, these findings were confirmed in patient and patient-derived xenograft samples. Taken together, our study uncovers previously unknown biologically relevant molecular properties of gemcitabine-resistant tumors, revealing an undescribed function of STIM1 as a rheostat directing the effects of calcium signaling and controlling epigenetic cell fate determination. It further reveals the potential benefit of targeting STIM1-controlled calcium signaling and its downstream effectors in PDAC. SIGNIFICANCE: Gemcitabine-resistant and some naïve tumors coamplify RRM1 and STIM1, which elicit gemcitabine resistance and induce a calcium signaling shift, promoting ER stress resistance and activation of NFAT signaling.

Functional mechanisms of MYRF DNA-binding domain mutations implicated in birth defects

Myrf is a pleiotropic membrane-bound transcription factor that plays critical roles in diverse organisms, including in oligodendrocyte differentiation, embryonic development, molting, and synaptic plasticity. Upon autolytic cleavage, the Myrf N-terminal fragment enters the nucleus as a homo-trimer and functions as a transcription factor. Homo-trimerization is essential for this function because it imparts DNA-binding specificity and affinity. Recent exome sequencing studies have implicated four de novo MYRF DNA-binding domain (DBD) mutations (F387S, Q403H, G435R, and L479V) in novel syndromic birth defects involving the diaphragm, heart, and the urogenital tract. It remains unknown whether and how these four mutations alter the transcription factor function of MYRF. Here, we studied them by introducing homologous mutations to the mouse Myrf protein. We found that the four DBD mutations abolish the transcriptional activity of the Myrf N-terminal fragment by interfering with its homo-trimerization ability by perturbing the DBD structure. Since the Myrf N-terminal fragment strictly functions as a homo-trimer, any loss-of-function mutation has the potential to act as a dominant negative. We observed that one copy of Myrf-F387S, Myrf-Q403H, or Myrf-L479V, but not Myrf-G435R, was tolerated by the Myrf N-terminal homo-trimer for structural and functional integrity. These data suggest that F387S, Q403H, and L479V cause birth defects by haploinsufficiency, while G435R does so via dominant negative functionality.