Computational prediction of microRNAs in marine bacteria of the genus Thalassospira

Thi Hoang Yen Dang; Sonika Tyagi; Glenn D’Cunha; Mrinal Bhave; Russell Crawford; Elena P. Ivanova

doi:10.1371/journal.pone.0212996

Abstract

MicroRNAs (miRNAs) are key players in regulation of gene expression at post-transcription level in eukaryotic cells. MiRNAs have been intensively studied in plants, animals and viruses. The investigations of bacterial miRNAs have gained less attention, except for the recent studies on miRNAs derived from Streptococcus mutans ATCC 25175 and Escherichia coli DH10B. In this study, high-throughput sequencing approach was employed to investigate the miRNA population in bacteria of the genus Thalassospira using both the miRDeep2 and CID-miRNA methods. A total of 984 putative miRNAs were identified in 9 species of the genus Thalassospira using both miRDeep and CID-miRNA analyses. Fifty seven conserved putative miRNAs were found in different species of the genus Thalassospira, and up to 6 miRNAs were found to be present at different locations in the T. alkalitolerans JCM 18968^T, T. lucentensis QMT2^T and T. xianhensis P-4^T. None of the putative miRNAs was found to share sequence to the reported miRNAs in E. coli DH10B and S. mutans ATCC 25175. The findings provide a comprehensive list of computationally identified miRNAs in 9 bacterial species of the genus Thalassospira and addressed the existing knowledge gap on the presence of miRNAs in the Thalassospira genomes.

Citation: Dang THY, Tyagi S, D’Cunha G, Bhave M, Crawford R, Ivanova EP (2019) Computational prediction of microRNAs in marine bacteria of the genus Thalassospira. PLoS ONE 14(3): e0212996. https://doi.org/10.1371/journal.pone.0212996

Editor: Patrick Jon Biggs, Massey University, NEW ZEALAND

Received: August 23, 2017; Accepted: February 13, 2019; Published: March 12, 2019

Copyright: © 2019 Dang et al. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.

Data Availability: All relevant data are within the paper and its Supporting Information files.

Funding: Mekong Project 1000 provided support in the form of scholarship for author THYD. Australian Genome Research Facility Ltd. provided support in the form of salaries to ST and GD. The funders had no role in the study design, data collection and analysis, decision to publish, or preparation of the manuscript. The specific roles of these authors are articulated in the ‘author contributions’ section.

Competing interests: The authors have declared that no competing interests exist.

Introduction

MicroRNA (miRNA) is a class of small, non-coding RNA molecules containing 19–22 nucleotides. Since the discovery of miRNA in Caenorhabditis elegans [1], a large number of predicted miRNA molecules has been reported in animals, plants and viruses as key players in regulation of gene expression networks [2–5]. In bacteria, the small non-coding RNAs (sRNAs) have been demonstrated to have a similar function to eukaryotic miRNA, in modulating the target mRNAs in various ways at a post-transcriptional level [6]. A number of sRNAs have been identified in bacteria, some of which were identified in marine bacteria such as Vibrio spp. and Synechococcus spp. that are functional analogues to plant miRNAs in response to environmental changes [7, 8]. The investigations of bacterial miRNAs, however, have gained little attention, except for the recent studies on miRNAs derived from Streptococcus mutans ATCC 25175 and Escherichia coli DH10B [9, 10].

The genus Thalassospira includes Gram-negative, aerobic and halophilic bacteria dwelling in a marine environment [11]. Bacteria of this genus are involved in the biodegradation of a variety of hydrocarbons [12–14]. For example, T. tepidiphila 1-1B^T and T. povalilytica Zumi 95^T have the ability to degrade polycyclic aromatic hydrocarbons and polyvinyl alcohol [12, 13]; and T. australica NP 3b2^T is able to utilise poly (ethylene) terephthalate (PET) plastic as a carbon source [14]. Currently the genus Thalassospira is comprised of 10 validly named species, of which the whole genome sequences of four species (T. australica NP 3b2^T, T. lucentensis QMT2^T, T. profundimaris WP0211^T and T. xiamenensis M-5^T) have been assembled and deposited in public databases [15–17]. The presence of miRNAs in bacteria of this genus has, however, not been demonstrated. Here, the investigation of miRNA populations was carried out using 9 validly named Thalassospira species, i.e., T. alkalitolerans JCM 18968^T, T. lucentensis QMT2^T, T. mesophila JCM 18969^T, T. povalilytica Zumi 95^T, T. profundimaris WP0211^T, T. tepidiphila 1-1B^T, T. xiamenensis M-5^T, T. xianhensis P-4^T and T. australica NP 3b2^T. Since Thalassospira lohafexi 139Z-12^T was not validly published until 2015 [11], this species was not included in this study. A deep sequencing approach has been used to identify miRNAs in bacteria in other studies [10]. Thus, the aim of this work was to identify the potential miRNAs in bacteria of the genus Thalassospira using computational approaches from small RNA sequence dataset generated by high-throughput sequencing technology.

Materials and methods

Bacterial strains and growth conditions

Nine type strains of Thalassospira species were obtained from various culture collections and used in this study. T. australica JCM 31222^T, T. alkalitolerans JCM 18968^T, T. mesophila JCM 18969^T, T. povalilytica JCM 18746^T, T. tepidiphila JCM 14578^T and T. xianhensis JCM 14850^T were obtained from RIKEN BRC-Japan Collection of Microorganisms (JCM) and T. lucentensis DSM 14000^T, T. profundimaris DSM 17430^T and T. xiamenensis DSM 17429^T were obtained from German Collection of Microorganisms and Cell Cultures (DSMZ). Bacteria were aerobically grown at their stationary phase in marine broth 2216 (BD, U.S.A.) at 22 ^oC for 48 h [18, 19]. Bacterial cells were collected by centrifugation and used for RNA extraction or stored at -80°C in marine broth 2216 (BD, USA) supplemented with 20% (v/v) glycerol.

RNA isolation and construction of small RNA libraries

