Abstract
Transcriptomic analyses across large scales of evolutionary distance have great potential to shed light on regulatory evolution but are complicated by difficulties in establishing orthology and limited availability of accessible software. We introduce here a method and a graphical user interface wrapper, called Annotator-RNAtor, for performing interspecies transcriptomic analysis and studying intragenus evolution. The pipeline uses third-party software to infer homologous genes in various species and highlight differences in the expression of the core-genes. To illustrate the methodology and demonstrate its usefulness, we focus on the emergence of the highly virulent Leptospira subclade known as P1+, which includes the causative agents of leptospirosis. Here, we expand on the genomic study through the comparison of transcriptomes between species from P1+ and their related P1- counterparts (low-virulent pathogens). In doing so, we shed light on differentially expressed pathways and focused on describing a specific example of adaptation based on a differential expression of PerRA-controlled genes. We showed that P1+ species exhibit higher expression of the katE gene, a well-known virulence determinant in pathogenic Leptospira species correlated with greater tolerance to peroxide. Switching PerRA alleles between P1+ and P1- species demonstrated that the lower repression of katE and greater tolerance to peroxide in P1+ species was solely controlled by PerRA and partly caused by a PerRA amino-acid permutation. Overall, these results demonstrate the strategic fit of the methodology and its ability to decipher adaptive transcriptomic changes, not observable by comparative genome analysis, that may have been implicated in the emergence of these pathogens.