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.

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.

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.

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.

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.

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.

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.

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.

Evolution and Expression Divergence of the CYP78A Subfamily Genes in Soybean.

Gene expression divergence is an important evolutionary driving force for the retention of duplicate genes. In this study, we identified three CYP78A subfamily genes in soybean, GmCYP78A70, GmCYP78A57 and GmCYP78A72, which experienced different duplication events. GmCYP78A70 was mainly expressed in leaf tissue and the vegetative phase, whereas GmCYP78A57 was mainly expressed in floral tissue and seed, i.e., the reproductive phase. Expression of GmCYP78A72 could be detected in all the tissues and phases mentioned above. The expression levels of GmCYP78A70 and GmCYP78A57 in different soybean cultivars showed positive correlations with leaf size and 100-seed weight, respectively. The population genetics analysis indicated that the three genes had experienced different selective pressures during domestication and improved breeding of soybean. Deciphering the function of this subfamily of genes may well prove useful to breeders for improving soybean's agronomic traits.

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.

A peptide encoded by circular form of LINC-PINT suppresses oncogenic transcriptional elongation in glioblastoma.

Circular RNAs (circRNAs) are a large class of transcripts in the mammalian genome. Although the translation of circRNAs was reported, additional coding circRNAs and the functions of their translated products remain elusive. Here, we demonstrate that an endogenous circRNA generated from a long noncoding RNA encodes regulatory peptides. Through ribosome nascent-chain complex-bound RNA sequencing (RNC-seq), we discover several peptides potentially encoded by circRNAs. We identify an 87-amino-acid peptide encoded by the circular form of the long intergenic non-protein-coding RNA p53-induced transcript (LINC-PINT) that suppresses glioblastoma cell proliferation in vitro and in vivo. This peptide directly interacts with polymerase associated factor complex (PAF1c) and inhibits the transcriptional elongation of multiple oncogenes. The expression of this peptide and its corresponding circRNA are decreased in glioblastoma compared with the levels in normal tissues. Our results establish the existence of peptides encoded by circRNAs and demonstrate their potential functions in glioblastoma tumorigenesis.

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.

A novel protein encoded by the circular form of the SHPRH gene suppresses glioma tumorigenesis.

Circular RNAs (circRNAs) are recognized as functional non-coding transcripts in eukaryotic cells. Recent evidence has indicated that even though circRNAs are generally expressed at low levels, they may be involved in many physiological or pathological processes, such as gene regulation, tissue development and carcinogenesis. Although the 'microRNA sponge' function is well characterized, most circRNAs do not contain perfect trapping sites for microRNAs, which suggests the possibility that circRNAs have functions that have not yet been defined. In this study, we show that a circRNA containing an open reading frame (ORF) driven by the internal ribosome entry site (IRES) can translate a functional protein. The circular form of the SNF2 histone linker PHD RING helicase (SHPRH) gene encodes a novel protein that we termed SHPRH-146aa. Circular SHPRH (circ-SHPRH) uses overlapping genetic codes to generate a 'UGA' stop codon, which results in the translation of the 17 kDa SHPRH-146aa. Both circ-SHPRH and SHPRH-146aa are abundantly expressed in normal human brains and are down-regulated in glioblastoma. The overexpression of SHPRH-146aa in U251 and U373 glioblastoma cells reduces their malignant behavior and tumorigenicity in vitro and in vivo. Mechanistically, SHPRH-146aa protects full-length SHPRH from degradation by the ubiquitin proteasome. Stabilized SHPRH sequentially ubiquitinates proliferating cell nuclear antigen (PCNA) as an E3 ligase, leading to inhibited cell proliferation and tumorigenicity. Our findings provide a novel perspective regarding circRNA function in physiological and pathological processes. Specifically, SHPRH-146aa generated from overlapping genetic codes of circ-SHPRH is a tumor suppressor in human glioblastoma.