Primary transcriptome map of the hyperthermophilic archaeon Thermococcus kodakarensis

Evolutionary mining and functional characterization of TnpB nucleases identify efficient miniature genome editors

As the evolutionary ancestor of Cas12 nuclease, the transposon (IS200/IS605)-encoded TnpB proteins act as compact RNA-guided DNA endonucleases. To explore their evolutionary diversity and potential as genome editors, we screened TnpBs from 64 annotated IS605 members and identified 25 active in Escherichia coli, of which three are active in human cells. Further characterization of these 25 TnpBs enables prediction of the transposon-associated motif (TAM) and the right-end element RNA (reRNA) directly from genomic sequences. We established a framework for annotating TnpB systems in prokaryotic genomes and applied it to identify 14 additional candidates. Among these, ISAam1 (369 amino acids (aa)) and ISYmu1 (382 aa) TnpBs demonstrated robust editing activity across dozens of genomic loci in human cells. Both RNA-guided genome editors demonstrated similar editing efficiency as SaCas9 (1,053 aa) while being substantially smaller. The enormous diversity of TnpBs holds potential for the discovery of additional valuable genome editors.

Transposons and CRISPR: Rewiring Gene Editing

CRISPR-Cas is driving a gene editing revolution because of its simple reprogramming. However, off-target effects and dependence on the double-strand break repair pathways impose important limitations. Because homology-directed repair acts primarily in actively dividing cells, many of the current gene correction/replacement approaches are restricted to a minority of cell types. Furthermore, current approaches display low efficiency upon insertion of large DNA cargos (e.g., sequences containing multiple gene circuits with tunable functionalities). Recent research has revealed new links between CRISPR-Cas systems and transposons providing new scaffolds that might overcome some of these limitations. Here, we comment on two new transposon-associated RNA-guided mechanisms considering their potential as new gene editing solutions. Initially, we focus on a group of small RNA-guided endonucleases of the IS200/IS605 family of transposons, which likely evolved into class 2 CRISPR effector nucleases (Cas9s and Cas12s). We explore the diversity of these nucleases (named OMEGA, obligate mobile element-guided activity) and analyze their similarities with class 2 gene editors. OMEGA nucleases can perform gene editing in human cells and constitute promising candidates for the design of new compact RNA-guided platforms. Then, we address the co-option of the RNA-guided activity of different CRISPR effector nucleases by a specialized group of Tn7-like transposons to target transposon integration. We describe the various mechanisms used by these RNA-guided transposons for target site selection and integration. Finally, we assess the potential of these new systems to circumvent some of the current gene editing challenges.

High-resolution temporal profiling of E. coli transcriptional response

Understanding how cells dynamically adapt to their environment is a primary focus of biology research. Temporal information about cellular behavior is often limited by both small numbers of data time-points and the methods used to analyze this data. Here, we apply unsupervised machine learning to a data set containing the activity of 1805 native promoters in E. coli measured every 10 minutes in a high-throughput microfluidic device via fluorescence time-lapse microscopy. Specifically, this data set reveals E. coli transcriptome dynamics when exposed to different heavy metal ions. We use a bioinformatics pipeline based on Independent Component Analysis (ICA) to generate insights and hypotheses from this data. We discovered three primary, time-dependent stages of promoter activation to heavy metal stress (fast, intermediate, and steady). Furthermore, we uncovered a global strategy E. coli uses to reallocate resources from stress-related promoters to growth-related promoters following exposure to heavy metal stress.

A robust genetic toolbox for fine-tuning gene expression in the CO2-Fixing methanogenic archaeon Methanococcus maripaludis

Libraries of well-characterized genetic elements for fine-tuning gene expression are essential for biological and biotechnological research and applications. The fast-growing and genetically tractable methanogen, Methanococcus maripaludis, is a promising host organism for biotechnological conversion of carbon dioxide and renewable hydrogen into fuels and value-added products, as well as fundamental biological studies of archaea. However, the lack of molecular tools for gene expression has hindered its application as a workhorse to fine-tune gene and metabolic pathway expressions. In this study, we developed a genetic toolbox, including libraries of promoters, ribosome binding sites (RBS), and neutral sites for chromosomal integration, to facilitate precise gene expression in M. maripaludis. We generated a promoter library consisting of 81 constitutive promoters with expression strengths spanning a ∼104-fold dynamic range. Importantly, we identified a base composition rule for strong archaeal promoters and successfully remodeled weak promoters, enhancing their activities by up to 120-fold. We also established an RBS library containing 42 diverse RBS sequences with translation strengths covering a ∼100-fold dynamic range. Additionally, we identified eight neutral sites and developed a one-step, Cas9-based marker-less knock-in approach for chromosomal integration. We successfully applied the characterized promoter and RBS elements to significantly improve recombinant protein expression by 41-fold and modulate essential gene expression to generate corresponding physiological changes in M. maripaludis. Therefore, this work establishes a solid foundation for utilizing this autotrophic methanogen as an ideal workhorse for archaeal biology and biotechnological studies and applications.

Internal transcription termination widely regulates differential expression of operon-organized genes including ribosomal protein and RNA polymerase genes in an archaeon

Genes organized within operons in prokaryotes benefit from coordinated expression. However, within many operons, genes are expressed at different levels, and the mechanisms for this remain obscure. By integrating PacBio-seq, dRNA-seq, Term-seq and Illumina-seq data of a representative archaeon Methanococcus maripaludis, internal transcription termination sites (ioTTSs) were identified within 38% of operons. Higher transcript and protein abundances were found for genes upstream than downstream of ioTTSs. For representative operons, these differences were confirmed by northern blotting, qRT-PCR and western blotting, demonstrating that these ioTTS terminations were functional. Of special interest, mutation of ioTTSs in ribosomal protein (RP)-RNA polymerase (RNAP) operons not only elevated expression of the downstream RNAP genes but also decreased production of the assembled RNAP complex, slowed whole cell transcription and translation, and inhibited growth. Overexpression of the RNAP subunits with a shuttle vector generated the similar physiological effects. Therefore, ioTTS termination is a general and physiologically significant regulatory mechanism of the operon gene expression. Because the RP-RNAP operons are found to be widely distributed in archaeal species, this regulatory mechanism could be commonly employed in archaea.

Molecular dynamics simulation-based trinucleotide and tetranucleotide level structural and energy characterization of the functional units of genomic DNA

Genomes of most organisms on earth are written in a universal language of life, made up of four units - adenine (A), thymine (T), guanine (G), and cytosine (C), and understanding the way they are put together has been a great challenge to date. Multiple efforts have been made to annotate this wonderfully engineered string of DNA using different methods but they lack a universal character. In this article, we have investigated the structural and energetic profiles of both prokaryotes and eukaryotes by considering two essential genomic sites, viz., the transcription start sites (TSS) and exon-intron boundaries. We have characterized these sites by mapping the structural and energy features of DNA obtained from molecular dynamics simulations, which considers all possible trinucleotide and tetranucleotide steps. For DNA, these physicochemical properties show distinct signatures at the TSS and intron-exon boundaries. Our results firmly convey the idea that DNA uses the same dialect for prokaryotes and eukaryotes and that it is worth going beyond sequence-level analyses to physicochemical space to determine the functional destiny of DNA sequences.

Histone variants in archaea - An undiscovered country

Exchanging core histones in the nucleosome for paralogous variants can have important functional ramifications. Many of these variants, and their physiological roles, have been characterized in exquisite detail in model eukaryotes, including humans. In comparison, our knowledge of histone biology in archaea remains rudimentary. This is true in particular for our knowledge of histone variants. Many archaea encode several histone genes that differ in sequence, but do these paralogs make distinct, adaptive contributions to genome organization and regulation in a manner comparable to eukaryotes? Below, we review what we know about histone variants in archaea at the level of structure, regulation, and evolution. In all areas, our knowledge pales when compared to the wealth of insight that has been gathered for eukaryotes. Recent findings, however, provide tantalizing glimpses into a rich and largely undiscovered country that is at times familiar and eukaryote-like and at times strange and uniquely archaeal. We sketch a preliminary roadmap for further exploration of this country; an undertaking that may ultimately shed light not only on chromatin biology in archaea but also on the origin of histone-based chromatin in eukaryotes.

