Transposable elements alter gene expression and may impact response to cisplatin therapy in ovarian cancer.

Cisplatin is widely employed for cancer treatment; therefore, understanding resistance to this drug is critical for therapeutic practice. While studies have delved into differential gene expression in the context of cisplatin resistance, findings remain somewhat scant. We performed a comprehensive investigation of Transposable Elements (TEs) expression and their impact in host genes in two cisplatin-treated ovarian cancer cell lines. RNA-seq, ATAC-seq, and in-depth bioinformatics analysis were used to compare cisplatin-sensitive and -resistant ovarian cancer cell lines. Our results reveal that cisplatin therapy alters the expression of protein-coding genes, but also key TEs, including LINE1, Alu, and endogenous retroviruses, in both cisplatin-sensitive and -resistant cell lines. By co-expressing with downstream genes or by creating chimeric transcripts with host genes at their insertion sites, these TEs seem to control the expression of protein-coding genes, including tumor-related genes. Our model uncovers TEs influencing the expression of cancer genes and cancer pathways. Collectively, our findings indicate that TEs alterations associated with cisplatin treatment occur in critical cancer genes and cellular pathways synergically. This research highlights the importance of considering the entire spectrum of transcribed elements in the genome, especially TE expression, for a complete understanding of complex models like cancer response to treatment.

Population Dynamics is a Cancer Driver.

Most tissues are continuously renovated through the division of stem cells and the death of old or damaged cells, which is known as cell turnover rate (CTOR). Despite being in steady state, tissues have different population dynamics and leading to diverse clonality levels. Here, we propose and test that cell population dynamics can be a cancer driver. We employed the evolutionary software esiCancer to show that CTOR, within a range comparable to what is observed in human tissues, can amplify the risk of a mutation due to ancestral selection (ANSEL). In a high CTOR tissue, a mutated ancestral cell is likely to be selected and persist over generations, which leads to a scenario of elevated ANSEL profile, characterized by few niches of large clones, which does not occur in low CTOR. We found that CTOR is significantly associated with the risk of developing cancer, even when correcting for mutation load, indicating that population dynamics per se is a cancer driver. This concept is central to understanding cancer risk and for the design of new therapeutic interventions that minimize the contribution of ANSEL in cancer growth.

Mutational signatures and increased retrotransposon insertions in xeroderma pigmentosum variant skin tumors.

Xeroderma pigmentosum variant (XP-V) is an autosomal recessive disease with an increased risk of developing cutaneous neoplasms in sunlight-exposed regions. These cells are deficient in the translesion synthesis (TLS) DNA polymerase eta, responsible for bypassing different types of DNA lesions. From the exome sequencing of 11 skin tumors of a genetic XP-V patients' cluster, classical mutational signatures related to sunlight exposure, such as C>T transitions targeted to pyrimidine dimers, were identified. However, basal cell carcinomas also showed distinct C>A mutation spectra reflecting a mutational signature possibly related to sunlight-induced oxidative stress. Moreover, four samples carry different mutational signatures, with C>A mutations associated with tobacco chewing or smoking usage. Thus, XP-V patients should be warned of the risk of these habits. Surprisingly, higher levels of retrotransposon somatic insertions were also detected when the tumors were compared with non-XP skin tumors, revealing other possible causes for XP-V tumors and novel functions for the TLS polymerase eta in suppressing retrotransposition. Finally, the expected high mutation burden found in most of these tumors renders these XP patients good candidates for checkpoint blockade immunotherapy.

Three-dimensional nuclear architecture distinguishes thyroid cancer histotypes.

Molecular markers can serve as diagnostic tools to support pathological analysis in thyroid neoplasms. However, because the same markers can be observed in some benign thyroid lesions, additional approaches are necessary to differentiate thyroid tumor subtypes, prevent overtreatment and tailor specific clinical management. This applies particularly to the recently described variant of thyroid cancer referred to as noninvasive follicular thyroid neoplasm with papillary-like nuclear features (NIFTP). This variant has an estimated prevalence of 4.4% to 9.1% of all papillary thyroid carcinomas worldwide. We studied 60 thyroid lesions: 20 classical papillary thyroid carcinoma (CPTC), 20 follicular variant of PTC (FVPTC) and 20 NIFTP. We examined morphological and molecular features to identify parameters that can differentiate NIFTP from the other PTC subtypes. When blindly investigating the nuclear architecture of thyroid neoplasms, we observed that NIFTP has significantly longer telomeres than CPTC and FVPTC. Super-resolved 3D-structured illumination microscopy demonstrated that NIFTP is heterogeneous and that its nuclei contain more densely packed DNA and smaller interchromatin spaces than CPTC and FVPTC, a pattern that resembles normal thyroid tissue. These data are consistent with the observed indolent biological behavior and favorable prognosis associated with NIFTP, which lacks BRAFV600E mutations. Of note, next-generation thyroid oncopanel sequencing was unable to distinguish the thyroid cancer histotypes in our study cohort. In summary, our data suggest that 3D nuclear architecture can be a powerful analytical tool to diagnose and guide clinical management of NIFTP.

