Large-scale benchmarking of circRNA detection tools reveals large differences in sensitivity but not in precision.

The detection of circular RNA molecules (circRNAs) is typically based on short-read RNA sequencing data processed using computational tools. Numerous such tools have been developed, but a systematic comparison with orthogonal validation is missing. Here, we set up a circRNA detection tool benchmarking study, in which 16 tools detected more than 315,000 unique circRNAs in three deeply sequenced human cell types. Next, 1,516 predicted circRNAs were validated using three orthogonal methods. Generally, tool-specific precision is high and similar (median of 98.8%, 96.3% and 95.5% for qPCR, RNase R and amplicon sequencing, respectively) whereas the sensitivity and number of predicted circRNAs (ranging from 1,372 to 58,032) are the most significant differentiators. Of note, precision values are lower when evaluating low-abundance circRNAs. We also show that the tools can be used complementarily to increase detection sensitivity. Finally, we offer recommendations for future circRNA detection and validation.

Synthesizing genome regulation data with vote-counting.

The increasing availability of high-throughput datasets allows amalgamating research information across a large body of genome regulation studies. Given the recent success of meta-analyses on transcriptional regulators, epigenetic marks, and enhancer:gene associations, we expect that such surveys will continue to provide novel and reproducible insights. However, meta-analyses are severely hampered by the diversity of available data, concurring protocols, an eclectic amount of bioinformatics tools, and myriads of conceivable parameter combinations. Such factors can easily bar life scientists from synthesizing omics data and substantially curb their interpretability. Despite statistical challenges of the method, we would like to emphasize the advantages of joining data from different sources through vote-counting and showcase examples that achieve a simple but highly intuitive data integration.

Locus-specific expression analysis of transposable elements.

Transposable elements (TEs) have been associated with many, frequently detrimental, biological roles. Consequently, the regulations of TEs, e.g. via DNA-methylation and histone modifications, are considered critical for maintaining genomic integrity and other functions. Still, the high-throughput study of TEs is usually limited to the family or consensus-sequence level because of alignment problems prompted by high-sequence similarities and short read lengths. To entirely comprehend the effects and reasons of TE expression, however, it is necessary to assess the TE expression at the level of individual instances. Our simulation study demonstrates that sequence similarities and short read lengths do not rule out the accurate assessment of (differential) expression of TEs at the instance-level. With only slight modifications to existing methods, TE expression analysis works surprisingly well for conventional paired-end sequencing data. We find that SalmonTE and Telescope can accurately tally a considerable amount of TE instances, allowing for differential expression recovery in model and non-model organisms.

Implementing FAIR data management within the German Network for Bioinformatics Infrastructure (de.NBI) exemplified by selected use cases.

This article describes some use case studies and self-assessments of FAIR status of de.NBI services to illustrate the challenges and requirements for the definition of the needs of adhering to the FAIR (findable, accessible, interoperable and reusable) data principles in a large distributed bioinformatics infrastructure. We address the challenge of heterogeneity of wet lab technologies, data, metadata, software, computational workflows and the levels of implementation and monitoring of FAIR principles within the different bioinformatics sub-disciplines joint in de.NBI. On the one hand, this broad service landscape and the excellent network of experts are a strong basis for the development of useful research data management plans. On the other hand, the large number of tools and techniques maintained by distributed teams renders FAIR compliance challenging.

Focal structural variants revealed by whole genome sequencing disrupt the histone demethylase KDM4C in B-cell lymphomas.

Histone methylation-modifiers, such as EZH2 and KMT2D, are recurrently altered in B-cell lymphomas. To comprehensively describe the landscape of alterations affecting genes encoding histone methylation-modifiers in lymphomagenesis we investigated whole genome and transcriptome data of 186 mature B-cell lymphomas sequenced in the ICGC MMML-Seq project. Besides confirming common alterations of KMT2D (47% of cases), EZH2 (17%), SETD1B (5%), PRDM9 (4%), KMT2C (4%), and SETD2 (4%), also identified by prior exome or RNA-sequencing studies, we here found recurrent alterations to KDM4C in chromosome 9p24, encoding a histone demethylase. Focal structural variation was the main mechanism of KDM4C alterations, and was independent from 9p24 amplification. We also identified KDM4C alterations in lymphoma cell lines including a focal homozygous deletion in a classical Hodgkin lymphoma cell line. By integrating RNA-sequencing and genome sequencing data we predict that KDM4C structural variants result in loss-offunction. By functional reconstitution studies in cell lines, we provide evidence that KDM4C can act as a tumor suppressor. Thus, we show that identification of structural variants in whole genome sequencing data adds to the comprehensive description of the mutational landscape of lymphomas and, moreover, establish KDM4C as a putative tumor suppressive gene recurrently altered in subsets of B-cell derived lymphomas.

