RIP-seq reveals RNAs that interact with RNA polymerase and primary sigma factors in bacteria.

Bacteria have evolved structured RNAs that can associate with RNA polymerase (RNAP). Two of them have been known so far-6S RNA and Ms1 RNA but it is unclear if any other types of RNAs binding to RNAP exist in bacteria. To identify all RNAs interacting with RNAP and the primary σ factors, we have established and performed native RIP-seq in Bacillus subtilis, Corynebacterium glutamicum, Streptomyces coelicolor, Mycobacterium smegmatis and the pathogenic Mycobacterium tuberculosis. Besides known 6S RNAs in B. subtilis and Ms1 in M. smegmatis, we detected MTS2823, a homologue of Ms1, on RNAP in M. tuberculosis. In C. glutamicum, we discovered novel types of structured RNAs that associate with RNAP. Furthermore, we identified other species-specific RNAs including full-length mRNAs, revealing a previously unknown landscape of RNAs interacting with the bacterial transcription machinery.

Clinical experience of reoperative right ventricular outflow tract reconstruction with valved conduits: risk factors for conduit failure in long-term follow-up.

Reconstruction of right ventricular outflow tract in patients with congenital heart disease in various age groups remains a controversial issue. Currently, a little is known about the fate of secondary and subsequent conduit. The aim of the study was to determine risk factors of conduit failure, evaluate long-term conduit survival, find out which type of conduit should be preferred in case of reoperations. We performed a retrospective analysis of a total of 249 records of valved conduit secondary and subsequent replacement in right ventricular outflow tract in 197 patients. Median follow-up was 5.7 years. The study endpoints were defined as conduit explants; balloon dilatation of the graft (excluding balloon dilatation of left/right pulmonary artery), transcatheter pulmonary valve implantation; heart transplantation or death of the patient. There were total of 21 deaths (11% mortality) among 197 patients during the follow-up, 2 patients underwent heart transplant, in 23 implanted conduits pulmonary angioplasty or/including transcatheter pulmonary valve implantation was afterwards performed due to graft failure, conduit had to be explanted in 46 cases. After 28 years follow-up, freedom from graft failure after 5 years was 77%, 48% after 10 years and 21% after 15 years. Reoperative right ventricular outflow tract reconstruction demonstrates good mid-term and acceptable long-term outcomes regardless of the type of conduit implanted. Worse long-term graft survival of secondary and further conduits is associated with younger age of the recipient at implantation, small size of the conduit, younger age of donor and male donor in case of allograft implantation.

Retinitis pigmentosa-associated mutations in mouse Prpf8 cause misexpression of circRNAs and degeneration of cerebellar granule cells.

A subset of patients with retinitis pigmentosa (RP) carry mutations in several spliceosomal components including the PRPF8 protein. Here, we established two alleles of murine Prpf8 that genocopy or mimic aberrant PRPF8 found in RP patients-the substitution p.Tyr2334Asn and an extended protein variant p.Glu2331ValfsX15. Homozygous mice expressing the aberrant Prpf8 variants developed within the first 2 mo progressive atrophy of the cerebellum because of extensive granule cell loss, whereas other cerebellar cells remained unaffected. We further show that a subset of circRNAs were deregulated in the cerebellum of both Prpf8-RP mouse strains. To identify potential risk factors that sensitize the cerebellum for Prpf8 mutations, we monitored the expression of several splicing proteins during the first 8 wk. We observed down-regulation of all selected splicing proteins in the WT cerebellum, which coincided with neurodegeneration onset. The decrease in splicing protein expression was further pronounced in mouse strains expressing mutated Prpf8. Collectively, we propose a model where physiological reduction in spliceosomal components during postnatal tissue maturation sensitizes cells to the expression of aberrant Prpf8 and the subsequent deregulation of circRNAs triggers neuronal death.

Clinical prediction models for mortality in patients with covid-19: external validation and individual participant data meta-analysis.