Tissue-resident glial cells associate with tumoral vasculature and promote cancer progression.

Cancer cells are embedded within the tissue and interact dynamically with its components during cancer progression. Understanding the contribution of cellular components within the tumor microenvironment is crucial for the success of therapeutic applications. Here, we reveal the presence of perivascular GFAP+/Plp1+ cells within the tumor microenvironment. Using in vivo inducible Cre/loxP mediated systems, we demonstrated that these cells derive from tissue-resident Schwann cells. Genetic ablation of endogenous Schwann cells slowed down tumor growth and angiogenesis. Schwann cell-specific depletion also induced a boost in the immune surveillance by increasing tumor-infiltrating anti-tumor lymphocytes, while reducing immune-suppressor cells. In humans, a retrospective in silico analysis of tumor biopsies revealed that increased expression of Schwann cell-related genes within melanoma was associated with improved survival. Collectively, our study suggests that Schwann cells regulate tumor progression, indicating that manipulation of Schwann cells may provide a valuable tool to improve cancer patients' outcomes.

Suspected Germline TP53 Variants and Clonal Hematopoiesis of Indeterminate Potential: Lessons Learned From a Molecular Tumor Board.

OBJECTIVE: Li-Fraumeni syndrome (LFS) is a pan-cancer predisposition syndrome caused by germline pathogenic variants in the gene TP53. The interpretation of TP53 variants in clinical scenarios outside the classic LFS criteria may be challenging. Here, we report a patient affected by 2 primary cancers at later ages, who harbored a likely pathogenic TP53 at low allele frequency detected in a blood sample. METHODS: The Molecular Tumor Board committee at our institution revisited the case of a patient who was enrolled in a research protocol for the investigation of genetic conditions associated with neuroendocrine tumors. Clinical, familial, and molecular data were reviewed. The patient received germline testing using a next generation sequencing multi-gene panel and was incidentally found to harbor a TP53 likely pathogenic variant, with 22% of variant allele fraction. Additional samples, including a second blood sample, oral swab, and saliva, were collected for DNA analysis. A new TP53 sequencing round was performed with the attempt to distinguish between a true constitutional germline variant and a somatically acquired variant due to aberrant clonal expansion of bone marrow precursors. RESULTS: Patient's personal and familial history of cancer did not meet classic nor Chompret LFS criteria. Environmental risk factors for cancer were identified, such as alcohol abuse and tobacco exposure. The TP53 variant initially found in the next-generation sequencing was confirmed by Sanger sequencing in the previous DNA sample extracted from blood for the first analysis and in a second blood sample collected 6 years later. The TP53 variant was not detected in the DNA extracted from the oral swab and saliva samples. CONCLUSION: Considering the low TP53 variant allele fraction in blood, absence of variant detection in oral swab and saliva samples, the lack of LFS clinical criteria, and history of exposure to environmental risk factors for cancer, the main hypothesis for this case was aberrant clonal expansion due to clonal hematopoiesis. Oncologists should interpret TP53 findings during germline testing with caution.

The paralogues MAGOH and MAGOHB are oncogenic factors in high-grade gliomas and safeguard the splicing of cell division and cell cycle genes.

