Next Article in Journal
Genome-Wide Identification, Characterization, and Expression Profiling of Glutathione S-Transferase (GST) Family in Pumpkin Reveals Likely Role in Cold-Stress Tolerance
Next Article in Special Issue
Biofilm Formation by the Acidophile Bacterium Acidithiobacillus thiooxidans Involves c-di-GMP Pathway and Pel exopolysaccharide
Previous Article in Journal
Re-Arrangements in the Cytoplasmic Distribution of Small RNAs Following the Maternal-to-Zygotic Transition in Drosophila Embryos
Previous Article in Special Issue
Development of an Effective 6-Methylpurine Counterselection Marker for Genetic Manipulation in Thermococcus barophilus
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Small RNAs of Haloferax mediterranei: Identification and Potential Involvement in Nitrogen Metabolism

by
Gloria Payá
1,
Vanesa Bautista
1,
Mónica Camacho
1,
Natalia Castejón-Fernández
2,
Luís A. Alcaraz
1,2,
María-José Bonete
1,* and
Julia Esclapez
1,*
1
Agrochemistry and Biochemistry Department, Biochemistry and Molecular Biology Division, Faculty of Science, University of Alicante, Ap 99, E-03080 Alicante, Spain
2
Bioarray, S.L., 03202 Alicante, Spain
*
Authors to whom correspondence should be addressed.
Genes 2018, 9(2), 83; https://doi.org/10.3390/genes9020083
Submission received: 30 December 2017 / Revised: 26 January 2018 / Accepted: 2 February 2018 / Published: 10 February 2018
(This article belongs to the Special Issue Genetics and Genomics of Extremophiles)

Abstract

:
Small RNAs have been studied in detail in domains Bacteria and Eukarya but, in the case of the domain Archaea, the knowledge is scarce and the physiological function of these small RNAs (sRNAs) is still uncertain. To extend the knowledge of sRNAs in the domain Archaea and their possible role in the regulation of the nitrogen assimilation metabolism in haloarchaea, Haloferax mediterranei has been used as a model microorganism. The bioinformatic approach has allowed for the prediction of 295 putative sRNAs genes in the genome of H. mediterranei, 88 of which have been verified by means of RNA-Sequencing (RNA-Seq). The secondary structure of these sRNAs and their possible targets have been identified. Curiously, some of them present as possible target genes relating to nitrogen assimilation, such as glutamate dehydrogenase and the nitrogen regulatory PII protein. Analysis of RNA-Seq data has also revealed differences in the expression pattern of 16 sRNAs according to the nitrogen source. Consequently, RNomic and bioinformatic approaches used in this work have allowed for the identification of new sRNAs in H. mediterranei, some of which show different expression patterns depending on the nitrogen source. This suggests that these sRNAs could be involved in the regulation of nitrogen assimilation and can constitute an important gene regulatory network.

1. Introduction

Small RNAs (sRNAs) play an essential role in the post-transcriptional regulation of many cellular processes in all domains of life, i.e., Eukarya, Bacteria, and Archaea. In eukaryotes, there are different classes of these, the best studied being microRNA (miRNAs), small interference RNA (siRNAs), and piwi-interacting RNAs (piRNAs). These are approximately 20–30 nucleotides (nt) in length and are involved in development, cellular activities, and different physiology processes [1,2,3]. In Archaea and Bacteria, sRNAs are much longer than eukaryotic small non-coding RNAs, which range in length from 50 to 500 nt [4,5,6,7]. Different mechanisms of action of sRNAs have been described, most of which affect the translation of the target messenger RNA and/or its stability [4]. Hence, they seem to be involved in the post-transcriptional regulation of metabolism, stress response, virulence processes, and so on. Depending on the location of their targets, sRNAs are classified into two groups: trans-encoded and cis-encoded sRNAs [8]. Trans-encoded sRNAs are those that are encoded within intergenic regions of the genome: they show a stable secondary structure and act on target sequences located at different positions in the genome. The complementarity between the sRNA and its target sequence is not complete and, for this reason, they require the presence of RNA chaperones to facilitate nucleotide binding [9]. In contrast, cis-encoded sRNAs originate in the nonsense strand of an open reading frame (ORF). Generally, that ORF corresponds to the sRNA target which presents complete complementarity to the cis-encoded sRNA. In addition, there have also been other types of RNA identified in Archaea: small nucleolar RNAs (snoRNAs) involved in the modifications of ribosomal RNA, whose presence was originally believed to be restricted to eukaryotic organisms [10]; sRNAs involved in the clustered regularly interspaced short palindromic repeats (CRISPR)/Cas prokaryotes immune system called crRNAs [11]; and sRNAs that derive from transfer RNAs called tRFs [12].
In Eukarya and Bacteria, many sRNAs have been characterised in detail. However, in domain Archaea, though there has been considerable progress in recent years, the number of sRNAs characterised is considerably lower than those in the other two domains [5,13,14]. The application of bioinformatic approaches and high-throughput sequencing systems for the analysis of complementary DNA (cDNA) libraries, RNA-Sequencing (RNA-Seq), have made it easier to understand the transcriptome and discover new sRNAs. However, to date, the identification of sRNAs using RNA-Seq analysis has only been carried out on seven species of Archaea under specific conditions: Archaeoglobus fulgidus [15], Sulfolobus solfataricus [16], Pyrococcus abyssi [10], Methanosarcina mazei [17], Haloferax volcanii [13], Pyrobaculum sp. [18], and Thermococcus kodakaraensis [19]. In these studies, it has been shown that the majority of sRNAs identified in Archaea are not conserved, even within species of the same genus [5]. This fact has also been observed in Bacteria, so it appears that the evolution of sRNA genes in prokaryotic organisms is greater than that of genes which encode proteins. Despite the large number of sRNAs identified in Archaea in recent years, very little is known about their biological functions and mechanisms of interaction, with their possible targets still an unknown and unexplored area of research.
Haloferax mediterranei is an extremely halophilic archaeon which belongs to the lineage of Euryarchaeota. This microorganism grows optimally at 2.5 M NaCl [20] in a defined medium with glucose as a carbon source and nitrate, nitrite, amino acids, or ammonium as sole nitrogen sources under aerobic conditions [21,22]. Most of the studies focused on nitrogen metabolism in halophilic archaea have been conducted using H. mediterranei as a haloarchaeal model. Specifically, biochemical, physiological, and transcriptomic studies of the assimilatory pathway in the presence of different nitrogen sources have been previously performed [22,23,24,25], as well as the development of molecular biology tools in H. mediterranei [26]. These results have revealed that the assimilatory pathway is highly regulated at a transcriptional level. However, little is known about the global regulatory networks that allow for the survival of this microorganism under stress or nitrogen starvation conditions. Besides transcriptional regulators analyses, the number of studies which show that sRNA is involved in a great variety of adaptive cellular responses to different stresses have surprisingly increased in the recent years.
In this work, we present an analysis of putative sRNAs in H. mediterranei expressed under two different nitrogen sources, nitrate and ammonium, using a combination of RNA-Seq and bioinformatic approaches. The main aim of this work is to extend the understanding of the global regulation network of the assimilatory pathway in H. mediterranei, specifically, and in domain Archaea, generally.

2. Materials and Methods

2.1. Strains and Growth Conditions

H. mediterranei strain R4 (ATCC 33500T) was grown at 42 °C with aeration at 225 rpm, contained in a 25% (w/v) mixture of inorganic salts (25% salt water) [27]. The pH value was adjusted to 7.3. H. mediterranei was grown in two different nitrogen sources, in a defined medium which contained 40 mM KNO3 or 40mM NH4Cl and supplemented with 5 g/L glucose, 0.0005 g/L FeCl3, and 0.5 g/L KH2PO4. Three independent biological replicates of each condition were employed.

2.2. Bioinformatic Prediction of small RNAs

The library of candidate sRNAs was constructed from putative sRNAs obtained in four different species of Archaea: H. volcanii [28], M. mazei [17], S. solfataricus [16], and A. fulgidus [15]. All putative sRNA sequences predicted from these four species were compared with the genome of H. mediterranei strain R4 (Genbank numbers: CP001868.2, CP001869.1, CP001870.1, and CP001871.1) using BLASTn [25]. The sequences obtained using the candidate sRNAs of other species (E-value < 0.05, p-value < 0.05) were prioritised using a function of identity, lower number of gaps (maximum 3), and lower number of mismatches.

