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

Computational analysis of genes with lethal knockout phenotype and prediction of essential genes in archaea

The identification of microbial genes essential for survival as those with lethal knockout phenotype (LKP) is a common strategy for functional interrogation of genomes. However, interpretation of the LKP is complicated because a substantial fraction of the genes with this phenotype remains poorly functionally characterized. Furthermore, many genes can exhibit LKP not because their products perform essential cellular functions but because their knockout activates the toxicity of other genes (conditionally essential genes). We analyzed the sets of LKP genes for two archaea, Methanococcus maripaludis and Sulfolobus islandicus, using a variety of computational approaches aiming to differentiate between essential and conditionally essential genes and to predict at least a general function for as many of the proteins encoded by these genes as possible. This analysis allowed us to predict the functions of several LKP genes including previously uncharacterized subunit of the GINS protein complex with an essential function in genome replication and of the KEOPS complex that is responsible for an essential tRNA modification as well as GRP protease implicated in protein quality control. Additionally, several novel antitoxins (conditionally essential genes) were predicted, and this prediction was experimentally validated by showing that the deletion of these genes together with the adjacent genes apparently encoding the cognate toxins caused no growth defect. We applied principal component analysis based on sequence and comparative genomic features showing that this approach can separate essential genes from conditionally essential ones and used it to predict essential genes in other archaeal genomes.IMPORTANCEOnly a relatively small fraction of the genes in any bacterium or archaeon is essential for survival as demonstrated by the lethal effect of their disruption. The identification of essential genes and their functions is crucial for understanding fundamental cell biology. However, many of the genes with a lethal knockout phenotype remain poorly functionally characterized, and furthermore, many genes can exhibit this phenotype not because their products perform essential cellular functions but because their knockout activates the toxicity of other genes. We applied state-of-the-art computational methods to predict the functions of a number of uncharacterized genes with the lethal knockout phenotype in two archaeal species and developed a computational approach to predict genes involved in essential functions. These findings advance the current understanding of key functionalities of archaeal cells.

The transcriptional regulator EarA and intergenic terminator sequences modulate archaellation in Pyrococcus furiosus

The regulation of archaellation, the formation of archaeal-specific cell appendages called archaella, is crucial for the motility, adhesion, and survival of archaeal organisms. Although the heavily archaellated and highly motile Pyrococcus furiosus is a key model organism for understanding the production and function of archaella in Euryarchaea, the transcriptional regulation of archaellum assembly is so far unknown. Here we show that the transcription factor EarA is the master regulator of the archaellum (arl) operon transcription, which is further modulated by intergenic transcription termination signals. EarA deletion or overexpression strains demonstrate that EarA is essential for archaellation in P. furiosus and governs the degree of archaellation. Providing a single-molecule update on the transcriptional landscape of the arl operon in P. furiosus, we identify sequence motifs for EarA binding upstream of the arl operon and intergenic terminator sequences as critical elements for fine-tuning the expression of the multicistronic arl cluster. Furthermore, transcriptome re-analysis across different Thermococcales species demonstrated a heterogeneous production of major archaellins, suggesting a more diverse composition of archaella than previously recognized. Overall, our study provides novel insights into the transcriptional regulation of archaellation and highlights the essential role of EarA in Pyrococcus furiosus. These findings advance our understanding of the mechanisms governing archaellation and have implications for the functional diversity of archaella.

The archaeal Lsm protein from Pyrococcus furiosus binds co-transcriptionally to poly(U)-rich target RNAs

Posttranscriptional processes in Bacteria include the association of small regulatory RNAs (sRNA) with a target mRNA. The sRNA/mRNA annealing process is often mediated by an RNA chaperone called Hfq. The functional role of bacterial and eukaryotic Lsm proteins is partially understood, whereas knowledge about archaeal Lsm proteins is scarce. Here, we used the genetically tractable archaeal hyperthermophile Pyrococcus furiosus to identify the protein interaction partners of the archaeal Sm-like proteins (PfuSmAP1) using mass spectrometry and performed a transcriptome-wide binding site analysis of PfuSmAP1. Most of the protein interaction partners we found are part of the RNA homoeostasis network in Archaea including ribosomal proteins, the exosome, RNA-modifying enzymes, but also RNA polymerase subunits, and transcription factors. We show that PfuSmAP1 preferentially binds messenger RNAs and antisense RNAs recognizing a gapped poly(U) sequence with high affinity. Furthermore, we found that SmAP1 co-transcriptionally associates with target RNAs. Our study reveals that in contrast to bacterial Hfq, PfuSmAP1 does not affect the transcriptional activity or the pausing behaviour of archaeal RNA polymerases. We propose that PfuSmAP1 recruits antisense RNAs to target mRNAs and thereby executes its putative regulatory function on the posttranscriptional level.

Uncovering the temporal dynamics and regulatory networks of thermal stress response in a hyperthermophile using transcriptomics and proteomics