The exon junction complex (EJC) plays key roles throughout the lifespan of RNA and is particularly relevant in the nervous system. We investigated the roles of two EJC members, the paralogs MAGOH and MAGOHB, with respect to brain tumour development. High MAGOH/MAGOHB expression was observed in 14 tumour types; glioblastoma (GBM) showed the greatest difference compared to normal tissue. Increased MAGOH/MAGOHB expression was associated with poor prognosis in glioma patients, while knockdown of MAGOH/MAGOHB affected different cancer phenotypes. Reduced MAGOH/MAGOHB expression in GBM cells caused alterations in the splicing profile, including re-splicing and skipping of multiple exons. The binding profiles of EJC proteins indicated that exons affected by MAGOH/MAGOHB knockdown accumulated fewer complexes on average, providing a possible explanation for their sensitivity to MAGOH/MAGOHB knockdown. Transcripts (genes) showing alterations in the splicing profile are mainly implicated in cell division, cell cycle, splicing, and translation. We propose that high MAGOH/MAGOHB levels are required to safeguard the splicing of genes in high demand in scenarios requiring increased cell proliferation (brain development and GBM growth), ensuring efficient cell division, cell cycle regulation, and gene expression (splicing and translation). Since differentiated neuronal cells do not require increased MAGOH/MAGOHB expression, targeting these paralogs is a potential option for treating GBM.

Depletion of Ric-8B leads to reduced mTORC2 activity.

mTOR, a serine/threonine protein kinase that is involved in a series of critical cellular processes, can be found in two functionally distinct complexes, mTORC1 and mTORC2. In contrast to mTORC1, little is known about the mechanisms that regulate mTORC2. Here we show that mTORC2 activity is reduced in mice with a hypomorphic mutation of the Ric-8B gene. Ric-8B is a highly conserved protein that acts as a non-canonical guanine nucleotide exchange factor (GEF) for heterotrimeric Gαs/olf type subunits. We found that Ric-8B hypomorph embryos are smaller than their wild type littermates, fail to close the neural tube in the cephalic region and die during mid-embryogenesis. Comparative transcriptome analysis revealed that signaling pathways involving GPCRs and G proteins are dysregulated in the Ric-8B mutant embryos. Interestingly, this analysis also revealed an unexpected impairment of the mTOR signaling pathway. Phosphorylation of Akt at Ser473 is downregulated in the Ric-8B mutant embryos, indicating a decreased activity of mTORC2. Knockdown of the endogenous Ric-8B gene in cultured cell lines leads to reduced phosphorylation levels of Akt (Ser473), further supporting the involvement of Ric-8B in mTORC2 activity. Our results reveal a crucial role for Ric-8B in development and provide novel insights into the signals that regulate mTORC2.

The Clinical and Molecular Profile of Lung Cancer Patients Harboring the TP53 R337H Germline Variant in a Brazilian Cancer Center: The Possible Mechanism of Carcinogenesis.

In southern and southeastern Brazil, the TP53 founder variant c.1010G>A (R337H) has been previously documented with a prevalence of 0.3% within the general population and linked to a heightened incidence of lung adenocarcinomas (LUADs). In the present investigation, we cover clinical and molecular characterizations of lung cancer patients from the Brazilian Li-Fraumeni Syndrome Study (BLISS) database. Among the 175 diagnosed malignant neoplasms, 28 (16%) were classified as LUADs, predominantly occurring in females (68%), aged above 50 years, and never-smokers (78.6%). Significantly, LUADs manifested as the initial clinical presentation of Li-Fraumeni Syndrome in 78.6% of cases. Molecular profiling was available for 20 patients, with 14 (70%) revealing EGFR family alterations. In total, 23 alterations in cancer driver genes were identified, comprising 7 actionable mutations and 4 linked to resistance against systemic treatments. In conclusion, the carriers of TP53 R337H demonstrate a predisposition to LUAD development. Furthermore, our results indicate that environmental pollution potentially impacts the carcinogenesis of lung tumors in the carriers of TP53 R337H.

sideRETRO: a pipeline for identifying somatic and polymorphic insertions of processed pseudogenes or retrocopies.

MOTIVATION: Retrocopies or processed pseudogenes are gene copies resulting from mRNA retrotransposition. These gene duplicates can be fixed, somatically inserted or polymorphic in the genome. However, knowledge regarding unfixed retrocopies (retroCNVs) is still limited, and the development of computational tools for effectively identifying and genotyping them is an urgent need. RESULTS: Here, we present sideRETRO, a pipeline dedicated not only to detecting retroCNVs in whole-genome or whole-exome sequencing data but also to revealing their insertion sites, zygosity and genomic context and classifying them as somatic or polymorphic events. We show that sideRETRO can identify novel retroCNVs and genotype them, in addition to finding polymorphic retroCNVs in whole-genome and whole-exome data. Therefore, sideRETRO fills a gap in the literature and presents an efficient and straightforward algorithm to accelerate the study of bona fide retroCNVs. AVAILABILITY AND IMPLEMENTATION: sideRETRO is available at https://github.com/galantelab/sideRETRO. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Antagonism between the RNA-binding protein Musashi1 and miR-137 and its potential impact on neurogenesis and glioblastoma development.

