Search

Abstract A large proportion of HIV-coinfected visceral leishmaniasis (VL-HIV) patients exhibit chronic disease with frequent VL recurrence. However, knowledge on immunological determinants underlying the disease course is scarce. We longitudinally profiled the circulatory cellular immunity of an Ethiopian HIV cohort that included VL developers. We show that chronic VL-HIV patients exhibit high and persistent levels of TIGIT and PD-1 on CD8+/CD8- T cells, in addition to a lower frequency of IFN-γ+ TIGIT- CD8+/CD8- T cells, suggestive of impaired T cell functionality. At single T cell transcriptome and clonal resolution, the patients show CD4+ T cell anergy, characterised by a lack of T cell activation and lymphoproliferative response. These findings suggest that PD-1 and TIGIT play a pivotal role in VL-HIV chronicity, and may be further explored for patient risk stratification. Our findings provide a strong rationale for adjunctive immunotherapy for the treatment of chronic VL-HIV patients to break the recurrent disease cycle.

Abstract Background Plasmodium vivax is the second most important cause of human malaria worldwide, and accounts for the majority of malaria cases in South America. A high-quality reference genome exists for Papua Indonesia (PvP01) and Thailand (PvW1), but is lacking for South America. A reference genome specifically for South America would be beneficial though, as P. vivax is a genetically diverse parasite with geographical clustering. Results This study presents a new high-quality assembly of a South American P. vivax isolate, referred to as PvPAM (P. vivax Peruvian AMazon). The genome was obtained from a low input patient sample from the Peruvian Amazon and sequenced using PacBio technology, resulting in a highly complete assembly with 6497 functional genes. Telomeric ends were present in 17 out of 28 chromosomal ends, and additional (sub)telomeric regions are present in 12 unassigned contigs. A comparison of multigene families between PvPAM and the PvP01 genome revealed remarkable variation in vir genes, and the presence of merozoite surface proteins (MSP) 3.6 and 3.7. Three dhfr and dhps drug resistance associated mutations are present in PvPAM, similar to those found in other Peruvian isolates. Mapping of publicly available South American whole genome sequencing (WGS) data to PvPAM resulted in significantly fewer variants and truncated reads compared to the use of PvP01 or PvW1 as reference genomes. To minimize the number of core genome variants in non-South American samples, PvW1 is most suited for Southeast Asian isolates, both PvPAM and PvW1 are suited for South Asian isolates, and PvPAM is recommended for African isolates. Interestingly, non-South American samples still contained the least subtelomeric variants when mapped to PvPAM, indicating high quality of the PvPAM subtelomeric regions. Conclusions Our findings show that the PvPAM reference genome more accurately represents South American P. vivax isolates in comparison to PvP01 and PvW1. In addition, PvPAM has a high level of completeness, and contains a similar number of annotated genes as PvP01 or PvW1. The PvPAM genome therefore will be a valuable resource to improve future genomic analyses on P. vivax isolates from the South American continent.

ABSTRACTWhole-genome sequencing of Plasmodium is becoming an increasingly important tool for genomic surveillance of malaria. Due to the predominance of human DNA in a patient blood sample, time-consuming laboratory procedures are required to deplete human DNA or enrich Plasmodium DNA. Here, we investigated the potential of nanopore adaptive sampling to enrich Plasmodium falciparum reads while sequencing unenriched patient blood samples. To compare adaptive sampling versus regular sequencing on a MinION device, a dilution series consisting of 0%–84% P. falciparum DNA in human DNA was sequenced. Half of the flow cell channels were run in adaptive sampling mode, enriching for the P. falciparum reference genome, resulting in a three- to five-fold enrichment of P. falciparum bases in samples containing 0.1%–8.4% P. falciparum DNA. This finding was confirmed by sequencing three P. falciparum patient blood samples with common levels of parasitemia, that is, 0.1%, 0.2%, and 0.6% in adaptive mode. Their estimated enrichment was 5.8, 3.9, and 2.7, respectively, which was sufficient to cover at least 97% of the P. falciparum reference genome at a median depth of 5 (lowest parasitemia) or 355 (highest parasitemia). In all, 38 drug resistance loci were compared to Sanger sequencing results, showing high concordance, which suggests that the obtained sequencing data are of sufficient quality to address common clinical research questions for patients with parasitemias of 0.1% and higher. Overall, our results indicate that adaptive nanopore sequencing has the potential to replace more time-consuming Plasmodium enrichment protocols in the future.IMPORTANCEMalaria is caused by parasites of the genus Plasmodium, and reached a global disease burden of 247 million cases in 2021. To study drug resistance mutations and parasite population dynamics, whole-genome sequencing of patient blood samples is commonly performed. However, the predominance of human DNA in these samples imposes the need for time-consuming laboratory procedures to enrich Plasmodium DNA. We used the Oxford Nanopore Technologies’ adaptive sampling feature to circumvent this problem and enrich Plasmodium reads directly during the sequencing run. We demonstrate that adaptive nanopore sequencing efficiently enriches Plasmodium reads, which simplifies and shortens the timeline from blood collection to parasite sequencing. In addition, we show that the obtained data can be used for monitoring genetic markers, or to generate nearly complete genomes. Finally, owing to its inherent mobility, this technology can be easily applied on-site in endemic areas where patients would benefit the most from genomic surveillance.

