Transcriptomic signatures reveal a shift towards an anti-inflammatory gene expression profile but also the induction of type I and type II interferon signaling networks through aryl hydrocarbon receptor activation in murine macrophages.

Introduction: The aryl hydrocarbon receptor (AhR) is a ligand-activated transcription factor that regulates a broad range of target genes involved in the xenobiotic response, cell cycle control and circadian rhythm. AhR is constitutively expressed in macrophages (Mϕ), acting as key regulator of cytokine production. While proinflammatory cytokines, i.e., IL-1ß, IL-6, IL-12, are suppressed through AhR activation, anti-inflammatory IL-10 is induced. However, the underlying mechanisms of those effects and the importance of the specific ligand structure are not yet completely understood. Methods: Therefore, we have compared the global gene expression pattern in activated murine bone marrow-derived macrophages (BMMs) subsequently to exposure with either benzo[a]pyrene (BaP) or indole-3-carbinol (I3C), representing high-affinity vs. low-affinity AhR ligands, respectively, by means of mRNA sequencing. AhR dependency of observed effects was proved using BMMs from AhR-knockout (Ahr-/-) mice. Results and discussion: In total, more than 1,000 differentially expressed genes (DEGs) could be mapped, covering a plethora of AhR-modulated effects on basal cellular processes, i.e., transcription and translation, but also immune functions, i.e., antigen presentation, cytokine production, and phagocytosis. Among DEGs were genes that are already known to be regulated by AhR, i.e., Irf1, Ido2, and Cd84. However, we identified DEGs not yet described to be AhR-regulated in Mϕ so far, i.e., Slpi, Il12rb1, and Il21r. All six genes likely contribute to shifting the Mϕ phenotype from proinflammatory to anti-inflammatory. The majority of DEGs induced through BaP were not affected through I3C exposure, probably due to higher AhR affinity of BaP in comparison to I3C. Mapping of known aryl hydrocarbon response element (AHRE) sequence motifs in identified DEGs revealed more than 200 genes not possessing any AHRE, and therefore being not eligible for canonical regulation. Bioinformatic approaches modeled a central role of type I and type II interferons in the regulation of those genes. Additionally, RT-qPCR and ELISA confirmed a AhR-dependent expressional induction and AhR-dependent secretion of IFN-γ in response to BaP exposure, suggesting an auto- or paracrine activation pathway of Mϕ.

Atypical fibroxanthoma and pleomorphic dermal sarcoma - gene expression analysis compared with undifferentiated cutaneous squamous cell carcinoma.

BACKGROUND: The histogenetic origin of atypical fibroxanthoma (AFX) and pleomorphic dermal sarcoma (PDS) has not been definitively elucidated. In addition to a fibroblastic origin, a keratinocytic differentiation is discussed due to strong clinical, histomorphological and molecular genetic similarities with undifferentiated cutaneous squamous cell carcinoma (cSCC). PATIENTS AND METHODS: 56 cases (36 AFXs, 8 PDSs, 12 undifferentiated cSCCs) were evaluated for their clinical, histomorphological, and immunohistochemical characteristics. RNA transcriptome analysis was performed on 18 cases (6 AFXs/PDSs, 6 undifferentiated cSCCs, 6 differentiated cSCCs). RESULTS: Clinically, the strong similarities in age, gender and tumor location were confirmed. Without further immunohistochemical staining, histomorphological differentiation between AFX/PDS and undifferentiated cSCC is often impossible. Principal component analysis of the RNA transcriptome analysis showed that AFX/PDS and differentiated cSCC each formed their own cluster, while the undifferentiated cSCCs fall in between these two groups, but without forming a cluster of their own. When examining differentially expressed genes (DEGs), the heat maps showed that there were cases within the undifferentiated cSCC that were more likely to be AFX/PDS than differentiated cSCC based on their expression profile. CONCLUSIONS: The results provide evidence of molecular similarities between AFX/PDS and undifferentiated cSCC and suggest a common histogenetic origin.

Low Energy Electron Irradiation Is a Potent Alternative to Gamma Irradiation for the Inactivation of (CAR-)NK-92 Cells in ATMP Manufacturing.

Background: With increasing clinical use of NK-92 cells and their CAR-modified derivatives in cancer immunotherapy, there is a growing demand for efficient production processes of these "off-the-shelf" therapeutics. In order to ensure safety and prevent the occurrence of secondary tumors, (CAR-)NK-92 cell proliferation has to be inactivated before transfusion. This is commonly achieved by gamma irradiation. Recently, we showed proof of concept that low energy electron irradiation (LEEI) is a new method for NK-92 inactivation. LEEI has several advantages over gamma irradiation, including a faster reaction time, a more reproducible dose rate and much less requirements on radiation shielding. Here, LEEI was further evaluated as a promising alternative to gamma irradiation yielding cells with highly maintained cytotoxic effector function. Methods: Effectiveness and efficiency of LEEI and gamma irradiation were analyzed using NK-92 and CD123-directed CAR-NK-92 cells. LEE-irradiated cells were extensively characterized and compared to gamma-irradiated cells via flow cytometry, cytotoxicity assays, and comet assays, amongst others. Results: Our results show that both irradiation methods caused a progressive decrease in cell viability and are, therefore, suitable for inhibition of cell proliferation. Notably, the NK-mediated specific lysis of tumor cells was maintained at stable levels for three days post-irradiation, with a trend towards higher activities after LEEI treatment as compared to gamma irradiation. Both gamma irradiation as well as LEEI led to substantial DNA damage and an accumulation of irradiated cells in the G2/M cell cycle phases. In addition, transcriptomic analysis of irradiated cells revealed approximately 12-fold more differentially expressed genes two hours after gamma irradiation, compared to LEEI. Analysis of surface molecules revealed an irradiation-induced decrease in surface expression of CD56, but no changes in the levels of the activating receptors NKp46, NKG2D, or NKp30. Conclusions: The presented data show that LEEI inactivates (CAR-)NK-92 cells as efficiently as gamma irradiation, but with less impact on the overall gene expression. Due to logistic advantages, LEEI might provide a superior alternative for the manufacture of (CAR-)NK-92 cells for clinical application.