Endonucleolytic processing plays a critical role in the maturation of ribosomal RNA in Methanococcus maripaludis

Ribosomal RNA (rRNA) processing and maturation are fundamentally important for ribosome biogenesis, but the mechanisms in archaea, the third form of life, remains largely elusive. This study aimed to investigate the rRNA maturation process in Methanococcus maripaludis, a representative archaeon lacking known 3'-5' exonucleases. Through cleavage site identification and enzymatic assays, the splicing endonuclease EndA was determined to process the bulge-helix-bulge (BHB) motifs in 16S and 23S rRNA precursors. After splicing, the circular processing intermediates were formed and this was confirmed by quantitative RT-PCR and Northern blot. Ribonuclease assay revealed a specific cleavage at a 10-nt A/U-rich motif at the mature 5' end of pre-16S rRNA, which linearized circular pre-16S rRNA intermediate. Further 3'-RACE and ribonuclease assays determined that the endonuclease Nob1 cleaved the 3' extension of pre-16S rRNA, and so generated the mature 3' end. Circularized RT-PCR (cRT-PCR) and 5'-RACE identified two cleavage sites near helix 1 at the 5' end of 23S rRNA, indicating that an RNA structure-based endonucleolytic processing linearized the circular pre-23S rRNA intermediate. In the maturation of pre-5S rRNA, multiple endonucleolytic processing sites were determined at the 10-nt A/U-rich motif in the leader and trailer sequence. This study demonstrates that endonucleolytic processing, particularly at the 10-nt A/U-rich motifs play an essential role in the pre-rRNA maturation of M. maripaludis, indicating diverse pathways of rRNA maturation in archaeal species.

Machine learning and statistics shape a novel path in archaeal promoter annotation

Background

Archaea are a vast and unexplored domain. Bioinformatic techniques might enlighten the path to a higher quality genome annotation in varied organisms. Promoter sequences of archaea have the action of a plethora of proteins upon it. The conservation found in a structural level of the binding site of proteins such as TBP, TFB, and TFE aids RNAP-DNA stabilization and makes the archaeal promoter prone to be explored by statistical and machine learning techniques.

Results and discussions

In this study, experimentally verified promoter sequences of the organisms Haloferax volcanii, Sulfolobus solfataricus, and Thermococcus kodakarensis were converted into DNA duplex stability attributes (i.e. numerical variables) and were classified through Artificial Neural Networks and an in-house statistical method of classification, being tested with three forms of controls. The recognition of these promoters enabled its use to validate unannotated promoter sequences in other organisms. As a result, the binding site of basal transcription factors was located through a DNA duplex stability codification. Additionally, the classification presented satisfactory results (above 90%) among varied levels of control.

Concluding remarks

The classification models were employed to perform genomic annotation into the archaea Aciduliprofundum boonei and Thermofilum pendens, from which potential promoters have been identified and uploaded into public repositories.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12859-022-04714-x.

Genetic and Functional Analyses of Archaeal ATP-Dependent RNA Ligase in C/D Box sRNA Circularization and Ribosomal RNA Processing

RNA ligases play important roles in repairing and circularizing RNAs post-transcriptionally. In this study, we generated an allelic knockout of ATP-dependent RNA ligase (Rnl) in the hyperthermophilic archaeon Thermococcus kodakarensis to identify its biological targets. A comparative analysis of circular RNA reveals that the Rnl-knockout strain represses circularization of C/D box sRNAs without affecting the circularization of tRNA and rRNA processing intermediates. Recombinant archaeal Rnl could circularize C/D box sRNAs with a mutation in the conserved C/D box sequence element but not when the terminal stem structures were disrupted, suggesting that proximity of the two ends could be critical for intramolecular ligation. Furthermore, T. kodakarensis accumulates aberrant RNA fragments derived from ribosomal RNA in the absence of Rnl. These results suggest that Rnl is responsible for C/D box sRNA circularization and may also play a role in ribosomal RNA processing.

Deep conservation of histone variants in Thermococcales archaea

Histones are ubiquitous in eukaryotes where they assemble into nucleosomes, binding and wrapping DNA to form chromatin. One process to modify chromatin and regulate DNA accessibility is the replacement of histones in the nucleosome with paralogous variants. Histones are also present in archaea but whether and how histone variants contribute to the generation of different physiologically relevant chromatin states in these organisms remains largely unknown. Conservation of paralogs with distinct properties can provide prima facie evidence for defined functional roles. We recently revealed deep conservation of histone paralogs with different properties in the Methanobacteriales, but little is known experimentally about these histones. In contrast, the two histones of the model archaeon Thermococcus kodakarensis, HTkA and HTkB, have been examined in some depth, both in vitro and in vivo. HTkA and HTkB exhibit distinct DNA-binding behaviors and elicit unique transcriptional responses when deleted. Here, we consider the evolution of HTkA/B and their orthologs across the order Thermococcales. We find histones with signature HTkA- and HTkB-like properties to be present in almost all Thermococcales genomes. Phylogenetic analysis indicates the presence of one HTkA- and one HTkB-like histone in the ancestor of Thermococcales and long-term maintenance of these two paralogs throughout Thermococcales diversification. Our results support the notion that archaea and eukaryotes have convergently evolved histone variants that carry out distinct adaptive functions. Intriguingly, we also detect more highly diverged histone-fold proteins, related to those found in some bacteria, in several Thermococcales genomes. The functions of these bacteria-type histones remain unknown, but structural modeling suggests that they can form heterodimers with HTkA/B-like histones.

The Piezo-Hyperthermophilic Archaeon Thermococcus piezophilus Regulates Its Energy Efficiency System to Cope With Large Hydrostatic Pressure Variations

Deep-sea ecosystems share a common physical parameter, namely high hydrostatic pressure (HHP). Some of the microorganisms isolated at great depths have a high physiological plasticity to face pressure variations. The adaptive strategies by which deep-sea microorganisms cope with HHP variations remain to be elucidated, especially considering the extent of their biotopes on Earth. Herein, we investigated the gene expression patterns of Thermococcus piezophilus, a piezohyperthermophilic archaeon isolated from the deepest hydrothermal vent known to date, under sub-optimal, optimal and supra-optimal pressures (0.1, 50, and 90 MPa, respectively). At stressful pressures [sub-optimal (0.1 MPa) and supra-optimal (90 MPa) conditions], no classical stress response was observed. Instead, we observed an unexpected transcriptional modulation of more than a hundred gene clusters, under the putative control of the master transcriptional regulator SurR, some of which are described as being involved in energy metabolism. This suggests a fine-tuning effect of HHP on the SurR regulon. Pressure could act on gene regulation, in addition to modulating their expression.

The hyperthermophilic archaeon Thermococcus kodakarensis is resistant to pervasive negative supercoiling activity of DNA gyrase

In all cells, DNA topoisomerases dynamically regulate DNA supercoiling allowing essential DNA processes such as transcription and replication to occur. How this complex system emerged in the course of evolution is poorly understood. Intriguingly, a single horizontal gene transfer event led to the successful establishment of bacterial gyrase in Archaea, but its emergent function remains a mystery. To better understand the challenges associated with the establishment of pervasive negative supercoiling activity, we expressed the gyrase of the bacterium Thermotoga maritima in a naïve archaeon Thermococcus kodakarensis which naturally has positively supercoiled DNA. We found that the gyrase was catalytically active in T. kodakarensis leading to strong negative supercoiling of plasmid DNA which was stably maintained over at least eighty generations. An increased sensitivity of gyrase-expressing T. kodakarensis to ciprofloxacin suggested that gyrase also modulated chromosomal topology. Accordingly, global transcriptome analyses revealed large scale gene expression deregulation and identified a subset of genes responding to the negative supercoiling activity of gyrase. Surprisingly, the artificially introduced dominant negative supercoiling activity did not have a measurable effect on T. kodakarensis growth rate. Our data suggest that gyrase can become established in Thermococcales archaea without critically interfering with DNA transaction processes.

