IFITM3 functions as a PIP3 scaffold to amplify PI3K signalling in B cells.

Interferon-induced transmembrane protein 3 (IFITM3) has previously been identified as an endosomal protein that blocks viral infection1-3. Here we studied clinical cohorts of patients with B cell leukaemia and lymphoma, and identified IFITM3 as a strong predictor of poor outcome. In normal resting B cells, IFITM3 was minimally expressed and mainly localized in endosomes. However, engagement of the B cell receptor (BCR) induced both expression of IFITM3 and phosphorylation of this protein at Tyr20, which resulted in the accumulation of IFITM3 at the cell surface. In B cell leukaemia, oncogenic kinases phosphorylate IFITM3 at Tyr20, which causes constitutive localization of this protein at the plasma membrane. In a mouse model, Ifitm3-/- naive B cells developed in normal numbers; however, the formation of germinal centres and the production of antigen-specific antibodies were compromised. Oncogenes that induce the development of leukaemia and lymphoma did not transform Ifitm3-/- B cells. Conversely, the phosphomimetic IFITM3(Y20E) mutant induced oncogenic PI3K signalling and initiated the transformation of premalignant B cells. Mechanistic experiments revealed that IFITM3 functions as a PIP3 scaffold and central amplifier of PI3K signalling. The amplification of PI3K signals depends on IFITM3 using two lysine residues (Lys83 and Lys104) in its conserved intracellular loop as a scaffold for the accumulation of PIP3. In Ifitm3-/- B cells, lipid rafts were depleted of PIP3, which resulted in the defective expression of over 60 lipid-raft-associated surface receptors, and impaired BCR signalling and cellular adhesion. We conclude that the phosphorylation of IFITM3 that occurs after B cells encounter antigen induces a dynamic switch from antiviral effector functions in endosomes to a PI3K amplification loop at the cell surface. IFITM3-dependent amplification of PI3K signalling, which in part acts downstream of the BCR, is critical for the rapid expansion of B cells with high affinity to antigen. In addition, multiple oncogenes depend on IFITM3 to assemble PIP3-dependent signalling complexes and amplify PI3K signalling for malignant transformation.

MAJIQlopedia: an encyclopedia of RNA splicing variations in human tissues and cancer.

Quantification of RNA splicing variations based on RNA-Sequencing can reveal tissue- and disease-specific splicing patterns. To study such splicing variations, we introduce MAJIQlopedia, an encyclopedia of splicing variations that encompasses 86 human tissues and 41 cancer datasets. MAJIQlopedia reports annotated and unannotated splicing events for a total of 486 175 alternative splice junctions in normal tissues and 338 317 alternative splice junctions in cancer. This database, available at https://majiq.biociphers.org/majiqlopedia/, includes a user-friendly interface that provides graphical representations of junction usage quantification for each junction across all tissue or cancer types. To demonstrate case usage of MAJIQlopedia, we review splicing variations in genes WT1, MAPT and BIN1, which all have known tissue or cancer-specific splicing variations. We also use MAJIQlopedia to highlight novel splicing variations in FDX1 and MEGF9 in normal tissues, and we uncover a novel exon inclusion event in RPS6KA6 that only occurs in two cancer types. Users can download the database, request the addition of data to the webtool, or install a MAJIQlopedia server to integrate proprietary data. MAJIQlopedia can serve as a reference database for researchers seeking to understand what splicing variations exist in genes of interest, and those looking to understand tissue- or cancer-specific splice isoform usage.

Alternative splicing of its 5'-UTR limits CD20 mRNA translation and enables resistance to CD20-directed immunotherapies.

