TransCRISPR-sgRNA design tool for CRISPR/Cas9 experiments targeting specific sequence motifs.

Eukaryotic genomes contain several types of recurrent sequence motifs, e.g. transcription factor motifs, miRNA binding sites, repetitive elements. CRISPR/Cas9 can facilitate identification and study of crucial motifs. We present transCRISPR, the first online tool dedicated to search for sequence motifs in the user-provided genomic regions and design optimal sgRNAs targeting them. Users can obtain sgRNAs for chosen motifs, for up to tens of thousands of target regions in 30 genomes, either for the Cas9 or dCas9 system. TransCRISPR provides user-friendly tables and visualizations, summarizing features of identified motifs and designed sgRNAs such as genomic localization, quality scores, closest transcription start sites and others. Experimental validation of sgRNAs for MYC binding sites designed with transCRISPR confirmed efficient disruption of the targeted motifs and effect on expression of MYC-regulated genes. TransCRISPR is available from https://transcrispr.igcz.poznan.pl/transcrispr/.

The lncRNA KTN1-AS1 co-regulates a variety of Myc-target genes and enhances proliferation of Burkitt lymphoma cells.

Long non-coding RNAs (lncRNAs) are involved in many normal and oncogenic pathways through a diverse repertoire of transcriptional and posttranscriptional regulatory mechanisms. LncRNAs that are under tight regulation of well-known oncogenic transcription factors such as c-Myc (Myc) are likely to be functionally involved in their disease-promoting mechanisms. Myc is a major driver of many subsets of B cell lymphoma and to date remains an undruggable target. We identified three Myc-induced and four Myc-repressed lncRNAs by use of multiple in vitro models of Myc-driven Burkitt lymphoma and detailed analysis of Myc binding profiles. We show that the top Myc-induced lncRNA KTN1-AS1 is strongly upregulated in different types of B cell lymphoma compared with their normal counterparts. We used CRISPR-mediated genome editing to confirm that the direct induction of KTN1-AS1 by Myc is dependent on the presence of a Myc E-box-binding motif. Knockdown of KTN1-AS1 revealed a strong negative effect on the growth of three BL cell lines. Global gene expression analysis upon KTN1-AS1 depletion shows a strong enrichment of key genes in the cholesterol biosynthesis pathway as well as co-regulation of many Myc-target genes, including a moderate negative effect on the levels of Myc itself. Our study suggests a critical role for KTN1-AS1 in supporting BL cell growth by mediating co-regulation of a variety of Myc-target genes and co-activating key genes involved in cholesterol biosynthesis. Therefore, KTN1-AS1 may represent a putative novel therapeutic target in lymphoma.

Core regions in immunoglobulin heavy chain enhancers essential for survival of non-Hodgkin lymphoma cells are identified by a CRISPR interference screen.

Chromosomal translocations in non-Hodgkin lymphoma (NHL) result in activation of oncogenes by placing them under the regulation of immunoglobulin heavy chain (IGH) super-enhancers. Aberrant expression of translocated oncogenes induced by enhancer activity can contribute to lymphomagenesis. The role of the IGH enhancers in normal B-cell development is well established, but knowledge regarding the precise mechanisms of their involvement in control of the translocated oncogenes is limited. The goal of this project was to define the critical regions in the IGH regulatory elements and identify enhancer RNAs (eRNA). We designed a sgRNA library densely covering the IGH enhancers and performed tiling CRISPR interference screens in three NHL cell lines. This revealed three regions crucial for NHL cell growth. With chromatin-enriched RNA-Seq we showed transcription from the core enhancer regions and subsequently validated expression of the eRNAs in a panel of NHL cell lines and tissue samples. Inhibition of the essential IGH enhancer regions decreased expression of eRNAs and translocated oncogenes in several NHL cell lines. The observed expression and growth patterns were consistent with the breakpoints in the IGH locus. Moreover, targeting the Eµ enhancer resulted in loss of B-cell receptor expression. In a Burkitt lymphoma cell line, MYC overexpression partially rescued the phenotype induced by IGH enhancer inhibition. Our results indicated the most critical regions in the IGH enhancers and provided new insights into the current understanding of the role of IGH enhancers in B-cell NHL. As such, this study forms a basis for development of potential therapeutic approaches.

Inhibition of the glutamate-cysteine ligase catalytic subunit with buthionine sulfoximine enhances the cytotoxic effect of doxorubicin and cyclophosphamide in Burkitt lymphoma cells.

