Transcriptional Responses to IFN-γ Require Mediator Kinase-Dependent Pause Release and Mechanistically Distinct CDK8 and CDK19 Functions.

Transcriptional responses to external stimuli remain poorly understood. Using global nuclear run-on followed by sequencing (GRO-seq) and precision nuclear run-on sequencing (PRO-seq), we show that CDK8 kinase activity promotes RNA polymerase II pause release in response to interferon-γ (IFN-γ), a universal cytokine involved in immunity and tumor surveillance. The Mediator kinase module contains CDK8 or CDK19, which are presumed to be functionally redundant. We implemented cortistatin A, chemical genetics, transcriptomics, and other methods to decouple their function while assessing enzymatic versus structural roles. Unexpectedly, CDK8 and CDK19 regulated different gene sets via distinct mechanisms. CDK8-dependent regulation required its kinase activity, whereas CDK19 governed IFN-γ responses through its scaffolding function (i.e., it was kinase independent). Accordingly, CDK8, not CDK19, phosphorylates the STAT1 transcription factor (TF) during IFN-γ stimulation, and CDK8 kinase inhibition blocked activation of JAK-STAT pathway TFs. Cytokines such as IFN-γ rapidly mobilize TFs to "reprogram" cellular transcription; our results implicate CDK8 and CDK19 as essential for this transcriptional reprogramming.

WAVES (Web-based tool for Analysis and Visualization of Environmental Samples)-a web application for visualization of wastewater pathogen sequencing results.

MOTIVATION: Environmental monitoring of pathogens provides an accurate and timely source of information for public health authorities and policymakers. In the last 2 years, wastewater sequencing proved to be an effective way of detection and quantification of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) variants circulating in population. Wastewater sequencing produces substantial amounts of geographical and genomic data. Proper visualization of spatial and temporal patterns in these data is crucial for the assessment of the epidemiological situation and forecasting. Here, we present a web-based dashboard application for the visualization and analysis of data obtained from sequencing of environmental samples. The dashboard provides multi-layered visualization of geographical and genomic data. It allows to display frequencies of detected pathogen variants as well as individual mutation frequencies. The features of WAVES (Web-based tool for Analysis and Visualization of Environmental Samples) for early tracking and detection of novel variants in the wastewater are demonstrated in an example of BA.1 variant and the signature Spike mutation S: E484A. WAVES dashboard is easily customized through the editable configuration file and can be used for different types of pathogens and environmental samples. AVAILABILITY AND IMPLEMENTATION: WAVES source code is freely available at https://github.com/ptriska/WavesDash under MIT license. A demo version of this application can be accessed at: https://wavesdashboard.azurewebsites.net/.

Genomic surveillance of SARS-CoV-2 evolution by a centralised pipeline and weekly focused sequencing, Austria, January 2021 to March 2023.

BackgroundThe COVID-19 pandemic was largely driven by genetic mutations of SARS-CoV-2, leading in some instances to enhanced infectiousness of the virus or its capacity to evade the host immune system. To closely monitor SARS-CoV-2 evolution and resulting variants at genomic-level, an innovative pipeline termed SARSeq was developed in Austria.AimWe discuss technical aspects of the SARSeq pipeline, describe its performance and present noteworthy results it enabled during the pandemic in Austria.MethodsThe SARSeq pipeline was set up as a collaboration between private and public clinical diagnostic laboratories, a public health agency, and an academic institution. Representative SARS-CoV-2 positive specimens from each of the nine Austrian provinces were obtained from SARS-CoV-2 testing laboratories and processed centrally in an academic setting for S-gene sequencing and analysis.ResultsSARS-CoV-2 sequences from up to 2,880 cases weekly resulted in 222,784 characterised case samples in January 2021-March 2023. Consequently, Austria delivered the fourth densest genomic surveillance worldwide in a very resource-efficient manner. While most SARS-CoV-2 variants during the study showed comparable kinetic behaviour in all of Austria, some, like Beta, had a more focused spread. This highlighted multifaceted aspects of local population-level acquired immunity. The nationwide surveillance system enabled reliable nowcasting. Measured early growth kinetics of variants were predictive of later incidence peaks.ConclusionWith low automation, labour, and cost requirements, SARSeq is adaptable to monitor other pathogens and advantageous even for resource-limited countries. This multiplexed genomic surveillance system has potential as a rapid response tool for future emerging threats.

