GmNAP1 is essential for trichome and leaf epidermal cell development in soybean.

KEY MESSAGE: Map-based cloning revealed that two novel soybean distorted trichome mutants were due to loss function of GmNAP1 gene, which affected the trichome morphology and pavement cell ploidy by regulating actin filament assembly. Trichomes increase both biotic and abiotic stress resistance in soybean. In this study, Gmdtm1-1 and Gmdtm1-2 mutants with shorter trichomes and bigger epidermal pavement cells were isolated from an ethyl methylsulfonate mutagenized population. Both of them had reduced plant height and smaller seeds. Map-based cloning and bulked segregant analysis identified that a G-A transition at the 3' boundary of the sixth intron of Glyma.20G019300 in the Gmdtm1-1 mutant and another G-A transition mutation at the 5' boundary of the fourteenth intron of Glyma.20G019300 in Gmdtm1-2; these mutations disrupted spliceosome recognition sites creating truncated proteins. Glyma.20G019300 encodes a Glycine max NCK-associated protein 1 homolog (GmNAP1) in soybean. Further analysis revealed that the GmNAP1 involved in actin filament assembling and genetic information processing pathways during trichome and pavement cell development. This study shows that GmNAP1 plays an important role in soybean growth and development and agronomic traits.

Circular HER2 RNA positive triple negative breast cancer is sensitive to Pertuzumab.

BACKGROUND: Triple negative breast cancer (TNBC) remains the most challenging breast cancer subtype so far. Specific therapeutic approaches have rarely achieved clinical improvements in treatment of TNBC patients and effective molecular biomarkers are largely unknown. METHODS: We used paired TNBC samples and high throughput RNA sequencing to identify differentially expressed circRNAs. Sucrose gradient polysome fractionation assay, antibody and Mass spectra were used to validate active circRNA translation. The novel protein function was validated in vitro and in vivo by gain or loss of function assays. Mechanistic results were concluded by immunoprecipitation analyses and kinase activity assay. RESULTS: Circular HER2 RNA (circ-HER2) encoded a novel protein, HER2-103. Unexpectedly, while HER2 mRNA and protein were barely detected, circ-HER2/HER2-103 was expressed in ~ 30% TNBC clinical samples. Circ-HER2/HER2-103 positive TNBC patients harbored worse overall prognosis than circ-HER2/HER2-103 negative patients. Knockdown circ-HER2 inhibited TNBC cells proliferation, invasion and tumorigenesis in vitro and in vivo, suggesting the critical role of circ-HER2/HER2-103 in TNBC tumorigenicity. Mechanistically, HER2-103 promoted homo/hetero dimerization of epidermal growth factor receptor (EGFR)/HER3, sustained AKT phosphorylation and downstream malignant phenotypes. Furthermore, HER2-103 shared most of the same amino acid sequences as HER2 CR1 domain which could be antagonized by Pertuzumab, a clinical used HER2 antibody. Pertuzumab markedly attenuated in vivo tumorigenicity of circ-HER2/HER2-103 expressing TNBC cells but showed no effects in circ-HER2/HER2-103 negative TNBC cells. CONCLUSION: Our results not only demonstrated that certain TNBCs were not truly 'HER2 negative' but also highlighted the clinical implications of Pertuzumab in circ-HER2/HER2-103 expressing TNBC patients.

Differential expression networks and inheritance patterns of long non-coding RNAs in castor bean seeds.