Facing rapid fluctuations in their natural environment, extremophiles, like the hyperthermophilic archaeon Pyrococcus furiosus, exhibit remarkable adaptability to extreme conditions. However, our understanding of their dynamic cellular responses remains limited. This study integrates RNA-sequencing and mass spectrometry data, thereby elucidating transcriptomic and proteomic responses to heat and cold shock stress in P. furiosus. Our results reveal rapid and dynamic changes in gene and protein expression following these stress responses. Heat shock triggers extensive transcriptome reprogramming, orchestrated by the transcriptional regulator Phr, targeting a broader gene repertoire than previously demonstrated. For heat shock signature genes, RNA levels swiftly return to baseline upon recovery, while protein levels remain persistently upregulated, reflecting a rapid but sustained response. Intriguingly, cold shock at 4°C elicits distinct short- and long-term responses at both RNA and protein levels. Cluster analysis identified gene sets with either congruent or contrasting trends in RNA and protein changes, representing well-separated arCOG groups tailored to their individual cellular responses. Particularly, upregulation of ribosomal proteins and significant enrichment of 5'-leadered sequences in cold-shock responsive genes suggest that translation regulation is important during cold shock adaption. Further investigating transcriptomic features, we reveal that thermal stress genes are equipped with basal sequence elements, such as strong promoter and poly(U)-terminators, facilitating a regulated response of the respective transcription units. Our study provides a comprehensive overview of the cellular response to temperature stress, advancing our understanding of stress response mechanisms in hyperthermophilic archaea and providing valuable insights into the molecular adaptations that facilitate life in extreme environments.IMPORTANCEExtreme environments provide unique challenges for life, and the study of extremophiles can shed light on the mechanisms of adaptation to such conditions. Pyrococcus furiosus, a hyperthermophilic archaeon, is a model organism for studying thermal stress response mechanisms. In this study, we used an integrated analysis of RNA-sequencing and mass spectrometry data to investigate the transcriptomic and proteomic responses of P. furiosus to heat and cold shock stress and recovery. Our results reveal the rapid and dynamic changes in gene and protein expression patterns associated with these stress responses, as well as the coordinated regulation of different gene sets in response to different stressors. These findings provide valuable insights into the molecular adaptations that facilitate life in extreme environments and advance our understanding of stress response mechanisms in hyperthermophilic archaea.

Current and emerging tools of computational biology to improve the detoxification of mycotoxins

Biological organisms carry a rich potential for removing toxins from our environment, but identifying suitable candidates and improving them remain challenging. We explore the use of computational tools to discover strains and enzymes that detoxify harmful compounds. In particular, we will focus on mycotoxins-fungi-produced toxins that contaminate food and feed-and biological enzymes that are capable of rendering them less harmful. We discuss the use of established and novel computational tools to complement existing empirical data in three directions: discovering the prospect of detoxification among underexplored organisms, finding important cellular processes that contribute to detoxification, and improving the performance of detoxifying enzymes. We hope to create a synergistic conversation between researchers in computational biology and those in the bioremediation field. We showcase open bioremediation questions where computational researchers can contribute and highlight relevant existing and emerging computational tools that could benefit bioremediation researchers.

Chromosome Genome Sequencing and Comparative Transcriptome-Based Analyses of Kloeckera apiculata 34-9 Unveil the Potential Biocontrol Mechanisms Against Citrus Green Mold

Biological control is an environmentally friendly, safe, and replaceable strategy for disease management. Genome sequences of a certain biocontrol agent could lay a solid foundation for the research of molecular biology, and the more refined the reference genome, the more information it provides. In the present study, a higher resolution genome of Kloeckera apiculata 34-9 was assembled using high-throughput chromosome conformation capture (Hi-C) technology. A total of 8.07 M sequences of K. apiculata 34-9 genome was anchored onto 7 pesudochromosomes, which accounting for about 99.51% of the whole assembled sequences, and 4,014 protein-coding genes were annotated. Meanwhile, the detailed gene expression changes of K. apiculata 34-9 were obtained under low temperature and co-incubation with Penicillium digitatum treatments, respectively. Totally 254 differentially expressed genes (DEGs) were detected with low temperature treatment, of which 184 and 70 genes were upregulated and downregulated, respectively. Some candidate genes were significantly enriched in ribosome biosynthesis in eukaryotes and ABC transporters. The expression of gene Kap003732 and Kap001595 remained upregulated and downregulated through the entire time-points, respectively, indicating that they might be core genes for positive and negative response to low temperature stress. When co-incubation with P. digitatum, a total of 2,364 DEGs were found, and there were 1,247 upregulated and 1,117 downregulated genes, respectively. Biosynthesis of lysine and arginine, and phenylalanine metabolism were the highest enrichment of the cluster and KEGG analyses of the co-DEGs, the results showed that they might be involved in the positive regulation of K. apiculata 34-9 response to P. digitatum. The completeness of K. apiculata 34-9 genome and the transcriptome data presented here are essential for providing a high-quality genomic resource and it might serve as valuable molecular properties for further studies on yeast genome, expression pattern of biocontrol system, and postharvest citrus storage and preservation.

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.

CopR, a Global Regulator of Transcription to Maintain Copper Homeostasis in Pyrococcus furiosus

Although copper is in many cases an essential micronutrient for cellular life, higher concentrations are toxic. Therefore, all living cells have developed strategies to maintain copper homeostasis. In this manuscript, we have analyzed the transcriptome-wide response of Pyrococcus furiosus to increased copper concentrations and described the essential role of the putative copper-sensing metalloregulator CopR in the detoxification process. To this end, we employed biochemical and biophysical methods to characterize the role of CopR. Additionally, a copR knockout strain revealed an amplified sensitivity in comparison to the parental strain towards increased copper levels, which designates an essential role of CopR for copper homeostasis. To learn more about the CopR-regulated gene network, we performed differential gene expression and ChIP-seq analysis under normal and 20 μM copper-shock conditions. By integrating the transcriptome and genome-wide binding data, we found that CopR binds to the upstream regions of many copper-induced genes. Negative-stain transmission electron microscopy and 2D class averaging revealed an octameric assembly formed from a tetramer of dimers for CopR, similar to published crystal structures from the Lrp family. In conclusion, we propose a model for CopR-regulated transcription and highlight the regulatory network that enables Pyrococcus to respond to increased copper concentrations.

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.