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.

Gene expression profile of human T cells following a single stimulation of peripheral blood mononuclear cells with anti-CD3 antibodies.

BACKGROUND: Anti-CD3 immunotherapy was initially approved for clinical use for renal transplantation rejection prevention. Subsequently, new generations of anti-CD3 antibodies have entered clinical trials for a broader spectrum of therapeutic applications, including cancer and autoimmune diseases. Despite their extensive use, little is known about the exact mechanism of these molecules, except that they are able to activate T cells, inducing an overall immunoregulatory and tolerogenic behavior. To better understand the effects of anti-CD3 antibodies on human T cells, PBMCs were stimulated, and then, we performed RNA-seq assays of enriched T cells to assess changes in their gene expression profiles. In this study, three different anti-CD3 antibodies were used for the stimulation: two recombinant antibody fragments, namely, a humanized and a chimeric FvFc molecule, and the prototype mouse mAb OKT3. RESULTS: Gene Ontology categories and individual immunoregulatory markers were compared, suggesting a similarity in modulated gene sets, mainly those for immunoregulatory and inflammatory terms. Upregulation of interleukin receptors, such as IL2RA, IL1R, IL12RB2, IL18R1, IL21R and IL23R, and of inhibitory molecules, such as FOXP3, CTLA4, TNFRSF18, LAG3 and PDCD1, were also observed, suggesting an inhibitory and exhausted phenotype. CONCLUSIONS: We used a deep transcriptome sequencing method for comparing three anti-CD3 antibodies in terms of Gene Ontology enrichment and immunological marker expression. The present data showed that both recombinant antibodies induced a compatible expression profile, suggesting that they might be candidates for a closer evaluation with respect to their therapeutic value. Moreover, the proposed methodology is amenable to be more generally applied for molecular comparison of cell receptor dependent antibody therapy.

DIEGO: detection of differential alternative splicing using Aitchison's geometry.

Motivation: Alternative splicing is a biological process of fundamental importance in most eukaryotes. It plays a pivotal role in cell differentiation and gene regulation and has been associated with a number of different diseases. The widespread availability of RNA-Sequencing capacities allows an ever closer investigation of differentially expressed isoforms. However, most tools for differential alternative splicing (DAS) analysis do not take split reads, i.e. the most direct evidence for a splice event, into account. Here, we present DIEGO, a compositional data analysis method able to detect DAS between two sets of RNA-Seq samples based on split reads. Results: The python tool DIEGO works without isoform annotations and is fast enough to analyze large experiments while being robust and accurate. We provide python and perl parsers for common formats. Availability and implementation: The software is available at: www.bioinf.uni-leipzig.de/Software/DIEGO. Contact: steve@bioinf.uni-leipzig.de. Supplementary information: Supplementary data are available at Bioinformatics online.

MINCR is a MYC-induced lncRNA able to modulate MYC's transcriptional network in Burkitt lymphoma cells.

Despite the established role of the transcription factor MYC in cancer, little is known about the impact of a new class of transcriptional regulators, the long noncoding RNAs (lncRNAs), on MYC ability to influence the cellular transcriptome. Here, we have intersected RNA-sequencing data from two MYC-inducible cell lines and a cohort of 91 B-cell lymphomas with or without genetic variants resulting in MYC overexpression. We identified 13 lncRNAs differentially expressed in IG-MYC-positive Burkitt lymphoma and regulated in the same direction by MYC in the model cell lines. Among them, we focused on a lncRNA that we named MYC-induced long noncoding RNA (MINCR), showing a strong correlation with MYC expression in MYC-positive lymphomas. To understand its cellular role, we performed RNAi and found that MINCR knockdown is associated with an impairment in cell cycle progression. Differential gene expression analysis after RNAi showed a significant enrichment of cell cycle genes among the genes down-regulated after MINCR knockdown. Interestingly, these genes are enriched in MYC binding sites in their promoters, suggesting that MINCR acts as a modulator of the MYC transcriptional program. Accordingly, MINCR knockdown was associated with a reduction in MYC binding to the promoters of selected cell cycle genes. Finally, we show that down-regulation of Aurora kinases A and B and chromatin licensing and DNA replication factor 1 may explain the reduction in cellular proliferation observed on MINCR knockdown. We, therefore, suggest that MINCR is a newly identified player in the MYC transcriptional network able to control the expression of cell cycle genes.

Atypical RNAs in the coelacanth transcriptome.