Long non-coding RNAs (lncRNAs) serve as versatile regulators of plant growth and development. The potential functions and inheritance patterns of lncRNAs, as well as the epigenetic regulation of lncRNA itself, remain largely uncharacterized in plant seeds, especially in the persistent endosperm of the dicotyledons. In this study, we investigated diverse RNA-seq data and catalogued 5356 lncRNAs in castor bean seeds. A small fraction of lncRNAs were transcribed from the same direction as the promoters of protein-coding genes (PCgenes) and exhibited strongly coordinated expression with the nearby PCgene. Co-expression analysis with weighted gene co-expression network analysis (WGCNA) showed these lncRNAs to be involved in differential transcription networks between the embryo and endosperm in the early developing seed. Genomic DNA methylation analyses revealed that the expression level of lncRNAs was tightly linked to DNA methylation and that endosperm hypomethylation could promote the expression of linked lncRNAs. Intriguingly, upon hybridization, most lncRNAs with divergent genome sequences between two parents could be reconciled and were expressed according to their parental genome contribution; however, some deviation in the expression of allelic lncRNAs was observed and found to be partially dependent on parental effects. In triploid endosperm, the expression of most lncRNAs was not dosage sensitive, as only 20 lncRNAs had balanced dosage. Our findings not only demonstrate that lncRNAs play potential roles in regulating the development of castor bean endosperm and embryo, but also provide novel insights into the parental effects, allelic expression and epigenetic regulation of lncRNAs in dicotyledonous seeds.

Correction to: A novel tumor suppressor protein encoded by circular AKT3 RNA inhibits glioblastoma tumorigenicity by competing with active phosphoinositide-dependent Kinase-1.

In the published article [1], an error was noticed in Fig. 6B. The western blot results were reversed between the overexpression group and the knockdown group of circ-AKT3. The corrected and updated Fig. 6 is provided below. This error does not affect the findings or conclusions of the article.

A novel tumor suppressor protein encoded by circular AKT3 RNA inhibits glioblastoma tumorigenicity by competing with active phosphoinositide-dependent Kinase-1.

BACKGROUND: The RTK/PI3K/AKT pathway plays key roles in the development and progression of many cancers, including GBM. As a regulatory molecule and a potential drug target, the oncogenic role of AKT has been substantially studied. Three isoforms of AKT have been identified, including AKT1, AKT2 and AKT3, but their individual functions in GBM remain controversial. Moreover, it is not known if there are more AKT alternative splicing variants. METHODS: High-throughput RNA sequencing and quantitative reverse transcription-PCR were used to identify the differentially expressed circRNAs in GBM samples and in paired normal tissues. High throughput RNA sequencing was used to identify circ-AKT3 regulated signaling pathways. Mass spectrometry, western blotting and immunofluorescence staining analyses were used to validate AKT3-174aa expression. The tumor suppressive role of AKT3-174aa was validated in vitro and in vivo. The competing interaction between AKT3-174aa and p-PDK1 was investigated by mass spectrometry and immunoprecipitation analyses. RESULTS: Circ-AKT3 is a previously uncharacterized AKT transcript variant. Circ-AKT3 is expressed at low levels in GBM tissues compared with the expression in paired adjacent normal brain tissues. Circ-AKT3 encodes a 174 amino acid (aa) novel protein, which we named AKT3-174aa, by utilizing overlapping start-stop codons. AKT3-174aa overexpression decreased the cell proliferation, radiation resistance and in vivo tumorigenicity of GBM cells, while the knockdown of circ-AKT3 enhanced the malignant phenotypes of astrocytoma cells. AKT3-174aa competitively interacts with phosphorylated PDK1, reduces AKT-thr308 phosphorylation, and plays a negative regulatory role in modulating the PI3K/AKT signal intensity. CONCLUSIONS: Our data indicate that the impaired circRNA expression of the AKT3 gene contributes to GBM tumorigenesis, and our data corroborate the hypothesis that restoring AKT3-174aa while inhibiting activated AKT may provide more benefits for certain GBM patients.

Chromosome level high-density integrated genetic maps improve the Pyrus bretschneideri 'DangshanSuli' v1.0 genome.