Burkitt lymphoma (BL) is a highly aggressive lymphoma that mainly affects children and young adults. Chemotherapy is effective in young BL patients but the outcome in adults is less satisfactory. Therefore, there is a need to enhance the cytotoxic effect of drugs used in BL treatment. Glutathione (GSH) is an important antioxidant involved in processes such as regulation of oxidative stress and drug detoxification. Elevated GSH levels have been observed in many cancers and were associated with chemoresistance. We previously identified GCLC, encoding an enzyme involved in GSH biosynthesis, as an essential gene in BL. We now confirm that knockout of GCLC decreases viability of BL cells and that the GCLC protein is overexpressed in BL tissues. Moreover, we demonstrate that buthionine sulfoximine (BSO), a known inhibitor of GCLC, decreases growth of BL cells but does not affect control B cells. Furthermore, we show for the first time that BSO enhances the cytotoxicity of compounds commonly used in BL treatment, doxorubicin, and cyclophosphamide. Given the fact that BSO itself was not toxic to control cells and well-tolerated in clinical trials, combination of chemotherapy with BSO may allow reduction of the doses of cytotoxic drugs required to obtain effective responses in BL patients.

Identification of ATM-dependent long non-coding RNAs induced in response to DNA damage.

DNA damage response (DDR) is a complex process, essential for cell survival. Especially deleterious type of DNA damage are DNA double-strand breaks (DSB), which can lead to genomic instability and malignant transformation if not repaired correctly. The central player in DSB detection and repair is the ATM kinase which orchestrates the action of several downstream factors. Recent studies have suggested that long non-coding RNAs (lncRNAs) are involved in DDR. Here, we aimed to identify lncRNAs induced upon DNA damage in an ATM-dependent manner. DNA damage was induced by ionizing radiation (IR) in immortalized lymphoblastoid cell lines derived from 4 patients with ataxia-telangiectasia (AT) and 4 healthy donors. RNA-seq revealed 10 lncRNAs significantly induced 1 h after IR in healthy donors, whereas none in AT patients. 149 lncRNAs were induced 8 h after IR in the control group, while only three in AT patients. Among IR-induced mRNAs, we found several genes with well-known functions in DDR. Gene Set Enrichment Analysis and Gene Ontology revealed delayed induction of key DDR pathways in AT patients compared to controls. The induction and dynamics of selected 9 lncRNAs were confirmed by RT-qPCR. Moreover, using a specific ATM inhibitor we proved that the induction of those lncRNAs is dependent on ATM. Some of the detected lncRNA genes are localized next to protein-coding genes involved in DDR. We observed that induction of lncRNAs after IR preceded changes in expression of adjacent genes. This indicates that IR-induced lncRNAs may regulate the transcription of nearby genes. Subcellular fractionation into chromatin, nuclear, and cytoplasmic fractions revealed that the majority of studied lncRNAs are localized in chromatin. In summary, our study revealed several lncRNAs induced by IR in an ATM-dependent manner. Their genomic co-localization and co-expression with genes involved in DDR suggest that those lncRNAs may be important players in cellular response to DNA damage.

CRISPR/Cas9 screen for genome-wide interrogation of essential MYC-bound E-boxes in cancer cells.

The transcription factor MYC is a proto-oncogene with a well-documented essential role in the pathogenesis and maintenance of several types of cancer. MYC binds to specific E-box sequences in the genome to regulate gene expression in a cell-type- and developmental-stage-specific manner. To date, a combined analysis of essential MYC-bound E-boxes and their downstream target genes important for growth of different types of cancer is missing. In this study, we designed a CRISPR/Cas9 library to destroy E-box sequences in a genome-wide fashion. In parallel, we used the Brunello library to knock out protein-coding genes. We performed high-throughput screens with these libraries in four MYC-dependent cancer cell lines-K562, ST486, HepG2, and MCF7-which revealed several essential E-boxes and genes. Among them, we pinpointed crucial common and cell-type-specific MYC-regulated genes involved in pathways associated with cancer development. Extensive validation of our approach confirmed that E-box disruption affects MYC binding, target-gene expression, and cell proliferation in vitro as well as tumor growth in vivo. Our unique, well-validated tool opens new possibilities to gain novel insights into MYC-dependent vulnerabilities in cancer cells.

7-[[(4-methyl-2-pyridinyl)amino](2-pyridinyl)methyl]-8-quinolinol (compound 30666) inhibits enhancer activity and reduces B-cell lymphoma growth - A question of specificity.