2.3. RNA Isolation

For RNA isolation, H. mediterranei was grown in the presence of two different nitrogen sources, nitrate and ammonium, to mid-exponential phase. The mid-exponential growth phase was reached at different times and values of OD600 nm depending on the nitrogen source [23]. RNA was isolated with the mirVanaTM miRNA isolation kit (Ambion, Thermo Fisher Scientific, Waltham, MA, USA) following product specifications. Afterwards, the RNA samples were treated with Turbo DNase (Ambion, Thermo Fisher Scientific, Waltham, MA, USA). RNA concentration was analysed by means of a Nanodrop ND-100 Spectrophotometer (Thermo Fisher Scientific, Waltham, MA, USA), and the quality was analysed using the Small RNA Analysis Kit on Agilent 2200 Tapestation (Agilent Technologies, Santa Clara, CA, USA), respectively.

2.4. Library of Complementary DNA Preparation and Sequencing

Library preparation and sequencing were performed by the Bioarray, S.L. company (Alicante, Spain). sRNA libraries were constructed using TruSeq Small RNA Library Prep (Illumina, San Diego, CA, USA) and sequenced on an llumina HiSeq 2500 system using a 50-base pair (bp) read length. The corresponding FASTQ files were obtained as a result.

2.5. RNA-Sequencing Bioinformatic Analysis

Raw reads from each FASTQ file were aligned into BAM files—compressed binary file used to represent aligned sequences up to 128Mb [29] for each sample of each condition—through bowtie2 2.3.0 [30] using the H. mediterranei ATCC 33500 genome (Genbank numbers: GCF_000306765.2_ASM30676v2; Assembly: GCA_000306765.2) as a reference. Seven ammonium BAM files and eight nitrogen BAM files were viewed in the Integrative Genomics Viewer (IGV) program [31,32], alongside the annotated H. mediterranei genome.
The raw readings, which aligned at the positions of the library of candidate sRNAs obtained from the different species against the H. mediterranei genome, were analysed manually using the Integrative Genomics Viewer (IGV). Only sequences with a score of at least 20 (frequency of reads in high-throughput sequencing) and that were present in the six samples of each condition were selected as sRNAs in H. mediterranei. The candidate sRNAs were classified as cis-encoded sRNA if they were encoded in the reverse direction of ORF, as trans-encoded sRNA if they were in intergenic regions, and as crRNA if they were in a CRISPR array. Since the size of the readings was generally greater than the homology shown in the different species analysed, the BAM files of each condition were pooled, thus, increasing the coverage of the readings, allowing more precise limiting of the sRNAs location in the H. mediterranei genome. Once the positions were obtained, the approximate sizes of the sRNAs were calculated. In addition, the genetic environment of each sRNA was analysed manually.
The counts of the sRNAs of H. mediterranei were obtained using the intersection nonempty feature from the HT-Seq program [33]. In this way, for each position i in the read, a set S(i) is defined as the set of all features overlapping the i position. Then, if S contains precisely one feature, the read is counted for this feature. If S is empty, the read is counted as no feature. Finally, if S contains more than one feature, the read is considered as ambiguous and is not counted for any features.
Differential expression analysis in the function of the nitrogen source (nitrate/ammonium) was performed using the DESeq2 library from the Bioconductor 3.5 package [34] obtaining the fold-change of each sRNA (p-value < 0.01 and p-adj < 0.05). Through this method, the variance-mean dependence in count data was estimated and tested for differential expression; this was based on a model using the negative binomial distribution between nitrate and ammonium group samples.
Both raw (FASTAQ files) and processed data (normalised counts) are available on the Gene Expression Omnibus (GEO) database (Series entry number: GSE108616) [35].
Mfold was used to predict the secondary structure of sRNAs obtained [36]. Potential gene targets for each sRNA were identified using TargetRNA2 [37]. Moreover, the cis-encoded sRNAs antisense of characterised ORF were analysed using IntaRNA [38]. BLASTn [25] was used to search homology regions in another organism.

2.6. Validation of sRNAs Using Reverse Transcription Polymerase Chain Reaction

The validation of 20 sRNAs was performed by reverse transcription polymerase chain reaction (RT-PCR), prioritising those with differential expression based on the nitrogen source and/or target genes of known functions. The RNA was isolated as described before. Between 0.5 and 0.8 µg of DNA-free RNA was used for the synthesis of cDNA using M-MuLV Reverse Transcriptase (Thermo Fisher Scientific, Waltham, MA, USA) and random hexamer primer (Thermo Fisher Scientific, Waltham, MA, USA) according to the manufacturer’s instructions. Negative controls for reverse transcription polymerase chain reaction (RT-PCR) were prepared by omitting reverse transcriptase and cDNA. The oligonucleotides used to carry out the RT-PCR were designed based on the sequence of sRNA candidates (Table S1). Amplified products were analysed using 3% agarose gel electrophoresis with GeneRulerTM Low Range DNA Ladder (Thermo Fisher Scientific, Waltham, MA, USA) run in parallel. PCR products were purified with the IllustraTM GFXTM PCR DNA and Gel Band Purification Kit (GE Healthcare, Little Chalfont, UK) and confirmed by Sanger sequencing (Stabvida, Caparica, Portugal).

3. Results

3.1. Identification of Putative small RNAs in Haloferax mediterranei

The library of candidate sRNAs was generated from putative sRNAs obtained in four different species of domain Archaea: H. volcanii [28], M. mazei [17], S. solfataricus [16], and A. fulgidus [15]. Since the candidate sRNAs of H. volcanii aligned at many positions in the genome of H. mediterranei (E-value < 0.05, p-value < 0.05), the results obtained were prioritised by a function of identity, lower number of gaps (maximum 3), and lower number of mismatches. Given the evolutionary distance, the candidate sRNAs of the remaining species aligned in fewer positions (E-value < 1) than H. volcanii in the genome of H. mediterranei, so it was not necessary to prioritise in the same way. From the bioinformatic analysis, a library of 295 candidate sRNAs in H. mediterranei were obtained (Table 1). The number of mismatches allowed depends on the length of the sequence, permitting a greater number of mismatches in larger sequences with a high identity percentage. The results obtained after alignment with the candidate sRNAs of other species are shown in detail in the supplementary material (Tables S2–S5).

3.2. Verification of Predicted 295 sRNAs Using RNA-Sequencing

The 295 sRNAs that were bioinformatically predicted were verified with the RNA-Seq results of the H. mediterranei ATCC33500. The raw aligned readings obtained from the analysis of RNA-Seq results and the 295 candidate sRNA sequences obtained by means of the bioinformatic approach were manually verified using the high-performance visualisation tool, IGV [31,32]. The transcripts expressed from intergenic regions, antisense to the characterised ORFs, and in the CRISPR array were classified in trans-encoded sRNA, cis-encoded sRNA, and crRNA, respectively. Due to the size of the readings being generally greater than the homology shown in different species being analysed, the BAM files of each condition were pooled, thus increasing the coverage of the readings and allowing for more accurate localization of the sRNAs in the H. mediterranei genome. Once the positions were obtained, the approximate size of the sRNAs was calculated. They were in the range of 20–500 nt. The results of the genetic environment analysis of the 88 verified sRNAs in H. mediterranei and their sequences are shown in Table 2 and Table S6, respectively.

3.3. Bioinformatic Analysis of Putative Small RNAs in H. mediterranei

3.3.1. Identification and Classification of small RNAs Verified by RNA-Sequencing

