Differential RNA sequencing (dRNA-seq) has been proven to be a valuable tool for studying bacterial and archaeal genomes by identifying transcription start sites (TSS). With the help of statistical analysis, researchers can analyze dRNA-seq reads from treated and untreated libraries based on the activity of the 5-monophosphate dependent terminator RNA exonuclease (TEX). Our study focuses on Halobacterium salinarum NRC-1, a type of extremophilic archaeon that is commonly found in highly saline environments and is a model organism for studying molecular biology and genetics in extreme environments. The objective of our study is to identify and map Alternative Transcription Start Sites (alTSS), Transcript Process Sites (TPS), and Transcription Termination Sites (TTS) using the dRNA-seq data in H. salinarum NRC-1. We modified the TSSAR tool to detect TSSs from dRNA-seq data, assuming that sequencing reads start at the exact position during transcription and follow a Gamma-Poisson distribution. We annotated alTSSs into four types based on the number of dRNA-seq library reads and differences between two main TSS locations. Alternative TSSs have lower RNA reads than primary ones and can have upstream open reading frames, leading to changes in gene regulation output, 5UTR isoform, and gene transcription pausing. Mapping alTSSs allowed us to explain changes in cell response to growth conditions and gene expression across different growth stages. Our findings revealed a significant number of falsely annotated internal transcription start sites (iTSSs) that were redefined as alternative TSSs in previous genome annotations. These alternative TSSs produced different protein isoforms depending on the length of the amino acid chain and the open reading frame. Additionally, alternative TSSs were identified not only in H. salinarum but also in other organisms, suggesting a crucial role in regulating gene expression across various species. Furthermore, we conducted a re-analysis of dRNA-seq data, focusing on non-primary transcripts (monophosphorylated RNAs) instead of the traditional method of enriching for primary transcripts (triphosphorylated RNAs) to identify genome-wide transcript processing sites (TPS) in H. salinarum NRC-1. We also applied this approach to Haloferax volcanii for comparative analysis. Lastly, we used dRNA-seq data to identify hairpin structures and mapped them onto the genome, providing insights into the potential role of transcription termination sites (TTS).

O sequenciamento diferencial de RNA (dRNA-seq) tem se mostrado uma ferramenta valiosa para o estudo de genomas bacterianos e arqueanos por meio da identificação de sítios de início de transcrição (TSS, do inglês Transcription Start Sites). Com o auxílio da análise estatística, os pesquisadores podem analisar as leituras de dRNA-seq de bibliotecas tratadas e não tratadas com base na atividade da exonucleasse de RNA terminador dependente de 5'-monofosfato (TEX). Nosso estudo se concentra em Halobacterium salinarum NRC-1, um tipo de arqueia extremófila comumente encontrada em ambientes altamente salinos e é um organismo modelo para estudar a biologia molecular e genética em ambientes extremos. O objetivo de nosso estudo é identificar e mapear sítios alternativos de início de transcrição (alTSS, do inglês Alternative Transcription Start Sites), sítios de processamento de transcrição (TPS, do inglês Transcript Proces Sites) e sítios de terminação de transcrição (TTS, do inglês Transcription Termination Sites) usando os dados dRNA-seq em H. salinarum NRC-1. Modificamos a ferramenta TSSAR para detectar TSSs a partir de dados dRNA-seq, assumindo que as leituras de sequenciamento começam na posição exata durante a transcrição e seguem uma distribuição Gamma-Poisson. Anotamos os alTSSs em quatro tipos com base no número de leituras da biblioteca dRNA-seq e diferenças entre duas localizações principais de TSS. Os TSSs alternativos possuem leituras de RNA mais baixas do que as primárias e podem ter fase de leitura aberta upstream, levando a mudanças na produção de regulação gênica, isoformas 5'UTR e pausas na transcrição gênica. O mapeamento de alTSSs nos permitiu explicar as mudanças na resposta celular às condições de crescimento e expressão gênica em diferentes estágios de crescimento. Nossas descobertas revelaram um número significativo de sítios de início de transcrição internos (iTSSs, do inglês Internal Transcription Start Sites) erroneamente anotados que foram redefinidos como TSSs alternativos nas anotações genômicas anteriores. Esses TSSs alternativos produzem diferentes isoformas de proteína dependendo do comprimento da cadeia de aminoácidos e do quadro de leitura aberto. Além disso, foram identificados TSSs alternativos não apenas em H. salinarum, mas também em outros organismos, sugerindo um papel crucial na regulação da expressão gênica em várias espécies. Ademais, executamos uma reanálise dos dados de dRNA-seq, com ênfase em transcritos não primários (RNAs monofosforilados) ao invés do método tradicional de enriquecimento em transcritos primários (RNAs trifosforilados) para identificar sítios de processamento de transcrição (TPS, do inglês Transcript Proces Sites em todo o genoma em H. salinarum NRC-1. Também aplicamos essa abordagem ao Haloferax volcanii para análise comparativa. Por fim, utilizamos dados de dRNA-seq para identificar estruturas em hairpin e mapeá-las no genoma, fornecendo informações sobre o potencial papel dos locais de terminação da transcrição (TTS).