A multiomic characterization of the leukemia cell line REH using short- and long-read sequencing.

The B-cell acute lymphoblastic leukemia (ALL) cell line REH, with the t(12;21) ETV6::RUNX1 translocation, is known to have a complex karyotype defined by a series of large-scale chromosomal rearrangements. Taken from a 15-yr-old at relapse, the cell line offers a practical model for the study of pediatric B-ALL. In recent years, short- and long-read DNA and RNA sequencing have emerged as a complement to karyotyping techniques in the resolution of structural variants in an oncological context. Here, we explore the integration of long-read PacBio and Oxford Nanopore whole-genome sequencing, IsoSeq RNA sequencing, and short-read Illumina sequencing to create a detailed genomic and transcriptomic characterization of the REH cell line. Whole-genome sequencing clarified the molecular traits of disrupted ALL-associated genes including CDKN2A, PAX5, BTG1, VPREB1, and TBL1XR1, as well as the glucocorticoid receptor NR3C1 Meanwhile, transcriptome sequencing identified seven fusion genes within the genomic breakpoints. Together, our extensive whole-genome investigation makes high-quality open-source data available to the leukemia genomics community.

Multimodal Single-Cell Sequencing of B Cells in Primary Sjögren's Syndrome.

OBJECTIVE: B cells are important in the pathogenesis of primary Sjögren's syndrome (pSS). Patients positive for Sjögren's syndrome antigen A/Sjögren syndrome antigen B (SSA/SSB) autoantibodies are more prone to systemic disease manifestations and adverse outcomes. We aimed to determine the role of B cell composition, gene expression, and B cell receptor usage in pSS subgroups stratified for SSA/SSB antibodies. METHODS: Over 230,000 B cells were isolated from peripheral blood of patients with pSS (n = 6 SSA-, n = 8 SSA+ single positive and n = 10 SSA/SSB+ double positive) and four healthy controls and processed for single-cell RNA sequencing (scRNA-seq) and single-cell variable, diversity, and joining (VDJ) gene sequencing (scVDJ-seq). RESULTS: We show that SSA/SSB+ patients present the highest and lowest proportion of naïve and memory B cells, respectively, and the highest up-regulation of interferon-induced genes across all B cell subtypes. Differential usage of IGHV showed that IGHV1-69 and IGHV4-30-4 were more often used in all pSS subgroups compared with controls. Memory B cells from SSA/SSB+ patients displayed a higher proportion of cells with unmutated VDJ transcripts compared with other pSS patient groups and controls, indicating altered somatic hypermutation processes. Comparison with previous studies revealed heterogeneous clonotype pools, with little overlap in CDR3 sequences. Joint analysis using scRNA-seq and scVDJ-seq data allowed unsupervised stratification of patients with pSS and identified novel parameters that correlated to disease manifestations and antibody status. CONCLUSION: We describe heterogeneity and molecular characteristics in B cells from patients with pSS, providing clues to intrinsic differences in B cells that affect the phenotype and outcome and allowing stratification of patients with pSS at improved resolution.

Multimodal classification of molecular subtypes in pediatric acute lymphoblastic leukemia.

Genomic analyses have redefined the molecular subgrouping of pediatric acute lymphoblastic leukemia (ALL). Molecular subgroups guide risk-stratification and targeted therapies, but outcomes of recently identified subtypes are often unclear, owing to limited cases with comprehensive profiling and cross-protocol studies. We developed a machine learning tool (ALLIUM) for the molecular subclassification of ALL in retrospective cohorts as well as for up-front diagnostics. ALLIUM uses DNA methylation and gene expression data from 1131 Nordic ALL patients to predict 17 ALL subtypes with high accuracy. ALLIUM was used to revise and verify the molecular subtype of 281 B-cell precursor ALL (BCP-ALL) cases with previously undefined molecular phenotype, resulting in a single revised subtype for 81.5% of these cases. Our study shows the power of combining DNA methylation and gene expression data for resolving ALL subtypes and provides a comprehensive population-based retrospective cohort study of molecular subtype frequencies in the Nordic countries.

Increased MYB alternative promoter usage is associated with relapse in acute lymphoblastic leukemia.

