Meta-Analysis of Mechano-Sensitive Ion Channels in Human Hearts: Chamber- and Disease-Preferential mRNA Expression.

The cardiac cell mechanical environment changes on a beat-by-beat basis as well as in the course of various cardiac diseases. Cells sense and respond to mechanical cues via specialized mechano-sensors initiating adaptive signaling cascades. With the aim of revealing new candidates underlying mechano-transduction relevant to cardiac diseases, we investigated mechano-sensitive ion channels (MSC) in human hearts for their chamber- and disease-preferential mRNA expression. Based on a meta-analysis of RNA sequencing studies, we compared the mRNA expression levels of MSC in human atrial and ventricular tissue samples from transplant donor hearts (no cardiac disease), and from patients in sinus rhythm (underlying diseases: heart failure, coronary artery disease, heart valve disease) or with atrial fibrillation. Our results suggest that a number of MSC genes are expressed chamber preferentially, e.g., CHRNE in the atria (compared to the ventricles), TRPV4 in the right atrium (compared to the left atrium), CACNA1B and KCNMB1 in the left atrium (compared to the right atrium), as well as KCNK2 and KCNJ2 in ventricles (compared to the atria). Furthermore, 15 MSC genes are differentially expressed in cardiac disease, out of which SCN9A (lower expressed in heart failure compared to donor tissue) and KCNQ5 (lower expressed in atrial fibrillation compared to sinus rhythm) show a more than twofold difference, indicative of possible functional relevance. Thus, we provide an overview of cardiac MSC mRNA expression in the four cardiac chambers from patients with different cardiac diseases. We suggest that the observed differences in MSC mRNA expression may identify candidates involved in altered mechano-transduction in the respective diseases.

A SARS-CoV-2 sequence submission tool for the European Nucleotide Archive.

SUMMARY: Many aspects of the global response to the COVID-19 pandemic are enabled by the fast and open publication of SARS-CoV-2 genetic sequence data. The European Nucleotide Archive (ENA) is the European recommended open repository for genetic sequences. In this work, we present a tool for submitting raw sequencing reads of SARS-CoV-2 to ENA. The tool features a single-step submission process, a graphical user interface, tabular-formatted metadata and the possibility to remove human reads prior to submission. A Galaxy wrap of the tool allows users with little or no bioinformatics knowledge to do bulk sequencing read submissions. The tool is also packed in a Docker container to ease deployment. AVAILABILITY AND IMPLEMENTATION: CLI ENA upload tool is available at github.com/usegalaxy-eu/ena-upload-cli (DOI 10.5281/zenodo.4537621); Galaxy ENA upload tool at toolshed.g2.bx.psu.edu/view/iuc/ena_upload/382518f24d6d and github.com/galaxyproject/tools-iuc/tree/master/tools/ena_upload (development); and ENA upload Galaxy container at github.com/ELIXIR-Belgium/ena-upload-container (DOI 10.5281/zenodo.4730785).

Deciphering human ribonucleoprotein regulatory networks.

RNA-binding proteins (RBPs) control and coordinate each stage in the life cycle of RNAs. Although in vivo binding sites of RBPs can now be determined genome-wide, most studies typically focused on individual RBPs. Here, we examined a large compendium of 114 high-quality transcriptome-wide in vivo RBP-RNA cross-linking interaction datasets generated by the same protocol in the same cell line and representing 64 distinct RBPs. Comparative analysis of categories of target RNA binding preference, sequence preference, and transcript region specificity was performed, and identified potential posttranscriptional regulatory modules, i.e. specific combinations of RBPs that bind to specific sets of RNAs and targeted regions. These regulatory modules represented functionally related proteins and exhibited distinct differences in RNA metabolism, expression variance, as well as subcellular localization. This integrative investigation of experimental RBP-RNA interaction evidence and RBP regulatory function in a human cell line will be a valuable resource for understanding the complexity of post-transcriptional regulation.

Temporospatial distribution and transcriptional profile of retinal microglia in the oxygen-induced retinopathy mouse model.