ADAR-deficiency perturbs the global splicing landscape in mouse tissues.

Adenosine-to-inosine RNA editing and pre-mRNA splicing largely occur cotranscriptionally and influence each other. Here, we use mice deficient in either one of the two editing enzymes ADAR (ADAR1) or ADARB1 (ADAR2) to determine the transcriptome-wide impact of RNA editing on splicing across different tissues. We find that ADAR has a 100× higher impact on splicing than ADARB1, although both enzymes target a similar number of substrates with a large common overlap. Consistently, differentially spliced regions frequently harbor ADAR editing sites. Moreover, catalytically dead ADAR also impacts splicing, demonstrating that RNA binding of ADAR affects splicing. In contrast, ADARB1 editing sites are found enriched 5' of differentially spliced regions. Several of these ADARB1-mediated editing events change splice consensus sequences, therefore strongly influencing splicing of some mRNAs. A significant overlap between differentially edited and differentially spliced sites suggests evolutionary selection toward splicing being regulated by editing in a tissue-specific manner.

VirPool: model-based estimation of SARS-CoV-2 variant proportions in wastewater samples.

BACKGROUND: The genomes of SARS-CoV-2 are classified into variants, some of which are monitored as variants of concern (e.g. the Delta variant B.1.617.2 or Omicron variant B.1.1.529). Proportions of these variants circulating in a human population are typically estimated by large-scale sequencing of individual patient samples. Sequencing a mixture of SARS-CoV-2 RNA molecules from wastewater provides a cost-effective alternative, but requires methods for estimating variant proportions in a mixed sample. RESULTS: We propose a new method based on a probabilistic model of sequencing reads, capturing sequence diversity present within individual variants, as well as sequencing errors. The algorithm is implemented in an open source Python program called VirPool. We evaluate the accuracy of VirPool on several simulated and real sequencing data sets from both Illumina and nanopore sequencing platforms, including wastewater samples from Austria and France monitoring the onset of the Alpha variant. CONCLUSIONS: VirPool is a versatile tool for wastewater and other mixed-sample analysis that can handle both short- and long-read sequencing data. Our approach does not require pre-selection of characteristic mutations for variant profiles, it is able to use the entire length of reads instead of just the most informative positions, and can also capture haplotype dependencies within a single read.

A high resolution A-to-I editing map in the mouse identifies editing events controlled by pre-mRNA splicing.

Pre-mRNA-splicing and adenosine to inosine (A-to-I) RNA-editing occur mostly cotranscriptionally. During A-to-I editing, a genomically encoded adenosine is deaminated to inosine by adenosine deaminases acting on RNA (ADARs). Editing-competent stems are frequently formed between exons and introns. Consistently, studies using reporter assays have shown that splicing efficiency can affect editing levels. Here, we use Nascent-seq and identify ∼90,000 novel A-to-I editing events in the mouse brain transcriptome. Most novel sites are located in intronic regions. Unlike previously assumed, we show that both ADAR (ADAR1) and ADARB1 (ADAR2) can edit repeat elements and regular transcripts to the same extent. We find that inhibition of splicing primarily increases editing levels at hundreds of sites, suggesting that reduced splicing efficiency extends the exposure of intronic and exonic sequences to ADAR enzymes. Lack of splicing factors NOVA1 or NOVA2 changes global editing levels, demonstrating that alternative splicing factors can modulate RNA editing. Finally, we show that intron retention rates correlate with editing levels across different brain tissues. We therefore demonstrate that splicing efficiency is a major factor controlling tissue-specific differences in editing levels.

Adaptor protein RapZ activates endoribonuclease RNase E by protein-protein interaction to cleave a small regulatory RNA.