Total RNAs were extracted from bacterial cells of Thalassospira species using TRIsure reagent (Bioline, Australia) according to the manufacturer’s instructions. The RNA concentration and purity were determined by measuring the absorbance at 260 nm and A260/A280 ratio using Nanodrop (ThermoFisher, Australia). The integrity of extracted RNAs was examined by a Bioanalyzer (Agilent Technology, USA). A RIN value ≥ 5 was applied for selection of good RNA quality with RIN (RNA integrity number) generated for each sample based on the ratio of ribosomal bands and the presence or absence of degradation products on the electrophoretic image. Small RNA-Seq libraries were prepared using the TruSeq Small RNA library preparation kit for Illumina (New England Biolabs, USA). Briefly, rRNA was depleted with Ribo Zero. RNAs were fragmented by heat and divalent cations and the 1^st strand cDNA was synthesised with SuperScript II Reverse Transcriptase (Invitrogen). The DNA fragments were 3’adenylated, ligated with 3’ and 5’ adapters and amplified via PCR (13 cycles). The amplified products were loaded on 6% polyacrylamide for small RNA (18–35 nucleotides) selection. The amplicons of small RNA were purified and sequenced using Illumina Hiseq2500 sequencing platform with 50 bp single end chemistry. The raw sequencing data generated from this study was submitted to the NCBI Sequence Read Archive (http://www.ncbi.nlm.nih.gov/sra) under accession number PRJNA505357.

Bioinformatics analysis of bacterial miRNAs

The raw reads (18–35 nucleotides) generated from the Illumina sequencing were first filtered by removing any contaminations including rDNA, Illumina small RNA adapter sequences or low-quality reads with < 18 nt length using Cutadapt [20]. The remaining reads after filtering were then aligned against the Thalassospira genome sequences available in the NCBI database using Bowtie (version 1.1.2) with parameters setting (-n 0 -e 80 -l 18 -a -m 5—best–strata) [21]. The whole genome sequences used as references for these 9 Thalassospira species were listed in S1 Table. The mapped sequences to genomes were further processed to remove redundancy and normalised expression was computed as CPM (Reads counts per Million) [22]. After normalisation, these mapped sequences were used to predict miRNAs. Two different miRNA prediction algorithms, miRDeep2 and CID-miRNA were used to detect miRNA in this study. The miRDeep2 software (version 2.0.0.7, http://www.mdc-berlin.de/rajewsky/miRDeep) was utilised to predict miRNAs based on an investigation of the secondary structure of the miRNA precursor sequences and integration of miRNA precursors with Dicer, providing the precursors and mature miRNA sequences [23]. CID-miRNA [24] is a method to identify miRNA precursors based on secondary structure filter and an algorithm of stochastic context free grammar (SCFG) [24]. This method only predicts the miRNA precursors. Therefore, MatureBayes (http://mirna.imbb.forth.gr/MatureBayes.html) was employed to identify mature miRNAs based on the sequence and structure of the miRNA precursors provided by CID-miRNA [25]. The small sequences mapped to the genome were processed through miRDeep2 ver2.0.0.7 as well as CID-miRNA and MatureBayes for identifying mature miRNAs with default parameters. The detected mature miRNAs were further investigated by searching any shared miRNA sequences generated in both methods using BLAST algorithm (version 4) with default setting. The potential precursor sequences predicted by both methods were also BLAST-searched to find any shared precursor sequences. Any shared miRNA sequences were then confirmed using ClustalW2 (https://www.ebi.ac.uk/Tools/msa/clustalw2/) with default setting. In order to find any conserved miRNA across the genus, the predicted mature miRNAs were aligned to those reported in Escherichia coli DH10B and Streptococcus mutans ATCC 25175 [9, 10]. The candidates were also aligned to miRNAs predicted within the genus Thalassospira to identify any conserved miRNAs across all bacteria within the genus.

Phylogenetic analysis and genome comparison

The 16S rRNA gene sequences of validly described Thalassospira species were collected from NCBI (https://www.ncbi.nlm.nih.gov/) and compared using the ClustalW2 (https://www.ebi.ac.uk/Tools/msa/clustalw2/) with default setting. Evolutionary phylogenetic trees were then generated using the neighbour-joining (NJ) algorithm [26]. Genetic distances for the NJ tree were calculated using Kimura’s two-parameter model [27], with the robustness of 1,000 replications, using MEGA 7 software [28].

Results and discussion

Generation of small RNA libraries and data evaluation

Small RNA-seq libraries were constructed for nine Thalassospira species, including T. alkalitolerans JCM 18968^T, T. lucentensis QMT2^T, T. mesophila JCM 18969^T, T. povalilytica Zumi 95^T, T. profundimaris WP0211^T, T. tepidiphila 1-1B^T, T. xiamenensis M-5^T, T. xianhensis P-4^T and T. australica NP 3b2^T. Next-generation small RNA-seq was used for detection of the presence of small RNAs. The results yielded 242,204,737 sequence reads of 18–33 nucleotides (nt) length, ranging from 22.99 to 31.43 million reads for 9 libraries (Table 1). The reads were curated for any contamination or artefacts such as rDNA or Illumina small RNA adaptor sequences using Cutadapt [20]. The cleaned reads were then aligned to the Thalassospira genome sequences available in the NCBI database using Bowtie software [21].

Download:

Table 1. Bacterial small RNA-Seq data yield from Hiseq2500 sequencing.

https://doi.org/10.1371/journal.pone.0212996.t001

Since the whole genome sequences of T. australica NP 3b2^T, T. lucentensis QMT2^T, T. profundimaris WP0211^T and T. xiamenensis M-5^T were available in the NCBI database, the small RNA sequence reads were aligned to these genome sequences. The remaining 5 Thalassospira species that did not have the whole genome sequence were aligned to these references based on their phylogenetically closest species [14] (S1 Table, S1 Fig). By aligning 34.23 million collapsed reads from nine species of the genus Thalassospira to the genome sequences available in NCBI, it was found that a proportion of cleaned reads was mappable to the reference genomes, ranging from 38.42% to 71.94% for 4 species with available genome sequences and from 2.55% to 28.16% for the remaining 5 species without whole genome sequenced (Table 1). The aligned reads were further analysed for miRNA identification.

Identification of miRNAs using the mirDeep2 method

Identification of bacterial miRNAs is in early stages in comparison to that in animals and plants. Therefore, no specific method has been developed for investigation of miRNAs in bacteria. miRDeep2 software (http://www.mdc-berlin.de/rajewsky/miRDeep) was first presented as an algorithm to predict miRNA in human and animals [23]. This software was later used to successfully identify miRNAs in plant [29]. Based on the applicability in different organisms, miRDeep2 software was, therefore, employed to detect putative bacterial miRNAs from the millions of short sequences generated from the dataset. The miRDeep2 software was designed to detect putative miRNAs based on an investigation of the secondary structure of the miRNA precursor sequences. The potential miRNA precursors were then integrated with a model for miRNA precursor processing using Dicer, releasing the mature miRNA sequences, star sequences (also called miRNA5p and miRNA3p sequences) and the loop [23]. By using miRDeep2, 86 potential miRNA precursor sequences could be identified from the dataset. The mature miRNA, star sequence and the loop were also identified within each precursor (Fig 1).

Download:

Fig 1. An example of miRDeep output for T.prof_5p_16873.

(A) Predicted secondary structure of miRNA precursor with the mature, star and loop sequences (highlighted in red, yellow and purple, respectively). The scores associated to the miRNA, the read counts for mature, star and loop sequences and the total read count are also given. (B) Sequences of the predicted miRNA precursors mapped to the mature, star and loop sequences on the genome (obs line) and the experimental sequence reported in miRBase (exp line). The frequency (reads column) and mismatches of the read with the genomic sequence (mm column) are also given with the mismatches shown in capital letters.

https://doi.org/10.1371/journal.pone.0212996.g001

The identified 86 putative mature miRNAs were located in either the 5’ arm or the 3’ arm of the precursors, with 11 sequences from T. australica NP 3b2^T, 11 from T. lucentensis QMT2^T, 14 from T. profundimaris WP0211^T, 6 from T. xiamenensis M-5^T, 11 from T. alkalitolerans JCM 18968^T, 5 from T. mesophila JCM 18969^T, 4 from T. povalilytica Zumi 95^T, 9 from T. tepidiphila 1-1B^T and 15 from T. xianhensis P-4^T (Table 2) being obtained.

Download:

Table 2. miRNAs of bacteria of the genus Thalassospira identified using miRDeep2.

https://doi.org/10.1371/journal.pone.0212996.t002

The identified putative miRNAs are 18–25 nucleotides in length (Table 2), which is in the range previously reported for animal and plant miRNAs [30]. These results are also consistent with the length distribution of miRNAs recently detected in E. coli DH10B and Streptococcus mutans ATCC 25175 [8, 9]. Based on the number of read counts, these putative miRNAs appear to have various expression levels among nine libraries, ranging from one to thousands of reads in each library, under the normal growth condition. Among the 86 detected putative miRNAs, the highest expression was found in T.prof_5p_16873, with 50145 reads, while the lowest expression was shared between T.alka_5p_4136 and T.alka_5p_1889 with only one read each being found in the libraries. Of the 86 detected putative miRNAs, thirteen were in high abundance, with over a thousand reads [10] with three miRNAs from T. australica NP 3b2^T, four from T. lucentensis QMT2^T, four from T. profundimaris WP0211^T and one each from T. tepidiphila 1-1B^T and T. xiamenensis M-5^T being identified. The high degree of expression of these putative miRNAs suggested that they may play specific roles involving in the growth and development processes of these bacteria (Fig 2).

Download:

Fig 2. Expression of miRNAs detected in bacteria of the genus Thalassospira.

Based on number of reads, the high degree of expression of miRNAs was found in T. australica NP 3b2^T, T. lucentensis QMT2^T, T. profundimaris WP0211^T, T. tepidiphila 1-1B^T and T. xiamenensis M-5^T.

https://doi.org/10.1371/journal.pone.0212996.g002

Identification of miRNAs using CID-miRNA method

The 86 putative miRNAs detected in nine Thalassospira species using miRDeep method appeared to be rather low in number compared to 400 miRNAs identified in E. coli DH10B [10] and 900 miRNAs in Streptococcus mutans ATCC 25175 [9]. Therefore, an alternative method, CID-miRNA[24] was employed in the expectation of identifying further potential miRNAs from the dataset. CID-miRNA is a web-server developed for the identification of the miRNA precursors based on the secondary structure filter and an algorithm of stochastic context free grammar (SCFG) [24]. The server was first used to predict the potential miRNA precursors in human genome [24] which was later applied in animals such as mouse, zebrafish and sea squirt [31]. Using this method, 449 potential miRNA precursors were identified from over 242 million reads of T. australica NP 3b2^T (77 potential miRNA precursors), T. lucentensis QMT2^T (161), T. profundimaris WP0211^T (78), T. xiamenensis M-5^T (16), T. alkalitolerans JCM 18968^T (21), T. mesophila JCM 18969^T (1), T. povalilytica Zumi 95^T (2), T. tepidiphila 1-1B^T (61) and T. xianhensis P-4^T (32).

In order to identify the mature miRNAs, these precursor sequences were further analysed by employing another web tool, MatureBayes (http://mirna.imbb.forth.gr/MatureBayes.html). MatureBayes is designed to incorporate to the Naïve Bayes classifier for identifying the putative mature miRNA molecules based on the sequence and structure of the miRNA precursors [25]. This tool provided putative mature miRNA positions in the 5’ and 3’ direction of the miRNA precursors, resulting in 898 putative mature miRNAs in the 9 Thalassospira species, comprised of T. australica NP 3b2^T (154 miRNAs), T. lucentensis QMT2^T (322 miRNAs), T. profundimaris WP0211^T (156 miRNAs), T. xiamenensis M-5^T (32 miRNAs), T. alkalitolerans JCM 18968^T (42 miRNAs), T. mesophila JCM 18969^T (2 miRNAs), T. povalilytica Zumi 95^T (4 miRNAs), T. tepidiphila 1-1B^T (122 miRNAs) and T. xianhensis P-4^T (64 miRNAs) (S2 Table). These identified putative miRNAs were used for further analysis of miRNA conservation.

Comparison of bacterial miRNAs identified using miRDeep2 and CID-miRNA

In order to determine whether any miRNAs were detected by both methods, miRNA data generated from miRDeep were BLAST-searched to those produced by CID-miRNA. The results showed that there were no shared putative mature miRNA sequences. miRDeep2 predicts putative miRNAs based on distribution of minimum free energy (MFE) and stability of secondary structure established for nematode, C. elegans, while CID-miRNA utilised secondary structure filter and stochastic context free grammar trained on human miRNAs [23, 24]. The differences of parameter setting may cause no overlap of putative mature miRNA sequences detected in both methods. Thakur et al [32] also pointed out that application of new parameters improved the accuracy of plant miRNA prediction using miRDeep in compared to the default setting for animals.

Study on miRNA biogenesis showed that one or more mature miRNAs can be produced from one pre-miRNA molecule [33]. These identified putative miRNAs, therefore, were aligned to the potential precursor miRNA sequences generated by both methods. BLAST-search between putative mature miRNA sequences identified from miRDeep2 and the potential precursor miRNAs detected by CID-miRNA was carried out and vice versa. 5 putative mature miRNA sequences identified by CID-miRNA (T.luce_5p_228121, T.luce_3p_228121, T.xian_5p_2738, T.xian_5p_2740 and T.xian_3p_2740) were found to locate in 3 potential precursor miRNAs producing T.luce_5p_11956, T.xian_5p_4477 and T.xian_3p_4479 from miRDeep2. The precursor miRNA sequences producing 5 putative mature miRNAs from CID-miRNA were then aligned to 3 potential precursor miRNAs from miRDeep2 using ClustalW2. The result showed the overlap of 3 precursor miRNAs predicted in both methods in which the putative mature miRNAs identified by each method located at different positions within the precursor sequences (Fig 3).

Download:

Fig 3. Sharing sequences of potential miRNA precursors identified in both miRDeep2 and CID-miRNA methods.

(A) T.luce_228121 precursor from CID-miRNA and T.luce_5p_11956 precursor from miRDeep2. (B) T.xian_2738 precursor from CID-miRNA and T.xian_5p_4477 precursor from miRDeep2. (C) T.xian_2740 precursor from CID-miRNA and T.xian_3p_4479 precursor from miRDeep2. The asterisk (*) indicate that the nucleotides are identical to the top sequence. The putative mature miRNAs are highlighted in grey. (CID-miRNA) indicated the potential precursor detected by CID-miRNA method. (miRDeep2) indicated the potential precursor detected by miRDeep2 method.

https://doi.org/10.1371/journal.pone.0212996.g003

The shared sequences from the same Thalassospira species proved the probability of putative miRNAs finding in both methods. Thus, a total of 984 putative miRNA candidates were preliminarily obtained from both methods with default setting together for 9 species (Table 3). Further work, such as real time PCR or Northern blot analysis, is needed to verify these putative miRNAs.

Download:

Table 3. Identification of miRNAs using miRDeep2 and CID-miRNA methods.

https://doi.org/10.1371/journal.pone.0212996.t003

Conservation of miRNAs in the genus Thalassospira and miRNAs previously described in Streptococcus mutans ATCC 25175 and E. coli DH10B

Previous studies showed the conservation of some miRNAs across animals or the plant kingdom [34, 35]; however, novel sequences can only be found in a particular species. The conserved and novel miRNAs in bacteria still remain largely unknown, in comparison to intensive studies of miRNAs in eukaryotic organisms and viruses. In this study, 984 putative miRNA candidates of the genus Thalassospira were BLAST-searched against those reported for E. coli DH10B (400 miRNAs) and S. mutans ATCC 25175 (900 miRNAs) [9, 10] in order to identify any conserved miRNA. However, no conserved sequences could be found without any mismatch or with three and fewer nucleotide substitution. E. coli is an enteric bacterium commonly found in low intestine of humans and animals [36], S. mutans is an oral pathogen that causes human dental caries [37], while bacteria of the Thalassospira genus are environmental bacteria [14]. The differences in their characteristics may influence the low degree of conservation of miRNAs among these bacteria. The results obtained in this study is in agreement with previous studies that reported the lack of significant sequence similarity in non-coding RNA homologues of different bacterial species [38, 39].

In order to identify any conserved miRNAs among species of the genus Thalassospira, the putative miRNA sequences obtained were compared. It appears that the nine species shared 57 common putative miRNAs. Among these, T. profundimaris WP0211^T and T. tepidiphila 1-1B^T had the highest number of conserved miRNAs, with 45 sequences presenting in both species. Five putative miRNAs in T. alkalitolerans JCM 18968^T were also found in the T. lucentensis QMT2^T, while T. xianhensis P-4^T and T. xiamenensis M-5^T shared 4 conserved sequences. One sequence was shared between T. australica NP 3b2^T and T. tepidiphila 1-1B^T, T. lucentensis QMT2^T and T. tepidiphila 1-1B^T and T. lucentensis QMT2^T and T. mesophila JCM 18969^T (Table 4). As seen from the data, T. lucentensis QMT2^T has identical miRNA sequences as T. alkalitolerans JCM 18968^T (5 sequences), T. tepidiphila 1-1B^T (1) and T. mesophila JCM 18969^T (1), while T. tepidiphila 1-1B^T also shared sequences with T. profundimaris WP0211^T (45) and T. australica NP 3b2^T (1).

Download:

Table 4. Conserved miRNAs in different species of bacteria of the genus Thalassospira.

https://doi.org/10.1371/journal.pone.0212996.t004

A comparative analysis of 16S rRNA sequence similarities revealed that the highest number of common miRNAs shared by the phylogenetically close species; e.g., T. tepidiphila 1-1B^T were found to be phylogenetically closely related to T. profundimaris WP0211^T (99.3% sequence similarity); T. xianhensis P-4^T and T. xiamenensis M-5^T also shared 99.3% 16S rRNA similarity, while T. mesophila JCM 18969^T and T. alkalitolerans JCM 18968^T shared 95.3% and 94.9% of 16S rRNA sequence similarity, respectively, with T. lucentensis QMT2^T [14]. Conserved miRNA sequences found in these species may indicate a close genetic relationship and that these miRNAs have a similar role in the regulation of the growth and development of bacteria. It is interesting to note that two species, T. australica NP 3b2^T and T. tepidiphila 1-1B^T, which are capable of hydrocarbon degradation [12, 14], shared one conserved miRNAs and the potential targets of this putative miRNA needs to be investigated for any role related to this process. In addition, T. alkalitolerans JCM 18968^T, T. lucentensis QMT2^T and T. xianhensis P-4^T were also found to have the same sequence presented in different genomic locations (Table 5). These miRNAs can have an influence on their expression and function at different genomic locations [40]. It will be of great interest to identify the target mRNA of these miRNAs and investigate their roles and the mechanisms of gene regulation in the physiology of these unique environmental bacteria.

Download:

Table 5. miRNAs present in different locations of the same Thalassospira species.

https://doi.org/10.1371/journal.pone.0212996.t005

Conclusions

Over 242 million reads of 18 to 33 nucleotides length were generated from nine bacterial species using high-throughput sequencing technology. Using miRDeep2 and CID-miRNA analyses, a total of 984 putative miRNAs were eventually identified, with typical miRNA length of 19–25 nucleotides. Compared to other species, these detected putative miRNAs were not conserved to those reported in E. coli DH10B and S. mutans ATCC 25175. This study presents the first comprehensive list of computationally identified miRNAs in 9 bacterial species of the genus Thalassospira without experimental verification. The further work, however, is needed to validate these candidates experimentally. In addition, further identification of miRNA targets will provide insights into the fundamental functions of miRNAs in the physiology of these bacteria.

Supporting information

S1 Fig. Phylogenetic tree of valid species of the genus Thalassospira.

https://doi.org/10.1371/journal.pone.0212996.s001

(DOCX)

S1 Table. Whole genome sequences of bacteria of the genus Thalassospira available in GenBank.

https://doi.org/10.1371/journal.pone.0212996.s002

(DOCX)

S2 Table. miRNAs retrieved from bacteria of the genus Thalassospira as identified by CID-miRNA analysis.

https://doi.org/10.1371/journal.pone.0212996.s003

(DOCX)

References

1. Lee RC, Feinbaum RL, Ambros V. The C. elegans heterochronic gene lin-4 encodes small RNAs with antisense complementarity to lin-14. Cell. 1993;75(5):843–54. pmid:8252621
- View Article
- PubMed/NCBI
- Google Scholar
2. Cullen BR. Viruses and microRNAs. Nature Genet. 2006;38:S25–S30. pmid:16736021
- View Article
- PubMed/NCBI
- Google Scholar
3. Takeda A, Watanabe Y. Small RNA world in plants. Tanpakushitsu Kakusan Koso. 2006;51(16 Suppl):2463–70. pmid:17471964
- View Article
- PubMed/NCBI
- Google Scholar
4. Gottesman S. Micros for microbes: Non-coding regulatory RNAs in bacteria. Trends Genet. 2005;21(7):399–404. pmid:15913835
- View Article
- PubMed/NCBI
- Google Scholar
5. Plasterk RHA. Micro RNAs in animal development. Cell. 2006;124(5):877–81. pmid:16530032
- View Article
- PubMed/NCBI
- Google Scholar
6. Gottesman S, Storz G. Bacterial small RNA regulators: versatile roles and rapidly evolving variations. Cold Spring Harb Perspect Biol. 2011;3(12).
- View Article
- Google Scholar
7. Gierga G, Voss B, Hess WR. Non-coding RNAs in marine Synechococcus and their regulation under environmentally relevant stress conditions. ISME Journal. 2012;6(8):1544–57. pmid:22258101
- View Article
- PubMed/NCBI
- Google Scholar
8. Nguyen A, Jacq A. Small RNAs in the Vibrionaceae: An ocean still to be explored. Wiley Interdiscip Rev RNA. 2014;5(3):381–92. pmid:24458378
- View Article
- PubMed/NCBI
- Google Scholar
9. Lee HJ, Hong SH. Analysis of microRNA-size, small RNAs in Streptococcus mutans by deep sequencing. FEMS Microbiol Lett. 2012;326(2):131–6. pmid:22092283
- View Article
- PubMed/NCBI
- Google Scholar
10. Kang SM, Choi JW, Lee Y, Hong SH, Lee HJ. Identification of microRNA-Size, Small RNAs in Escherichia coli. Curr Microbiol. 2013;67(5):609–13. pmid:23783561
- View Article
- PubMed/NCBI
- Google Scholar
11. Shivaji S, Sathyanarayana Reddy G, Sundareswaran VR, Thomas C. Description of Thalassospira lohafexi sp. nov., isolated from Southern Ocean, Antarctica. Arch Microbiol. 2015;197(5):627–37. pmid:25702315
- View Article
- PubMed/NCBI
- Google Scholar
12. Kodama Y, Stiknowati LI, Ueki A, Ueki K, Watanabe K. Thalassospira tepidiphila sp. nov., a polycyclic aromatic hydrocarbon-degrading bacterium isolated from seawater. Int J Syst Evol Microbiol. 2008;58(3):711–5.
- View Article
- Google Scholar
13. Nogi Y, Yoshizumi M, Miyazaki M. Thalassospira povalilytica sp. nov., a polyvinyl-alcohol-degrading marine bacterium. Int J Syst Evol Microbiol. 2014;64(PART 4):1149–53.
- View Article
- Google Scholar
14. Ivanova EP, López-Pérez M, Webb HK, Ng HJ, Dang THY, Zhukova NV, et al. Thalassospira australica sp. nov. isolated from sea water. Antonie van Leeuwenhoek. 2016;109(8):1091–100. pmid:27180097
- View Article
- PubMed/NCBI
- Google Scholar
15. Lai Q, Shao Z. Genome sequence of Thalassospira xiamenensis type strain M-5. J Bacteriol. 2012;194(24):6957–. PMC3510571. pmid:23209216
- View Article
- PubMed/NCBI
- Google Scholar
16. Lai Q, Shao Z. Genome sequence of Thalassospira profundimaris type strain WP0211. J Bacteriol. 2012;194(24):6956–. PMC3510614. pmid:23209215
- View Article
- PubMed/NCBI
- Google Scholar
17. López-Pérez M, Rodriguez-Valera F, Webb HK, Crawford RJ, Ivanova EP. Genome Sequence of “Thalassospira australica” NP3b2T Isolated from St. Kilda Beach, Tasman Sea. Genome Announc. 2014;2(6):e01139–14. PMC4241657. pmid:25395631
- View Article
- PubMed/NCBI
- Google Scholar
18. Hütz A, Schubert K, Overmann J. Thalassospira sp. isolated from the oligotrophic eastern mediterranean sea exhibits chemotaxis toward inorganic phosphate during starvation. Appl Environ Microbiol. 2011;77(13):4412–21. pmid:21602377
- View Article
- PubMed/NCBI
- Google Scholar
19. Guo J, Zhang L, Lu X, Zeng S, Zhang Y, Xu H, et al. Medium optimization and fermentation kinetics for κ-carrageenase production by Thalassospira sp. Fjfst-332. Molecules. 2016;21(11):1479. pmid:27827964
- View Article
- PubMed/NCBI
- Google Scholar
20. Martin M. Cutadapt removes adapter sequences from high-throughput sequencing reads. EMBnetjournal. 2011;17(1):1–3. pp. 10–12.
- View Article
- Google Scholar
21. Langmead B, Trapnell C, Pop M, Salzberg SL. Ultrafast and memory-efficient alignment of short DNA sequences to the human genome. Genome Biol. 2009;10(3):1–10. pmid:19261174
- View Article
- PubMed/NCBI
- Google Scholar
22. Ozhuner E, Eldem V, Ipek A, Okay S, Sakcali S, Zhang B, et al. Boron stress responsive micrornas and their targets in barley. PLoS ONE. 2013;8(3). PubMed Central PMCID: PMCPMC3608689. pmid:23555702
- View Article
- PubMed/NCBI
- Google Scholar
23. Friedlander MR, Chen W, Adamidi C, Maaskola J, Einspanier R, Knespel S, et al. Discovering microRNAs from deep sequencing data using miRDeep. Nat Biotech. 2008;26(4):407–15. pmid:18392026
- View Article
- PubMed/NCBI
- Google Scholar
24. Tyagi S, Vaz C, Gupta V, Bhatia R, Maheshwari S, Srinivasan A, et al. CID-miRNA: A web server for prediction of novel miRNA precursors in human genome. Biochem Biophys Res Commun. 2008;372(4):831–4. pmid:18522801
- View Article
- PubMed/NCBI
- Google Scholar
25. Gkirtzou K, Tsamardinos I, Tsakalides P, Poirazi P. MatureBayes: a probabilistic algorithm for identifying the mature miRNA within novel precursors. PLoS ONE. 2010;5(8):1–14. pmid:20700506
- View Article
- PubMed/NCBI
- Google Scholar
26. Saitou N, Nei M. The neighbor-joining method: a new method for reconstructing phylogenetic trees. Mol Biol Evol. 1987;4(4):406–25. pmid:3447015
- View Article
- PubMed/NCBI
- Google Scholar
27. Kimura M. A simple method for estimating evolutionary rates of base substitutions through comparative studies of nucleotide sequences. J Mol Evol. 1980;16(2):111–20. pmid:7463489
- View Article
- PubMed/NCBI
- Google Scholar
28. Hall BG. Building phylogenetic trees from molecular data with MEGA. Mol Biol Evol. 2013;30(5):1229–35. pmid:23486614
- View Article
- PubMed/NCBI
- Google Scholar
29. Hackenberg M, Huang PJ, Huang CY, Shi BJ, Gustafson P, Langridge P. A Comprehensive expression profile of micrornas and other classes of non-coding small RNAs in barley under phosphorous-deficient and-sufficient conditions. DNA Research. 2013;20(2):109–25. pmid:23266877
- View Article
- PubMed/NCBI
- Google Scholar
30. Millar AA, Waterhouse PM. Plant and animal microRNAs: Similarities and differences. Funct Integr Genomics. 2005;5(3):129–35. pmid:15875226
- View Article
- PubMed/NCBI
- Google Scholar
31. Tempel S, Tahi F. A fast ab-initio method for predicting miRNA precursors in genomes. Nucleic Acids Research. 2012;40(11):1–9. pmid:22362754
- View Article
- PubMed/NCBI
- Google Scholar
32. Thakur V, Wanchana S, Xu M, Bruskiewich R, Quick WP, Mosig A, et al. Characterization of statistical features for plant microRNA prediction. BMC Genomics. 2011;12:108. pmid:21324149
- View Article
- PubMed/NCBI
- Google Scholar
33. Bologna NG, Schapire AL, Palatnik JF. Processing of plant microRNA precursors. Briefings in Functional Genomics. 2013;12(1):37–45. pmid:23148323
- View Article
- PubMed/NCBI
- Google Scholar
34. Wheeler BM, Heimberg AM, Moy VN, Sperling EA, Holstein TW, Heber S, et al. The deep evolution of metazoan microRNAs. Evol Dev. 2009;11(1):50–68. pmid:19196333
- View Article
- PubMed/NCBI
- Google Scholar
35. Schreiber AW, Shi B-J, Huang C-Y, Langridge P, Baumann U. Discovery of barley miRNAs through deep sequencing of short reads. BMC Genomics. 2011;12:1–21. WOS:000288489600001. pmid:21352554
- View Article
- PubMed/NCBI
- Google Scholar
36. Kaper JB, Nataro JP, Mobley HLT. Pathogenic Escherichia coli. Nature Rev Microbiol. 2004;2(2):123–40. pmid:15040260
- View Article
- PubMed/NCBI
- Google Scholar
37. Ajdić D, McShan WM, McLaughlin RE, Savić G, Chang J, Carson MB, et al. Genome sequence of Streptococcus mutans UA159, a cariogenic dental pathogen. Proc Natl Acad Sci USA. 2002;99(22):14434–9. pmid:12397186
- View Article
- PubMed/NCBI
- Google Scholar
38. Zhao T, Zhang R, Wang M. Prediction of candidate small non-coding RNAs in by Agrobacterium by computational analysis. Journal of Biomedical Research. 2010;24(1):33–42. pmid:23554609
- View Article
- PubMed/NCBI
- Google Scholar
39. Livny J, Waldor MK. Identification of small RNAs in diverse bacterial species. Current Opinion in Microbiology. 2007;10(2):96–101. pmid:17383222
- View Article
- PubMed/NCBI
- Google Scholar
40. Paczynska P, Grzemski A, Szydlowski M. Distribution of miRNA genes in the pig genome. BMC Genetics. 2015;16(1):1–12. PMC4318388. pmid:25632794
- View Article
- PubMed/NCBI
- Google Scholar

[ref1] 1. Lee RC, Feinbaum RL, Ambros V. The C. elegans heterochronic gene lin-4 encodes small RNAs with antisense complementarity to lin-14. Cell. 1993;75(5):843–54. pmid:8252621
View Article
PubMed/NCBI
Google Scholar

[2] View Article

[3] PubMed/NCBI

[4] Google Scholar

[ref2] 2. Cullen BR. Viruses and microRNAs. Nature Genet. 2006;38:S25–S30. pmid:16736021
View Article
PubMed/NCBI
Google Scholar

[6] View Article

[7] PubMed/NCBI

[8] Google Scholar

[ref3] 3. Takeda A, Watanabe Y. Small RNA world in plants. Tanpakushitsu Kakusan Koso. 2006;51(16 Suppl):2463–70. pmid:17471964
View Article
PubMed/NCBI
Google Scholar

[10] View Article

[11] PubMed/NCBI

[12] Google Scholar

[ref4] 4. Gottesman S. Micros for microbes: Non-coding regulatory RNAs in bacteria. Trends Genet. 2005;21(7):399–404. pmid:15913835
View Article
PubMed/NCBI
Google Scholar

[14] View Article

[15] PubMed/NCBI

[16] Google Scholar

[ref5] 5. Plasterk RHA. Micro RNAs in animal development. Cell. 2006;124(5):877–81. pmid:16530032
View Article
PubMed/NCBI
Google Scholar

[18] View Article

[19] PubMed/NCBI

[20] Google Scholar

[ref6] 6. Gottesman S, Storz G. Bacterial small RNA regulators: versatile roles and rapidly evolving variations. Cold Spring Harb Perspect Biol. 2011;3(12).
View Article
Google Scholar

[22] View Article

[23] Google Scholar

[ref7] 7. Gierga G, Voss B, Hess WR. Non-coding RNAs in marine Synechococcus and their regulation under environmentally relevant stress conditions. ISME Journal. 2012;6(8):1544–57. pmid:22258101
View Article
PubMed/NCBI
Google Scholar

[25] View Article

[26] PubMed/NCBI

[27] Google Scholar

[ref8] 8. Nguyen A, Jacq A. Small RNAs in the Vibrionaceae: An ocean still to be explored. Wiley Interdiscip Rev RNA. 2014;5(3):381–92. pmid:24458378
View Article
PubMed/NCBI
Google Scholar

[29] View Article

[30] PubMed/NCBI

[31] Google Scholar

[ref9] 9. Lee HJ, Hong SH. Analysis of microRNA-size, small RNAs in Streptococcus mutans by deep sequencing. FEMS Microbiol Lett. 2012;326(2):131–6. pmid:22092283
View Article
PubMed/NCBI
Google Scholar

[33] View Article

[34] PubMed/NCBI

[35] Google Scholar

[ref10] 10. Kang SM, Choi JW, Lee Y, Hong SH, Lee HJ. Identification of microRNA-Size, Small RNAs in Escherichia coli. Curr Microbiol. 2013;67(5):609–13. pmid:23783561
View Article
PubMed/NCBI
Google Scholar

[37] View Article

[38] PubMed/NCBI

[39] Google Scholar

[ref11] 11. Shivaji S, Sathyanarayana Reddy G, Sundareswaran VR, Thomas C. Description of Thalassospira lohafexi sp. nov., isolated from Southern Ocean, Antarctica. Arch Microbiol. 2015;197(5):627–37. pmid:25702315
View Article
PubMed/NCBI
Google Scholar

[41] View Article

[42] PubMed/NCBI

[43] Google Scholar

[ref12] 12. Kodama Y, Stiknowati LI, Ueki A, Ueki K, Watanabe K. Thalassospira tepidiphila sp. nov., a polycyclic aromatic hydrocarbon-degrading bacterium isolated from seawater. Int J Syst Evol Microbiol. 2008;58(3):711–5.
View Article
Google Scholar

[45] View Article

[46] Google Scholar

[ref13] 13. Nogi Y, Yoshizumi M, Miyazaki M. Thalassospira povalilytica sp. nov., a polyvinyl-alcohol-degrading marine bacterium. Int J Syst Evol Microbiol. 2014;64(PART 4):1149–53.
View Article
Google Scholar

[48] View Article

[49] Google Scholar

[ref14] 14. Ivanova EP, López-Pérez M, Webb HK, Ng HJ, Dang THY, Zhukova NV, et al. Thalassospira australica sp. nov. isolated from sea water. Antonie van Leeuwenhoek. 2016;109(8):1091–100. pmid:27180097
View Article
PubMed/NCBI
Google Scholar

[51] View Article

[52] PubMed/NCBI

[53] Google Scholar

[ref15] 15. Lai Q, Shao Z. Genome sequence of Thalassospira xiamenensis type strain M-5. J Bacteriol. 2012;194(24):6957–. PMC3510571. pmid:23209216
View Article
PubMed/NCBI
Google Scholar

[55] View Article

[56] PubMed/NCBI

[57] Google Scholar

[ref16] 16. Lai Q, Shao Z. Genome sequence of Thalassospira profundimaris type strain WP0211. J Bacteriol. 2012;194(24):6956–. PMC3510614. pmid:23209215
View Article
PubMed/NCBI
Google Scholar

[59] View Article

[60] PubMed/NCBI

[61] Google Scholar

[ref17] 17. López-Pérez M, Rodriguez-Valera F, Webb HK, Crawford RJ, Ivanova EP. Genome Sequence of “Thalassospira australica” NP3b2T Isolated from St. Kilda Beach, Tasman Sea. Genome Announc. 2014;2(6):e01139–14. PMC4241657. pmid:25395631
View Article
PubMed/NCBI
Google Scholar

[63] View Article

[64] PubMed/NCBI

[65] Google Scholar

[ref18] 18. Hütz A, Schubert K, Overmann J. Thalassospira sp. isolated from the oligotrophic eastern mediterranean sea exhibits chemotaxis toward inorganic phosphate during starvation. Appl Environ Microbiol. 2011;77(13):4412–21. pmid:21602377
View Article
PubMed/NCBI
Google Scholar

[67] View Article

[68] PubMed/NCBI

[69] Google Scholar

[ref19] 19. Guo J, Zhang L, Lu X, Zeng S, Zhang Y, Xu H, et al. Medium optimization and fermentation kinetics for κ-carrageenase production by Thalassospira sp. Fjfst-332. Molecules. 2016;21(11):1479. pmid:27827964
View Article
PubMed/NCBI
Google Scholar

[71] View Article

[72] PubMed/NCBI

[73] Google Scholar

[ref20] 20. Martin M. Cutadapt removes adapter sequences from high-throughput sequencing reads. EMBnetjournal. 2011;17(1):1–3. pp. 10–12.
View Article
Google Scholar

[75] View Article

[76] Google Scholar

[ref21] 21. Langmead B, Trapnell C, Pop M, Salzberg SL. Ultrafast and memory-efficient alignment of short DNA sequences to the human genome. Genome Biol. 2009;10(3):1–10. pmid:19261174
View Article
PubMed/NCBI
Google Scholar

[78] View Article

[79] PubMed/NCBI

[80] Google Scholar

[ref22] 22. Ozhuner E, Eldem V, Ipek A, Okay S, Sakcali S, Zhang B, et al. Boron stress responsive micrornas and their targets in barley. PLoS ONE. 2013;8(3). PubMed Central PMCID: PMCPMC3608689. pmid:23555702
View Article
PubMed/NCBI
Google Scholar

[82] View Article

[83] PubMed/NCBI

[84] Google Scholar

[ref23] 23. Friedlander MR, Chen W, Adamidi C, Maaskola J, Einspanier R, Knespel S, et al. Discovering microRNAs from deep sequencing data using miRDeep. Nat Biotech. 2008;26(4):407–15. pmid:18392026
View Article
PubMed/NCBI
Google Scholar

[86] View Article

[87] PubMed/NCBI

[88] Google Scholar

[ref24] 24. Tyagi S, Vaz C, Gupta V, Bhatia R, Maheshwari S, Srinivasan A, et al. CID-miRNA: A web server for prediction of novel miRNA precursors in human genome. Biochem Biophys Res Commun. 2008;372(4):831–4. pmid:18522801
View Article
PubMed/NCBI
Google Scholar

[90] View Article

[91] PubMed/NCBI

[92] Google Scholar

[ref25] 25. Gkirtzou K, Tsamardinos I, Tsakalides P, Poirazi P. MatureBayes: a probabilistic algorithm for identifying the mature miRNA within novel precursors. PLoS ONE. 2010;5(8):1–14. pmid:20700506
View Article
PubMed/NCBI
Google Scholar

[94] View Article

[95] PubMed/NCBI

[96] Google Scholar

[ref26] 26. Saitou N, Nei M. The neighbor-joining method: a new method for reconstructing phylogenetic trees. Mol Biol Evol. 1987;4(4):406–25. pmid:3447015
View Article
PubMed/NCBI
Google Scholar

[98] View Article

[99] PubMed/NCBI

[100] Google Scholar

[ref27] 27. Kimura M. A simple method for estimating evolutionary rates of base substitutions through comparative studies of nucleotide sequences. J Mol Evol. 1980;16(2):111–20. pmid:7463489
View Article
PubMed/NCBI
Google Scholar

[102] View Article

[103] PubMed/NCBI

[104] Google Scholar

[ref28] 28. Hall BG. Building phylogenetic trees from molecular data with MEGA. Mol Biol Evol. 2013;30(5):1229–35. pmid:23486614
View Article
PubMed/NCBI
Google Scholar

[106] View Article

[107] PubMed/NCBI

[108] Google Scholar

[ref29] 29. Hackenberg M, Huang PJ, Huang CY, Shi BJ, Gustafson P, Langridge P. A Comprehensive expression profile of micrornas and other classes of non-coding small RNAs in barley under phosphorous-deficient and-sufficient conditions. DNA Research. 2013;20(2):109–25. pmid:23266877
View Article
PubMed/NCBI
Google Scholar

[110] View Article

[111] PubMed/NCBI

[112] Google Scholar

[ref30] 30. Millar AA, Waterhouse PM. Plant and animal microRNAs: Similarities and differences. Funct Integr Genomics. 2005;5(3):129–35. pmid:15875226
View Article
PubMed/NCBI
Google Scholar

[114] View Article

[115] PubMed/NCBI

[116] Google Scholar

[ref31] 31. Tempel S, Tahi F. A fast ab-initio method for predicting miRNA precursors in genomes. Nucleic Acids Research. 2012;40(11):1–9. pmid:22362754
View Article
PubMed/NCBI
Google Scholar

[118] View Article

[119] PubMed/NCBI

[120] Google Scholar

[ref32] 32. Thakur V, Wanchana S, Xu M, Bruskiewich R, Quick WP, Mosig A, et al. Characterization of statistical features for plant microRNA prediction. BMC Genomics. 2011;12:108. pmid:21324149
View Article
PubMed/NCBI
Google Scholar

[122] View Article

[123] PubMed/NCBI

[124] Google Scholar

[ref33] 33. Bologna NG, Schapire AL, Palatnik JF. Processing of plant microRNA precursors. Briefings in Functional Genomics. 2013;12(1):37–45. pmid:23148323
View Article
PubMed/NCBI
Google Scholar

[126] View Article

[127] PubMed/NCBI

[128] Google Scholar

[ref34] 34. Wheeler BM, Heimberg AM, Moy VN, Sperling EA, Holstein TW, Heber S, et al. The deep evolution of metazoan microRNAs. Evol Dev. 2009;11(1):50–68. pmid:19196333
View Article
PubMed/NCBI
Google Scholar

[130] View Article

[131] PubMed/NCBI

[132] Google Scholar

[ref35] 35. Schreiber AW, Shi B-J, Huang C-Y, Langridge P, Baumann U. Discovery of barley miRNAs through deep sequencing of short reads. BMC Genomics. 2011;12:1–21. WOS:000288489600001. pmid:21352554
View Article
PubMed/NCBI
Google Scholar

[134] View Article

[135] PubMed/NCBI

[136] Google Scholar

[ref36] 36. Kaper JB, Nataro JP, Mobley HLT. Pathogenic Escherichia coli. Nature Rev Microbiol. 2004;2(2):123–40. pmid:15040260
View Article
PubMed/NCBI
Google Scholar

[138] View Article

[139] PubMed/NCBI

[140] Google Scholar

[ref37] 37. Ajdić D, McShan WM, McLaughlin RE, Savić G, Chang J, Carson MB, et al. Genome sequence of Streptococcus mutans UA159, a cariogenic dental pathogen. Proc Natl Acad Sci USA. 2002;99(22):14434–9. pmid:12397186
View Article
PubMed/NCBI
Google Scholar

[142] View Article

[143] PubMed/NCBI

[144] Google Scholar

[ref38] 38. Zhao T, Zhang R, Wang M. Prediction of candidate small non-coding RNAs in by Agrobacterium by computational analysis. Journal of Biomedical Research. 2010;24(1):33–42. pmid:23554609
View Article
PubMed/NCBI
Google Scholar

[146] View Article

[147] PubMed/NCBI

[148] Google Scholar

[ref39] 39. Livny J, Waldor MK. Identification of small RNAs in diverse bacterial species. Current Opinion in Microbiology. 2007;10(2):96–101. pmid:17383222
View Article
PubMed/NCBI
Google Scholar

[150] View Article

[151] PubMed/NCBI

[152] Google Scholar

[ref40] 40. Paczynska P, Grzemski A, Szydlowski M. Distribution of miRNA genes in the pig genome. BMC Genetics. 2015;16(1):1–12. PMC4318388. pmid:25632794
View Article
PubMed/NCBI
Google Scholar

[154] View Article

[155] PubMed/NCBI

[156] Google Scholar

Figures

Abstract

Introduction

Materials and methods

Bacterial strains and growth conditions

RNA isolation and construction of small RNA libraries

Bioinformatics analysis of bacterial miRNAs

Phylogenetic analysis and genome comparison

Results and discussion

Generation of small RNA libraries and data evaluation

Identification of miRNAs using the mirDeep2 method

Identification of miRNAs using CID-miRNA method

Comparison of bacterial miRNAs identified using miRDeep2 and CID-miRNA

Conservation of miRNAs in the genus Thalassospira and miRNAs previously described in Streptococcus mutans ATCC 25175 and E. coli DH10B

Conclusions

Supporting information

S1 Fig. Phylogenetic tree of valid species of the genus Thalassospira.

S1 Table. Whole genome sequences of bacteria of the genus Thalassospira available in GenBank.

S2 Table. miRNAs retrieved from bacteria of the genus Thalassospira as identified by CID-miRNA analysis.

References