BACKGROUND: Chromosomal level reference genomes provide a crucial foundation for genomics research such as genome-wide association studies (GWAS) and whole genome selection. The chromosomal-level sequences of both the European (Pyrus communis) and Chinese (P. bretschneideri) pear genomes have not been published in public databases so far. RESULTS: To anchor the scaffolds of P. bretschneideri 'DangshanSuli' (DS) v1.0 genome into pseudo-chromosomes, two genetic maps (MH and YM maps) were constructed using half sibling populations of Chinese pear crosses, 'Mantianhong' (MTH) × 'Hongxiangsu' (HXS) and 'Yuluxiang' (YLX) × MTH, from 345 and 162 seedlings, respectively, which were prepared for SNP discovery using genotyping-by-sequencing (GBS) technology. The MH and YM maps, each with 17 linkage groups (LGs), were constructed from 2606 and 2489 SNP markers and spanned 1847 and 1668 cM, respectively, with average marker intervals of 0.7. The two maps were further merged with a previously published genetic map (BD) based on the cross 'Bayuehong' (BYH) × 'Dangshansuli' (DS) to build a new integrated MH-YM-BD map. By using 7757 markers located on the integrated MH-YM-BD map, 898 scaffolds (400.57 Mb) of the DS v1.0 assembly were successfully anchored into 17 pseudo-chromosomes, accounting for 78.8% of the assembled genome size. About 88.31% of them (793 scaffolds) were directionally anchored with two or more markers on the pseudo-chromosomes. Furthermore, the errors in each pseudo-chromosome (especially 1, 5, 7 and 11) were manually corrected and pseudo-chromosomes 1, 5 and 7 were extended by adding 19, 12 and 14 scaffolds respectively in the newly constructed DS v1.1 genome. Synteny analyses revealed that the DS v1.1 genome had high collinearity with the apple genome, and the homologous fragments between pseudo-chromosomes were similar to those found in previous studies. Moreover, the red-skin trait of Asian pear was mapped to an identical locus as identified previously. CONCLUSIONS: The accuracy of DS v1.1 genome was improved by using larger mapping populations and merged genetic map. With more than 400 MB anchored to 17 pseudo-chromosomes, the new DS v1.1 genome provides a critical tool that is essential for studies of pear genetics, genomics and molecular breeding.

Whole-genome sequencing of cultivated and wild peppers provides insights into Capsicum domestication and specialization.

As an economic crop, pepper satisfies people's spicy taste and has medicinal uses worldwide. To gain a better understanding of Capsicum evolution, domestication, and specialization, we present here the genome sequence of the cultivated pepper Zunla-1 (C. annuum L.) and its wild progenitor Chiltepin (C. annuum var. glabriusculum). We estimate that the pepper genome expanded ∼0.3 Mya (with respect to the genome of other Solanaceae) by a rapid amplification of retrotransposons elements, resulting in a genome comprised of ∼81% repetitive sequences. Approximately 79% of 3.48-Gb scaffolds containing 34,476 protein-coding genes were anchored to chromosomes by a high-density genetic map. Comparison of cultivated and wild pepper genomes with 20 resequencing accessions revealed molecular footprints of artificial selection, providing us with a list of candidate domestication genes. We also found that dosage compensation effect of tandem duplication genes probably contributed to the pungent diversification in pepper. The Capsicum reference genome provides crucial information for the study of not only the evolution of the pepper genome but also, the Solanaceae family, and it will facilitate the establishment of more effective pepper breeding programs.

Genome-wide comparative analysis of digital gene expression tag profiles during maize ear development.

The present study profiled and analyzed gene expression of the maize ear at four key developmental stages. Based on genome-wide profile analysis, we detected differential mRNA of maize genes. Some of the differentially expressed genes (DEGs) were predicted to be potential candidates of maize ear development. Several well-known genes were found with reported mutant analyses, such as, compact plant2 (ct2), zea AGAMOUS homolog1 (zag1), bearded ear (bde), and silky1 (si1). MicroRNAs such as microRNA156 were predicted to target genes involved in maize ear development. Antisense transcripts were widespread throughout all the four stages, and are suspected to play important roles in maize ear development. Thus, identification and characterization of important genes and regulators at all the four developmental stages will contribute to an improved understanding of the molecular mechanisms responsible for maize ear development.