In Escherichia coli, endoribonuclease RNase E initiates degradation of many RNAs and represents a hub for post-transcriptional regulation. The tetrameric adaptor protein RapZ targets the small regulatory RNA GlmZ to degradation by RNase E. RapZ binds GlmZ through a domain located at the carboxyl terminus and interacts with RNase E, promoting GlmZ cleavage in the base-pairing region. When necessary, cleavage of GlmZ is counteracted by the homologous small RNA GlmY, which sequesters RapZ through molecular mimicry. In the current study, we addressed the molecular mechanism employed by RapZ. We show that RapZ mutants impaired in RNA-binding but proficient in binding RNase E are able to stimulate GlmZ cleavage in vivo and in vitro when provided at increased concentrations. In contrast, a truncated RapZ variant retaining RNA-binding activity but incapable of contacting RNase E lacks this activity. In agreement, we find that tetrameric RapZ binds the likewise tetrameric RNase E through direct interaction with its large globular domain within the catalytic amino terminus, independent of RNA. Although RapZ stimulates cleavage of at least one non-cognate RNA by RNase E in vitro, its activity is restricted to GlmZ in vivo as revealed by RNA sequencing, suggesting that certain features within the RNA substrate are also required for cleavage. In conclusion, RapZ boosts RNase E activity through interaction with its catalytic domain, which represents a novel mechanism of RNase E activation. In contrast, RNA-binding has a recruiting role, increasing the likelihood that productive RapZ/GlmZ/RNase E complexes form.

Dual RNA-seq unveils noncoding RNA functions in host-pathogen interactions.

Bacteria express many small RNAs for which the regulatory roles in pathogenesis have remained poorly understood due to a paucity of robust phenotypes in standard virulence assays. Here we use a generic 'dual RNA-seq' approach to profile RNA expression simultaneously in pathogen and host during Salmonella enterica serovar Typhimurium infection and reveal the molecular impact of bacterial riboregulators. We identify a PhoP-activated small RNA, PinT, which upon bacterial internalization temporally controls the expression of both invasion-associated effectors and virulence genes required for intracellular survival. This riboregulatory activity causes pervasive changes in coding and noncoding transcripts of the host. Interspecies correlation analysis links PinT to host cell JAK-STAT signalling, and we identify infection-specific alterations in multiple long noncoding RNAs. Our study provides a paradigm for a sensitive RNA-based analysis of intracellular bacterial pathogens and their hosts without physical separation, as well as a new discovery route for hidden functions of pathogen genes.

Functional characterization of the human tRNA methyltransferases TRMT10A and TRMT10B.

The TRM10 family of methyltransferases is responsible for the N1-methylation of purines at position 9 of tRNAs in Archaea and Eukarya. The human genome encodes three TRM10-type enzymes, of which only the mitochondrial TRMT10C was previously characterized in detail, whereas the functional significance of the two presumably nuclear enzymes TRMT10A and TRMT10B remained unexplained. Here we show that TRMT10A is m1G9-specific and methylates a subset of nuclear-encoded tRNAs, whilst TRMT10B is the first m1A9-specific tRNA methyltransferase found in eukaryotes and is responsible for the modification of a single nuclear-encoded tRNA. Furthermore, we show that the lack of G9 methylation causes a decrease in the steady-state levels of the initiator tRNAiMet-CAT and an alteration in its further post-transcriptional modification. Our work finally clarifies the function of TRMT10A and TRMT10B in vivo and provides evidence that the loss of TRMT10A affects the pool of cytosolic tRNAs required for protein synthesis.

An internal deletion of ADAR rescued by MAVS deficiency leads to a minute phenotype.

The RNA-editing protein ADAR is essential for early development in the mouse. Genetic evidence suggests that A to I editing marks endogenous RNAs as 'self'. Today, different Adar knockout alleles have been generated that show a common phenotype of apoptosis, liver disintegration, elevated immune response and lethality at E12.5. All the Adar knockout alleles can be rescued by a concomitant deletion of the innate immunity genes Mavs or Ifih1 (MDA5), albeit to different extents. This suggests multiple functions of ADAR. We analyze AdarΔ7-9 mice that show a unique growth defect phenotype when rescued by Mavs. We show that AdarΔ7-9 can form a truncated, unstable, editing deficient protein that is mislocalized. Histological and hematologic analysis of these mice indicate multiple tissue- and hematopoietic defects. Gene expression profiling shows dysregulation of Rps3a1 and Rps3a3 in rescued AdarΔ7-9. Consistently, a distortion in 40S and 60S ribosome ratios is observed in liver cells. This dysregulation is also seen in AdarΔ2-13; Mavs-/- but not in AdarE861A/E861A; Ifih1-/- mice, suggesting editing-independent functions of ADAR in regulating expression levels of Rps3a1 and Rps3a3. In conclusion, our study demonstrates the importance of ADAR in post-natal development which cannot be compensated by ADARB1.

