Evolutionary divergent clusters of transcribed extinct truncated retroposons drive low mRNA expression and developmental regulation in the protozoan Leishmania

BACKGROUND

The Leishmania genome harbors formerly active short interspersed degenerated retroposons (SIDERs) representing the largest family of repetitive elements among trypanosomatids. Their substantial expansion in Leishmania is a strong predictor of important biological functions. In this study, we combined multilevel bioinformatic predictions with high-throughput genomic and transcriptomic analyses to gain novel insights into the diversified roles retroposons of the SIDER2 subfamily play in Leishmania genome evolution and expression.

RESULTS

We show that SIDER2 retroposons form various evolutionary divergent clusters, each harboring homologous SIDER2 sequences usually located nearby in the linear sequence of chromosomes. This intriguing genomic organization underscores the importance of SIDER2 proximity in shaping chromosome dynamics and co-regulation. Accordingly, we show that transcripts belonging to the same SIDER2 cluster can display similar levels of expression. SIDER2 retroposons are mostly transcribed as part of 3'UTRs and account for 13% of the Leishmania transcriptome. Genome-wide expression profiling studies underscore SIDER2 association generally with low mRNA expression. The remarkable link of SIDER2 retroposons with downregulation of gene expression supports their co-option as major regulators of mRNA abundance. SIDER2 sequences also add to the diversification of the Leishmania gene expression repertoire since ~ 35% of SIDER2-containing transcripts can be differentially regulated throughout the parasite development, with a few encoding key virulence factors. In addition, we provide evidence for a functional bias of SIDER2-containing transcripts with protein kinase and transmembrane transporter activities being most represented.

CONCLUSIONS

Altogether, these findings provide important conceptual advances into evolutionary innovations of transcribed extinct retroposons acting as major RNA cis-regulators.

In silico analysis of Leishmania proteomes and protein-protein interaction network: Prioritizing therapeutic targets and drugs for repurposing to treat leishmaniasis

Leishmaniasis is a serious world health problem and its current therapies have several limitations demanding to develop novel therapeutics for this disease. The present study aims to prioritize novel broad-spectrum targets using proteomics and protein-protein interaction network (PPIN) data for 11 Leishmania species. Proteome comparison and host non-homology analysis resulted in 3605 pathogen-specific conserved core proteins. Gene ontology analysis indicated their involvement in major molecular functions like DNA binding, transportation, dioxygenase, and catalytic activity. PPIN analysis of these core proteins identified eight hub proteins (viz., vesicle-trafficking protein (LBRM2903_190011800), ribosomal proteins S17 (LBRM2903_34004790) and L2 (LBRM2903_080008100), eukaryotic translation initiation factor 3 (LBRM2903_350086700), replication factor A (LBRM2903_150008000), U3 small nucleolar RNA-associated protein (LBRM2903_340025600), exonuclease (LBRM2903_200021800), and mitochondrial RNA ligase (LBRM2903_200074100)). Among the hub proteins, six were classified as drug targets and two as vaccine candidates. Further, druggability analysis indicated three hub proteins, namely eukaryotic translation initiation factor 3, ribosomal proteins S17 and L2 as druggable. Their three-dimensional structures were modelled and docked with the identified ligands (2-methylthio-N6-isopentenyl-adenosine-5'-monophosphate, artenimol and omacetaxine mepesuccinate). These ligands could be experimentally validated (in vitro and in vivo) and repurposed for the development of novel antileishmanial agents.

Unraveling structural insights of ribokinase from Leishmania donovani

Ribokinase (RK) is an ATP dependent sugar kinase that enables the entry of ribose in the metabolism. Leishmania accumulates ribose into the cytosol through hydrolysis of nucleosides and by transport from the extracellular environment. Activation by RK is critical to mobilize the ribose into the metabolism of Leishmania. To understand the catalytic role, the crystal structure of RK (LdRK) from L. donovani was determined in the apo and complex forms with several nucleotides (ATP, AMPPCP and ADP) in the presence of Na⁺ ion. The dual insertion of five amino acid stretches makes LdRK structurally unique from other reported structures of RKs. The structure of LdRK-ATP provided the basis for positioning of γ-phosphate of ATP by conserved -GAGD- motif. Liganded and unliganded structures of LdRK exists in similar conformation, which suggests binding of nucleotides does not make any significant conformational changes in nucleotide-bound structures. Substitution of a conserved asparagine with phenylalanine in ribose binding pocket differentiates the LdRK from other RKs. Glycerol molecule bound in the substrate binding pocket mimics the enzyme-substrate interactions but in turn, hampers the binding of ribose to LdRK. Comparative structural analysis revealed the flexibility of γ-phosphate, which adopts multiple conformations in the absence of divalent metal ion and ribose. Similar to other RKs, LdRK is also dependent on monovalent as well as divalent cations for its catalytic activity.