Skip to main content

Main menu

  • Online first
    • Online first
  • Current issue
    • Current issue
  • Archive
    • Archive
  • Submit a paper
    • Online submission site
    • Information for authors
  • About the journal
    • About the journal
    • Editorial board
    • Information for authors
    • FAQs
    • Thank you to our reviewers
      • Thank you to our reviewers
    • American Federation for Medical Research
  • Help
    • Contact us
    • Feedback form
    • Reprints
    • Permissions
    • Advertising
  • BMJ Journals

User menu

  • Login

Search

  • Advanced search
  • BMJ Journals
  • Login
  • Facebook
  • Twitter
JIM

Advanced Search

  • Online first
    • Online first
  • Current issue
    • Current issue
  • Archive
    • Archive
  • Submit a paper
    • Online submission site
    • Information for authors
  • About the journal
    • About the journal
    • Editorial board
    • Information for authors
    • FAQs
    • Thank you to our reviewers
    • American Federation for Medical Research
  • Help
    • Contact us
    • Feedback form
    • Reprints
    • Permissions
    • Advertising
Open Access

Bioinformatic and experimental methods to identify and validate bacterial RNA-human RNA interactions

Kylie I Krohmaly, Robert J Freishtat, Andrea L Hahn
DOI: 10.1136/jim-2022-002509 Published 26 September 2022
Kylie I Krohmaly
1Center for Genetic Medicine Research, Children's National Research Institute, Washington, District of Columbia, USA
2Institute for Biomedical Sciences, The George Washington University School of Medicine and Health Sciences, Washington, District of Columbia, USA
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Robert J Freishtat
1Center for Genetic Medicine Research, Children's National Research Institute, Washington, District of Columbia, USA
3Division of Emergency Medicine, Children's National Hospital, Washington, District of Columbia, USA
4Department of Pediatrics, The George Washington University School of Medicine and Health Sciences, Washington, District of Columbia, USA
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Andrea L Hahn
1Center for Genetic Medicine Research, Children's National Research Institute, Washington, District of Columbia, USA
4Department of Pediatrics, The George Washington University School of Medicine and Health Sciences, Washington, District of Columbia, USA
5Division of Infectious Diseases, Children's National Hospital, Washington, District of Columbia, USA
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
  • ORCID record for Andrea L Hahn
  • Article
  • Figures & Data
  • eLetters
  • Info & Metrics
  • PDF
Loading

Abstract

Ample evidence supports the importance of the microbiota on human health and disease. Recent studies suggest that extracellular vesicles are an important means of bacterial-host communication, in part via the transport of small RNAs (sRNAs). Bacterial sRNAs have been shown to co-precipitate with human and mouse RNA-induced silencing complex, hinting that some may regulate gene expression as eukaryotic microRNAs do. Bioinformatic tools, including those that can incorporate an sRNA’s secondary structure, can be used to predict interactions between bacterial sRNAs and human messenger RNAs (mRNAs). Validation of these potential interactions using reproducible experimental methods is essential to move the field forward. This review will cover the evidence of interspecies communication via sRNAs, bioinformatic tools currently available to identify potential bacterial sRNA-host (specifically, human) mRNA interactions, and experimental methods to identify and validate those interactions.

Section 1: background and evidence of bacteria-to-human communication via small RNAs

The function and dysfunction of microbiota have been implicated in numerous human conditions, including inflammatory bowel diseases, asthma, and even infertility.1–3 However, the mechanism by which the community of microbiota that reside on and in humans contribute to human disease remains poorly understood. More recent evidence supports that microbiota impact the host in part via molecular mechanisms induced by bacterial extracellular vesicles (EVs) and their contents.4

Bacterial EVs5 are lipid bilayer-bound products deriving from the outer and/or inner membranes of bacteria6 and can be between 10 and 500 nm depending on the species.7 EVs can be isolated through various methods including filtration, ultracentrifugation, precipitation kits, or size-exclusion chromatography, and characterized by size, protein content, lipid content, function, and more.8 There is no gold standard for bacterial EV isolation and characterization, however the minimal information for studies of EV guidelines updated in 2018 provide recommended protocols and reporting criteria to enhance reproducibility of methods.8

Bacterial small RNAs are among the RNA cargo of bacterial EVs

Bacterial EVs carry a range of cargo, including RNA, DNA, and proteins.6 9 10 RNA sizes and types found within EVs include protein coding genes, transfer RNAs (tRNAs), ribosomal RNAs (rRNAs), and short non-coding RNAs.9 11–13 Short non-coding RNAs are regulatory molecules that are able to modulate gene expression via base pairing with a potential target transcript.14 Eukaryotic short non-coding RNAs often take the form of microRNAs (miRNA) which are processed from miRNA genes,15 while in bacteria, short non-coding RNAs are called small RNAs (sRNAs) and are generally derived from intergenic regions.16,17 tRNA fragments are considered a type of sRNA,18,19 and are found in bacterial EVs.9 Additionally, both rRNA and tRNA fragments have been implicated in inducing apoptosis in human cells.20

Whereas eukaryotic miRNAs are around 20 nucleotides (nts) long, bacterial sRNAs can be as long as 500+ nt,21 but are generally between 50 and 400 nt. This extra length allows bacterial sRNAs to form complex secondary structures that impact the stability of the sRNA, how it regulates its target transcript, and how and if it binds its target.22–2829 30 Bacteria also produce miRNA-sized RNAs,18 29–33 which may function similar to human miRNAs. sRNA is a consistent and often large portion of EV cargo.9 11–13 Interestingly, bacterial EV RNA cargo is distinct from the intracellular RNA population,9 18 34 and from the non-EV extracellular RNA population,9 indicating that RNA is strategically loaded into EVs for their target cell.

Bacterial sRNAs circulate in the bloodstream and affect mammalian gene expression

Bacterial sRNAs35–37 and particularly those encapsulated in bacterial EVs38 39 have been detected in the plasma of both healthy individuals and those that are affected by various illnesses. Mammalian cells internalize bacterial EVs11 40 and furthermore, EV cargo—specifically, their RNA cargo—localizes to both the cytoplasm33 40 and the nucleus.11 33 40 Exogenous bacterial EVs also traverse the blood-brain barrier.33 40

Bacterial sRNAs have been bioinformatically predicted and experimentally validated to affect certain pathways in human and murine cells. Using three periodontal pathogens’ EVs, Choi et al confirmed the presence of miRNA-sized sRNAs.32 Using bioinformatic tools, they predicted human immune system-related mRNA targets of several sRNAs.32 Jurkat T cells exposed to synthetic oligos of the most highly expressed sRNAs showed that these sRNAs downregulated several cytokines (interleukin (IL)-5, IL-13, and IL-15).32

Additionally, Han et al investigated the potential of EV-derived RNAs from Aggregatibacter actinomycetemcomitans (Aa), a periodontal pathogen, to induce an immune response.29 U937 cells exposed to both intact Aa EVs and lysed Aa EVs (both exposed to RNase) showed that intact EVs induced a greater tumor necrosis factor-alpha response than lysed EVs.33 Compared with RNase-treated lysed EVs, cells exposed to DNase-treated lysed EVs also upregulated nuclear factor kappa B (NF-κB) activation and toll-like receptor 8 expression (both of which sense single-stranded RNAs); the authors attributed the difference in outcome to the RNAs of the DNase-treated lysed EVs.35 The authors also found that BV2 (ie, murine microglial) cells activated NF-κB and upregulated IL-6 after exposure to DNase-treated lysed Aa EVs or intact EVs more so than RNase-treated lysed EVs.41

Bacterial sRNAs interact with mammalian RISC

RNA-induced silencing complex (RISC) in mammalian cells comprises Argonaute (Ago) and a miRNA, which then can bind the miRNA’s target mRNA to result in its degradation.15 Han et al described bacterial sRNAs that interact with host RISC after Aa EV exposure to U397 cells using total RNA sequencing (RNA-Seq) and RNA immunoprecipitation and sequencing (RIP-Seq).33 While Aa-derived sRNAs only accounted for 0.02% of reads via RNA-Seq, 0.35% of reads originated from Aa when looking at RISC-associated sRNA through RIP-Seq.33

In a similar study, Furuse et al investigated five pathogenic intracellular bacteria which produced RISC-enriched RNAs (compared with total RNA levels) that were close to the same length as authentic eukaryotic miRNAs.29 Additionally, one bacterial species, Mycobacterium marinum (Mm), produced an RNA, MM-H, the expected size (23 nt) with similar characteristics to canonical miRNAs (ie, stem-loop structure prediction, discrete 5’ end consistent with Dicer trimming).29

These studies demonstrate that bacterial sRNAs have the potential to operate as human regulatory RNAs via RISC. Bacterial sRNAs—specifically, bacterial miRNA-sized RNAs—interacting with RISC are a likely mechanism by which some sRNAs exert their impact on host gene expression.29 33 While these findings taken together are very compelling, more studies are needed to confirm this potential mechanism and investigate any off-target effects (perhaps relating to the machinery itself), or the possibility that host miRNA may be targeting the bacterial sRNAs. Nonetheless, discovering that sRNAs interact with miRNA machinery is encouraging evidence that the interactions and resulting effects may be sequence-specific rather than due to foreign single-stranded RNA.

