A novel entity of massive multifocal osteolyses in the elderly.

Osteolyses are common findings in elderly patients and most frequently represent malignant or locally aggressive bone tumors, infection, inflammatory and endocrine disorders, histiocytoses, and rare diseases such as Gorham-Stout syndrome. We here report on a novel entity of massive multifocal osteolyses in both shoulders, the right hip and left knee joint and the dens of an 83-year-old patient not relatable to any previously known etiopathology of bone disorders. The soft tissue mass is of myxoid stroma with an unspecific granulomatous inflammatory process, aggressively destroying extensive cortical and cancellous bone segments and encroaching on articulating bones in diarthrodial large joints. Radiological, nuclear medical, serological, histological, and immunohistochemical analyses were incapable of further classifying the disease pattern within the existing scheme of pathology. Quantitative polymerase chain reaction and next generation sequencing revealed that mutations are not suggestive of any known hereditary or acquired bone disease. Possible treatment options include radionuclide therapy for pain palliation and percutaneous radiation to arrest bone resorption while surgical treatment is inevitable for pathological fractures. This case study shall increase the awareness of the musculoskeletal community and motivate to collect further information on this rare but mutilating disorder.

Cas9-mediated knockout of Ndrg2 enhances the regenerative potential of dendritic cells for wound healing.

Chronic wounds impose a significant healthcare burden to a broad patient population. Cell-based therapies, while having shown benefits for the treatment of chronic wounds, have not yet achieved widespread adoption into clinical practice. We developed a CRISPR/Cas9 approach to precisely edit murine dendritic cells to enhance their therapeutic potential for healing chronic wounds. Using single-cell RNA sequencing of tolerogenic dendritic cells, we identified N-myc downregulated gene 2 (Ndrg2), which marks a specific population of dendritic cell progenitors, as a promising target for CRISPR knockout. Ndrg2-knockout alters the transcriptomic profile of dendritic cells and preserves an immature cell state with a strong pro-angiogenic and regenerative capacity. We then incorporated our CRISPR-based cell engineering within a therapeutic hydrogel for in vivo cell delivery and developed an effective translational approach for dendritic cell-based immunotherapy that accelerated healing of full-thickness wounds in both non-diabetic and diabetic mouse models. These findings could open the door to future clinical trials using safe gene editing in dendritic cells for treating various types of chronic wounds.

Time-resolved RNA signatures of CD4+ T cells in Parkinson's disease.

Parkinson's disease (PD) emerges as a complex, multifactorial disease. While there is increasing evidence that dysregulated T cells play a central role in PD pathogenesis, elucidation of the pathomechanical changes in related signaling is still in its beginnings. We employed time-resolved RNA expression upon the activation of peripheral CD4+ T cells to track and functionally relate changes on cellular signaling in representative cases of patients at different stages of PD. While only few miRNAs showed time-course related expression changes in PD, we identified groups of genes with significantly altered expression for each different time window. Towards a further understanding of the functional consequences, we highlighted pathways with decreased or increased activity in PD, including the most prominent altered IL-17 pathway. Flow cytometric analyses showed not only an increased prevalence of Th17 cells but also a specific subtype of IL-17 producing γδ-T cells, indicating a previously unknown role in PD pathogenesis.

Changes of Protein Expression after CRISPR/Cas9 Knockout of miRNA-142 in Cell Lines Derived from Diffuse Large B-Cell Lymphoma.

BACKGROUND: As microRNA-142 (miR-142) is the only human microRNA gene where mutations have consistently been found in about 20% of all cases of diffuse large B-cell lymphoma (DLBCL), we wanted to determine the impact of miR-142 inactivation on protein expression of DLBCL cell lines. METHODS: miR-142 was deleted by CRISPR/Cas9 knockout in cell lines from DLBCL. RESULTS: By proteome analyses, miR-142 knockout resulted in a consistent up-regulation of 52 but also down-regulation of 41 proteins in GC-DLBCL lines BJAB and SUDHL4. Various mitochondrial ribosomal proteins were up-regulated in line with their pro-tumorigenic properties, while proteins necessary for MHC-I presentation were down-regulated in accordance with the finding that miR-142 knockout mice have a defective immune response. CFL2, CLIC4, STAU1, and TWF1 are known targets of miR-142, and we could additionally confirm AKT1S1, CCNB1, LIMA1, and TFRC as new targets of miR-142-3p or -5p. CONCLUSIONS: Seed-sequence mutants of miR-142 confirmed potential targets and novel targets of miRNAs can be identified in miRNA knockout cell lines. Due to the complex contribution of miRNAs within cellular regulatory networks, in particular when miRNAs highly present in RISC complexes are replaced by other miRNAs, primary effects on gene expression may be covered by secondary layers of regulation.

