We recently implemented a bioinformatics pipeline that can uncover novel, but rare, riboswitch candidates as well as other noncoding RNA structures in bacteria. A prominent candidate revealed by our initial search efforts was called the ‘thiS motif’ because of its frequent association with a gene coding for the ThiS protein, which delivers sulfur to form the thiazole moiety of the thiamin precursor HET-P. In the current report, we describe biochemical and genetic data demonstrating that thiS motif RNAs function as sensors of the thiamin precursor HMP-PP, which is fused with HET-P ultimately to form the final active coenzyme thiamin pyrophosphate (TPP). HMP-PP riboswitches exhibit a distinctive architecture wherein an unusually small ligand-sensing aptamer is almost entirely embedded within an otherwise classic intrinsic transcription terminator stem. This arrangement yields remarkably compact genetic switches that bacteria use to tune the levels of thiamin precursors during the biosynthesis of this universally distributed coenzyme.https://doi.org/10.7554/eLife.45210.001
Many bacteria use small genetic devices called riboswitches to sense molecules that are essential for life and regulate the genes necessary to make, break or move these molecules. Riboswitches are made of molecules of RNA and appear to have ancient origins that predate the evolution of bacteria and other lifeforms made of cells. Inside modern bacteria, chunks of DNA in the genome provide the instructions to make riboswitches and around 40 different types of riboswitch have been identified so far. However, it has been proposed that the instructions for thousands more riboswitches may remain hidden in the DNA of bacteria.
All of the currently known riboswitches contain a region called an aptamer that binds to a target molecule. This binding causes another structure in the riboswitch RNA to switch a specific gene on or off. For example, the aptamer binding might cause a hairpin-like structure called a terminator to form, which stops a gene being used to make new RNA molecules.
In 2019 a team of researchers reported using a computational approach to identify new riboswitches in bacteria. This approach identified many different chunks of DNA that might code for a riboswitch, including one known as the thiS motif. This potential new riboswitch appeared to be associated with a gene that encodes a protein required to make a vitamin called thiamin (also known as vitamin B1).
To test the new computational approach, Atilho et al. including several of the researchers involved in the earlier work used genetic and biochemical techniques to study the thiS motif. The experiments revealed that the motif binds to a molecule called HMP-PP, which bacteria use to make thiamin. Unexpectedly, the aptamer of the riboswitch was nested within a terminator, rather than being a separate entity.
The findings of Atilho et al. reveal that riboswitches can be even more compact than previously thought. Furthermore, these findings reveal new insights into how bacteria use riboswitches to manage their vitamin levels. In the future it may be possible to develop drugs that target such riboswitches to starve bacteria of these essential molecules.https://doi.org/10.7554/eLife.45210.002
Approximately 40 distinct riboswitch classes that regulate gene expression in various bacterial species have been experimentally validated to date (McCown et al., 2017; Serganov and Nudler, 2013; Sherwood and Henkin, 2016; Breaker, 2011). Based on the abundances and distributions of these known riboswitch classes, it has been proposed that many thousands of additional riboswitch classes remain to be discovered in the eubacterial domain of life (Ames and Breaker, 2010; Breaker, 2012; McCown et al., 2017). The collection of known riboswitch classes largely sense compounds or ions that are of fundamental importance to organisms from all three domains of life, and these ligands also exhibit a bias in favor of compounds (enzyme cofactors, RNA nucleotides and their precursors or derivatives) that are predicted to be of ancient origin (Breaker, 2012; McCown et al., 2017; Nelson and Breaker, 2017). If these trends hold, it seems likely that numerous additional riboswitch classes that regulate fundamental biological processes remain to be discovered. Unfortunately, the vast majority of these undiscovered riboswitch classes are predicted to be exceedingly rare, and this characteristic is likely to cause difficulties for researchers who seek to identify them.
To address this challenge, we developed a computational pipeline that first identifies the regions of bacterial genomes that are most likely to serve as transcription templates for structured noncoding RNAs (ncRNAs), and then uses comparative sequence and structural analyses to identify novel candidate RNA motifs (Stav et al., 2019). Specifically, this approach examines only the putative noncoding regions of a given sequenced bacterial genome, and evaluates each intergenic region (IGR) based on two parameters: (i) percent guanosine and cytidine (GC) nucleotide content and (ii) length in nucleotides. For many bacteria, structured ncRNAs are GC-rich compared to other regions of the genome (Klein et al., 2002; Schattner, 2002), and the IGRs that serve as synthesis templates for these ncRNAs tend to be much longer than typical bacterial IGRs that contain only an RNA polymerase promoter and/or a protein-specific regulatory domain. Our current bioinformatics pipeline is based on earlier implementations of this search strategy that were used to discover several novel structured ncRNA motifs (Meyer et al., 2009), including the SAM-V riboswitch class (Poiata et al., 2009).
Our updated computational pipeline was employed to comprehensively examine the genomes of five bacterial species, which revealed the existence of as many as 70 novel genetic elements, including 30 candidate ncRNA motifs (Stav et al., 2019). Of the eight candidate riboswitch classes uncovered in this search, the most promising was called the ‘thiS motif’ (Figure 1A) because representatives most commonly reside immediately upstream of thiS genes, which code for a protein that delivers sulfur to the pathway for the production of the thiamin biosynthetic intermediate 5-(2-hydroxyethyl)−4-methylthiazole phosphate (HET-P) (Begley et al., 2012). Other genes associated with thiS motif RNAs appear to code for proteins that participate in the production of HET-P or its fusion to 4-amino-5-hydroxymethyl-2-methylpyrimidine diphosphate (HMP-PP), to ultimately produce the bioactive coenzyme thiamin pyrophosphate (TPP) (Jurgenson et al., 2009).
An additional clue regarding the function of thiS motif RNAs was derived from the fact that this novel RNA structure occasionally resides in tandem with TPP riboswitches (Stav et al., 2019). Tandem arrangements of other riboswitch classes have been shown to function as two-input Boolean logic gates (Sudarsan et al., 2006; Stoddard and Batey, 2006; Lee et al., 2010; Sherlock et al., 2018), suggesting that the genes associated with tandem arrangements of TPP riboswitches and thiS motif RNAs likely respond to concentration changes of two distinct ligands. Together, these observations strongly indicate that thiS motif RNAs function as riboswitches that respond to a biochemical intermediate of the TPP biosynthetic pathway.
Previously (Stav et al., 2019), we created a riboswitch-reporter fusion construct by joining a thiS motif RNA representative to a β-galactosidase gene. Using this construct, we observed robust gene expression in host Bacillus subtilis cells with a deleted thiS gene (ΔthiS). This same reporter-fusion construct yields no reporter gene expression in host B. subtilis cells that naturally carry a thiS gene. These findings were consistent with our hypothesis that thiS motif RNAs function as riboswitches, but the precise ligand sensed by the unusual architecture of this RNA remained unknown.
In the current report, we describe a series of bioinformatic, genetic and biochemical analyses that provide conclusive evidence that the ligand for thiS motif RNAs is the TPP biosynthetic intermediate HMP-PP. Furthermore, our findings demonstrate that this unusually small riboswitch employs a distinct architecture to regulate RNA transcription termination. This regulatory RNA provides cells with an efficient mechanism to balance the production of two key biosynthetic intermediates, HMP-PP and HET-P, which are then fused to make the essential coenzyme TPP.
Most experimentally validated riboswitch classes are composed of distinct, but partially overlapping aptamer and expression platform domains (Barrick and Breaker, 2007; Breaker, 2012). In contrast, the thiS motif exhibits an unusual arrangement wherein the predominant secondary structure, derived by thermodynamic modeling, is an extended hairpin structure that maximizes conventional Watson/Crick base pairing (Figure 1A, left). This structure exhibits all the features characteristic of bacterial intrinsic terminator stems, including an uninterrupted and strong base-paired stem followed by a run of six or more uridine (U) residues (Wilson and von Hippel, 1995; Yarnell and Roberts, 1999). As a result, we concluded that transcription termination was certain to be a major function of thiS motif RNAs.
We frequently encounter simple terminator stems when using bioinformatics search algorithms to identify novel ncRNA candidates, and we have now begun to assign predicted terminator function to such motifs and then quickly move on to examine other more promising ncRNA candidates. When evaluating the thiS candidate, however, we noted four features that suggested this terminator stem was peculiar. First, the loop of the terminator hairpin is abnormally well conserved compared to the loop sequences of more typical terminator stems, which are usually irrelevant to the mechanism of transcription termination. This unique terminator element is found in species from two phyla and from several classes within Firmicutes (Supplementary file 1, Supplementary file 2). Second, another consensus structural model was also consistent with the comparative sequence analysis data (Figure 1A, right). This architecture, including a possible pseudoknot and two major base-paired stems, disrupts the contiguous terminator stem near the run of U residues, suggesting that mutually exclusive and competing structures and functions might exist for thiS motif RNAs. Third, these apparently specialized terminator stems associate exclusively with genes related to the biosynthesis and utilization of HET-P, which is a precursor of the coenzyme TPP (Figure 1B). These sequence, structure, and genomic distribution characteristics suggested to us that each thiS motif RNA might function as a compact ligand sensor and regulator of TPP coenzyme biosynthesis.
A fourth feature of thiS motif RNAs is that approximately 30% of the known representatives reside immediately downstream of riboswitches that sense and respond to TPP (Stav et al., 2019). Each associated TPP riboswitch appears to use a terminator stem as an expression platform. The tandem arrangement of two terminator stems for a single riboswitch aptamer would be unprecedented, and so this observation also supported our hypothesis that thiS motif RNAs represent an unusual form of riboswitch. Due to these tandem arrangements, we speculated that the natural ligand for this riboswitch candidate would not be TPP. Rather, it seemed more likely that the ligand would be a major precursor of TPP (either HET-P or HMP-PP), and that each tandem TPP-thiS motif system would function as a two-input Boolean logic gate (Sudarsan et al., 2006; Stoddard and Batey, 2006; Lee et al., 2010; Sherlock et al., 2018) to regulate HET-P production in response to the cellular concentrations of both TPP and one of its biosynthetic precursors.
Furthermore, if the two secondary structure states called the ‘terminator configuration’ and the ‘aptamer configuration’ (Figure 1A) represent the riboswitch ‘OFF’ and ‘ON’ configurations, respectively, then ligand binding by the aptamer in thiS motif RNAs is expected to activate gene expression. The vast majority of riboswitches that sense metabolites and control biosynthetic pathways are OFF switches, and thus the accumulation of the metabolite ligand results in reduced expression of the protein products that otherwise would make (or import) more of the desired metabolite. Although thiS motif representatives associate with genes for the production of HET-P (Figure 1B), it does not make sense that HET-P would turn on its own production when it is already abundant. Finally, we recognized that the thiE gene (coding for thiamin-phosphate synthase) is occasionally associated with thiS motif representatives, suggesting that the RNA motif determines if cellular conditions are suitable to couple the two key precursors of thiamin monophosphate (TMP) to eventually yield TPP. Taken together, all these bioinformatic observations are consistent with the hypothesis that thiS motif RNAs function as ON riboswitches for the TPP precursor HMP-PP.
As an initial test of our hypothesis that HMP-PP is the ligand for thiS motif RNAs, we employed a chromosomally-integrated reporter construct (Stav et al., 2019) carrying a β-galactosidase gene fused downstream of nucleotides encompassing the thiS motif from Clostridium species Maddingley (Rosewarne et al., 2013) (Figure 2A). Transformed B. subtilis cells carrying a transcriptional fusion of the wild-type (WT) riboswitch to a lacZ reporter gene, which was integrated into the amyE locus, exhibited no β-galactosidase activity in response to HMP added to rich (LB) liquid culture media (Figure 2B). Perhaps the natural suppression of TPP biosynthesis when cells have sufficient amounts of this coenzyme precludes the formation of excess HMP-PP in this surrogate organism. Specifically, if translation of the thiD gene coding for HMP/HMP-P kinase is largely suppressed under normal cellular conditions, the externally supplied HMP cannot be phosphorylated to generate HMP-PP.
In contrast, robust β-galactosidase activity is evident in cells co-transformed with the reporter construct and a vector overexpressing the ThiD protein, but only when growth media is supplemented with HMP (Figure 2B). Reporter gene expression levels are dependent on the concentration of HMP added to the culture medium (Figure 2C), suggesting that cells take up HMP and use the thiD gene product to add two phosphates to generate excess HMP-PP, which activates lacZ reporter expression by binding to the riboswitch.
Furthermore, the characteristics of a series of mutant riboswitch-reporter constructs examined in B. subtilis cells likewise indicates that thiS motif RNAs function as riboswitch aptamers for HMP-PP. Constructs carrying mutations M1 (G26C) or M2 (C35G) (Figure 2A), which alter strictly-conserved nucleotides that are predicted to be part of the pseudoknot structure of the aptamer (Figure 1A), are not activated by HMP addition (Figure 2C). Importantly, construct M3 (U46A, U47A, U48A) exhibits a detectable level of reporter gene expression in LB media alone, and this expression is further enhanced by the addition of HMP. This result suggests that these nucleotides are involved in forming a strong terminator stem, but that the identities of the nucleotides at these positions are not critical for ligand binding by the putative HMP-PP aptamer.
The WT riboswitch-reporter construct (Figure 2A) was also used to assess gene regulation in response to differences in growth media, and to several genetic disruptions of the TPP biosynthetic pathway (Figure 2—figure supplement 1). B. subtilis cells grown in LB medium are expected to suppress the genes needed to produce both the pyrimidine moiety HMP-PP and the thiazole moiety HET-P. Thus, as expected, expression of the reporter gene fused to the putative HMP-PP riboswitch is off, regardless of the genetic background tested. Likewise, WT cells grown in minimal (GMM) medium also do not express the reporter gene, presumably because HMP-PP does not accumulate due to its rapid and efficient conversion to the final product, TPP. By contrast, ΔthiS and ΔthiE cells (carrying genetic disruptions of the thiS and thiE genes, respectively) exhibit high levels of reporter gene expression (Figure 2—figure supplement 1). These two genetic knock-out strains are predicted to accumulate HMP-PP because they either lack the protein (ThiS) that initiates the production of HET-P, or they lack the protein (ThiE) that fuses HMP-PP to HET-P even when HET-P is available.
Additional support for the hypothesis that HMP-PP is the ligand, rather than HMP-P, is provided by the level of expression of the reporter construct when present in host cells carrying a deletion of the gene (ΔthiD) required to phosphorylate HMP-P to make HMP-PP. Reporter gene expression is not observed in ΔthiD cells (Figure 2—figure supplement 1C), despite the fact that these cells should accumulate HMP-P. Overall, these genetic results strongly indicate that HMP-PP is the ligand for a riboswitch class represented by thiS motif RNAs that turns on gene expression when this ligand is abundant.
The two possible architectures of the thiS motif RNA from C. sp. Maddingley (Figure 3A), also observed in other members of this candidate riboswitch class (Figure 1A), suggest a mechanism for gene regulation involving the mutually exclusive formation of an intrinsic terminator stem and its competing ligand-bound aptamer state. Such transcription control mechanisms for riboswitches have been experimentally validated in the past (e.g. McDaniel et al., 2003; Mironov et al., 2002; Sudarsan et al., 2003; Wickiser et al., 2005a) by using single-round in vitro transcription assays (Landick et al., 1996) to reveal ligand-dependent modulation of RNA transcript lengths. Therefore, we sought to examine the proposed riboswitch mechanism by conducting in vitro transcription reactions with the expected natural ligand, HMP-PP.
Unfortunately, our use of such transcription assays was made more difficult, because HMP-PP is not commercially available and is known to be relatively unstable (Hanes et al., 2007). Therefore, we freshly generated HMP-PP from HMP and ATP enzymatically (Figure 3B) by using recombinantly produced ThiD protein. HMP-PP production as described previously (Hanes et al., 2007) was confirmed by mass spectrum analysis (Figure 3C and Figure 3—figure supplement 1). The enzymatic reactions were deproteinized by filtration, and the resulting mixtures containing HMP-PP were added, without further purification, to various assays as described for each experiment.
The DNA template for the WT thiS motif RNA construct (Figure 3A) yields only ~14% full-length transcript when transcription reactions are conducted in the presence of HMP alone. In contrast, nearly 40% of the transcripts are full-length (FL) when transcription reactions are conducted in the presence of enzyme-prepared HMP-PP (Figure 3D). Furthermore, mutations M1 and M2, which disrupt possible pseudoknot formation by the aptamer and which eliminate gene expression in vivo (Figure 2C), fail to produce more FL transcripts when HMP-PP is present. These results suggest that HMP-PP binding induces RNA polymerase to transcribe past the intrinsic terminator element to yield fewer terminated (T) products and a greater proportion of FL RNAs.
From the assay data presented, this switching effect is incomplete in vitro, as evident by the fact that construct M4, which carries three nucleotide changes that disrupt the otherwise perfect base pairing of the terminator stem, yields nearly 100% ‘FL’ transcripts. Construct M5 yields transcripts that approximate the length of ‘T’ RNAs, which are generated naturally by termination due to the action of the intrinsic terminator stem. Incomplete switching by terminator-regulating riboswitches is typical for such assays conducted in vitro (e.g. McDaniel et al., 2003; Mironov et al., 2002; Sherlock et al., 2018; Sudarsan et al., 2003; Wickiser et al., 2005a). However, it seems likely that a larger dynamic range for transcription control is exploited by this riboswitch in cells, as evident by the robust differences in gene expression for the nearly identical reporter constructs used in this study (Figure 2 and Figure 2—figure supplement 1).
The unique architecture of thiS motif RNAs also served as an obstacle for evaluating the ability of these RNAs to directly bind a ligand. We frequently employ in-line probing assays (Soukup and Breaker, 1999; Regulski and Breaker, 2008) to determine if RNAs undergo structural changes in response to ligand binding. For thiS motif RNAs, the terminator stem is expected to dominate over the aptamer configuration (Figure 1A) simply due to the different number of base-pairs present in each structural state. Therefore, under thermodynamic equilibrium conditions typically experienced by RNAs subjected to in-line probing reactions, a full-length construct is not expected to bind HMP-PP because the RNA will always favor the terminator configuration (‘OFF’ state).
To address this problem, we reasoned that a shorter construct that weakens the terminator stem might permit the aptamer configuration to be adopted. Such constructs are also likely to better represent the structures naturally adopted by the riboswitch during transcription. An RNA polymerase paused at the run of U nucleotides at the end of the intrinsic terminator element will sequester ~12 nucleotides of the transcript within the protein structure (Monforte et al., 1990; Komissarova and Kashlev, 1998; Vassylyev et al., 2002). Therefore, only the first 50 to 54 nucleotides of the nascent transcript of the natural C. sp. Maddingley thiS motif should be exposed if the RNA polymerase complex is paused within the run of U nucleotides of the intrinsic terminator stem.
A series of four thiS motif RNA constructs was prepared to investigate the hypothesis that shorter constructs might permit aptamer formation by weakening the terminator stem. The longest RNA construct carries the full terminator stem (66 thiS), and three progressively truncated variants (54 thiS, 53 thiS and 52 thiS) represent only the RNA regions that would be exposed when RNA polymerase is stalled at various locations within the run of U nucleotides (Figure 4A). In-line probing assays (Figure 4B) reveal that 66 thiS indeed exclusively forms the terminator stem, as evident by robust spontaneous RNA cleavage in the unstructured regions near the 5′ terminus through nucleotide position 9, and by the unstructured loop region including nucleotides 35 through 39. This loop region remains unstructured with all constructs tested, demonstrating that the upper base-paired portion is common to all structural configurations.
Construct 54 thiS appears to retain all base-pairs of the terminator stem, except those that are disrupted by the deletion of the 3′-terminal nucleotides from position 54 and beyond (Figure 4B). In stark contrast, construct 52 thiS adopts a configuration that matches the consensus model for the aptamer configuration, including the formation of stems P1 and P2 (Figure 1A). Intriguingly, this shortest construct exhibits modest evidence (Figure 4B, asterisk) of a weakening of base-pairing within the series of G-U wobble interactions formed between nucleotides 24 to 27 and nucleotides 46 to 49. This structural flexibility might permit the formation of a pseudoknot or some other tertiary interactions between these highly-conserved G nucleotides and pyrimidines in the loop region upon ligand binding. Also noteworthy is the fact that the intermediate length construct, 53 thiS, exhibits an in-line probing pattern consistent with the formation of both the terminator configuration and the aptamer configuration. These results suggest that the genetic decision whether to form the terminator stem configuration or the ligand-bound aptamer configuration is made within a very narrow window of transcription progression.
Given the ability of the 52 thiS construct to exclusively adopt the aptamer configuration, we used this RNA to seek additional evidence for direct binding of HMP-PP by the riboswitch. However, this construct still failed to exhibit ligand-induced changes in the banding pattern resulting from in-line probing reactions. We speculate that nucleotide positions 24–27 and positions 46–49 form base-pairs that trap this portion of the motif in its terminator configuration. To further favor the formation of the desired aptamer configuration, the M3 version of the 52 thiS construct was made (Figure 5A). The same M3 mutations, made in the context of the full-length construct, earlier were observed to retain ligand responsiveness in riboswitch-reporter assays in vivo (Figure 2B).
The 52 thiS M3 construct indeed exhibits structural modulation upon introduction of HMP-PP (Figure 5B). Although the main base-paired regions P1 and P2 remain unchanged by HMP-PP addition, nucleotides involved in forming the putative pseudoknot (sites 1 and 2), and other nucleotides in the loop of P2 (site 3) appear to become more structured. Additional mutations were introduced into construct 52 thiS M3 to evaluate the effects of mutations known or expected to disrupt gene control function. Specifically, previous mutations M1 or M2 (Figure 2A), or M6 (Figure 5—figure supplement 1A) were introduced to create constructs, termed M7, M8 and M9, which carry alterations to highly-conserved nucleotides that presumably disrupt pseudoknot formation. These additional mutations eliminate structural modulation by HMP-PP (Figure 5—figure supplement 1B), as would be expected if these highly-conserved sequence and structural features are critical for riboswitch aptamer function.
We suspect that the ligand concentration is insufficient to saturate all RNAs in the sample. Regardless, by quantifying the band intensities at the three sites of structural modulation and by assuming that these values could become zero (or background) upon ligand binding, we generated a partial HMP-PP binding curve for 52 thiS M3 (Figure 5C), which is consistent with a 1-to-1 interaction between the ligand and the RNA aptamer. Given our use of enzymatically prepared HMP-PP samples, we cannot use this data to precisely determine the dissociation constant (KD) for this interaction. However, if we assume that all HMP was enzymatically converted into HMP-PP, this product was fully recovered after removal of protein by filtration, and that there was no degradation over the time frame of the assays, the KD value cannot be greater than 500 μM for this variant aptamer construct.
A series of shortened constructs carrying the M3 mutations (Figure 5—figure supplement 2A) recapitulate the transition between the terminator and aptamer configurations as originally observed for the series of WT truncated RNAs (Figure 4). Importantly, the 53 thiS M3 construct exhibits evidence of both binding to HMP-PP and of structural switching to favor the aptamer configuration (‘ON’ state) (Figure 5—figure supplement 2B). This finding suggests that it is necessary to prevent the formation of G-U wobble base-pairs between nucleotide positions 25–27 and positions 46–48 for ligand binding to be observed by using in-line probing assays. The time scale of in-line probing reactions (24 to 48 hr) is much longer than the time scale for the natural genetic decision to take place (probably a few seconds). Therefore our biochemical assays permit the RNAs to reach thermodynamic equilibrium (Wickiser et al., 2005a; Wickiser et al., 2005b), which should favor the terminator configuration. By mutating the U nucleotides at positions 46–48, we prevent constructs from becoming thermodynamically trapped in the terminator configuration during in-line probing, which permits HMP-PP binding to be observed. This same effect could be achieved with natural thiS motif sequences by having the U nucleotides at positions 46–48 remain sequestered in the RNA-exit channel (Vassylyev et al., 2002; Hein et al., 2014) of an RNA polymerase that is paused within the run of U nucleotides of the intrinsic terminator stem. We speculate that momentary ligand binding during this paused state will be sufficient to prevent terminator formation on the vastly shorter time scale that is relevant to the genetic decision process in cells.
Overall, our biochemical data are strongly consistent with the hypothesis that thiS motif RNAs function as riboswitches that directly bind to HMP-PP. As a result, we favor renaming thiS motif RNAs as HMP-PP riboswitches. The pursuit of more precise biochemical and biophysical characteristics, and further evidence for the mechanism proposed above, will require both the judicious use of RNA constructs that form the aptamer configuration, and either the preparation of pure samples of HMP-PP with consideration for its relative instability or the use of more stable analogs that can trigger riboswitch function.
As noted above, approximately 30% of the representatives of HMP-PP riboswitches reside in tandem with TPP riboswitches. In these tandem systems, the TPP riboswitch always occurs first, and the TPP aptamer is always associated with its own expression platform, which is routinely a readily recognizable intrinsic terminator stem. The HMP-PP riboswitch follows each complete TPP riboswitch, and by its unique architecture, carries its own intrinsic terminator expression platform.
This arrangement, as observed in a representative from the bacterium Clostridium lundense (Figure 6A), makes apparent the genetic decisions made by the host. Specifically, abundant TPP should repress expression of the associated gene, which usually codes for a protein needed to biosynthesize HET-P (Figure 1B). This repression, regardless of the concentration of HMP-PP, makes sense because cells do not need to make more TPP when this enzyme cofactor is already abundant. However, when TPP is in short supply, transcripts should read through the first intrinsic terminator stem to begin transcription of the HMP-PP riboswitch. The HMP-PP aptamer only permits transcription read-through of its co-resident intrinsic terminator stem if the relative abundance of its ligand is high. If HMP-PP is abundant, then the resulting full-length mRNAs will produce proteins that biosynthesize more HET-P. This two-input tandem riboswitch system has a truth table that matches a ‘converse nonimplication’ Boolean logic function (Figure 6B).
To confirm that a tandem TPP and HMP-PP riboswitch system functions with this Boolean logic, we created a tandem riboswitch-reporter fusion construct based on the natural C. lundense RNA (Figure 6A), and examined its function in B. subtilis cells. Unfortunately, we cannot easily create cellular conditions to fully examine all four possible ligand states represented in the truth table (Figure 6B). Most obviously, TPP is an essential cofactor, and thus we cannot deplete its concentration to zero while maintaining cell viability. Also, the various TPP biosynthetic intermediates are likely to be present in low concentrations if the cell is continuously making more TPP, even at low flux. As a result, a cell might experience conditions where an individual nascent RNA transcript of the tandem riboswitch system could respond to any of the four possible states, but these responses when summed over multiple nascent transcripts will yield a hybrid gene expression response that generally reflects the state of ligands in the cell.
To fully examine the gene expression outputs for all four possible states, despite the complications described above, we only partially manipulated the state of ligands, while also creating artificial representations of the missing states by using riboswitch aptamer mutations. Specifically, by adding HMP to the culture medium, we could change the ligand state from state 2 (+TPP, ‒ HMP-PP) to state 4 (+TPP,+HMP-PP). As predicted by the truth table, ligand states 2 and 4 do not yield expression of the reporter gene (Figure 6C). We also artificially created states 1 (‒ TPP, ‒ HMP-PP) and 3 (‒ TPP,+HMP-PP) by employing disruptive mutations in the aptamers for TPP (M10), HMP-PP (M1), or both (M11). Only the artificial state 3 condition, as created by mutating the TPP aptamer and growing cells in the presence of HMP, yields robust reporter gene expression (Figure 6C, M10). This gene expression effect is lost when disabling mutations are placed in both aptamers (M11) to simulate state 1. These findings are consistent with our hypothesis that each tandem TPP and HMP-PP riboswitch system functions as a Boolean converse nonimplication logic gate.
All of our findings derived from bioinformatic, genetic and biochemical analyses indicate that thiS motif RNAs represent an unusual class of riboswitches that sense and respond to the thiamin pyrophosphate biosynthetic intermediate HMP-PP. Its architecture is unusual because the majority of the nucleotides and structures required for ligand recognition appear to be encompassed by an intrinsic transcription terminator, which forms in a mutually-exclusive manner to regulate gene expression. All other riboswitch classes reported to date (McCown et al., 2017) that employ an intrinsic terminator structure as an expression platform carry the most-highly conserved nucleotides and structural features largely upstream of, and sometimes only partly overlapping, the terminator structure. This unique overlapping architecture of HMP-PP riboswitches provides bacteria with a highly-compact, yet effective RNA-based gene control system for measuring HMP-PP concentrations and tuning the production of HET-P to more efficiently produce TPP.
The embedded nature of the aptamer and expression platform domains, however, complicated the pursuit of experiments typically used to biochemically validate a newly discovered riboswitch class. Most importantly, our experiments required the use of non-native constructs to successfully reveal HMP-PP binding in vitro. Therefore, it is notable that modifications to certain key constructs maintain function in vitro or in vivo. For example, the trimming of nucleotides at the 3′ region of the terminator stem to create RNA construct 52 thiS yields a molecule that is expected to represent the exposed portion of a polymerase-paused nascent RNA transcript (Figure 4A). This truncated construct adopts the predicted P1 stem of the aptamer configuration (Figure 4B). Also, the addition of nucleotide changes to create the M3 variant of 52 thiS yields a construct that binds HMP-PP, whereas this same M3 mutation in a related construct retains ligand-dependent gene-control function in cells (Figure 2C).
Our results also reveal a mechanism for gene control involving the stalling of RNA polymerase at a precise location to yield an exposed nascent RNA transcript that is transiently capable of forming an HMP-PP binding pocket. The 52 thiS construct predominantly exists in the aptamer configuration, which is distinct from even modestly longer forms that partially (53 thiS) or exclusively (54 thiS) adopt the terminator configuration (Figure 4). If HMP-PP is present and docks in the aptamer when the nascent mRNA transcript is exposed up to nucleotide position 52, terminator stem formation is precluded and transcription proceeds to generate the full-length RNA. Forward movement of RNA polymerase by even a single nucleotide yields a construct that begins to commit to the terminator stem structure, and at this point transcription would be terminated regardless of the presence of HMP-PP.
The thiS motif was originally identified as a riboswitch candidate by using an updated computational pipeline (Stav et al., 2019) that focuses attention on relatively long, GC-rich IGRs. This method of riboswitch discovery was developed to improve our ability to uncover rare and/or small riboswitch candidates. Thousands of riboswitch classes are predicted to remain undiscovered in bacteria (Ames and Breaker, 2010; McCown et al., 2017), and many of these are likely to be exceedingly rare and/or difficult to find because of their structural simplicity. Experimental demonstration that thiS motif RNAs are representatives of a compact riboswitch class that responds to HMP-PP provides initial confirmation that the updated computational pipeline can uncover riboswitch classes that have resisted discovery by other bioinformatic or genetic methods. We are hopeful that many more riboswitch classes will be revealed by employing such bioinformatics search algorithms on additional sequenced bacterial genomes.
Finally, we note that the discovery of an HMP-PP riboswitch provides yet another example of a modern RNA class that selectively binds a molecule that was likely present in an RNA World (Benner et al., 1989; Gilbert, 1986). Many of the common enzyme cofactors are derived from RNA nucleotides or their precursors, including TPP and its precursor HMP-PP, which is a characteristic previously used to support the hypothesis that these molecules predate proteins (White, 1976). TPP and many other RNA-derived coenzymes might have participated in a complex metabolic state run entirely by RNA enzymes and receptors, and therefore many riboswitches for coenzymes likewise might be molecular descendants from RNAs that utilized these coenzymes long ago (Breaker, 2012; McCown et al., 2017). If true, HMP-PP riboswitches provide another opportunity to learn how simple RNA sequences and structures could have selectively bound important nucleotide-derived metabolites in primitive organisms of the RNA World.
Chemicals were purchased from Sigma-Aldrich with the exception of 4-amino-5-hydroxymethyl-2-methylpyrimidine, also known as (4-amino-2-methylpyrimidin-5-yl) methanol, which was purchased from Enamine Ltd. The radiolabeled molecules [γ-32P]-ATP and [α-32P]-ATP were purchased from PerkinElmer. All enzymes were purchased from New England BioLabs unless otherwise specified. Synthetic DNA oligonucleotides were purchased from Sigma-Aldrich or Integrated DNA Technologies. A list of oligonucleotides used in this study can be found in Supplementary file 3.
BL21 (DE3) E. coli cells were purchased from New England Biolabs and transformed with the appropriate plasmid for overproduction of the ThiD protein, whose enzymatic function ultimately was confirmed by monitoring the production of HMP-PP by mass spec analysis. The parent B. subtilis 168 (BGSC 1A1) strain, and the corresponding mutant strains ΔthiS (BGSC 11680), ΔthiE (BGSC 38290) and ΔthiD (BGSC 11710) were obtained from the Bacillus Genetic Stock Center (BGSC) at The Ohio State University, and genetically modified as described herein. All bacterial strains were verified by testing for the expected growth, antibiotic resistance, and reporter gene expression phenotypes.
Representative thiS motif RNAs used to create the consensus sequence and structure models (Figure 1A) were identified using Infernal 1.1 (Nawrocki and Eddy, 2013) from RefSeq version 80 and certain environmental microbial databases as described previously (Weinberg et al., 2017). A total of 400 unique and complete representatives were used to generate an updated consensus model relative to that published previously (Stav et al., 2019). The consensus sequence and structural models were derived by using the R2R algorithm (Weinberg and Breaker, 2011), which employs weighting. To prevent irrelevant differences between the aptamer and terminator configurations caused by weighting, the consensus model of the aptamer configuration was used to annotate the consensus sequence for the terminator configuration for the overlapping region.
Riboswitch-reporter constructs, integrating riboswitch representatives from either C. sp Maddingley or C. lundense DSM 17049, were prepared as synthetic oligonucleotides, amplified by PCR and cloned into vector pDG1661 upstream of the E. coli lacZ gene as described previously (Sudarsan et al., 2003; Nelson et al., 2017). Transcription initiation of the constructs is driven by the B. subtilis lysC gene promoter. The resulting WT and mutant reporter constructs were integrated into the amyE locus of WT (1A1 strain 168 Δtrp) or thiamin biosynthetic knockout strains (ΔthiS, ΔthiE or ΔthiD) as indicated. The resulting transformed strains were verified as previously described (Sherlock et al., 2018).
The thiD gene construct was generated by amplifying this gene from B. subtilis genomic DNA by PCR and inserted into the StuI site of a modified pDG148 vector using ligation-independent cloning as described previously (Joseph et al., 2001). The lacI gene in this vector has been mutated so that the thiD gene is expected to give constitutive expression. The resulting protein expression vector was then transformed into B. subtilis strains containing WT or mutant riboswitch reporter constructs, as indicated for each experiment.
Riboswitch-reporter assays were performed by inoculating various B. subtilis strains into Lysogeny Broth (LB) with appropriate antibiotics and growing overnight at 37°C. For liquid-culture reporter assays with thiamin biosynthetic knock-out strains, overnight cultures grown in LB were then diluted 1/20 into Spizizen glucose minimal medium (GMM) (Anagnostopoulos and Spizizen, 1961) and grown overnight at 37°C. The residual thiamin from LB is sufficient for growth in GMM over the duration of the assay. For riboswitch reporter experiments with ThiD-producing strains, bacteria were diluted directly into LB and grown overnight with or without supplementation with HMP. Liquid media (LB or GMM) was supplemented with X-gal (200 μg mL−1) to allow visual detection of reporter gene expression. Similarly, reporter expression analysis using 4-methylumbelliferyl β-D-galactopyranoside was conducted as described previously (Nelson et al., 2015; Atilho et al., 2019) to establish fluorescence units.
An N-terminal 6xHis-tagged thiD gene from Salmonella typhimurium was cloned into a pETDuet vector, which was then transformed into E. coli strain BL21(DE3). Transformed cells were grown in Terrific Broth medium until the OD600 reached 0.8. The resulting culture was incubated overnight at 16°C for protein expression that was induced by the addition of 0.3 mM isopropyl β-D-1-thiogalactopyranoside. Cells were pelleted, resuspended in Buffer A [50 mM Tris (pH 8 at 23°C), 400 mM NaCl, 10 mM imidazole, 5% glycerol, 0.1 mM tris(2-carboxyethyl)phosphine (TCEP)], and lysed with a microfluidizer. The resulting lysate was clarified via centrifugation, applied to NiNTA resin, and washed with 10 column volumes of Buffer A. Tagged protein was eluted from the column with three column volumes of Buffer B [50 mM Tris (pH 8 at 23°C), 400 mM NaCl, 400 mM imidazole, 5% glycerol, 0.1 mM TCEP] and applied to a size-exclusion column equilibrated in Buffer C [50 mM Tris (pH 8 at 23°C), 150 mM NaCl, 0.1 mM TCEP]. 10% Glycerol (v/v, final concentration) was added to the protein sample before storage at −80°C.
Pyrophosphorylation of HMP was carried out enzymatically as described previously (Hanes et al., 2007). Briefly, the reaction was initiated by adding 10 μL of Salmonella typhimurium HMP-P kinase (ThiD, 20 mg/mL) to a 90 μL reaction preparation to yield a final concentration of 5 mM HMP, 20 mM ATP, 50 mM Tris-HCl buffer (pH 7.5 at 23°C), 2 mM TCEP, and 5 mM MgCl2. This 100 μL reaction mixture was allowed to incubate at 23°C overnight, at which time the protein was removed by ultrafiltration using an Amicon Ultra-0.5 centrifugal filter unit with a 3 kDa cutoff membrane. The HMP-PP generated in this manner was used without further purification, typically on the same day but no later than 2 weeks after production. HMP-PP stock solutions were stored at −20°C. Control assays for transcription termination and in-line probing conducted without the addition of enzymatically prepared HMP-PP contained an equivalent amount of the enzymatic reaction preparation wherein the ThiD protein was excluded. These control assays therefore contain HMP, whereas test reactions contain HMP-PP.
Enzymatically prepared HMP-PP samples were sent to the MS and Proteomics Resource at Yale University for analysis. The presence of HMP-PP was confirmed by data (Figure 3C—figure supplement 1) generated using a Thermo Scientific LTQ Orbitrap ELITE mass spectrometer. Data was acquired and analyzed with Xcalibur (v2.1). Peaks were considered to have the same mass-to-charge ratio as HMP-PP if they were within 10 ppm of the calculated ratio.
The protocol used for single-round in vitro transcription assays was adapted from that described previously (Landick et al., 1996). DNA constructs were designed to include the promoter sequence of the lysC gene from B. subtilis, the riboswitch aptamer, and the expression platform of the thiS gene from C. sp Maddignly to 33 nucleotides following the terminator stem. Additional non-native nucleotides were added to the 5′ region upstream of the HMP-PP aptamer to increase the amount of [α-32P]-ATP incorporation.
To assemble each in vitro transcription reaction, approximately 2 pmol of the purified, PCR amplified DNA template was added to a transcription initiation mixture [final concentration of 20 mM Tris (pH 8.0 at 23°C), 20 mM NaCl, 14 mM MgCl2, 100 µM EDTA, 10 µg mL−1 bovine serum albumin, 130 µM ApA dinucleotide, 1% glycerol, 0.04 U µL−1 E. coli RNA polymerase holoenzyme, 2.5 µM GTP, 2.5 µM UTP, and 1 µM ATP]. Approximately 8 µCi [α-32P]-ATP was added to the 90 µL transcription reaction and transcription was allowed to proceed at 37°C for 30 min, leading to formation of a stalled polymerase complex before the first cytidine of each transcript. The reaction mixture was then distributed in 8 µL aliquots into separate microfuge tubes, which contained 1 µL of an HMP or HMP-PP solution plus 1 µL of 10x elongation buffer (200 mM Tris [pH 8.0 at 23°C], 200 mM NaCl, 140 mM MgCl2, 1 mM EDTA, 1 mg mL−1 heparin, 1.5 mM each of ATP, GTP, and CTP, and 0.5 mM UTP). For transcription termination assays, the maximum concentration of HMP-PP is 1.5 mM, assuming 100% enzymatic conversion of HMP and no loss to instability. However, HMP-PP concentration is likely substantially lower due to incomplete conversion of HMP to HMP-PP (Figure 3—figure supplement 1). Transcription elongation was allowed to proceed for 45 min at 37°C.
The transcription products were separated by denaturing (8 M urea) 10% polyacrylamide gel electrophoresis (PAGE) then imaged and quantified using a Typhoon Phosphorimager and ImageQuaNT software. The fraction of full length (FL) and terminated (T) RNA transcripts was calculated by measuring band intensity values and using the equation Fraction FL = (FL intensity)/(FL intensity +T intensity). The differences in specific activities between the FL and T products due to [α-32P]-ATP incorporation were considered negligible.
RNAs were prepared by in vitro transcription using DNA oligonucleotides containing a T7 RNA polymerase promoter sequence upstream of the desired template sequence. The resulting desired RNA transcripts were purified, enzymatically 5′ 32P-labeled, and repurified as previously described (Mirihana Arachchilage et al., 2018; Atilho et al., 2019).
In-line probing assays (Soukup and Breaker, 1999; Regulski and Breaker, 2008) were performed precisely as described previously (Atilho et al., 2019; Mirihana Arachchilage et al., 2018). The maximum HMP-PP concentration (max) is achieved by using a 3/10 dilution of enzymatically prepared HMP-PP as described above. The control reaction was performed using the 1/10 dilution of the control sample that lacks HMP-PP.
The Chemical Biology of Nucleic AcidsBacterial riboswitch discovery and analysis, The Chemical Biology of Nucleic Acids, Chichester, UK, Wiley, 10.1002/9780470664001.ch20.
Thiamin biosynthesis: still yielding fascinating biological chemistryBiochemical Society Transactions 40:555–560.https://doi.org/10.1042/BST20120084
Thiamin phosphate synthase: the rate of pyrimidine carbocation formationJournal of the American Chemical Society 129:4860–4861.https://doi.org/10.1021/ja0679634
RNA polymerase pausing and nascent-RNA structure formation are linked through clamp-domain movementNature Structural & Molecular Biology 21:794–802.https://doi.org/10.1038/nsmb.2867
The structural and biochemical foundations of thiamin biosynthesisAnnual Review of Biochemistry 78:569–603.https://doi.org/10.1146/annurev.biochem.78.072407.102340
Infernal 1.1: 100-fold faster RNA homology searchesBioinformatics 29:2933–2935.https://doi.org/10.1093/bioinformatics/btt509
Searching for RNA genes using base-composition statisticsNucleic Acids Research 30:2076–2082.https://doi.org/10.1093/nar/30.9.2076
Riboswitch-Mediated gene regulation: novel RNA architectures dictate gene expression responsesAnnual Review of Microbiology 70:361–374.https://doi.org/10.1146/annurev-micro-091014-104306
An mRNA structure in bacteria that controls gene expression by binding lysineGenes & Development 17:2688–2697.https://doi.org/10.1101/gad.1140003
Coenzymes as fossils of an earlier metabolic stateJournal of Molecular Evolution 7:101–104.https://doi.org/10.1007/BF01732468
The kinetics of ligand binding by an adenine-sensing riboswitchBiochemistry 44:13404–13414.https://doi.org/10.1021/bi051008u
Ailong KeReviewing Editor; Cornell University
Gisela StorzSenior Editor; National Institute of Child Health and Human Development, United States
Jinwei ZhangReviewer; NIDDK, United States
In the interests of transparency, eLife includes the editorial decision letter and accompanying author responses. A lightly edited version of the letter sent to the authors after peer review is shown, indicating the most substantive concerns; minor comments are not usually included.
Thank you for submitting your article "A bacterial riboswitch class for the thiamin precursor HMP-PP employs a terminator-embedded aptamer" for consideration by eLife. Your article has been reviewed by four peer reviewers, one of whom is a guest Reviewing Editor, and the evaluation has been overseen Gisela Storz as the Senior Editor. The following individual involved in review of your submission has agreed to reveal their identity: Jinwei Zhang (Reviewer #4).
The reviewers have discussed the reviews with one another and the Reviewing Editor has drafted this decision to help you prepare a revised submission.
The manuscript by Atilho et al. presents a new twist in the riboswitch field: an unexpected finding of a ligand-sensing domain embedded in the expression platform, a transcriptional terminator. It begs the question: how many of these mini-systems are hidden out there? The authors did excellent work proving that thiS motif is a functional riboswitch that responds to an important precursor of thiamine pyrophosphate, the active form of vitamin B1. The size of this riboswitch is remarkably small, and yet not only can this riboswitch recognizes the ligand very specifically, it can also function in conjunction with another riboswitch to precisely balance the expression of relevant genes. Undoubtedly, this work will attract attention of many researchers in different fields and will force researchers to re-think composition and role of one of the most widely-spread molecular elements – transcription terminator. The results will be of great interest to broad readership of eLife. Collectively, the reviewers have a few major points for clarification.
(The authors do not necessarily need to address every question raised here, however, an effort to provide a more completely picture of the thiS-thiM regulated bacterial physiology would certainly further improve the impact of the work.)
1) There is a lack of sufficient information on the distribution and genetic arrangements of thiS riboswitches. The authors cited an unpublished paper by Stav et al., 2019 but did not provide this manuscript and did not reveal details relevant to the current work. It is unclear when the Stav et al. work will be published and it is recommended to give a better summary of the findings and perhaps a supplementary figure with genetic arrangements, which would nicely complement a supplementary figure with riboswitch sequences. What are "uncommon" genes that reside downstream of thiS motif? Are all of them mentioned in the third paragraph of the Introduction? Can the assigned specificity of this motif explain regulation of these uncommon genes? If no, could it be possible for the riboswitch to recognize more than one ligand, especially given the fact that HMP-PP binding may be weak? Can this motif be a family of distinct riboswitches? Were there any attempts made to separate the motif on submotifs based on the genetic association? What are transporter genes controlled by the motif? ThiW appears to be an essential thiazole transporter and although the presence of the thiS motif controlling this protein suggests the regulation of coupling of the two TPP moieties, one can imagine direct regulation of a subtype of this motif by precursors of HET-P (Figure 1B). Figure 1B shows several genes associated with this motif and it would be helpful to identify and discuss these genes. Is each of these genes preceded by its own copy of the thiS motif? Are these genes controlled by a single thiS motif and not by a tandem this-thiM motif? The third paragraph of the subsection “The unusual architecture and genetic distribution of the thiS riboswitch candidate”, says that 30% of the known representatives reside near TPP riboswitches. Is there a genetic bias for the tandem riboswitch arrangement? Why the consensus motif (Figure 1A) is presented differently in the terminator and aptamer configurations? In the aptamer configuration, the loop in P1 is shown as variable while no extra nts are indicated for this region in the terminator configuration. The variable loop size contradicts the statement: "This structure exhibits all the features characteristic of bacterial intrinsic terminator stems, including an uninterrupted and strong base-paired stem". Actually the entire 5' end is different in the two configurations. Answers to these questions will help to convince this reviewer that the thiS motif is indeed a HMP-PP binder and will also help microbiologists to better understanding regulatory networks associated with vitamin production in bacteria.
2) The novel computational pipeline the authors describe is from unpublished work. Though not the subject of this manuscript, it is hard to see how this new approach differs from approaches used by this lab in the past. What makes this new approach more sensitive to rare riboswitch elements? Does this sensitivity arise from looking at smaller groups of bacteria (5 species) instead of across numerous species? The novelty of the approach is touted throughout, but not well explained.https://doi.org/10.7554/eLife.45210.019
1) There is a lack of sufficient information on the distribution and genetic arrangements of thiS riboswitches. The authors cited an unpublished paper by Stav et al., 2019 but did not provide this manuscript and did not reveal details relevant to the current work. It is unclear when the Stav et al. work will be published and it is recommended to give a better summary of the findings and perhaps a supplementary figure with genetic arrangements, which would nicely complement a supplementary figure with riboswitch sequences.
What are "uncommon" genes that reside downstream of thiS motif? Are all of them mentioned in the third paragraph of the Introduction? Can the assigned specificity of this motif explain regulation of these uncommon genes? If no, could it be possible for the riboswitch to recognize more than one ligand, especially given the fact that HMP-PP binding may be weak? Can this motif be a family of distinct riboswitches? Were there any attempts made to separate the motif on submotifs based on the genetic association? What are transporter genes controlled by the motif? ThiW appears to be an essential thiazole transporter and although the presence of the thiS motif controlling this protein suggests the regulation of coupling of the two TPP moieties, one can imagine direct regulation of a subtype of this motif by precursors of HET-P (Figure 1B). Figure 1B shows several genes associated with this motif and it would be helpful to identify and discuss these genes. Is each of these genes preceded by its own copy of the thiS motif? Are these genes controlled by a single thiS motif and not by a tandem this-thiM motif? The third paragraph of the subsection “The unusual architecture and genetic distribution of the thiS riboswitch candidate”, says that 30% of the known representatives reside near TPP riboswitches. Is there a genetic bias for the tandem riboswitch arrangement? Why the consensus motif (Figure 1A) is presented differently in the terminator and aptamer configurations? In the aptamer configuration, the loop in P1 is shown as variable while no extra nts are indicated for this region in the terminator configuration. The variable loop size contradicts the statement: "This structure exhibits all the features characteristic of bacterial intrinsic terminator stems, including an uninterrupted and strong base-paired stem". Actually the entire 5' end is different in the two configurations. Answers to these questions will help to convince this reviewer that the thiS motif is indeed a HMP-PP binder and will also help microbiologists to better understanding regulatory networks associated with vitamin production in bacteria.
Point 1A: The reviewers note a lack of bioinformatics information on the thiS motif that is otherwise typically supplied for novel riboswitches.
As the reviewers correctly speculate, this comprehensive information set was included for publication in an earlier manuscript (Stav et al.) reporting the bioinformatics discovery of the thiS riboswitch candidate (and numerous other RNA motifs). Unfortunately, this manuscript had experienced a remarkably long delay in the initial review process (~9 months). On March 7, 2019 we received notice that this original discovery manuscript has been accepted for publication by BMC Microbiology. Regardless, we have also decided to include an updated and detailed bioinformatics dataset as a Supplementary file 1 to the current revised manuscript for eLife. The original two supplementary files have been renamed accordingly. See also our response to Major Point 2 below.
Point 1B: The reviewers encourage us to comment on uncommon gene associations for the thiS motif and to discuss the possibility that the riboswitch aptamer might recognize different ligands to regulate these uncommon genes.
This question from the reviewers is an important one, particularly given the recent and striking examples of variant “ykkC orphan riboswitch” RNAs that have adapted to respond to at least five distinct ligands. However, nearly all the genes associated with thiS motif RNAs appear to code for proteins that participate in the production of HET-P or its coupling to HMP-PP. In some instances, associated genes from organisms whose genomes have not been annotated appear as ‘proteins of unknown function’. However, we have examined these genes and nearly all of these can easily be classified as thiS. Therefore, we do not have any indication that thiS motif RNAs have changed their specificity to regulate genes unrelated to thiamin pyrophosphate biosynthesis in general, or HET-P biosynthesis in particular. Also, the sequence alignments are quite uniform (matching the consensus), and so we do not believe that the RNAs have specialized based on any particular gene association. To clarify these consistent characteristics, we have edited the following text to the paragraph noted by the reviewers:
“Most other genes associated with thiS motif RNAs appear to code for proteins that participate in the production of HET-P or its fusion to 4-amino-5-hydroxymethyl-2-methylpyrimidine diphosphate (HMP-PP), to ultimately produce the bioactive coenzyme thiamin pyrophosphate (TPP) (Jurgenson et al., 2009).”
Point 1C: The reviewers ask for specific clarification on the thiS motif association with a putative thiazole transporter.
We have found that transporter annotations can be uncertain (and are frequently incorrect) and so we must speculate with considerable caution. However, the transporter genes associated with thiS motif RNAs would make most sense if they import molecules related to HET-P production. Specifically, we believe that the thiW gene association is a logical association for an HMP-PP “ON” switch because excess HMP-PP would trigger the import of the thiazole derivative HET, which can be used to biosynthesize more HET-P. In other words, we see no reason to doubt the current annotation of this transporter gene, and thus we have no evidence of a ligand specificity switch based on the current transporter gene annotations associated with thiS motif representatives.
Point 1D: The reviewers encourage us to discuss each of the genes associated with thiS motif RNAs and to describe the arrangements of their associated riboswitches.
Generally, the thiS gene is the gene most proximal to the RNA motif, whereas the other genes appear to be part of an operon regulated by HMP-PP riboswitches. Based on these gene names, and these operon arrangements, we believe the gene products catalyze the reaction steps for HET-P or TPP biosynthesis as depicted in Figure 1B. This list of gene functions has now been included in the new Supplementary file 1 that reports the various thiS representatives and their associated genes. Unfortunately we do not have a basis for discussing these predicted functions beyond that implied by the figure.
Point 1E: The reviewers ask if there is any special bias in the distribution of tandem TPP and HMP-PP (thiS) riboswitch systems.
We do not observe any special bias in the distribution of single versus tandem riboswitches, and therefore we cannot make any profound conclusions regarding the use of these arrangements. Presumably some species are exploiting the advantage of using tandem TPP and HMP-PP riboswitches, whereas others fail to take advantage of this arrangement. Such imperfect exploitations of tandem riboswitch arrangements have been observed for other riboswitch classes in the past.
Point 1F: The reviewers ask why the consensus models between the two structural states depicted in Figure 1A are different.
The vast majority of the nucleotide annotations in the original Figure 1A were identical between the two models. However, sequence alignments that match the proposed distinct structural states (specifically the complete terminator stem versus P1) caused slight differences in the consensus sequence model near the 5’ terminus (due to alignment alterations). To avoid depicting these differences, we have revised Figure 1A so that we depict the two structural states using only the consensus nucleotide sequence model for the aptamer configuration alignment. This updated figure is included in the revised manuscript file.
Point 1G: The reviewers state that our claim that the motif has the classic characteristics of an intrinsic terminator stem is incorrect.
We disagree. As we indicate in the original manuscript, the structural model as drawn on the left side of Figure 1A (terminator configuration) has all the key characteristics expected for a classic intrinsic terminator. Specifically, this part of the figure depicts a long hairpin with strong and continuous base-pairing, followed by a run of U nucleotides. Even as originally drafted, we do not believe that readers will be unduly confused by the right side of Figure 1A (aptamer configuration), as we carefully label these two structural states. Thus, given the clarity of the text and the figure annotations, we feel that the original text and the revised graphics should serve the readers well, and so no additional edits beyond those noted in Point 1F have been made.
2) The novel computational pipeline the authors describe is from unpublished work. Though not the subject of this manuscript, it is hard to see how this new approach differs from approaches used by this lab in the past. What makes this new approach more sensitive to rare riboswitch elements? Does this sensitivity arise from looking at smaller groups of bacteria (5 species) instead of across numerous species? The novelty of the approach is touted throughout, but not well explained.
The reviewers ask us to clarify the differences between the current discovery pipeline used to identify thiS RNAs compared to our previously used methods. Specifically, how is the current pipeline able to uncover rare riboswitch candidates?
These questions on the bioinformatics pipeline are fully addressed in the preceding manuscript (Stav et al.) that is mentioned in the response to Major Point 1A above. Briefly, our approach is a substantially updated version of a search strategy we first published in 2009. Specifically, this approach enriches for stretches of bacterial genomes that are likely to carry novel noncoding RNA motifs like riboswitches by collecting long, GC-rich intergenic regions (IGRs). This massively reduces the number of false positive candidates, and allows us to focus on motifs that are more likely to have the biochemical functions we seek. Although the delay in the initial review of the Stav et al. manuscript has been unusually long, as noted in the response to Point 1A, we have recently received word that our revised manuscript has been accepted for publication.https://doi.org/10.7554/eLife.45210.020
- Ronald R Breaker
- Ruben M Atilho
- Ronald R Breaker
The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.
We thank Adam Roth, Narasimhan Sudarsan, Shira Stav, and other members of the Breaker laboratory for helpful discussions. We are grateful to Rob Bjornson for assisting our use of the Yale Life Sciences High Performance Computing Center (NIH grant RR19895-02), Dr. Yong Xiong for assistance in the production and purification of the ThiD protein, and WeiWei Wang for assistance with mass spectrometry. RMA was supported by the National Science Foundation Graduate Research Fellowship Program (DGE1122492). This work was also supported by NIH grants to RRB (GM022778). In addition, RRB is supported by the Howard Hughes Medical Institute.
- Gisela Storz, National Institute of Child Health and Human Development, United States
- Ailong Ke, Cornell University
- Jinwei Zhang, NIDDK, United States
© 2019, Atilho et al.
This article is distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use and redistribution provided that the original author and source are credited.