Increased prostaglandin-D2 in male STAT3-deficient hearts shifts cardiac progenitor cells from endothelial to white adipocyte differentiation.

Cardiac levels of the signal transducer and activator of transcription factor-3 (STAT3) decline with age, and male but not female mice with a cardiomyocyte-specific STAT3 deficiency conditional knockout (CKO) display premature age-related heart failure associated with reduced cardiac capillary density. In the present study, isolated male and female CKO-cardiomyocytes exhibit increased prostaglandin (PG)-generating cyclooxygenase-2 (COX-2) expression. The PG-degrading hydroxyprostaglandin-dehydrogenase-15 (HPGD) expression is only reduced in male cardiomyocytes, which is associated with increased prostaglandin D2 (PGD2) secretion from isolated male but not female CKO-cardiomyocytes. Reduced HPGD expression in male cardiomyocytes derive from impaired androgen receptor (AR)-signaling due to loss of its cofactor STAT3. Elevated PGD2 secretion in males is associated with increased white adipocyte accumulation in aged male but not female hearts. Adipocyte differentiation is enhanced in isolated stem cell antigen-1 (SCA-1)+ cardiac progenitor cells (CPC) from young male CKO-mice compared with the adipocyte differentiation of male wild-type (WT)-CPC and CPC isolated from female mice. Epigenetic analysis in freshly isolated male CKO-CPC display hypermethylation in pro-angiogenic genes (Fgfr2, Epas1) and hypomethylation in the white adipocyte differentiation gene Zfp423 associated with up-regulated ZFP423 expression and a shift from endothelial to white adipocyte differentiation compared with WT-CPC. The expression of the histone-methyltransferase EZH2 is reduced in male CKO-CPC compared with male WT-CPC, whereas no differences in the EZH2 expression in female CPC were observed. Clonally expanded CPC can differentiate into endothelial cells or into adipocytes depending on the differentiation conditions. ZFP423 overexpression is sufficient to induce white adipocyte differentiation of clonal CPC. In isolated WT-CPC, PGD2 stimulation reduces the expression of EZH2, thereby up-regulating ZFP423 expression and promoting white adipocyte differentiation. The treatment of young male CKO mice with the COX inhibitor Ibuprofen or the PGD2 receptor (DP)2 receptor antagonist BAY-u 3405 in vivo increased EZH2 expression and reduced ZFP423 expression and adipocyte differentiation in CKO-CPC. Thus, cardiomyocyte STAT3 deficiency leads to age-related and sex-specific cardiac remodeling and failure in part due to sex-specific alterations in PGD2 secretion and subsequent epigenetic impairment of the differentiation potential of CPC. Causally involved is the impaired AR signaling in absence of STAT3, which reduces the expression of the PG-degrading enzyme HPGD.

Transcription factor RFX7 governs a tumor suppressor network in response to p53 and stress.

Despite its prominence, the mechanisms through which the tumor suppressor p53 regulates most genes remain unclear. Recently, the regulatory factor X 7 (RFX7) emerged as a suppressor of lymphoid neoplasms, but its regulation and target genes mediating tumor suppression remain unknown. Here, we identify a novel p53-RFX7 signaling axis. Integrative analysis of the RFX7 DNA binding landscape and the RFX7-regulated transcriptome in three distinct cell systems reveals that RFX7 directly controls multiple established tumor suppressors, including PDCD4, PIK3IP1, MXD4, and PNRC1, across cell types and is the missing link for their activation in response to p53 and stress. RFX7 target gene expression correlates with cell differentiation and better prognosis in numerous cancer types. Interestingly, we find that RFX7 sensitizes cells to Doxorubicin by promoting apoptosis. Together, our work establishes RFX7's role as a ubiquitous regulator of cell growth and fate determination and a key node in the p53 transcriptional program.