OBJECTIVE: To externally validate various prognostic models and scoring rules for predicting short term mortality in patients admitted to hospital for covid-19. DESIGN: Two stage individual participant data meta-analysis. SETTING: Secondary and tertiary care. PARTICIPANTS: 46 914 patients across 18 countries, admitted to a hospital with polymerase chain reaction confirmed covid-19 from November 2019 to April 2021. DATA SOURCES: Multiple (clustered) cohorts in Brazil, Belgium, China, Czech Republic, Egypt, France, Iran, Israel, Italy, Mexico, Netherlands, Portugal, Russia, Saudi Arabia, Spain, Sweden, United Kingdom, and United States previously identified by a living systematic review of covid-19 prediction models published in The BMJ, and through PROSPERO, reference checking, and expert knowledge. MODEL SELECTION AND ELIGIBILITY CRITERIA: Prognostic models identified by the living systematic review and through contacting experts. A priori models were excluded that had a high risk of bias in the participant domain of PROBAST (prediction model study risk of bias assessment tool) or for which the applicability was deemed poor. METHODS: Eight prognostic models with diverse predictors were identified and validated. A two stage individual participant data meta-analysis was performed of the estimated model concordance (C) statistic, calibration slope, calibration-in-the-large, and observed to expected ratio (O:E) across the included clusters. MAIN OUTCOME MEASURES: 30 day mortality or in-hospital mortality. RESULTS: Datasets included 27 clusters from 18 different countries and contained data on 46 914patients. The pooled estimates ranged from 0.67 to 0.80 (C statistic), 0.22 to 1.22 (calibration slope), and 0.18 to 2.59 (O:E ratio) and were prone to substantial between study heterogeneity. The 4C Mortality Score by Knight et al (pooled C statistic 0.80, 95% confidence interval 0.75 to 0.84, 95% prediction interval 0.72 to 0.86) and clinical model by Wang et al (0.77, 0.73 to 0.80, 0.63 to 0.87) had the highest discriminative ability. On average, 29% fewer deaths were observed than predicted by the 4C Mortality Score (pooled O:E 0.71, 95% confidence interval 0.45 to 1.11, 95% prediction interval 0.21 to 2.39), 35% fewer than predicted by the Wang clinical model (0.65, 0.52 to 0.82, 0.23 to 1.89), and 4% fewer than predicted by Xie et al's model (0.96, 0.59 to 1.55, 0.21 to 4.28). CONCLUSION: The prognostic value of the included models varied greatly between the data sources. Although the Knight 4C Mortality Score and Wang clinical model appeared most promising, recalibration (intercept and slope updates) is needed before implementation in routine care.

Nicotinic Acetylcholine Receptors Expressed by Striatal Interneurons Inhibit Striatal Activity and Control Striatal-Dependent Behaviors.

Acetylcholine is an important modulator of striatal activity, and it is vital to controlling striatal-dependent behaviors, including motor and cognitive functions. Despite this significance, the mechanisms determining how acetylcholine impacts striatal signaling are still not fully understood. In particular, little is known about the role of nAChRs expressed by striatal interneurons. In the present study, we used FISH to determine which neuronal types express the most prevalent beta2 nicotinic subunit in the mouse striatum. Our data support a common view that nAChR expression is mostly restricted to striatal interneurons. Surprisingly though, cholinergic interneurons were identified as a population with the highest expression of beta2 nicotinic subunit. To investigate the functional significance of beta2-containing nAChRs in striatal interneurons, we deleted them by injecting the AAV-Cre vector into the striatum of beta2-flox/flox male mice. The deletion led to alterations in several behavioral domains, namely, to an increased anxiety-like behavior, decrease in sociability ratio, deficit in discrimination learning, and increased amphetamine-induced hyperlocomotion and c-Fos expression in mice with beta2 deletion. Further colocalization analysis showed that the increased c-Fos expression was present in both medium spiny neurons and presumed striatal interneurons. The present study concludes that, despite being relatively rare, beta2-containing nAChRs are primarily expressed in striatal neurons by cholinergic interneurons and play a significant role in behavior.SIGNIFICANCE STATEMENT A large variety of nAChRs are expressed in the striatum, a brain region that is crucial in the control of behavior. The complexity of receptors with different functions is hindering our understanding of mechanisms through which striatal acetylcholine modulates behavior. We focused on the role of a small population of beta2-containing nAChRs. We identified neuronal types expressing these receptors and determined their impact in the control of explorative behavior, anxiety-like behavior, learning, and sensitivity to stimulants. Additional experiments showed that these alterations were associated with an overall increased activity of striatal neurons. Thus, the small population of nicotinic receptors represents an interesting target for a modulation of response to stimulant drugs and other striatal-based behavior.

Response to elexacaftor/tezacaftor/ivacaftor in intestinal organoids derived from people with cystic fibrosis.