The results of the localisation and classification of the sRNAs verified by RNA-Seq reveals that 58 sRNAs of the 116 sRNAs predicted by homology with H. volcanii (Table 1 and Table S2) could be assigned as sRNAs in H. mediterranei. There were 56 sRNAs on the chromosome and two on the plasmid pHM500. According to the classification established above, 51 sRNAs are intergenic sRNAs, six are antisense sRNAs, and one is crRNA. The homology analysis performed with M. mazei results in 55 predicted sRNAs (Table 1 and Table S3), of which 16 could be sRNAs in H. mediterranei. Twelve sRNAs are located on the chromosome, three on the plasmid pHM500, and one on the plasmid pHM100. Related to the classification according to the location of the target, 11 sRNAs are intergenic sRNAs and five antisense. Of the 47 sRNAs predicted by homology with S. solfataricus (Table 1 and Table S4), only eight sRNAs were identified in the results of RNA-Seq. These are located on the chromosome, and are classified as four intergenic sRNAs and four antisense sRNAs. From the analysis with A. fulgidus (Table 1 and Table S5), only six candidate sRNAs from a total of 74 sRNAs predicted by homology could be assigned as sRNAs. As in the analysis of S. solfataricus, all of them are located on the chromosome, with four being intergenic sRNAs and two being antisense sRNAs. From a total of 88 sRNAs verified in H. mediterranei by RNA-Seq, 93.18% were located on the chromosome, with most of them (79.54%) being intergenic (trans-encoded sRNA).

3.3.2. Structure and Targets of 88 small RNAs

Usually, sRNAs are characterised by stable secondary structures. In order to analyse the secondary structure of the 88 putative sRNAs, Mfold software [36] was used with default parameters set and with the temperature modified to 42 °C. The core algorithm predicts a minimum free energy, ΔG°, as well as minimum free energies for folding that must contain any particular base pair. All of the sRNAs were shown to have significant predicted secondary structures since they presented with ΔG < 0. However, only the secondary structures of 16 sRNAs with differential expression according to the nitrogen source (nitrate/ammonium) are shown in Figure 1. Furthermore, most secondary structures of sRNAs revealed highly structured molecules, including more than one hairpin loop and high stability (ΔG < 0). It is noteworthy that the structures of HM39_V, HM54_V, HM38_V, and HM8_S sRNAs present predominant structures with two or three hairpin loops and are similar to the structures of other characterised sRNAs [39,40].
The identification of potential sRNA targets was carried out via a bioinformatic approach using TargetRNA2 [37]. Gene targets were found for 71.59% of the sRNA verified, of which the first five are most likely due to their energy values (Table S8). Taking the five gene targets for each sRNA analysed into account, the most commonly predicted targets correspond to hypothetical proteins with unknown functions. Many examples of transcriptional regulators such as the ArsR family, nitrogen regulatory PII protein, transporters (including metal and ATP-binding cassette (ABC) transporters), proteins related to RNA (i.e., H/ACA RNA-protein complex component Gar1, ribonuclease P subunit p30, 50S ribosomal protein L24e, 50S ribosomal protein L37Ae, DNA-directed RNA polymerase subunit B″, DNA-directed RNA polymerase subunit F, 30S ribosomal protein S8P, etc.), and proteins related to DNA metabolism (i.e., DNA lyase, DNA ligase, DNA double-strand break repair protein mre11, DNA primase large subunit) were also identified. To go in depth with this analysis, the cis-encoded sRNAs antisense of characterised ORF were analysed using IntaRNA [38] (Table 3). The results show that nine of the 17 cis-encoded sRNAs analysed present high interaction energy and full complementarity with the ORF where they were identified. The combination of these software programs has allowed the identification of possible gene targets for more than 70% of the 88 possible sRNAs identified by RNA-Seq. However, there are still 26 sRNAs whose possible targets could not be found with the software employed.

3.3.3. Conservation of small RNAs Verified in H. mediterranei

The sequence conservation of the sRNAs predicted in H. mediterranei was analysed using BLASTn [25], comparing each sequence to all sequenced archaeal genomes (E value 10 × 10−6). Only hits with a nucleotide identity higher than 60% combined with a coverage between the query and subject sequence higher than 80% were considered as conserved. The results of homology are shown in Figure 2 and Table S7. It was found that 27% of verified sRNAs present conserved sequences in other Halobacteria, mainly in the Haloferacales order (86%) and Haloferacaceae family (79%). The remaining 63 sRNAs showed no sequence homology with any other Archaea.

3.4. Expression Analysis According to the Nitrogen Source

The expression analysis of the sRNAs in H. mediterranei according to the nitrogen source was performed using the HT-Seq program, which counts the readings of the sRNAs verified by RNA-Seq and DESeq2, which carries out the differential expression analysis (nitrate/ammonium). From the 88 candidate sRNAs analysed, 16 sRNAs met the statistical criteria with a p-value and p-adj lower than 0.02 and 0.05, respectively (Table 4). These sRNAs show differences in their expression pattern. Eight of them are overexpressed in the presence of nitrate as a sole nitrogen source (log2-fold-change between 0.519 and 8.699), whereas the other eight present a decrease in their transcriptional level in presence of ammonium (log2-fold-change between −0.478 and 1.325).

3.5. Validation of Expression of Small RNAs Using Reverse Transcription Polymerase Chain Reaction

To confirm the presence of 20 sRNAs predicted by RNA-Seq we used RT-PCR [41,42,43,44]. Of 20 predicted sRNAs tested, nine (40%) gave positive results of the expected size as shown in Figure 3. These PCR products were subsequently verified using Sanger sequencing. This proportion of successfully validated sRNAs (45%) agrees with the results obtained in other works, where validation is often successful around 40–50% [39,42].
The majority of sRNAs (HM6_S, HM37_V, HM7_S, HM1_S, HM3_M, and HM32_V) are expressed in both conditions, in the presence of nitrate and ammonium. However, three sRNAs show a different expression pattern according to the nitrogen source; HM8_S and HM1_M sRNAs are expressed exclusively in the presence of nitrate and HM16_M sRNA is expressed in the presence of ammonium.

4. Discussion

