Barbara Papadopoulou

Leishmania relies on posttranscriptional control to regulate gene expression. Protein-coding genes are synthesised as polycistronic precursors that are processed into individual mRNAs by trans-splicing adding the spliced leader (SL) RNA to the 5’-end and 3’ cleavage-polyadenylation. Here, we employ Nanopore direct RNA sequencing (DRS) combined with Illumina RNA-Seq to comprehensively interrogate the transcriptomes of Leishmania infantum developmental stages at single-molecule resolution. Analysis of DRS full-length reads of poly(A)+-enriched RNA from L. infantum developmental stages enabled us to precisely determine the primary SL and poly(A) sites for 52% of the protein-coding transcripts and to accurately define their 5’- and 3’-end and the length of UTRs. In addition, our analysis confirmed the motifs ‘[C/A/T] A|G’ being associated with 94.8% of the SL cleavage sites and better defined the genomic context for cleavage and polyadenylation. Overall, we observed more diversity for poly(A) than SL sites per transcript. The frequency of the primary SL and poly(A) sites was 64.2% and 24%, respectively, with most transcripts having additional poly(A) sites nearby. Alternative polyadenylation was detected in 11-13% of transcripts with ~ 20% of these having different primary poly(A) sites between promastigote and amastigote developmental stages. Furthermore, DRS uncovered multiple processing events occurring mostly within 3’UTRs, leading to the formation of long non-coding RNAs (lncRNAs). The L. infantum transcriptome expresses a rich repertoire of 1,825 lncRNAs, of which 98% were not previously annotated in L. infantum and only 21.5% were found in L. major. These lncRNAs exhibit generally distinct expression patterns from the 3’UTRs they derived and several are developmentally regulated, representing ~ 27% of the L. infantum stage-regulated transcriptome. Their expression was generally higher in amastigotes than in promastigotes, highlighting their importance in parasite intracellular development. Protein prediction tools combined to mass-spectrometry revealed that 7.6% of these lncRNAs have a limited protein-coding potential. This is the first comprehensive transcriptomic analysis of L. infantum developmental stages using single-molecule Nanopore DRS. Our findings advance knowledge on existing Leishmania expression datasets and provide new insights into the transcriptome complexity and dynamics of both protein-coding and non-coding sequences throughout the parasite development.

Background: Leishmania relies on posttranscriptional control to regulate gene expression. Protein-coding genes are synthesised as polycistronic precursors that are processed into individual mRNAs by coupled trans-splicing adding the spliced leader (SL) RNA to the 5'-end of all transcripts and 3' cleavage-polyadenylation. Here, we employ Nanopore direct RNA sequencing (DRS) combined with Illumina RNA-Seq to comprehensively interrogate the transcriptomes of Leishmania infantum developmental stages at single-molecule resolution. Results: Analysis of DRS full-length reads of poly(A)+-enriched RNA from L. infantum developmental stages enabled us to precisely determine the primary SL and poly(A) sites for 52% of the protein-coding transcripts and to accurately define their 5'- and 3'-end boundaries and the length of UTRs. We observed a higher diversity for the poly(A) sites per transcript than the SL sites. The frequency of the primary SL site was ~ 64.2%, whereas the incidence of the primary poly(A) site was 24%, with most transcripts having additional poly(A) sites nearby. Alternative polyadenylation was detected in 11–13% of transcripts, of which 28% were developmentally regulated. Similarly, 26.5% of transcripts with different primary SL sites between developmental stages were differentially expressed. Our analysis confirmed the motifs ‘[C/A/T] A|G’ being associated with 94.8% of the SL cleavage sites and better defined the genomic context for cleavage and polyadenylation. Further, DRS uncovered multiple processing events occurring mostly within long 3'UTRs, leading to the formation of long non-coding RNAs (lncRNAs). The L. infantum transcriptome expresses a rich repertoire of 1,825 lncRNAs, of which 98% were not previously annotated in L. infantum and most are absent from L. major and may be specific to L. infantum. Generally, these lncRNAs exhibit distinct expression patterns from the 3'UTRs they derived, and several are developmentally regulated, representing 27% of the L. infantum stage-regulated transcriptome. Protein prediction tools combined to mass-spectrometry revealed that 7.6% of the lncRNAs have a limited protein-coding potential. Conclusions: This is the first comprehensive transcriptomic analysis of L. infantum developmental stages using single-molecule Nanopore DRS. Our findings advance knowledge on existing Leishmania expression datasets and provide new insights into the transcriptome complexity and dynamics of protein-coding and non-coding sequences throughout the parasite development.

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.