Although Furuse et al investigated if MM-H interacted with a perfectly complementary target using a reporter assay,29 no human target prediction was performed for the sRNA. Future studies will need to investigate the potential specific targets of the discovered sRNAs more deeply. There is a need for identification and validation to determine sRNA candidates for microbial-host communication in conditions for which host microbiota have been implicated (eg, inflammatory bowel disease or asthma).1 2 In the following two sections, we will discuss bioinformatic methods (table 1) and experimental methods (table 2) that could be employed to discover these interactions and identify and validate potential human targets of bacterial sRNAs.

View this table:
  • View inline
  • View popup
Table 1

Summary of bioinformatic tools that can be used to predict bacterial RNA-human RNA targets

View this table:
  • View inline
  • View popup
Table 2

Summary of experimental approaches that can be used to identify and validate bacterial RNA-human RNA interactions

Section 2: bioinformatic approaches for predicting bacterial sRNA-human mRNA interactions

Homology and alignment tools for sRNA target prediction in silico

As RNA-RNA interactions are mediated via base pairing, sequence complementarity for at least the seed region (ie, the perfectly complementary region that ‘seeds’ the interaction) is required. This region of the sRNA may be conserved across species, and a novel sRNA’s function in one species may be similar to that of a homologous or identical sRNA (ie, a perfect match after alignment) in another species.

Koeppen et al identified human targets for 10 Pseudomonas aeruginosa sRNAs in EVs via BLAST (and other tools discussed below), which they then used to guide their functional validation.18 Jorth and Whitely, after identifying and characterizing a lysine riboswitch controlling a lysine transporter in Aa, identified a homologous transporter gene in Haemophilus influenzae (Hi) through BLASTp.42 They confirmed the riboswitch’s effect on the transporter in Hi, demonstrating conservation of the system across the two species.42 In other words, they aligned one species’ sRNA target in BLAST in order to determine if the latter produced a similar target. This type of approach can be ued in identifying sRNA targets in different genomes; after identifying the transcript and homologous transcripts in other genomes, and one can then go backwards from there to identify sRNAs that target the transcript. BLAST (and other alignment tools combined with appropriate databases such as bowtie41) can also be useful in identifying sRNAs themselves through alignment to sRNA databases—with whole sequences, and with portions of sequences that may be conserved. For more tools on sRNA prediction and detection, refer to table 3.

View this table:
  • View inline
  • View popup
Table 3

Summary of bioinformatic tools that detect bacterial small RNAs and their genes

The use of alignment and homology to annotate sRNAs with targets in different species (and even within a species) is currently limited in part because there may be species-specific sRNAs, and because there are still many sRNAs that need functional validation and many species left to investigate. For example, one group assigned functions to detected bacterial sRNAs using Rfam43 and other databases by aligning the detected sRNAs to the annotated sRNA sequences within the databases.44 However, only 19 of 422 sRNAs were assigned a function.44 These databases can contain predicted and experimentally validated sRNAs. sRNA alignment and homology may be more widely used after more sRNA targets have been identified and fully characterized.

Importantly, alignment to directly determine sRNA targets (ie, a reverse-complement match to an sRNA) may be less reliable, however, as the secondary structure of the sRNA and mRNA are not taken into account in alignment tools. This is important to consider because although a sequence match may be found between a pair, the sequences may not be able to interact—for example, the seed region may be occluded.

Energy-based and accessibility-based models consider secondary structure

Besides homology, energy-based models are RNA-RNA interaction prediction algorithms and pipelines that calculate interaction or hybridization energy; that is, the energy required for two RNAs to interact via base pairing.22–24 An interaction is favorable when the energy is negative, which means that it is a spontaneous interaction,45 and is used as a proxy to interaction stability.27 One such model is RNAhybrid.46 However, RNAhybrid does not consider the secondary structure of the two RNAs before they interact.46 Importantly, the binding site that is predicted in RNAhybrid may be concealed and unable to interact in its natural secondary structure, generating a false positive ‘hit’. Other tools take secondary structure into account by considering the ‘accessible energy’ of the interaction, which is the energy needed to make the binding site accessible or open,22 23 27 47 or the probability that the binding site is unpaired in the secondary structure.24 In these models, the reported energy is typically the sum of accessible energy and hybridization or interaction energy,22 23 27 with the lowest energy being a more favorable or plausible interaction.

RNAup,23 RNAplex,24 and IntaRNA22 27 are designed for short non-coding RNA interaction prediction and consider both hybridization and accessible energy. Generally, at least one sRNA and putative target are required. These targets may be a few mRNAs to a whole transcriptome or genome, as in IntaRNA.22 27 RNAplex and IntaRNA are able to be applied to both bacterial sRNA and miRNA target prediction.24 27 Umu and Gardner performed a benchmark study of 15 RNA-RNA interaction prediction tools comparing datasets of eukaryotic, bacterial, and archaeal short RNAs along with their RNA targets.48 Overall, energy-based and accessibility-based tools performed the best in terms of true positive rates (sensitivity) and positive predictive values (precision).48 Specifically, RNAup overall performed the best, followed by IntaRNA, and RNAplex.48

From the developers of GLASSgo, IntaRNA is an energy-based and accessibility-based RNA-RNA interaction prediction tool that predicts both prokaryotic49 and eukaryotic interactions.50 Recently, IntaRNA has been incorporated into ViRBase V.3.0, a virus-host interaction database, to aid in interspecies RNA-RNA interaction binding site prediction.51 Uniquely, in addition to outputting the potential binding sites and energy scores as RNAplex and RNAup do, IntaRNA provides minimal energy profiles that allow visual exploration of interaction patterns.27

RIsearch247 facilitates exploration of predicted sRNA targets by integrating experimental transcriptomic data. This pipeline assumes that there are target (intended) interactions and off-target interactions, and calculates targeting probabilities for each potential target of an sRNA using interaction probabilities derived from both energies.47 As output, RIsearch2 reports the predicted interaction target, the start and stop positions with a visualized binding site, and the energy of the interaction.52

Similarly, RNAnue53 is a bioinformatic pipeline that analyzes experimental RNA-RNA interaction data derived from crosslinking and ligation methods (discussed in section 3). RNAnue considers hybridization and accessible energy, and complementarity of RNA partners.53 The main outputs of RNAnue are alignments of the RNA partners (chimeras, discussed further in section 3) to the genome and complementarity and energy scores.53 Overall, when benchmarked using experimental data produced by the methods, the RNAnue pipeline outperformed the original bioinformatic pipelines.53

Bioinformatic tools are able to predict interkingdom RNA-RNA interactions

Dual RNA sequencing (dRNA-Seq) is simply the simultaneous sequencing of two or more organisms without enriching or depleting any RNA type.54 After sequencing, reads are aligned to human and bacterial genomes in parallel, and subsequent differential expression, pathway, or network analysis, and more can be performed to detect sRNAs in the dataset and determine their targets.55 56 Westermann et al used dRNA-Seq and differential gene expression (DGE) and correlation analyses to identify a novel Salmonella enterica sRNA, PinT, after infecting HeLa cells.55 56 Through this strategy, they also determined that PinT regulates bacterial virulence factors which impact host signaling.55 56 Although this sRNA was not shown to have a direct human RNA target, this technique is still applicable to identifying direct interspecies RNA-RNA interactions.

Ingenuity pathway analysis (IPA)57 is a bioinformatic tool that allows transcriptome integration and pathway visualization facilitating sRNA target exploration and has been employed to predict bacterial sRNA-induced human gene expression changes.18 IPA was developed for eukaryotic transcriptome data (such as that generated from dRNA-Seq58), as well as metabolome and proteome data analysis. IPA integrates experimentally validated interactions from miRNA-mRNA interaction databases, and TargetScan, a miRNA-mRNA interaction prediction tool.59 Using IPA, Mika-Gospodorz et al identified differentially expressed genes and affected pathways in host cells infected with an intracellular bacterium,60 while Koeppen et al analyzed which human proteins were impacted by the addition of their synthetic sRNA of interest in human bronchial epithelial cells.18 Using a platform such as IPA to analyze dual RNA sequencing data54,58 61 would facilitate target prediction and contextualization within the biological system. TargetScan (currently at release V.8.0) predicts miRNA targets by scanning 6–8 nt long sites on mRNAs complementary to the miRNA.62 63 It requires the selection of a species supported by TargetScan in addition to either a gene or transcript ID and/or select a conserved miRNA family, miRBase annotation, or that the user enter the name of a known miRNA.62 63

Koeppen et al also predicted interkingdom RNA-RNA interactions between bacterial sRNAs and human mRNAs using miRanda.18 miRanda computes miRNA targets by assessing complementarity of the two input sequences and then calculating free energy of the interaction; the lowest (most favorable) energy of the resulting linked sequence structure is computed64 miRanda also considers conservation of the miRNA by scanning it against 3’ UTR datasets.64 Additionally, Choi et al used DIANA-microT for the same purpose.32 Using alignment of the driver sequence (ie, the first 9 nt of the miRNA, to databases of targets’ 3’ UTRs65 and coding sequences66), the algorithm identifies miRNA recognition elements, which are 7–9 consecutive and complementary nucleotides to the driver sequence.65 DIANA-microT uses RNAhybrid to compute hybridization energy for sites with fewer than seven consecutive complementary nucleotides.65 Pairings are ranked based on conservation of the miRNA recognition elements in several species.65 Additionally, DIANA-microT contains integrated pipelines that can, for example, analyze expression data or perform enrichment analysis of miRNA-target pairings.66 While all three tools were designed to predict human miRNA targets, Koeppen et al and Choi et al have established their use in predicting interspecies—even more, interkingdom—RNA-RNA interactions (ie, human targets of bacterial sRNAs).18 32