Transposon-associated TnpB is a programmable RNA-guided DNA endonuclease

Transposition has a key role in reshaping genomes of all living organisms1. Insertion sequences of IS200/IS605 and IS607 families2 are among the simplest mobile genetic elements and contain only the genes that are required for their transposition and its regulation. These elements encode tnpA transposase, which is essential for mobilization, and often carry an accessory tnpB gene, which is dispensable for transposition. Although the role of TnpA in transposon mobilization of IS200/IS605 is well documented, the function of TnpB has remained largely unknown. It had been suggested that TnpB has a role in the regulation of transposition, although no mechanism for this has been established3-5. A bioinformatic analysis indicated that TnpB might be a predecessor of the CRISPR-Cas9/Cas12 nucleases6-8. However, no biochemical activities have been ascribed to TnpB. Here we show that TnpB of Deinococcus radiodurans ISDra2 is an RNA-directed nuclease that is guided by an RNA, derived from the right-end element of a transposon, to cleave DNA next to the 5'-TTGAT transposon-associated motif. We also show that TnpB could be reprogrammed to cleave DNA target sites in human cells. Together, this study expands our understanding of transposition mechanisms by highlighting the role of TnpB in transposition, experimentally confirms that TnpB is a functional progenitor of CRISPR-Cas nucleases and establishes TnpB as a prototype of a new system for genome editing.

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.

Archaeal tyrosine recombinases

The integration of mobile genetic elements into their host chromosome influences the immediate fate of cellular organisms and gradually shapes their evolution. Site-specific recombinases catalyzing this integration have been extensively characterized both in bacteria and eukarya. More recently, a number of reports provided the in-depth characterization of archaeal tyrosine recombinases and highlighted new particular features not observed in the other two domains. In addition to being active in extreme environments, archaeal integrases catalyze reactions beyond site-specific recombination. Some of these integrases can catalyze low-sequence specificity recombination reactions with the same outcome as homologous recombination events generating deep rearrangements of their host genome. A large proportion of archaeal integrases are termed suicidal due to the presence of a specific recombination target within their own gene. The paradoxical maintenance of integrases that disrupt their gene upon integration implies novel mechanisms for their evolution. In this review, we assess the diversity of the archaeal tyrosine recombinases using a phylogenomic analysis based on an exhaustive similarity network. We outline the biochemical, ecological and evolutionary properties of these enzymes in the context of the families we identified and emphasize similarities and differences between archaeal recombinases and their bacterial and eukaryal counterparts.

Widespread divergent transcription from bacterial and archaeal promoters is a consequence of DNA-sequence symmetry

Transcription initiates at promoters, DNA regions recognized by a DNA-dependent RNA polymerase. We previously identified horizontally acquired Escherichia coli promoters from which the direction of transcription was unclear. In the present study, we show that more than half of these promoters are bidirectional and drive divergent transcription. Using genome-scale approaches, we demonstrate that 19% of all transcription start sites detected in E. coli are associated with a bidirectional promoter. Bidirectional promoters are similarly common in diverse bacteria and archaea, and have inherent symmetry: specific bases required for transcription initiation are reciprocally co-located on opposite DNA strands. Bidirectional promoters enable co-regulation of divergent genes and are enriched in both intergenic and horizontally acquired regions. Divergent transcription is conserved among bacteria, archaea and eukaryotes, but the underlying mechanisms for bidirectionality are different.

Functional Nucleic Acids Under Unusual Conditions

Functional nucleic acids (FNAs), including naturally occurring ribozymes and riboswitches as well as artificially created DNAzymes and aptamers, have been popular molecular toolboxes for diverse applications. Given the high chemical stability of nucleic acids and their ability to fold into diverse sequence-dependent structures, FNAs are suggested to be highly functional under unusual reaction conditions. This review will examine the progress of research on FNAs under conditions of low pH, high temperature, freezing conditions, and the inclusion of organic solvents and denaturants that are known to disrupt nucleic acid structures. The FNA species to be discussed include ribozymes, riboswitches, G-quadruplex-based peroxidase mimicking DNAzymes, RNA-cleaving DNAzymes, and aptamers. Research within this space has not only revealed the hidden talents of FNAs but has also laid important groundwork for pursuing these intriguing functional macromolecules for unique applications.

Extended Archaeal Histone-Based Chromatin Structure Regulates Global Gene Expression in Thermococcus kodakarensis

Histone proteins compact and organize DNA resulting in a dynamic chromatin architecture impacting DNA accessibility and ultimately gene expression. Eukaryotic chromatin landscapes are structured through histone protein variants, epigenetic marks, the activities of chromatin-remodeling complexes, and post-translational modification of histone proteins. In most Archaea, histone-based chromatin structure is dominated by the helical polymerization of histone proteins wrapping DNA into a repetitive and closely gyred configuration. The formation of the archaeal-histone chromatin-superhelix is a regulatory force of adaptive gene expression and is likely critical for regulation of gene expression in all histone-encoding Archaea. Single amino acid substitutions in archaeal histones that block formation of tightly packed chromatin structures have profound effects on cellular fitness, but the underlying gene expression changes resultant from an altered chromatin landscape have not been resolved. Using the model organism Thermococcus kodakarensis, we genetically alter the chromatin landscape and quantify the resultant changes in gene expression, including unanticipated and significant impacts on provirus transcription. Global transcriptome changes resultant from varying chromatin landscapes reveal the regulatory importance of higher-order histone-based chromatin architectures in regulating archaeal gene expression.

Histone variants in archaea and the evolution of combinatorial chromatin complexity

Nucleosomes in eukaryotes act as platforms for the dynamic integration of epigenetic information. Posttranslational modifications are reversibly added or removed and core histones exchanged for paralogous variants, in concert with changing demands on transcription and genome accessibility. Histones are also common in archaea. Their role in genome regulation, however, and the capacity of individual paralogs to assemble into histone-DNA complexes with distinct properties remain poorly understood. Here, we combine structural modeling with phylogenetic analysis to shed light on archaeal histone paralogs, their evolutionary history, and capacity to generate combinatorial chromatin states through hetero-oligomeric assembly. Focusing on the human commensal Methanosphaera stadtmanae as a model archaeal system, we show that the heteromeric complexes that can be assembled from its seven histone paralogs vary substantially in DNA binding affinity and tetramer stability. Using molecular dynamics simulations, we go on to identify unique paralogs in M. stadtmanae and Methanobrevibacter smithii that are characterized by unstable interfaces between dimers. We propose that these paralogs act as capstones that prevent stable tetramer formation and extension into longer oligomers characteristic of model archaeal histones. Importantly, we provide evidence from phylogeny and genome architecture that these capstones, as well as other paralogs in the Methanobacteriales, have been maintained for hundreds of millions of years following ancient duplication events. Taken together, our findings indicate that at least some archaeal histone paralogs have evolved to play distinct and conserved functional roles, reminiscent of eukaryotic histone variants. We conclude that combinatorially complex histone-based chromatin is not restricted to eukaryotes and likely predates their emergence.

A novel method SEProm for prokaryotic promoter prediction based on DNA structure and energetics

MOTIVATION

Despite conservation in general architecture of promoters and protein-DNA interaction interface of RNA polymerases among various prokaryotes, identification of promoter regions in the whole genome sequences remains a daunting challenge. The available tools for promoter prediction do not seem to address the problem satisfactorily, apparently because the biochemical nature of promoter signals is yet to be understood fully. Using 28 structural and 3 energetic parameters, we found that prokaryotic promoter regions have a unique structural and energy state, quite distinct from that of coding regions and the information for this signature state is in-built in their sequences. We developed a novel promoter prediction tool from these 31 parameters using various statistical techniques.