Two different approaches, bioinformatic and RNomic, have allowed the identification of 88 sRNAs in H. mediterranei. Most of them are located on the chromosome, as in the case of H. volcanii. However, there is a difference between these two halophilic microorganisms. H. mediterranei sRNAs are predominantly intergenic (trans-encoded sRNA), whereas H. volcanii sRNAs are mostly antisense instead [28,45].
The 88 sRNAs verified by RNA-Seq were analysed using different strategies to obtain further information about their characteristics and their possible involvement in different cellular processes. Some of the predicted secondary structures of these sRNAs remain similar to other sRNA structures which have been analysed in detail (Figure 1) [39,40]. According to Gaimster et al. [39], these results are significant because they imply that many of these sRNAs could have the potential to form complex conformations like those commonly associated with many other directly-acting RNA transcripts, including known bacterial sRNAs. Moreover, additional molecular mechanisms of sRNA functions in bacteria have been described, including the destabilisation or stabilisation of the target mRNAs. These actions are conducted via specific binding to a protein or full complementarity to their target mRNA [5]. Thus, putative gene targets of all the sRNAs candidates in H. mediterranei were predicted using TargetRNA2 [37]. As it has been observed in other works [39], the results obtained reveal that the most commonly predicted targets matched to hypothetical proteins. This information could potentially be useful in eventually assigning a function to these unknown proteins. It has also been identified to target some transcriptional regulators such as the ArsR family, as well as the nitrogen regulatory PII protein. Interestingly, the expression of the ArsR and PII proteins in H. mediterranei is closely related to the nitrogen source [24], so the existence of a regulatory network where the sRNAs act by activating or repressing the expression of these proteins is highly likely. In other bacterial studies, it has been observed that global regulators are subject to regulation by multiple Hfq-dependent sRNAs [46]. In the case of domain Archaea, similar interactions have been described between sRNAs and Like Sm (Lsm) homologous proteins, instead of the bacterial Hfq [47]. H. mediterranei contains one gene encoding the Lsm1 protein (HFX_2733) in its genome that shares 99% sequence identity with the H. volcanii RNA-binding Lsm protein (E value 6 × 10−51). These results suggest that the interaction between Lsm proteins and sRNAs and their participation in the RNA metabolism are also possible in H. mediterranei. Transporters are also included in the different targets predicted, specifically metal and ABC transporters, which have been observed in other analyses of bacterial sRNAs such as Paracoccus denitrificans [39] and Ruegeria pomeroyi [48], as well as in Archaeal sRNA such as S. solfataricus [49] and H. volcanii [50]. Based on these results, it could be possible that the regulation of transporters is not a characteristic only of the Bacteria domain. This suggests that antisense sRNA might be a common regulatory mechanism for such genes. It has also been identified as a target of different proteins related to RNA and DNA metabolism, which support the concept that sRNAs play an essential role in transcriptional and post-transcriptional regulation. Furthermore, analysis of co-purification with ribosomal protein L7Ae in S. solfataricus allowed for the identification of one sRNA in that microorganism which suggests a direct interaction of these proteins with sRNAs [51]. In parallel, the cis-encoded sRNAs antisense of ORF characterised in H. mediterranei were also analysed using IntaRNA [38]. This analysis showed positive results in nine of the 17 cis-encoded sRNAs. Four of them were present in only one target and the other five could interact with other possible targets. Most sRNAs that did not exhibit interaction with RNA using IntaRNA did exhibit interaction with other mRNAs using TargetRNA2. Therefore, it is likely that some sRNAs are specific to a gene, while others have a broader range of action, being involved in several metabolic pathways (Table S8). Only 26 of the 88 sRNAs analysed did not show positive results using either IntaRNA or TargetRNA2. Hence, more than 70% of the 88 sRNAs verified by RNA-Seq present possible target genes using such software. Nevertheless, it is clear that the only way to confirm the interaction between RNA and its target is through experimental validation. Therefore, in future works we will try to characterise the targets of the sRNAs that present a differential expression in the function of the nitrogen source through the interaction with the protein Lsm1, the generation of deletion sRNA mutants, and other analyses (such as overexpression mutants) to demonstrate the influence of these sRNAs on the respective target.
BLASTn [25] analysis reveals that 27% of sRNAs identified in H. mediterranei present conserved sequences in other Halobacteria. Consequently, it is likely that these conserved sRNAs may play a conserved role in closely related species. However, despite the fact that the sequences of the sRNAs are slightly conserved in Archaea, this does not imply that they are all true sRNA candidates. The remaining 72% of sRNAs identified showed no sequence homology to any other archaea. According to Babski et al. [5], many sRNA genes of Archaea are not even shared by species of the same genus. Therefore, this 72% of the sRNA could be specific to the halophilic archaea H. mediterranei.
The transcript pattern analysis of H. mediterranei sRNA according to the nitrogen source demonstrates that 16 sRNA showed higher or lower transcript levels with statistically significant parameters (p-value < 0.02 and p-adj < 0.05). Of these, eight show overexpression when H. mediterranei grows with nitrate as the sole nitrogen source, whereas the other eight show a decrease in their transcriptional level in the presence of ammonium. The majority of sRNAs validated by RT-PCR are expressed in the presence of nitrogen sources nitrate and ammonium. However, HM8_S and HM16_M sRNAs are expressed exclusively in the presence of nitrate and ammonium, respectively. Curiously, although the HM1_M sRNA transcript level does not show differences in the RNA-Seq analysis according to the nitrogen source, RT-PCR results clearly reveal that this sRNA is only expressed in the presence of nitrate (Figure 3).
Interestingly, some of the sRNAs which show differences in their expression pattern according to the nitrogen source possibly have target genes whose expression also depends on the nitrogen source. HM8_S, which is overexpressed in nitrate as the nitrogen source, has the protein glutamate dehydrogenase as a possible target (Table 3). Glutamate dehydrogenase is underexpressed in nitrate’s presence. Therefore, this sRNA could be involved in the repression of this key enzyme in the metabolism of nitrogen. Moreover, two sRNAs (i.e., HM7_S and HM54_V) are slightly more expressed in the presence of nitrate. Curiously, their possible targets are transcriptional regulators belonging to the ArsR family. Previous results revealed that the expression of different ArsR proteins depends on the nitrogen source [24]. Therefore, sRNAs could also be involved in the regulation of the expression of these transcriptional regulators, and it could be an example of a gene regulatory network related to nitrogen assimilation. Finally, HM1_A sRNA is expressed in the presence of ammonium as a nitrogen source, and it has been found to be a possible target of the amt1 gene (Table 3), which encodes an ammonium transporter. Different studies performed with H. mediterranei have confirmed that ammonium transporters are expressed in the presence of nitrate or under nitrogen starvation [24]. Hence, this sRNA could also be related to the regulation of Amt transporters expression and, consequently, be adjusting the uptake of ammonium from the medium. More work is needed to confirm these hypotheses and to find out the role of these sRNAs in nitrogen metabolism.

5. Conclusions

This work, focused on identifying sRNAs involved in nitrogen assimilation, has also increased the knowledge about sRNAs in the domain Archaea. Specifically, 88 sRNAs have been identified in H. mediterranei using bioinformatic and RNomic approaches, some of which show different expression patterns depending on the nitrogen source and/or present genes involved in the nitrogen assimilation as a potential gene target. This data suggests that some of these sRNAs could be related to the regulation of nitrogen assimilation, being able to constitute an important gene regulatory network which involves enzymes, transporters, and transcriptional regulators in this metabolism. Undoubtedly, this work constitutes an excellent starting point to elucidate the role of these sRNAs in the nitrogen metabolism of haloarchaea.

Supplementary Materials

The following are available online at https://www.mdpi.com/2073-4425/9/2/83/s1. Table S1: Nucleotide sequences of PCR primers for amplifying sRNA genes. Table S2: Library of sRNAs obtained after alignment of H. volcanii sRNAs [28] against the H. mediterranei genome (E-value < 0.05). Table S3: Library of sRNAs obtained after alignment of M. mazei sRNAs [17] against the H. mediterranei genome (E-value < 1). Table S4: Library of sRNAs obtained after alignment of S. solfataricus sRNAs [16] against the H. mediterranei genome (E-value < 1). Table S5: Library of sRNAs obtained after alignment of A. fulgidus sRNAs [15] against the H. mediterranei genome (E-value < 1). Table S6: Sequences of 88 candidates sRNAs identified in H. mediterranei by RNA-Seq. Table S7: Conservation of the sequences of sRNAs across different orders, classes, and species. BLASTn, homology (E value 10 × 10−6, query cover at least 80% and sequence identity at least 60%). Table S8: Prediction of putative gene targets of 88 sRNA in H. mediterranei using TargetRNA2 [36].

Acknowledgments

We would like to thank Jörg Soppa (Goethe-University, Frankfurt, Germany) for useful comments and assistance in this work. This work was supported by grant BIO2013-42921-P from the Spanish Ministry of Economy and Competitiveness. Funds from “Programa Propio para el fomento de la I+D+i del Vicerrectorado de Investigación y Transferencia de Conocimiento” de la Universidad de Alicante are also acknowledged.

Author Contributions