Section 3: experimental approaches for identifying and validating bacterial sRNA-human mRNA interactions

Use of two known interacting RNA partners to validate an interaction

Originally used to investigate nucleic acid-protein interactions, an electrophoretic mobility shift assay (EMSA) can validate the interaction of two suspected RNA partners. This first involves in vitro transcription of the two RNAs and labeling of one of them for detection, such as radiolabel49 67 68 or digoxigenin.69 The two RNAs are then incubated and shifted on a gel together; if the two RNAs interact and form a duplex, a second radiolabeled band will form and be visible.49 67–69

Reporter systems are also helpful in determining the interaction of two known RNAs. Furuse et al used a Renilla luciferase (Rluc) reporter assay to determine if MM-H would bind and repress two perfectly complementary sequences.29 One modification to this approach is to clone a potential human target transcript into the Rluc vector, providing a clearer picture of if MM-H can repress a human transcript. Also, in lacZ reporter system by Mai et al, derepression of the sRNA will decrease beta-galactosidase activity if it interacts and binds with the candidate transcript’s 5’ UTR.49 Mai et al observed this decrease in beta-galactosidase activity after 6C sRNA was derepressed, indicating its binding to its target, the lacZ-fused panD transcript.49

Alteration of sequence and expression can similarly determine if two known RNAs interact. Zhang et al, interested in cotton-fungal interactions, identified several cotton miRNAs in infecting fungal hyphae through deep sequencing of sRNA.70 Zhang et al mutated fungal transcripts to eliminate base pairing with the identified miRNAs.70 After transformation into the fungus, expression and cotton infection, the resistant transcripts accumulated while the wild-type transcripts were silenced, demonstrating direct targeting by the cotton miRNAs.70 Similarly, Jiao and Peng overexpressed and knocked-down a particular wheat miRNA (using a viral-transformed bacterial infiltration method); miRNA overexpression resulted in downregulation of target fungal transcripts, while knocked-down miRNA allowed an accumulation of target transcripts.71

RNA pull-down assays can demonstrate direct RNA-RNA interaction between known RNAs.72 Here, both interacting partners must be known, as transfection of one known RNA occurs and the end point is reverse transcription-quantitative PCR (RT qPCR), which requires primers of both sequences for detection. In this technique, a biotin-labelled known RNA of interest is transfected into cells and after lysis, the cells are added to magnetic streptavidin beads.72 The known biotinylated RNA base pairing to its suspected partner are pulled down using the beads, extracted, then amplified with RT qPCR using primers of both sequences.72 For Yang et al, base-pairing without crosslinking was enough to pull-down the pair and allow detection of the target RNA via RT-qPCR.72 This method is helpful in validating a few RNA-RNA interactions at a time, but again the interacting pairs need to be known beforehand.

Targeted identification to enable identification of one or both RNA partners

‘MS2-affinity purification coupled with RNA sequencing’ (MAPS) requires knowledge of one of the two sequences.73–75 First, an MS2 (a coat protein from a bacteriophage) aptamer is fused to a sequence of interest and expressed in vivo.73–75 Lysates undergo affinity chromatography using the MS2 maltose binding protein, which binds the MS2 aptamer with high affinity and specificity,75 to enrich RNAs fused to the MS2 aptamer and those base pairing with fused RNAs.73–75 RNA is then extracted and sequenced.73–75 This method involves manipulation of the native RNA of interest, which could potentially impact the secondary structure of the RNA or its binding. In contrast to RNA pull-down assays, only one interacting partner needs to be known as MAPS uses sequencing of the associated RNAs instead of RT qPCR.

Furuse et al and Han et al used RIP-Seq to identify microRNA-sized bacterial sRNAs associated with mammalian RISC.29 33 After infection with intracellular bacteria29 or dosage with EVs,33 cells are lysed and incubated with beads coated with anti-Ago antibodies. RNA is then extracted and subjected to deep sRNA-Seq.29 33 This technique allows identification of mammalian RISC-associated bacterial sRNAs, and previously discussed bioinformatic tools may be employed to help with this.

Crosslinking and/or ligation methods to enable global or near-global identification of interacting RNA partners

As some RNA-RNA interactions are facilitated by certain proteins, crosslinking RNA to proteins and/or ligation to the linked RNA with subsequent pull-down of specific anchor proteins known to facilitate or be involved with RNA-RNA interactions can aid in identifying targets of sRNAs that use that protein. Generally, techniques that involve both crosslinking and ligation follow these steps: (1) cells are grown and nearby RNAs or proteins are crosslinked; (2) cells are lysed and RNA is extracted; (3) crosslinked RNAs are ligated at one or both ends, forming a chimera; (4) crosslinks are reversed; (5) RNA is purified, prepared for sequencing and sequenced (figure 1). Timing of RNA pull-down depends on the technique. Crosslinking can be achieved by several molecules, but psoralen-derivative 4’-aminomethyltrioxsalen (AMT) is commonly used.76–80 AMT inserts itself in nucleic acids and on irradiation with 365 nm UV light, AMT will covalently bind two nearby RNAs together (intermolecularly and intramolecularly), with a preference for thymine or uridine,81 82 which may introduce bias for regions containing more uridine.83 Formaldehyde crosslinks RNA and protein to protein, while disuccinimidyl glutarate crosslinks protein to protein.80 Crosslinking RNA with protein, and protein to protein, enables inclusion of RNA whose direct interactions are facilitated by proteins, and also RNA indirectly involved through peripheral proteins.80 Conveniently, these crosslinks are reversed on irradiation with 254 nm UV light. To our knowledge, the methods discussed in this section have not been applied to identify direct RNA-RNA interactions between bacteria and humans.

Figure 1
  • Download figure
  • Open in new tab
  • Download powerpoint
Figure 1

Summary and comparison of crosslinking and ligation methods. Crosslinking, RNA extraction, ligation, crosslink reversal, and RNA purification and sequencing represent a core set of steps shared by many discussed methods and are indicated with images. For each method, the protocol is listed stepwise in descending order, starting at crosslinking. Differences between the methods are indicated with text as they appear in the protocol. RNA purification and sequencing (not depicted) is always the final step in the protocol. 2D, two-dimensional; AMT, 4’-aminomethyltrioxsalen; DSG, disuccinimidyl glutarate; FA, formaldehyde; LIGR-Seq, ligation of interacting RNA followed by high throughput sequencing; PARIS, psoralen analysis of RNA interactions and structures; RAP-RNA, RNA antisense purification to systemically map RNA-RNA interactions; RIC-Seq, RNA in situ conformation; RIL-Seq, RNA interaction by ligation and sequencing; rRNA, ribosomal RNA.

Gay et al debuted a ‘quick cross-linking and sequencing of hybrids’ (qCLASH) that enables the pull-down of Ago-crosslinked RNAs.84 In this technique, only UV irradiation crosslinks RNAs to proteins. Specifically, Ago proteins, along with the crosslinked RNAs, are then immunoprecipitated on beads.84 RNA pairs are ligated while on the beads and Ago is digested, while RNA is purified and sequenced.84 In a similar CLASH, Hfq, a bacterial chaperone that interacts with sRNAs, and associated RNAs are pulled down.85

‘RNA interaction by ligation and sequencing’ (RIL-Seq) similarly targets Hfq and has very similar general steps to qCLASH and the Hfq-CLASH.86 RIL-Seq, however, requires Hfq to express a Flag tag to facilitate its precipitation.86 The computational portion of RIL-Seq (and CLASH) uses the ends of sequence fragments to map to the genome, as some of these reads are chimeras made of two RNA molecules that could be from different loci in the genome (or different genomes altogether, if looking at interspecies RNA-RNA interactions), and identifies statistically significant over-represented chimeras.86

‘RNA antisense purification to systemically map RNA-RNA interactions’ (RAP-RNA) uses biotinylated DNA probes after utilizing only crosslinking.80 These probes are antisense to an RNA of interest. There are three versions of RAP-RNA that differ mainly in type of crosslinking molecule: RAP-RNA (FA), (FA+DSG) or (AMT).80 In RAP-RNA (AMT), RNA-RNA complexes are crosslinked with AMT, and RNA antisense to the biotinylated probes are sequenced.80 In this technique, the crosslinks are not reversed, and reverse transcription is halted at the crosslinked RNA which enables mapping of RNA-RNA interaction sites, as opposed to just the identification of interaction partners.

To positively select crosslinked interaction pairs, ‘sequencing of psoralen crosslinked, ligated, and selected hybrids’ (SPLASH) uses biotinylated psoralen to crosslink interacting RNAs, enabling pull-down of the RNA complexes by streptavidin-coated beads before ligation.79 Similarly, ‘RNA in situ conformation’ (RIC-Seq) labels interaction pairs with biotinylated cytidine (bis) phosphate (pCp) at the 3’ ends of transcripts before ligation of both ends of the pair, which are then positively selected for using biotin-pCp.87