RESULTS

Here, we introduce SEProm, a novel tool that is developed by studying and utilizing the in-built structural and energy information of DNA sequences, which is applicable to all prokaryotes including archaea. Compared to five most recent, diverged and current best available tools, SEProm performs much better, predicting promoters with an 'F-value' of 82.04 and 'Precision' of 81.08. The next best 'F-value' was obtained with PromPredict (72.14) followed by BProm (68.37). On the basis of 'Precision' value, the next best 'Precision' was observed for Pepper (75.39) followed by PromPredict (72.01). SEProm maintained the lead even when comparison was done on two test organisms (not involved in training for SEProm).

AVAILABILITY AND IMPLEMENTATION

The software is freely available with easy to follow instructions (www.scfbio-iitd.res.in/software/TSS_Predict.jsp).

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Recent Advances in Archaeal Translation Initiation

Translation initiation (TI) allows accurate selection of the initiation codon on a messenger RNA (mRNA) and defines the reading frame. In all domains of life, translation initiation generally occurs within a macromolecular complex made up of the small ribosomal subunit, the mRNA, a specialized methionylated initiator tRNA, and translation initiation factors (IFs). Once the start codon is selected at the P site of the ribosome and the large subunit is associated, the IFs are released and a ribosome competent for elongation is formed. However, even if the general principles are the same in the three domains of life, the molecular mechanisms are different in bacteria, eukaryotes, and archaea and may also vary depending on the mRNA. Because TI mechanisms have evolved lately, their studies bring important information about the evolutionary relationships between extant organisms. In this context, recent structural data on ribosomal complexes and genome-wide studies are particularly valuable. This review focuses on archaeal translation initiation highlighting its relationships with either the eukaryotic or the bacterial world. Eukaryotic features of the archaeal small ribosomal subunit are presented. Ribosome evolution and TI mechanisms diversity in archaeal branches are discussed. Next, the use of leaderless mRNAs and that of leadered mRNAs having Shine-Dalgarno sequences is analyzed. Finally, the current knowledge on TI mechanisms of SD-leadered and leaderless mRNAs is detailed.

Genome-Wide Transcription Start Site Mapping and Promoter Assignments to a Sigma Factor in the Human Enteropathogen Clostridioides difficile

The emerging human enteropathogen Clostridioides difficile is the main cause of diarrhea associated with antibiotherapy. Regulatory pathways underlying the adaptive responses remain understudied and the global view of C. difficile promoter structure is still missing. In the genome of C. difficile 630, 22 genes encoding sigma factors are present suggesting a complex pattern of transcription in this bacterium. We present here the first transcriptional map of the C. difficile genome resulting from the identification of transcriptional start sites (TSS), promoter motifs and operon structures. By 5′-end RNA-seq approach, we mapped more than 1000 TSS upstream of genes. In addition to these primary TSS, this analysis revealed complex structure of transcriptional units such as alternative and internal promoters, potential RNA processing events and 5′ untranslated regions. By following an in silico iterative strategy that used as an input previously published consensus sequences and transcriptomic analysis, we identified candidate promoters upstream of most of protein-coding and non-coding RNAs genes. This strategy also led to refine consensus sequences of promoters recognized by major sigma factors of C. difficile. Detailed analysis focuses on the transcription in the pathogenicity locus and regulatory genes, as well as regulons of transition phase and sporulation sigma factors as important components of C. difficile regulatory network governing toxin gene expression and spore formation. Among the still uncharacterized regulons of the major sigma factors of C. difficile, we defined the SigL regulon by combining transcriptome and in silico analyses. We showed that the SigL regulon is largely involved in amino-acid degradation, a metabolism crucial for C. difficile gut colonization. Finally, we combined our TSS mapping, in silico identification of promoters and RNA-seq data to improve gene annotation and to suggest operon organization in C. difficile. These data will considerably improve our knowledge of global regulatory circuits controlling gene expression in C. difficile and will serve as a useful rich resource for scientific community both for the detailed analysis of specific genes and systems biology studies.

FttA is a CPSF73 homologue that terminates transcription in Archaea

Regulated gene expression is largely achieved by controlling the activities of essential, multisubunit RNA polymerase transcription elongation complexes (TECs). The extreme stability required of TECs to processively transcribe large genomic regions necessitates robust mechanisms to terminate transcription. Efficient transcription termination is particularly critical for gene-dense bacterial and archaeal genomes1-3 in which continued transcription would necessarily transcribe immediately adjacent genes and result in conflicts between the transcription and replication apparatuses4-6; the coupling of transcription and translation7,8 would permit the loading of ribosomes onto aberrant transcripts. Only select sequences or transcription termination factors can disrupt the otherwise extremely stable TEC and we demonstrate that one of the last universally conserved archaeal proteins with unknown biological function is the Factor that terminates transcription in Archaea (FttA). FttA resolves the dichotomy of a prokaryotic gene structure (operons and polarity) and eukaryotic molecular homology (general transcription apparatus) that is observed in Archaea. This missing link between prokaryotic and eukaryotic transcription regulation provides the most parsimonious link to the evolution of the processing activities involved in RNA 3'-end formation in Eukarya.

The DNA-binding protein HTa from Thermoplasma acidophilum is an archaeal histone analog

Histones are a principal constituent of chromatin in eukaryotes and fundamental to our understanding of eukaryotic gene regulation. In archaea, histones are widespread but not universal: several lineages have lost histone genes. What prompted or facilitated these losses and how archaea without histones organize their chromatin remains largely unknown. Here, we elucidate primary chromatin architecture in an archaeon without histones, Thermoplasma acidophilum, which harbors a HU family protein (HTa) that protects part of the genome from micrococcal nuclease digestion. Charting HTa-based chromatin architecture in vitro, in vivo and in an HTa-expressing E. coli strain, we present evidence that HTa is an archaeal histone analog. HTa preferentially binds to GC-rich sequences, exhibits invariant positioning throughout the growth cycle, and shows archaeal histone-like oligomerization behavior. Our results suggest that HTa, a DNA-binding protein of bacterial origin, has converged onto an architectural role filled by histones in other archaea.

Chromatinization of Escherichia coli with archaeal histones

Nucleosomes restrict DNA accessibility throughout eukaryotic genomes, with repercussions for replication, transcription, and other DNA-templated processes. How this globally restrictive organization emerged during evolution remains poorly understood. Here, to better understand the challenges associated with establishing globally restrictive chromatin, we express histones in a naive system that has not evolved to deal with nucleosomal structures: Escherichia coli. We find that histone proteins from the archaeon Methanothermus fervidus assemble on the E. coli chromosome in vivo and protect DNA from micrococcal nuclease digestion, allowing us to map binding footprints genome-wide. We show that higher nucleosome occupancy at promoters is associated with lower transcript levels, consistent with local repressive effects. Surprisingly, however, this sudden enforced chromatinization has only mild repercussions for growth unless cells experience topological stress. Our results suggest that histones can become established as ubiquitous chromatin proteins without interfering critically with key DNA-templated processes.

Identification of Dephospho-Coenzyme A (Dephospho-CoA) Kinase in Thermococcus kodakarensis and Elucidation of the Entire CoA Biosynthesis Pathway in Archaea

CoA is utilized in a wide range of metabolic pathways, and its biosynthesis is essential for all life. Pathways for CoA biosynthesis in bacteria and eukaryotes have been established. In archaea, however, the enzyme that catalyzes the final step in CoA biosynthesis, dephospho-CoA kinase (DPCK), had not been identified. In the present study, bioinformatic analyses identified a candidate for the DPCK in archaea, which was biochemically and genetically confirmed in the hyperthermophilic archaeon Thermococcus kodakarensis. Genetic analyses on genes presumed to encode bifunctional phosphopantothenoylcysteine synthetase-phosphopantothenoylcysteine decarboxylase and phosphopantetheine adenylyltransferase confirmed their involvement in CoA biosynthesis. Taken together with previous studies, the results reveal the entire pathway for CoA biosynthesis in a single archaeon and provide insight into the different mechanisms of CoA biosynthesis and their distribution in nature.