Comparative Omics Analysis of Historic and Recent Isolates of Bordetella pertussis and Effects of Genome Rearrangements on Evolution.

Despite high vaccination coverage, pertussis is increasing in many industrialized countries, including the Czech Republic. To better understand Bordetella pertussis resurgence, we analyzed historic strains and recent clinical isolates by using a comparative omics approach. Whole-genome sequencing showed that historic and recent isolates of B. pertussis have substantial variation in genome organization and form separate phylogenetic clusters. Subsequent RNA sequence analysis and liquid chromatography with mass tandem spectrometry analyses showed that these variations translated into discretely separated transcriptomic and proteomic profiles. When compared with historic strains, recent isolates showed increased expression of flagellar genes and genes involved in lipopolysaccharide biosynthesis and decreased expression of polysaccharide capsule genes. Compared with reference strain Tohama I, all strains had increased expression and production of the type III secretion system apparatus. We detected the potential link between observed effects and insertion sequence element-induced changes in gene context only for a few genes.

Inosine induces context-dependent recoding and translational stalling.

RNA modifications are present in all classes of RNAs. They control the fate of mRNAs by affecting their processing, translation, or stability. Inosine is a particularly widespread modification in metazoan mRNA arising from deamination of adenosine catalyzed by the RNA-targeting adenosine deaminases ADAR1 or ADAR2. Inosine is commonly thought to be interpreted as guanosine by cellular machines and during translation. Here, we systematically test ribosomal decoding using mass spectrometry. We show that while inosine is primarily interpreted as guanosine it can also be decoded as adenosine, and rarely even as uracil. Decoding of inosine as adenosine and uracil is context-dependent. In addition, mass spectrometry analysis indicates that inosine causes ribosome stalling especially when multiple inosines are present in the codon. Indeed, ribosome profiling data from human tissues confirm inosine-dependent ribosome stalling in vivo. To our knowledge this is the first study where decoding of inosine is tested in a comprehensive and unbiased way. Thus, our study shows novel, unanticipated functions for inosines in mRNAs, further expanding coding potential and affecting translational efficiency.

iRAPs curb antisense transcription in E. coli.

RNA polymerase-binding RNA aptamers (RAPs) are natural RNA elements that control transcription in cis by directly contacting RNA polymerase. Many RAPs inhibit transcription by inducing Rho-dependent termination in Escherichia coli. Here, we studied the role of inhibitory RAPs (iRAPs) in modulation of antisense transcription (AT) using in silico and in vivo approaches. We revisited the antisense transcriptome in cells with impaired AT regulators (Rho, H-NS and RNaseIII) and searched for the presence of RAPs within antisense RNAs. Many of these RAPs were found at key genomic positions where they terminate AT. By exploring the activity of several RAPs both in a reporter system and in their natural genomic context, we confirmed their significant role in AT regulation. RAPs coordinate Rho activity at the antisense strand and terminate antisense transcripts. In some cases, they stimulated sense expression by alleviating ongoing transcriptional interference. Essentially, our data postulate RAPs as key determinants of Rho-mediated AT regulation in E. coli.

Signal transduction-dependent small regulatory RNA is involved in glutamate metabolism of the human pathogen Bordetella pertussis.