Circular and apparently trans-spliced RNAs have recently been reported as abundant types of transcripts in mammalian transcriptome data. Both types of non-colinear RNAs are also abundant in RNA-seq of different tissue from both the African and the Indonesian coelacanth. We observe more than 8,000 lincRNAs with normal gene structure and several thousands of circularized and trans-spliced products, showing that such atypical RNAs form a substantial contribution to the transcriptome. Surprisingly, the majority of the circularizing and trans-connecting splice junctions are unique to atypical forms, that is, are not used in normal isoforms.

Bayesian sampling of evolutionarily conserved RNA secondary structures with pseudoknots.

MOTIVATION: Today many non-coding RNAs are known to play an active role in various important biological processes. Since RNA's functionality is correlated with specific structural motifs that are often conserved in phylogenetically related molecules, computational prediction of RNA structure should ideally be based on a set of homologous primary structures. But many available RNA secondary structure prediction programs that use sequence alignments do not consider pseudoknots or their estimations consist on a single structure without information on uncertainty. RESULTS: In this article we present a method that takes advantage of the evolutionary history of a group of aligned RNA sequences for sampling consensus secondary structures, including pseudoknots, according to their approximate posterior probability. We investigate the benefit of using evolutionary history and demonstrate the competitiveness of our method compared with similar methods based on RNase P RNA sequences and simulated data. AVAILABILITY: PhyloQFold, a C + + implementation of our method, is freely available from http://evol.bio.lmu.de/_statgen/software/phyloqfold/.

Mapping the RNA-Seq trash bin: unusual transcripts in prokaryotic transcriptome sequencing data.

Prokaryotic transcripts constitute almost always uninterrupted intervals when mapped back to the genome. Split reads, i.e., RNA-seq reads consisting of parts that only map to discontiguous loci, are thus disregarded in most analysis pipelines. There are, however, some well-known exceptions, in particular, tRNA splicing and circularized small RNAs in Archaea as well as self-splicing introns. Here, we reanalyze a series of published RNA-seq data sets, screening them specifically for non-contiguously mapping reads. We recover most of the known cases together with several novel archaeal ncRNAs associated with circularized products. In Eubacteria, only a handful of interesting candidates were obtained beyond a few previously described group I and group II introns. Most of the atypically mapping reads do not appear to correspond to well-defined, specifically processed products. Whether this diffuse background is, at least in part, an incidental by-product of prokaryotic RNA processing or whether it consists entirely of technical artifacts of reverse transcription or amplification remains unknown.

LINC00892 Is an lncRNA Induced by T Cell Activation and Expressed by Follicular Lymphoma-Resident T Helper Cells.

Successful immunotherapy in both solid tumors and in hematological malignancies relies on the ability of T lymphocytes to infiltrate the cancer tissue and mount an immune response against the tumor. Biomarkers able to discern the amount and the types of T lymphocytes infiltrating a given tumor therefore have high diagnostic and prognostic value. Given that lncRNAs are known to have a highly cell-type-specific expression pattern, we searched for lncRNAs specifically expressed by activated T cells and at the same time in a kind of lymphoma, follicular lymphoma, where the microenvironment is known to play a critical role in the regulation of antitumor immunity. We focused on a non-coding transcript, annotated as LINC00892, which reaches extremely high expression levels following cell activation in Jurkat cells. Interestingly LINC00892 has an expression pattern resembling that of genes involved in T cell memory. Accordingly, LINC00892 is mostly expressed by the effector memory and helper CD4+ T cell sub-types but not by naïve T cells. In situ analyses of LINC00892 expression in normal lymph nodes and in follicular lymphoma biopsies show that its expression is limited to CD4+ PD1hi T cells, with a subcellular localization within the germinal center matching that of follicular helper T cells. Our analysis therefore suggests that the previously uncharacterized lncRNA LINC00892 could be a useful biomarker for the detection of CD4+ memory T cells in both normal and tumor tissues.

Detection of Atherosclerosis by Small RNA-Sequencing Analysis of Extracellular Vesicle Enriched Serum Samples.