RNA-binding proteins (RBPs) and miRNAs are critical gene expression regulators that interact with one another in cooperative and antagonistic fashions. We identified Musashi1 (Msi1) and miR-137 as regulators of a molecular switch between self-renewal and differentiation. Msi1 and miR-137 have opposite expression patterns and functions, and Msi1 is repressed by miR-137. Msi1 is a stem-cell protein implicated in self-renewal while miR-137 functions as a proneuronal differentiation miRNA. In gliomas, miR-137 functions as a tumor suppressor while Msi1 is a prooncogenic factor. We suggest that the balance between Msi1 and miR-137 is a key determinant in cell fate decisions and disruption of this balance could contribute to neurodegenerative diseases and glioma development. Genomic analyses revealed that Msi1 and miR-137 share 141 target genes associated with differentiation, development, and morphogenesis. Initial results pointed out that these two regulators have an opposite impact on the expression of their target genes. Therefore, we propose an antagonistic model in which this network of shared targets could be either repressed by miR-137 or activated by Msi1, leading to different outcomes (self-renewal, proliferation, tumorigenesis).

Gene copy-number polymorphism caused by retrotransposition in humans.

The era of whole-genome sequencing has revealed that gene copy-number changes caused by duplication and deletion events have important evolutionary, functional, and phenotypic consequences. Recent studies have therefore focused on revealing the extent of variation in copy-number within natural populations of humans and other species. These studies have found a large number of copy-number variants (CNVs) in humans, many of which have been shown to have clinical or evolutionary importance. For the most part, these studies have failed to detect an important class of gene copy-number polymorphism: gene duplications caused by retrotransposition, which result in a new intron-less copy of the parental gene being inserted into a random location in the genome. Here we describe a computational approach leveraging next-generation sequence data to detect gene copy-number variants caused by retrotransposition (retroCNVs), and we report the first genome-wide analysis of these variants in humans. We find that retroCNVs account for a substantial fraction of gene copy-number differences between any two individuals. Moreover, we show that these variants may often result in expressed chimeric transcripts, underscoring their potential for the evolution of novel gene functions. By locating the insertion sites of these duplicates, we are able to show that retroCNVs have had an important role in recent human adaptation, and we also uncover evidence that positive selection may currently be driving multiple retroCNVs toward fixation. Together these findings imply that retroCNVs are an especially important class of polymorphism, and that future studies of copy-number variation should search for these variants in order to illuminate their potential evolutionary and functional relevance.

ICRmax: an optimized approach to detect tumor-specific interchromosomal rearrangements for clinical application.

Somatically acquired chromosomal rearrangements occur at early stages during tumorigenesis and can be used to indirectly detect tumor cells, serving as highly sensitive and tumor-specific biomarkers. Advances in high-throughput sequencing have allowed the genome-wide identification of patient-specific chromosomal rearrangements to be used as personalized biomarkers to efficiently assess response to treatment, detect residual disease and monitor disease recurrence. However, sequencing and data processing costs still represent major obstacles for the widespread application of personalized biomarkers in oncology. We developed a computational pipeline (ICRmax) for the cost-effective identification of a minimal set of tumor-specific interchromosomal rearrangements (ICRs). We examined ICRmax performance on sequencing data from rectal tumors and simulated data achieving an average accuracy of 68% for ICR identification. ICRmax identifies ICRs from low-coverage sequenced tumors, eliminates the need to sequence a matched normal tissue and significantly reduces the costs that limit the utilization of personalized biomarkers in the clinical setting.

RCPedia: a database of retrocopied genes.

MOTIVATION: Retrocopies are copies of mature RNAs that are usually devoid of regulatory sequences and introns. They have routinely been classified as processed pseudo-genes with little or no biological relevance. However, recent findings have revealed functional roles for retrocopies, as well as their high frequency in some organisms, such as primates. Despite their increasing importance, there is no user-friendly and publicly available resource for the study of retrocopies. RESULTS: Here, we present RCPedia, an integrative and user-friendly database designed for the study of retrocopied genes. RCPedia contains a complete catalogue of the retrocopies that are known to be present in human and five other primate genomes, their genomic context, inter-species conservation and gene expression data. RCPedia also offers a streamlined data representation and an efficient query system. AVAILABILITY AND IMPLEMENTATION: RCPedia is available at http://www.bioinfo.mochsl.org.br/rcpedia.

