Bart Cuypers

Metabolic adaptations of Leishmania donovani in relation to differentiation, drug resistance, and drug pressure

Antimonial (sodium stibogluconate, SSG) resistance and differentiation have been shown to be closely linked in Leishmania donovani, with SSG‐resistant strains showing an increased capacity to generate infectious (metacyclic) forms. This is the first untargeted LC‐MS metabolomics study which integrated both phenomena in one experimental design and provided insights into metabolic differences between three clinical L. donovani strains with a similar genetic background but different SSG‐susceptibilities. We performed this analysis at different stages during promastigote growth and in the absence or presence of drug pressure. When comparing SSG‐resistant and SSG‐sensitive strains, a number of metabolic changes appeared to be constitutively present in all growth stages, pointing towards a clear link with SSG‐resistance, whereas most metabolic changes were only detected in the stationary stage. These changes reflect the close intertwinement between SSG‐resistance and an increased metacyclogenesis in resistant parasites. The metabolic changes suggest that SSG‐resistant parasites have (i) an increased capacity for protection against oxidative stress; (ii) a higher fluidity of the plasma membrane; and (iii) a metabolic survival kit to better endure infection. These changes were even more pronounced in a resistant strain kept under SbIII drug pressure.

Genomic and Molecular Characterization of Miltefosine Resistance in Leishmania infantum Strains with Either Natural or Acquired Resistance through Experimental Selection of Intracellular Amastigotes.

During the last decade miltefosine (MIL) has been used as first-line treatment for visceral leishmaniasis in endemic areas with antimonial resistance, but a decline in clinical effectiveness is now being reported. While only two MIL-resistant Leishmania infantum strains from HIV co-infected patients have been documented, phenotypic MIL-resistance for L. donovani has not yet been identified in the laboratory. Hence, a better understanding of the factors contributing to increased MIL-treatment failure is necessary. Given the paucity of defined MIL-resistant L. donovani clinical isolates, this study used an experimental amastigote-selected MIL-resistant L. infantum isolate (LEM3323). In-depth exploration of the MIL-resistant phenotype was performed by coupling genomic with phenotypic data to gain insight into gene function and the mutant phenotype. A naturally MIL-resistant L. infantum clinical isolate (LEM5159) was included to compare both datasets. Phenotypically, resistance was evaluated by determining intracellular amastigote susceptibility in vitro and actual MIL-uptake. Genomic analysis provided supportive evidence that the resistance selection model on intracellular amastigotes can be a good proxy for the in vivo field situation since both resistant strains showed mutations in the same inward transporter system responsible for the acquired MIL-resistant phenotype. In line with previous literature findings in promastigotes, our data confirm a defective import machinery through inactivation of the LiMT/LiRos3 protein complex as the main mechanism for MIL-resistance also in intracellular amastigotes. Whole genome sequencing analysis of LEM3323 revealed a 2 base pair deletion in the LiMT gene that led to the formation an early stop codon and a truncation of the LiMT protein. Interestingly, LEM5159 revealed mutations in both the LiMT and LiRos3 genes, resulting in an aberrant expression of the LiMT protein. To verify that these mutations were indeed accountable for the acquired resistance, transfection experiments were performed to re-establish MIL-susceptibility. In LEM3323, susceptibility was restored upon expression of a LiMT wild-type gene, whereas the MIL-susceptibility of LEM5159 could be reversed after expression of the LiRos3 wild-type gene. The aberrant expression profile of the LiMT protein could be restored upon rescue of the LiRos3 gene both in the LEM5159 clinical isolate and a ΔLiRos3 strain, showing that expression of LdMT is dependent on LdRos3 expression. The present findings clearly corroborate the pivotal role of the LiMT/LiRos3 complex in resistance towards MIL.

Experimental Resistance to Drug Combinations in Leishmania donovani: Metabolic and Phenotypic Adaptations