Poly(A) Binding Proteins (PABPs) are major eukaryotic RNA-binding proteins (RBPs) with multiple roles associated with mRNA stability and translation and characterized mainly from multicellular organisms and yeasts. A variable number of PABP homologues are seen in different organisms however the biological reasons for multiple PABPs are generally not well understood. In the unicellular Leishmania, dependent on post-transcriptional mechanisms for the control of its gene expression, three distinct PABPs are found, with yet undefined functional distinctions. Here, using RNA-immunoprecipitation sequencing analysis we show that the Leishmania PABP1 preferentially associates with mRNAs encoding ribosomal proteins, while PABP2 and PABP3 bind to an overlapping set of mRNAs distinct to those enriched in PABP1. Immunoprecipitation studies combined to mass-spectrometry analysis identified RBPs differentially associated with PABP1 or PABP2, including RBP23 and DRBD2, respectively, that were investigated further. Both RBP23 and DRBD2 bind directly to the three PABPs in vitro, but reciprocal experiments confirmed preferential co-immunoprecipitation of PABP1, as well as the EIF4E4/EIF4G3 based translation initiation complex, with RBP23. Other RBP23 binding partners also imply a direct role in translation. DRBD2, in contrast, co-immunoprecipitated with PABP2, PABP3 and with RBPs unrelated to translation. Over 90% of the RBP23-bound mRNAs code for ribosomal proteins, mainly absent from the transcripts co-precipitated with DRBD2. These experiments suggest a novel and specific route for translation of the ribosomal protein mRNAs, mediated by RBP23, PABP1 and the associated EIF4E4/EIF4G3 complex. They also highlight the unique roles that different PABP homologues may have in eukaryotic cells associated with mRNA translation.

DDX3 is a multifaceted RNA helicase of the DEAD-box family that plays central roles in all aspects of RNA metabolism including translation initiation. Here, we provide evidence that the Leishmania DDX3 ortholog functions in post-initiation steps of translation. We show that genetic depletion of DDX3 slows down ribosome movement resulting in elongation-stalled ribosomes, impaired translation elongation and decreased de novo protein synthesis. We also demonstrate that the essential ribosome recycling factor Rli1/ABCE1 and termination factors eRF3 and GTPBP1 are less recruited to ribosomes upon DDX3 loss, suggesting that arrested ribosomes may be inefficiently dissociated and recycled. Furthermore, we show that prolonged ribosome stalling triggers co-translational ubiquitination of nascent polypeptide chains and a higher recruitment of E3 ubiquitin ligases and proteasome components to ribosomes of DDX3 knockout cells, which further supports that ribosomes are not elongating optimally. Impaired elongation of translating ribosomes also results in the accumulation of cytoplasmic protein aggregates, which implies that defects in translation overwhelm the normal quality controls. The partial recovery of translation by overexpressing Hsp70 supports this possibility. Collectively, these results suggest an important novel contribution of DDX3 to optimal elongation of translating ribosomes by preventing prolonged translation stalls and stimulating recycling of arrested ribosomes.

Valosin‐containing protein (VCP)/p97/Cdc48 is an AAA + ATPase associated with many ubiquitin-dependent cellular pathways that are central to protein quality control. VCP binds various cofactors, which determine pathway selectivity and substrate processing. Here, we used co-immunoprecipitation and mass spectrometry studies coupled to in silico analyses to identify the Leishmania infantum VCP (LiVCP) interactome and to predict molecular interactions between LiVCP and its major cofactors. Our data support a largely conserved VCP protein network in Leishmania including known but also novel interaction partners. Network proteomics analysis confirmed LiVCP-cofactor interactions and provided novel insights into cofactor-specific partners and the diversity of LiVCP complexes, including the well-characterized VCP-UFD1-NPL4 complex. Gene Ontology analysis coupled with digitonin fractionation and immunofluorescence studies support cofactor subcellular compartmentalization with either cytoplasmic or organellar or vacuolar localization. Furthermore, in silico models based on 3D homology modeling and protein–protein docking indicated that the conserved binding modules of LiVCP cofactors, except for NPL4, interact with specific binding sites in the hexameric LiVCP protein, similarly to their eukaryotic orthologs. Altogether, these results allowed us to build the first VCP protein interaction network in parasitic protozoa through the identification of known and novel interacting partners potentially associated with distinct VCP complexes.

RNA-binding proteins (RBPs) play key roles in many aspects of RNA metabolism. In Leishmania, a unicellular eukaryote that favors the posttranscriptional mode of regulation for controlling gene expression levels, the function of RBPs becomes even more critical. However, due largely to limited in vivo approaches available for identifying RBPs in these parasites, there have been no significant advances to our understanding of the role these proteins play in posttranscriptional control through binding to cis-acting elements in the 3′ untranslated region (3′UTR) of mRNAs. Here we describe an optimized in vivo RNA tethering approach using the bacteriophage MS2 coat protein combined to immunoprecipitation and mass spectrometry analysis to identify RBPs specifically interacting with 3′UTR short interspersed degenerated retroposon elements (SIDERs). Members of the SIDER2 subfamily were shown previously to promote mRNA degradation through a novel mechanism of mRNA decay. Using this modified MS2 tethering approach, we have identified candidate RBPs specifically interacting with SIDER2 elements and contributing to the decay mechanism.