Impact of the thyroid hormone T3 and its nuclear receptor TRα1 on colon cancer stem cell phenotypes and response to chemotherapies.

Colorectal cancers (CRCs) are highly heterogeneous and show a hierarchical organization, with cancer stem cells (CSCs) responsible for tumor development, maintenance, and drug resistance. Our previous studies showed the importance of thyroid hormone-dependent signaling on intestinal tumor development and progression through action on stem cells. These results have a translational value, given that the thyroid hormone nuclear receptor TRα1 is upregulated in human CRCs, including in the molecular subtypes associated with CSC features. We used an established spheroid model generated from the human colon adenocarcinoma cell line Caco2 to study the effects of T3 and TRα1 on spheroid formation, growth, and response to conventional chemotherapies. Our results show that T3 treatment and/or increased TRα1 expression in spheroids impaired the response to FOLFIRI and conferred a survival advantage. This was achieved by stimulating drug detoxification pathways and increasing ALDH1A1-expressing cells, including CSCs, within spheroids. These results suggest that clinical evaluation of the thyroid axis and assessing TRα1 levels in CRCs could help to select optimal therapeutic regimens for patients with CRC. Proposed mechanism of action of T3/TRα1 in colon cancer spheroids. In the control condition, TRα1 participates in maintaining homeostatic cell conditions. The presence of T3 in the culture medium activates TRα1 action on target genes, including the drug efflux pumps ABCG2 and ABCB1. In the case of chemotherapy FOLFIRI, the increased expression of ABC transcripts and proteins induced by T3 treatment is responsible for the augmented efflux of 5-FU and Irinotecan from the cancer cells. Taken together, these mechanisms contribute to the decreased efficacy of the chemotherapy and allow cells to escape the treatment. Created with BioRender.com .

SERBP1 interacts with PARP1 and is present in PARylation-dependent protein complexes regulating splicing, cell division, and ribosome biogenesis.

RNA binding proteins (RBPs) containing intrinsically disordered regions (IDRs) are present in diverse molecular complexes where they function as dynamic regulators. Their characteristics promote liquid-liquid phase separation (LLPS) and the formation of membraneless organelles such as stress granules and nucleoli. IDR-RBPs are particularly relevant in the nervous system and their dysfunction is associated with neurodegenerative diseases and brain tumor development. SERBP1 is a unique member of this group, being mostly disordered and lacking canonical RNA-binding domains. Using a proteomics approach followed by functional analysis, we defined SERBP1's interactome. We uncovered novel SERBP1 roles in splicing, cell division, and ribosomal biogenesis and showed its participation in pathological stress granules and Tau aggregates in Alzheimer's disease brains. SERBP1 preferentially interacts with other G-quadruplex (G4) binders, implicated in different stages of gene expression, suggesting that G4 binding is a critical component of SERBP1 function in different settings. Similarly, we identified important associations between SERBP1 and PARP1/polyADP-ribosylation (PARylation). SERBP1 interacts with PARP1 and its associated factors and influences PARylation. Moreover, protein complexes in which SERBP1 participates contain mostly PARylated proteins and PAR binders. Based on these results, we propose a feedback regulatory model in which SERBP1 influences PARP1 function and PARylation, while PARylation modulates SERBP1 functions and participation in regulatory complexes.

Alternative splicing and genetic diversity: silencers are more frequently modified by SNVs associated with alternative exon/intron borders.

With the availability of a large amount of genomic data it is expected that the influence of single nucleotide variations (SNVs) in many biological phenomena will be elucidated. Here, we approached the problem of how SNVs affect alternative splicing. First, we observed that SNVs and exonic splicing regulators (ESRs) independently show a biased distribution in alternative exons. More importantly, SNVs map more frequently in ESRs located in alternative exons than in ESRs located in constitutive exons. By looking at SNVs associated with alternative exon/intron borders (by their common presence in the same cDNA molecule), we observed that a specific type of ESR, the exonic splicing silencers (ESSs), are more frequently modified by SNVs. Our results establish a clear association between genetic diversity and alternative splicing involving ESSs.

Distinct patterns of somatic alterations in a lymphoblastoid and a tumor genome derived from the same individual.