Simultaneous regression and classification for drug sensitivity prediction using an advanced random forest method.

Machine learning methods trained on cancer cell line panels are intensively studied for the prediction of optimal anti-cancer therapies. While classification approaches distinguish effective from ineffective drugs, regression approaches aim to quantify the degree of drug effectiveness. However, the high specificity of most anti-cancer drugs induces a skewed distribution of drug response values in favor of the more drug-resistant cell lines, negatively affecting the classification performance (class imbalance) and regression performance (regression imbalance) for the sensitive cell lines. Here, we present a novel approach called SimultAneoUs Regression and classificatiON Random Forests (SAURON-RF) based on the idea of performing a joint regression and classification analysis. We demonstrate that SAURON-RF improves the classification and regression performance for the sensitive cell lines at the expense of a moderate loss for the resistant ones. Furthermore, our results show that simultaneous classification and regression can be superior to regression or classification alone.

The DGCR8 E518K mutation found in Wilms tumors leads to a partial miRNA processing defect that alters gene expression patterns and biological processes.

Wilms tumor (WT) is the most common renal tumor in childhood. We and others have previously identified oncogenic driver mutations affecting the microprocessor genes DROSHA and DGCR8 that lead to altered miRNA expression patterns. In the case of DGCR8, a single recurrent hotspot mutation (E518K) was found in the RNA binding domain. To functionally assess this mutation in vitro, we generated mouse Dgcr8-KO embryonic stem cell (mESC) lines with an inducible expression of wild-type or mutant DGCR8, mirroring the hemizygous mutant expression seen in WT. RNA-seq analysis revealed significant differences of miRNA expression profiles in DGCR8-E518K compared with DGCR8-wild-type mESCs. The E518K mutation only led to a partial rescue of the reported miRNA processing defect in Dgcr8-KO, with selectively reduced expression of numerous canonical miRNAs. Nevertheless, DGCR8-E518K retained significant activity given its ability to still process many miRNAs. Subsequent to altered miRNA levels, the expression of mRNA targets was likewise changed. Functional assays showed that DGCR8-E518K cells still have a partial proliferation and differentiation defect but were able to rescue critical biological processes in embryoid body development. The stem cell program could be shut down and all three germ layers were formed. These findings suggest that the E518K mutation leads to a partial reduction of microprocessor activity and altered specificity with selective impairment only in certain developmental contexts, apparently including nephrogenesis.

GeneTrail: A Framework for the Analysis of High-Throughput Profiles.

Experimental high-throughput techniques, like next-generation sequencing or microarrays, are nowadays routinely applied to create detailed molecular profiles of cells. In general, these platforms generate high-dimensional and noisy data sets. For their analysis, powerful bioinformatics tools are required to gain novel insights into the biological processes under investigation. Here, we present an overview of the GeneTrail tool suite that offers rich functionality for the analysis and visualization of (epi-)genomic, transcriptomic, miRNomic, and proteomic profiles. Our framework enables the analysis of standard bulk, time-series, and single-cell measurements and includes various state-of-the-art methods to identify potentially deregulated biological processes and to detect driving factors within those deregulated processes. We highlight the capabilities of our web service with an analysis of a single-cell COVID-19 data set that demonstrates its potential for uncovering complex molecular mechanisms. GeneTrail can be accessed freely and without login requirements at http://genetrail.bioinf.uni-sb.de.

MERIDA: a novel Boolean logic-based integer linear program for personalized cancer therapy.