An Analysis of Methylome Evolution in Primates.

Although the investigation of the epigenome becomes increasingly important, still little is known about the long-term evolution of epigenetic marks and systematic investigation strategies are still lacking. Here, we systematically demonstrate the transfer of classic phylogenetic methods such as maximum likelihood based on substitution models, parsimony, and distance-based to interval-scaled epigenetic data. Using a great apes blood data set, we demonstrate that DNA methylation is evolutionarily conserved at the level of individual CpGs in promotors, enhancers, and genic regions. Our analysis also reveals that this epigenomic conservation is significantly correlated with its transcription factor binding density. Binding sites for transcription factors involved in neuron differentiation and components of AP-1 evolve at a significantly higher rate at methylation than at the nucleotide level. Moreover, our models suggest an accelerated epigenomic evolution at binding sites of BRCA1, chromobox homolog protein 2, and factors of the polycomb repressor 2 complex in humans. For most genomic regions, the methylation-based reconstruction of phylogenetic trees is at par with sequence-based reconstruction. Most strikingly, phylogenetic reconstruction using methylation rates in enhancer regions was ineffective independently of the chosen model. We identify a set of phylogenetically uninformative CpG sites enriched in enhancers controlling immune-related genes.

DREAM and RB cooperate to induce gene repression and cell-cycle arrest in response to p53 activation.

Most human cancers acquire mutations causing defects in the p53 signaling pathway. The tumor suppressor p53 becomes activated in response to genotoxic stress and is essential for arresting the cell cycle to facilitate DNA repair or to initiate apoptosis. p53-induced cell cycle-arrest is mediated by expression of the CDK inhibitor p21WAF1/Cip1, which prevents phosphorylation and inactivation of the pocket proteins RB, p130, and p107. In a hypophosphorylated state, pocket proteins bind to E2F factors forming RB-E2F and DREAM transcriptional repressor complexes. Here, we analyze the influence of RB and DREAM on p53-induced gene repression and cell-cycle arrest. We show that abrogation of DREAM function by knockout of the DREAM component LIN37 results in a reduced repression of cell-cycle genes. We identify the genes repressed by the p53-DREAM pathway and describe a set of genes that is downregulated by p53 independent of LIN37/DREAM. Most strikingly, p53-dependent repression of cell-cycle genes is completely abrogated in LIN37-/-;RB-/- cells leading to a loss of the G1/S checkpoint. Taken together, we show that DREAM and RB are key factors in the p53 signaling pathway to downregulate a large number of cell-cycle genes and to arrest the cell cycle at the G1/S transition.

Evidence-Based Study Design for Hospital-Acquired Bacterial Pneumonia and Ventilator-Associated Bacterial Pneumonia.

BACKGROUND: The US Food and Drug Administration solicited evidence-based recommendations to improve guidance for studies of hospital-acquired bacterial pneumonia (HABP) and ventilator-associated bacterial pneumonia (VABP). METHODS: We analyzed 7 HABP/VABP datasets to explore novel noninferiority study endpoints and designs, focusing on alternatives to all-cause mortality (ACM). RESULTS: ACM at day 28 differed for ventilated HABP (27.8%), VABP (18.0%), and nonventilated HABP (14.5%). A "mortality-plus" (ACM+) composite endpoint was constructed by combining ACM with patient-relevant, infection-related adverse events from the Medical Dictionary for Regulatory Activities toxic/septic shock standardized query. The ACM+ rate was 3-10 percentage points above that of ACM across the studies and treatment groups. Predictors of higher ACM/ACM+ rates included older age and elevated acute physiology and chronic health evaluation (APACHE) II score. Only patients in the nonventilated HABP group were able to report pneumonia symptom changes. CONCLUSIONS: If disease groups and patient characteristics in future studies produce an ACM rate so low (<10%-15%) that a fixed noninferiority margin of 10% cannot be justified (requiring an odds ratio analysis), an ACM+ endpoint could lower sample size. Enrichment of studies with patients with a higher severity of illness would increase ACM. Data on symptom resolution in nonventilated HABP support development of a patient-reported outcome instrument.

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.