Dephospho-coenzyme A (dephospho-CoA) kinase (DPCK) catalyzes the ATP-dependent phosphorylation of dephospho-CoA, the final step in coenzyme A (CoA) biosynthesis. DPCK has been identified and characterized in bacteria and eukaryotes but not in archaea. The hyperthermophilic archaeon Thermococcus kodakarensis encodes two homologs of bacterial DPCK and the DPCK domain of eukaryotic CoA synthase, TK1334 and TK2192. We purified the recombinant TK1334 and TK2192 proteins and found that they lacked DPCK activity. Bioinformatic analyses showed that, in several archaea, the uncharacterized gene from arCOG04076 protein is fused with the gene for phosphopantetheine adenylyltransferase (PPAT), which catalyzes the reaction upstream of the DPCK reaction in CoA biosynthesis. This observation suggested that members of arCOG04076, both fused to PPAT and standalone, could be the missing archaeal DPCKs. We purified the recombinant TK1697 protein, a standalone member of arCOG04076 from T. kodakarensis, and demonstrated its GTP-dependent DPCK activity. Disruption of the TK1697 resulted in CoA auxotrophy, indicating that TK1697 encodes a DPCK that contributes to CoA biosynthesis in T. kodakarensis. TK1697 homologs are widely distributed in archaea, suggesting that the arCOG04076 protein represents a novel family of DPCK that is not homologous to bacterial and eukaryotic DPCKs but is distantly related to bacterial and eukaryotic thiamine pyrophosphokinases. We also constructed and characterized gene disruption strains of TK0517 and TK2128, homologs of bifunctional phosphopantothenoylcysteine synthetase-phosphopantothenoylcysteine decarboxylase and PPAT, respectively. Both strains displayed CoA auxotrophy, indicating their contribution to CoA biosynthesis. Taken together with previous studies, the results experimentally validate the entire CoA biosynthesis pathway in T. kodakarensis.

Next Generation DNA-Seq and Differential RNA-Seq Allow Re-annotation of the Pyrococcus furiosus DSM 3638 Genome and Provide Insights Into Archaeal Antisense Transcription

Pyrococcus furiosus DSM 3638 is a model organism for hyperthermophilic archaea with an optimal growth temperature near 100°C. The genome was sequenced about 18 years ago. However, some publications suggest that in contrast to other Pyrococcus species, the genome of P. furiosus DSM 3638 is prone to genomic rearrangements. Therefore, we re-sequenced the genome using third generation sequencing techniques. The new de novo assembled genome is 1,889,914 bp in size and exhibits high sequence identity to the published sequence. However, two major deviations were detected: (1) The genome is 18,342 bp smaller than the NCBI reference genome due to a recently described deletion. (2) The region between PF0349 and PF0388 is inverted most likely due an assembly problem for the original sequence. In addition, numerous minor variations, ranging from single nucleotide exchanges, deletions or insertions were identified. The total number of insertion sequence (IS) elements is also reduced from 30 to 24 in the new sequence. Re-sequencing of a 2-year-old “lab culture” using Nanopore sequencing confirmed the overall stability of the P. furiosus DSM 3638 genome even under normal lab conditions without taking any special care. To improve genome annotation, the updated DNA sequence was combined with an RNA sequencing approach. Here, RNAs from eight different growth conditions were pooled to increase the number of detected transcripts. Furthermore, a differential RNA-Seq approach was employed for the identification of transcription start sites (TSSs). In total, 2515 TSSs were detected and classified into 834 primary (pTSS), 797 antisense (aTSS), 739 internal and 145 secondary TSSs. Our analysis of the upstream regions revealed a well conserved archaeal promoter structure. Interrogation of the distances between pTSSs and aTSSs revealed a significant number of antisense transcripts, which are a result of bidirectional transcription from the same TATA box. This mechanism of antisense transcript production could be further confirmed by in vitro transcription experiments. We assume that bidirectional transcription gives rise to non-functional antisense RNAs and that this is a widespread phenomenon in archaea due to the architecture of the TATA element and the symmetric structure of the TATA-binding protein.

A Haloarchaeal Small Regulatory RNA (sRNA) Is Essential for Rapid Adaptation to Phosphate Starvation Conditions

The haloarchaeon Haloferax volcanii contains nearly 2800 small non-coding RNAs (sRNAs). One intergenic sRNA, sRNA132, was chosen for a detailed characterization. A deletion mutant had a growth defect and thus underscored the importance of sRNA132. A microarray analysis identified the transcript of an operon for a phosphate-specific ABC transporter as a putative target of sRNA132. Both the sRNA132 and the operon transcript accumulated under low phosphate concentrations, indicating a positive regulatory role of sRNA132. A kinetic analysis revealed that sRNA132 is essential shortly after the onset of phosphate starvation, while other regulatory processes take over after several hours. Comparison of the transcriptomes of wild-type and the sRNA132 gene deletion mutant 30 min after the onset of phosphate starvation revealed that sRNA132 controls a regulon of about 40 genes. Remarkably, the regulon included a second operon for a phosphate-specific ABC transporter, which also depended on sRNA132 for rapid induction in the absence of phosphate. Competitive growth experiments of the wild-type and ABC transporter operon deletion mutants underscored the importance of both transporters for growth at low phosphate concentrations. Northern blot analyses of four additional members of the sRNA132 regulon verified that all four transcripts depended on sRNA132 for rapid regulation after the onset of phosphate starvation. Importantly, this is the first example for the transient importance of a sRNA for any archaeal and bacterial species. In addition, this study unraveled the first sRNA regulon for haloarchaea.

Characterization of the transcriptome of Haloferax volcanii, grown under four different conditions, with mixed RNA-Seq

Haloferax volcanii is a well-established model species for haloarchaea. Small scale RNomics and bioinformatics predictions were used to identify small non-coding RNAs (sRNAs), and deletion mutants revealed that sRNAs have important regulatory functions. A recent dRNA-Seq study was used to characterize the primary transcriptome. Unexpectedly, it was revealed that, under optimal conditions, H. volcanii contains more non-coding sRNAs than protein-encoding mRNAs. However, the dRNA-Seq approach did not contain any length information. Therefore, a mixed RNA-Seq approach was used to determine transcript length and to identify additional transcripts, which are not present under optimal conditions. In total, 50 million paired end reads of 150 nt length were obtained. 1861 protein-coding RNAs (cdRNAs) were detected, which encoded 3092 proteins. This nearly doubled the coverage of cdRNAs, compared to the previous dRNA-Seq study. About 2/3 of the cdRNAs were monocistronic, and 1/3 covered more than one gene. In addition, 1635 non-coding sRNAs were identified. The highest fraction of non-coding RNAs were cis antisense RNAs (asRNAs). Analysis of the length distribution revealed that sRNAs have a median length of about 150 nt. Based on the RNA-Seq and dRNA-Seq results, genes were chosen to exemplify characteristics of the H. volcanii transcriptome by Northern blot analyses, e.g. 1) the transcript patterns of gene clusters can be straightforward, but also very complex, 2) many transcripts differ in expression level under the four analyzed conditions, 3) some genes are transcribed into RNA isoforms of different length, which can be differentially regulated, 4) transcripts with very long 5'-UTRs and with very long 3'-UTRs exist, and 5) about 30% of all cdRNAs have overlapping 3'-ends, which indicates, together with the asRNAs, that H. volcanii makes ample use of sense-antisense interactions. Taken together, this RNA-Seq study, together with a previous dRNA-Seq study, enabled an unprecedented view on the H. volcanii transcriptome.