Radiotherapy is the standard of care for breast cancer. However, surviving radioresistant cells can repopulate following treatment and provoke relapse. Better understanding of the molecular mechanisms of radiation resistance may help to improve treatment of radioresistant tumours. To emulate radiation therapy at the cellular level, we exposed MCF7 breast cancer cells to daily radiation doses of 2 Gy up to an accumulated dose of 20 Gy. Fractionally irradiated cells (FIR20) displayed increased clonogenic survival and population doubling time as compared with age‐matched sham‐irradiated cells and untreated parental MCF7 cells. RNA‐sequencing revealed a core signature of 229 mRNAs and 7 circular RNAs of which the expression was significantly altered in FIR20 cells. Dysregulation of several top genes was mirrored at the protein level. The FIR20 cell transcriptome overlapped significantly with canonical radiation response signatures and demonstrated a remarkable commonality with radiation and endocrine therapy resistance expression profiles, suggesting crosstalk between both acquired resistance pathways, as indicated by reduced sensitivity to tamoxifen cytotoxicity of FIR20 cells. Using predictive analyses and functional enrichment, we identified a gene‐regulatory network that promotes stemness and inflammatory signalling in FIR20 cells. We propose that these phenotypic traits render breast cancer cells more radioresistant but may at the same time serve as potential targets for combination therapies.

Aneuploidy causes system-wide disruptions in the stochiometric balances of transcripts, proteins, and metabolites, often resulting in detrimental effects for the organism. The protozoan parasite Leishmania has an unusually high tolerance for aneuploidy, but the molecular and functional consequences for the pathogen remain poorly understood. Here, we addressed this question in vitro and present the first integrated analysis of the genome, transcriptome, proteome, and metabolome of highly aneuploid Leishmania donovani strains. Our analyses unambiguously establish that aneuploidy in Leishmania proportionally impacts the average transcript- and protein abundance levels of affected chromosomes, ultimately correlating with the degree of metabolic differences between closely related aneuploid strains. This proportionality was present in both proliferative and non-proliferative in vitro promastigotes. However, as in other Eukaryotes, we observed attenuation of dosage effects for protein complex subunits and in addition, non-cytoplasmic proteins. Differentially expressed transcripts and proteins between aneuploid Leishmania strains also originated from non-aneuploid chromosomes. At protein level, these were enriched for proteins involved in protein metabolism, such as chaperones and chaperonins, peptidases, and heat-shock proteins. In conclusion, our results further support the view that aneuploidy in Leishmania can be adaptive. Additionally, we believe that the high karyotype diversity in vitro and absence of classical transcriptional regulation make Leishmania an attractive model to study processes of protein homeostasis in the context of aneuploidy and beyond.