Superior efficacy of elexacaftor/tezacaftor/ivacaftor (ELX/TEZ/IVA) over tezacaftor/ivacaftor (TEZ/IVA) in people with cystic fibrosis (CF) and Phe508del/Phe508del genotype was shown in clinical trials. We utilized intestinal organoid approach to compare in vitro responses to these 2 CFTR modulator drug combinations and to check potential inter-individual variability in therapeutic response to the triple combination. Organoids from 17 subjects with Phe508del/Phe508del were screened with forskolin induced swelling assay. Significantly larger swelling, when exposed to ELX/TEZ/IVA as compared to TEZ/IVA, was observed in 16 of them. However, 1 sample showed no additional effect of ELX. The finding of unique CFTR variants in this sample indicates that genetic traits other than CF-causing CFTR mutation are worth exploring as they may have an impact on the definitive modulator drug response.

Fecal Microbiome Changes and Specific Anti-Bacterial Response in Patients with IBD during Anti-TNF Therapy.

Inflammatory bowel diseases (IBD) are chronic disorders of the gastrointestinal tract that have been linked to microbiome dysbiosis and immune system dysregulation. We investigated the longitudinal effect of anti-TNF therapy on gut microbiota composition and specific immune response to commensals in IBD patients. The study included 52 patients tracked over 38 weeks of therapy and 37 healthy controls (HC). To characterize the diversity and composition of the gut microbiota, we used amplicon sequencing of the V3V4 region of 16S rRNA for the bacterial community and of the ITS1 region for the fungal community. We measured total antibody levels as well as specific antibodies against assorted gut commensals by ELISA. We found diversity differences between HC, Crohn's disease, and ulcerative colitis patients. The bacterial community of patients with IBD was more similar to HC at the study endpoint, suggesting a beneficial shift in the microbiome in response to treatment. We identified factors such as disease severity, localization, and surgical intervention that significantly contribute to the observed changes in the gut bacteriome. Furthermore, we revealed increased IgM levels against specific gut commensals after anti-TNF treatment. In summary, this study, with its longitudinal design, brings insights into the course of anti-TNF therapy in patients with IBD and correlates the bacterial diversity with disease severity in patients with ulcerative colitis (UC).

Disease progression of 213 patients hospitalized with Covid-19 in the Czech Republic in March-October 2020: An exploratory analysis.

We collected a multi-centric retrospective dataset of patients (N = 213) who were admitted to ten hospitals in Czech Republic and tested positive for SARS-CoV-2 during the early phases of the pandemic in March-October 2020. The dataset contains baseline patient characteristics, breathing support required, pharmacological treatment received and multiple markers on daily resolution. Patients in the dataset were treated with hydroxychloroquine (N = 108), azithromycin (N = 72), favipiravir (N = 9), convalescent plasma (N = 7), dexamethasone (N = 4) and remdesivir (N = 3), often in combination. To explore association between treatments and patient outcomes we performed multiverse analysis, observing how the conclusions change between defensible choices of statistical model, predictors included in the model and other analytical degrees of freedom. Weak evidence to constrain the potential efficacy of azithromycin and favipiravir can be extracted from the data. Additionally, we performed external validation of several proposed prognostic models for Covid-19 severity showing that they mostly perform unsatisfactorily on our dataset.

Transcriptome sequencing and multi-plex imaging of prostate cancer microenvironment reveals a dominant role for monocytic cells in progression.

BACKGROUND: Prostate cancer is caused by genomic aberrations in normal epithelial cells, however clinical translation of findings from analyses of cancer cells alone has been very limited. A deeper understanding of the tumour microenvironment is needed to identify the key drivers of disease progression and reveal novel therapeutic opportunities. RESULTS: In this study, the experimental enrichment of selected cell-types, the development of a Bayesian inference model for continuous differential transcript abundance, and multiplex immunohistochemistry permitted us to define the transcriptional landscape of the prostate cancer microenvironment along the disease progression axis. An important role of monocytes and macrophages in prostate cancer progression and disease recurrence was uncovered, supported by both transcriptional landscape findings and by differential tissue composition analyses. These findings were corroborated and validated by spatial analyses at the single-cell level using multiplex immunohistochemistry. CONCLUSIONS: This study advances our knowledge concerning the role of monocyte-derived recruitment in primary prostate cancer, and supports their key role in disease progression, patient survival and prostate microenvironment immune modulation.

Probabilistic outlier identification for RNA sequencing generalized linear models.