Myeloid cells such as resident retinal microglia (MG) or infiltrating blood-derived macrophages (Mϕ) accumulate in areas of retinal ischemia and neovascularization (RNV) and modulate neovascular eye disease. Their temporospatial distribution and biological function in this process, however, remain unclarified. Using state-of-the-art methods, including cell-specific reporter mice and high-throughput RNA sequencing (RNA Seq), this study determined the extent of MG proliferation and Mϕ infiltration in areas with retinal ischemia and RNV in Cx3cr1CreERT2 :Rosa26-tdTomato mice and examined the transcriptional profile of MG in the mouse model of oxygen-induced retinopathy (OIR). For RNA Seq, tdTomato-positive retinal MG were sorted by flow cytometry followed by Gene ontology (GO) cluster analysis. Furthermore, intraperitoneal injections of the cell proliferation marker 5-ethynyl-2'-deoxyuridine (EdU) were performed from postnatal day (p) 12 to p16. We found that MG is the predominant myeloid cell population while Mϕ rarely appears in areas of RNV. Thirty percent of retinal MG in areas of RNV were EdU-positive indicating a considerable local MG cell expansion. GO cluster analysis revealed an enrichment of clusters related to cell division, tubulin binding, ATPase activity, protein kinase regulatory activity, and chemokine receptor binding in MG in the OIR model compared to untreated controls. In conclusion, activated retinal MG alter their transcriptional profile, exhibit considerable proliferative ability and are by far the most frequent myeloid cell population in areas of ischemia and RNV in the OIR model thus presenting a potential target for future therapeutic approaches.

RCAS: an RNA centric annotation system for transcriptome-wide regions of interest.

In the field of RNA, the technologies for studying the transcriptome have created a tremendous potential for deciphering the puzzles of the RNA biology. Along with the excitement, the unprecedented volume of RNA related omics data is creating great challenges in bioinformatics analyses. Here, we present the RNA Centric Annotation System (RCAS), an R package, which is designed to ease the process of creating gene-centric annotations and analysis for the genomic regions of interest obtained from various RNA-based omics technologies. The design of RCAS is modular, which enables flexible usage and convenient integration with other bioinformatics workflows. RCAS is an R/Bioconductor package but we also created graphical user interfaces including a Galaxy wrapper and a stand-alone web service. The application of RCAS on published datasets shows that RCAS is not only able to reproduce published findings but also helps generate novel knowledge and hypotheses. The meta-gene profiles, gene-centric annotation, motif analysis and gene-set analysis provided by RCAS provide contextual knowledge which is necessary for understanding the functional aspects of different biological events that involve RNAs. In addition, the array of different interfaces and deployment options adds the convenience of use for different levels of users. RCAS is available at http://bioconductor.org/packages/release/bioc/html/RCAS.html and http://rcas.mdc-berlin.de.

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.

Small RNA profiling of low biomass samples: identification and removal of contaminants.

BACKGROUND: Sequencing-based analyses of low-biomass samples are known to be prone to misinterpretation due to the potential presence of contaminating molecules derived from laboratory reagents and environments. DNA contamination has been previously reported, yet contamination with RNA is usually considered to be very unlikely due to its inherent instability. Small RNAs (sRNAs) identified in tissues and bodily fluids, such as blood plasma, have implications for physiology and pathology, and therefore the potential to act as disease biomarkers. Thus, the possibility for RNA contaminants demands careful evaluation. RESULTS: Herein, we report on the presence of small RNA (sRNA) contaminants in widely used microRNA extraction kits and propose an approach for their depletion. We sequenced sRNAs extracted from human plasma samples and detected important levels of non-human (exogenous) sequences whose source could be traced to the microRNA extraction columns through a careful qPCR-based analysis of several laboratory reagents. Furthermore, we also detected the presence of artefactual sequences related to these contaminants in a range of published datasets, thereby arguing in particular for a re-evaluation of reports suggesting the presence of exogenous RNAs of microbial and dietary origin in blood plasma. To avoid artefacts in future experiments, we also devise several protocols for the removal of contaminant RNAs, define minimal amounts of starting material for artefact-free analyses, and confirm the reduction of contaminant levels for identification of bona fide sequences using 'ultra-clean' extraction kits. CONCLUSION: This is the first report on the presence of RNA molecules as contaminants in RNA extraction kits. The described protocols should be applied in the future to avoid confounding sRNA studies.

Animal snoRNAs and scaRNAs with exceptional structures.

The overwhelming majority of small nucleolar RNAs (snoRNAs) fall into two clearly defined classes characterized by distinctive secondary structures and sequence motifs. A small group of diverse ncRNAs, however, shares the hallmarks of one or both classes of snoRNAs but differs substantially from the norm in some respects. Here, we compile the available information on these exceptional cases, conduct a thorough homology search throughout the available metazoan genomes, provide improved and expanded alignments, and investigate the evolutionary histories of these ncRNA families as well as their mutual relationships.

Bcheck: a wrapper tool for detecting RNase P RNA genes.