Comparison between transcriptomic responses to short-term stress exposures of a common Holarctic and endemic Lake Baikal amphipods.

BACKGROUND: Lake Baikal is one of the oldest freshwater lakes and has constituted a stable environment for millions of years, in stark contrast to small, transient bodies of water in its immediate vicinity. A highly diverse endemic endemic amphipod fauna is found in one, but not the other habitat. We ask here whether differences in stress response can explain the immiscibility barrier between Lake Baikal and non-Baikal faunas. To this end, we conducted exposure experiments to increased temperature and the toxic heavy metal cadmium as stressors. RESULTS: Here we obtained high-quality de novo transcriptome assemblies, covering mutiple conditions, of three amphipod species, and compared their transcriptomic stress responses. Two of these species, Eulimnogammarus verrucosus and E. cyaneus, are endemic to Lake Baikal, while the Holarctic Gammarus lacustris is a potential invader. CONCLUSIONS: Both Baikal species possess intact stress response systems and respond to elevated temperature with relatively similar changes in their expression profiles. G. lacustris reacts less strongly to the same stressors, possibly because its transcriptome is already perturbed by acclimation conditions.

metilene: fast and sensitive calling of differentially methylated regions from bisulfite sequencing data.

The detection of differentially methylated regions (DMRs) is a necessary prerequisite for characterizing different epigenetic states. We present a novel program, metilene, to identify DMRs within whole-genome and targeted data with unrivaled specificity and sensitivity. A binary segmentation algorithm combined with a two-dimensional statistical test allows the detection of DMRs in large methylation experiments with multiple groups of samples in minutes rather than days using off-the-shelf hardware. metilene outperforms other state-of-the-art tools for low coverage data and can estimate missing data. Hence, metilene is a versatile tool to study the effect of epigenetic modifications in differentiation/development, tumorigenesis, and systems biology on a global, genome-wide level. Whether in the framework of international consortia with dozens of samples per group, or even without biological replicates, it produces highly significant and reliable results.

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.

The African coelacanth genome provides insights into tetrapod evolution.

The discovery of a living coelacanth specimen in 1938 was remarkable, as this lineage of lobe-finned fish was thought to have become extinct 70 million years ago. The modern coelacanth looks remarkably similar to many of its ancient relatives, and its evolutionary proximity to our own fish ancestors provides a glimpse of the fish that first walked on land. Here we report the genome sequence of the African coelacanth, Latimeria chalumnae. Through a phylogenomic analysis, we conclude that the lungfish, and not the coelacanth, is the closest living relative of tetrapods. Coelacanth protein-coding genes are significantly more slowly evolving than those of tetrapods, unlike other genomic features. Analyses of changes in genes and regulatory elements during the vertebrate adaptation to land highlight genes involved in immunity, nitrogen excretion and the development of fins, tail, ear, eye, brain and olfaction. Functional assays of enhancers involved in the fin-to-limb transition and in the emergence of extra-embryonic tissues show the importance of the coelacanth genome as a blueprint for understanding tetrapod evolution.

The RNA workbench: best practices for RNA and high-throughput sequencing bioinformatics in Galaxy.

RNA-based regulation has become a major research topic in molecular biology. The analysis of epigenetic and expression data is therefore incomplete if RNA-based regulation is not taken into account. Thus, it is increasingly important but not yet standard to combine RNA-centric data and analysis tools with other types of experimental data such as RNA-seq or ChIP-seq. Here, we present the RNA workbench, a comprehensive set of analysis tools and consolidated workflows that enable the researcher to combine these two worlds. Based on the Galaxy framework the workbench guarantees simple access, easy extension, flexible adaption to personal and security needs, and sophisticated analyses that are independent of command-line knowledge. Currently, it includes more than 50 bioinformatics tools that are dedicated to different research areas of RNA biology including RNA structure analysis, RNA alignment, RNA annotation, RNA-protein interaction, ribosome profiling, RNA-seq analysis and RNA target prediction. The workbench is developed and maintained by experts in RNA bioinformatics and the Galaxy framework. Together with the growing community evolving around this workbench, we are committed to keep the workbench up-to-date for future standards and needs, providing researchers with a reliable and robust framework for RNA data analysis. AVAILABILITY: The RNA workbench is available at https://github.com/bgruening/galaxy-rna-workbench.