J.E. and M.-J.B. conceived and designed the experiments; G.P., J.E., V.B., and M.C. performed most of the experiments; N.C.-F. and L.A.A. carried out specific experimental parts; G.P., J.E., N.C.-F., and M.-J.B. analysed the data; M.-J.B. supervised the experimental work; G.P., J.E., and M.-J.B. wrote the paper. V.B. and M.C. contributed to the writing of the manuscript.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Chu, C.Y.; Rana, T.M. Small RNAs: Regulators and guardians of the genome. J. Cell. Physiol. 2007, 213, 412–419. [Google Scholar] [CrossRef] [PubMed]
  2. Choudhuri, S. Lesser known relatives of miRNA. Biochem. Biophys. Res. Commun. 2009, 388, 177–180. [Google Scholar] [CrossRef] [PubMed]
  3. Liu, Q.; Paroo, Z. Biochemical principles of small RNA pathways. Annu. Rev. Biochem. 2010, 79, 295–319. [Google Scholar] [CrossRef] [PubMed]
  4. Storz, G.; Vogel, J.; Wassarman, K.M. Regulation by small RNAs in bacteria: Expanding frontiers. Mol. Cell. 2011, 43, 880–891. [Google Scholar] [CrossRef] [PubMed]
  5. Babski, J.; Maier, L.K.; Heyer, R.; Jaschinski, K.; Prasse, D.; Jäger, D.; Randau, L.; Schmitz, R.A.; Marchfelder, A.; Soppa, J. Small regulatory RNAs in Archaea. RNA Biol. 2014, 11, 484–493. [Google Scholar] [CrossRef] [PubMed]
  6. Michaux, C.; Hartke, A.; Martini, C.; Reiss, S.; Albrecht, D.; Budin-Verneuil, A.; Sanguinetti, M.; Engelmann, S.; Hain, T.; Verneuil, N.; et al. Involvement of Enterococcus faecalis small RNAs in stress response and virulence. Infect. Immun. 2014, 82, 3599–3611. [Google Scholar] [CrossRef] [PubMed]
  7. Durand, S.; Tomasini, A.; Braun, F.; Condon, C.; Romby, P. sRNA and mRNA turnover in Gram-positive bacteria. FEMS Microbiol. Rev. 2015, 39, 316–330. [Google Scholar] [CrossRef] [PubMed]
  8. Gottesman, S.; Storz, G. Bacterial small RNA regulators: Versatile roles and rapidly evolving variations. Cold Spring Harb. Perspect. Biol. 2011, 3. [Google Scholar] [CrossRef] [PubMed]
  9. Moll, I.; Leitsch, D.; Steinhauser, T.; Bläsi, U. RNA chaperone activity of the Sm-like Hfq protein. EMBO Rep. 2003, 4, 284–289. [Google Scholar] [CrossRef] [PubMed]
  10. Muller, S.; Leclerc, F.; Behm-Ansmant, I.; Fourmann, J.B.; Charpentier, B.; Branlant, C. Combined in silico and experimental identification of the Pyrococcus abyssi H/ACA sRNAs and their target sites in ribosomal RNAs. Nucleic Acids Res. 2008, 36, 2459–2475. [Google Scholar] [CrossRef] [PubMed]
  11. Phok, K.; Moisan, A.; Rinaldi, D.; Brucato, N.; Carpousis, A.J.; Gaspin, C.; Clouet-d’Orval, B. Identification of CRISPR and riboswitch related RNAs among novel noncoding RNAs of the euryarchaeon Pyrococcus abyssi. BMC Genom. 2011, 12, 312. [Google Scholar] [CrossRef] [PubMed]
  12. Heyer, R.; Dörr, M.; Jellen-Ritter, A.; Späth, B.; Babski, J.; Jaschinski, K.; Soppa, J.; Marchfelder, A. High throughput sequencing reveals a plethora of small RNAs including tRNA derived fragments in Haloferax volcanii. RNA Biol. 2012, 9, 1011–1018. [Google Scholar] [CrossRef] [PubMed]
  13. Soppa, J.; Straub, J.; Brenneis, M.; Jellen-Ritter, A.; Heyer, R.; Fischer, S.; Granzow, M.; Voss, B.; Hess, W.R.; Tjaden, B.; et al. Small RNAs of the halophilic archaeon Haloferax volcanii. Biochem. Soc. Trans. 2009, 37, 133–136. [Google Scholar] [CrossRef] [PubMed]
  14. Wagner, E.G.; Romby, P. Small RNAs in bacteria and archaea: Who they are, what they do, and how they do it. Adv. Genet. 2015, 90, 133–208. [Google Scholar] [CrossRef] [PubMed]
  15. Tang, T.H.; Bachellerie, J.P.; Rozhdestvensky, T.; Bortolin, M.L.; Huber, H.; Drungowski, M.; Elge, T.; Brosius, J.; Hüttenhofer, A. Identification of 86 candidates for small non-messenger RNAs from the archaeon Archaeoglobus fulgidus. Proc. Natl. Acad. Sci. USA 2002, 99, 7536–7541. [Google Scholar] [CrossRef] [PubMed]
  16. Tang, T.-H.; Polacek, N.; Zywicki, M.; Huber, H.; Brugger, K.; Garrett, R.; Bachellerie, J.P.; Hüttenhofer, A. Identification of novel non-coding RNAs as potential antisense regulators in the archaeon Sulfolobus solfataricus. Mol. Microbiol. 2005, 55, 469–481. [Google Scholar] [CrossRef] [PubMed]
  17. Jäger, D.; Sharma, C.M.; Thomsen, J.; Ehlers, C.; Vogel, J.; Schmitz, R.A. Deep sequencing analysis of the Methanosarcina mazei Gö1 transcriptome in response to nitrogen availability. Proc. Natl. Acad. Sci. USA 2009, 106, 21878–21882. [Google Scholar] [CrossRef] [PubMed]
  18. Bernick, D.L.; Dennis, P.P.; Höchsmann, M.; Lowe, T.M. Discovery of Pyrobaculum small RNA families with atypical pseudouridine guide RNA features. RNA 2012, 18, 402–411. [Google Scholar] [CrossRef] [PubMed]
  19. Jäger, D.; Förstner, K.U.; Sharma, C.M.; Santangelo, T.J.; John, N.; Reeve, J.N. Primary transcriptome map of the hyperthermophilic archaeon Thermococcus kodakarensis. BMC Genom. 2014, 15, 684. [Google Scholar] [CrossRef] [PubMed]
  20. Torreblanca, M.; Rodriguez-Valera, F.; Juez, G.; Ventosa, A.; Masahiro, K.; Morris, K. Classification of non-alkaliphilic halobacteria based on numerical taxonomy and polar lipid composition, and description of Haloarcula gen. nov. and Haloferax gen. nov. Syst. Appl. Microbiol. 1986, 8, 89–99. [Google Scholar] [CrossRef]
  21. Rodríguez-Valera, F.; Juez, G.; Kushner, D.J. Halobacterium mediterranei spec, nov., a new carbohydrate-utilizing extreme halophile. Syst. Appl. Microbiol. 1983, 4, 369–381. [Google Scholar] [CrossRef]
  22. Bonete, M.J.; Martínez-Espinosa, R.M.; Pire, C.; Zafrilla, B.; Richardson, D.J. Nitrogen metabolism in haloarchaea. Saline Syst. 2008, 4, 9. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  23. Esclapez, J.; Bravo-Barrales, G.; Bautista, V.; Pire, C.; Camacho, M.; Bonete, M.J. Effects of nitrogen sources on the nitrate assimilation in Haloferax mediterranei: Growth kinetics and transcriptomic analysis. FEMS Microbiol. Lett. 2014, 350, 168–174. [Google Scholar] [CrossRef] [PubMed]
  24. Esclapez, J.; Pire, C.; Camacho, M.; Bautista, V.; Martínez-Espinosa, R.M.; Zafrilla, B.; Vegara, A.; Alcaraz, L.A.; Bonete, M.J. Transcriptional profiles of Haloferax mediterranei based on nitrogen availability. J. Biotechnol. 2015, 193, 100–107. [Google Scholar] [CrossRef] [PubMed]
  25. Esclapez, J.; Camacho, M.; Pire, C.; Bautista, V.; Vegara, A.; Pedro-Roig, L.; Pérez-Pomares, F.; Martínez-Espinosa, R.M.; Bonete, M.J. Recent Advances in The Nitrogen Metabolism in Haloarchaea and Its Biotechnological Applications. In Biotechnology of Extremophiles: Advances and Challenges; Rampelotto, P., Ed.; Springer: Cham, Switzerland, 2016; Volume 1, pp. 273–301. ISBN 978-3-319-13520-5. [Google Scholar]
  26. Pedro-Roig, L.; Lange, C.; Bonete, M.J.; Soppa, J.; Maupin-Furlow, J. Nitrogen regulation of protein-protein interactions and transcript levels of GlnK PII regulator and AmtB ammonium transporter homologs in Archaea. Microbiologyopen 2013, 2, 826–840. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  27. Rodríguez-Valera, F.; Ruiz-Berraquero, F.; Ramos-Cormenzana, A. Behaviour of mixed populations of halophilic bacteria in continuous cultures. Can. J. Microbiol. 1980, 26, 1259–1263. [Google Scholar] [CrossRef] [PubMed]
  28. Babski, J.; Tjaden, B.; Voss, B.; Jellen-Ritter, A.; Marchfelder, A.; Hess, W.R.; Soppa, J. Bioinformatic prediction and experimental verification of sRNAs in the haloarchaeon Haloferax volcanii. RNA Biol. 2011, 8, 806–816. [Google Scholar] [CrossRef] [PubMed]
  29. Altschul, S.F.; Gish, W.; Miller, W.; Myers, E.W.; Lipman, D.J. Basic local alignment search tool. J. Mol. Biol. 1990, 215, 403–410. [Google Scholar] [CrossRef]
  30. Langmead, B.; Salzberg, S. Fast gapped-read alignment with Bowtie 2. Nat. Methods 2012, 9, 357–359. [Google Scholar] [CrossRef] [PubMed]
  31. Robinson, J.T.; Thorvaldsdóttir, H.; Winckler, W.; Guttman, M.; Lander, E.S.; Getz, G.; Mesirov, J.P. Integrative genomics viewer. Nat. Biotechnol. 2011, 29, 24–26. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  32. Thorvaldsdóttir, H.; Robinson, J.T.; Mesirov, J.P. Integrative genomics viewer (IGV): High-performance genomics. Brief. Bioinform. 2013, 14, 178–192. [Google Scholar] [CrossRef] [PubMed]
  33. Anders, S.; Pyl, P.T.; Huber, W. HTSeq—A Python framework to work with high-throughput sequencing data. Bioinformatics 2015, 31, 166–169. [Google Scholar] [CrossRef] [PubMed]
  34. Love, M.I.; Huber, W.; Anders, S. Moderated estimation of fold change and dispersion for RNA-Seq data with DESeq2. Genome Biol. 2014, 15, 550. [Google Scholar] [CrossRef] [PubMed]
  35. Edgar, R.; Domrachev, M.; Lash, A.E. Gene Expression Omnibus: NCBI gene expression and hybridization array data repository. Nucleic Acids Res. 2002 30, 207–210. [CrossRef]
  36. Zuker, M. Mfold web server for nucleic acid folding and hybridization prediction. Nucleic Acids Res. 2003, 31, 3406–3415. [Google Scholar] [CrossRef] [PubMed]
  37. Kery, M.B.F.M.; Livny, J.; Tjaden, B. TargetRNA2: Identifying targets of small regulatory RNAs in bacteria. Nucleic Acids Res. 2014, 42, 124–129. [Google Scholar] [CrossRef] [PubMed]
  38. Mann, M.; Wright, P.R.; Backofen, R. IntaRNA 2.0: Enhanced and customizable prediction of RNA–RNA interactions. Nucleic Acids Res. 2017, 45, 435–439. [Google Scholar] [CrossRef] [PubMed]
  39. Gaimster, H.; Chalklen, L.; Alston, M.; Munnoch, J.T.; Richardson, D.J.; Gates, A.J.; Rowley, G. Genome-wide discovery of putative sRNAs in Paracoccus denitrificans expressed under nitrous oxide emitting conditions. Front. Microbiol. 2016, 7, 1806. [Google Scholar] [CrossRef] [PubMed]
  40. Prasse, D.; Förstner, K.U.; Jäger, D.; Backofen, R.; Schmitz, R.A. sRNA154 a newly identified regulator of nitrogen fixation in Methanosarcina mazei strain Gö1. RNA Biol. 2017, 14, 1544–1558. [Google Scholar] [CrossRef] [PubMed]
  41. Pánek, J.; Bobek, J.; Mikulík, K.; Basler, M.; Vohradský, J. Biocomputational prediction of small non-coding RNAs in Streptomyces. BMC Genom. 2008, 9, 217. [Google Scholar] [CrossRef] [PubMed]
  42. Khoo, J.S.; Chai, S.F.; Mohamed, R.; Nathan, S.; Firdaus-Raih, M. Computational discovery and RT-PCR validation of novel Burkholderia conserved and Burkholderia pseudomallei unique sRNAs. BMC Genom. 2012, 13. [Google Scholar] [CrossRef]
  43. Panda, G.; Tanwer, P.; Ansari, S.; Khare, D.; Bhatnagar, R. Regulation and RNA-binding properties of Hfq-like RNA chaperones in Bacillus anthracis. Biochim. Biophys. Acta 2015, 1850, 1661–1668. [Google Scholar] [CrossRef] [PubMed]
  44. Kwenda, S.; Gorshkov, V.; Ramesh, A.M.; Naidoo, S.; Rubagotti, E.; Birch, P.R.; Moleleki, L.N. Discovery and profiling of small RNAs responsive to stress conditions in the plant pathogen Pectobacterium atrosepticum. BMC Genom. 2016, 12, 47. [Google Scholar] [CrossRef] [PubMed]
  45. Babski, J.; Haas, K.A.; Näther-Schindler, D.; Pfeiffer, F.; Förstner, K.U.; Hammelmann, M.; Hilker, R.; Becker, A.; Sharma, C.M.; Marchfelder, A.; Soppa, J. Genome-wide identification of transcriptional start sites in the haloarchaeon Haloferax volcanii based on differential RNA-Seq (dRNA-Seq). BMC Genom. 2016, 17, 629. [Google Scholar] [CrossRef] [PubMed]
  46. Lee, H.J.; Gottesman, S. sRNA roles in regulating transcriptional regulators: Lrp and SoxS regulation by sRNAs. Nucleic Acids Res. 2016, 44, 6907–6923. [Google Scholar] [CrossRef] [PubMed]
  47. Fischer, S.; Benz, J.; Späth, B.; Maier, L.K.; Straub, J.; Granzow, M.; Raabe, M.; Urlaub, H.; Hoffmann, J.; Brutschy, B.; et al. The archaeal Lsm protein binds to small RNAs. J. Biol. Chem. 2010, 285, 45, 34429–34438. [Google Scholar] [CrossRef] [PubMed]
  48. Rivers, A.R.; Burns, A.S.; Chan, L.K.; Moran, M.A. Experimental identification of small non-coding RNAs in the model marine bacterium Ruegeria pomeroyi DSS-3. Front. Microbiol. 2016, 29, 380. [Google Scholar] [CrossRef] [PubMed]
  49. Wurtzel, O.; Sapra, R.; Chen, F.; Zhu, Y.; Simmons, B.A.; Sorek, R. A single-base resolution map of an archaeal transcriptome. Genome Res. 2010, 20, 133–141. [Google Scholar] [CrossRef] [PubMed]
  50. Straub, J.; Brenneis, M.; Jellen-Ritter, A.; Heyer, R.; Soppa, J.; Marchfelder, A. Small RNAs in haloarchaea: Identification, differential expression and biological function. RNA Biol. 2009, 6, 281–292. [Google Scholar] [CrossRef] [PubMed]
  51. Zago, M.A.; Dennis, P.P.; Omer, A.D. The expanding world of small RNAs in the hyperthermophilic archaeon Sulfolobus solfataricus. Mol. Microbiol. 2005, 55, 1812–1828. [Google Scholar] [CrossRef] [PubMed]
