Ten quick tips for building FAIR workflows.

Research data is accumulating rapidly and with it the challenge of fully reproducible science. As a consequence, implementation of high-quality management of scientific data has become a global priority. The FAIR (Findable, Accesible, Interoperable and Reusable) principles provide practical guidelines for maximizing the value of research data; however, processing data using workflows-systematic executions of a series of computational tools-is equally important for good data management. The FAIR principles have recently been adapted to Research Software (FAIR4RS Principles) to promote the reproducibility and reusability of any type of research software. Here, we propose a set of 10 quick tips, drafted by experienced workflow developers that will help researchers to apply FAIR4RS principles to workflows. The tips have been arranged according to the FAIR acronym, clarifying the purpose of each tip with respect to the FAIR4RS principles. Altogether, these tips can be seen as practical guidelines for workflow developers who aim to contribute to more reproducible and sustainable computational science, aiming to positively impact the open science and FAIR community.

Comprehensive diagnostics of acute myeloid leukemia by whole transcriptome RNA sequencing.

Acute myeloid leukemia (AML) is caused by genetic aberrations that also govern the prognosis of patients and guide risk-adapted and targeted therapy. Genetic aberrations in AML are structurally diverse and currently detected by different diagnostic assays. This study sought to establish whole transcriptome RNA sequencing as single, comprehensive, and flexible platform for AML diagnostics. We developed HAMLET (Human AML Expedited Transcriptomics) as bioinformatics pipeline for simultaneous detection of fusion genes, small variants, tandem duplications, and gene expression with all information assembled in an annotated, user-friendly output file. Whole transcriptome RNA sequencing was performed on 100 AML cases and HAMLET results were validated by reference assays and targeted resequencing. The data showed that HAMLET accurately detected all fusion genes and overexpression of EVI1 irrespective of 3q26 aberrations. In addition, small variants in 13 genes that are often mutated in AML were called with 99.2% sensitivity and 100% specificity, and tandem duplications in FLT3 and KMT2A were detected by a novel algorithm based on soft-clipped reads with 100% sensitivity and 97.1% specificity. In conclusion, HAMLET has the potential to provide accurate comprehensive diagnostic information relevant for AML classification, risk assessment and targeted therapy on a single technology platform.

Repeated FcεRI triggering reveals modified mast cell function related to chronic allergic responses in tissue.

BACKGROUND: Activation of mast cells through FcεRI plays an important role in acute allergic reactions. However, little is known about the function of mast cells in patients with chronic allergic inflammation or the effect of repeated FcεRI triggering occurring in such responses. OBJECTIVE: We aimed to identify changes in mast cell function after repeated FcεRI triggering and to correlate these changes to chronic allergic responses in tissue. METHODS: Human cord blood-derived mast cells were treated for 2 weeks with anti-IgE. The function of naive or treated mast cells was analyzed by means of RNA sequencing, quantitative RT-PCR, flow cytometry, and functional assays. Protein secretion was measured with ELISAs and multiplex assays. RESULTS: We observed several changes in mast cell function after repeated anti-IgE triggering. Although the acute response was dampened, we identified 289 genes significantly upregulated after repeated anti-IgE. Most of these genes (84%) were not upregulated after a single anti-IgE stimulus, indicating a significantly different response mode characterized by increased antigen presentation, response to bacteria, and chemotaxis. Changes in mast cell function were related to changes in expression of the transcription factors RXRA and BATF and others. Importantly, we found a substantial overlap between genes upregulated after repeated anti-IgE triggering and genes upregulated in tissue from patients with chronic allergy, in particular those of patients with chronic rhinosinusitis. CONCLUSION: Our study provides evidence for intrinsic modulation of mast cell function on repeated FcεRI-mediated activation. The overlap with gene expression in tissues is suggestive of a direct link between repeated IgE-mediated activation of mast cells and chronic allergy.

Integrated Whole Genome and Transcriptome Analysis Identified a Therapeutic Minor Histocompatibility Antigen in a Splice Variant of ITGB2.

PURPOSE: In HLA-matched allogeneic hematopoietic stem cell transplantation (alloSCT), donor T cells recognizing minor histocompatibility antigens (MiHAs) can mediate desired antitumor immunity as well as undesired side effects. MiHAs with hematopoiesis-restricted expression are relevant targets to augment antitumor immunity after alloSCT without side effects. To identify therapeutic MiHAs, we analyzed the in vivo immune response in a patient with strong antitumor immunity after alloSCT. EXPERIMENTAL DESIGN: T-cell clones recognizing patient, but not donor, hematopoietic cells were selected for MiHA discovery by whole genome association scanning. RNA-sequence data from the GEUVADIS project were analyzed to investigate alternative transcripts, and expression patterns were determined by microarray analysis and qPCR. T-cell reactivity was measured by cytokine release and cytotoxicity. RESULTS: T-cell clones were isolated for two HLA-B*15:01-restricted MiHA. LB-GLE1-1V is encoded by a nonsynonymous SNP in exon 6 of GLE1 For the other MiHAs, an associating SNP in intron 3 of ITGB2 was found, but no SNP disparity was present in the normal gene transcript between patient and donor. RNA-sequence analysis identified an alternative ITGB2 transcript containing part of intron 3. qPCR demonstrated that this transcript is restricted to hematopoietic cells and SNP-positive individuals. In silico translation revealed LB-ITGB2-1 as HLA-B*15:01-binding peptide, which was validated as hematopoietic MiHA by T-cell experiments. CONCLUSIONS: Whole genome and transcriptome analysis identified LB-ITGB2-1 as MiHAs encoded by an alternative transcript. Our data support the therapeutic relevance of LB-ITGB2-1 and illustrate the value of RNA-sequence analysis for discovery of immune targets encoded by alternative transcripts. Clin Cancer Res; 22(16); 4185-96. ©2016 AACR.