MOTIVATION: A major goal of personalized medicine in oncology is the optimization of treatment strategies given measurements of the genetic and molecular profiles of cancer cells. To further our knowledge on drug sensitivity, machine learning techniques are commonly applied to cancer cell line panels. RESULTS: We present a novel integer linear programming formulation, called MEthod for Rule Identification with multi-omics DAta (MERIDA), for predicting the drug sensitivity of cancer cells. The method represents a modified version of the LOBICO method and yields easily interpretable models amenable to a Boolean logic-based interpretation. Since the proposed altered logical rules lead to an enormous acceleration of the running times of MERIDA compared to LOBICO, we cannot only consider larger input feature sets integrated from genetic and molecular omics data but also build more comprehensive models that mirror the complexity of cancer initiation and progression. Moreover, we enable the inclusion of a priori knowledge that can either stem from biomarker databases or can also be newly acquired knowledge gathered iteratively by previous runs of MERIDA. Our results show that this approach does not only lead to an improved predictive performance but also identifies a variety of putative sensitivity and resistance biomarkers. We also compare our approach to state-of-the-art machine learning methods and demonstrate the superior performance of our method. Hence, MERIDA has great potential to deepen our understanding of the molecular mechanisms causing drug sensitivity or resistance. AVAILABILITY AND IMPLEMENTATION: The corresponding code is available on github (https://github.com/unisb-bioinf/MERIDA.git). SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

miRTargetLink 2.0-interactive miRNA target gene and target pathway networks.

Which genes, gene sets or pathways are regulated by certain miRNAs? Which miRNAs regulate a particular target gene or target pathway in a certain physiological context? Answering such common research questions can be time consuming and labor intensive. Especially for researchers without computational experience, the integration of different data sources, selection of the right parameters and concise visualization can be demanding. A comprehensive analysis should be central to present adequate answers to complex biological questions. With miRTargetLink 2.0, we develop an all-in-one solution for human, mouse and rat miRNA networks. Users input in the unidirectional search mode either a single gene, gene set or gene pathway, alternatively a single miRNA, a set of miRNAs or an miRNA pathway. Moreover, genes and miRNAs can jointly be provided to the tool in the bidirectional search mode. For the selected entities, interaction graphs are generated from different data sources and dynamically presented. Connected application programming interfaces (APIs) to the tailored enrichment tools miEAA and GeneTrail facilitate downstream analysis of pathways and context-annotated categories of network nodes. MiRTargetLink 2.0 is freely accessible at https://www.ccb.uni-saarland.de/mirtargetlink2.

Validation of human microRNA target pathways enables evaluation of target prediction tools.

MicroRNAs are regulators of gene expression. A wide-spread, yet not validated, assumption is that the targetome of miRNAs is non-randomly distributed across the transcriptome and that targets share functional pathways. We developed a computational and experimental strategy termed high-throughput miRNA interaction reporter assay (HiTmIR) to facilitate the validation of target pathways. First, targets and target pathways are predicted and prioritized by computational means to increase the specificity and positive predictive value. Second, the novel webtool miRTaH facilitates guided designs of reporter assay constructs at scale. Third, automated and standardized reporter assays are performed. We evaluated HiTmIR using miR-34a-5p, for which TNF- and TGFB-signaling, and Parkinson's Disease (PD)-related categories were identified and repeated the pipeline for miR-7-5p. HiTmIR validated 58.9% of the target genes for miR-34a-5p and 46.7% for miR-7-5p. We confirmed the targeting by measuring the endogenous protein levels of targets in a neuronal cell model. The standardized positive and negative targets are collected in the new miRATBase database, representing a resource for training, or benchmarking new target predictors. Applied to 88 target predictors with different confidence scores, TargetScan 7.2 and miRanda outperformed other tools. Our experiments demonstrate the efficiency of HiTmIR and provide evidence for an orchestrated miRNA-gene targeting.

Wrinkle in the plan: miR-34a-5p impacts chemokine signaling by modulating CXCL10/CXCL11/CXCR3-axis in CD4+, CD8+ T cells, and M1 macrophages.

BACKGROUND: In 2016 the first-in-human phase I study of a miRNA-based cancer therapy with a liposomal mimic of microRNA-34a-5p (miR-34a-5p) was closed due to five immune related serious adverse events (SAEs) resulting in four patient deaths. For future applications of miRNA mimics in cancer therapy it is mandatory to unravel the miRNA effects both on the tumor tissue and on immune cells. Here, we set out to analyze the impact of miR-34a-5p over-expression on the CXCL10/CXCL11/CXCR3 axis, which is central for the development of an effective cancer control. METHODS: We performed a whole genome expression analysis of miR-34a-5p transfected M1 macrophages followed by an over-representation and a protein-protein network analysis. In-silico miRNA target prediction and dual luciferase assays were used for target identification and verification. Target genes involved in chemokine signaling were functionally analyzed in M1 macrophages, CD4+ and CD8+ T cells. RESULTS: A whole genome expression analysis of M1 macrophages with induced miR-34a-5p over-expression revealed an interaction network of downregulated target mRNAs including CXCL10 and CXCL11. In-silico target prediction in combination with dual luciferase assays identified direct binding of miR-34a-5p to the 3'UTRs of CXCL10 and CXCL11. Decreased CXCL10 and CXCL11 secretion was shown on the endogenous protein level and in the supernatant of miR-34a-5p transfected and activated M1 macrophages. To complete the analysis of the CXCL10/CXCL11/CXCR3 axis, we activated miR-34a-5p transfected CD4+ and CD8+ T cells by PMA/Ionomycin and found reduced levels of endogenous CXCR3 and CXCR3 on the cell surface. CONCLUSIONS: MiR-34a-5p mimic administered by intravenous administration will likely not only be up-taken by the tumor cells but also by the immune cells. Our results indicate that miR-34a-5p over-expression leads in M1 macrophages to a reduced secretion of CXCL10 and CXCL11 chemokines and in CD4+ and CD8+ T cells to a reduced expression of CXCR3. As a result, less immune cells will be attracted to the tumor site. Furthermore, high levels of miR-34a-5p in naive CD4+ T cells can in turn hinder Th1 cell polarization through the downregulation of CXCR3 leading to a less pronounced activation of cytotoxic T lymphocytes, natural killer, and natural killer T cells and possibly contributing to lymphocytopenia.

Quantitative and time-resolved miRNA pattern of early human T cell activation.

T cells are central to the immune response against various pathogens and cancer cells. Complex networks of transcriptional and post-transcriptional regulators, including microRNAs (miRNAs), coordinate the T cell activation process. Available miRNA datasets, however, do not sufficiently dissolve the dynamic changes of miRNA controlled networks upon T cell activation. Here, we established a quantitative and time-resolved expression pattern for the entire miRNome over a period of 24 h upon human T-cell activation. Based on our time-resolved datasets, we identified central miRNAs and specified common miRNA expression profiles. We found the most prominent quantitative expression changes for miR-155-5p with a range from initially 40 molecules/cell to 1600 molecules/cell upon T-cell activation. We established a comprehensive dynamic regulatory network of both the up- and downstream regulation of miR-155. Upstream, we highlight IRF4 and its complexes with SPI1 and BATF as central for the transcriptional regulation of miR-155. Downstream of miR-155-5p, we verified 17 of its target genes by the time-resolved data recorded after T cell activation. Our data provide comprehensive insights into the range of stimulus induced miRNA abundance changes and lay the ground to identify efficient points of intervention for modifying the T cell response.

GeneTrail 3: advanced high-throughput enrichment analysis.

We present GeneTrail 3, a major extension of our web service GeneTrail that offers rich functionality for the identification, analysis, and visualization of deregulated biological processes. Our web service provides a comprehensive collection of biological processes and signaling pathways for 12 model organisms that can be analyzed with a powerful framework for enrichment and network analysis of transcriptomic, miRNomic, proteomic, and genomic data sets. Moreover, GeneTrail offers novel workflows for the analysis of epigenetic marks, time series experiments, and single cell data. We demonstrate the capabilities of our web service in two case-studies, which highlight that GeneTrail is well equipped for uncovering complex molecular mechanisms. GeneTrail is freely accessible at: http://genetrail.bioinf.uni-sb.de.

miRPathDB 2.0: a novel release of the miRNA Pathway Dictionary Database.

Since the initial release of miRPathDB, tremendous progress has been made in the field of microRNA (miRNA) research. New miRNA reference databases have emerged, a vast amount of new miRNA candidates has been discovered and the number of experimentally validated target genes has increased considerably. Hence, the demand for a major upgrade of miRPathDB, including extended analysis functionality and intuitive visualizations of query results has emerged. Here, we present the novel release 2.0 of the miRNA Pathway Dictionary Database (miRPathDB) that is freely accessible at https://mpd.bioinf.uni-sb.de/. miRPathDB 2.0 comes with a ten-fold increase of pre-processed data. In total, the updated database provides putative associations between 27 452 (candidate) miRNAs, 28 352 targets and 16 833 pathways for Homo sapiens, as well as interactions of 1978 miRNAs, 24 898 targets and 6511 functional categories for Mus musculus. Additionally, we analyzed publications citing miRPathDB to identify common use-cases and further extensions. Based on this evaluation, we added new functionality for interactive visualizations and down-stream analyses of bulk queries. In summary, the updated version of miRPathDB, with its new custom-tailored features, is one of the most comprehensive and advanced resources for miRNAs and their target pathways.

miR-34a as hub of T cell regulation networks.

BACKGROUND: Micro(mi)RNAs are increasingly recognized as central regulators of immune cell function. While it has been predicted that miRNAs have multiple targets, the majority of these predictions still await experimental confirmation. Here, miR-34a, a well-known tumor suppressor, is analyzed for targeting genes involved in immune system processes of leucocytes. METHODS: Using an in-silico approach, we combined miRNA target prediction with GeneTrail2, a web tool for Multi-omics enrichment analysis, to identify miR-34a target genes, which are involved in the immune system process subcategory of Gene Ontology. RESULTS: Out of the 193 predicted target genes in this subcategory we experimentally tested 22 target genes and confirmed binding of miR-34a to 14 target genes including VAMP2, IKBKE, MYH9, MARCH8, KLRK1, CD11A, TRAFD1, CCR1, PYDC1, PRF1, PIK3R2, PIK3CD, AP1B1, and ADAM10 by dual luciferase assays. By transfecting Jurkat, primary CD4+ and CD8+ T cells with miR-34a, we demonstrated that ectopic expression of miR-34a leads to reduced levels of endogenous VAMP2 and CD11A, which are central to the analyzed subcategories. Functional downstream analysis of miR-34a over-expression in activated CD8+ T cells exhibits a distinct decrease of PRF1 secretion. CONCLUSIONS: By simultaneous targeting of 14 mRNAs miR-34a acts as major hub of T cell regulatory networks suggesting to utilize miR-34a as target of intervention towards a modulation of the immune responsiveness of T-cells in a broad tumor context.

ClinOmicsTrailbc: a visual analytics tool for breast cancer treatment stratification.

MOTIVATION: Breast cancer is the second leading cause of cancer death among women. Tumors, even of the same histopathological subtype, exhibit a high genotypic diversity that impedes therapy stratification and that hence must be accounted for in the treatment decision-making process. RESULTS: Here, we present ClinOmicsTrailbc, a comprehensive visual analytics tool for breast cancer decision support that provides a holistic assessment of standard-of-care targeted drugs, candidates for drug repositioning and immunotherapeutic approaches. To this end, our tool analyzes and visualizes clinical markers and (epi-)genomics and transcriptomics datasets to identify and evaluate the tumor's main driver mutations, the tumor mutational burden, activity patterns of core cancer-relevant pathways, drug-specific biomarkers, the status of molecular drug targets and pharmacogenomic influences. In order to demonstrate ClinOmicsTrailbc's rich functionality, we present three case studies highlighting various ways in which ClinOmicsTrailbc can support breast cancer precision medicine. ClinOmicsTrailbc is a powerful integrated visual analytics tool for breast cancer research in general and for therapy stratification in particular, assisting oncologists to find the best possible treatment options for their breast cancer patients based on actionable, evidence-based results. AVAILABILITY AND IMPLEMENTATION: ClinOmicsTrailbc can be freely accessed at https://clinomicstrail.bioinf.uni-sb.de. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

The role of TCF3 as potential master regulator in blastemal Wilms tumors.

Wilms tumors are the most common type of pediatric kidney tumors. While the overall prognosis for patients is favorable, especially tumors that exhibit a blastemal subtype after preoperative chemotherapy have a poor prognosis. For an improved risk assessment and therapy stratification, it is essential to identify the driving factors that are distinctive for this aggressive subtype. In our study, we compared gene expression profiles of 33 tumor biopsies (17 blastemal and 16 other tumors) after neoadjuvant chemotherapy. The analysis of this dataset using the Regulator Gene Association Enrichment algorithm successfully identified several biomarkers and associated molecular mechanisms that distinguish between blastemal and nonblastemal Wilms tumors. Specifically, regulators involved in embryonic development and epigenetic processes like chromatin remodeling and histone modification play an essential role in blastemal tumors. In this context, we especially identified TCF3 as the central regulatory element. Furthermore, the comparison of ChIP-Seq data of Wilms tumor cell cultures from a blastemal mouse xenograft and a stromal tumor provided further evidence that the chromatin states of blastemal cells share characteristics with embryonic stem cells that are not present in the stromal tumor cell line. These stem-cell like characteristics could potentially add to the increased malignancy and chemoresistance of the blastemal subtype. Along with TCF3, we detected several additional biomarkers that are distinctive for blastemal Wilms tumors after neoadjuvant chemotherapy and that may provide leads for new therapeutic regimens.

Large-scale validation of miRNAs by disease association, evolutionary conservation and pathway activity.

The validation of microRNAs (miRNAs) identified by next generation sequencing involves amplification-free and hybridization-based detection of transcripts as criteria for confirming valid miRNAs. Since respective validation is frequently not performed, miRNA repositories likely still contain a substantial fraction of false positive candidates while true miRNAs are not stored in the repositories yet. Especially if downstream analyses are performed with these candidates (e.g. target or pathway prediction), the results may be misleading. In the present study, we evaluated 558 mature miRNAs from miRBase and 1,709 miRNA candidates from next generation sequencing experiments by amplification-free hybridization and investigated their distributions in patients with various disease conditions. Notably, the most significant miRNAs in diseases are often not contained in the miRBase. However, these candidates are evolutionary highly conserved. From the expression patterns, target gene and pathway analyses and evolutionary conservation analyses, we were able to shed light on the complexity of miRNAs in humans. Our data also highlight that a more thorough validation of miRNAs identified by next generation sequencing is required. The results are available in miRCarta ( https://mircarta.cs.uni-saarland.de ).

REGGAE: a novel approach for the identification of key transcriptional regulators.

Motivation: Transcriptional regulators play a major role in most biological processes. Alterations in their activities are associated with a variety of diseases and in particular with tumor development and progression. Hence, it is important to assess the effects of deregulated regulators on pathological processes. Results: Here, we present REGulator-Gene Association Enrichment (REGGAE), a novel method for the identification of key transcriptional regulators that have a significant effect on the expression of a given set of genes, e.g. genes that are differentially expressed between two sample groups. REGGAE uses a Kolmogorov-Smirnov-like test statistic that implicitly combines associations between regulators and their target genes with an enrichment approach to prioritize the influence of transcriptional regulators. We evaluated our method in two different application scenarios, which demonstrate that REGGAE is well suited for uncovering the influence of transcriptional regulators and is a valuable tool for the elucidation of complex regulatory mechanisms. Availability and implementation: REGGAE is freely available at https://regulatortrail.bioinf.uni-sb.de. Supplementary information: Supplementary data are available at Bioinformatics online.

miRCarta: a central repository for collecting miRNA candidates.

The continuous increase of available biological data as consequence of modern high-throughput technologies poses new challenges for analysis techniques and database applications. Especially for miRNAs, one class of small non-coding RNAs, many algorithms have been developed to predict new candidates from next-generation sequencing data. While the amount of publications describing novel miRNA candidates keeps steadily increasing, the current gold standard database for miRNAs - miRBase - has not been updated since June 2014. As a result, publications describing new miRNA candidates in the last three to five years might have a substantial overlap of candidates without noticing. With miRCarta we implemented a database to collect novel miRNA candidates and augment the information provided by miRBase. In the first stage, miRCarta is thought to be a highly sensitive collection of potential miRNA candidates with a high degree of analysis functionality, annotations and details on each miRNA. We added-besides the full content of the miRBase-12,857 human miRNA precursors to miRCarta. Users can match their own predictions to the entries of miRCarta to reduce potential redundancies in their studies. miRCarta provides the most comprehensive collection of human miRNAs and miRNA candidates to form a basis for further refinement and validation studies. The database is freely accessible at https://mircarta.cs.uni-saarland.de/.