The Primary Antisense Transcriptome of Halobacterium salinarum NRC-1

Antisense RNAs (asRNAs) are present in diverse organisms and play important roles in gene regulation. In this work, we mapped the primary antisense transcriptome in the halophilic archaeon Halobacterium salinarum NRC-1. By reanalyzing publicly available data, we mapped antisense transcription start sites (aTSSs) and inferred the probable 3′ ends of these transcripts. We analyzed the resulting asRNAs according to the size, location, function of genes on the opposite strand, expression levels and conservation. We show that at least 21% of the genes contain asRNAs in H. salinarum. Most of these asRNAs are expressed at low levels. They are located antisense to genes related to distinctive characteristics of H. salinarum, such as bacteriorhodopsin, gas vesicles, transposases and other important biological processes such as translation. We provide evidence to support asRNAs in type II toxin–antitoxin systems in archaea. We also analyzed public Ribosome profiling (Ribo-seq) data and found that ~10% of the asRNAs are ribosome-associated non-coding RNAs (rancRNAs), with asRNAs from transposases overrepresented. Using a comparative transcriptomics approach, we found that ~19% of the asRNAs annotated in H. salinarum belong to genes with an ortholog in Haloferax volcanii, in which an aTSS could be identified with positional equivalence. This shows that most asRNAs are not conserved between these halophilic archaea.

Genome-wide screening of potential RNase Y-processed mRNAs in the M49 serotype Streptococcus pyogenes NZ131

RNase Y is a major endoribonuclease in Group A streptococcus (GAS) and other Gram-positive bacteria. Our previous study showed that RNase Y was involved in mRNA degradation and processing in GAS. We hypothesized that mRNA processing regulated the expression of important GAS virulence factors via altering their mRNA stabilities and that RNase Y mediated at least some of the mRNA-processing events. The aims of this study were to (1) identify mRNAs that were processed by RNase Y and (2) confirm the mRNA-processing events. The transcriptomes of Streptococcus pyogenes NZ131 wild type and its RNase Y mutant (Δrny) were examined with RNA-seq. The data were further analyzed to define GAS operons. The mRNA stabilities of the wild type and Δrny at subgene level were determined with tiling array analysis. Operons displaying segmental stability in the wild type but not in the Δrny were predicted to be RNase Y processed. Overall 865 operons were defined and their boundaries predicted. Further analysis narrowed down 15 mRNAs potentially processed by RNase Y. A selection of four candidates including folC1 (folylpolyglutamate synthetase), prtF (fibronectin-binding protein), speG (streptococcal exotoxin G), ropB (transcriptional regulator of speB), and ypaA (riboflavin transporter) mRNAs was examined with Northern blot analysis. However, only folC1 was confirmed to be processed, but it is unlikely that RNase Y is responsible. We conclude that GAS use RNase Y to selectively process mRNA, but the overall impact is confined to selected virulence factors.

The yeast C/D box snoRNA U14 adopts a "weak" K-turn like conformation recognized by the Snu13 core protein in solution

Non-coding RNAs associate with proteins to form ribonucleoproteins (RNPs), such as ribosome, box C/D snoRNPs, H/ACA snoRNPs, ribonuclease P, telomerase and spliceosome to ensure cell viability. The assembly of these RNA-protein complexes relies on the ability of the RNA to adopt the correct bound conformation. K-turn motifs represent ubiquitous binding platform for proteins found in several cellular environment. This structural motif has an internal three-nucleotide bulge flanked on its 3' side by a G•A/A•G tandem pairs followed by one or two non-Watson-Crick pairs, and on its 5' side by a classical RNA helix. This peculiar arrangement induces a strong curvature of the phosphodiester backbone, which makes it conducive to multiple tertiary interactions. SNU13/Snu13p (Human/Yeast) binds specifically the U14 C/D box snoRNA K-turn sequence motif. This event is the prerequisite to promote the assembly of the RNP, which contains NOP58/Nop58 and NOP56/Nop56 core proteins and the 2'-O-methyl-transferase, Fibrillarin/Nop1p. The U14 small nucleolar RNA is a conserved non-coding RNA found in yeast and vertebrates required for the pre-rRNA maturation and ribose methylation. Here, we report the solution structure of the free U14 snoRNA K-turn motif (kt-U14) as determined by Nuclear Magnetic Resonance. We demonstrate that a major fraction of free kt-U14 adopts a pre-folded conformation similar to protein bound K-turn, even in the absence of divalent ions. In contrast to the kt-U4 or tyrS RNA, kt-U14 displays a sharp bent in the phosphodiester backbone. The U•U and G•A tandem base pairs are formed through weak hydrogen bonds. Finally, we show that the structure of kt-U14 is stabilized upon Snu13p binding. The structure of the free U14 RNA is the first reference example for the canonical motifs of the C/D box snoRNA family.

Distinct Physiological Roles of the Three Ferredoxins Encoded in the Hyperthermophilic Archaeon Thermococcus kodakarensis

High-energy electrons liberated during catabolic processes can be exploited for energy-conserving mechanisms. Maximal energy gains demand these valuable electrons be accurately shuttled from electron donor to appropriate electron acceptor. Proteinaceous electron carriers such as ferredoxins offer opportunities to exploit specific ferredoxin partnerships to ensure that electron flux to critical physiological pathways is aligned with maximal energy gains. Most species encode many ferredoxin isoforms, but very little is known about the role of individual ferredoxins in most systems. Our results detail that ferredoxin isoforms make largely unique and distinct protein interactions in vivo and that flux through one ferredoxin often cannot be recovered by flux through a different ferredoxin isoform. The results obtained more broadly suggest that ferredoxin isoforms throughout biological life have evolved not as generic electron shuttles, but rather serve as selective couriers of valuable low-potential electrons from select electron donors to desirable electron acceptors.

Control of electron flux is critical in both natural and bioengineered systems to maximize energy gains. Both small molecules and proteins shuttle high-energy, low-potential electrons liberated during catabolism through diverse metabolic landscapes. Ferredoxin (Fd) proteins—an abundant class of Fe-S-containing small proteins—are essential in many species for energy conservation and ATP production strategies. It remains difficult to model electron flow through complicated metabolisms and in systems in which multiple Fd proteins are present. The overlap of activity and/or limitations of electron flux through each Fd can limit physiology and metabolic engineering strategies. Here we establish the interplay, reactivity, and physiological role(s) of the three ferredoxin proteins in the model hyperthermophile Thermococcus kodakarensis. We demonstrate that the three loci encoding known Fds are subject to distinct regulatory mechanisms and that specific Fds are utilized to shuttle electrons to separate respiratory and energy production complexes during different physiological states. The results obtained argue that unique physiological roles have been established for each Fd and that continued use of T. kodakarensis and related hydrogen-evolving species as bioengineering platforms must account for the distinct Fd partnerships that limit flux to desired electron acceptors. Extrapolating our results more broadly, the retention of multiple Fd isoforms in most species argues that specialized Fd partnerships are likely to influence electron flux throughout biology.

TFS and Spt4/5 accelerate transcription through archaeal histone-based chromatin

RNA polymerase must surmount translocation barriers for continued transcription. In Eukarya and most Archaea, DNA-bound histone proteins represent the most common and troublesome barrier to transcription elongation. Eukaryotes encode a plethora of chromatin-remodeling complexes, histone-modification enzymes and transcription elongation factors to aid transcription through nucleosomes, while archaea seemingly lack machinery to remodel/modify histone-based chromatin and thus must rely on elongation factors to accelerate transcription through chromatin-barriers. TFS (TFIIS in Eukarya) and the Spt4-Spt5 complex are universally encoded in archaeal genomes, and here we demonstrate that both elongation factors, via different mechanisms, can accelerate transcription through archaeal histone-based chromatin. Histone proteins in Thermococcus kodakarensis are sufficiently abundant to completely wrap all genomic DNA, resulting in a consistent protein barrier to transcription elongation. TFS-enhanced cleavage of RNAs in backtracked transcription complexes reactivates stalled RNAPs and dramatically accelerates transcription through histone-barriers, while Spt4-Spt5 changes to clamp-domain dynamics play a lesser-role in stabilizing transcription. Repeated attempts to delete TFS, Spt4 and Spt5 from the T. kodakarensis genome were not successful, and the essentiality of both conserved transcription elongation factors suggests that both conserved elongation factors play important roles in transcription regulation in vivo, including mechanisms to accelerate transcription through downstream protein barriers.