Tumor cell migration screen identifies SRPK1 as breast cancer metastasis determinant.

Tumor cell migration is a key process for cancer cell dissemination and metastasis that is controlled by signal-mediated cytoskeletal and cell matrix adhesion remodeling. Using a phagokinetic track assay with migratory H1299 cells, we performed an siRNA screen of almost 1,500 genes encoding kinases/phosphatases and adhesome- and migration-related proteins to identify genes that affect tumor cell migration speed and persistence. Thirty candidate genes that altered cell migration were validated in live tumor cell migration assays. Eight were associated with metastasis-free survival in breast cancer patients, with integrin ß3-binding protein (ITGB3BP), MAP3K8, NIMA-related kinase (NEK2), and SHC-transforming protein 1 (SHC1) being the most predictive. Examination of genes that modulate migration indicated that SRPK1, encoding the splicing factor kinase SRSF protein kinase 1, is relevant to breast cancer outcomes, as it was highly expressed in basal breast cancer. Furthermore, high SRPK1 expression correlated with poor breast cancer disease outcome and preferential metastasis to the lungs and brain. In 2 independent murine models of breast tumor metastasis, stable shRNA-based SRPK1 knockdown suppressed metastasis to distant organs, including lung, liver, and spleen, and inhibited focal adhesion reorganization. Our study provides comprehensive information on the molecular determinants of tumor cell migration and suggests that SRPK1 has potential as a drug target for limiting breast cancer metastasis.

RNA expression profiling in brains of familial hemiplegic migraine type 1 knock-in mice.

BACKGROUND: Various CACNA1A missense mutations cause familial hemiplegic migraine type 1 (FHM1), a rare monogenic subtype of migraine with aura. FHM1 mutation R192Q is associated with pure hemiplegic migraine, whereas the S218L mutation causes hemiplegic migraine, cerebellar ataxia, seizures, and mild head trauma-induced brain edema. Transgenic knock-in (KI) migraine mouse models were generated that carried either the FHM1 R192Q or the S218L mutation and were shown to exhibit increased CaV2.1 channel activity. Here we investigated their cerebellar and caudal cortical transcriptome. METHODS: Caudal cortical and cerebellar RNA expression profiles from mutant and wild-type mice were studied using microarrays. Respective brain regions were selected based on their relevance to migraine aura and ataxia. Relevant expression changes were further investigated at RNA and protein level by quantitative polymerase chain reaction (qPCR) and/or immunohistochemistry, respectively. RESULTS: Expression differences in the cerebellum were most pronounced in S218L mice. Particularly, tyrosine hydroxylase, a marker of delayed cerebellar maturation, appeared strongly upregulated in S218L cerebella. In contrast, only minimal expression differences were observed in the caudal cortex of either mutant mice strain. CONCLUSION: Despite pronounced consequences of migraine gene mutations at the neurobiological level, changes in cortical RNA expression in FHM1 migraine mice compared to wild-type are modest. In contrast, pronounced RNA expression changes are seen in the cerebellum of S218L mice and may explain their cerebellar ataxia phenotype.

Systematic identification of trans eQTLs as putative drivers of known disease associations.

Identifying the downstream effects of disease-associated SNPs is challenging. To help overcome this problem, we performed expression quantitative trait locus (eQTL) meta-analysis in non-transformed peripheral blood samples from 5,311 individuals with replication in 2,775 individuals. We identified and replicated trans eQTLs for 233 SNPs (reflecting 103 independent loci) that were previously associated with complex traits at genome-wide significance. Some of these SNPs affect multiple genes in trans that are known to be altered in individuals with disease: rs4917014, previously associated with systemic lupus erythematosus (SLE), altered gene expression of C1QB and five type I interferon response genes, both hallmarks of SLE. DeepSAGE RNA sequencing showed that rs4917014 strongly alters the 3' UTR levels of IKZF1 in cis, and chromatin immunoprecipitation and sequencing analysis of the trans-regulated genes implicated IKZF1 as the causal gene. Variants associated with cholesterol metabolism and type 1 diabetes showed similar phenomena, indicating that large-scale eQTL mapping provides insight into the downstream effects of many trait-associated variants.