ABSTRACT Together with vector control, chemotherapy is an essential tool for the control of visceral leishmaniasis (VL), but its efficacy is jeopardized by growing resistance and treatment failure against first-line drugs. To delay the emergence of resistance, the use of drug combinations of existing antileishmanial agents has been tested systematically in clinical trials for the treatment of visceral leishmaniasis (VL). In vitro, Leishmania donovani promastigotes are able to develop experimental resistance to several combinations of different antileishmanial drugs after 10 weeks of drug pressure. Using an untargeted liquid chromatography-mass spectrometry (LC-MS) metabolomics approach, we identified metabolic changes in lines that were experimentally resistant to drug combinations and their respective single-resistant lines. This highlighted both collective metabolic changes (found in all combination therapy-resistant [CTR] lines) and specific ones (found in certain CTR lines). We demonstrated that single-resistant and CTR parasite cell lines show distinct metabolic adaptations, which all converge on the same defensive mechanisms that were experimentally validated: protection against drug-induced and external oxidative stress and changes in membrane fluidity. The membrane fluidity changes were accompanied by changes in drug uptake only in the lines that were resistant against drug combinations with antimonials, and surprisingly, drug accumulation was higher in these lines. Together, these results highlight the importance and the central role of protection against oxidative stress in the different resistant lines. Ultimately, these phenotypic changes might interfere with the mode of action of all drugs that are currently used for the treatment of VL and should be taken into account in drug development.

LC-MS metabolomics from study design to data-analysis-using a versatile pathogen as a test case

Thanks to significant improvements in LC-MS technology, metabolomics is increasingly used as a tool to discriminate the responses of organisms to various stimuli or drugs. In this minireview we discuss all aspects of the LC-MS metabolomics pipeline, using a complex and versatile model organism, Leishmania donovani, as an illustrative example. The benefits of a hyphenated mass spectrometry platform and a detailed overview of the entire experimental pipeline from sampling, sample storage and sample list set-up to LC-MS measurements and the generation of meaningful results with state-of-the-art data-analysis software will be thoroughly discussed. Finally, we also highlight important pitfalls in the processing of LC-MS data and comment on the benefits of implementing metabolomics in a systems biology approach.

Modulation of Aneuploidy in Leishmania donovani during Adaptation to Different In Vitro and In Vivo Environments and Its Impact on Gene Expression

ABSTRACT Aneuploidy is usually deleterious in multicellular organisms but appears to be tolerated and potentially beneficial in unicellular organisms, including pathogens. Leishmania, a major protozoan parasite, is emerging as a new model for aneuploidy, since in vitro-cultivated strains are highly aneuploid, with interstrain diversity and intrastrain mosaicism. The alternation of two life stages in different environments (extracellular promastigotes and intracellular amastigotes) offers a unique opportunity to study the impact of environment on aneuploidy and gene expression. We sequenced the whole genomes and transcriptomes of Leishmania donovani strains throughout their adaptation to in vivo conditions mimicking natural vertebrate and invertebrate host environments. The nucleotide sequences were almost unchanged within a strain, in contrast to highly variable aneuploidy. Although high in promastigotes in vitro, aneuploidy dropped significantly in hamster amastigotes, in a progressive and strain-specific manner, accompanied by the emergence of new polysomies. After a passage through a sand fly, smaller yet consistent karyotype changes were detected. Changes in chromosome copy numbers were correlated with the corresponding transcript levels, but additional aneuploidy-independent regulation of gene expression was observed. This affected stage-specific gene expression, downregulation of the entire chromosome 31, and upregulation of gene arrays on chromosomes 5 and 8. Aneuploidy changes in Leishmania are probably adaptive and exploited to modulate the dosage and expression of specific genes; they are well tolerated, but additional mechanisms may exist to regulate the transcript levels of other genes located on aneuploid chromosomes. Our model should allow studies of the impact of aneuploidy on molecular adaptations and cellular fitness. IMPORTANCE Aneuploidy is usually detrimental in multicellular organisms, but in several microorganisms, it can be tolerated and even beneficial. Leishmania—a protozoan parasite that kills more than 30,000 people each year—is emerging as a new model for aneuploidy studies, as unexpectedly high levels of aneuploidy are found in clinical isolates. Leishmania lacks classical regulation of transcription at initiation through promoters, so aneuploidy could represent a major adaptive strategy of this parasite to modulate gene dosage in response to stressful environments. For the first time, we document the dynamics of aneuploidy throughout the life cycle of the parasite, in vitro and in vivo. We show its adaptive impact on transcription and its interaction with regulation. Besides offering a new model for aneuploidy studies, we show that further genomic studies should be done directly in clinical samples without parasite isolation and that adequate methods should be developed for this. IMPORTANCE Aneuploidy is usually detrimental in multicellular organisms, but in several microorganisms, it can be tolerated and even beneficial. Leishmania—a protozoan parasite that kills more than 30,000 people each year—is emerging as a new model for aneuploidy studies, as unexpectedly high levels of aneuploidy are found in clinical isolates. Leishmania lacks classical regulation of transcription at initiation through promoters, so aneuploidy could represent a major adaptive strategy of this parasite to modulate gene dosage in response to stressful environments. For the first time, we document the dynamics of aneuploidy throughout the life cycle of the parasite, in vitro and in vivo. We show its adaptive impact on transcription and its interaction with regulation. Besides offering a new model for aneuploidy studies, we show that further genomic studies should be done directly in clinical samples without parasite isolation and that adequate methods should be developed for this.