A linear pathway for mevalonate production supports growth of Thermococcus kodakarensis

The sole unifying feature of Archaea is the use of isoprenoid-based glycerol lipid ethers to compose cellular membranes. The branched hydrocarbon tails of archaeal lipids are synthesized via the polymerization of isopentenyl diphosphate (IPP) and dimethylallyl diphosphate (DMAPP), but many questions still surround the pathway(s) that result in production of IPP and DMAPP in archaeal species. Isotopic-labeling strategies argue for multiple biological routes for production of mevalonate, but biochemical and bioinformatic studies support only a linear pathway for mevalonate production. Here, we use a combination of genetic and biochemical assays to detail the production of mevalonate in the model archaeon Thermococcus kodakarensis. We demonstrate that a single, linear pathway to mevalonate biosynthesis is essential and that alternative routes of mevalonate production, if present, are not biologically sufficient to support growth in the absence of the classical mevalonate pathway resulting in IPP production from acetyl-CoA. Archaeal species provide an ideal platform for production of high-value isoprenoids in large quantities, and the results obtained provide avenues to further increase the production of mevalonate to drive isoprenoid production in archaeal hosts.

The role of RNA structure in translational regulation by L7Ae protein in archaea

A recent study has shown that archaeal L7Ae binds to a putative k-turn structure in the 5'-leader of the mRNA of its structural gene to regulate translation. To function as a regulator, the RNA should be unstructured in the absence of protein, but it should adopt a k-turn-containing stem-loop on binding L7Ae. Sequence analysis of UTR sequences indicates that their k-turn elements will be unable to fold in the absence of L7Ae, and we have demonstrated this experimentally in solution using FRET for the Archaeoglobus fulgidus sequence. We have solved the X-ray crystal structure of the complex of the A. fulgidus RNA bound to its cognate L7Ae protein. The RNA adopts a standard k-turn conformation that is specifically recognized by the L7Ae protein, so stabilizing the stem-loop. In-line probing of the natural-sequence UTR shows that the RNA is unstructured in the absence of L7Ae binding, but folds on binding the protein such that the ribosome binding site is occluded. Thus, L7Ae regulates its own translation by switching the conformation of the RNA to alter accessibility.

Insights into RNA-processing pathways and associated RNA-degrading enzymes in Archaea

RNA-processing pathways are at the centre of regulation of gene expression. All RNA transcripts undergo multiple maturation steps in addition to covalent chemical modifications to become functional in the cell. This includes destroying unnecessary or defective cellular RNAs. In Archaea, information on mechanisms by which RNA species reach their mature forms and associated RNA-modifying enzymes are still fragmentary. To date, most archaeal actors and pathways have been proposed in light of information gathered from Bacteria and Eukarya. In this context, this review provides a state of the art overview of archaeal endoribonucleases and exoribonucleases that cleave and trim RNA species and also of the key small archaeal proteins that bind RNAs. Furthermore, synthetic up-to-date views of processing and biogenesis pathways of archaeal transfer and ribosomal RNAs as well as of maturation of stable small non-coding RNAs such as CRISPR RNAs, small C/D and H/ACA box guide RNAs, and other emerging classes of small RNAs are described. Finally, prospective post-transcriptional mechanisms to control archaeal messenger RNA quality and quantity are discussed.

Archaeal DNA replication and repair: new genetic, biophysical and molecular tools for discovering and characterizing enzymes, pathways and mechanisms

DNA replication and repair are essential biological processes needed for the survival of all organisms. Although these processes are fundamentally conserved in the three domains, archaea, bacteria and eukarya, the proteins and complexes involved differ. The genetic and biophysical tools developed for archaea in the last several years have accelerated the study of DNA replication and repair in this domain. In this review, the current knowledge of DNA replication and repair processes in archaea will be summarized, with emphasis on the contribution of genetics and other recently developed biophysical and molecular tools, including capillary gel electrophoresis, next-generation sequencing and single-molecule approaches. How these new tools will continue to drive archaeal DNA replication and repair research will also be discussed.

Unique Archaeal Small RNAs

Advances in genome-wide sequence technologies allow for detailed insights into the complexity of RNA landscapes of organisms from all three domains of life. Recent analyses of archaeal transcriptomes identified interaction and regulation networks of noncoding RNAs in this understudied domain. Here, we review current knowledge of small, noncoding RNAs with important functions for the archaeal lifestyle, which often requires adaptation to extreme environments. One focus is RNA metabolism at elevated temperatures in hyperthermophilic archaea, which reveals elevated amounts of RNA-guided RNA modification and virus defense strategies. Genome rearrangement events result in unique fragmentation patterns of noncoding RNA genes that require elaborate maturation pathways to yield functional transcripts. RNA-binding proteins, e.g., L7Ae and LSm, are important for many posttranscriptional control functions of RNA molecules in archaeal cells. We also discuss recent insights into the regulatory potential of their noncoding RNA partners.

Internal RNAs overlapping coding sequences can drive the production of alternative proteins in archaea

Prokaryotic genomes show a high level of information compaction often with different molecules transcribed from the same locus. Although antisense RNAs have been relatively well studied, RNAs in the same strand, internal RNAs (intraRNAs), are still poorly understood. The question of how common is the translation of overlapping reading frames remains open. We address this question in the model archaeon Halobacterium salinarum. In the present work we used differential RNA-seq (dRNA-seq) in H. salinarum NRC-1 to locate intraRNA signals in subsets of internal transcription start sites (iTSS) and establish the open reading frames associated to them (intraORFs). Using C-terminally flagged proteins, we experimentally observed isoforms accurately predicted by intraRNA translation for kef1, acs3 and orc4 genes. We also recovered from the literature and mass spectrometry databases several instances of protein isoforms consistent with intraRNA translation such as the gas vesicle protein gene gvpC1. We found evidence for intraRNAs in horizontally transferred genes such as the chaperone dnaK and the aerobic respiration related cydA in both H. salinarum and Escherichia coli. Also, intraRNA translation evidence in H. salinarum, E. coli and yeast of a universal elongation factor (aEF-2, fusA and eEF-2) suggests that this is an ancient phenomenon present in all domains of life.

Fine-Tuned Protein Production in Methanosarcina acetivorans C2A

Methanogenic archaea can be integrated into a sustainable, carbon-neutral cycle for producing organic chemicals from C₁ compounds if the rate, yield, and titer of product synthesis can be improved using metabolic engineering. However, metabolic engineering techniques are limited in methanogens by insufficient methods for controlling cellular protein levels. We conducted a systematic approach to tune protein levels in Methanosarcina acetivorans C2A, a model methanogen, by regulating transcription and translation initiation. Rationally designed core promoter and ribosome binding site mutations in M. acetivorans C2A resulted in a predicable change in protein levels over a 60 fold range. The overall range of protein levels was increased an additional 3 fold by introducing the 5' untranslated region of the mcrB transcript. This work demonstrates a wide range of precisely controlled protein levels in M. acetivorans C2A, which will help facilitate systematic metabolic engineering efforts in methanogens.

Toward a Universal Structural and Energetic Model for Prokaryotic Promoters