BACKGROUND: Effective bioinformatics solutions are needed to tackle challenges posed by industrial-scale genome annotation. We present Bcheck, a wrapper tool which predicts RNase P RNA genes by combining the speed of pattern matching and sensitivity of covariance models. The core of Bcheck is a library of subfamily specific descriptor models and covariance models. RESULTS: Scanning all microbial genomes in GenBank identifies RNase P RNA genes in 98% of 1024 microbial chromosomal sequences within just 4 hours on single CPU. Comparing to existing annotations found in 387 of the GenBank files, Bcheck predictions have more intact structure and are automatically classified by subfamily membership. For eukaryotic chromosomes Bcheck could identify the known RNase P RNA genes in 84 out of 85 metazoan genomes and 19 out of 21 fungi genomes. Bcheck predicted 37 novel eukaryotic RNase P RNA genes, 32 of which are from fungi. Gene duplication events are observed in at least 20 metazoan organisms. Scanning of meta-genomic data from the Global Ocean Sampling Expedition, comprising over 10 million sample sequences (18 Gigabases), predicted 2909 unique genes, 98% of which fall into ancestral bacteria A type of RNase P RNA and 66% of which have no close homolog to known prokaryotic RNase P RNA. CONCLUSIONS: The combination of efficient filtering by means of a descriptor-based search and subsequent construction of a high-quality gene model by means of a covariance model provides an efficient method for the detection of RNase P RNA genes in large-scale sequencing data. Bcheck is implemented as webserver and can also be downloaded for local use from http://rna.tbi.univie.ac.at/bcheck.

Generation of pure monocultures of human microglia-like cells from induced pluripotent stem cells.

Microglia are resident tissue macrophages of the central nervous system (CNS) that arise from erythromyeloid progenitors during embryonic development. They play essential roles in CNS development, homeostasis and response to disease. Since microglia are difficult to procure from the human brain, several protocols have been developed to generate microglia-like cells from human induced pluripotent stem cells (hiPSCs). However, some concerns remain over the purity and quality of in vitro generated microglia. Here, we describe a new protocol that does not require co-culture with neural cells and yields cultures of 100% P2Y12+ 95% TMEM119+ ramified human microglia-like cells (hiPSC-MG). In the presence of neural precursor cell-conditioned media, hiPSC-MG expressed high levels of human microglia signature genes, including SALL1, CSF1R, P2RY12, TMEM119, TREM2, HEXB and SIGLEC11, as revealed by whole-transcriptome analysis. Stimulation of hiPSC-MG with lipopolysaccharide resulted in downregulation of P2Y12 expression, induction of IL1B mRNA expression and increase in cell capacitance. HiPSC-MG were phagocytically active and maintained their cell identity after transplantation into murine brain slices and human brain spheroids. Together, our new protocol for the generation of microglia-like cells from human iPSCs will facilitate the study of human microglial function in health and disease.

Reproducible inference of transcription factor footprints in ATAC-seq and DNase-seq datasets using protocol-specific bias modeling.

BACKGROUND: DNase-seq and ATAC-seq are broadly used methods to assay open chromatin regions genome-wide. The single nucleotide resolution of DNase-seq has been further exploited to infer transcription factor binding sites (TFBSs) in regulatory regions through footprinting. Recent studies have demonstrated the sequence bias of DNase I and its adverse effects on footprinting efficiency. However, footprinting and the impact of sequence bias have not been extensively studied for ATAC-seq. RESULTS: Here, we undertake a systematic comparison of the two methods and show that a modification to the ATAC-seq protocol increases its yield and its agreement with DNase-seq data from the same cell line. We demonstrate that the two methods have distinct sequence biases and correct for these protocol-specific biases when performing footprinting. Despite the differences in footprint shapes, the locations of the inferred footprints in ATAC-seq and DNase-seq are largely concordant. However, the protocol-specific sequence biases in conjunction with the sequence content of TFBSs impact the discrimination of footprint from the background, which leads to one method outperforming the other for some TFs. Finally, we address the depth required for reproducible identification of open chromatin regions and TF footprints. CONCLUSIONS: We demonstrate that the impact of bias correction on footprinting performance is greater for DNase-seq than for ATAC-seq and that DNase-seq footprinting leads to better performance. It is possible to infer concordant footprints by using replicates, highlighting the importance of reproducibility assessment. The results presented here provide an overview of the advantages and limitations of footprinting analyses using ATAC-seq and DNase-seq.

Community-Driven Data Analysis Training for Biology.