Aberrant skipping of coding exons in CD19 and CD22 compromises the response to immunotherapy in B-cell malignancies. Here, we showed that the MS4A1 gene encoding human CD20 also produces several messenger RNA (mRNA) isoforms with distinct 5' untranslated regions. Four variants (V1-4) were detected using RNA sequencing (RNA-seq) at distinct stages of normal B-cell differentiation and B-lymphoid malignancies, with V1 and V3 being the most abundant. During B-cell activation and Epstein-Barr virus infection, redirection of splicing from V1 to V3 coincided with increased CD20 positivity. Similarly, in diffuse large B-cell lymphoma, only V3, but not V1, correlated with CD20 protein levels, suggesting that V1 might be translation-deficient. Indeed, the longer V1 isoform contained upstream open reading frames and a stem-loop structure, which cooperatively inhibited polysome recruitment. By modulating CD20 isoforms with splice-switching morpholino oligomers, we enhanced CD20 expression and anti-CD20 antibody rituximab-mediated cytotoxicity in a panel of B-cell lines. Furthermore, reconstitution of CD20-knockout cells with V3 mRNA led to the recovery of CD20 positivity, whereas V1-reconstituted cells had undetectable levels of CD20 protein. Surprisingly, in vitro CD20-directed chimeric antigen receptor T cells were able to kill both V3- and V1-expressing cells, but the bispecific T-cell engager mosunetuzumab was only effective against V3-expressing cells. To determine whether CD20 splicing is involved in immunotherapy resistance, we performed RNA-seq on 4 postmosunetuzumab follicular lymphoma relapses and discovered that in 2 of them, the downregulation of CD20 was accompanied by a V3-to-V1 shift. Thus, splicing-mediated mechanisms of epitope loss extend to CD20-directed immunotherapies.

Aberrant splicing in B-cell acute lymphoblastic leukemia.

Aberrant splicing is a hallmark of leukemias with mutations in splicing factor (SF)-encoding genes. Here we investigated its prevalence in pediatric B-cell acute lymphoblastic leukemias (B-ALL), where SFs are not mutated. By comparing these samples to normal pro-B cells, we found thousands of aberrant local splice variations (LSVs) per sample, with 279 LSVs in 241 genes present in every comparison. These genes were enriched in RNA processing pathways and encoded ∼100 SFs, e.g. hnRNPA1. HNRNPA1 3'UTR was most pervasively mis-spliced, yielding the transcript subject to nonsense-mediated decay. To mimic this event, we knocked it down in B-lymphoblastoid cells and identified 213 hnRNPA1-regulated exon usage events comprising the hnRNPA1 splicing signature in pediatric leukemia. Some of its elements were LSVs in DICER1 and NT5C2, known cancer drivers. We searched for LSVs in other leukemia and lymphoma drivers and discovered 81 LSVs in 41 additional genes. Seventy-seven LSVs out of 81 were confirmed using two large independent B-ALL RNA-seq datasets, and the twenty most common B-ALL drivers, including NT5C2, showed higher prevalence of aberrant splicing than of somatic mutations. Thus, post-transcriptional deregulation of SF can drive widespread changes in B-ALL splicing and likely contributes to disease pathogenesis.

The miR-17-92 microRNA cluster regulates multiple components of the TGF-ß pathway in neuroblastoma.

The miR-17-92 microRNA cluster is often activated in cancer cells, but the identity of its targets remains elusive. Using SILAC and quantitative mass spectrometry, we examined the effects of activation of the miR-17-92 cluster on global protein expression in neuroblastoma (NB) cells. Our results reveal cooperation between individual miR-17-92 miRNAs and implicate miR-17-92 in multiple hallmarks of cancer, including proliferation and cell adhesion. Most importantly, we show that miR-17-92 is a potent inhibitor of TGF-ß signaling. By functioning both upstream and downstream of pSMAD2, miR-17-92 activation triggers downregulation of multiple key effectors along the TGF-ß signaling cascade as well as direct inhibition of TGF-ß-responsive genes.

miR-17-92 cluster components analysis in Burkitt lymphoma: overexpression of miR-17 is associated with poor prognosis.

Burkitt lymphoma (BL) is an aggressive B cell lymphoma characterized by the reciprocal translocation of the c-Myc gene with immunoglobulin genes. Recently, MYC has been shown to maintain the neoplastic state via the miR-17-92 microRNA cluster that suppresses chromatin regulatory genes and the apoptosis regulator Bim. However, the expression and prognostic impact of miR-17-92 members in pediatric BL (pBL) are unknown. Therefore, we investigated miR-17, miR-19a, miR-19b, miR-20, and miR-92a expression and prognostic impact in a series of 41 pBL samples. In addition, Bim protein expression was evaluated and compared to miR-17, miR-19a, miR-19b, miR-20, and miR-92a levels and patient outcomes. The expression of miR-17-92 members was evaluated by qPCR and Bim protein by immunohistochemistry. Log-rank test was employed to assess prognostic impact. We found that upregulated expression of miR-17 and miR-20a correlates with lack of pro-apoptotic Bim expression. Patients bearing tumors with upregulated miR-17 displayed decreased overall survival (OS), and multivariate analysis revealed that miR-17 was a significant predictor of shortened OS. Using hairpin inhibitors, we showed that inhibition of miR-17 resulted in enhanced Bim expression in a BL cell line overexpressing the miR-17-92 cluster. Our results describe for the first time miR-17, miR-19a, miR-19b, miR-20a, and miR-92a expression profiles in pBL. The prognostic impact of miR-17 should be validated in a larger series, and may provide new therapeutic avenues in the era of anti-miRNA therapy research. Additional functional studies are further required to understand the specific role of miR-17-92 cluster members in BL.