With almost no consensus promoter sequence in prokaryotes, recruitment of RNA polymerase (RNAP) to precise transcriptional start sites (TSSs) has remained an unsolved puzzle. Uncovering the underlying mechanism is critical for understanding the principle of gene regulation. We attempted to search the hidden code in ∼16,500 promoters of 12 prokaryotes representing two kingdoms in their structure and energetics. Twenty-eight fundamental parameters of DNA structure including backbone angles, basepair axis, and interbasepair and intrabasepair parameters were used, and information was extracted from x-ray crystallography data. Three parameters (solvation energy, hydrogen-bond energy, and stacking energy) were selected for creating energetics profiles using in-house programs. DNA of promoter regions was found to be inherently designed to undergo a change in every parameter undertaken for the study, in all prokaryotes. The change starts from some distance upstream of TSSs and continues past some distance from TSS, hence giving a signature state to promoter regions. These signature states might be the universal hidden codes recognized by RNAP. This observation was reiterated when randomly selected promoter sequences (with little sequence conservation) were subjected to structure generation; all developed into very similar three-dimensional structures quite distinct from those of conventional B-DNA and coding sequences. Fine structural details at important motifs (viz. -11, -35, and -75 positions relative to TSS) of promoters reveal novel to our knowledge and pointed insights for RNAP interaction at these locations; it could be correlated with how some particular structural changes at the -11 region may allow insertion of RNAP amino acids in interbasepair space as well as facilitate the flipping out of bases from the DNA duplex.

Modeling leaderless transcription and atypical genes results in more accurate gene prediction in prokaryotes

In a conventional view of the prokaryotic genome organization, promoters precede operons and ribosome binding sites (RBSs) with Shine-Dalgarno consensus precede genes. However, recent experimental research suggesting a more diverse view motivated us to develop an algorithm with improved gene-finding accuracy. We describe GeneMarkS-2, an ab initio algorithm that uses a model derived by self-training for finding species-specific (native) genes, along with an array of precomputed “heuristic” models designed to identify harder-to-detect genes (likely horizontally transferred). Importantly, we designed GeneMarkS-2 to identify several types of distinct sequence patterns (signals) involved in gene expression control, among them the patterns characteristic for leaderless transcription as well as noncanonical RBS patterns. To assess the accuracy of GeneMarkS-2, we used genes validated by COG (Clusters of Orthologous Groups) annotation, proteomics experiments, and N-terminal protein sequencing. We observed that GeneMarkS-2 performed better on average in all accuracy measures when compared with the current state-of-the-art gene prediction tools. Furthermore, the screening of ∼5000 representative prokaryotic genomes made by GeneMarkS-2 predicted frequent leaderless transcription in both archaea and bacteria. We also observed that the RBS sites in some species with leadered transcription did not necessarily exhibit the Shine-Dalgarno consensus. The modeling of different types of sequence motifs regulating gene expression prompted a division of prokaryotic genomes into five categories with distinct sequence patterns around the gene starts.

The Transcriptional Regulator TFB-RF1 Activates Transcription of a Putative ABC Transporter in Pyrococcus furiosus

Transcription factor B recruiting factor 1 (TFB-RF1; PF1088) is a transcription regulator which activates transcription on archaeal promoters containing weak TFB recognition elements (BRE) by recruiting TFB to the promoter. The mechanism of activation is described in detail, but nothing is known about the biological function of this protein in Pyrococcus furiosus. The protein is located in an operon structure together with the hypothetical gene pf1089 and western blot as well as end-point RT-PCR experiments revealed an extremely low expression rate of both proteins. Furthermore, conditions to induce the expression of the operon are not known. By introducing an additional copy of tfb-RF1 using a Pyrococcus shuttle vector we could circumvent the lacking expression of both proteins under standard growth conditions as indicated by western blot as well as end-point RT-PCR experiments. A ChIP-seq experiment revealed an additional binding site of TFB-RF1 in the upstream region of the pf1011/1012 operon, beside the expected target of the pf1089/tfb-RF1 region. This operon codes for a putative ABC transporter which is most-related to a multidrug export system and in vitro analysis using gel shift assays, DNase I footprinting and in vitro transcription confirmed the activator function of TFB-RF1 on the corresponding promoter. These findings are also in agreement with in vivo data, as RT-qPCR experiments also indicate transcriptional activation of both operons. Taken together, the overexpression strategy of tfb-RF1 enabled the identification of an additional operon of the TFB-RF1 regulon which indicates a transport-related function and provides a promising starting position to decipher the physiological function of the TFB-RF1 gene regulatory network in P. furiosus.

An Archaeal Fluoride-Responsive Riboswitch Provides an Inducible Expression System for Hyperthermophiles

Robust genetic systems for the hyperthermophilic Thermococcales have facilitated the overexpression of native genes, enabled the addition of sequences encoding secretion signals, epitope, and affinity tags to coding regions, and aided the introduction of sequences encoding new proteins in these fast-growing fermentative heterotrophs. However, tightly controlled and easily manipulated systems facilitating regulated gene expression are limited for these hosts. Here, we describe an alternative method for regulatory control reliant on a cis-encoded functional riboswitch in the model archaeon Thermococcus kodakarensis Despite the hyperthermophilic growth temperatures, the proposed structure of the riboswitch conforms to a fluoride-responsive riboswitch encoded in many bacteria and similarly functions to regulate a component-conserved fluoride export pathway. Deleting components of the fluoride export pathway generates T. kodakarensis strains with increased fluoride sensitivity. The mechanism underlying regulated expression suggested that the riboswitch-encoding sequences could be utilized as a tunable expression cassette. When appended to a reporter gene, the riboswitch-mediated control system provides fluoride-dependent tunable regulatory potential, offering an alternative system for regulating gene expression. Riboswitch-regulated expression is thus ubiquitous in extant life and can be exploited to generate regulated expression systems for hyperthermophiles.IMPORTANCE Gene expression is controlled by a myriad of interconnected mechanisms that interpret metabolic states and environmental cues to balance cell physiology. Transcription regulation in Archaea is known to employ both typical repressors-operators and transcription activators to regulate transcription initiation in addition to the regulation afforded by chromatin structure. It was perhaps surprising that the presumed ancient mechanism of riboswitch-mediated regulation is found in Bacteria and Eukarya, but seemingly absent in Archaea We demonstrate here that a fluoride-responsive riboswitch functions to regulate a detoxification pathway in the hyperthermophilic archaeon Thermococcus kodakarensis The results obtained define a universal role for riboswitch-mediated regulation, adumbrate the presence of several riboswitch-regulated genes in Thermococcus kodakarensis, demonstrate the utility of RNA-based regulation at high temperatures, and provide a novel riboswitch-regulated expression system to employ in hyperthermophiles.

The Non-Coding Regulatory RNA Revolution in Archaea

Small non-coding RNAs (sRNAs) are ubiquitously found in the three domains of life playing large-scale roles in gene regulation, transposable element silencing and defense against foreign elements. While a substantial body of experimental work has been done to uncover function of sRNAs in Bacteria and Eukarya, the functional roles of sRNAs in Archaea are still poorly understood. Recently, high throughput studies using RNA-sequencing revealed that sRNAs are broadly expressed in the Archaea, comprising thousands of transcripts within the transcriptome during non-challenged and stressed conditions. Antisense sRNAs, which overlap a portion of a gene on the opposite strand (cis-acting), are the most abundantly expressed non-coding RNAs and they can be classified based on their binding patterns to mRNAs (3' untranslated region (UTR), 5' UTR, CDS-binding). These antisense sRNAs target many genes and pathways, suggesting extensive roles in gene regulation. Intergenic sRNAs are less abundantly expressed and their targets are difficult to find because of a lack of complete overlap between sRNAs and target mRNAs (trans-acting). While many sRNAs have been validated experimentally, a regulatory role has only been reported for very few of them. Further work is needed to elucidate sRNA-RNA binding mechanisms, the molecular determinants of sRNA-mediated regulation, whether protein components are involved and how sRNAs integrate with complex regulatory networks.

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

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.