Therapy-resistant disease is a major cause of death in patients with acute lymphoblastic leukemia (ALL). Activation of the MYB oncogene is associated with ALL and leads to uncontrolled neoplastic cell proliferation and blocked differentiation. Here, we used RNA-seq to study the clinical significance of MYB expression and MYB alternative promoter (TSS2) usage in 133 pediatric ALLs. RNA-seq revealed that all cases analyzed overexpressed MYB and demonstrated MYB TSS2 activity. qPCR analyses confirmed the expression of the alternative MYB promoter also in seven ALL cell lines. Notably, high MYB TSS2 activity was significantly associated with relapse (p = 0.007). Moreover, cases with high MYB TSS2 usage showed evidence of therapy-resistant disease with increased expression of ABC multidrug resistance transporter genes (e.g., ABCA2, ABCB5, and ABCC10) and enzymes catalyzing drug degradation (e.g., CYP1A2, CYP2C9, and CYP3A5). Elevated MYB TSS2 activity was further associated with augmented KRAS signaling (p < 0.05) and decreased methylation of the conventional MYB promoter (p < 0.01). Taken together, our results suggest that MYB alternative promoter usage is a novel potential prognostic biomarker for relapse and therapy resistance in pediatric ALL.

Differential Transcriptional Reprogramming by Wild Type and Lymphoma-Associated Mutant MYC Proteins as B-Cells Convert to a Lymphoma Phenotype.

The MYC transcription factor regulates a vast number of genes and is implicated in many human malignancies. In some hematological malignancies, MYC is frequently subject to missense mutations that enhance its transformation activity. Here, we use a novel murine cell system to (i) characterize the transcriptional effects of progressively increasing MYC levels as normal primary B-cells transform to lymphoma cells and (ii) determine how this gene regulation program is modified by lymphoma-associated MYC mutations (T58A and T58I) that enhance its transformation activity. Unlike many previous studies, the cell system exploits primary B-cells that are transduced to allow regulated MYC expression under circumstances where apoptosis and senescence pathways are abrogated by the over-expression of the Bcl-xL and BMI1 proteins. In such cells, transition from a normal to a lymphoma phenotype is directly dependent on the MYC expression level, without a requirement for secondary events that are normally required during MYC-driven oncogenic transformation. A generalized linear model approach allowed an integrated analysis of RNA sequencing data to identify regulated genes in relation to both progressively increasing MYC level and wild type or mutant status. Using this design, a total of 7569 regulated genes were identified, of which the majority (n = 7263) were regulated in response to progressively increased levels of wild type MYC, while a smaller number of genes (n = 917) were differentially regulated, compared to wild type MYC, in T58A MYC- and/or T58I MYC-expressing cells. Unlike most genes that are similarly regulated by both wild type and mutant MYC genes, the set of 917 genes did not significantly overlap with known lipopolysaccharide regulated genes, which represent genes regulated by MYC in normal B cells. The genes that were differently regulated in cells expressing mutant MYC proteins were significantly enriched in DNA replication and G2 phase to mitosis transition genes. Thus, mutants affecting MYC proteins may augment quantitative oncogenic effects on the expression of normal MYC-target genes with qualitative oncogenic effects, by which sets of cell cycle genes are abnormally targeted by MYC as B cells transition into lymphoma cells. The T58A and T58I mutations augment MYC-driven transformation by distinct mechanisms.

3D heterospecies spheroids of pancreatic stroma and cancer cells demonstrate key phenotypes of pancreatic ductal adenocarcinoma.

Pancreatic ductal adenocarcinoma (PDAC) is one of the most lethal malignancies, partly due to the dense desmoplasia and a lack of suitable model systems to study. In the present work, we developed a 3D heterospecies spheroid model to study the microenvironmental interactions between tumor cells and stellate cells which can also be employed to test therapeutic regimens. We set up monospheroids and heterospheroids made up from murine pancreatic stellate cells (mPSCs) and human PDAC cells (Panc1), which allowed for direct isolation of mRNA from a mixed cell population followed by an in silico separation of the RNA-seq reads. Global transcript level changes for cells in heterospheroids versus monospheroids were calculated, followed by gene set enrichment analysis and molecular subtype analysis. We observed an apparent shift of Panc1 from the classical to the squamous/basal-like phenotype upon co-culture with mPSCs. Moreover, mPSCs acquired a different cancer-associated fibroblast-related phenotype upon co-culture with Panc1. We analyzed the tumor cell-specific chemosensitivities towards gemcitabine, paclitaxel and SN38 and compared these to published pharmacotranscriptomic signatures. In conclusion, our heterospecies spheroid model reflected key aspects of PDAC and facilitated the study of intercellular interactions between tumor and stroma while additionally proving to be a good model for studying therapeutic responses.

Differential B-Cell Receptor Signaling Requirement for Adhesion of Mantle Cell Lymphoma Cells to Stromal Cells.