Genome-wide SNP analysis reveals distinct origins of Trypanosoma evansi and Trypanosoma equiperdum

Abstract Trypanosomes cause a variety of diseases in man and domestic animals in Africa, Latin America, and Asia. In the Trypanozoon subgenus, Trypanosoma brucei gambiense and Trypanosoma brucei rhodesiense cause human African trypanosomiasis, whereas Trypanosoma brucei brucei, Trypanosoma evansi, and Trypanosoma equiperdum are responsible for nagana, surra, and dourine in domestic animals, respectively. The genetic relationships between T. evansi and T. equiperdum and other Trypanozoon species remain unclear because the majority of phylogenetic analyses has been based on only a few genes. In this study, we have conducted a phylogenetic analysis based on genome-wide SNP analysis comprising 56 genomes from the Trypanozoon subgenus. Our data reveal that T. equiperdum has emerged at least once in Eastern Africa and T. evansi at two independent occasions in Western Africa. The genomes within the T. equiperdum and T. evansi monophyletic clusters show extremely little variation, probably due to the clonal spread linked to the independence from tsetse flies for their transmission.

Apolipoprotein L1 variant associated with increased susceptibility to trypanosome infection

ABSTRACT African trypanosomes, except Trypanosoma brucei gambiense and Trypanosoma brucei rhodesiense, which cause human African trypanosomiasis, are lysed by the human serum protein apolipoprotein L1 (ApoL1). These two subspecies can resist human ApoL1 because they express the serum resistance proteins T. b. gambiense glycoprotein (TgsGP) and serum resistance-associated protein (SRA), respectively. Whereas in T. b. rhodesiense, SRA is necessary and sufficient to inhibit ApoL1, in T. b. gambiense, TgsGP cannot protect against high ApoL1 uptake, so different additional mechanisms contribute to limit this uptake. Here we report a complex interplay between trypanosomes and an ApoL1 variant, revealing important insights into innate human immunity against these parasites. Using whole-genome sequencing, we characterized an atypical T. b. gambiense infection in a patient in Ghana. We show that the infecting trypanosome has diverged from the classical T. b. gambiense strains and lacks the TgsGP defense mechanism against human serum. By sequencing the ApoL1 gene of the patient and subsequent in vitro mutagenesis experiments, we demonstrate that a homozygous missense substitution (N264K) in the membrane-addressing domain of this ApoL1 variant knocks down the trypanolytic activity, allowing the trypanosome to avoid ApoL1-mediated immunity. IMPORTANCE Most African trypanosomes are lysed by the ApoL1 protein in human serum. Only the subspecies Trypanosoma b. gambiense and T. b. rhodesiense can resist lysis by ApoL1 because they express specific serum resistance proteins. We here report a complex interplay between trypanosomes and an ApoL1 variant characterized by a homozygous missense substitution (N264K) in the domain that we hypothesize interacts with the endolysosomal membranes of trypanosomes. The N264K substitution knocks down the lytic activity of ApoL1 against T. b. gambiense strains lacking the TgsGP defense mechanism and against T. b. rhodesiense if N264K is accompanied by additional substitutions in the SRA-interacting domain. Our data suggest that populations with high frequencies of the homozygous N264K ApoL1 variant may be at increased risk of contracting human African trypanosomiasis. Most African trypanosomes are lysed by the ApoL1 protein in human serum. Only the subspecies Trypanosoma b. gambiense and T. b. rhodesiense can resist lysis by ApoL1 because they express specific serum resistance proteins. We here report a complex interplay between trypanosomes and an ApoL1 variant characterized by a homozygous missense substitution (N264K) in the domain that we hypothesize interacts with the endolysosomal membranes of trypanosomes. The N264K substitution knocks down the lytic activity of ApoL1 against T. b. gambiense strains lacking the TgsGP defense mechanism and against T. b. rhodesiense if N264K is accompanied by additional substitutions in the SRA-interacting domain. Our data suggest that populations with high frequencies of the homozygous N264K ApoL1 variant may be at increased risk of contracting human African trypanosomiasis.