Bordetella pertussis is the causative agent of human whooping cough, a highly contagious respiratory disease which despite vaccination programs remains the major cause of infant morbidity and mortality. The requirement of the RNA chaperone Hfq for virulence of B. pertussis suggested that Hfq-dependent small regulatory RNAs are involved in the modulation of gene expression. High-throughput RNA sequencing revealed hundreds of putative noncoding RNAs including the RgtA sRNA. Abundance of RgtA is strongly decreased in the absence of the Hfq protein and its expression is modulated by the activities of the two-component regulatory system BvgAS and another response regulator RisA. Whereas RgtA levels were elevated under modulatory conditions or in the absence of bvg genes, deletion of the risA gene completely abolished RgtA expression. Profiling of the ΔrgtA mutant in the ΔbvgA genetic background identified the BP3831 gene encoding a periplasmic amino acid-binding protein of an ABC transporter as a possible target gene. The results of site-directed mutagenesis and in silico analysis indicate that RgtA base-pairs with the region upstream of the start codon of the BP3831 mRNA and thereby weakens the BP3831 protein production. Furthermore, our data suggest that the function of the BP3831 protein is related to transport of glutamate, an important metabolite in the B. pertussis physiology. We propose that the BvgAS/RisA interplay regulates the expression of RgtA which upon infection, when glutamate might be scarce, attenuates translation of the glutamate transporter and thereby assists in adaptation of the pathogen to other sources of energy.

Transcriptional profiling of human macrophages during infection with Bordetella pertussis.

Bordetella pertussis, a strictly human re-emerging pathogen and the causative agent of whooping cough, exploits a broad variety of virulence factors to establish efficient infection. Here, we used RNA sequencing to analyse the changes in gene expression profiles of human THP-1 macrophages resulting from B. pertussis infection. In parallel, we attempted to determine the changes in intracellular B. pertussis-specific transcriptomic profiles resulting from interaction with macrophages. Our analysis revealed that global gene expression profiles in THP-1 macrophages are extensively rewired 6 h post-infection. Among the highly expressed genes, we identified those encoding cytokines, chemokines, and transcription regulators involved in the induction of the M1 and M2 macrophage polarization programmes. Notably, several host genes involved in the control of apoptosis and inflammation which are known to be hijacked by intracellular bacterial pathogens were overexpressed upon infection. Furthermore, in silico analyses identified large temporal changes in expression of specific gene subsets involved in signalling and metabolic pathways. Despite limited numbers of the bacterial reads, we observed reduced expression of majority of virulence factors and upregulation of several transcriptional regulators during infection suggesting that intracellular B. pertussis cells switch from virulent to avirulent phase and actively adapt to intracellular environment, respectively.

Identification of RNA 3´ ends and termination sites in Haloferax volcanii.

Archaeal genomes are densely packed; thus, correct transcription termination is an important factor for orchestrated gene expression. A systematic analysis of RNA 3´ termini, to identify transcription termination sites (TTS) using RNAseq data has hitherto only been performed in two archaea, Methanosarcina mazei and Sulfolobus acidocaldarius. In this study, only regions directly downstream of annotated genes were analysed, and thus, only part of the genome had been investigated. Here, we developed a novel algorithm (Internal Enrichment-Peak Calling) that allows an unbiased, genome-wide identification of RNA 3´ termini independent of annotation. In an RNA fraction enriched for primary transcripts by terminator exonuclease (TEX) treatment we identified 1,543 RNA 3´ termini. Approximately half of these were located in intergenic regions, and the remainder were found in coding regions. A strong sequence signature consistent with known termination events at intergenic loci indicates a clear enrichment for native TTS among them. Using these data we determined distinct putative termination motifs for intergenic (a T stretch) and coding regions (AGATC). In vivo reporter gene tests of selected TTS confirmed termination at these sites, which exemplify the different motifs. For several genes, more than one termination site was detected, resulting in transcripts with different lengths of the 3´ untranslated region (3´ UTR).

Accurate mapping of tRNA reads.