VCP/p97/Cdc48 is one of the best‐characterized type II cytosolic AAA+ ATPases most known for their role in ubiquitin‐dependent protein quality control. Here, we provide functional insights into the role of the Leishmania VCP/p97 homolog (LiVCP) in the parasite intracellular development. We demonstrate that although LiVCP is an essential gene, L. infantum promastigotes can grow with less VCP. In contrast, growth of axenic and intracellular amastigotes is dramatically affected upon decreased LiVCP levels in heterozygous and temperature sensitive LiVCP mutants or the expression of dominant negative mutants known to specifically target the second conserved VCP ATPase domain, a major contributor of the VCP overall ATPase activity. Interestingly, these VCP mutants are also unable to survive heat stress and a temperature sensitive VCP mutant is defective in amastigote growth. Consistent with LiVCP's essential function in amastigotes, LiVCP mRNA undergoes 3'UTR‐mediated developmental regulation, resulting in higher VCP expression in amastigotes. Furthermore, we show that parasite mutant lines expressing lower VCP levels or dominant negative VCP forms exhibit high accumulation of polyubiquitinated proteins and increased sensitivity to proteotoxic stress, supporting the ubiquitin‐selective chaperone function of LiVCP. Together, these results emphasize the crucial role LiVCP plays under heat stress and during the parasite intracellular development.

The Poly-A Binding Protein (PABP) is a conserved eukaryotic polypeptide involved in many aspects of mRNA metabolism. During translation initiation, PABP interacts with the translation initiation complex eIF4F and enhances the translation of polyadenylated mRNAs. Schematically, most PABPs can be divided into an N-terminal RNA-binding region, a non-conserved linker segment and the C-terminal MLLE domain. In pathogenic Leishmania protozoans, three PABP homologues have been identified, with the first one (PABP1) targeted by phosphorylation and shown to co-immunoprecipitate with an eIF4F-like complex (EIF4E4/EIF4G3) implicated in translation initiation. Here, PABP1 phosphorylation was shown to be linked to logarithmic cell growth, reminiscent of EIF4E4 phosphorylation, and coincides with polysomal association. Phosphorylation targets multiple serine-proline (SP) or threonine-proline (TP) residues within the PABP1 linker region. This is an essential protein, but phosphorylation is not needed for its association with polysomes or cell viability. Mutations which do impair PABP1 polysomal association and are required for viability do not prevent phosphorylation, although further mutations lead to a presumed inactive protein largely lacking phosphorylated isoforms. Co-immunoprecipitation experiments were carried out to investigate PABP1 function further, identifying several novel protein partners and the EIF4E4/EIF4G3 complex, but no other eIF4F-like complex or subunit. A novel, direct interaction between PABP1 and EIF4E4 was also investigated and found to be mediated by the PABP1 MLLE binding to PABP Interacting Motifs (PAM2) within the EIF4E4 N-terminus. The results shown here are consistent with phosphorylation of PABP1 being part of a novel pathway controlling its function and possibly translation in Leishmania.

Reevaluation of treatment guidelines for Old and New World leishmaniasis is urgently needed on a global basis because treatment failure is an increasing problem. Drug resistance is a fundamental determinant of treatment failure, although other factors also contribute to this phenomenon, including the global HIV/AIDS epidemic with its accompanying impact on the immune system. Pentavalent antimonials have been used successfully worldwide for the treatment of leishmaniasis since the first half of the 20th century, but the last 10 to 20 years have witnessed an increase in clinical resistance, e.g., in North Bihar in India. In this review, we discuss the meaning of “resistance” related to leishmaniasis and discuss its molecular epidemiology, particularly for Leishmania donovani that causes visceral leishmaniasis. We also discuss how resistance can affect drug combination therapies. Molecular mechanisms known to contribute to resistance to antimonials, amphotericin B, and miltefosine are also outlined. Chemotherapy is central to the control and management of leishmaniasis. Antimonials remain the primary drugs against different forms of leishmaniasis in several regions. However, resistance to antimony has necessitated the use of alternative medications, especially in the Indian subcontinent (ISC). Compounds, notably the orally available miltefosine (MIL), parenteral paromomycin, and amphotericin B (AmB), are increasingly used to treat leishmaniasis. Although treatment failure (TF) has been observed in patients treated with most anti-leishmanials, its frequency of appearance may be important in patients treated with MIL, which has replaced antimonials within the kala-azar elimination program in the ISC. AmB is highly efficacious, and the associated toxic effects—when administered in its free deoxycholate form—are somewhat ameliorated in its liposomal formulation. Regrettably, laboratory experimentation has demonstrated a risk of resistance towards AmB as well. The rise of drug resistance impacts treatment outcome, and understanding its causes, spread, and impact will help us manage the risks it imposes. Here, we review the problem of TF in leishmaniasis and the contribution of drug resistance to the problem. Molecular mechanisms causing resistance to anti-leishmanials are discussed along with the appropriate use of additional available drugs, as well as the urgent need to consolidate strategies to monitor drug efficacy, epidemiological surveillance, and local policies. Coordination of these activities in national and international programs against leishmaniasis might represent a successful guide to further research and prevention activities.