Targeting of TGFß signature and its essential component CTGF by miR-18 correlates with improved survival in glioblastoma.

The miR-17∼92 cluster is thought to be an oncogene, yet its expression is low in glioblastoma multiforme (GBM) cell lines. This could allow unfettered expression of miR-17∼92 target genes such as connective tissue growth factor (CTGF; or CCN2), which is known to contribute to GBM pathogenesis. Indeed, microRNA-18a (but not other miR-17∼92 members) has a functional site in the CTGF 3' UTR, and its forced reexpression sharply reduces CTGF protein and mRNA levels. Interestingly, it also reduces the levels of CTGF primary transcript. The unexpected effects of miR-18a on CTGF transcription are mediated in part by direct targeting of Smad3 and ensuing weakening of TGFß signaling. Having defined the TGFß signature in GBM cells, we demonstrate a significant anti-correlation between miR-18 and TGFß signaling in primary GBM samples from The Cancer Genome Atlas. Most importantly, high levels of miR-18 combined with low levels of the TGFß metagene correlate with prolonged patient survival. Thus, low expression of the miR-17∼92 cluster, and specifically miR-18a, could significantly contribute to GBM pathogenesis.

The Myc-miR-17-92 axis amplifies B-cell receptor signaling via inhibition of ITIM proteins: a novel lymphomagenic feed-forward loop.

The c-Myc oncoprotein regulates >15% of the human transcriptome and a limited number of microRNAs (miRNAs). Here, we establish that in a human B-lymphoid cell line, Myc-repressed, but not Myc-stimulated, genes are significantly enriched for predicted binding sites of Myc-regulated miRNAs, primarily those comprising the Myc-activated miR-17~92 cluster. Notably, gene set enrichment analysis demonstrates that miR-17∼92 is a major regulator of B-cell receptor (BCR) pathway components. Many of them are immunoreceptor tyrosine inhibitory motif (ITIM)-containing proteins, and ITIM proteins CD22 and FCGR2B were found to be direct targets of miR-17∼92. Consistent with the propensity of ITIM proteins to recruit phosphatases, either MYC or miR-17~92 expression was necessary to sustain phosphorylation of spleen tyrosine kinase (SYK) and the B-cell linker protein (BLNK) upon ligation of the BCR. Further downstream, stimulation of the BCR response by miR-17-92 resulted in the enhanced calcium flux and elevated levels of Myc itself. Notably, inhibition of the miR-17~92 cluster in diffuse large B-cell lymphoma (DLBCL) cell lines diminished the BCR response as measured by SYK and BLNK phosphorylation. Conversely, human DLBCLs of the BCR subtype express higher Myc and mir17hg transcript levels than other subtypes. Hence, the Myc-miR-17-92-BCR axis, frequently affected by genomic rearrangements, constitutes a novel lymphomagenic feed-forward loop.

Augmentation of tumor angiogenesis by a Myc-activated microRNA cluster.

Human adenocarcinomas commonly harbor mutations in the KRAS and MYC proto-oncogenes and the TP53 tumor suppressor gene. All three genetic lesions are potentially pro-angiogenic, as they sustain production of vascular endothelial growth factor (VEGF). Yet Kras-transformed mouse colonocytes lacking p53 formed indolent, poorly vascularized tumors, whereas additional transduction with a Myc-encoding retrovirus promoted vigorous vascularization and growth. In addition, VEGF levels were unaffected by Myc, but enhanced neovascularization correlated with downregulation of anti-angiogenic thrombospondin-1 (Tsp1) and related proteins, such as connective tissue growth factor (CTGF). Both Tsp1 and CTGF are predicted targets for repression by the miR-17-92 microRNA cluster, which was upregulated in colonocytes coexpressing K-Ras and c-Myc. Indeed, miR-17-92 knockdown with antisense 2'-O-methyl oligoribonucleotides partly restored Tsp1 and CTGF expression; in addition, transduction of Ras-only cells with a miR-17-92-encoding retrovirus reduced Tsp1 and CTGF levels. Notably, miR-17-92-transduced cells formed larger, better-perfused tumors. These findings establish a role for microRNAs in non-cell-autonomous Myc-induced tumor phenotypes.

