A key aspect of RNA secondary structure prediction is the identification of novel functional elements. This is a challenging task because these elements typically are embedded in longer transcripts where the borders between the element and flanking regions have to be defined. The flanking sequences impact the folding of the functional elements both at the level of computational analyses and when the element is extracted as a transcript for experimental analysis. Here, we analyze how different flanking region lengths impact folding into a constrained structure by computing probabilities of folding for different sizes of flanking regions. Our method, RNAcop (RNA context optimization by probability), is tested on known and de novo predicted structures. In vitro experiments support the computational analysis and suggest that for a number of structures, choosing proper lengths of flanking regions is critical. RNAcop is available as web server and stand-alone software via http://rth.dk/resources/rnacop.

Fibrillarin (FBL) is a dual-function nucleolar protein that catalyzes 2'-O methylation of pre-rRNA and methylation of histone H2A at glutamine 104 (H2AQ104me). The mechanisms that regulate FBL activity are unexplored. Here, we show that FBL is acetylated at several lysine residues by the acetyltransferase CBP and deacetylated by SIRT7. While reversible acetylation does not impact FBL-mediated pre-rRNA methylation, hyperacetylation impairs the interaction of FBL with histone H2A and chromatin, thereby compromising H2AQ104 methylation (H2AQ104me) and rDNA transcription. SIRT7-dependent deacetylation of FBL ensures H2AQ104me and high levels of rRNA synthesis during interphase. At the onset of mitosis, nucleolar disassembly is accompanied by hyperacetylation of FBL, loss of H2AQ104me, and repression of polymerase I (Pol I) transcription. Overexpression of an acetylation-deficient, but not an acetylation-mimicking, FBL mutant restores H2AQ104me and transcriptional activity. The results reveal that SIRT7-dependent deacetylation impacts nucleolar activity by an FBL-driven circuitry that mediates cell-cycle-dependent fluctuation of rDNA transcription.

Twin-ribozyme introns contain a branching ribozyme (GIR1) followed by a homing endonuclease (HE) encoding sequence embedded in a peripheral domain of a group I splicing ribozyme (GIR2). GIR1 catalyses the formation of a lariat with 3 nt in the loop, which caps the HE mRNA. GIR1 is structurally related to group I ribozymes raising the question about how two closely related ribozymes can carry out very different reactions. Modelling of GIR1 based on new biochemical and mutational data shows an extended substrate domain containing a GoU pair distinct from the nucleophilic residue that dock onto a catalytic core showing a different topology from that of group I ribozymes. The differences include a core J8/7 region that has been reduced and is complemented by residues from the pre-lariat fold. These findings provide the basis for an evolutionary mechanism that accounts for the change from group I splicing ribozyme to the branching GIR1 architecture. Such an evolutionary mechanism can be applied to other large RNAs such as the ribonuclease P.

A knowledge of the positions of introns in eukaryotic genes is important for understanding the evolution of introns. Despite this, there has been relatively little focus on the distribution of intron positions in genes.

The 5' cap structure of eukaryotic mRNA is critical for its processing, transport, translation, and stability. The many functions of the cap and the fact that most, if not all, mRNA carries the same type of cap makes it difficult to analyze cap function in vivo at individual steps of gene expression. We have used the lariat capping ribozyme (LCrz) from the myxomycete Didymium to replace the mRNA m7G cap of a single reporter mRNA species with a tiny lariat in which the first and the third nucleotide are joined by a 2', 5' phosphodiester bond. We show that the ribozyme functions in vivo in the budding yeast Saccharomyces cerevisiae presumably without cofactors and that lariat capping occurs cotranscriptionally. The lariat-capped reporter mRNA is efficiently exported to the cytoplasm where it is found to be oligoadenylated and evenly distributed. Both the oligoadenylated form and a lariat-capped mRNA with a templated poly(A) tail translates poorly, underlining the critical importance of the m7G cap in translation. Finally, the lariat-capped RNA exhibits a threefold longer half-life compared to its m7G-capped counterpart, consistent with a key role for the m7G cap in mRNA turnover. Our study emphasizes important activities of the m7G cap and suggests new utilities of lariat capping as a molecular tool in vivo.

During zebrafish development, an early type of rRNA is gradually replaced by a late type that is substantially different in sequence. We applied RiboMeth-seq to rRNA from developmental stages for profiling of 2'-O-Me, to learn if changes in methylation pattern were a component of the shift. We compiled a catalog of 2'-O-Me sites and cognate box C/D guide RNAs comprising 98 high-confidence sites, including 10 sites that were not known from other vertebrates, one of which was specific to late-type rRNA. We identified a subset of sites that changed in methylation status during development and found that some of these could be explained by availability of their cognate SNORDs. Sites that changed during development were enriched in the novel sites revealed in zebrafish. We propose that the early type of rRNA is a specialized form and that its structure and ribose methylation pattern may be an adaptation to features of development, including translation of specific maternal mRNAs.