An ultra-high-density map as a community resource for discerning the genetic basis of quantitative traits in maize.

BACKGROUND: To safeguard the food supply for the growing human population, it is important to understand and exploit the genetic basis of quantitative traits. Next-generation sequencing technology performs advantageously and effectively in genetic mapping and genome analysis of diverse genetic resources. Hence, we combined re-sequencing technology and a bin map strategy to construct an ultra-high-density bin map with thousands of bin markers to precisely map a quantitative trait locus. RESULTS: In this study, we generated a linkage map containing 1,151,856 high quality SNPs between Mo17 and B73, which were verified in the maize intermated B73 × Mo17 (IBM) Syn10 population. This resource is an excellent complement to existing maize genetic maps available in an online database  (iPlant, http://data.maizecode.org/maize/qtl/syn10/ ). Moreover, in this population combined with the IBM Syn4 RIL population, we detected 135 QTLs for flowering time and plant height traits across the two populations. Eighteen known functional genes and twenty-five candidate genes for flowering time and plant height trait were fine-mapped into a 2.21-4.96 Mb interval. Map expansion and segregation distortion were also analyzed, and evidence for inadvertent selection of early flowering time in the process of mapping population development was observed. Furthermore, an updated integrated map with 1,151,856 high-quality SNPs, 2,916 traditional markers and 6,618 bin markers was constructed. The data were deposited into the iPlant Discovery Environment (DE), which provides a fundamental resource of genetic data for the maize genetic research community. CONCLUSIONS: Our findings provide basic essential genetic data for the maize genetic research community. An updated IBM Syn10 population and a reliable, verified high-quality SNP set between Mo17 and B73 will aid in future molecular breeding efforts.

Identification of miRNAs and their target genes in developing maize ears by combined small RNA and degradome sequencing.

BACKGROUND: In plants, microRNAs (miRNAs) are endogenous ~22 nt RNAs that play important regulatory roles in many aspects of plant biology, including metabolism, hormone response, epigenetic control of transposable elements, and stress response. Extensive studies of miRNAs have been performed in model plants such as rice and Arabidopsis thaliana. In maize, most miRNAs and their target genes were analyzed and identified by clearly different treatments, such as response to low nitrate, salt and drought stress. However, little is known about miRNAs involved in maize ear development. The objective of this study is to identify conserved and novel miRNAs and their target genes by combined small RNA and degradome sequencing at four inflorescence developmental stages. RESULTS: We used deep-sequencing, miRNA microarray assays and computational methods to identify, profile, and describe conserved and non-conserved miRNAs at four ear developmental stages, which resulted in identification of 22 conserved and 21-maize-specific miRNA families together with their corresponding miRNA*. Comparison of miRNA expression in these developmental stages revealed 18 differentially expressed miRNA families. Finally, a total of 141 genes (251 transcripts) targeted by 102 small RNAs including 98 miRNAs and 4 ta-siRNAs were identified by genomic-scale high-throughput sequencing of miRNA cleaved mRNAs. Moreover, the differentially expressed miRNAs-mediated pathways that regulate the development of ears were discussed. CONCLUSIONS: This study confirmed 22 conserved miRNA families and discovered 26 novel miRNAs in maize. Moreover, we identified 141 target genes of known and new miRNAs and ta-siRNAs. Of these, 72 genes (117 transcripts) targeted by 62 differentially expressed miRNAs may attribute to the development of maize ears. Identification and characterization of these important classes of regulatory genes in maize may improve our understanding of molecular mechanisms controlling ear development.

Circular RNA encoded MET variant promotes glioblastoma tumorigenesis.

Activated by its single ligand, hepatocyte growth factor (HGF), the receptor tyrosine kinase MET is pivotal in promoting glioblastoma (GBM) stem cell self-renewal, invasiveness and tumorigenicity. Nevertheless, HGF/MET-targeted therapy has shown limited clinical benefits in GBM patients, suggesting hidden mechanisms of MET signalling in GBM. Here, we show that circular MET RNA (circMET) encodes a 404-amino-acid MET variant (MET404) facilitated by the N6-methyladenosine (m6A) reader YTHDF2. Genetic ablation of circMET inhibits MET404 expression in mice and attenuates MET signalling. Conversely, MET404 knock-in (KI) plus P53 knock-out (KO) in mouse astrocytes initiates GBM tumorigenesis and shortens the overall survival. MET404 directly interacts with the MET ß subunit and forms a constitutively activated MET receptor whose activity does not require HGF stimulation. High MET404 expression predicts poor prognosis in GBM patients, indicating its clinical relevance. Targeting MET404 through a neutralizing antibody or genetic ablation reduces GBM tumorigenicity in vitro and in vivo, and combinatorial benefits are obtained with the addition of a traditional MET inhibitor. Overall, we identify a MET variant that promotes GBM tumorigenicity, offering a potential therapeutic strategy for GBM patients, especially those with MET hyperactivation.

A novel protein RASON encoded by a lncRNA controls oncogenic RAS signaling in KRAS mutant cancers.

Mutations of the RAS oncogene are found in around 30% of all human cancers yet direct targeting of RAS is still considered clinically impractical except for the KRASG12C mutant. Here we report that RAS-ON (RASON), a novel protein encoded by the long intergenic non-protein coding RNA 00673 (LINC00673), is a positive regulator of oncogenic RAS signaling. RASON is aberrantly overexpressed in pancreatic ductal adenocarcinoma (PDAC) patients, and it promotes proliferation of human PDAC cell lines in vitro and tumor growth in vivo. CRISPR/Cas9-mediated knockout of Rason in mouse embryonic fibroblasts inhibits KRAS-mediated tumor transformation. Genetic deletion of Rason abolishes oncogenic KRAS-driven pancreatic and lung cancer tumorigenesis in LSL-KrasG12D; Trp53R172H/+ mice. Mechanistically, RASON directly binds to KRASG12D/V and inhibits both intrinsic and GTPase activating protein (GAP)-mediated GTP hydrolysis, thus sustaining KRASG12D/V in the GTP-bound hyperactive state. Therapeutically, deprivation of RASON sensitizes KRAS mutant pancreatic cancer cells and patient-derived organoids to EGFR inhibitors. Our findings identify RASON as a critical regulator of oncogenic KRAS signaling and a promising therapeutic target for KRAS mutant cancers.

Circular EZH2-encoded EZH2-92aa mediates immune evasion in glioblastoma via inhibition of surface NKG2D ligands.

Glioblastoma (GBM) is a highly aggressive primary brain tumour and is resistant to nearly all available treatments, including natural killer (NK) cell immunotherapy. However, the factors mediating NK cell evasion in GBM remain largely unclear. Here, we report that EZH2-92aa, a protein encoded by circular EZH2, is overexpressed in GBM and induces the immune evasion of GBM stem cells (GSCs) from NK cells. Positively regulated by DEAD-box helicase 3 (DDX3), EZH2-92aa directly binds the major histocompatibility complex class I polypeptide-related sequence A/B (MICA/B) promoters and represses their transcription; it also indirectly represses UL16-binding protein (ULBP) transcription by stabilizing EZH2. The downregulation of NK group 2D ligands (NKG2DLs, including MICA/B and ULBPs) in GSCs mediates NK cell resistance. Moreover, stable EZH2-92aa knockdown enhances NK cell-mediated GSC eradication in vitro and in vivo and synergizes with anti-PD1 therapy. Our results highlight the immunosuppressive function of EZH2-92aa in inhibiting the NK cell response in GBM and the clinical potential of targeting EZH2-92aa for NK-cell-directed immune therapy.

Eight soybean reference genome resources from varying latitudes and agronomic traits.

Comparative analysis of multiple reference genomes representing diverse genetic backgrounds is critical for understanding the role of key alleles important in domestication and genetic breeding of important crops such as soybean. To enrich the genetic resources for soybean, we describe the generation, technical assessment, and preliminary genomic variation analysis of eight de novo reference-grade soybean genome assemblies from wild and cultivated accessions. These resources represent soybeans cultured at different latitudes and exhibiting different agronomical traits. Of these eight soybeans, five are from new accessions that have not been sequenced before. We demonstrate the usage of these genomes to identify small and large genomic variations affecting known genes as well as screening for genic PAV regions for identifying candidates for further functional studies.

A Robust and Rapid Candidate Gene Mapping Pipeline Based on M2 Populations.

The whole-genome sequencing-based bulked segregant analysis (WGS-BSA) has facilitated the mapping candidate causal variations for cloning target plant genes. Here, we report an improved WGS-BSA method termed as M2-seq to expedite the mapping candidate mutant loci by studying just M2 generation. It is an efficient mutant gene mapping tool, rapid, and comparable to the previously reported approaches, such as Mutmap and Mutmap+ that require studying M3 or advanced selfed generations. In M2-seq, background variations among the M2 populations can be removed efficiently without knowledge of the variations of the wild-type progenitor plant. Furthermore, the use of absolute delta single-nucleotide polymorphism (SNP) index values can effectively remove the background variation caused by repulsion phase linkages of adjacent mutant alleles; and thereby facilitating the identification of the causal mutation in target genes. Here, we demonstrated the application of M2-seq in successfully mapping the genomic regions harboring causal mutations for mutant phenotypes among 10 independent M2 populations of soybean. The mapping candidate mutant genes just in M2 generation with the aid of the M2-seq method should be particularly useful in expediting gene cloning especially among the plant species with long generation time.

Rolling-translated EGFR variants sustain EGFR signaling and promote glioblastoma tumorigenicity.

BACKGROUND: Aberrant epidermal growth factor receptor (EGFR) activation is observed in over 50% of cases of adult glioblastoma (GBM). Nevertheless, EGFR antibodies are ineffective in clinical GBM treatment, suggesting the existence of redundant EGFR activation mechanisms. Whether circular RNA (circRNA) encodes a protein involved in EGFR-driven GBM remains unclear. We reported an unexpected mechanism in which circular EGFR RNA (circ-EGFR) encodes a novel EGFR variant to sustained EGFR activation. METHOD: We used RNA-seq, Northern blot, and Sanger sequencing to confirm the existence of circ-EGFR. Antibodies and a liquid chromatograph tandem mass spectrometer were used to identify circ-EGFR protein products. Lentivirus-transfected stable cell lines were used to assess the biological functions of the novel protein in vitro and in vivo. Clinical implications of circ-EGFR were assessed using 97 pathologically diagnosed GBM patient samples. RESULTS: The infinite open reading frame (iORF) in circ-EGFR translated repeating amino acid sequences via rolling translation and programmed -1 ribosomal frameshifting (-1PRF) induced out-of-frame stop codon (OSC), forming a polymetric novel protein-complex, which we termed rolling-translated EGFR (rtEGFR). rtEGFR directly interacted with EGFR, maintained EGFR membrane localization and attenuated EGFR endocytosis and degradation. Importantly, circ-EGFR levels correlated with the EGFR signature and predicted the poor prognosis of GBM patients. Deprivation of rtEGFR in brain tumor-initiating cells (BTICs) attenuated tumorigenicity and enhanced the anti-GBM effect. CONCLUSION: Our findings identified the endogenous rolling-translated protein and provided strong clinical evidence that targeting rtEGFR could improve the efficiency of EGFR-targeting therapies in GBM.

An Upstream Open Reading Frame in Phosphatase and Tensin Homolog Encodes a Circuit Breaker of Lactate Metabolism.

The metabolic role of micropeptides generated from untranslated regions remains unclear. Here we describe MP31, a micropeptide encoded by the upstream open reading frame (uORF) of phosphatase and tensin homolog (PTEN) acting as a "circuit breaker" that limits lactate-pyruvate conversion in mitochondria by competing with mitochondrial lactate dehydrogenase (mLDH) for nicotinamide adenine dinucleotide (NAD+). Knocking out the MP31 homolog in mice enhanced global lactate metabolism, manifesting as accelerated oxidative phosphorylation (OXPHOS) and increased lactate consumption and production. Conditional knockout (cKO) of MP31 homolog in mouse astrocytes initiated gliomagenesis and shortened the overall survival of the animals, establishing a tumor-suppressing role for MP31. Recombinant MP31 administered intraperitoneally penetrated the blood-brain barrier and inhibited mice GBM xenografts without neurological toxicity, suggesting the clinical implication and application of this micropeptide. Our findings reveal a novel mode of MP31-orchestrated lactate metabolism reprogramming in glioblastoma.

A novel protein encoded by circular SMO RNA is essential for Hedgehog signaling activation and glioblastoma tumorigenicity.

BACKGROUND: Aberrant activation of the Hedgehog pathway drives tumorigenesis of many cancers, including glioblastoma. However, the sensitization mechanism of the G protein-coupled-like receptor smoothened (SMO), a key component of Hedgehog signaling, remains largely unknown. RESULTS: In this study, we describe a novel protein SMO-193a.a. that is essential for Hedgehog signaling activation in glioblastoma. Encoded by circular SMO (circ-SMO), SMO-193a.a. is required for sonic hedgehog (Shh) induced SMO activation, via interacting with SMO, enhancing SMO cholesterol modification, and releasing SMO from the inhibition of patched transmembrane receptors. Deprivation of SMO-193a.a. in brain cancer stem cells attenuates Hedgehog signaling intensity and suppresses self-renewal, proliferation in vitro, and tumorigenicity in vivo. Moreover, circ-SMO/SMO-193a.a. is positively regulated by FUS, a direct transcriptional target of Gli1. Shh/Gli1/FUS/SMO-193a.a. form a positive feedback loop to sustain Hedgehog signaling activation in glioblastoma. Clinically, SMO-193a.a. is more specifically expressed in glioblastoma than SMO and is relevant to Gli1 expression. Higher expression of SMO-193a.a. predicts worse overall survival of glioblastoma patients, indicating its prognostic value. CONCLUSIONS: Our study reveals that SMO-193a.a., a novel protein encoded by circular SMO, is critical for Hedgehog signaling, drives glioblastoma tumorigenesis and is a novel target for glioblastoma treatment.

Circular RNA-encoded oncogenic E-cadherin variant promotes glioblastoma tumorigenicity through activation of EGFR-STAT3 signalling.

Activated EGFR signalling drives tumorigenicity in 50% of glioblastoma (GBM). However, EGFR-targeting therapy has proven ineffective in treating patients with GBM, indicating that there is redundant EGFR activation. Circular RNAs are covalently closed RNA transcripts that are involved in various physiological and pathological processes. Herein, we report an additional activation mechanism of EGFR signalling in GBM by an undescribed secretory E-cadherin protein variant (C-E-Cad) encoded by a circular E-cadherin (circ-E-Cad) RNA through multiple-round open reading frame translation. C-E-Cad is overexpressed in GBM and promotes glioma stem cell tumorigenicity. C-E-Cad activates EGFR independent of EGF through association with the EGFR CR2 domain using a unique 14-amino-acid carboxy terminus, thereby maintaining glioma stem cell tumorigenicity. Notably, inhibition of C-E-Cad markedly enhances the antitumour activity of therapeutic anti-EGFR strategies in GBM. Our results uncover a critical role of C-E-Cad in stimulating EGFR signalling and provide a promising approach for treating EGFR-driven GBM.