Macromolecular biosynthetic parameters and metabolic profile in different life stages of Leishmania braziliensis: Amastigotes as a functionally less active stage

It was recently hypothesized that Leishmania amastigotes could constitute a semi-quiescent stage characterized by low replication and reduced metabolic activity. This concept developed with Leishmania (Leishmania) mexicana and Leishmania (Leishmania) major models might explain numerous clinical and sub-clinical features of Leishmania (Viannia) braziliensis infections, like reactivation of the disease, non-response to chemotherapy or asymptomatic infections. We compared here in vitro the proliferative capability of L. (V.) braziliensis amastigotes and promastigotes, assessed the expression of key molecular parameters and performed metabolomic analysis. We found that contrary to the highly proliferative promastigotes, amastigotes (axenic and intracellular) do not show evidence of extensive proliferation. In parallel, amastigotes showed a significant decrease of (i) the kDNA mini-circle abundance, (ii) the intracellular ATP level, (iii) the ribosomal components: rRNA subunits 18S and 28S α and ribosomal proteins RPS15 and RPL19, (iv) total RNA and protein levels. An untargeted metabolomic study identified clear differences between the different life stages: in comparison to logarithmic promastigotes, axenic amastigotes showed (a) a strong decrease of 14 essential and non-essential amino acids and eight metabolites involved in polyamine synthesis, (b) extensive changes in the phospholipids composition and (c) increased levels of several endogenous and exogenous sterols. Altogether, our results show that L. (V.) braziliensis amastigotes can show a phenotype with negligible rate of proliferation, a lower capacity of biosynthesis, a reduced bio-energetic level and a strongly altered metabolism. Our results pave the way for further exploration of quiescence among amastigotes of this species.

On the feasibility of mining CD8+ T cell receptor patterns underlying immunogenic peptide recognition

Current T cell epitope prediction tools are a valuable resource in designing targeted immunogenicity experiments. They typically focus on, and are able to, accurately predict peptide binding and presentation by major histocompatibility complex (MHC) molecules on the surface of antigen-presenting cells. However, recognition of the peptide-MHC complex by a T cell receptor (TCR) is often not included in these tools. We developed a classification approach based on random forest classifiers to predict recognition of a peptide by a T cell receptor and discover patterns that contribute to recognition. We considered two approaches to solve this problem: (1) distinguishing between two sets of TCRs that each bind to a known peptide and (2) retrieving TCRs that bind to a given peptide from a large pool of TCRs. Evaluation of the models on two HIV-1, B*08-restricted epitopes reveals good performance and hints towards structural CDR3 features that can determine peptide immunogenicity. These results are of particular importance as they show that prediction of T cell epitope and T cell epitope recognition based on sequence data is a feasible approach. In addition, the validity of our models not only serves as a proof of concept for the prediction of immunogenic T cell epitopes but also paves the way for more general and high-performing models.

Toxicogenomics of the flame retardant tris (2-butoxyethyl) phosphate in HepG2 cells using RNA-seq

Tris (2-butoxyethyl) phosphate (TBOEP) is a compound produced at high volume that is used as both a flame retardant and a plasticizer. It is persistent and bioaccumulative, yet little is known of its toxicological modes of action. Such insight may aid risk assessment in a weight-of-evidence approach supplementing current testing strategies. We used an RNA sequencing approach as an unbiased and sensitive tool to explore potential negative health effects of sub-cytotoxic concentrations of TBOEP on the transcriptome of the human liver hepatocellular carcinoma cell line, HepG2, with the lowest concentration used potentially holding relevance to human physiological levels. Over-representation and gene set enrichment analysis corresponded well and revealed that TBOEP treatments resulted in an upregulation of genes involved in protein and energy metabolism, along with DNA replication. Such increases in cell and macromolecule metabolism could explain the increase in mitochondrial activity at lower TBOEP concentrations. In addition, TBOEP affected a wide variety of biological processes, the most notable one being the general stress response, wound healing. Finally, TBOEP showed effects on steroid hormone biosynthesis and activation, regulation, and potentiation of immune responses, in agreement with other studies. As such, this study is the first study investigating genome-wide changes in gene transcription in response to TBOEP in human cells.