Convalescent plasma has been widely used to treat COVID-19 and is under investigation in numerous randomized clinical trials, but results are publicly available only for a small number of trials. The objective of this study was to assess the benefits of convalescent plasma treatment compared to placebo or no treatment and all-cause mortality in patients with COVID-19, using data from all available randomized clinical trials, including unpublished and ongoing trials (Open Science Framework, https://doi.org/10.17605/OSF.IO/GEHFX ).

Cancer cells are addicted to ribosome biogenesis and high levels of translation. Thus, differential inhibition of cancer cells can be achieved by targeting aspects of ribosome biogenesis or ribosome function. Using RiboMeth-seq for profiling of the ∼112 2'-O-Me sites in human ribosomal RNA, we demonstrated pronounced hypomethylation at several sites in patient-derived diffuse large B-cell lymphoma (DLBCL) cell lines with a more severe perturbation in ABC-DLBCL compared to GBC-DLBCL. We extended our analysis to tumor samples from patients and demonstrated significant changes to the ribosomal modification pattern that appeared to consist of cell growth-related as well as tumor-specific changes. Sites of hypomethylation in patient samples are discussed as potential drug targets, using as an example a site in the small subunit (SSU-C1440) located in a ribosomal substructure that can be linked to DLBCL pathogenesis.

The aim of this study was to provide information about immunity against COVID-19 along with risk factors and behavior among employees in day care facilities and preschools (DCS) in Denmark. In collaboration with the Danish Union of Pedagogues, during February and March 2021, 47,810 members were offered a point-of-care rapid SARS-CoV-2 antibody test (POCT) at work and were invited to fill in an electronic questionnaire covering COVID-19 exposure. Seroprevalence data from Danish blood donors (total Ig enzyme-linked immunosorbent assay [ELISA]) were used as a proxy for the Danish population. A total of 21,018 (45%) DCS employees completed the questionnaire and reported their POCT result {median age, 44.3 years (interquartile range [IQR], [32.7 to 53.6]); females, 84.1%}, of which 20,267 (96.4%) were unvaccinated and included in analysis. A total of 1,857 (9.2%) participants tested seropositive, significantly higher than a seroprevalence at 7.6% (risk ratio [RR], 1.2; 95% confidence interval [CI], 1.14 to 1.27) among 40,541 healthy blood donors (median age, 42 years [IQR, 28 to 53]; males, 51.3%). Exposure at work (RR, 2.9; 95% CI, 2.3 to 3.6) was less of a risk factor than exposure within the household (RR, 12.7; 95% CI, 10.2 to 15.8). Less than 25% of participants reported wearing face protection at work. Most of the participants expressed some degree of fear of contracting COVID-19 both at work and outside work. SARS-CoV-2 seroprevalence was slightly higher in DCS staff than in blood donors, but possible exposure at home was associated with a higher risk than at work. DCS staff expressed fear of contracting COVID-19, though there was limited use of face protection at work. IMPORTANCE Identifying at-risk groups and evaluating preventive interventions in at-risk groups is imperative for the ongoing pandemic as well as for the control of future epidemics. Although DCS staff have a much higher risk of being infected within their own household than at their workplace, most are fearful of being infected with COVID-19 or bringing COVID-19 to work. This represents an interesting dilemma and an important issue which should be addressed by public health authorities for risk communication and pandemic planning. This study design can be used in a strategy for ongoing surveillance of COVID-19 immunity or other infections in the population. The findings of this study can be used to assess the need for future preventive interventions in DCS, such as the use of personal protective equipment.

The emerging pathogen Campylobacter concisus has been isolated from patients with gastrointestinal diseases; however, it is also present in the gut of healthy individuals. The aim of this study was to compare IL-8 production in HT-29 cells after infection with C. concisus from different gastrointestinal disease phenotypes. Additionally, to investigate whether differentiation of isolates in genomospecies (GS1 and GS2) or presence of the zot gene, encoding the Zot toxin, affects IL-8 production. A total of 37 C. concisus isolates from patients with microscopic colitis (n = 20), ulcerative colitis (n = 5), Crohn's disease (n = 5), diarrhoea (n = 2) and from healthy controls (n = 5) were used. Intestinal HT-29 cells were infected and incubated for 24 h. Supernatants were subsequently removed and analysed for IL-8 by MILLIPLEX. All isolates were able to stimulate IL-8 production and IL-8 levels were higher than in non-infected HT-29 cells. No difference was observed between disease phenotypes or GS1 and GS2, whereas presence of the zot gene showed a tendency towards higher IL-8 production. Further investigations in other inflammatory and physiological models are needed to conclude whether C. concisus strains from different gastrointestinal disease phenotypes differ in pathogenic potential and play a part in gastrointestinal disease.

Campylobacter showae a bacterium historically linked to gingivitis and periodontitis, has recently been associated with inflammatory bowel disease and colorectal cancer. Our aim was to generate genome sequences for new clinical C. showae strains and identify functional properties explaining their pathogenic potential. Eight C. showae genomes were assessed, four strains isolated from inflamed gut tissues from paediatric Crohn's disease patients, three strains from colonic adenomas, and one from a gastroenteritis patient stool. Genome assemblies were analyzed alongside the only 3 deposited C. showae genomes. The pangenome from these 11 strains consisted of 4686 unique protein families, and the core genome size was estimated at 1050 ± 15 genes with each new genome contributing an additional 206 ± 16 genes. Functional assays indicated that colonic strains segregated into 2 groups: adherent/invasive vs. non-adherent/non-invasive strains. The former possessed Type IV secretion machinery and S-layer proteins, while the latter contained Cas genes and other CRISPR associated proteins. Comparison of gene profiles with strains in Human Microbiome Project metagenomes showed that gut-derived isolates share genes specific to tongue dorsum and supragingival plaque counterparts. Our findings indicate that C. showae strains are phenotypically and genetically diverse and suggest that secretion systems may play an important role in virulence potential.

Impairment of ribosome function activates the MAPKKK ZAK, leading to activation of mitogen-activated protein (MAP) kinases p38 and JNK and inflammatory signaling. The mechanistic basis for activation of this ribotoxic stress response (RSR) remains completely obscure. We show that the long isoform of ZAK (ZAKα) directly associates with ribosomes by inserting its flexible C terminus into the ribosomal intersubunit space. Here, ZAKα binds helix 14 of 18S ribosomal RNA (rRNA). An adjacent domain in ZAKα also probes the ribosome, and together, these sensor domains are critically required for RSR activation after inhibition of both the E-site, the peptidyl transferase center (PTC), and ribotoxin action. Finally, we show that ablation of the RSR response leads to organismal phenotypes and decreased lifespan in the nematode Caenorhabditis elegans (C. elegans). Our findings yield mechanistic insight into how cells detect ribotoxic stress and provide experimental in vivo evidence for its physiological importance.

The prediction of protein subcellular localization is of great relevance for proteomics research. Here, we propose an update to the popular tool DeepLoc with multi-localization prediction and improvements in both performance and interpretability. For training and validation, we curate eukaryotic and human multi-location protein datasets with stringent homology partitioning and enriched with sorting signal information compiled from the literature. We achieve state-of-the-art performance in DeepLoc 2.0 by using a pre-trained protein language model. It has the further advantage that it uses sequence input rather than relying on slower protein profiles. We provide two means of better interpretability: an attention output along the sequence and highly accurate prediction of nine different types of protein sorting signals. We find that the attention output correlates well with the position of sorting signals. The webserver is available at services.healthtech.dtu.dk/service.php?DeepLoc-2.0.

In bioinformatics, machine learning methods have been used to predict features embedded in the sequences. In contrast to what is generally assumed, machine learning approaches can also provide new insights into the underlying biology. Here, we demonstrate this by presenting TargetP 2.0, a novel state-of-the-art method to identify N-terminal sorting signals, which direct proteins to the secretory pathway, mitochondria, and chloroplasts or other plastids. By examining the strongest signals from the attention layer in the network, we find that the second residue in the protein, that is, the one following the initial methionine, has a strong influence on the classification. We observe that two-thirds of chloroplast and thylakoid transit peptides have an alanine in position 2, compared with 20% in other plant proteins. We also note that in fungi and single-celled eukaryotes, less than 30% of the targeting peptides have an amino acid that allows the removal of the N-terminal methionine compared with 60% for the proteins without targeting peptide. The importance of this feature for predictions has not been highlighted before.

The prediction of protein sub-cellular localization is an important step toward elucidating protein function. For each query protein sequence, LocTree2 applies machine learning (profile kernel SVM) to predict the native sub-cellular localization in 18 classes for eukaryotes, in six for bacteria and in three for archaea. The method outputs a score that reflects the reliability of each prediction. LocTree2 has performed on par with or better than any other state-of-the-art method. Here, we report the availability of LocTree3 as a public web server. The server includes the machine learning-based LocTree2 and improves over it through the addition of homology-based inference. Assessed on sequence-unique data, LocTree3 reached an 18-state accuracy Q18=80±3% for eukaryotes and a six-state accuracy Q6=89±4% for bacteria. The server accepts submissions ranging from single protein sequences to entire proteomes. Response time of the unloaded server is about 90 s for a 300-residue eukaryotic protein and a few hours for an entire eukaryotic proteome not considering the generation of the alignments. For over 1000 entirely sequenced organisms, the predictions are directly available as downloads. The web server is available at http://www.rostlab.org/services/loctree3.

To evaluate the positive predictive value (PPV) of ICD-10 diagnosis codes for brain abscess in the Danish National Patient Registry (DNPR).

Ribose methylation is one of the two most abundant modifications in human ribosomal RNA and is believed to be important for ribosome biogenesis, mRNA selectivity and translational fidelity. We have applied RiboMeth-seq to rRNA from HeLa cells for ribosome-wide, quantitative mapping of 2'-O-Me sites and obtained a comprehensive set of 106 sites, including two novel sites, and with plausible box C/D guide RNAs assigned to all but three sites. We find approximately two-thirds of the sites to be fully methylated and the remainder to be fractionally modified in support of ribosome heterogeneity at the level of RNA modifications. A comparison to HCT116 cells reveals similar 2'-O-Me profiles with distinct differences at several sites. This study constitutes the first comprehensive mapping of 2'-O-Me sites in human rRNA using a high throughput sequencing approach. It establishes the existence of a core of constitutively methylated positions and a subset of variable, potentially regulatory positions, and paves the way for experimental analyses of the role of variations in rRNA methylation under different physiological or pathological settings.

Recent experimental and computational progress has revealed a large potential for RNA structure in the genome. This has been driven by computational strategies that exploit multiple genomes of related organisms to identify common sequences and secondary structures. However, these computational approaches have two main challenges: they are computationally expensive and they have a relatively high false discovery rate (FDR). Simultaneously, RNA 3D structure analysis has revealed modules composed of non-canonical base pairs which occur in non-homologous positions, apparently by independent evolution. These modules can, for example, occur inside structural elements which in RNA 2D predictions appear as internal loops. Hence one question is if the use of such RNA 3D information can improve the prediction accuracy of RNA secondary structure at a genome-wide level. Here, we use RNAz in combination with 3D module prediction tools and apply them on a 13-way vertebrate sequence-based alignment. We find that RNA 3D modules predicted by metaRNAmodules and JAR3D are significantly enriched in the screened windows compared to their shuffled counterparts. The initially estimated FDR of 47.0% is lowered to below 25% when certain 3D module predictions are present in the window of the 2D prediction. We discuss the implications and prospects for further development of computational strategies for detection of RNA 2D structure in genomic sequence.

In all free-living organisms a late-stage checkpoint in the biogenesis of the small ribosomal subunit involves rRNA modification by an RsmA/Dim1 methyltransferase. The hyperthermophilic archaeon Nanoarchaeum equitans, whose existence is confined to the surface of a second archaeon, Ignicoccus hospitalis, lacks an RsmA/Dim1 homolog. We demonstrate here that the I. hospitalis host possesses the homolog Igni_1059, which dimethylates the N6-positions of two invariant adenosines within helix 45 of 16S rRNA in a manner identical to other RsmA/Dim1 enzymes. However, Igni_1059 is not transferred from I. hospitalis to N. equitans across their fused cell membrane structures and the corresponding nucleotides in N. equitans 16S rRNA remain unmethylated. An alternative mechanism for ribosomal subunit maturation in N. equitans is suggested by sRNA interactions that span the redundant RsmA/Dim1 site to introduce 2΄-O-ribose methylations within helices 44 and 45 of the rRNA.

Primary microcephaly (MCPH) is characterized by reduced brain size and intellectual disability. The exact pathophysiological mechanism underlying MCPH remains to be elucidated, but dysfunction of neuronal progenitors in the developing neocortex plays a major role. We identified a homozygous missense mutation (p.W155C) in Ribosomal RNA Processing 7 Homolog A, RRP7A, segregating with MCPH in a consanguineous family with 10 affected individuals. RRP7A is highly expressed in neural stem cells in developing human forebrain, and targeted mutation of Rrp7a leads to defects in neurogenesis and proliferation in a mouse stem cell model. RRP7A localizes to centrosomes, cilia and nucleoli, and patient-derived fibroblasts display defects in ribosomal RNA processing, primary cilia resorption, and cell cycle progression. Analysis of zebrafish embryos supported that the patient mutation in RRP7A causes reduced brain size, impaired neurogenesis and cell proliferation, and defective ribosomal RNA processing. These findings provide novel insight into human brain development and MCPH.