Contrasting and Combining Transcriptome Complexity Captured by Short and Long RNA Sequencing Reads.

Mapping transcriptomic variations using either short or long reads RNA sequencing is a staple of genomic research. Long reads are able to capture entire isoforms and overcome repetitive regions, while short reads still provides improved coverage and error rates. Yet how to quantitatively compare the technologies, can we combine those, and what may be the benefit of such a combined view remain open questions. We tackle these questions by first creating a pipeline to assess matched long and short reads data using a variety of transcriptome statistics. We find that across datasets, algorithms and technologies, matched short reads data detects roughly 50% more splice junctions, with 10-30% of the splice junctions included at 20% or more are missed by long reads. In contrast, long reads detect many more intron retention events, pointing to the benefit of combining the technologies. We introduce MAJIQ-L, an extension of the MAJIQ software to enable a unified view of transcriptome variations from both technologies and demonstrate its benefits. Our software can be used to assess any future long reads technology or algorithm, and combine it with short reads data for improved transcriptome analysis.

A Bayesian model for unsupervised detection of RNA splicing based subtypes in cancers.

Identification of cancer sub-types is a pivotal step for developing personalized treatment. Specifically, sub-typing based on changes in RNA splicing has been motivated by several recent studies. We thus develop CHESSBOARD, an unsupervised algorithm tailored for RNA splicing data that captures "tiles" in the data, defined by a subset of unique splicing changes in a subset of patients. CHESSBOARD allows for a flexible number of tiles, accounts for uncertainty of splicing quantification, and is able to model missing values as additional signals. We first apply CHESSBOARD to synthetic data to assess its domain specific modeling advantages, followed by analysis of several leukemia datasets. We show detected tiles are reproducible in independent studies, investigate their possible regulatory drivers and probe their relation to known AML mutations. Finally, we demonstrate the potential clinical utility of CHESSBOARD by supplementing mutation based diagnostic assays with discovered splicing profiles to improve drug response correlation.

Discovery of antibodies targeting multipass transmembrane proteins using a suspension cell-based evolutionary approach.

Due to their critical functions in cell sensing and signal processing, membrane proteins are highly preferred as pharmacological targets, and antibody drugs constitute the fastest growing category of therapeutic agents on the pharmaceutical market. However, major limitations exist in developing antibodies that recognize complex, multipass transmembrane proteins, such as G-protein-coupled receptors (GPCRs). These challenges, largely due to difficulties with recombinant expression of multipass transmembrane proteins, can be overcome using whole-cell screening techniques, which enable presentation of the functional antigen in its native conformation. Here, we developed suspension cell-based whole-cell panning methodologies to screen for specific binders against GPCRs within a naive yeast-displayed antibody library. We implemented our strategy to discover high-affinity antibodies against four distinct GPCR target proteins, demonstrating the potential for our cell-based screening workflow to advance the discovery of antibody therapeutics targeting membrane proteins.

FBXW7ß isoform drives transcriptional activation of the proinflammatory TNF cluster in human pro-B cells.

Noncanonical exon usage plays many important roles in cellular phenotypes, but its contribution to human B-cell development remains sketchily understood. To fill this gap, we collected various B-cell fractions from bone marrow (BM) and tonsil donors, performed RNA sequencing, and examined transcript variants. We identified 150 genes that harbor local splicing variations in all pairwise comparisons. One of them encodes FBXW7, an E3 ubiquitin ligase implicated as a driver in several blood cancers. Surprisingly, we discovered that in normal human pro-B cells, the predominant transcript used an alternative first exon to produce the poorly characterized FBXW7ß isoform, previously thought to be restricted to neural tissues. The FBXW7ß transcript was also abundant in cell lines and primary samples of pediatric B-cell acute lymphoblastic leukemia (B-ALL), which originates in the BM. When overexpressed in a heterologous cell system, this transcript yielded the expected protein product, as judged by anti-FLAG immunoblotting and mass spectrometry. Furthermore, in REH B-ALL cells, FBXW7ß mRNA was the only FBXW7 isoform enriched in the polyribosome fraction. To shed light on possible functions of FBXW7ß, we used gain- and loss-of-function approaches and identified an FBXW7-dependent inflammatory gene signature, apparent in a subset of B-ALL with high FBXW7ß expression. This signature contained several members of the tumor necrosis factor superfamily, including those comprising the HLA Class III cluster (LTB, LST1, NCR3, LTA, and NFKBIL1). Our findings suggest that FBXW7ß expression drives proinflammatory responses, which could contribute to normal B-cell development, leukemogenesis, and responses to anticancer therapies.