Figure 1. Prediction of the secondary structure of 16 sRNAs with differential expression (nitrate/ammonium) in Haloferax mediterranei using Mfold [36]. (A) sRNAs overexpressed in nitrate. (B) sRNAs overexpressed in ammonium.
Figure 1. Prediction of the secondary structure of 16 sRNAs with differential expression (nitrate/ammonium) in Haloferax mediterranei using Mfold [36]. (A) sRNAs overexpressed in nitrate. (B) sRNAs overexpressed in ammonium.
Genes 09 00083 g001
Figure 2. Conservation percentage of the sequences of sRNAs across different (A) classes, (B) orders, and (C) families. BLASTn, homology (E value 10 × 10−6); query covered at least 80% and sequence identity covered at least 60%.
Figure 2. Conservation percentage of the sequences of sRNAs across different (A) classes, (B) orders, and (C) families. BLASTn, homology (E value 10 × 10−6); query covered at least 80% and sequence identity covered at least 60%.
Genes 09 00083 g002
Figure 3. Reverse transcription polymerase chain reaction (RT-PCR) validation of 20 sRNAs. Electrophoresis of PCR products of 20 sRNAs on 3% agarose gel. Lanes 1 and 12: GeneRuler Low Range DNA Ladder (Thermo Scientific, Waltham, MA, USA). Lanes 2–11: PCR products of sRNAs.
Figure 3. Reverse transcription polymerase chain reaction (RT-PCR) validation of 20 sRNAs. Electrophoresis of PCR products of 20 sRNAs on 3% agarose gel. Lanes 1 and 12: GeneRuler Low Range DNA Ladder (Thermo Scientific, Waltham, MA, USA). Lanes 2–11: PCR products of sRNAs.
Genes 09 00083 g003
Table 1. Preliminary library of 295 sRNAs in H. mediterranei obtained by a bioinformatic approach based on candidate small RNAs of other species.
Table 1. Preliminary library of 295 sRNAs in H. mediterranei obtained by a bioinformatic approach based on candidate small RNAs of other species.
Reference SpeciesNumber of Candidate sRNAsIdentity (%)Number of MismatchesNumber of GapsE-Value
H. volcanii11674.19–100<30<11<0.05
M. mazei5582.86–100<5<2<1
S. solfataricus4795.83–100<3<3<1
A. fulgidus7473.28–100<30<4<1
Table 2. Compilation of positive predicted sRNAs in H. mediterranei using RNA-Sequencing.
Table 2. Compilation of positive predicted sRNAs in H. mediterranei using RNA-Sequencing.
sRNA NameSize sRNA (pb)Position (Start-Stop)Localisation H. MedStrandClassificationGene Environment
HM1_V1312947735-2947604CHR-IntergenicHFX_3032 (gbp3)/HFX_0001(cdc6A-1)
HM3_V64602004-601940CHR-IntergenicHFX_0639/HFX_640 (aspS)
HM4_V1432374110-2374253CHR-IntergenicHFX_2271(htr18-1)/HFX_2270
HM5_V46576412-576366CHR-IntergenicHFX_0609(cdc6C)/HFX_0610
HM6_V117401997-401880CHR-IntergenicHFX_0430(ilvA-1)/HFX_0431
HM8_V18310608-310590CHR-IntergenicHFX_0331(tRNA)/HFX_0332 (rpoH)
HM10_V421913258-1913216CHR-IntergenicHFX_2748(rpl7AE)/HFX_2747
HM12_V1301928092-1927962CHR-IntergenicHFX_2734/HFX_2733 (lsm1)
HM15_V1132841808-2841695CHR-IntergenicHFX_1820(rRNA)/HFX_2934
HM17_V81136329-136248CHR-IntergenicHFX_0140 (tRNA)/HFX_0141
HM19_V1112175910-2175799CHR-IntergenicHFX_2478/HFX_2477(tRNA)
HM20_V1041737798-1737694CHR-IntergenicHFX_1813/HFX_1814
HM24_V96440309-440213CHR-IntergenicHFX_0475/HFX_0476
HM25_V902584538-2584448CHR-IntergenicHFX_2084(ydjP-1)/HFX_2083
HM26_V1021737680-1737578CHR-IntergenicHFX_1813/HFX_1814
HM28_V771298728-1298651CHR-IntergenicHFX_1375(aroA)/HFX_1377
HM30_V160248344-248184CHR-IntergenicHFX_0263/HFX_0264
HM35_V1682385624-2385456CHR-AntisenseHFX_2256
HM36_V26227244-26982CHR-IntergenicHFX_0025/HFX_0026
HM38_V912836827-2836736CHR-IntergenicHFX_1823/HFX_1822
HM39_V902685508-2685418CHR-AntisenseHFX_1980 (abc22A)
HM42_V581769424-1769366CHR-IntergenicHFX_2905/HFX_2904(folP)
HM46_V2002756370-2756170CHR-IntergenicHFX_1903/HFX_1902
HM47_V1241816376-1816252CHR-IntergenicHFX_2862/HFX_2861
HM49_V38359173-359135pHM500-IntergenicHFX_6336/HFX_6337
HM50_V781116972-1116894CHR-IntergenicHFX_1173 (tRNA)/HFX_1174(tatAE)
HM52_V4572020941-2020484CHR-IntergenicHFX_2638(dkgB)/HFX_2637(pepC)
HM53_V156397174-397018CHR-IntergenicHFX_0425/HFX_0426
HM54_V722591733-2591661CHR-IntergenicHFX_2076/2075 (tnp4)
HM55_V33486185-486152CHR-IntergenicHFX_0525(mutS)/HFX_0526(livK)
HM56_V1702852306-2852136CHR-IntergenicHFX_2941(pheS)/HFX_2942
HM2_V52681041-681093CHR+IntergenicHFX_0721/HFX_0722
HM7_V24310803-310827CHR+IntergenicHFX_0331 (tRNA)/HFX_0332 (rpoH)
HM9_V271913243-1913270CHR+IntergenicHFX_2748(rpl7AE)/HFX_2747
HM11_V291927979-1928008CHR+IntergenicHFX_2734/HFX_2733 (lsm1)
HM13_V721930742-1930814CHR+IntergenicHFX_2731(purF)/HFX_2730
HM14_V491931964-1932013CHR+IntergenicHFX_2729/HFX_2728
HM16_V246136148-136394CHR+IntergenicHFX_0140 (tRNA)/HFX_0141
HM18_V1142175756-2175870CHR+IntergenicHFX_2478/HFX_2477(tRNA)
HM21_V2627000891-27000917CHR+IntergenicHFX_1962/HFX_1961
HM22_V212249860-2249881CHR+IntergenicHFX_2397 (tRNA)/HFX_2396
HM23_V111257489-257600CHR+IntergenicHFX_0274/HFX_0275 (tRNA)
HM27_V147581307-581454CHR+IntergenicHFX_0614(tRNA)/HFX_0615
HM29_V258248240-248498CHR+IntergenicHFX_0263/HFX_0264
HM31_V1251824730-1824855CHR+IntergenicHFX_2852/HFX_2851
HM32_V1022342182-2342284CHR+IntergenicHFX_2304/HFX_2303
HM33_V161221020-221181CHR+AntisenseHFX_0231 (tfb1-1)
HM34_V951259873-1259968CHR+crRNAHFX_1335/HFX_1336
HM37_V14995312-95461CHR+AntisenseHFX_0090 (HemL)
HM40_V351734001-1734036CHR+AntisenseHFX_1808 (ygcJ)
HM41_V532246499-2246552CHR+AntisenseHFX_2401 (xnuC-1)
HM43_V1001104687-1104787CHR+IntergenicHFX_1163 (tRNA)/HFX_1164
HM44_V38979793-979831CHR+IntergenicHFX_1024/HFX_1026
HM45_V90809546-809636CHR+IntergenicHFX_0853/HFX_0854
HM48_V195359045-359240pHM500+IntergenicHFX_6336/HFX_6337
HM51_V4572020484-2020941CHR+IntergenicHFX_2638(dkgB)/HFX_2637(pepC)
HM54_V351299991-1300026CHR+IntergenicHFX_1375(aroA)/HFX_1377
HM55_V441300050-1300094CHR+IntergenicHFX_1375(aroA)/HFX_1377
HM1_M922342183-2342275CHR+IntergenicHFX_2304/HFX_2303 (lactoyglutathione lyase)
HM2_M3456050-56016pHM300-IntergenicHFX_5046 (cyp1_citocromo P450)/ HFX_5047
HM3_M2922357250-2357542CHR+AntisenseHFX_2287 (selR)
HM4_M35201705-201740pHM500+AntisenseHFX_6177 (nahC2_Antiporter Na/H)
HM5_M1814274-14256pHM500-IntergenicHFX_6013 (nrd-1)/HFX_6014 (sodA)
HM6_M51338357-338306pHM500-AntisenseHFX_6316 (csh2-1_CRISPR-associated csh2 family protein)
HM7_M832448865-2448948CHR+IntergenicHFX_2200/2199 (imd3)
HM8_M176804495-804671CHR+IntergenicHFX_0847/HFX_0848 (epf2-1)
HM9_M531087812-1087865CHR+AntisenseHFX_1142
HM10_M971087972-1087875CHR-IntergenicHFX_1142/1143
HM11_M257709963-709706CHR-IntergenicHFX_0754 (prot.membr)/HFX_0755
HM12_M78862151-862073CHR-IntergenicHFX_0905/HFX_0906
HM13_M411287240-1287281CHR+IntergenicHFX_1366/HFX_1367 (htlD-1)
HM14_M411287281-1287240CHR-IntergenicHFX_1366/HFX_1367 (htlD-1)
HM15_M1532508794-2508641CHR-AntisenseHFX_2148 (polysaccharide biosynthesis transporter)
HM16_M572661554-2661497CHR-IntergenicHFX_2006 (paaJ-1)/HFX_2005 (nasD)
HM1_S2082948115-2948323CHR+IntergenicHFX_3032 (gbp3-1, GTP-binding protein)/HFX_0001 (cell division control protein)
HM2_S80412151-412231CHR+IntergenicHFX_0441 (pyrrolo-quinoline quinone)/HFX_0442 (serine/treonine protein kinase)
HM3_S981691336-1691238CHR-AntisenseHFX_1779 (integrase-like protein)
HM4_S81310404-310323CHR-IntergenicHFX_0329/HFX_0330 (tRNA)
HM5_S311104534-1104503CHR-AntisenseHFX_1161 (tRNA)
HM6_S60547749-547689CHR-AntisenseHFX_0587 (ark-6-1, putative signal-transducing histadine kinase)
HM7_S1852620188-2620373CHR+IntergenicHFX_2046/HFX_2045
HM8_S481432982-1432982CHR-AntisenseHFX_1518 (gdhA1, glutamate dehydrogenase NAD (P))
HM1_A14199868-100009CHR+AntisenseHM_0095 (amt1-1, ammonium transporter)
HM2_A482792936-793418CHR+AntisenseHFX_0839
HM3_A481277917-1277965CHR+IntergenicHFX_1356/HFX_1357
HM4_A191778748-1778729CHR-IntergenicHFX_2897 (fumC-1, fumarate hydratase)/HFX_2896 (carbohydrate ATP-binding cassette (ABC) transporter substrate-binding protein)
HM5_A222113156-2113178CHR+IntergenicHFX_2528/HFX_2529 (gpm, phosphoglycerate mutase)
HM6_A1272520322-2520449CHR+IntergenicHFX_2140 (cdc6H, cell division control protein cdc 6-like protein)/HFX_2139(galE, nucleoside-diphosphate sugar epimerase)
Table 3. Analysis of possible targets of antisense sRNAs using IntaRNA [38].
Table 3. Analysis of possible targets of antisense sRNAs using IntaRNA [38].
sRNAsRNA Size (pb)Gene TargetRNA-Seq Ammonium (+/−)RNA-Seq Nitrate (+/−)Target PositionsRNA PositionEnergy (Kcal/mol)
HM33_V161tfb1-1++513–6621–150−247.734
HM35_V168HFX_2256+414–46810–88−19.201
HM37_V149HemL++971–11201–149−238.470
HM39_V90abc22A++375–41228–62−15.561
HM40_V35ygcJ++260–2951–35−55.516
HM41_V53xnuC-1++583–6361–53−86.827
HM3_M292selR++316–423186–292−162.400
HM4_M35nahC2, antiporter Na+/H+++1151–11861–35−55.151
HM6_M51csh2-1_CRISPR-associated csh2 family protein++
HM9_M53HFX_1142++1–241–24−33.421
HM15_M153Polysaccharide biosynthesis transporter++93–110115–129−14.191
HM1_A141amt1-1, ammonium transporter+1330–14711–141−229.279
HM2_A482HFX_0839++513–662301–450−245.770
HM3_S98Integrase-like protein++17–2921–33−9.201
HM5_S31tRNA++47–5425–32−1.222
HM6_S60ark-6-1, putative signal-transducing histadine kinase++1546–155626–37−12.128
HM8_S48gdhA1, glutamate dehydrogenase NAD (P)+230–2362–8−6.608
Table 4. Summary of the characteristics of sRNAs with differential expression (nitrate/ammonium).
Table 4. Summary of the characteristics of sRNAs with differential expression (nitrate/ammonium).
sRNA NameLocalisationStrandPositionSize (nt)ClassificationGene Target 1log2 FoldChangep-Valuep-adj
HM16_MCHR-2661497-266155457Intergenic HFX_1537, pgk, HFX_1505, trkA, HFX_17898.69855.05 × 10−2598.09 × 10−258
HM7_SCHR+2620188-2620373185Intergenic HFX_0627, acnA, HFX_1005, arsR, HFX_30061.75471.29 × 10−241.03 × 10−23
HM52_VCHR-2020941-2020484457Intergenic hflC, HFX_0589, HFX_03551.36331.06 × 10−215.96 × 10−20
HM46_VCHR-2756370-2756170200Intergenic HFX_0329, psmA, HFX_2148, HFX_1154, graD51.19662.59 × 10−67.25 × 10−5
HM39_VCHR-2685508-268541890Antisense ilvB, cbs_4, HFX_0847, gldA, ligA0.84980.0060.0413
HM54_VCHR+2591733-259166135Intergenic HFX_0906, HFX_1268, HFX_0562, arsR, nfi0.67080.0040.0338
HM12_VCHR-1928092-1927962130Intergenic HFX_1575, mutT, HFX_2074, coxB40.62630.0030.0333
HM1_SCHR+2948115-2948323208IntergeniclivG, HFX_1497, HFX_2088, HFX_2739, nce10.51920.00030.00086
HM37_VCHR+95312-95461149AntisenseHFX_0565, nifS, HFX_1510, HFX_0591−0.47830.00110.0158
HM36_VCHR-27244-26982262Intergenic HFX_2059, rpoF, apa−0.48450.00360.0332
HM38_VCHR-2836827-283673691IntergenicmenD, hcpE, HFX_1197, HFX_0402−0.61200.0080.0451
HM3_MCHR+2357250-2357542292AntisenseHFX_2088, rnhB, petE, HFX_2672, HFX_2025−0.62470.01120.0894
HM8_VCHR-310608-31059018Intergenic no targets−0.73870.00790.0451
HM6_SCHR-547689-54774960AntisensexthA, dap2, HFX_0771, HFX_0201, yfmJ1−0.83730.00030.00086
HM8_SCHR-1432982-143303048Antisenseatpl, HFX_0366, gatD, atpF, gvpJ−0.87910.00210.00413
HM9_VCHR+1913243-191327027IntergenicppiB, HFX_1294, HFX_2540, ispA, HFX_2807−1.32530.00070.0128
1 Values of free energy and p-value are available in Table S8.