Relative transcript abundance has proven to be a valuable tool for understanding the function of genes in biological systems. For the differential analysis of transcript abundance using RNA sequencing data, the negative binomial model is by far the most frequently adopted. However, common methods that are based on a negative binomial model are not robust to extreme outliers, which we found to be abundant in public datasets. So far, no rigorous and probabilistic methods for detection of outliers have been developed for RNA sequencing data, leaving the identification mostly to visual inspection. Recent advances in Bayesian computation allow large-scale comparison of observed data against its theoretical distribution given in a statistical model. Here we propose ppcseq, a key quality-control tool for identifying transcripts that include outlier data points in differential expression analysis, which do not follow a negative binomial distribution. Applying ppcseq to analyse several publicly available datasets using popular tools, we show that from 3 to 10 percent of differentially abundant transcripts across algorithms and datasets had statistics inflated by the presence of outliers.

Comparing quality of reporting between preprints and peer-reviewed articles in the biomedical literature.

BACKGROUND: Preprint usage is growing rapidly in the life sciences; however, questions remain on the relative quality of preprints when compared to published articles. An objective dimension of quality that is readily measurable is completeness of reporting, as transparency can improve the reader's ability to independently interpret data and reproduce findings. METHODS: In this observational study, we initially compared independent samples of articles published in bioRxiv and in PubMed-indexed journals in 2016 using a quality of reporting questionnaire. After that, we performed paired comparisons between preprints from bioRxiv to their own peer-reviewed versions in journals. RESULTS: Peer-reviewed articles had, on average, higher quality of reporting than preprints, although the difference was small, with absolute differences of 5.0% [95% CI 1.4, 8.6] and 4.7% [95% CI 2.4, 7.0] of reported items in the independent samples and paired sample comparison, respectively. There were larger differences favoring peer-reviewed articles in subjective ratings of how clearly titles and abstracts presented the main findings and how easy it was to locate relevant reporting information. Changes in reporting from preprints to peer-reviewed versions did not correlate with the impact factor of the publication venue or with the time lag from bioRxiv to journal publication. CONCLUSIONS: Our results suggest that, on average, publication in a peer-reviewed journal is associated with improvement in quality of reporting. They also show that quality of reporting in preprints in the life sciences is within a similar range as that of peer-reviewed articles, albeit slightly lower on average, supporting the idea that preprints should be considered valid scientific contributions.

Author Correction: Bacterial nanotubes as a manifestation of cell death.

An amendment to this paper has been published and can be accessed via a link at the top of the paper.

Bacterial nanotubes as a manifestation of cell death.

Bacterial nanotubes are membranous structures that have been reported to function as conduits between cells to exchange DNA, proteins, and nutrients. Here, we investigate the morphology and formation of bacterial nanotubes using Bacillus subtilis. We show that nanotube formation is associated with stress conditions, and is highly sensitive to the cells' genetic background, growth phase, and sample preparation methods. Remarkably, nanotubes appear to be extruded exclusively from dying cells, likely as a result of biophysical forces. Their emergence is extremely fast, occurring within seconds by cannibalizing the cell membrane. Subsequent experiments reveal that cell-to-cell transfer of non-conjugative plasmids depends strictly on the competence system of the cell, and not on nanotube formation. Our study thus supports the notion that bacterial nanotubes are a post mortem phenomenon involved in cell disintegration, and are unlikely to be involved in cytoplasmic content exchange between live cells.

rboAnalyzer: A Software to Improve Characterization of Non-coding RNAs From Sequence Database Search Output.

Searching for similar sequences in a database via BLAST or a similar tool is one of the most common bioinformatics tasks applied in general, and to non-coding RNAs in particular. However, the results of the search might be difficult to interpret due to the presence of partial matches to the database subject sequences. Here, we present rboAnalyzer - a tool that helps with interpreting sequence search result by (1) extending partial matches into plausible full-length subject sequences, (2) predicting homology of RNAs represented by full-length subject sequences to the query RNA, (3) pooling information across homologous RNAs found in the search results and public databases such as Rfam to predict more reliable secondary structures for all matches, and (4) contextualizing the matches by providing the prediction results and other relevant information in a rich graphical output. Using predicted full-length matches improves secondary structure prediction and makes rboAnalyzer robust with regards to identification of homology. The output of the tool should help the user to reliably characterize non-coding RNAs in BLAST output. The usefulness of the rboAnalyzer and its ability to correctly extend partial matches to full-length is demonstrated on known homologous RNAs. To allow the user to use custom databases and search options, rboAnalyzer accepts any search results as a text file in the BLAST format. The main output is an interactive HTML page displaying the computed characteristics and other context of the matches. The output can also be exported in an appropriate sequence and/or secondary structure formats.