Alternative splicing of its 5'-UTR limits CD20 mRNA translation and enables resistance to CD20-directed immunotherapies.

Aberrant skipping of coding exons in CD19 and CD22 compromises responses to immunotherapy for B-cell malignancies. Here, we show that the MS4A1 gene encoding human CD20 also produces several mRNA isoforms with distinct 5' untranslated regions (5'-UTR). Four variants (V1-4) were detectable by RNA-seq in distinct stages of normal B-cell differentiation and B-lymphoid malignancies, with V1 and V3 being the most abundant by far. During B-cell activation and Epstein-Barr virus infection, redirection of splicing from V1 to V3 coincided with increased CD20 positivity. Similarly, in diffuse large B-cell lymphoma only V3, but not V1, correlated with CD20 protein levels, suggesting that V1 might be translation-deficient. Indeed, the longer V1 isoform was found to contain upstream open reading frames (uORFs) and a stem-loop structure, which cooperatively inhibited polysome recruitment. By modulating CD20 isoforms with splice-switching Morpholino oligomers, we enhanced CD20 expression and anti-CD20 antibody rituximab-mediated cytotoxicity in a panel of B-cell lines. Furthermore, reconstitution of CD20-knockout cells with V3 mRNA led to the recovery of CD20 positivity, while V1-reconstituted cells had undetectable levels of CD20 protein. Surprisingly, in vitro CD20-directed CAR T cells were able to kill both V3- and V1-expressing cells, but the bispecific T cell engager mosunetuzumab was only effective against V3-expressing cells. To determine whether CD20 splicing is involved in immunotherapy resistance, we performed RNA-seq on four post-mosunetuzumab follicular lymphoma relapses and discovered that in two of them downregulation of CD20 was accompanied by the V3-to-V1 shift. Thus, splicing-mediated mechanisms of epitope loss extend to CD20-directed immunotherapies. Key Points: In normal & malignant human B cells, CD20 mRNA is alternatively spliced into four 5'-UTR isoforms, some of which are translation-deficient.The balance between translation-deficient and -competent isoforms modulates CD20 protein levels & responses to CD20-directed immunotherapies. Explanation of Novelty: We discovered that in normal and malignant B-cells, CD20 mRNA is alternatively spliced to generate four distinct 5'-UTRs, including the longer translation-deficient V1 variant. Cells predominantly expressing V1 were still sensitive to CD20-targeting chimeric antigen receptor T-cells. However, they were resistant to the bispecific anti-CD3/CD20 antibody mosunetuzumab, and the shift to V1 were observed in CD20-negative post-mosunetuzumab relapses of follicular lymphoma.

Mapping and modeling human colorectal carcinoma interactions with the tumor microenvironment.

The initiation and progression of cancer are intricately linked to the tumor microenvironment (TME). Understanding the function of specific cancer-TME interactions poses a major challenge due in part to the complexity of the in vivo microenvironment. Here we predict cancer-TME interactions from single cell transcriptomic maps of both human colorectal cancers (CRCs) and mouse CRC models, ask how these interactions are altered in human tumor organoid (tumoroid) cultures, and functionally recapitulate human myeloid-carcinoma interactions in vitro. Tumoroid cultures suppress gene expression programs involved in inflammation and immune cell migration, providing a reductive platform for re-establishing carcinoma-immune cell interactions in vitro. Introduction of human monocyte-derived macrophages into tumoroid cultures instructs macrophages to acquire immunosuppressive and pro-tumorigenic gene expression programs similar to those observed in vivo. This includes hallmark induction of SPP1, encoding Osteopontin, an extracellular CD44 ligand with established oncogenic effects. Taken together, these findings offer a framework for understanding CRC-TME interactions and provide a reductionist tool for modeling specific aspects of these interactions.