For clear delineation between and enrichment for ligated interaction pairs after sequencing, ‘mapping RNA interactome in vivo’ (MARIO) ligates unambiguous biotin-tagged RNA linkers between the RNA pairs.88 First, all RNA-protein and protein-protein complexes are crosslinked via UV irradiation for capture of RNA-RNA interactions facilitated by any protein, and the proteins are denatured while the RNA is fragmented.88 Biotinylation occurs at cysteine residues and is used to immobilize the RNA-protein complexes on streptavidin beads. Next, a biotin-tagged RNA linker is ligated to the 5’ ends of RNAs, and then the interacting RNA pairs are ligated to form the chimeric product ‘RNA1-linker-RNA2’.88 Nguyen et al developed ‘MARIO tools’ (available at: https://mariotools.ucsd.edu), which is a bioinformatic pipeline to aid in the analysis of MARIO output and determine enriched interacting pairs from different genomic locations.88

Some techniques omit protein crosslinking and opt for only RNA-RNA crosslinking. Some techniques use the same general crosslinking and ligation procedure summarized above without pull-down of any kind and some minor additions. In modified CLASH by Liu et al, extracted RNAs are size-selected before ligation.76 In ‘ligation of interacting RNA followed by high throughput sequencing’ (LIGR-Seq), RNase R digests uncrosslinked (likely non-interacting) RNA and thus enriches crosslinked duplexes after ligation.77 Additionally, in ‘psoralen analysis of RNA interactions and structures’ (PARIS) proteins are totally digested, RNA is partially digested, and crosslinked RNA is purified through gel electrophoresis.78 RNAnue, mentioned in the previous section, outperformed the original analyzation pipelines for SPLASH-generated, LIGR-Seq-generated, and PARIS-generated datasets.53 This or a program like it would be ideal for analyzation of data from crosslinking/ligation methods that generate chimeras.

Various techniques use both photoactive crosslinking and subsequent ligation of RNA partners to determine direct RNA-RNA interactions without the use of anchors such as RNA-binding proteins or complementary oligonucleotide baits. Crosslinking RNAs to each other, instead of anchoring proteins, and ligation of those partners enables investigation into a global picture of RNA interaction, including protein-independent interaction.76 In contrast to some methods mentioned above, most of these methods do not require the use of genetically modified proteins or systems in order to investigate RNA-RNA interactions. Furthermore, these techniques do not require an RNA of interest, as they are able to map global or near-global RNA-RNA interactions. However, as the RNA ligation step is inefficient and ligation of non-interacting RNA occurs,83 the number of chimeric reads is typically low. Schönberger et al83 offer several improvements to these methods, including introducing non-complementary overhangs to improve ligation efficiency and a combinatorial approach to enriching crosslinked RNA partners.

These crosslinking and ligation techniques could be used to identify interspecies or interkingdom RNA-RNA interactions—in particular, bacterial sRNA-human mRNA interactions—if employed either after exposure of a singular sRNA of interest or a mixed population of sRNAs (such as that inside bacterial EVs), as these techniques identify RNA-RNA interactions globally. While the main advantage of these techniques is the global or near-global identification of RNA pairing partners, one disadvantage is the possibility of a high number of pairing partners that need to be evaluated. However, excluding interactions that do not involve interspecies pairings lowers that number considerably, especially since a small proportion of such interactions are expected.29 33

Conclusion

Bacterial sRNAs are frequently found within EVs.9 11–13 These EVs can affect immune system-related gene expression,32 33 40 and bacterial sRNAs associate with eukaryotic miRNA machinery.29 33 While not a complete list of bioinformatic and experimental approaches that exist, this review summarizes a wide range of techniques that can be employed or modified for use to identify and validate the human targets of newly identified bacterial sRNAs (and more broadly, interspecies RNA-RNA interactions).

We discussed bioinformatic tools that allow annotation of sRNAs via aligning to known sRNAs (alignment/homology tools), and those that consider energy and accessibility to assess interaction favorability (summarized in table 1). Although the latter are more thorough and consider secondary structure, using a combination of bioinformatic techniques to identify sRNA targets is the most robust approach. Even when a particular target prediction tool is designed for a certain biological system, it may be applicable elsewhere—such as IntaRNA, which was originally intended for sRNA target prediction,22 but has been used for eukaryotic target prediction as well.50 Similarly, miRNA target prediction tools, such as TargetScan (IPA), DIANA microT, and miRanda, have been successfully used to predict human targets of bacterial sRNAs.18 32

Experimental techniques to identify RNA-RNA interactions that require the use of an anchor (eg, a biotinylated probe) and those that employ crosslinking and ligation of interacting RNA pairs were also discussed (summarized in table 2). EMSA, reporter systems, altering expression, and RNA pull-down assays are useful for validation of interactions, whereas crosslinking and ligation methods are more global and will be useful for exploration and identification of RNA-RNA interaction networks. While experimental methods to verify RNA-RNA interactions have often been described for use in one cell type or species, they can be adapted in experimental design to allow identification and validation of interkingdom RNA-RNA interactions. For example, the aforementioned crosslinking and ligation methods could be employed after cell exposure to bacterial sRNAs or EVs.

To our knowledge, the vast majority of the bioinformatic and experimental methods discussed here have not been used to analyze interspecies interactions, particularly bacteria-human. Employing these methods, particularly the experimental methods, to determine bacterial-human RNA-RNA interactions will facilitate the important research that is needed to determine the potential impact microbiota exert on the host through these interspecies RNA-RNA interactions. It is only with improved understanding of these interactions that we can begin to effectively alter or adapt microbiota and positively impact human disease.

Ethics statements

Patient consent for publication

Not applicable.

Ethics approval

Not applicable.

Footnotes

  • Contributors KIK, RJF, and ALH were all involved in the concept development. KIK wrote the first draft of the manuscript. RJF and ALH provided critical review for intellectual content. All authors approve of the final version as written.

  • Funding The authors have not declared a specific grant for this research from any funding agency in the public, commercial or not-for-profit sectors.

  • Competing interests None declared.

  • Provenance and peer review Commissioned; externally peer reviewed.

http://creativecommons.org/licenses/by-nc/4.0/

This is an open access article distributed in accordance with the Creative Commons Attribution Non Commercial (CC BY-NC 4.0) license, which permits others to distribute, remix, adapt, build upon this work non-commercially, and license their derivative works on different terms, provided the original work is properly cited, an indication of whether changes were made, and the use is non-commercial. See: http://creativecommons.org/licenses/by-nc/4.0/.

References

  1. ↵
    1. Sultan S,
    2. El-Mowafy M,
    3. Elgaml A, et al
    . Metabolic influences of gut microbiota dysbiosis on inflammatory bowel disease. Front Physiol 2021;12:715506.doi:10.3389/fphys.2021.715506pmid:http://www.ncbi.nlm.nih.gov/pubmed/34646151
    OpenUrlPubMed
  2. ↵
    1. Hufnagl K,
    2. Pali-Schöll I,
    3. Roth-Walter F, et al
    . Dysbiosis of the gut and lung microbiome has a role in asthma. Semin Immunopathol 2020;42:75–93.doi:10.1007/s00281-019-00775-ypmid:http://www.ncbi.nlm.nih.gov/pubmed/32072252
    OpenUrlPubMed
  3. ↵
    1. Saraf VS,
    2. Sheikh SA,
    3. Ahmad A, et al
    . Vaginal microbiome: normalcy vs dysbiosis. Arch Microbiol 2021;203:3793–802.doi:10.1007/s00203-021-02414-3pmid:http://www.ncbi.nlm.nih.gov/pubmed/34120200
    OpenUrlPubMed
  4. ↵
    1. Macia L,
    2. Nanan R,
    3. Hosseini-Beheshti E, et al
    . Host- and Microbiota-Derived extracellular vesicles, immune function, and disease development. Int J Mol Sci 2019;21:201.doi:10.3390/ijms21010107pmid:http://www.ncbi.nlm.nih.gov/pubmed/31877909
    OpenUrlPubMed
  5. ↵
    1. Ñahui Palomino RA,
    2. Vanpouille C,
    3. Costantini PE, et al
    . Microbiota-host communications: bacterial extracellular vesicles as a common language. PLoS Pathog 2021;17:e1009508.doi:10.1371/journal.ppat.1009508pmid:http://www.ncbi.nlm.nih.gov/pubmed/33984071
    OpenUrlPubMed
  6. ↵
    1. Pérez-Cruz C,
    2. Delgado L,
    3. López-Iglesias C, et al
    . Outer-inner membrane vesicles naturally secreted by gram-negative pathogenic bacteria. PLoS One 2015;10:e0116896.doi:10.1371/journal.pone.0116896pmid:http://www.ncbi.nlm.nih.gov/pubmed/25581302
    OpenUrlCrossRefPubMed
  7. ↵
    1. Bose S,
    2. Aggarwal S,
    3. Singh DV, et al
    . Extracellular vesicles: an emerging platform in Gram-positive bacteria. Microb Cell 2020;7:312–22.doi:10.15698/mic2020.12.737pmid:http://www.ncbi.nlm.nih.gov/pubmed/33335921
    OpenUrlPubMed
  8. ↵
    1. Théry C,
    2. Witwer KW,
    3. Aikawa E, et al
    . Minimal information for studies of extracellular vesicles 2018 (MISEV2018): a position statement of the International Society for extracellular vesicles and update of the MISEV2014 guidelines. J Extracell Vesicles 2018;7:1535750.doi:10.1080/20013078.2018.1535750pmid:http://www.ncbi.nlm.nih.gov/pubmed/30637094
    OpenUrlPubMed
  9. ↵
    1. Ghosal A,
    2. Upadhyaya BB,
    3. Fritz JV, et al
    . The extracellular RNA complement of Escherichia coli. Microbiologyopen 2015;4:252–66.doi:10.1002/mbo3.235pmid:http://www.ncbi.nlm.nih.gov/pubmed/25611733
    OpenUrlCrossRefPubMed
  10. ↵
    1. Chamberlain CA,
    2. Hatch M,
    3. Garrett TJ
    . Extracellular vesicle analysis by paper spray ionization mass spectrometry. Metabolites 2021;11:308.doi:10.3390/metabo11050308pmid:http://www.ncbi.nlm.nih.gov/pubmed/34065030
    OpenUrlPubMed
  11. ↵
    1. Blenkiron C,
    2. Simonov D,
    3. Muthukaruppan A, et al
    . Uropathogenic Escherichia coli releases extracellular vesicles that are associated with RNA. PLoS One 2016;11:e0160440.doi:10.1371/journal.pone.0160440pmid:http://www.ncbi.nlm.nih.gov/pubmed/27500956
    OpenUrlPubMed
  12. ↵
    1. Joshi B,
    2. Singh B,
    3. Nadeem A, et al
    . Transcriptome Profiling of Staphylococcus aureus Associated Extracellular Vesicles Reveals Presence of Small RNA-Cargo. Front Mol Biosci 2020;7:566207.doi:10.3389/fmolb.2020.566207pmid:http://www.ncbi.nlm.nih.gov/pubmed/33521050
    OpenUrlPubMed
  13. ↵
    1. Tsatsaronis JA,
    2. Franch-Arroyo S,
    3. Resch U, et al
    . Extracellular vesicle RNA: a universal mediator of microbial communication? Trends Microbiol 2018;26:401–10.doi:10.1016/j.tim.2018.02.009pmid:http://www.ncbi.nlm.nih.gov/pubmed/29548832
    OpenUrlCrossRefPubMed
  14. ↵
    1. Beisel CL,
    2. Storz G
    . Base pairing small RNAs and their roles in global regulatory networks. FEMS Microbiol Rev 2010;34:866–82.doi:10.1111/j.1574-6976.2010.00241.xpmid:http://www.ncbi.nlm.nih.gov/pubmed/20662934
    OpenUrlCrossRefPubMedWeb of Science
  15. ↵
    1. Nakanishi K
    . Anatomy of RISC: how do small RNAs and chaperones activate Argonaute proteins? Wiley Interdiscip Rev RNA 2016;7:637–60.doi:10.1002/wrna.1356pmid:http://www.ncbi.nlm.nih.gov/pubmed/27184117
    OpenUrlPubMed
  16. ↵
    1. Fritsch TE,
    2. Siqueira FM,
    3. Schrank IS
    . Global analysis of sRNA target genes in Mycoplasma hyopneumoniae. BMC Genomics 2018;19:767.doi:10.1186/s12864-018-5136-5pmid:http://www.ncbi.nlm.nih.gov/pubmed/30352553
    OpenUrlPubMed
  17. ↵
    1. Miravet-Verde S,
    2. Lloréns-Rico V,
    3. Serrano L
    . Alternative transcriptional regulation in genome-reduced bacteria. Curr Opin Microbiol 2017;39:89–95.doi:10.1016/j.mib.2017.10.022pmid:http://www.ncbi.nlm.nih.gov/pubmed/29154025
    OpenUrlPubMed
  18. ↵
    1. Koeppen K,
    2. Hampton TH,
    3. Jarek M, et al
    . A novel mechanism of host-pathogen interaction through sRNA in bacterial outer membrane vesicles. PLoS Pathog 2016;12:e1005672.doi:10.1371/journal.ppat.1005672pmid:http://www.ncbi.nlm.nih.gov/pubmed/27295279
    OpenUrlCrossRefPubMed
  19. ↵
    1. Carrier M-C,
    2. Lalaouna D,
    3. Massé E
    . Broadening the definition of bacterial small RNAs: characteristics and mechanisms of action. Annu Rev Microbiol 2018;72:141–61.doi:10.1146/annurev-micro-090817-062607pmid:http://www.ncbi.nlm.nih.gov/pubmed/30200848
    OpenUrlCrossRefPubMed
  20. ↵
    1. Obregón-Henao A,
    2. Duque-Correa MA,
    3. Rojas M, et al
    . Stable extracellular RNA fragments of Mycobacterium tuberculosis induce early apoptosis in human monocytes via a caspase-8 dependent mechanism. PLoS One 2012;7:e29970.doi:10.1371/journal.pone.0029970pmid:http://www.ncbi.nlm.nih.gov/pubmed/22253841
    OpenUrlCrossRefPubMed
  21. ↵
    1. RKr G,
    2. Luong TT,
    3. Lee CY
    . RNAIII of the Staphylococcus aureus agr system activates global regulator MgrA by stabilizing mRNA. Proceedings of the National Academy of Sciences 2015;112:14036–41.
    OpenUrlAbstract/FREE Full Text
  22. ↵
    1. Busch A,
    2. Richter AS,
    3. Backofen R
    . IntaRNA: efficient prediction of bacterial sRNA targets incorporating target site accessibility and seed regions. Bioinformatics 2008;24:2849–56.doi:10.1093/bioinformatics/btn544pmid:http://www.ncbi.nlm.nih.gov/pubmed/18940824
    OpenUrlCrossRefPubMedWeb of Science
  23. ↵
    1. Mückstein U,
    2. Tafer H,
    3. Hackermüller J, et al
    . Thermodynamics of RNA-RNA binding. Bioinformatics 2006;22:1177–82.doi:10.1093/bioinformatics/btl024pmid:http://www.ncbi.nlm.nih.gov/pubmed/16446276
    OpenUrlCrossRefPubMedWeb of Science
  24. ↵
    1. Tafer H,
    2. Amman F,
    3. Eggenhofer F, et al
    . Fast accessibility-based prediction of RNA-RNA interactions. Bioinformatics 2011;27:1934–40.doi:10.1093/bioinformatics/btr281pmid:http://www.ncbi.nlm.nih.gov/pubmed/21593134
    OpenUrlCrossRefPubMedWeb of Science
  25. ↵
    1. Wenzel A,
    2. Akbasli E,
    3. Gorodkin J
    . RIsearch: fast RNA-RNA interaction search using a simplified nearest-neighbor energy model. Bioinformatics 2012;28:2738–46.doi:10.1093/bioinformatics/bts519pmid:http://www.ncbi.nlm.nih.gov/pubmed/22923300
    OpenUrlCrossRefPubMedWeb of Science
  26. ↵
    1. Fukunaga T,
    2. Hamada M
    . RIblast: an ultrafast RNA-RNA interaction prediction system based on a seed-and-extension approach. Bioinformatics 2017;33:2666–74.doi:10.1093/bioinformatics/btx287pmid:http://www.ncbi.nlm.nih.gov/pubmed/28459942
    OpenUrlPubMed
  27. ↵
    1. Mann M,
    2. Wright PR,
    3. Backofen R
    . IntaRNA 2.0: enhanced and customizable prediction of RNA-RNA interactions. Nucleic Acids Res 2017;45:W435–9.doi:10.1093/nar/gkx279pmid:http://www.ncbi.nlm.nih.gov/pubmed/28472523
    OpenUrlCrossRefPubMed
  28. ↵
    1. Andrade JM,
    2. Pobre V,
    3. Arraiano CM
    . Small RNA modules confer different stabilities and interact differently with multiple targets. PLoS One 2013;8:e52866.doi:10.1371/journal.pone.0052866pmid:http://www.ncbi.nlm.nih.gov/pubmed/23349691
    OpenUrlCrossRefPubMed
  29. ↵
    1. Furuse Y,
    2. Finethy R,
    3. Saka HA, et al
    . Search for microRNAs expressed by intracellular bacterial pathogens in infected mammalian cells. PLoS One 2014;9:e106434.doi:10.1371/journal.pone.0106434pmid:http://www.ncbi.nlm.nih.gov/pubmed/25184567
    OpenUrlCrossRefPubMed
  30. ↵
    1. Choi J-W,
    2. Kwon T-Y,
    3. Hong S-H, et al
    . Isolation and characterization of a microRNA-size secretable small RNA in Streptococcus sanguinis. Cell Biochem Biophys 2018;76:293–301.doi:10.1007/s12013-016-0770-5pmid:http://www.ncbi.nlm.nih.gov/pubmed/27796789
    OpenUrlPubMed
  31. ↵
    1. Kang S-M,
    2. Choi J-W,
    3. Lee Y, et al
    . Identification of microRNA-size, small RNAs in Escherichia coli. Curr Microbiol 2013;67:609–13.doi:10.1007/s00284-013-0411-9pmid:http://www.ncbi.nlm.nih.gov/pubmed/23783561
    OpenUrlCrossRefPubMed
  32. ↵
    1. Choi J-W,
    2. Kim S-C,
    3. Hong S-H, et al
    . Secretable small RNAs via outer membrane vesicles in periodontal pathogens. J Dent Res 2017;96:458–66.doi:10.1177/0022034516685071pmid:http://www.ncbi.nlm.nih.gov/pubmed/28068479
    OpenUrlCrossRefPubMed
  33. ↵
    1. Han E-C,
    2. Choi S-Y,
    3. Lee Y, et al
    . Extracellular RNAs in periodontopathogenic outer membrane vesicles promote TNF-α production in human macrophages and cross the blood-brain barrier in mice. Faseb J 2019;33:13412–22.doi:10.1096/fj.201901575Rpmid:http://www.ncbi.nlm.nih.gov/pubmed/31545910
    OpenUrlCrossRefPubMed
  34. ↵
    1. Malabirade A,
    2. Habier J,
    3. Heintz-Buschart A, et al
    . The RNA Complement of Outer Membrane Vesicles From Salmonella enterica Serovar Typhimurium Under Distinct Culture Conditions. Front Microbiol 2018;9:2015.doi:10.3389/fmicb.2018.02015pmid:http://www.ncbi.nlm.nih.gov/pubmed/30214435
    OpenUrlPubMed
  35. ↵
    1. Ormseth MJ,
    2. Wu Q,
    3. Zhao S, et al
    . Circulating microbial small RNAs are altered in patients with rheumatoid arthritis. Ann Rheum Dis 2020;79:1557–64.doi:10.1136/annrheumdis-2020-217589pmid:http://www.ncbi.nlm.nih.gov/pubmed/32958509
    OpenUrlAbstract/FREE Full Text
  36. ↵
    1. Beatty M,
    2. Guduric-Fuchs J,
    3. Brown E, et al
    . Small RNAs from plants, bacteria and fungi within the order Hypocreales are ubiquitous in human plasma. BMC Genomics 2014;15:933.doi:10.1186/1471-2164-15-933pmid:http://www.ncbi.nlm.nih.gov/pubmed/25344700
    OpenUrlPubMed
  37. ↵
    1. Wang K,
    2. Li H,
    3. Yuan Y, et al
    . The complex exogenous RNA spectra in human plasma: an interface with human gut biota? PLoS One 2012;7:e51009.doi:10.1371/journal.pone.0051009pmid:http://www.ncbi.nlm.nih.gov/pubmed/23251414
    OpenUrlCrossRefPubMed
  38. ↵
    1. Tulkens J,
    2. De Wever O,
    3. Hendrix A
    . Analyzing bacterial extracellular vesicles in human body fluids by orthogonal biophysical separation and biochemical characterization. Nat Protoc 2020;15:40–67.doi:10.1038/s41596-019-0236-5pmid:http://www.ncbi.nlm.nih.gov/pubmed/31776460
    OpenUrlPubMed
  39. ↵
    1. Krohmaly K,
    2. Hoptay C,
    3. Hahn A
    . Detection of bacterial extracellular vesicles in blood from healthy volunteers. Poster presented at: ISEV2020 Annual Meeting, 2020.
  40. ↵
    1. Ha JY,
    2. Choi S-Y,
    3. Lee JH, et al
    . Delivery of periodontopathogenic extracellular vesicles to brain monocytes and microglial IL-6 promotion by RNA cargo. Front Mol Biosci 2020;7:596366.doi:10.3389/fmolb.2020.596366pmid:http://www.ncbi.nlm.nih.gov/pubmed/33330627
    OpenUrlPubMed
  41. ↵
    1. Langmead B,
    2. Trapnell C,
    3. Pop M, et al
    . Ultrafast and memory-efficient alignment of short DNA sequences to the human genome. Genome Biol 2009;10:R25.doi:10.1186/gb-2009-10-3-r25pmid:http://www.ncbi.nlm.nih.gov/pubmed/19261174
    OpenUrlCrossRefPubMed
  42. ↵
    1. Jorth P,
    2. Whiteley M
    . Characterization of a novel riboswitch-regulated lysine transporter in Aggregatibacter actinomycetemcomitans. J Bacteriol 2010;192:6240–50.doi:10.1128/JB.00935-10pmid:http://www.ncbi.nlm.nih.gov/pubmed/20889741
    OpenUrlAbstract/FREE Full Text
  43. ↵
    1. Kalvari I,
    2. Nawrocki EP,
    3. Ontiveros-Palacios N, et al
    . Rfam 14: expanded coverage of metagenomic, viral and microRNA families. Nucleic Acids Res 2021;49:D192–200.doi:10.1093/nar/gkaa1047pmid:http://www.ncbi.nlm.nih.gov/pubmed/33211869
    OpenUrlCrossRefPubMed
  44. ↵
    1. Grüll MP,
    2. Peña-Castillo L,
    3. Mulligan ME, et al
    . Genome-wide identification and characterization of small RNAs in Rhodobacter capsulatus and identification of small RNAs affected by loss of the response regulator CtrA. RNA Biol 2017;14:914–25.doi:10.1080/15476286.2017.1306175pmid:http://www.ncbi.nlm.nih.gov/pubmed/28296577
    OpenUrlPubMed
  45. ↵
    1. Ragan C,
    2. Zuker M,
    3. Ragan MA
    . Quantitative prediction of miRNA-mRNA interaction based on equilibrium concentrations. PLoS Comput Biol 2011;7:e1001090.doi:10.1371/journal.pcbi.1001090pmid:http://www.ncbi.nlm.nih.gov/pubmed/21390282
    OpenUrlCrossRefPubMed
  46. ↵
    1. Rehmsmeier M,
    2. Steffen P,
    3. Höchsmann M, et al
    . Fast and effective prediction of microRNA/target duplexes. RNA 2004;10:1507–17.doi:10.1261/rna.5248604pmid:http://www.ncbi.nlm.nih.gov/pubmed/15383676
    OpenUrlAbstract/FREE Full Text
  47. ↵
    1. Alkan F,
    2. Wenzel A,
    3. Palasca O, et al
    . RIsearch2: suffix array-based large-scale prediction of RNA-RNA interactions and siRNA off-targets. Nucleic Acids Res 2017;45:e60.doi:10.1093/nar/gkw1325pmid:http://www.ncbi.nlm.nih.gov/pubmed/28108657
    OpenUrlPubMed
  48. ↵
    1. Umu SU,
    2. Gardner PP
    . A comprehensive benchmark of RNA–RNA interaction prediction tools for all domains of life. Bioinformatics 2016;7:btw728–96.doi:10.1093/bioinformatics/btw728
    OpenUrl
  49. ↵
    1. Mai J,
    2. Rao C,
    3. Watt J, et al
    . Mycobacterium tuberculosis 6C sRNA binds multiple mRNA targets via C-rich loops independent of RNA chaperones. Nucleic Acids Res 2019;47:4292–307.doi:10.1093/nar/gkz149pmid:http://www.ncbi.nlm.nih.gov/pubmed/30820540
    OpenUrlCrossRefPubMed
  50. ↵
    1. Zhao P,
    2. Li T,
    3. Wang Y, et al
    . Lncrna MYCNOS promotes glioblastoma cell proliferation by regulating miR-216b/FOXM1 axis. Metab Brain Dis 2021;36:1185–9.doi:10.1007/s11011-021-00729-0pmid:http://www.ncbi.nlm.nih.gov/pubmed/33871770
    OpenUrlPubMed
  51. ↵
    1. Cheng J,
    2. Lin Y,
    3. Xu L, et al
    . ViRBase v3.0: a virus and host ncRNA-associated interaction Repository with increased coverage and annotation. Nucleic Acids Res 2022;50:D928–33.doi:10.1093/nar/gkab1029pmid:http://www.ncbi.nlm.nih.gov/pubmed/34723320
    OpenUrlPubMed
  52. ↵
    1. Aklan F,
    2. Wenzel A,
    3. Palasca O, et al
    . RIsearch2: user manual, 2018. Available: https://rth.dk/resources/risearch/download/Manual_RIsearch2.pdf [Accessed 3 May 2022].
  53. ↵
    1. Schäfer RA,
    2. Voß B
    . RNAnue: efficient data analysis for RNA-RNA interactomics. Nucleic Acids Res 2021;49:5493–501.doi:10.1093/nar/gkab340pmid:http://www.ncbi.nlm.nih.gov/pubmed/34019662
    OpenUrlPubMed
  54. ↵
    1. Westermann AJ,
    2. Gorski SA,
    3. Vogel J
    . Dual RNA-seq of pathogen and host. Nat Rev Microbiol 2012;10:618–30.doi:10.1038/nrmicro2852pmid:http://www.ncbi.nlm.nih.gov/pubmed/22890146
    OpenUrlCrossRefPubMed
  55. ↵
    1. Westermann AJ,
    2. Förstner KU,
    3. Amman F, et al
    . Dual RNA-seq unveils noncoding RNA functions in host-pathogen interactions. Nature 2016;529:496–501.doi:10.1038/nature16547pmid:http://www.ncbi.nlm.nih.gov/pubmed/26789254
    OpenUrlCrossRefPubMed
  56. ↵
    1. Westermann AJ,
    2. Barquist L,
    3. Vogel J
    . Resolving host-pathogen interactions by dual RNA-seq. PLoS Pathog 2017;13:e1006033.doi:10.1371/journal.ppat.1006033pmid:http://www.ncbi.nlm.nih.gov/pubmed/28207848
    OpenUrlCrossRefPubMed
  57. ↵
    1. QIAGEN Ingenuity Pathway Analysis (QIAGEN IPA)
    . Qiagen.com. Available: https://digitalinsights.qiagen.com/products-overview/discovery-insights-portfolio/analysis-and-visualization/qiagen-ipa/ [Accessed 4 May 2022].
  58. ↵
    1. Chung M,
    2. Bruno VM,
    3. Rasko DA, et al
    . Best practices on the differential expression analysis of multi-species RNA-seq. Genome Biol 2021;22:121.doi:10.1186/s13059-021-02337-8pmid:http://www.ncbi.nlm.nih.gov/pubmed/33926528
    OpenUrlPubMed
  59. ↵
    1. Qiagen.com
    . Microrna target filter. Available: https://digitalinsights.qiagen.com/products-overview/discovery-insights-portfolio/analysis-and-visualization/qiagen-ipa/features/microrna-target-filter/ [Accessed 3 May 2022].
  60. ↵
    1. Mika-Gospodorz B,
    2. Giengkam S,
    3. Westermann AJ, et al
    . Dual RNA-seq of Orientia tsutsugamushi informs on host-pathogen interactions for this neglected intracellular human pathogen. Nat Commun 2020;11:3363.doi:10.1038/s41467-020-17094-8pmid:http://www.ncbi.nlm.nih.gov/pubmed/32620750
    OpenUrlCrossRefPubMed
  61. ↵
    1. Schulze S,
    2. Schleicher J,
    3. Guthke R, et al
    . How to predict molecular interactions between species? Front Microbiol 2016;7:442.doi:10.3389/fmicb.2016.00442pmid:http://www.ncbi.nlm.nih.gov/pubmed/27065992
    OpenUrlPubMed
  62. ↵
    1. Lewis BP,
    2. Burge CB,
    3. Bartel DP
    . Conserved seed pairing, often flanked by adenosines, indicates that thousands of human genes are microRNA targets. Cell 2005;120:15–20.doi:10.1016/j.cell.2004.12.035pmid:http://www.ncbi.nlm.nih.gov/pubmed/15652477
    OpenUrlCrossRefPubMedWeb of Science
  63. ↵
    TargetScan. Available: https://www.targetscan.org/vert_80/ [Accessed 4 May 2022].
  64. ↵
    1. Enright AJ,
    2. John B,
    3. Gaul U, et al
    . MicroRNA targets in Drosophila. Genome Biol 2003;5:R1.doi:10.1186/gb-2003-5-1-r1pmid:http://www.ncbi.nlm.nih.gov/pubmed/14709173
    OpenUrlCrossRefPubMed
  65. ↵
    1. Maragkakis M,
    2. Alexiou P,
    3. Papadopoulos GL, et al
    . Accurate microRNA target prediction correlates with protein repression levels. BMC Bioinformatics 2009;10:295.doi:10.1186/1471-2105-10-295pmid:http://www.ncbi.nlm.nih.gov/pubmed/19765283
    OpenUrlCrossRefPubMed
  66. ↵
    1. Paraskevopoulou MD,
    2. Georgakilas G,
    3. Kostoulas N, et al
    . DIANA-microT web server v5.0: service integration into miRNA functional analysis workflows. Nucleic Acids Res 2013;41:W169–73.doi:10.1093/nar/gkt393pmid:http://www.ncbi.nlm.nih.gov/pubmed/23680784
    OpenUrlCrossRefPubMedWeb of Science
  67. ↵
    1. Henderson CA,
    2. Vincent HA,
    3. Stone CM, et al
    . Characterization of MICA interactions suggests a potential novel means of gene regulation by small non-coding RNAs. Nucleic Acids Res 2013;41:3386–97.doi:10.1093/nar/gkt008pmid:http://www.ncbi.nlm.nih.gov/pubmed/23361466
    OpenUrlCrossRefPubMed
  68. ↵
    1. Silva IJ,
    2. Ortega AD,
    3. Viegas SC, et al
    . An rpoS-dependent sRNA regulates the expression of a chaperone involved in protein folding. RNA 2013;19:1253–65.doi:10.1261/rna.039537.113pmid:http://www.ncbi.nlm.nih.gov/pubmed/23893734
    OpenUrlAbstract/FREE Full Text
  69. ↵
    1. Ota C,
    2. Morisaki H,
    3. Nakata M, et al
    . Streptococcus sanguinis Noncoding cia-Dependent Small RNAs Negatively Regulate Expression of Type IV Pilus Retraction ATPase PilT and Biofilm Formation. Infect Immun 2018;86:e00894–17.doi:10.1128/IAI.00894-17pmid:http://www.ncbi.nlm.nih.gov/pubmed/29263111
    OpenUrlPubMed
  70. ↵
    1. Zhang T,
    2. Zhao Y-L,
    3. Zhao J-H, et al
    . Cotton plants export microRNAs to inhibit virulence gene expression in a fungal pathogen. Nat Plants 2016;2:16153.doi:10.1038/nplants.2016.153pmid:http://www.ncbi.nlm.nih.gov/pubmed/27668926
    OpenUrlPubMed
  71. ↵
    1. Jiao J,
    2. Peng D
    . Wheat microRNA1023 suppresses invasion of Fusarium graminearum via targeting and silencing FGSG_03101. Journal of Plant Interactions 2018;13:514–21.
    OpenUrl
  72. ↵
    1. Yang M,
    2. Chen W,
    3. Liu H, et al
    . Long non-coding RNA CBR3 antisense RNA 1 is downregulated in colorectal cancer and inhibits miR-29a-Mediated cell migration and invasion. Mol Biotechnol 2022;64:773–9.doi:10.1007/s12033-021-00444-2pmid:http://www.ncbi.nlm.nih.gov/pubmed/35107754
    OpenUrlCrossRefPubMed
  73. ↵
    1. Carrier M-C,
    2. Lalaouna D,
    3. Massé E
    . A game of tag: MAPS catches up on RNA interactomes. RNA Biol 2016;13:473–6.doi:10.1080/15476286.2016.1156830pmid:http://www.ncbi.nlm.nih.gov/pubmed/26967018
    OpenUrlPubMed
  74. ↵
    1. Mercier N,
    2. Prévost K,
    3. Massé E, et al
    . MS2-Affinity purification coupled with RNA sequencing in Gram-positive bacteria. J Vis Exp 2021. doi:doi:10.3791/61731. [Epub ahead of print: 23 Feb 2021].pmid:http://www.ncbi.nlm.nih.gov/pubmed/33720114
    OpenUrlPubMed
  75. ↵
    1. Lalaouna D,
    2. Desgranges E,
    3. Caldelari I
    . MS2-Affinity purification coupled with RNA sequencing approach in the human pathogen Staphylococcus aureus 2018:393–411.
  76. ↵
    1. Liu T,
    2. Zhang K,
    3. Xu S, et al
    . Detecting RNA-RNA interactions in E. coli using a modified clash method. BMC Genomics 2017;18:343.doi:10.1186/s12864-017-3725-3pmid:http://www.ncbi.nlm.nih.gov/pubmed/28468647
    OpenUrlPubMed
  77. ↵
    1. Sharma E,
    2. Sterne-Weiler T,
    3. O'Hanlon D, et al
    . Global mapping of human RNA-RNA interactions. Mol Cell 2016;62:618–26.doi:10.1016/j.molcel.2016.04.030pmid:http://www.ncbi.nlm.nih.gov/pubmed/27184080
    OpenUrlCrossRefPubMed
  78. ↵
    1. Lu Z,
    2. Zhang QC,
    3. Lee B, et al
    . RNA duplex map in living cells reveals higher-order transcriptome structure. Cell 2016;165:1267–79.doi:10.1016/j.cell.2016.04.028pmid:http://www.ncbi.nlm.nih.gov/pubmed/27180905
    OpenUrlCrossRefPubMed
  79. ↵
    1. Aw JGA,
    2. Shen Y,
    3. Wilm A, et al
    . In Vivo Mapping of Eukaryotic RNA Interactomes Reveals Principles of Higher-Order Organization and Regulation. Mol Cell 2016;62:603–17.doi:10.1016/j.molcel.2016.04.028pmid:http://www.ncbi.nlm.nih.gov/pubmed/27184079
    OpenUrlCrossRefPubMed
  80. ↵
    1. Engreitz JM,
    2. Sirokman K,
    3. McDonel P, et al
    . RNA-RNA interactions enable specific targeting of noncoding RNAs to nascent pre-mRNAs and chromatin sites. Cell 2014;159:188–99.doi:10.1016/j.cell.2014.08.018pmid:http://www.ncbi.nlm.nih.gov/pubmed/25259926
    OpenUrlCrossRefPubMed
  81. ↵
    1. Zhen WP,
    2. Buchardt O,
    3. Nielsen H, et al
    . Site specificity of psoralen-DNA interstrand cross-linking determined by nuclease Bal31 digestion. Biochemistry 1986;25:6598–603.doi:10.1021/bi00369a039pmid:http://www.ncbi.nlm.nih.gov/pubmed/3024713
    OpenUrlPubMed
  82. ↵
    1. Cimino GD,
    2. Gamper HB,
    3. Isaacs ST, et al
    . Psoralens as photoactive probes of nucleic acid structure and function: organic chemistry, photochemistry, and biochemistry. Annu Rev Biochem 1985;54:1151–93.doi:10.1146/annurev.bi.54.070185.005443pmid:http://www.ncbi.nlm.nih.gov/pubmed/2411210
    OpenUrlCrossRefPubMedWeb of Science
  83. ↵
    1. Schönberger B,
    2. Schaal C,
    3. Schäfer R, et al
    . Rna interactomics: recent advances and remaining challenges. F1000Res 2018;7:1824.doi:10.12688/f1000research.16146.1pmid:http://www.ncbi.nlm.nih.gov/pubmed/30519453
    OpenUrlPubMed
  84. ↵
    1. Gay LA,
    2. Turner PC,
    3. Renne R
    . Modified cross-linking, ligation, and sequencing of hybrids (qCLASH) to identify microRNA targets. Curr Protoc 2021;1:e257.doi:10.1002/cpz1.257pmid:http://www.ncbi.nlm.nih.gov/pubmed/34610213
    OpenUrlPubMed
  85. ↵
    1. Iosub IA,
    2. van Nues RW,
    3. McKellar SW, et al
    . Hfq clash uncovers sRNA-target interaction networks linked to nutrient availability adaptation. Elife 2020;9:e54655.doi:10.7554/eLife.54655pmid:http://www.ncbi.nlm.nih.gov/pubmed/32356726
    OpenUrlPubMed
  86. ↵
    1. Melamed S,
    2. Faigenbaum-Romm R,
    3. Peer A, et al
    . Mapping the small RNA interactome in bacteria using RIL-seq. Nat Protoc 2018;13:1–33.doi:10.1038/nprot.2017.115pmid:http://www.ncbi.nlm.nih.gov/pubmed/29215635
    OpenUrlCrossRefPubMed
  87. ↵
    1. Cai Z,
    2. Cao C,
    3. Ji L, et al
    . RIC-seq for global in situ profiling of RNA-RNA spatial interactions. Nature 2020;582:432–7.doi:10.1038/s41586-020-2249-1pmid:http://www.ncbi.nlm.nih.gov/pubmed/32499643
    OpenUrlCrossRefPubMed
  88. ↵
    1. Nguyen TC,
    2. Cao X,
    3. Yu P, et al
    . Mapping RNA-RNA interactome and RNA structure in vivo by Mario. Nat Commun 2016;7:12023.doi:10.1038/ncomms12023pmid:http://www.ncbi.nlm.nih.gov/pubmed/27338251
    OpenUrlCrossRefPubMed
    1. Johnson M,
    2. Zaretskaya I,
    3. Raytselis Y, et al
    . NCBI BLAST: a better web interface. Nucleic Acids Res 2008;36:W5–9.doi:10.1093/nar/gkn201pmid:http://www.ncbi.nlm.nih.gov/pubmed/18440982
    OpenUrlCrossRefPubMedWeb of Science
    1. Ermolaeva MD,
    2. Khalak HG,
    3. White O, et al
    . Prediction of transcription terminators in bacterial genomes. J Mol Biol 2000;301:27–33.doi:10.1006/jmbi.2000.3836pmid:http://www.ncbi.nlm.nih.gov/pubmed/10926490
    OpenUrlCrossRefPubMedWeb of Science
    1. Gautheret D,
    2. Lambert A
    . Direct RNA motif definition and identification from multiple sequence alignments using secondary structure profiles. J Mol Biol 2001;313:1003–11.doi:10.1006/jmbi.2001.5102pmid:http://www.ncbi.nlm.nih.gov/pubmed/11700055
    OpenUrlCrossRefPubMedWeb of Science
    1. Gardner PP,
    2. Barquist L,
    3. Bateman A, et al
    . RNIE: genome-wide prediction of bacterial intrinsic terminators. Nucleic Acids Res 2011;39:5845–52.doi:10.1093/nar/gkr168pmid:http://www.ncbi.nlm.nih.gov/pubmed/21478170
    OpenUrlCrossRefPubMedWeb of Science
    1. McClure R,
    2. Balasubramanian D,
    3. Sun Y, et al
    . Computational analysis of bacterial RNA-Seq data. Nucleic Acids Res 2013;41:e140.doi:10.1093/nar/gkt444pmid:http://www.ncbi.nlm.nih.gov/pubmed/23716638
    OpenUrlCrossRefPubMed
    1. Leonard S,
    2. Meyer S,
    3. Lacour S, et al
    . APERO: a genome-wide approach for identifying bacterial small RNAs from RNA-Seq data. Nucleic Acids Res 2019;47:e88.doi:10.1093/nar/gkz485pmid:http://www.ncbi.nlm.nih.gov/pubmed/31147705
    OpenUrlPubMed
    1. Toffano-Nioche C,
    2. Luo Y,
    3. Kuchly C, et al
    . Detection of non-coding RNA in bacteria and archaea using the DETR'PROK Galaxy pipeline. Methods 2013;63:60–5.doi:10.1016/j.ymeth.2013.06.003pmid:http://www.ncbi.nlm.nih.gov/pubmed/23806640
    OpenUrlCrossRefPubMedWeb of Science
    1. Lott SC,
    2. Schäfer RA,
    3. Mann M, et al
    . GLASSgo - Automated and Reliable Detection of sRNA Homologs From a Single Input Sequence. Front Genet 2018;9:124.doi:10.3389/fgene.2018.00124pmid:http://www.ncbi.nlm.nih.gov/pubmed/29719549
    OpenUrlCrossRefPubMed
PreviousNext
Back to top
Vol 70 Issue 8 Table of Contents
Journal of Investigative Medicine: 70 (8)
  • Table of Contents
  • Table of Contents (PDF)
  • About the Cover
  • Index by author
  • Front Matter (PDF)
Email

Thank you for your interest in spreading the word on JIM.

NOTE: We only request your email address so that the person you are recommending the page to knows that you wanted them to see it, and that it is not junk mail. We do not capture any email address.

Enter multiple addresses on separate lines or separate them with commas.
Bioinformatic and experimental methods to identify and validate bacterial RNA-human RNA interactions
(Your Name) has sent you a message from JIM
(Your Name) thought you would like to see the JIM web site.
CAPTCHA
This question is for testing whether or not you are a human visitor and to prevent automated spam submissions.
Print
Alerts
Sign In to Email Alerts with your Email Address
Citation Tools
Bioinformatic and experimental methods to identify and validate bacterial RNA-human RNA interactions
Kylie I Krohmaly, Robert J Freishtat, Andrea L Hahn
Journal of Investigative Medicine Sep 2022, jim-2022-002509; DOI: 10.1136/jim-2022-002509

Citation Manager Formats

  • BibTeX
  • Bookends
  • EasyBib
  • EndNote (tagged)
  • EndNote 8 (xml)
  • Medlars
  • Mendeley
  • Papers
  • RefWorks Tagged
  • Ref Manager
  • RIS
  • Zotero
Cite This
  • APA
  • Chicago
  • Endnote
  • MLA
Loading
Bioinformatic and experimental methods to identify and validate bacterial RNA-human RNA interactions
Kylie I Krohmaly, Robert J Freishtat, Andrea L Hahn
Journal of Investigative Medicine Sep 2022, jim-2022-002509; DOI: 10.1136/jim-2022-002509
Download PDF

Share
Bioinformatic and experimental methods to identify and validate bacterial RNA-human RNA interactions
Kylie I Krohmaly, Robert J Freishtat, Andrea L Hahn
Journal of Investigative Medicine Sep 2022, jim-2022-002509; DOI: 10.1136/jim-2022-002509
Reddit logo Twitter logo Facebook logo Mendeley logo
Respond to this article
  • Tweet Widget
  • Facebook Like
  • Google Plus One
  • Article
    • Abstract
    • Section 1: background and evidence of bacteria-to-human communication via small RNAs
    • Section 2: bioinformatic approaches for predicting bacterial sRNA-human mRNA interactions
    • Section 3: experimental approaches for identifying and validating bacterial sRNA-human mRNA interactions
    • Conclusion
    • Ethics statements
    • Footnotes
    • References
  • Figures & Data
  • eLetters
  • Info & Metrics
  • PDF

Related Articles

Cited By...

More in this TOC Section

  • Glucocorticoid therapy in respiratory illness: bench to bedside
  • Remote patient monitoring in the management of chronic obstructive pulmonary disease
  • Expanded hemodialysis: a new concept of renal replacement therapy
Show more Review

Similar Articles

 

CONTENT

  • Latest content
  • Current issue
  • Archive
  • Sign up for email alerts
  • RSS

JOURNAL

  • About the journal
  • Editorial board
  • Subscribe
  • Thank you to our reviewers
  • American Federation for Medical Research

AUTHORS

  • Information for authors
  • Submit a paper
  • Track your article
  • Open Access at BMJ

HELP

  • Contact us
  • Reprints
  • Permissions
  • Advertising
  • Feedback form

© 2023 American Federation for Medical Research