MicroRNA dilution during oocyte growth disables the microRNA pathway in mammalian oocytes.

MicroRNAs (miRNAs) are ubiquitous small RNAs guiding post-transcriptional gene repression in countless biological processes. However, the miRNA pathway in mouse oocytes appears inactive and dispensable for development. We propose that marginalization of the miRNA pathway activity stems from the constraints and adaptations of RNA metabolism elicited by the diluting effects of oocyte growth. We report that miRNAs do not accumulate like mRNAs during the oocyte growth because miRNA turnover has not adapted to it. The most abundant miRNAs total tens of thousands of molecules in growing (∅ 40 µm) and fully grown (∅ 80 µm) oocytes, a number similar to that observed in much smaller fibroblasts. The lack of miRNA accumulation results in a 100-fold lower miRNA concentration in fully grown oocytes than in somatic cells. This brings a knock-down-like effect, where diluted miRNAs engage targets but are not abundant enough for significant repression. Low-miRNA concentrations were observed in rat, hamster, porcine and bovine oocytes, arguing that miRNA inactivity is not mouse-specific but a common mammalian oocyte feature. Injection of 250,000 miRNA molecules was sufficient to restore reporter repression in mouse and porcine oocytes, suggesting that miRNA inactivity comes from low-miRNA abundance and not from some suppressor of the pathway.

Meta-analysis of genotype-phenotype associations in Bardet-Biedl syndrome uncovers differences among causative genes.

Bardet-Biedl syndrome (BBS) is a recessive genetic disease causing multiple organ anomalies. Most patients carry mutations in genes encoding for the subunits of the BBSome, an octameric ciliary transport complex, or accessory proteins involved in the BBSome assembly or function. BBS proteins have been extensively studied using in vitro, cellular, and animal models. However, the molecular functions of particular BBS proteins and the etiology of the BBS symptoms are still largely elusive. In this study, we applied a meta-analysis approach to study the genotype-phenotype association in humans using our database of all reported BBS patients. The analysis revealed that the identity of the causative gene and the character of the mutation partially predict the clinical outcome of the disease. Besides their potential use for clinical prognosis, our analysis revealed functional differences of particular BBS genes in humans. Core BBSome subunits BBS2, BBS7, and BBS9 manifest as more critical for the function and development of kidneys than peripheral subunits BBS1, BBS4, and BBS8/TTC8, suggesting that incomplete BBSome retains residual function at least in the kidney.

Genexpi: a toolset for identifying regulons and validating gene regulatory networks using time-course expression data.

BACKGROUND: Identifying regulons of sigma factors is a vital subtask of gene network inference. Integrating multiple sources of data is essential for correct identification of regulons and complete gene regulatory networks. Time series of expression data measured with microarrays or RNA-seq combined with static binding experiments (e.g., ChIP-seq) or literature mining may be used for inference of sigma factor regulatory networks. RESULTS: We introduce Genexpi: a tool to identify sigma factors by combining candidates obtained from ChIP experiments or literature mining with time-course gene expression data. While Genexpi can be used to infer other types of regulatory interactions, it was designed and validated on real biological data from bacterial regulons. In this paper, we put primary focus on CyGenexpi: a plugin integrating Genexpi with the Cytoscape software for ease of use. As a part of this effort, a plugin for handling time series data in Cytoscape called CyDataseries has been developed and made available. Genexpi is also available as a standalone command line tool and an R package. CONCLUSIONS: Genexpi is a useful part of gene network inference toolbox. It provides meaningful information about the composition of regulons and delivers biologically interpretable results.

An Algorithm for Template-Based Prediction of Secondary Structures of Individual RNA Sequences.

While understanding the structure of RNA molecules is vital for deciphering their functions, determining RNA structures experimentally is exceptionally hard. At the same time, extant approaches to computational RNA structure prediction have limited applicability and reliability. In this paper we provide a method to solve a simpler yet still biologically relevant problem: prediction of secondary RNA structure using structure of different molecules as a template. Our method identifies conserved and unconserved subsequences within an RNA molecule. For conserved subsequences, the template structure is directly transferred into the generated structure and combined with de-novo predicted structure for the unconserved subsequences with low evolutionary conservation. The method also determines, when the generated structure is unreliable. The method is validated using experimentally identified structures. The accuracy of the method exceeds that of classical prediction algorithms and constrained prediction methods. This is demonstrated by comparison using large number of heterogeneous RNAs. The presented method is fast and robust, and useful for various applications requiring knowledge of secondary structures of individual RNA sequences.