Genomic and transcriptomic heterogeneity of colorectal tumours arising in Lynch syndrome.

Colorectal cancer (CRC) arising in Lynch syndrome (LS) comprises tumours with constitutional mutations in DNA mismatch repair genes. There is still a lack of whole-genome and transcriptome studies of LS-CRC to address questions about similarities and differences in mutation and gene expression characteristics between LS-CRC and sporadic CRC, about the molecular heterogeneity of LS-CRC, and about specific mechanisms of LS-CRC genesis linked to dysfunctional mismatch repair in LS colonic mucosa and the possible role of immune editing. Here, we provide a first molecular characterization of LS tumours and of matched tumour-distant reference colonic mucosa based on whole-genome DNA-sequencing and RNA-sequencing analyses. Our data support two subgroups of LS-CRCs, G1 and G2, whereby G1 tumours show a higher number of somatic mutations, a higher amount of microsatellite slippage, and a different mutation spectrum. The gene expression phenotypes support this difference. Reference mucosa of G1 shows a strong immune response associated with the expression of HLA and immune checkpoint genes and the invasion of CD4+ T cells. Such an immune response is not observed in LS tumours, G2 reference and normal (non-Lynch) mucosa, and sporadic CRC. We hypothesize that G1 tumours are edited for escape from a highly immunogenic microenvironment via loss of HLA presentation and T-cell exhaustion. In contrast, G2 tumours seem to develop in a less immunogenic microenvironment where tumour-promoting inflammation parallels tumourigenesis. Larger studies on non-neoplastic mucosa tissue of mutation carriers are required to better understand the early phases of emerging tumours. Copyright © 2017 Pathological Society of Great Britain and Ireland. Published by John Wiley & Sons, Ltd.

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.

Patient-Reported Outcome Assessments as Endpoints in Studies in Infectious Diseases.

The goal of administering medical interventions is to help patients live longer or live better. In keeping with this goal, there has been increasing interest in taking the "voice" of the patient into account during the development process, specifically in the evaluation of treatment benefits of medical interventions, and use of patient-centered outcome data to justify reimbursement. Patient-reported outcomes (PROs) are outcome assessments (OAs) used to define endpoints that can provide direct evidence of treatment benefit on how patients feel or function. When PROs are appropriately developed, they can increase the efficiency and clinical relevance of clinical trials. Several PROs have been developed for OA in specific infectious diseases indications, and more are under development. PROs also hold promise for use in evaluating adherence, adverse effects, satisfaction with care, and routine clinical practice.

Alterations of microRNA and microRNA-regulated messenger RNA expression in germinal center B-cell lymphomas determined by integrative sequencing analysis.

MicroRNA are well-established players in post-transcriptional gene regulation. However, information on the effects of microRNA deregulation mainly relies on bioinformatic prediction of potential targets, whereas proof of the direct physical microRNA/target messenger RNA interaction is mostly lacking. Within the International Cancer Genome Consortium Project "Determining Molecular Mechanisms in Malignant Lymphoma by Sequencing", we performed miRnome sequencing from 16 Burkitt lymphomas, 19 diffuse large B-cell lymphomas, and 21 follicular lymphomas. Twenty-two miRNA separated Burkitt lymphomas from diffuse large B-cell lymphomas/follicular lymphomas, of which 13 have shown regulation by MYC. Moreover, we found expression of three hitherto unreported microRNA. Additionally, we detected recurrent mutations of hsa-miR-142 in diffuse large B-cell lymphomas and follicular lymphomas, and editing of the hsa-miR-376 cluster, providing evidence for microRNA editing in lymphomagenesis. To interrogate the direct physical interactions of microRNA with messenger RNA, we performed Argonaute-2 photoactivatable ribonucleoside-enhanced cross-linking and immunoprecipitation experiments. MicroRNA directly targeted 208 messsenger RNA in the Burkitt lymphomas and 328 messenger RNA in the non-Burkitt lymphoma models. This integrative analysis discovered several regulatory pathways of relevance in lymphomagenesis including Ras, PI3K-Akt and MAPK signaling pathways, also recurrently deregulated in lymphomas by mutations. Our dataset reveals that messenger RNA deregulation through microRNA is a highly relevant mechanism in lymphomagenesis.