Atherosclerosis can occur throughout the arterial vascular system and lead to various diseases. Early diagnosis of atherosclerotic processes and of individual disease patterns would be more likely to be successful if targeted therapies were available. For this, it is important to find reliable biomarkers that are easily accessible and with little inconvenience for patients. There are many cell culture, animal model or tissue studies that found biomarkers at the microRNA (miRNA) and mRNA level describing atherosclerotic processes. However, little is known about their potential as circulating and liquid biopsy markers in patients. In this study, we examined serum-derived miRNA - profiles from 129 patients and 28 volunteers to identify potential biomarkers. The patients had four different atherosclerotic manifestations: abdominal aneurysm (n = 35), coronary heart disease (n = 34), carotid artery stenosis (n = 24) and peripheral arterial disease (n = 36). The samples were processed with an extracellular vesicle enrichment protocol, total-RNA extraction and small RNA-sequencing were performed. A differential expression analysis was performed bioinformatically to find potentially regulated miRNA biomarkers. Resulting miRNA candidates served as a starting point for an overrepresentation analysis in which relevant target mRNAs were identified. The Gene Ontology database revealed relevant biological functions in relation to atherosclerotic processes. In patients, expression of specific miRNAs changed significantly compared to healthy volunteers; 27 differentially expressed miRNAs were identified. We were able to detect a group-specific miRNA fingerprint: miR-122-5p, miR-2110 and miR-483-5p for abdominal aortic aneurysm, miR-370-3p and miR-409-3p for coronary heart disease, miR-335-3p, miR-381-3p, miR493-5p and miR654-3p for carotid artery stenosis, miR-199a-5p, miR-215-5p, miR-3168, miR-582-3p and miR-769-5p for peripheral arterial disease. The results of the study show that some of the identified miRNAs have already been associated with atherosclerosis in previous studies. Overrepresentation analysis on this data detected biological processes that are clearly relevant for atherosclerosis, its development and progression showing the potential of these miRNAs as biomarker candidates. In a next step, the relevance of these findings on the mRNA level is to be investigated and substantiated.

Mutational mechanisms shaping the coding and noncoding genome of germinal center derived B-cell lymphomas.

B cells have the unique property to somatically alter their immunoglobulin (IG) genes by V(D)J recombination, somatic hypermutation (SHM) and class-switch recombination (CSR). Aberrant targeting of these mechanisms is implicated in lymphomagenesis, but the mutational processes are poorly understood. By performing whole genome and transcriptome sequencing of 181 germinal center derived B-cell lymphomas (gcBCL) we identified distinct mutational signatures linked to SHM and CSR. We show that not only SHM, but presumably also CSR causes off-target mutations in non-IG genes. Kataegis clusters with high mutational density mainly affected early replicating regions and were enriched for SHM- and CSR-mediated off-target mutations. Moreover, they often co-occurred in loci physically interacting in the nucleus, suggesting that mutation hotspots promote increased mutation targeting of spatially co-localized loci (termed hypermutation by proxy). Only around 1% of somatic small variants were in protein coding sequences, but in about half of the driver genes, a contribution of B-cell specific mutational processes to their mutations was found. The B-cell-specific mutational processes contribute to both lymphoma initiation and intratumoral heterogeneity. Overall, we demonstrate that mutational processes involved in the development of gcBCL are more complex than previously appreciated, and that B cell-specific mutational processes contribute via diverse mechanisms to lymphomagenesis.

Divergent evolution in the genomes of closely related lacertids, Lacerta viridis and L. bilineata, and implications for speciation.

BACKGROUND: Lacerta viridis and Lacerta bilineata are sister species of European green lizards (eastern and western clades, respectively) that, until recently, were grouped together as the L. viridis complex. Genetic incompatibilities were observed between lacertid populations through crossing experiments, which led to the delineation of two separate species within the L. viridis complex. The population history of these sister species and processes driving divergence are unknown. We constructed the first high-quality de novo genome assemblies for both L. viridis and L. bilineata through Illumina and PacBio sequencing, with annotation support provided from transcriptome sequencing of several tissues. To estimate gene flow between the two species and identify factors involved in reproductive isolation, we studied their evolutionary history, identified genomic rearrangements, detected signatures of selection on non-coding RNA, and on protein-coding genes. FINDINGS: Here we show that gene flow was primarily unidirectional from L. bilineata to L. viridis after their split at least 1.15 million years ago. We detected positive selection of the non-coding repertoire; mutations in transcription factors; accumulation of divergence through inversions; selection on genes involved in neural development, reproduction, and behavior, as well as in ultraviolet-response, possibly driven by sexual selection, whose contribution to reproductive isolation between these lacertid species needs to be further evaluated. CONCLUSION: The combination of short and long sequence reads resulted in one of the most complete lizard genome assemblies. The characterization of a diverse array of genomic features provided valuable insights into the demographic history of divergence among European green lizards, as well as key species differences, some of which are candidates that could have played a role in speciation. In addition, our study generated valuable genomic resources that can be used to address conservation-related issues in lacertids.

Genomic and transcriptomic changes complement each other in the pathogenesis of sporadic Burkitt lymphoma.