The Plasmodium vivax reticulocyte invasion process is still poorly understood, with only a few receptor-ligand interactions identified to date. Individuals with the Southeast Asian ovalocytosis (SAO) phenotype have a deletion in the band 3 protein on the surface of erythrocytes, and are reported to have a lower incidence of clinical P. vivax malaria. Based on this observation, band 3 has been put forward as a receptor for P. vivax invasion, although direct proof is still lacking. In this study, we combined functional ex vivo invasion assays and transcriptome sequencing to uncover a band 3–mediated invasion pathway in P. vivax and potential band 3 ligands. Invasion by P. vivax field isolates was 67%-71% lower in SAO reticulocytes compared with non-SAO reticulocytes. Reticulocyte invasion was decreased by 40% and 27%-31% when blocking with an anti-band 3 polyclonal antibody and a PvTRAg38 peptide, respectively. To identify new band 3 receptor candidates, we mRNA-sequenced schizont-stage isolates used in the invasion assays, and observed high transcriptional variability in multigene and invasion-related families. Transcriptomes of isolates with low or high dependency on band 3 for invasion were compared by differential expression analysis, which produced a list of band 3 ligand candidates with high representation of PvTRAg genes. Our ex vivo invasion assays have demonstrated that band 3 is a P. vivax invasion receptor and confirm previous in vitro studies showing binding between PvTRAg38 and band 3, although the lower and variable inhibition levels observed suggest the involvement of other ligands. By coupling transcriptomes and invasion phenotypes from the same isolates, we identified a list of band 3 ligand candidates, of which the overrepresented PvTRAg genes are the most promising for future research.

Abstract Melanoma remains a deadly cancer despite significant advances in immune checkpoint blockade and targeted therapies. The incidence of melanoma is also growing worldwide, which highlights the need for novel treatment options and strategic combination of therapies. Here, we investigate non‐thermal plasma (NTP), an ionized gas, as a promising, therapeutic option. In a melanoma mouse model, direct treatment of tumors with NTP results in reduced tumor burden and prolonged survival. Physical characterization of NTP treatment in situ reveals the deposited NTP energy and temperature associated with therapy response, and whole transcriptome analysis of the tumor identified several modulated pathways. NTP treatment also enhances the cancer‐immunity cycle, as immune cells in both the tumor and tumor‐draining lymph nodes appear more stimulated to perform their anti‐cancer functions. Thus, our data suggest that local NTP therapy stimulates systemic, anti‐cancer immunity. We discuss, in detail, how these fundamental insights will help direct the translation of NTP technology into the clinic and inform rational combination strategies to address the challenges in melanoma therapy.

Overexpression of the transmembrane mucin MUC13, as seen in inflammatory bowel diseases (IBD), could potentially impact barrier function. This study aimed to explore how inflammation-induced MUC13 disrupts epithelial barrier integrity by affecting junctional protein expression in IBD, thereby also considering the involvement of MUC1. RNA sequencing and permeability assays were performed using LS513 cells transfected with MUC1 and MUC13 siRNA and subsequently stimulated with IL-22. In vivo intestinal permeability and MUC13-related signaling pathways affecting barrier function were investigated in acute and chronic DSS-induced colitis wildtype and Muc13−/− mice. Finally, the expression of MUC13, its regulators and other barrier mediators were studied in IBD and control patients. Mucin knockdown in intestinal epithelial cells affected gene expression of several barrier mediators in the presence/absence of inflammation. IL-22-induced MUC13 expression impacted barrier function by modulating the JAK1/STAT3, SNAI1/ZEB1 and ROCK2/MAPK signaling pathways, with a cooperating role for MUC1. In response to DSS, MUC13 was protective during the acute phase whereas it caused more harm upon chronic colitis. The pathways accounting for the MUC13-mediated barrier dysfunction were also altered upon inflammation in IBD patients. These novel findings indicate an active role for aberrant MUC13 signaling inducing intestinal barrier dysfunction upon inflammation with MUC1 as collaborating partner.

ABSTRACT The implementation of prospective drug resistance (DR) studies in the research-and-development (R&D) pipeline is a common practice for many infectious diseases but not for neglected tropical diseases (NTDs). Here, we explored and demonstrated the importance of this approach using as paradigms Leishmania donovani, the etiological agent of visceral leishmaniasis (VL), and TCMDC-143345, a promising compound of the GlaxoSmithKline (GSK) “Leishbox” to treat VL. We experimentally selected resistance to TCMDC-143345 in vitro and characterized resistant parasites at the genomic and phenotypic levels. We found that it took more time to develop resistance to TCMDC-143345 than to other drugs in clinical use and that there was no cross-resistance to these drugs, suggesting a new and unique mechanism. By whole-genome sequencing, we found two mutations in the gene encoding the L. donovani dynamin-1-like protein (LdoDLP1) that were fixed at the highest drug pressure. Through phylogenetic analysis, we identified LdoDLP1 as a family member of the dynamin-related proteins, a group of proteins that impacts the shapes of biological membranes by mediating fusion and fission events, with a putative role in mitochondrial fission. We found that L. donovani lines genetically engineered to harbor the two identified LdoDLP1 mutations were resistant to TCMDC-143345 and displayed altered mitochondrial properties. By homology modeling, we showed how the two LdoDLP1 mutations may influence protein structure and function. Taken together, our data reveal a clear involvement of LdoDLP1 in the adaptation/reduced susceptibility of L. donovani to TCMDC-143345. IMPORTANCE Humans and their pathogens are continuously locked in a molecular arms race during which the eventual emergence of pathogen drug resistance (DR) seems inevitable. For neglected tropical diseases (NTDs), DR is generally studied retrospectively once it has already been established in clinical settings. We previously recommended to keep one step ahead in the host-pathogen arms race and implement prospective DR studies in the R&D pipeline, a common practice for many infectious diseases but not for NTDs. Here, using Leishmania donovani, the etiological agent of visceral leishmaniasis (VL), and TCMDC-143345, a promising compound of the GSK Leishbox to treat VL, as paradigms, we experimentally selected resistance to the compound and proceeded to genomic and phenotypic characterization of DR parasites. The results gathered in the present study suggest a new DR mechanism involving the L. donovani dynamin-1-like protein (LdoDLP1) and demonstrate the practical relevance of prospective DR studies.

Auranofin (AF) is an FDA-approved antirheumatic drug with anticancer properties that acts as a thioredoxin reductase 1 (TrxR) inhibitor. The exact mechanisms through which AF targets cancer cells remain elusive. To shed light on the mode of action, this study provides an in-depth analysis on the molecular mechanisms and immunogenicity of AF-mediated cytotoxicity in the non-small cell lung cancer (NSCLC) cell line NCI–H1299 (p53 Null) and its two isogenic derivates with mutant p53 R175H or R273H accumulation.TrxR is highly expressed in a panel of 72 NSCLC patients, making it a valid druggable target in NSCLC for AF. The presence of mutant p53 overexpression was identified as an important sensitizer for AF in (isogenic) NSCLC cells as it was correlated with reduced thioredoxin (Trx) levels in vitro. Transcriptome analysis revealed dysregulation of genes involved in oxidative stress response, DNA damage, granzyme A (GZMA) signaling and ferroptosis. Although functionally AF appeared a potent inhibitor of GPX4 in all NCI–H1299 cell lines, the induction of lipid peroxidation and consequently ferroptosis was limited to the p53 R273H expressing cells. In the p53 R175H cells, AF mainly induced large-scale DNA damage and replication stress, leading to the induction of apoptotic cell death rather than ferroptosis. Importantly, all cell death types were immunogenic since the release of danger signals (ecto-calreticulin, ATP and HMGB1) and dendritic cell maturation occurred irrespective of (mutant) p53 expression. Finally, we show that AF sensitized cancer cells to caspase-independent natural killer cell-mediated killing by downregulation of several key targets of GZMA.Our data provides novel insights on AF as a potent, clinically available, off-patent cancer drug by targeting mutant p53 cancer cells through distinct cell death mechanisms (apoptosis and ferroptosis). In addition, AF improves the innate immune response at both cytostatic (natural killer cell-mediated killing) and cytotoxic concentrations (dendritic cell maturation).

Background: Visceral leishmaniasis (VL) treatment in HIV patients very often fails and is followed by high relapse and case-fatality rates. Hence, treatment efficacy assessment is imperative but based on invasive organ aspiration for parasite detection. In the search of a less-invasive alternative and because the host immune response is pivotal for treatment outcome in immunocompromised VL patients, we studied changes in the whole blood transcriptional profile of VL-HIV patients during treatment. Methods: Embedded in a clinical trial in Northwest Ethiopia, RNA-Seq was performed on whole blood samples of 28 VL-HIV patients before and after completion of a 29-day treatment regimen of AmBisome or AmBisome/miltefosine. Pathway analyses were combined with a machine learning approach to establish a clinically-useful 4-gene set. Findings: Distinct signatures of differentially expressed genes between D0 and D29 were identified for patients who failed treatment and were successfully treated. Pathway analyses in the latter highlighted a downregulation of genes associated with host cellular activity and immunity, and upregulation of antimicrobial peptide activity in phagolysosomes. No signs of disease remission nor pathway enrichment were observed in treatment failure patients. Next, we identified a 4-gene pre-post signature (PRSS33, IL10, SLFN14, HRH4) that could accurately discriminate treatment outcome at end of treatment (D29), displaying an average area-under-the-ROC-curve of 0.95 (CI: 0.75–1.00). Interpretation: A simple blood-based signature thus holds significant promise to facilitate treatment efficacy monitoring and provide an alternative test-of-cure to guide patient management in VL-HIV patients. Funding: Project funding was provided by the AfricoLeish project, supported by the European Union Seventh Framework Programme (EU FP7).

Summary: The C-terminal variants G1 and G2 of apolipoprotein L1 (APOL1) confer human resistance to the sleeping sickness parasite Trypanosoma rhodesiense, but they also increase the risk of kidney disease. APOL1 and APOL3 are death-promoting proteins that are partially associated with the endoplasmic reticulum and Golgi membranes. We report that in podocytes, either APOL1 C-terminal helix truncation (APOL1Δ) or APOL3 deletion (APOL3KO) induces similar actomyosin reorganization linked to the inhibition of phosphatidylinositol-4-phosphate [PI(4)P] synthesis by the Golgi PI(4)-kinase IIIB (PI4KB). Both APOL1 and APOL3 can form K+ channels, but only APOL3 exhibits Ca2+-dependent binding of high affinity to neuronal calcium sensor-1 (NCS-1), promoting NCS-1-PI4KB interaction and stimulating PI4KB activity. Alteration of the APOL1 C-terminal helix triggers APOL1 unfolding and increased binding to APOL3, affecting APOL3-NCS-1 interaction. Since the podocytes of G1 and G2 patients exhibit an APOL1Δ or APOL3KO-like phenotype, APOL1 C-terminal variants may induce kidney disease by preventing APOL3 from activating PI4KB, with consecutive actomyosin reorganization of podocytes. : Uzureau et al. propose a molecular explanation for the linkage between resistance to sleeping sickness and high risk of kidney disease in African individuals expressing APOL1 variants. These variants resist neutralization by T. rhodesiense SRA and kill the parasite, but they also indirectly trigger podocyte actomyosin reorganization through interaction with APOL3, inactivating the NCS-1-PI4KB complex. Keywords: APOL1, APOL3, NCS-1, PI4KB, MYH9, kidney disease, sleeping sickness, actomyosin cytoskeleton, phosphoinositide control

Abstract This article describes the motivation, origin and evolution of the student symposia series organised by the ISCB Student Council. The meeting series started thirteen years ago in Madrid and has spread to four continents. The article concludes with the highlights of the most recent edition of annual Student Council Symposium held in conjunction with the 25th Conference on Intelligent Systems for Molecular Biology and the 16th European Conference on Computational Biology, in Prague, in July 2017.

Disease relapse and therapy resistance remain key challenges in treating multiple myeloma. Underlying (epi-)mutational events can promote myelomagenesis and contribute to multi-drug and apoptosis resistance. Therefore, compounds inducing ferroptosis, a form of iron and lipid peroxidation-regulated cell death, are appealing alternative treatment strategies for multiple myeloma and other malignancies. Both ferroptosis and the epigenetic machinery are heavily influenced by oxidative stress and iron metabolism changes. Yet, only a limited number of epigenetic enzymes and modifications have been identified as ferroptosis regulators. In this study, we found that MM1 multiple myeloma cells are sensitive to ferroptosis induction and epigenetic reprogramming by RSL3, irrespective of their glucocorticoid-sensitivity status. LC-MS/MS analysis revealed the formation of non-heme iron-histone complexes and altered expression of histone modifications associated with DNA repair and cellular senescence. In line with this observation, EPIC BeadChip measurements of significant DNA methylation changes in ferroptotic myeloma cells demonstrated an enrichment of CpG probes located in genes associated with cell cycle progression and senescence, such as Nuclear Receptor Subfamily 4 Group A member 2 (NR4A2). Overall, our data show that ferroptotic cell death is associated with an epigenomic stress response that might advance the therapeutic applicability of ferroptotic compounds.

Abstract High throughput sequencing techniques are poorly adapted for in vivo studies of parasites, which require prior in vitro culturing and purification. Trypanosomatids, a group of kinetoplastid protozoans, possess a distinctive feature in their transcriptional mechanism whereby a specific Spliced Leader (SL) sequence is added to the 5′end of each mRNA by trans-splicing. This allows to discriminate Trypansomatid RNA from mammalian RNA and forms the basis of our new multiplexed protocol for high-throughput, selective RNA-sequencing called SL-seq. We provided a proof-of-concept of SL-seq in Leishmania donovani, the main causative agent of visceral leishmaniasis in humans, and successfully applied the method to sequence Leishmania mRNA directly from infected macrophages and from highly diluted mixes with human RNA. mRNA profiles obtained with SL-seq corresponded largely to those obtained from conventional poly-A tail purification methods, indicating both enumerate the same mRNA pool. However, SL-seq offers additional advantages, including lower sequencing depth requirements, fast and simple library prep and high resolution splice site detection. SL-seq is therefore ideal for fast and massive parallel sequencing of parasite transcriptomes directly from host tissues. Since SLs are also present in Nematodes, Cnidaria and primitive chordates, this method could also have high potential for transcriptomics studies in other organisms.

Susceptibility for leishmaniasis is largely dependent on host genetic and immune factors. Despite the previously described association of human leukocyte antigen (HLA) gene cluster variants as genetic susceptibility factors for leishmaniasis, little is known regarding the mechanisms that underpin these associations. To better understand this underlying functionality, we first collected all known leishmaniasis-associated HLA variants in a thorough literature review. Next, we aligned and compared the protection- and risk-associated HLA-DRB1 allele sequences. This identified several amino acid polymorphisms that distinguish protection- from risk-associated HLA-DRB1 alleles. Subsequently, T cell epitope binding predictions were carried out across these alleles to map the impact of these polymorphisms on the epitope binding repertoires. For these predictions, we used epitopes derived from entire proteomes of multiple Leishmania species. Epitopes binding to protection-associated HLA-DRB1 alleles shared common binding core motifs, mapping to the identified HLA-DRB1 amino acid polymorphisms. These results strongly suggest that HLA polymorphism, resulting in differential antigen presentation, affects the association between HLA and leishmaniasis disease development. Finally, we established a valuable open-access resource of putative epitopes. A set of 14 HLA-unrestricted strong-binding epitopes, conserved across species, was prioritized for further epitope discovery in the search for novel subunit-based vaccines.

Cytosine C5 methylation is an important epigenetic control mechanism in a wide array of eukaryotic organisms and generally carried out by proteins of the C-5 DNA methyltransferase family (DNMTs). In several protozoans, the status of this mechanism remains elusive, such as in Leishmania, the causative agent of the disease leishmaniasis in humans and a wide array of vertebrate animals. In this work, we showed that the Leishmania donovani genome contains a C-5 DNA methyltransferase (DNMT) from the DNMT6 subfamily, whose function is still unclear, and verified its expression at the RNA level. We created viable overexpressor and knock-out lines of this enzyme and characterized their genome-wide methylation patterns using whole-genome bisulfite sequencing, together with promastigote and amastigote control lines. Interestingly, despite the DNMT6 presence, we found that methylation levels were equal to or lower than 0.0003% at CpG sites, 0.0005% at CHG sites, and 0.0126% at CHH sites at the genomic scale. As none of the methylated sites were retained after manual verification, we conclude that there is no evidence for DNA methylation in this species. We demonstrated that this difference in DNA methylation between the parasite (no detectable DNA methylation) and the vertebrate host (DNA methylation) allowed enrichment of parasite vs. host DNA using methyl-CpG-binding domain columns, readily available in commercial kits. As such, we depleted methylated DNA from mixes of Leishmania promastigote and amastigote DNA with human DNA, resulting in average Leishmania:human enrichments from 62× up to 263×. These results open a promising avenue for unmethylated DNA enrichment as a pre-enrichment step before sequencing Leishmania clinical samples.

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.

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.

Leishmania donovani causes visceral leishmaniasis (VL), the second most deadly vector-borne parasitic disease. A recent epidemic in the Indian subcontinent (ISC) caused up to 80% of global VL and over 30,000 deaths per year. Resistance against antimonial drugs has probably been a contributing factor in the persistence of this epidemic. Here we use whole genome sequences from 204 clinical isolates to track the evolution and epidemiology of L. donovani from the ISC. We identify independent radiations that have emerged since a bottleneck coincident with 1960s DDT spraying campaigns. A genetically distinct population frequently resistant to antimonials has a two base-pair insertion in the aquaglyceroporin gene LdAQP1 that prevents the transport of trivalent antimonials. We find evidence of genetic exchange between ISC populations, and show that the mutation in LdAQP1 has spread by recombination. Our results reveal the complexity of L. donovani evolution in the ISC in response to drug treatment.

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.

ABSTRACT Leishmania donovani is an intracellular protozoan parasite that causes leishmaniasis, which can range from a self-healing cutaneous disease to a fatal visceral disease depending on the infecting species. Miltefosine is currently the latest and only oral antileishmanial that came out of drug discovery pipelines in the past few decades, but recent reports indicate a significant decline in its efficacy against visceral leishmaniasis (also known as kala-azar) in the Indian subcontinent. This relapse rate of up to 20% within 12 months after treatment was shown not to be related to reinfection, drug quality, drug exposure, or drug-resistant parasites. We therefore aimed to assess other phenotypes of the parasite that may affect treatment outcome and found a significant association between the number of metacyclic parasites, parasite infectivity, and patient treatment outcome in the Indian subcontinent. Together with previous studies on resistance of L. donovani against pentavalent antimonials, these data suggest that the infectivity of the parasite, or related phenotypes, might be a more determinant factor for treatment failure in visceral leishmaniasis than drug susceptibility, warranting a reassessment of our current view on treatment failure and drug resistance in leishmaniasis and beyond. IMPORTANCE The high miltefosine relapse rate poses a major challenge for the current Kala-Azar Elimination Program in the Indian subcontinent and other leishmaniasis control programs worldwide. This relapse rate could not be related to reinfection, drug-resistant parasites, or reduced treatment quality. Here we report that an increased infectivity of the parasite is associated with miltefosine relapse of visceral leishmaniasis (VL) patients. These results supplement those obtained with antimonial-resistant L. donovani where an increased infectivity was also observed. This challenges the current view of Leishmania drug susceptibility being the biggest parasitic factor that contributes to treatment failure in leishmaniasis. These selected more infectious parasites may pose an additional burden to leishmaniasis control programs, highlighting the importance of multifaceted control measures to achieve leishmaniasis elimination in the Indian subcontinent and other regions where leishmaniasis is endemic.

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.

Since 2004, the ISCB Student Council has been organizing different symposia worldwide, gathering together the community of young computational biologists. Due to the coronavirus disease 2019 (COVID-19) pandemic situation, the world scientific community was forced to cancel in-person meetings for almost two years, imposing the adoption of virtual formats instead. After the successful editions of our continental symposia in 2020 in the USA, Latin America, and Europe, we organized our flagship global event, the Student Council Symposium (SCS) 2021, trying to apply all previous lessons learned and to exploit the advantages that virtuality has to offer.

In this meeting overview, we summarise the scientific program and organisation of the 16th International Society for Computational Biology Student Council Symposium in 2020 (ISCB SCS2020). This symposium was the first virtual edition in an uninterrupted series of symposia that has been going on for 15 years, aiming to unite computational biology students and early career researchers across the globe.

Regional Student Groups (RSGs) of the International Society for Computational Biology Student Council (ISCB-SC) have been instrumental to connect computational biologists globally and to create more awareness about bioinformatics education. This article highlights the initiatives carried out by the RSGs both nationally and internationally to strengthen the present and future of the bioinformatics community. Moreover, we discuss the future directions the organization will take and the challenges to advance further in the ISCB-SC main mission: “Nurture the new generation of computational biologists”.

This editorial provides a brief overview of the 12th International Society for Computational Biology (ISCB) Student Council Symposium and the 4th European Student Council Symposium held in Florida, USA and The Hague, Netherlands, respectively. Further, the role of the ISCB Student Council in promoting education and networking in the field of computational biology is also highlighted.