B-cell non-Hodgkin lymphoma (NHL) is among the ten most common malignancies. Survival rates range from very poor to over 90% and highly depend on the stage and subtype. Characteristic features of NHL are recurrent translocations juxtaposing an oncogene (e.g. MYC, BCL2) to the enhancers in the immunoglobulin heavy chain (IGH) locus. Survival and proliferation of many B-cell lymphomas depend on the expression of the translocated oncogene. Thus, targeting IGH enhancers as an anti-lymphoma treatment seems a promising strategy. Recently, a small molecule - 7-[[(4-methyl-2-pyridinyl)amino](2-pyridinyl)methyl]-8-quinolinol (compound 30666) was identified to decrease activity of the Eµ enhancer and reduce the expression of translocated oncogenes in multiple myeloma and some NHL cell lines (Dolloff, 2019). Here, we aimed to test the effect of compound 30666 in Burkitt lymphoma (BL) and diffuse large B-cell lymphoma (DLBCL) and shed light on its mechanism of action. We report that both IGH-translocation positive NHL cells as well as IGH-translocation negative B cells and non-B cell controls treated with compound 30666 exhibited consistent growth inhibition. A statistically significant increase in cell percentage in sub-G1 phase of cell cycle was observed, suggesting induction of apoptosis. Compound 30666 downregulated MYC levels in BL cell lines and altered IGH enhancer RNA expression. Moreover, a global decrease of H3K27ac and an increase of H3K4me1 was observed upon 30666 treatment, which suggests switching enhancers to a poised or primed state. Altogether, our findings indicate that 30666 inhibitor affects enhancer activity but might not be as specific for IGH enhancers as previously reported.

The miR-26b-5p/KPNA2 Axis Is an Important Regulator of Burkitt Lymphoma Cell Growth.

The expression of several microRNAs (miRNAs) is known to be changed in Burkitt lymphoma (BL), compared to its normal counterparts. Although for some miRNAs, a role in BL was demonstrated, for most of them, their function is unclear. In this study, we aimed to identify miRNAs that control BL cell growth. Two BL cell lines were infected with lentiviral pools containing either 58 miRNA inhibitors or 44 miRNA overexpression constructs. Eighteen constructs showed significant changes in abundance over time, indicating that they affected BL growth. The screening results were validated by individual green fluorescent protein (GFP) growth competition assays for fifteen of the eighteen constructs. For functional follow-up studies, we focused on miR-26b-5p, whose overexpression inhibited BL cell growth. Argonaute 2 RNA immunoprecipitation (Ago2-IP) in two BL cell lines revealed 47 potential target genes of miR-26b-5p. Overlapping the list of putative targets with genes showing a growth repression phenotype in a genome-wide CRISPR/Cas9 knockout screen, revealed eight genes. The top-5 candidates included EZH2, COPS2, KPNA2, MRPL15, and NOL12. EZH2 is a known target of miR-26b-5p, with oncogenic properties in BL. The relevance of the latter four targets was confirmed using sgRNAs targeting these genes in individual GFP growth competition assays. Luciferase reporter assay confirmed binding of miR-26b-5p to the predicted target site for KPNA2, but not to the other genes. In summary, we identified 18 miRNAs that affected BL cell growth in a loss- or gain-of-function screening. A tumor suppressor role was confirmed for miR-26b-5p, and this effect could at least in part be attributed to KPNA2, a known regulator of OCT4, c-jun, and MYC.

MAVS is not a Likely Susceptibility Locus for Addison's Disease and Type 1 Diabetes.

Mitochondrial antiviral signaling (MAVS) protein is an intracellular adaptor molecule, downstream of viral sensors, retinoid acid-inducible gene I (RIG-I)-like receptors (RLRs). Impaired antiviral cell signaling might contribute to autoimmunity. Studies have recently shown variations in genes encoding RLRs as risk factors for autoimmune diseases. We investigated whether MAVS coding polymorphisms are associated with Addison's disease (AD) and type 1 diabetes (T1D) in Polish population. We genotyped 140 AD, 532 T1D patients and 600 healthy controls for MAVS rs17857295, rs7262903, rs45437096 and rs7269320. Genotyping was performed by TaqMan assays. Distribution of the MAVS genotypes and alleles did not reveal significant differences between patients and controls (p > 0.05). This analysis did not indicate the association of the MAVS locus with susceptibility to AD and T1D.