|Home | About | Journals | Submit | Contact Us | Français|
OLE (Ornate, Large, Extremophilic) RNAs represent a recently discovered non-coding RNA class found in extremophilic anaerobic bacteria, including certain human pathogens. OLE RNAs exhibit several unusual characteristics that indicate a potentially novel function, including exceptionally high expression and localization to cell membranes via interaction with a protein partner called OLE-associated protein (OAP). In the current study, new genetic and phenotypic characteristics of OLE RNA from Bacillus halodurans C-125 were established. OLE RNA is transcribed at high levels from its own promoter under normal growth conditions and the transcript is exceptionally stable compared to most other RNAs. Expression is increased by ~7-fold when cells are exposed to near lethal concentrations of short-chain alcohols such as ethanol or methanol. Strains wherein the genes for OLE and/or OAP are deleted are more susceptible to growth inhibition by alcohol and also become more sensitive to cold. Normal growth characteristics can be restored by expressing the genes for OLE and OAP from plasmids or from elsewhere on the chromosome. Our findings confirm a functional link between OLE and OAP and reveal the importance of a large non-coding RNA in the response to alcohol-induced stress.
In addition to the fundamental tasks many RNAs perform in biological information transfer, non-coding RNAs are involved in many diverse biochemical functions in bacteria including RNA processing (1–4), gene regulation (5–8), defense against phage (9) and protein localization (10) among others. The number of discoveries of additional non-coding RNA classes has been increasing as more microbial DNA is sequenced and as bioinformatics search algorithms for structured RNA identification continue to improve. Intriguingly, novel classes of bacterial non-coding RNAs that rank among the largest and most complex known have been reported recently (11–14), suggesting that numerous additional RNAs with distinct biochemical functions important for many species remain to be discovered.
The functions of most newfound classes of large non-coding RNAs remain to be validated. RNAs whose functions were established long ago typically were discovered by investigating the molecular basis for specific cellular processes. For example, processing of precursor tRNAs to yield mature tRNAs necessarily involves RNase activity. Studies revealed the existence of a large non-coding RNA (15) that was determined to function as an RNA-cleaving ribozyme called RNase P (16). In the rare instances, when the specific non-coding RNA was discovered first (e.g. tmRNA and 6S RNA), decades passed between their initial identification and experimental validation of their biological and biochemical functions (17,18). As bioinformatics searches yield new non-coding RNA classes that are not yet linked to biological processes, additional work will be required to gather clues about their biochemical functions. Such efforts are particularly challenging for RNAs that are not involved in processes that are universal in biology or that are essential under typical laboratory growth conditions.
One of the recently discovered classes of large non-coding RNAs is called OLE (Ornate, Large, Extremophilic) RNA (11,19). OLE RNAs are similar in length and structural complexity to several large ribozymes, including Group I and II self-splicing introns and RNase P (13). OLE RNA representatives occur exclusively in extremophilic and/or anaerobic Firmicutes (Gram-positive bacteria) (11), including the causative agents for both botulism (Clostridium botulinum) and tetanus (Clostridium tetani). All organisms that contain OLE RNA also carry a gene for a predicted integral membrane protein of unknown function, called OLE-associated protein (OAP). We recently demonstrated that OLE and OAP form an RNA–protein complex in vitro and in vivo, and this binding is necessary for OLE RNA to localize to the membrane of a cell (19). Membrane localization by OLE RNA is almost unique among known non-coding RNAs. One prominent exception is the signal recognition particle (SRP) RNA from bacteria (10), which binds a soluble protein to form a complex that transiently localizes to the membrane (20,21).
Under normal growth conditions, OLE RNA is one of the most abundant transcripts in Bacillus halodurans C-125, a bacterial alkaliphile and facultative anaerobe (19). The number of copies of OLE RNA is nearly equal to that observed for RNase P or for SRP RNA. Due to its large size, complex structure, unusual phylogenetic distribution, unusual subcellular localization and high expression level, we hypothesized that OLE RNA might perform a biochemical function that has not yet been observed for non-coding RNAs. In an attempt to link OLE RNA to a biological process, we performed a series of bioinformatics analyses and genetic and molecular biology experiments using B. halodurans C-125, which is one of the few genetically tractable organisms that carries OLE RNA and that can grow in the presence of oxygen.
Bacillus halodurans C-125 was purchased from the ATCC (Catalog #BAA-125). Plasmids pMK4, pHCMC05, pE194 and pBGSC6 were obtained from the Bacillus Genetic Stock Center (The Ohio State University). Plasmid pUC19 was purchased from Invitrogen, and pGFPuv was purchased from Clontech. Unless otherwise specified, B. halodurans was grown in LB (Lysogeny Broth) broth (USB Corporation) that was prepared at 90% volume, autoclaved and adjusted to full volume and pH ~10.5 with 10% (w/v) filter-sterilized Na2CO3 [1% (w/v) final concentration]. Anaerobic cultures were grown in nutrient broth [NB; 5g peptone (Difco) and 3g beef extract (Difco) per liter] supplemented with 1% (w/v) glucose (J.T. Baker), 5mM KNO3 (Sigma-Aldrich) and 1μg ml−1 resazurin (Sigma-Aldrich), and adjusted to pH ~10.5 as described above. Media was solidified with 1.5% (w/v) agar (BD) as needed.Unless otherwise indicated, all cells were grown at 37°C and shaken at ~200rpm.
Mutual information (MI; a measure of the dependence of one variable on another, such as covariance of two nucleotides) analysis and other comparative sequence analyses (22,23) were conducted on OLE RNAs from 52 sequenced genomes and 156 unique sequences from environmental DNA samples. Sequences were aligned using CMfinder (http://bio.cs.washington.edu/yzizhen/CMfinder/ 20 April 2012, date last accessed) (24) and MI content was calculated as described previously (5). P-values were calculated by comparing MI scores with those from 100000 neutral-evolution simulations to determine the probability of a given MI score arising by chance. The OLE RNA secondary structure diagram was generated by first processing the alignment with the RNA structure drawing algorithm R2R (http://breaker.research.yale.edu/R2R/ 20 April 2012, date last accessed) (25) and then manually adjusting the output to improve clarity.
Total cellular RNA was extracted from cells that were pelleted from media and resuspended in TE buffer (10mM Tris–HCl, 1mM EDTA, pH 8.0 at 23°C) containing 3mg ml−1 lysozyme. The resuspended mixture was incubated either at 23°C or on ice for 5min. Cells were lysed by freeze–thawing three times, and total RNA was extracted with TRIzol reagent (Invitrogen) according to the manufacturer’s directions. Residual genomic DNA was removed with RQ1 DNase (Promega) according to the manufacturer’s directions.
To assess possible splicing of OLE RNA, total RNA (1μg) was reverse-transcribed with Superscript II reverse transcriptase (Invitrogen) according to the manufacturer’s protocol with 50ng random hexamer. cDNA was then amplified in 50μl PCR reactions containing 200nM primers (each), 200μM dNTPs (NEB), 2.5 U Taq polymerase (NEB), PCR buffer [1.5mM MgCl2, 50mM KCl, 10mM Tris–HCl (pH 8.3 at 23°C), 0.01% (w/v) gelatin] and 0.5μl cDNA. Amplification was performed for 30–50 cycles with an extension at 72°C for 2min, and products were separated and analyzed by agarose (0.8%) gel electrophoresis.
Plasmids for knocking out OLE and OAP were constructed using techniques similar to those reported previously (26). Briefly, the TaqI fragment of pE194 containing its temperature-sensitive oriC was amplified via PCR with primers that introduce BamHI restriction sites. The PCR product was cloned using a TOPO-TA cloning kit (Invitrogen), after which the insert was removed by digestion with BamHI (NEB), gel purified and ligated into pUC19 to generate plasmid pUCE194. Inserts to knockout-specific genes (as described in Figure 4a) were created by overlap-extension PCR and ligated into the EcoRI site of pUCE194. Plasmids were then transformed into B. halodurans C-125 via protoplast transformation as described previously (27) and the cells were recovered on medium containing 3μgml−1 chloramphenicol. Subsequent transformations were conducted using improvements in B. halodurans cloning as described elsewhere (28). All constructs were confirmed via sequencing by the W. M. Keck Biotechnology Resource Laboratory at Yale University.
Plasmid integration was achieved using published protocols (26,29) with minor modifications. In brief, overnight cultures grown at 30°C were diluted 1:100 in fresh media with antibiotic (chloramphenicol, 3μgml−1) and grown 2h at 30°C. Cultures were then shifted to 45°C and grown 3h to dilute plasmid concentration per cell. Serial 10-fold dilutions of each culture were then plated on LB pH ~10.5 agar with 3μgml−1 chloramphenicol and incubated at 45°C overnight to recover colonies that had integrated the plasmid. Double-crossover events were identified via PCR. The chloramphenicol resistance cassette was subsequently excised through Xer recombination (30) as follows. Cells were grown 24h at 37°C in LB (pH~10.5; no antibiotic) and plated on LB (pH~10.5) agar for overnight growth. The resulting colonies were replica-plated onto LB (pH~10.5) agar with 3μgml−1 chloramphenicol and incubated at 37°C overnight. Colonies that failed to grow on the chloramphenicol-containing plates were picked from the initial plate and screened for the desired genetic change by PCR and DNA sequencing.
Competitive co-culture experiments were conducted by inoculating liquid LB medium (pH~10.5) with equal amounts of exponentially growing wild-type (WT) and knock-out (KO) cells (0.001 initial OD600) and incubating them with shaking for 0, 24 or 48h in LB (pH~10.5) containing 6% (v/v) ethanol. Culture samples were serially diluted 10-fold and plated on LB (pH~10.5) agar, incubated overnight at 37°C, and individual colonies were picked and subjected to PCR using primers specific for the ole-oap locus. PCR products corresponding to amplification of the full-length WT locus or the shortened KO locus were identified by separation using 1% agarose gel electrophoresis and staining with ethidium bromide. Cell growth was estimated by counting the total number of colonies per plate and using the plate’s dilution factor and the PCR genotype proportions to determine the total number of cells of each genotype at each time point.
Transcriptome analyses were conducted in duplicate using wild-type and ΔoleΔoap B. halodurans C-125 strains grown to mid-exponential phase (OD600 0.2–0.3) in LB (pH~10.5). Triplicate subcultures of 3ml each were incubated in sealed tubes for 3h at 37°C with shaking after the addition of ethanol (95% aqueous-distilled, Pharmco-AAPER) to 5% (v/v) concentration. In all instances, ethanol concentrations noted are initial values and may change due to evaporation or bacterial metabolic activities. Cultures were pooled and RNA was extracted as described above, and rRNA was removed by using a RiboMinus Kit for Bacteria (Invitrogen). cDNA library preparation and sequencing was conducted by Otogenetics (Tucker, Georgia, USA). Reads were aligned to the B. halodurans reference genome with SOAP2 (http://soap.genomics.org.cn/ 20 April 2012, date last accessed) (31), and the expression of individual genes was tabulated by custom Perl scripts. Significant differences in gene expression were determined with edgeR (http://www.bioconductor.org/packages/release/bioc/html/edgeR.html 20 April 2012, date last accessed) (32), and relative expression levels were computed by calculating RPKM scores (Reads Per Kilobase per Million reads) (33). Raw and mapped data are available at the NCBI Gene Expression Omnibus database (http://www.ncbi.nlm.nih.gov/geo/ 20 April 2012, date last accessed), accession GSE33303.
Constructs for plasmid-based genetic rescues were each constructed using plasmid pHCMC05 carrying the appropriate ole-oap construct controlled by an IPTG-inducible promoter. Plasmids were transformed into the indicated B. halodurans C-125 strain as described earlier. Rescue transformants were grown in media containing 3μgml−1 chloramphenicol and induced with 1mM IPTG for at least 3h prior to the start of the assays.
Chromosomal rescue constructs were made by inserting the oriC-containing TaqI fragment of pE194 into pUC19 at the HindIII restriction site to create plasmid pUCE195. 5′ and 3′ homology segments of the lacZ gene (BH2723) of B. halodurans and a chloramphenicol resistance cassette from pBGSC6 were amplified by PCR with primers containing the necessary restriction sites for cloning into pUCE195. Integration of the final vector disrupts lacZ and places the chloramphenicol resistance cassette on the chromosome. For rescue, the full ole-oap region, including the promoter (see text), was amplified and inserted upstream of the chloramphenicol resistance cassette. Plasmids were transformed into B. halodurans C-125 and integrated as described for making KOs.
Stress assays were conducted using 3ml cultures (LB at pH~10.5) of B. halodurans C-125 cells in exponential growth phase. For ethanol stress, the initial OD600 was 0.001 for 48-h assays, or ~0.1 for 24-h assays. Ethanol was added to culture tubes to a final concentration of 5% (northern blots) or 6% (growth assays). Tube caps were sealed and the cultures shaken at 37°C for the indicated times. For cold stress, cells were diluted to OD600 0.01 and placed in a 15°C incubator for 7 days.
GFPuv mRNA reporter constructs were prepared by PCR amplification from pGFPuv (Clontech) and fused to the OLE RNA promoter via overlap-extension PCR. The resulting construct was cloned into plasmid pMK4 and transformed into B. halodurans as described above. mRNA levels were assayed via northern blot analysis.
RNA decay rates were established by growing wild-type B. halodurans to mid-exponential phase (OD600~0.3) in LB (pH~10.5). Rifampicin (Sigma-Aldrich) was added to 100μgml−1 final concentration, and aliquots of the culture were withdrawn at specific time points. Total RNA was extracted as described above and probed via northern blot analysis.
qRT-PCR was conducted on total RNA samples extracted as described above. Equal amounts of total RNA (0.5–1.0μg) were reverse-transcribed with Superscript II reverse transcriptase (Invitrogen) according to the manufacturer’s protocol with 100ng of random hexamers. The resulting cDNA was diluted by the addition of nine volumes dH2O and subjected to quantitative PCR on an Applied Biosystems 7500 Real-Time PCR system using reactions consisting of 10μl Power SYBR Green Master Mix (Applied Biosystems), 4μl water, 1μl primer mix (10mM each) and 5μl template. Expression levels were determined according to published protocols (34).
Northern blot analyses were performed according to published protocols (35). Hybridization was carried out in GE Rapid-Hyb buffer at 44°C according to the manufacturer’s protocol. Signals were normalized to 16S rRNA for quantification unless otherwise specified.
A MI analysis (22,23) was conducted using approximately 200 unique OLE RNA representatives derived from sequenced genomes and from metagenomic data (19). MI is a statistical measure of how predictable the identity of one particular nucleotide is given the identity of another, and this method has been used previously to predict and refine the structures of other conserved RNAs (e.g. 5). MI analysis is particularly well suited to reveal base-paired structures that can co-vary, although other tertiary-structure contacts also can be revealed.
MI scores confirm most of the secondary-structure elements predicted previously (Figure 1). In addition, this analysis revealed the likely existence of three previously unidentified stems: P4a, P14a and P15. P4a is formed in part from one of two short regions of conserved RNA sequence proposed to be important for OLE RNA binding to the protein partner OAP (19). Formation of this newfound P4a stem would make the local structure at this first protein-binding site similar to the structure of the second proposed OAP-binding site (stem P2). It is not known whether protein binding at each site requires stem formation as predicted by our MI analysis, but such a requirement could explain the selective binding of OLE RNA by OAP despite the relatively low information content of the binding site (19).
Formation of stems P14a and P15 further supports our original suggestion (11) that the nucleotides in this region form a sequence and structural domain with partial 2-fold symmetry. Three of four loop sequences conform to a consensus GNRA tetraloop, and analogous hairpins are known to stabilize RNA structures and can dock to the minor groove of distal helices of specific architectures (36–38). The remaining sequence is typically GAAAN, which closely approximates a GNRA tetraloop. The proximity of these loops to one another could indicate that several distal portions of OLE RNA may form tertiary contacts with this symmetrical domain. This entire region of OLE RNA (P12 through P15) could function as a stabilizing scaffold analogous to Domain I in Group II introns (39). However, there are numerous other feasible reasons for the clustering of these tetraloops, including the possibility that they might form contacts with another RNA molecule.
Our MI analysis did not reveal the presence of potential pseudoknots, which is surprising for an RNA of this size and structural complexity. Perhaps OLE RNAs do not rely on pseudoknots to form their tertiary structure, or perhaps pseudoknots are present but their component nucleotides are subject to other selective pressures that keep their identities too highly conserved for MI analysis to reveal co-variation. However, several distal interactions that indicate non-canonical tertiary contacts (Figure 1 and Inset) were revealed. Additional tertiary contacts are also likely to exist but there were insufficient occurrences of co-variation to declare their existence with confidence. The newfound structural elements further refine our model for OLE RNA and place it firmly among the most complex bacterial RNAs known. Only the two large ribosomal RNAs, Group II self-splicing introns, and the recently identified GOLLD RNA are known to conserve greater lengths and structural complexity (13).
In B. halodurans, OLE RNA is a highly abundant transcript whose expression level is within a few folds of that for the RNase P and SRP RNAs (19). At least some OLE RNA is co-transcribed with neighboring genes (11), but a putative promoter located immediately upstream of the RNA gene also has been identified (40) (Figure 2a). To determine if this proximal promoter is used, we cloned the OLE RNA gene (hereafter called ole) either with or without the putative promoter into a plasmid and transformed the constructs into a B. halodurans strain with an ole deletion (Δole, see description below).
Total RNA was extracted from cells isolated either during exponential growth or stationary phases. The RNA was subjected to quantitative reverse transcription and polymerase chain reaction amplification (qRT-PCR), and the levels of OLE RNA transcription were determined and normalized relative to the levels of mRNA expression of the chloramphenicol acetyltransferase (cat) antibiotic resistance gene also carried by the plasmids. Inclusion of the promoter increases OLE expression by >100-fold (Figure 2b), demonstrating that the promoter is indeed active under normal growth conditions.
To assess the stability of OLE RNA, we treated exponentially growing B. halodurans C-125 with rifampicin to halt RNA synthesis and extracted total RNA at various times. A northern blot analysis was conducted to detect several abundant transcripts, including OLE (Figure 2c). In contrast to 16S rRNA and the mRNA for the tuf gene (encoding elongation factor Tu, which is expressed at approximately the same level as OLE RNA) (19), OLE RNA exhibits almost no degradation even 4h after the addition of rifampicin. However, each lane of the gel was loaded with an equivalent amount of RNA, and SRP RNA (an exceptionally stable RNA) appears to be increasing in abundance because it is becoming a greater proportion of the total RNA remaining. Therefore, although OLE RNA appears to maintain a constant level on the gel, it is likely to be degrading slowly over the course of this assay.
When normalized to the amount of SRP RNA, the half-life of OLE RNA is estimated to be ~3h. Given that bacterial doubling time under these conditions is ~45min (data not shown), it appears that the decay of OLE RNA does not play an appreciable role in cells under these growth conditions. We also evaluated the stability of OLE RNA in B. halodurans cells grown at neutral pH or in cells lacking the OAP protein, but no variation in stability was observed (data not shown). Our findings indicate that OLE RNA could be considered as a ‘stable’ RNA (41), and might be degraded only in times of stress or damage.
Northern blot analyses reveal that there is only a single product band for OLE RNA (Supplementary Figure S1). Furthermore, previous 5′ RACE data (11) revealed that sites of OLE RNA processing are located only near the start of the P1 stem (Figure 2a). These sites are very close to the predicted transcription start site for the proximal promoter, indicating that the vast majority of independent OLE RNA transcripts and any OLE RNAs processed from larger RNA precursors likely have similar 5′ termini.
Self-splicing ribozymes are similar in size and complexity to OLE RNAs (13), which raises the possibility that OLE RNAs could represent a new type of self-splicing RNA. However, the northern blot data described above suggest that the RNA does not undergo splicing because there is no evidence of smaller processed products. To further evaluate this hypothesis, we performed RT-PCR on RNA extracts from B. halodurans growing under various conditions. No spliced transcripts were ever detected (Figure 2d and e), even after 50 cycles of PCR amplification (data not shown). Additionally, although we have not rigorously ruled out an RNA self-cleavage function, the size and complexity of OLE RNAs relative to known self-cleaving ribozymes (42) implies that if self-cleavage activity exists, this is unlikely to be its sole function.
Cell culture conditions were sought under which OLE RNA expression is altered relative to its established baseline. Due to the strong association of OLE and OAP genes with extremophilic organisms (11), we hypothesized that the RNA might play a role in protecting cells against stresses present in extreme environments. Bacillis. halodurans C-125 cells were subjected to a variety of stresses (Supplementary Table S1) and OLE expression changes were monitored via northern blot analyses. Conditions that produced a 2-fold or larger change in OLE RNA levels were reassessed with a reporter construct (gfpuv mRNA fused to the promoter for OLE RNA) (Figure 2a) to confirm changes in transcript levels. The only conditions that consistently altered OLE RNA expression were exposure to short-chain alcohols (Figure 3a). Treatment with heat (55°C, the limit of B. halodurans growth) or shift to anaerobic conditions also appeared to induce OLE (Supplementary Table S1 and Supplementary Figure S2), but these results were not consistently reproducible with the reporter construct.
All three short-chain alcohols examined (methanol, ethanol and isopropanol) induce OLE RNA production. Likewise, these alcohols induce expression of the gfpuv reporter gene. Each alcohol was tested at a concentration that is near equivalent in toxicity compared with the other alcohols, and the expression levels are similar in all three analyses. This result suggests that the increase in OLE expression is not due to any specific molecular recognition event, but rather to a common effect caused by these compounds.
Given the high-basal expression of OLE RNA (19), its biochemical function is probably also beneficial under normal growth conditions but becomes more beneficial when cells are stressed with alcohols. Alcohols cause many stresses on cells, but one of the strongest effects is to make cell membranes more permeable to ions and small molecules (43,44). Attempts to duplicate this membrane permeability stress with detergents or ionophoric antibiotics failed (Supplementary Table S1), which suggests that there may be a different aspect of alcohol stress that is responsible for triggering increased OLE production.
Since ethanol is presumably the short-chain alcohol that is most common in the natural environment, we performed a more detailed analysis of the response of B. halodurans to sublethal ethanol concentration (Figure 3b). Addition of ethanol causes an almost immediate increase in both OLE and reporter transcripts. Ethanol-mediated induction also occurs at neutral pH (~7.0) (Supplementary Figure S2), indicating that this response is independent of the pH of the medium.
A series of KO strains were generated in B. halodurans C-125 cells to assess the effects of the loss of OLE and OAP. These KO strains carry internal deletions of ole, oap or both, and were constructed in a manner that minimizes the possibility of unanticipated polar effects on neighboring genes (Figure 4a). The KO strains were subjected to a variety of phenotype assays (Supplementary Table S2), wherein the vast majority of conditions generally show no difference from wild type. Only under conditions of alcohol stress (Figure 4b) or cold stress (Supplementary Figure S3) are KO cells significantly affected. Notably, the growth defect under alcohol stress corresponds well with the observation that OLE RNA production is increased under these same conditions. This finding supports the hypothesis that OLE RNA is important for resisting alcohol-induced stress. Cold stress causes growth reduction for the KO cells, but does not induce OLE RNA expression (Supplementary Table S1). Even though OLE RNAs are present in many different types of extremophiles, they have not been found in psychrophiles (cold-loving organisms). Therefore, the relevance of OLE RNA to cold stress adaptation is less clear.
Ethanol causes the largest growth defect among the short-chain alcohols tested (Figure 4b; Supplementary Figure S4). For each alcohol type, all three KO strains (Δole, Δoap and ΔoleΔoap) are affected similarly, suggesting that OLE and OAP are part of the same pathway and affect the same biological function. To more precisely quantify growth differences, we established competitive co-cultures with WT and ΔoleΔoap KO cells in the presence of 6% ethanol and removed samples at 24 and 48h to establish the relative fraction of each cell type by plating, picking colonies and classifying their type by PCR (Figure 4c). The fraction of KO cells relative to WT cells decreases by >3.6-fold (Figure 4d), which is consistent with our observations of reduced growth of KO cells compared with WT when they are cultured independently (Figure 4b). We also performed quantitative PCR of co-culture genomic DNA to ensure that each colony reflected cell growth and not other effects (e.g. cell aggregation), with the same results (data not shown).
To determine if the absence of OLE or OAP is responsible for the observed growth defects caused by exposure to ethanol or to cold, we created a series of B. halodurans rescue strains to supply the missing gene(s) by either plasmid or by chromosomal integration at a site distal from the original locus. Expression of plasmid-borne ole and oap, controlled by an IPTG-inducible promoter, resulted in partial rescue of growth when cells were grown in 6% ethanol (Supplementary Figure S3a). Likewise, the rescue strain carrying ole and oap on the chromosome at a distal locus (lacZ, BH2723) under control of the native promoter partially restores resistance to ethanol (Figure 4b) and yields complete restoration of growth under cold stress conditions (Supplementary Figure S3b). These results indicate that the observed phenotypes are due to OLE and OAP rather than to unanticipated effects on neighboring genes.
In an attempt to link other pathways to the ethanol sensitivity phenotype, RNA-seq (45) was performed on both WT and ΔoleΔoap cells before and after growth for 3h with 5% ethanol. This concentration was chosen because it strongly induces OLE RNA production in WT cells (Figure 3b) but KO cells experience only minor growth defects (data not shown), implying they are largely able to compensate for lack of OLE RNA at this concentration. Transcript analyses reveal essentially no difference between WT and KO cells before ethanol stress (Figure 5a and data not shown), consistent with the lack of a phenotype under these conditions. After ethanol stress, the expression levels of only a small number of genes are significantly different between the strains (Figure 5b and Supplementary Table S3). However, we did not identify any obvious connections between their functions that would guide the design of additional experiments to directly reveal the biochemical function of OLE RNA. Furthermore, although OLE RNA is observed to increase substantially upon ethanol exposure, the mRNA for OAP increases only modestly (Supplementary Table S4). This suggests that if OAP levels increase, regulation may occur at the translation level and not at the level of mRNA abundance.
Since the genes with significantly different expression levels between strains represent only 1.5% of the approximately 1600 genes that change expression during the experiment (data not shown), the absence of OLE and OAP appears to have only a minor effect on the transcriptional response of cells to ethanol stress. It is curious that the deletion of such a highly expressed RNA would have so little effect on the transcriptome, especially since under ethanol stress OLE RNA is the fifth most expressed transcript (excluding rRNAs and tRNAs) (Figure 5b). Despite the fact that the point of maximum OLE RNA induction occurs at around 3h post-exposure to ethanol (Figure 3b), the cell’s transcriptional response to the absence of OLE RNA may occur at a different time. Alternatively, compensation may occur mostly at the post-transcriptional level, or perhaps the functions of OLE and OAP are at the end of a biochemical pathway that has little effect on other pathways.
Both the induction of OLE RNA in response to short-chain alcohols and the alcohol-sensitive phenotype of ole and oap deletion strains demonstrate that at least one function of OLE and OAP is to protect cells from alcohol-induced stress. Alcohols cause a wide variety of stresses on cells, but one major effect is permeabilization of membranes to ions and small molecules (43,44). In response to ethanol-mediated membrane stress, bacteria make changes in the protein and lipid composition of membranes (46). Since OLE and OAP localize to the cell membrane (19), this complex also may be important for fortifying the membrane against leakage. The large size and complexity of OLE RNA, which is similar to that of RNase P and self-splicing introns, implies that the OLE–OAP complex might function as a ribonucleoprotein particle with enzymatic activity. Other possibilities (e.g. gene regulation, signal transduction, etc.) are still possible, and therefore further experiments will be required to determine the biochemical function of the RNA and its protein partner.
It is unknown how accurately the conditions of our cell culture media approximate the natural environment of OLE-containing organisms. Little data exist on the presence of alcohols in the microenvironments of most organisms that carry OLE RNA (e.g. thermal vents, soda lakes and vertebrate gut), but, in many cases, it would be expected to be low. It is possible that local microenvironments could allow alcohols to build up to appreciable concentrations, or that natural stresses are simply less severe and cause a smaller but still significant fitness benefit for organisms with OLE. Even small fitness benefits can lead to gene fixation, and the extreme ethanol concentrations used in our assays may simply push the cells to the point where the fitness differences become large. Indeed, the 5–6% ethanol concentrations used in this study are at the very limit for viability of B. halodurans.
An interesting question is why cells would use a large and complex non-coding RNA to protect against alcohol stress instead of exclusively exploiting protein factors. Other widely conserved large non-coding RNAs in bacteria appear to be ancient RNAs that carry out fundamental biochemical tasks (e.g. RNAs in translation or RNAs processing RNA) that either are evolutionarily difficult to replace or are particularly well suited for RNA, such as specifically recognizing other nucleic acids. Given these precedents, it seems possible that the biochemical function of OLE may be ancient and fundamental, or involve other nucleic acids, or both.
Supplementary Data are available at NAR Online: Supplementary Tables 1–4, Supplementary Figures 1–5 and Supplementary Reference .
National Institutes of Health [GM 022778]; Howard Hughes Medical Institute. Funding for open access charge: Howard Hughes Medical Institute.
Conflict of interest statement. None declared.
We are grateful to members of the Breaker laboratory and to Dr Nick Ornston for helpful discussions.