Interactions between lymphoma cells and stromal cells play a key role in promoting tumor survival and development of drug resistance. We identified differences in key signaling pathways between the JeKo-1 and REC-1 mantle cell lymphoma (MCL) cell lines, displaying different patterns of stromal cell adhesion and chemotaxis towards stroma-conditioned medium. The identified adhesion-regulated genes reciprocated important aspects of microenvironment-mediated gene modulation in MCL patients. Five-hundred and ninety genes were differently regulated between the cell lines upon adhesion to stromal cells, while 32 genes were similarly regulated in both cell lines. Regulation of B-cell Receptor (BCR) signature genes in adherent cells was specific for JeKo-1. Inhibition of BCR using siRNA or clinically approved inhibitors, Ibrutinib and Acalabrutinib, decreased adhesion of JeKo-1, but not REC-1 cells. Cell surface levels of chemokine receptor CXCR4 were higher in JeKo-1, facilitating migration and adhesion of JeKo-1 but not REC-1 cells. Surface levels of ICAM1 adhesion protein differ for REC-1 and JeKo-1. While ICAM1 played a positive role in adherence of both cell lines to stromal cells, S1PR1 had an inhibitory effect. Our results provide a model framework for further investigation of mechanistic differences in patient-response to new pathway-specific drugs.

A Protein Intrinsic Disorder Approach for Characterising Differentially Expressed Genes in Transcriptome Data: Analysis of Cell-Adhesion Regulated Gene Expression in Lymphoma Cells.

Conformational protein properties are coupled to protein functionality and could provide a useful parameter for functional annotation of differentially expressed genes in transcriptome studies. The aim was to determine whether predicted intrinsic protein disorder was differentially associated with proteins encoded by genes that are differentially regulated in lymphoma cells upon interaction with stromal cells, an interaction that occurs in microenvironments, such as lymph nodes that are protective for lymphoma cells during chemotherapy. Intrinsic disorder protein properties were extracted from the Database of Disordered Protein Prediction (D²P²), which contains data from nine intrinsic disorder predictors. Proteins encoded by differentially regulated cell-adhesion regulated genes were enriched in intrinsically disordered regions (IDRs) compared to other genes both with regard to IDR number and length. The enrichment was further ascribed to down-regulated genes. Consistently, a higher proportion of proteins encoded by down-regulated genes contained at least one IDR or were completely disordered. We conclude that down-regulated genes in stromal cell-adherent lymphoma cells encode proteins that are characterized by elevated levels of intrinsically disordered conformation, indicating the importance of down-regulating functional mechanisms associated with intrinsically disordered proteins in these cells. Further, the approach provides a generally applicable and complementary alternative to classification of differentially regulated genes using gene ontology or pathway enrichment analysis.

Mixed-species RNAseq analysis of human lymphoma cells adhering to mouse stromal cells identifies a core gene set that is also differentially expressed in the lymph node microenvironment of mantle cell lymphoma and chronic lymphocytic leukemia patients.

A subset of hematologic cancer patients is refractory to treatment or suffers relapse, due in part to minimal residual disease, whereby some cancer cells survive treatment. Cell-adhesion-mediated drug resistance is an important mechanism, whereby cancer cells receive survival signals via interaction with e.g. stromal cells. No genome-wide studies of in vitro systems have yet been performed to compare gene expression in different cell subsets within a co-culture and cells grown separately. Using RNA sequencing and species-specific read mapping, we compared transcript levels in human Jeko-1 mantle cell lymphoma cells stably adhered to mouse MS-5 stromal cells or in suspension within a co-culture or cultured separately as well as in stromal cells in co-culture or in separate culture. From 1050 differentially expressed transcripts in adherent mantle cell lymphoma cells, we identified 24 functional categories that together represent four main functional themes, anti-apoptosis, B-cell signaling, cell adhesion/migration and early mitosis. A comparison with previous mantle cell lymphoma and chronic lymphocytic leukemia studies, of gene expression differences between lymph node and blood, identified 116 genes that are differentially expressed in all three studies. From these genes, we suggest a core set of genes (CCL3, CCL4, DUSP4, ETV5, ICAM1, IL15RA, IL21R, IL4I1, MFSD2A, NFKB1, NFKBIE, SEMA7A, TMEM2) characteristic of cells undergoing cell-adhesion-mediated microenvironment signaling in mantle cell lymphoma/chronic lymphocytic leukemia. The model system developed and characterized here together with the core gene set will be useful for future studies of pathways that mediate increased cancer cell survival and drug resistance mechanisms.