The primary problem with the explosion of biomedical datasets is not the data, not computational resources, and not the required storage space, but the general lack of trained and skilled researchers to manipulate and analyze these data. Eliminating this problem requires development of comprehensive educational resources. Here we present a community-driven framework that enables modern, interactive teaching of data analytics in life sciences and facilitates the development of training materials. The key feature of our system is that it is not a static but a continuously improved collection of tutorials. By coupling tutorials with a web-based analysis framework, biomedical researchers can learn by performing computation themselves through a web browser without the need to install software or search for example datasets. Our ultimate goal is to expand the breadth of training materials to include fundamental statistical and data science topics and to precipitate a complete re-engineering of undergraduate and graduate curricula in life sciences. This project is accessible at https://training.galaxyproject.org.

Strategies for analyzing bisulfite sequencing data.

DNA methylation is one of the main epigenetic modifications in the eukaryotic genome; it has been shown to play a role in cell-type specific regulation of gene expression, and therefore cell-type identity. Bisulfite sequencing is the gold-standard for measuring methylation over the genomes of interest. Here, we review several techniques used for the analysis of high-throughput bisulfite sequencing. We introduce specialized short-read alignment techniques as well as pre/post-alignment quality check methods to ensure data quality. Furthermore, we discuss subsequent analysis steps after alignment. We introduce various differential methylation methods and compare their performance using simulated and real bisulfite sequencing datasets. We also discuss the methods used to segment methylomes in order to pinpoint regulatory regions. We introduce annotation methods that can be used for further classification of regions returned by segmentation and differential methylation methods. Finally, we review software packages that implement strategies to efficiently deal with large bisulfite sequencing datasets locally and we discuss online analysis workflows that do not require any prior programming skills. The analysis strategies described in this review will guide researchers at any level to the best practices of bisulfite sequencing analysis.

The extracellular RNA complement of Escherichia coli.

The secretion of biomolecules into the extracellular milieu is a common and well-conserved phenomenon in biology. In bacteria, secreted biomolecules are not only involved in intra-species communication but they also play roles in inter-kingdom exchanges and pathogenicity. To date, released products, such as small molecules, DNA, peptides, and proteins, have been well studied in bacteria. However, the bacterial extracellular RNA complement has so far not been comprehensively characterized. Here, we have analyzed, using a combination of physical characterization and high-throughput sequencing, the extracellular RNA complement of both outer membrane vesicle (OMV)-associated and OMV-free RNA of the enteric Gram-negative model bacterium Escherichia coli K-12 substrain MG1655 and have compared it to its intracellular RNA complement. Our results demonstrate that a large part of the extracellular RNA complement is in the size range between 15 and 40 nucleotides and is derived from specific intracellular RNAs. Furthermore, RNA is associated with OMVs and the relative abundances of RNA biotypes in the intracellular, OMV and OMV-free fractions are distinct. Apart from rRNA fragments, a significant portion of the extracellular RNA complement is composed of specific cleavage products of functionally important structural noncoding RNAs, including tRNAs, 4.5S RNA, 6S RNA, and tmRNA. In addition, the extracellular RNA pool includes RNA biotypes from cryptic prophages, intergenic, and coding regions, of which some are so far uncharacterised, for example, transcripts mapping to the fimA-fimL and ves-spy intergenic regions. Our study provides the first detailed characterization of the extracellular RNA complement of the enteric model bacterium E. coli. Analogous to findings in eukaryotes, our results suggest the selective export of specific RNA biotypes by E. coli, which in turn indicates a potential role for extracellular bacterial RNAs in intercellular communication.

An alternative method for the evaluation of docking performance: RSR vs RMSD.

A new assessment criterion for docking poses is proposed in which experimental electron density is taken into account when evaluating the ability of docking programs to reproduce experimentally observed binding modes. Three docking programs (Gold, Glide, and Fred) were used to generate poses for a set of 88 protein-ligand complexes for which the crystal structure is known. The new criterion is based on the real space R-factor (RSR), which measures how well a group of atoms-in our case the ligand-fits the experimental electron density by comparing that density to the expected density, calculated from the model (i.e., the predicted ligand pose). The RSR-based measure is compared to the traditional criterion, the root-mean-square distance (RMSD) between the docking pose and the binding configuration in the crystallographic model. The results highlight several shortcomings of the RMSD criterion that do not affect the RSR-based measure. Examples illustrate that the RSR-derived approach allows a more meaningful a posteriori assessment of docking methods and results. Practical implications for docking evaluations and for methodological development work in this field are discussed.