Although patterns of somatic alterations have been reported for tumor genomes, little is known on how they compare with alterations present in non-tumor genomes. A comparison of the two would be crucial to better characterize the genetic alterations driving tumorigenesis. We sequenced the genomes of a lymphoblastoid (HCC1954BL) and a breast tumor (HCC1954) cell line derived from the same patient and compared the somatic alterations present in both. The lymphoblastoid genome presents a comparable number and similar spectrum of nucleotide substitutions to that found in the tumor genome. However, a significant difference in the ratio of non-synonymous to synonymous substitutions was observed between both genomes (P = 0.031). Protein-protein interaction analysis revealed that mutations in the tumor genome preferentially affect hub-genes (P = 0.0017) and are co-selected to present synergistic functions (P < 0.0001). KEGG analysis showed that in the tumor genome most mutated genes were organized into signaling pathways related to tumorigenesis. No such organization or synergy was observed in the lymphoblastoid genome. Our results indicate that endogenous mutagens and replication errors can generate the overall number of mutations required to drive tumorigenesis and that it is the combination rather than the frequency of mutations that is crucial to complete tumorigenic transformation.

Retro-miRs: novel and functional miRNAs originating from mRNA retrotransposition.

BACKGROUND: Reverse-transcribed gene copies (retrocopies) have emerged as major sources of evolutionary novelty. MicroRNAs (miRNAs) are small and highly conserved RNA molecules that serve as key post-transcriptional regulators of gene expression. The origin and subsequent evolution of miRNAs have been addressed but not fully elucidated. RESULTS: In this study, we performed a comprehensive investigation of miRNA origination through retroduplicated mRNA sequences (retro-miRs). We identified 17 retro-miRs that emerged from the mRNA retrocopies. Four of these retro-miRs had de novo origins within retrocopied sequences, while 13 retro-miRNAs were located within exon regions and duplicated along with their host mRNAs. We found that retro-miRs were primate-specific, including five retro-miRs conserved among all primates and two human-specific retro-miRs. All retro-miRs were expressed, with predicted and experimentally validated target genes except miR-10527. Notably, the target genes of retro-miRs are involved in key biological processes such as metabolic processes, cell signaling, and regulation of neurotransmitters in the central nervous system. Additionally, we found that these retro-miRs play a potential oncogenic role in cancer by targeting key cancer genes and are overexpressed in several cancer types, including liver hepatocellular carcinoma and stomach adenocarcinoma. CONCLUSIONS: Our findings demonstrated that mRNA retrotransposition is a key mechanism for the generation of novel miRNAs (retro-miRs) in primates. These retro-miRs are expressed, conserved, have target genes with important cellular functions, and play important roles in cancer.

ELF4 is a critical component of a miRNA-transcription factor network and is a bridge regulator of glioblastoma receptor signaling and lipid dynamics.

BACKGROUND: The loss of neurogenic tumor suppressor microRNAs miR-124, miR-128, and miR-137 is associated with glioblastoma's undifferentiated state. Most of their impact comes via the repression of a network of oncogenic transcription factors. We conducted a high-throughput functional siRNA screen in glioblastoma cells and identify E74 like ETS transcription factor 4 (ELF4) as the leading contributor to oncogenic phenotypes. METHODS: In vitro and in vivo assays were used to assess ELF4 impact on cancer phenotypes. We characterized ELF4's mechanism of action via genomic and lipidomic analyses. A MAPK reporter assay verified ELF4's impact on MAPK signaling, and qRT-PCR and western blotting were used to corroborate ELF4 regulatory role on most relevant target genes. RESULTS: ELF4 knockdown resulted in significant proliferation delay and apoptosis in GBM cells and long-term growth delay and morphological changes in glioma stem cells (GSCs). Transcriptomic analyses revealed that ELF4 controls two interlinked pathways: 1) Receptor tyrosine kinase signaling and 2) Lipid dynamics. ELF4 modulation directly affected receptor tyrosine kinase (RTK) signaling, as mitogen-activated protein kinase (MAPK) activity was dependent upon ELF4 levels. Furthermore, shotgun lipidomics revealed that ELF4 depletion disrupted several phospholipid classes, highlighting ELF4's importance in lipid homeostasis. CONCLUSIONS: We found that ELF4 is critical for the GBM cell identity by controlling genes of two dependent pathways: RTK signaling (SRC, PTK2B, and TNK2) and lipid dynamics (LRP1, APOE, ABCA7, PLA2G6, and PITPNM2). Our data suggest that targeting these two pathways simultaneously may be therapeutically beneficial to GBM patients.