Motivation: Many repetitive DNA elements are transcribed at appreciable expression levels. Mapping the corresponding RNA sequencing reads back to a reference genome is notoriously difficult and error-prone task, however. This is in particular true if chemical modifications introduce systematic mismatches, while at the same time the genomic loci are only approximately identical, as in the case of tRNAs. Results: We therefore developed a dedicated mapping strategy to handle RNA-seq reads that map to tRNAs relying on a modified target genome in which known tRNA loci are masked and instead intronless tRNA precursor sequences are appended as artificial 'chromosomes'. In a first pass, reads that overlap the boundaries of mature tRNAs are extracted. In the second pass, the remaining reads are mapped to a tRNA-masked target that is augmented by representative mature tRNA sequences. Using both simulated and real life data we show that our best-practice workflow removes most of the mapping artefacts introduced by simpler mapping schemes and makes it possible to reliably identify many of chemical tRNA modifications in generic small RNA-seq data. Using simulated data the FDR is only 2%. We find compelling evidence for tissue specific differences of tRNA modification patterns. Availability and implementation: The workflow is available both as a bash script and as a Galaxy workflow from https://github.com/AnneHoffmann/tRNA-read-mapping. Contact: fabian@tbi.univie.ac.at. Supplementary information: Supplementary data are available at Bioinformatics online.

The SmAP1/2 proteins of the crenarchaeon Sulfolobus solfataricus interact with the exosome and stimulate A-rich tailing of transcripts.

The conserved Sm and Sm-like proteins are involved in different aspects of RNA metabolism. Here, we explored the interactome of SmAP1 and SmAP2 of the crenarchaeon Sulfolobus solfataricus (Sso) to shed light on their physiological function(s). Both, SmAP1 and SmAP2 co-purified with several proteins involved in RNA-processing/modification, translation and protein turnover as well as with components of the exosome involved in 3΄ to 5΄ degradation of RNA. In follow-up studies a direct interaction with the poly(A) binding and accessory exosomal subunit DnaG was demonstrated. Moreover, elevated levels of both SmAPs resulted in increased abundance of the soluble exosome fraction, suggesting that they affect the subcellular localization of the exosome in the cell. The increased solubility of the exosome was accompanied by augmented levels of RNAs with A-rich tails that were further characterized using RNASeq. Hence, the observation that the Sso SmAPs impact on the activity of the exosome revealed a hitherto unrecognized function of SmAPs in archaea.

Comparative Integrated Omics Analysis of the Hfq Regulon in Bordetella pertussis.

Bordetella pertussis is a Gram-negative strictly human pathogen of the respiratory tract and the etiological agent of whooping cough (pertussis). Previously, we have shown that RNA chaperone Hfq is required for virulence of B. pertussis. Furthermore, microarray analysis revealed that a large number of genes are affected by the lack of Hfq. This study represents the first attempt to characterize the Hfq regulon in bacterial pathogen using an integrative omics approach. Gene expression profiles were analyzed by RNA-seq and protein amounts in cell-associated and cell-free fractions were determined by LC-MS/MS technique. Comparative analysis of transcriptomic and proteomic data revealed solid correlation (r2 = 0.4) considering the role of Hfq in post-transcriptional control of gene expression. Importantly, our study confirms and further enlightens the role of Hfq in pathogenicity of B. pertussis as it shows that Δhfq strain displays strongly impaired secretion of substrates of Type III secretion system (T3SS) and substantially reduced resistance to serum killing. On the other hand, significantly increased production of proteins implicated in transport of important metabolites and essential nutrients observed in the mutant seems to compensate for the physiological defect introduced by the deletion of the hfq gene.

Organ-wide profiling in mouse reveals high editing levels of Filamin B mRNA in the musculoskeletal system.

Adenosine to inosine RNA editing in protein-coding messenger RNAs (mRNAs) potentially leads to changes in the amino acid composition of the encoded proteins. The mRNAs encoding the ubiquitously expressed actin-crosslinking proteins Filamin A and Filamin B undergo RNA editing leading to a highly conserved glutamine to arginine exchange at the identical position in either protein. Here, by targeted amplicon sequencing we analysed the RNA editing of Filamin B across several mouse tissues during post-natal development. We find highest filamin B editing levels in skeletal muscles, cartilage and bones, tissues where Filamin B function seems most important. Through the analysis of Filamin B editing in mice deficient in either ADAR1 or 2, we identified ADAR2 as the enzyme responsible for Filamin B RNA editing. We show that in neuronal tissues Filamin B editing drops in spliced transcripts indicating regulated maturation of edited transcripts. We show further that the variability of Filamin B editing across several organs correlates with its mRNA expression.