Master and servant: LINC00152 - a STAT3-induced long noncoding RNA regulates STAT3 in a positive feedback in human multiple myeloma.

BACKGROUND: The survival of INA-6 human multiple myeloma cells is strictly dependent upon the Interleukin-6-activated transcription factor STAT3. Although transcriptional analyses have revealed many genes regulated by STAT3, to date no protein-coding STAT3 target gene is known to mediate survival in INA-6 cells. Therefore, the aim here was to identify and analyze non-protein-coding STAT3 target genes. In addition to the oncogenic microRNA-21, we previously described five long noncoding RNAs (lncRNAs) induced by STAT3, named STAiRs. Here, we focus on STAT3-induced RNA 18 (STAiR18), an mRNA-like, long ncRNA that is duplicated in the human lineage. One STAiR18 locus is annotated as the already well described LINC00152/CYTOR, however, the other harbors the MIR4435-2HG gene and is, up to now, barely described. METHODS: CAPTURE-RNA-sequencing was used to analyze STAiR18 transcript architecture. To identify the STAiR18 and STAT3 phenotype, siRNA-based knockdowns were performed and microarrays were applied to identify their target genes. RNA-binding partners of STAiR18 were determined by Chromatin-Isolation-by-RNA-Purification (ChIRP) and subsequent sequencing. STAT3 expression in dependence of STAiR18 was investigated by immunoblots, chromatin- and RNA-immunoprecipitations. RESULTS: As identified by CAPTURE-RNA sequencing, a complex splice pattern originates from both STAiR18 loci, generating different transcripts. Knockdown of the most abundant STAiR18 isoforms dramatically decreased INA-6 cell vitality, suggesting a functional role in myeloma cells. Additionally, STAiR18 and STAT3 knockdowns yielded overlapping changes of transcription patterns in INA-6 cells, suggesting a close functional interplay between the two factors. Moreover, Chromatin isolation by RNA purification (ChIRP), followed by genome-wide RNA sequencing showed that STAiR18 associates specifically with the STAT3 primary transcript. Furthermore, the knockdown of STAiR18 reduced STAT3 levels on both the RNA and protein levels, suggesting a positive feedback between both molecules. Furthermore, STAiR18 knockdown changes the histone methylation status of the STAT3 locus, which explains the positive feedback and indicates that STAiR18 is an epigenetic modulator. CONCLUSION: Hence, STAiR18 is an important regulator of myeloma cell survival and is strongly associated with the oncogenic function of STAT3. The close functional interplay between STAT3 and STAiR18 suggests a novel principle of regulatory interactions between long ncRNAs and signaling pathways.

uap: reproducible and robust HTS data analysis.

BACKGROUND: A lack of reproducibility has been repeatedly criticized in computational research. High throughput sequencing (HTS) data analysis is a complex multi-step process. For most of the steps a range of bioinformatic tools is available and for most tools manifold parameters need to be set. Due to this complexity, HTS data analysis is particularly prone to reproducibility and consistency issues. We have defined four criteria that in our opinion ensure a minimal degree of reproducible research for HTS data analysis. A series of workflow management systems is available for assisting complex multi-step data analyses. However, to the best of our knowledge, none of the currently available work flow management systems satisfies all four criteria for reproducible HTS analysis. RESULTS: Here we present uap, a workflow management system dedicated to robust, consistent, and reproducible HTS data analysis. uap is optimized for the application to omics data, but can be easily extended to other complex analyses. It is available under the GNU GPL v3 license at https://github.com/yigbt/uap. CONCLUSIONS: uap is a freely available tool that enables researchers to easily adhere to reproducible research principles for HTS data analyses.

STAT3-induced long noncoding RNAs in multiple myeloma cells display different properties in cancer.

Interleukin-6 (IL-6)-activated Signal Transducer and Activator of Transcription 3 (STAT3) facilitates survival in the multiple myeloma cell line INA-6 and therefore represents an oncogenic key player. However, the biological mechanisms are still not fully understood. In previous studies we identified microRNA-21 as a STAT3 target gene with strong anti-apoptotic potential, suggesting that noncoding RNAs have an impact on the pathogenesis of human multiple myeloma. Here, we describe five long noncoding RNAs (lncRNAs) induced by IL-6-activated STAT3, which we named STAiRs. While STAiRs 1, 2 and 6 remain unprocessed in the nucleus and show myeloma-specific expression, STAiRs 15 and 18 are spliced and broadly expressed. Especially STAiR2 and STAiR18 are promising candidates. STAiR2 originates from the first intron of a tumor suppressor gene. Our data support a mutually exclusive expression of either STAiR2 or the functional tumor suppressor in INA-6 cells and thus a contribution of STAiR2 to tumorigenesis. Furthermore, STAiR18 was shown to be overexpressed in every tested tumor entity, indicating its global role in tumor pathogenesis. Taken together, our study reveals a number of STAT3-induced lncRNAs suggesting that the interplay between the coding and noncoding worlds represents a fundamental principle of STAT3-driven cancer development in multiple myeloma and beyond.