Burkitt lymphoma (BL) is the most common B-cell lymphoma in children. Within the International Cancer Genome Consortium (ICGC), we performed whole genome and transcriptome sequencing of 39 sporadic BL. Here, we unravel interaction of structural, mutational, and transcriptional changes, which contribute to MYC oncogene dysregulation together with the pathognomonic IG-MYC translocation. Moreover, by mapping IGH translocation breakpoints, we provide evidence that the precursor of at least a subset of BL is a B-cell poised to express IGHA. We describe the landscape of mutations, structural variants, and mutational processes, and identified a series of driver genes in the pathogenesis of BL, which can be targeted by various mechanisms, including IG-non MYC translocations, germline and somatic mutations, fusion transcripts, and alternative splicing.

RNA-seq analysis identifies different transcriptomic types and developmental trajectories of primary melanomas.

Recent studies revealed trajectories of mutational events in early melanomagenesis, but the accompanying changes in gene expression are far less understood. Therefore, we performed a comprehensive RNA-seq analysis of laser-microdissected melanocytic nevi (n = 23) and primary melanoma samples (n = 57) and characterized the molecular mechanisms of early melanoma development. Using self-organizing maps, unsupervised clustering, and analysis of pseudotime (PT) dynamics to identify evolutionary trajectories, we describe here two transcriptomic types of melanocytic nevi (N1 and N2) and primary melanomas (M1 and M2). N1/M1 lesions are characterized by pigmentation-type and MITF gene signatures, and a high prevalence of NRAS mutations in M1 melanomas. N2/M2 lesions are characterized by inflammatory-type and AXL gene signatures with an equal distribution of wild-type and mutated BRAF and low prevalence of NRAS mutations in M2 melanomas. Interestingly, N1 nevi and M1 melanomas and N2 nevi and M2 melanomas, respectively, cluster together, but there is no clustering in a stage-dependent manner. Transcriptional signatures of M1 melanomas harbor signatures of BRAF/MEK inhibitor resistance and M2 melanomas harbor signatures of anti-PD-1 antibody treatment resistance. Pseudotime dynamics of nevus and melanoma samples are suggestive for a switch-like immune-escape mechanism in melanoma development with downregulation of immune genes paralleled by an increasing expression of a cell cycle signature in late-stage melanomas. Taken together, the transcriptome analysis identifies gene signatures and mechanisms underlying development of melanoma in early and late stages with relevance for diagnostics and therapy.

Rare Splice Variants in Long Non-Coding RNAs.

Long non-coding RNAs (lncRNAs) form a substantial component of the transcriptome and are involved in a wide variety of regulatory mechanisms. Compared to protein-coding genes, they are often expressed at low levels and are restricted to a narrow range of cell types or developmental stages. As a consequence, the diversity of their isoforms is still far from being recorded and catalogued in its entirety, and the debate is ongoing about what fraction of non-coding RNAs truly conveys biological function rather than being "junk". Here, using a collection of more than 100 transcriptomes from related B cell lymphoma, we show that lncRNA loci produce a very defined set of splice variants. While some of them are so rare that they become recognizable only in the superposition of dozens or hundreds of transcriptome datasets and not infrequently include introns or exons that have not been included in available genome annotation data, there is still a very limited number of processing products for any given locus. The combined depth of our sequencing data is large enough to effectively exhaust the isoform diversity: the overwhelming majority of splice junctions that are observed at all are represented by multiple junction-spanning reads. We conclude that the human transcriptome produces virtually no background of RNAs that are processed at effectively random positions, but is-under normal circumstances-confined to a well defined set of splice variants.

Changes of bivalent chromatin coincide with increased expression of developmental genes in cancer.

Bivalent (poised or paused) chromatin comprises activating and repressing histone modifications at the same location. This combination of epigenetic marks at promoter or enhancer regions keeps genes expressed at low levels but poised for rapid activation. Typically, DNA at bivalent promoters is only lowly methylated in normal cells, but frequently shows elevated methylation levels in cancer samples. Here, we developed a universal classifier built from chromatin data that can identify cancer samples solely from hypermethylation of bivalent chromatin. Tested on over 7,000 DNA methylation data sets from several cancer types, it reaches an AUC of 0.92. Although higher levels of DNA methylation are often associated with transcriptional silencing, counter-intuitive positive statistical dependencies between DNA methylation and expression levels have been recently reported for two cancer types. Here, we re-analyze combined expression and DNA methylation data sets, comprising over 5,000 samples, and demonstrate that the conjunction of hypermethylation of bivalent chromatin and up-regulation of the corresponding genes is a general phenomenon in cancer. This up-regulation affects many developmental genes and transcription factors, including dozens of homeobox genes and other genes implicated in cancer. Thus, we reason that the disturbance of bivalent chromatin may be intimately linked to tumorigenesis.