Share and Cite

MDPI and ACS Style

Payá, G.; Bautista, V.; Camacho, M.; Castejón-Fernández, N.; Alcaraz, L.A.; Bonete, M.-J.; Esclapez, J. Small RNAs of Haloferax mediterranei: Identification and Potential Involvement in Nitrogen Metabolism. Genes 2018, 9, 83. https://doi.org/10.3390/genes9020083

AMA Style

Payá G, Bautista V, Camacho M, Castejón-Fernández N, Alcaraz LA, Bonete M-J, Esclapez J. Small RNAs of Haloferax mediterranei: Identification and Potential Involvement in Nitrogen Metabolism. Genes. 2018; 9(2):83. https://doi.org/10.3390/genes9020083

Chicago/Turabian Style

Payá, Gloria, Vanesa Bautista, Mónica Camacho, Natalia Castejón-Fernández, Luís A. Alcaraz, María-José Bonete, and Julia Esclapez. 2018. "Small RNAs of Haloferax mediterranei: Identification and Potential Involvement in Nitrogen Metabolism" Genes 9, no. 2: 83. https://doi.org/10.3390/genes9020083

APA Style

Payá, G., Bautista, V., Camacho, M., Castejón-Fernández, N., Alcaraz, L. A., Bonete, M. -J., & Esclapez, J. (2018). Small RNAs of Haloferax mediterranei: Identification and Potential Involvement in Nitrogen Metabolism. Genes, 9(2), 83. https://doi.org/10.3390/genes9020083

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop