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

A conserved transcription factor controls gluconeogenesis via distinct targets in hypersaline-adapted archaea with diverse metabolic capabilities

Timely regulation of carbon metabolic pathways is essential for cellular processes and to prevent futile cycling of intracellular metabolites. In Halobacterium salinarum, a hypersaline adapted archaeon, a sugar-sensing TrmB family protein controls gluconeogenesis and other biosynthetic pathways. Notably, Hbt. salinarum does not utilize carbohydrates for energy, uncommon among Haloarchaea. We characterized a TrmB-family transcriptional regulator in a saccharolytic generalist, Haloarcula hispanica, to investigate whether the targets and function of TrmB, or its regulon, is conserved in related species with distinct metabolic capabilities. In Har. hispanica, TrmB binds to 15 sites in the genome and induces the expression of genes primarily involved in gluconeogenesis and tryptophan biosynthesis. An important regulatory control point in Hbt. salinarum, activation of ppsA and repression of pykA, is absent in Har. hispanica. Contrary to its role in Hbt. salinarum and saccharolytic hyperthermophiles, TrmB does not act as a global regulator: it does not directly repress the expression of glycolytic enzymes, peripheral pathways such as cofactor biosynthesis, or catabolism of other carbon sources in Har. hispanica. Cumulatively, these findings suggest rewiring of the TrmB regulon alongside metabolic network evolution in Haloarchaea.

Genome-wide transcriptional response to silver stress in extremely halophilic archaeon Haloferax alexandrinus DSM 27206 T

BACKGROUND

The extremely halophilic archaeon Haloferax (Hfx.) alexandrinus DSM 27206 T was previously documented for the ability to biosynthesize silver nanoparticles while mechanisms underlying its silver tolerance were overlooked. In the current study, we aimed to assess the transcriptional response of this haloarchaeon to varying concentrations of silver, seeking a comprehensive understanding of the molecular determinants underpinning its heavy metal tolerance.

RESULTS

The growth curves confirmed the capacity of Hfx. alexandrinus to surmount silver stress, while the SEM-EDS analysis illustrated the presence of silver nanoparticles in cultures exposed to 0.5 mM silver nitrate. The RNA-Seq based transcriptomic analysis of Hfx. alexandrinus cells exposed to 0.1, 0.25, and 0.5 mM silver nitrate revealed the differential expression of multiple sets of genes potentially employed in heavy-metal stress response, genes mostly related to metal transporters, basic metabolism, oxidative stress response and cellular motility. The RT-qPCR analysis of selected transcripts was conducted to verify and validate the generated RNA-Seq data.

CONCLUSIONS

Our results indicated that copA, encoding the copper ATPase, is essential for the survival of Hfx. alexandrinus cells in silver-containing saline media. The silver-exposed cultures underwent several metabolic adjustments that enabled the activation of enzymes involved in the oxidative stress response and impairment of the cellular movement capacity. To our knowledge, this study represents the first comprehensive analysis of gene expression in halophillic archaea facing increased levels of heavy metals.

Comparative Analysis of rRNA Removal Methods for RNA-Seq Differential Expression in Halophilic Archaea

Despite intense recent research interest in archaea, the scientific community has experienced a bottleneck in the study of genome-scale gene expression experiments by RNA-seq due to the lack of commercial and specifically designed rRNA depletion kits. The high rRNA:mRNA ratio (80–90%: ~10%) in prokaryotes hampers global transcriptomic analysis. Insufficient ribodepletion results in low sequence coverage of mRNA, and therefore, requires a substantially higher number of replicate samples and/or sequencing reads to achieve statistically reliable conclusions regarding the significance of differential gene expression between case and control samples. Here, we show that after the discontinuation of the previous version of RiboZero (Illumina, San Diego, CA, USA) that was useful in partially or completely depleting rRNA from archaea, archaeal transcriptomics studies have experienced a slowdown. To overcome this limitation, here, we analyze the efficiency for four different hybridization-based kits from three different commercial suppliers, each with two sets of sequence-specific probes to remove rRNA from four different species of halophilic archaea. We conclude that the key for transcriptomic success with the currently available tools is the probe-specificity for the rRNA sequence hybridization. With this paper, we provide insights into the archaeal community for selecting certain reagents and strategies over others depending on the archaeal species of interest. These methods yield improved RNA-seq sensitivity and enhanced detection of low abundance transcripts.

Post-transcriptional regulation of redox homeostasis by the small RNA SHOxi in haloarchaea

While haloarchaea are highly resistant to oxidative stress, a comprehensive understanding of the processes regulating this remarkable response is lacking. Oxidative stress-responsive small non-coding RNAs (sRNAs) have been reported in the model archaeon, Haloferax volc anii, but targets and mechanisms have not been elucidated. Using a combination of high throughput and reverse molecular genetic approaches, we elucidated the functional role of the most up-regulated intergenic sRNA during oxidative stress in H. volcanii, named Small RNA in Haloferax Oxidative Stress (SHOxi). SHOxi was predicted to form a stable secondary structure with a conserved stem-loop region as the potential binding site for trans-targets. NAD-dependent malic enzyme mRNA, identified as a putative target of SHOxi, interacted directly with a putative 'seed' region within the predicted stem loop of SHOxi. Malic enzyme catalyzes the oxidative decarboxylation of malate into pyruvate using NAD+ as a cofactor. The destabilization of malic enzyme mRNA, and the decrease in the NAD+/NADH ratio, resulting from the direct RNA-RNA interaction between SHOxi and its trans-target was essential for the survival of H. volcanii to oxidative stress. These findings indicate that SHOxi likely regulates redox homoeostasis during oxidative stress by the post-transcriptional destabilization of malic enzyme mRNA. SHOxi-mediated regulation provides evidence that the fine-tuning of metabolic cofactors could be a core strategy to mitigate damage from oxidative stress and confer resistance. This study is the first to establish the regulatory effects of sRNAs on mRNAs during the oxidative stress response in Archaea.

Novel Glutamate-Putrescine Ligase Activity in Haloferax mediterranei: A New Function for glnA-2 Gene

The genome of the halophilic archaea Haloferax mediterranei contains three ORFs that show homology with glutamine synthetase (GS) (glnA-1, glnA-2, and glnA-3). Previous studies have focused on the role of GlnA-1, suggesting that proteins GlnA-2 and GlnA-3 could play a different role to that of GS. Glutamine synthetase (EC 6.3.1.2) belongs to the class of ligases, including 20 subclasses of other different enzymes, such as aspartate–ammonia ligase (EC 6.3.1.1), glutamate–ethylamine ligase (EC 6.3.1.6), and glutamate–putrescine ligase (EC 6.3.1.11). The reaction catalyzed by glutamate–putrescine ligase is comparable to the reaction catalyzed by glutamine synthetase (GS). Both enzymes can bind a glutamate molecule to an amino group: ammonium (GS) or putrescine (glutamate–putrescine ligase). In addition, they present the characteristic catalytic domain of GS, showing significant similarities in their structure. Although these proteins are annotated as GS, the bioinformatics and experimental results obtained in this work indicate that the GlnA-2 protein (HFX_1688) is a glutamate–putrescine ligase, involved in polyamine catabolism. The most significant results are those related to glutamate–putrescine ligase’s activity and the analysis of the transcriptional and translational expression of the glnA-2 gene in the presence of different nitrogen sources. This work confirms a new metabolic pathway in the Archaea domain which extends the knowledge regarding the utilization of alternative nitrogen sources in this domain.

Functional analysis of Lsm protein under multiple stress conditions in the extreme haloarchaeon Haloferax mediterranei

The Sm, like-Sm, and Hfq proteins belonging to the Sm superfamily of proteins are represented in all domains of life. These proteins are involved in several RNA metabolism pathways. The functions of bacterial Hfq and eukaryotic Sm proteins have been described, but knowledge about the in vivo functions of archaeal Sm proteins remains limited. This study aims to improve the understanding of Lsm proteins and their role using the haloarchaeon Haloferax mediterranei as a model microorganism. The Haloferax mediterranei genome contains one lsm gene that overlaps with the rpl37e gene. To determine the expression of lsm and rpl37e genes and the co-transcription of both, reverse transcription-polymerase chain reaction (RT-PCR) analyses were performed under different standard and stress conditions. The results suggest that the expression of lsm and rpl37e is constitutive. Co-transcription occurs at sub-optimal salt concentrations and temperatures, depending on the growth phase. The halophilic Lsm protein contains two Sm motifs, Sm1 and Sm2, and the sequence encoding the Sm2 motif also constitutes the promoter of the rpl37e gene. To investigate their biological functions, the lsm deletion mutant and the Sm1 motif deletion mutant, where the Sm2 motif remained intact, were generated and characterised. Comparison of the lsm deletion mutant, Sm1 deletion mutant, and the parental strain HM26 under standard and stress growth conditions revealed growth differences. Finally, swarming assays in complex and defined media showed greater swarming capacity in the deletion mutants.

A Small RNA Is Linking CRISPR-Cas and Zinc Transport

The function and mode of action of small regulatory RNAs is currently still understudied in archaea. In the halophilic archaeon Haloferax volcanii, a plethora of sRNAs have been identified; however, in-depth functional analysis is missing for most of them. We selected a small RNA (s479) from Haloferax volcanii for detailed characterization. The sRNA gene is encoded between a CRISPR RNA locus and the Cas protein gene cluster, and the s479 deletion strain is viable and was characterized in detail. Transcriptome studies of wild-type Haloferax cells and the deletion mutant revealed upregulation of six genes in the deletion strain, showing that this sRNA has a clearly defined function. Three of the six upregulated genes encode potential zinc transporter proteins (ZnuA1, ZnuB1, and ZnuC1) suggesting the involvement of s479 in the regulation of zinc transport. Upregulation of these genes in the deletion strain was confirmed by northern blot and proteome analyses. Furthermore, electrophoretic mobility shift assays demonstrate a direct interaction of s479 with the target znuC1 mRNA. Proteome comparison of wild-type and deletion strains further expanded the regulon of s479 deeply rooting this sRNA within the metabolism of H. volcanii especially the regulation of transporter abundance. Interestingly, s479 is not only encoded next to CRISPR-cas genes, but the mature s479 contains a crRNA-like 5' handle, and experiments with Cas protein deletion strains indicate maturation by Cas6 and interaction with Cas proteins. Together, this might suggest that the CRISPR-Cas system is involved in s479 function.

Towards the Elucidation of Assimilative nasABC Operon Transcriptional Regulation in Haloferax mediterranei

The assimilatory pathway of the nitrogen cycle in the haloarchaeon Haloferax mediterranei has been well described and characterized in previous studies. However, the regulatory mechanisms involved in the gene expression of this pathway remain unknown in haloarchaea. This work focuses on elucidating the regulation at the transcriptional level of the assimilative nasABC operon (HFX_2002 to HFX_2004) through different approaches. Characterization of its promoter region using β-galactosidase as a reporter gene and site-directed mutagenesis has allowed us to identify possible candidate binding regions for a transcriptional factor. The identification of a potential transcriptional regulator related to nitrogen metabolism has become a real challenge due to the lack of information on haloarchaea. The investigation of protein–DNA binding by streptavidin bead pull-down analysis combined with mass spectrometry resulted in the in vitro identification of a transcriptional regulator belonging to the Lrp/AsnC family, which binds to the nasABC operon promoter (p.nasABC). To our knowledge, this study is the first report to suggest the AsnC transcriptional regulator as a powerful candidate to play a regulatory role in nasABC gene expression in Hfx. mediterranei and, in general, in the assimilatory nitrogen pathway.

Gene Expression of Haloferax volcanii on Intermediate and Abundant Sources of Fixed Nitrogen

Haloferax volcanii, a well-developed model archaeon for genomic, transcriptomic, and proteomic analyses, can grow on a defined medium of abundant and intermediate levels of fixed nitrogen. Here we report a global profiling of gene expression of H. volcanii grown on ammonium as an abundant source of fixed nitrogen compared to l-alanine, the latter of which exemplifies an intermediate source of nitrogen that can be obtained from dead cells in natural habitats. By comparing the two growth conditions, 30 genes were found to be differentially expressed, including 16 genes associated with amino acid metabolism and transport. The gene expression profiles contributed to mapping ammonium and l-alanine usage with respect to transporters and metabolic pathways. In addition, conserved DNA motifs were identified in the putative promoter regions and transcription factors were found to be in synteny with the differentially expressed genes, leading us to propose regulons of transcriptionally co-regulated operons. This study provides insight to how H. volcanii responds to and utilizes intermediate vs. abundant sources of fixed nitrogen for growth, with implications for conserved functions in related halophilic archaea.

Global Transcriptional Programs in Archaea Share Features with the Eukaryotic Environmental Stress Response

The environmental stress response (ESR), a global transcriptional program originally identified in yeast, is characterized by a rapid and transient transcriptional response composed of large, oppositely regulated gene clusters. Genes induced during the ESR encode core components of stress tolerance, macromolecular repair, and maintenance of homeostasis. In this review, we investigate the possibility for conservation of the ESR across the eukaryotic and archaeal domains of life. We first re-analyze existing transcriptomics data sets to illustrate that a similar transcriptional response is identifiable in Halobacterium salinarum, an archaeal model organism. To substantiate the archaeal ESR, we calculated gene-by-gene correlations, gene function enrichment, and comparison of temporal dynamics. We note reported examples of variation in the ESR across fungi, then synthesize high-level trends present in expression data of other archaeal species. In particular, we emphasize the need for additional high-throughput time series expression data to further characterize stress-responsive transcriptional programs in the Archaea. Together, this review explores an open question regarding features of global transcriptional stress response programs shared across domains of life.

Whole genome sequence of moderate halophilic marine bacterium Marinobacter litoralis SW-45: Abundance of non-coding RNAs

A report on the de novo Whole Genome Sequence (WGS) of Marinobacter litoralis SW-45, a moderately salt-tolerant bacterium isolated from the seawater in Malaysia is presented. The strain has a genome size of 3.45 Mb and is capable of producing halophilic lipase, protease and esterase enzymes. Computational prediction of non-coding RNA (ncRNA) genes in M. litoralis SW-45 was performed using standalone software known as the non-coding RNA characterization (nocoRNAc). In addition, a phylogenetic tree showing the evolutionary relationship between the strain and other members of the genus Marinobacter was constructed using 16SrRNA sequence information. A total of 385 ncRNA transcripts, 1124 terminator region, and 2350 Stress Induced Duplex Destabilization sites were predicted. The current WGS shotgun project has provided the relevant genetic information that may be useful for the strain's improvement studies. This manuscript gives the first description of M. litoralis with a complete genome.

Putative small RNAs controlling detoxification of industrial cyanide-containing wastewaters by Pseudomonas pseudoalcaligenes CECT5344

The alkaliphilic bacterium Pseudomonas pseudoalcaligenes CECT5344 uses free cyanide and several metal−cyanide complexes as the sole nitrogen source and tolerates high concentrations of metals like copper, zinc and iron, which are present in the jewelry wastewaters. To understand deeply the regulatory mechanisms involved in the transcriptional regulation of cyanide-containing wastewaters detoxification by P. pseudoalcaligenes CECT5344, RNA-Seq has been performed from cells cultured with a cyanide-containing jewelry wastewater, sodium cyanide or ammonium chloride as the sole nitrogen source. Small RNAs (sRNAs) that may have potential regulatory functions under cyanotrophic conditions were identified. In total 20 sRNAs were identified to be differentially expressed when compared the jewelry residue versus ammonium as nitrogen source, 16 of which could be amplified successfully by RT-PCR. As predicted targets of these 16 sRNAs were several components of the nit1C gene cluster encoding the nitrilase NitC essential for cyanide assimilation, the cioAB gene cluster that codes for the cyanide-insensitive cytochrome bd-type terminal oxidase, the medium length-polyhydroxyalkanoates (ml-PHAs) gene cluster, and gene clusters related with a global nitrogen limitation response like those coding for glutamine synthase and urease. Other targets were non-clustered genes (or their products) involved in metal resistance and iron acquisition, such as metal extrusion systems and the ferric uptake regulatory (Fur) protein, and a GntR-like regulatory family member probably involved in the regulation of the cyanide assimilation process in the strain CECT5344. Induction of genes targeted by sRNAs in the jewelry residue was demonstrated by qRT-PCR.

Characterization of the Microbial Population Inhabiting a Solar Saltern Pond of the Odiel Marshlands (SW Spain)

The solar salterns located in the Odiel marshlands, in southwest Spain, are an excellent example of a hypersaline environment inhabited by microbial populations specialized in thriving under conditions of high salinity, which remains poorly explored. Traditional culture-dependent taxonomic studies have usually under-estimated the biodiversity in saline environments due to the difficulties that many of these species have to grow at laboratory conditions. Here we compare two molecular methods to profile the microbial population present in the Odiel saltern hypersaline water ponds (33% salinity). On the one hand, the construction and characterization of two clone PCR amplified-16S rRNA libraries, and on the other, a high throughput 16S rRNA sequencing approach based on the Illumina MiSeq platform. The results reveal that both methods are comparable for the estimation of major genera, although massive sequencing provides more information about the less abundant ones. The obtained data indicate that Salinibacter ruber is the most abundant genus, followed by the archaea genera, Halorubrum and Haloquadratum. However, more than 100 additional species can be detected by Next Generation Sequencing (NGS). In addition, a preliminary study to test the biotechnological applications of this microbial population, based on its ability to produce and excrete haloenzymes, is shown.