Differential transcriptional responses to Ebola and Marburg virus infection in bat and human cells.

The unprecedented outbreak of Ebola in West Africa resulted in over 28,000 cases and 11,000 deaths, underlining the need for a better understanding of the biology of this highly pathogenic virus to develop specific counter strategies. Two filoviruses, the Ebola and Marburg viruses, result in a severe and often fatal infection in humans. However, bats are natural hosts and survive filovirus infections without obvious symptoms. The molecular basis of this striking difference in the response to filovirus infections is not well understood. We report a systematic overview of differentially expressed genes, activity motifs and pathways in human and bat cells infected with the Ebola and Marburg viruses, and we demonstrate that the replication of filoviruses is more rapid in human cells than in bat cells. We also found that the most strongly regulated genes upon filovirus infection are chemokine ligands and transcription factors. We observed a strong induction of the JAK/STAT pathway, of several genes encoding inhibitors of MAP kinases (DUSP genes) and of PPP1R15A, which is involved in ER stress-induced cell death. We used comparative transcriptomics to provide a data resource that can be used to identify cellular responses that might allow bats to survive filovirus infections.

DNA methylome analysis in Burkitt and follicular lymphomas identifies differentially methylated regions linked to somatic mutation and transcriptional control.

Although Burkitt lymphomas and follicular lymphomas both have features of germinal center B cells, they are biologically and clinically quite distinct. Here we performed whole-genome bisulfite, genome and transcriptome sequencing in 13 IG-MYC translocation-positive Burkitt lymphoma, nine BCL2 translocation-positive follicular lymphoma and four normal germinal center B cell samples. Comparison of Burkitt and follicular lymphoma samples showed differential methylation of intragenic regions that strongly correlated with expression of associated genes, for example, genes active in germinal center dark-zone and light-zone B cells. Integrative pathway analyses of regions differentially methylated in Burkitt and follicular lymphomas implicated DNA methylation as cooperating with somatic mutation of sphingosine phosphate signaling, as well as the TCF3-ID3 and SWI/SNF complexes, in a large fraction of Burkitt lymphomas. Taken together, our results demonstrate a tight connection between somatic mutation, DNA methylation and transcriptional control in key B cell pathways deregulated differentially in Burkitt lymphoma and other germinal center B cell lymphomas.

A multi-split mapping algorithm for circular RNA, splicing, trans-splicing and fusion detection.

Numerous high-throughput sequencing studies have focused on detecting conventionally spliced mRNAs in RNA-seq data. However, non-standard RNAs arising through gene fusion, circularization or trans-splicing are often neglected. We introduce a novel, unbiased algorithm to detect splice junctions from single-end cDNA sequences. In contrast to other methods, our approach accommodates multi-junction structures. Our method compares favorably with competing tools for conventionally spliced mRNAs and, with a gain of up to 40% of recall, systematically outperforms them on reads with multiple splits, trans-splicing and circular products. The algorithm is integrated into our mapping tool segemehl (http://www.bioinf.uni-leipzig.de/Software/segemehl/).

PROSO II--a new method for protein solubility prediction.

Many fields of science and industry depend on efficient production of active protein using heterologous expression in Escherichia coli. The solubility of proteins upon expression is dependent on their amino acid sequence. Prediction of solubility from sequence is therefore highly valuable. We present a novel machine-learning-based model called PROSO II which makes use of new classification methods and growth in experimental data to improve coverage and accuracy of solubility predictions. The classification algorithm is organized as a two-layered structure in which the output of a primary Parzen window model for sequence similarity and a logistic regression classifier of amino acid k-mer composition serve as input for a second-level logistic regression classifier. Compared with previously published research our model is trained on five times more data than used by any other method before (82 000 proteins). When tested on a separate holdout set not used at any point of method development our server attained the best results in comparison with other currently available methods: accuracy 75.4%, Matthew's correlation coefficient 0.39, sensitivity 0.731, specificity 0.759, gain (soluble) 2.263. In summary, due to utilization of cutting edge machine learning technologies combined with the largest currently available experimental data set the PROSO II server constitutes a substantial improvement in protein solubility predictions. PROSO II is available at http://mips.helmholtz-muenchen.de/prosoII.