Hypoxia-Inducible Factor-1 Alpha Stabilization in Human Macrophages during Leishmania major Infection Is Impaired by Parasite Virulence.

Hypoxia-inducible factor-1 alpha (HIF-1α) is one of the master regulators of immune and metabolic cellular functions. HIF-1α, a transcriptional factor whose activity is closely related to oxygen levels, is a target for understanding infectious disease control. Several studies have demonstrated that HIF-1α plays an important role during the infectious process, while its role in relation to parasite virulence has not been addressed. In this work, we studied the expression levels of HIF-1α and related angiogenic vascular endothelial growth factor A (VEGF-A) in human macrophages infected with promastigotes of hypo- or hyper-virulent Leishmania major human isolates. L. major parasites readily subverted host macrophage functions for their survival and induced local oxygen consumption at the site of infection. In contrast to hypo-virulent parasites that induce high HIF-1α expression levels, hyper-virulent L. major reduced HIF-1α expression in macrophages under normoxic or hypoxic conditions, and consequently impeded the expression of VEGF-A mRNA. HIF-1α may play a key role during control of disease chronicity, severity, or outcome.

Healed Lesions of Human Cutaneous Leishmaniasis Caused By Leishmania major Do Not Shelter Persistent Residual Parasites.

In human cutaneous leishmaniasis (HCL) caused by Leishmania (L.) major, the cutaneous lesions heal spontaneously and induce a Th1-type immunity that confers solid protection against reinfection. The same holds true for the experimental leishmaniasis induced by L. major in C57BL/6 mice where residual parasites persist after spontaneous clinical cure and induce sustainable memory immune responses and resistance to reinfection. Whether residual parasites also persist in scars of cured HCL caused by L. major is still unknown. Cutaneous scars from 53 volunteers with healed HCL caused by L. major were biopsied and the tissue sample homogenates were analyzed for residual parasites by four methods: i) microscope detection of amastigotes, ii) parasite culture by inoculation on biphasic medium, iii) inoculation of tissue exctracts to the footpad of BALB/c mice, an inbred strain highly susceptible to L. major, and iv) amplification of parasite kDNA by a highly sensitive real-time PCR (RT-PCR). Our results show that the scars of healed lesions of HCL caused by L. major do not contain detectable residual parasites, suggesting that this form likely induces a sterile cure at least within the scars. This feature contrasts with other Leishmania species causing chronic, diffuse, or recidivating forms of leishmaniasis where parasites do persist in healed lesions. The possibility that alternative mechanisms to parasite persistence are needed to boost and maintain long-term immunity to L. major, should be taken into consideration in vaccine development against L. major infection.

Molecular Epidemiology of SARS-CoV-2 in Tunisia (North Africa) through Several Successive Waves of COVID-19.

Documenting the circulation dynamics of SARS-CoV-2 variants in different regions of the world is crucial for monitoring virus transmission worldwide and contributing to global efforts towards combating the pandemic. Tunisia has experienced several waves of COVID-19 with a significant number of infections and deaths. The present study provides genetic information on the different lineages of SARS-CoV-2 that circulated in Tunisia over 17 months. Lineages were assigned for 1359 samples using whole-genome sequencing, partial S gene sequencing and variant-specific real-time RT-PCR tests. Forty-eight different lineages of SARS-CoV-2 were identified, including variants of concern (VOCs), variants of interest (VOIs) and variants under monitoring (VUMs), particularly Alpha, Beta, Delta, A.27, Zeta and Eta. The first wave, limited to imported and import-related cases, was characterized by a small number of positive samples and lineages. During the second wave, a large number of lineages were detected; the third wave was marked by the predominance of the Alpha VOC, and the fourth wave was characterized by the predominance of the Delta VOC. This study adds new genomic data to the global context of COVID-19, particularly from the North African region, and highlights the importance of the timely molecular characterization of circulating strains.

Designing and running an advanced Bioinformatics and genome analyses course in Tunisia.

Genome data, with underlying new knowledge, are accumulating at exponential rate thanks to ever-improving sequencing technologies and the parallel development of dedicated efficient Bioinformatics methods and tools. Advanced Education in Bioinformatics and Genome Analyses is to a large extent not accessible to students in developing countries where endeavors to set up Bioinformatics courses concern most often only basic levels. Here, we report a pioneering pilot experience concerning the design and implementation, from scratch, of a three-months advanced and extensive course in Bioinformatics and Genome Analyses in the Institut Pasteur de Tunis. Most significantly the outcome of the course was upgrading the participants' skills in Bioinformatics and Genome Analyses to recognized international standards. Here we detail the different steps involved in the implementation of this course as well as the topics covered in the program. The description of this pilot experience might be helpful for the implementation of other similar educational projects, notably in developing countries, aiming to go beyond basics and providing young researchers with high-level skills.

Leishmania Genome Dynamics during Environmental Adaptation Reveal Strain-Specific Differences in Gene Copy Number Variation, Karyotype Instability, and Telomeric Amplification.

Protozoan parasites of the genus Leishmania adapt to environmental change through chromosome and gene copy number variations. Only little is known about external or intrinsic factors that govern Leishmania genomic adaptation. Here, by conducting longitudinal genome analyses of 10 new Leishmania clinical isolates, we uncovered important differences in gene copy number among genetically highly related strains and revealed gain and loss of gene copies as potential drivers of long-term environmental adaptation in the field. In contrast, chromosome rather than gene amplification was associated with short-term environmental adaptation to in vitro culture. Karyotypic solutions were highly reproducible but unique for a given strain, suggesting that chromosome amplification is under positive selection and dependent on species- and strain-specific intrinsic factors. We revealed a progressive increase in read depth towards the chromosome ends for various Leishmania isolates, which may represent a nonclassical mechanism of telomere maintenance that can preserve integrity of chromosome ends during selection for fast in vitro growth. Together our data draw a complex picture of Leishmania genomic adaptation in the field and in culture, which is driven by a combination of intrinsic genetic factors that generate strain-specific phenotypic variations, which are under environmental selection and allow for fitness gain.IMPORTANCE Protozoan parasites of the genus Leishmania cause severe human and veterinary diseases worldwide, termed leishmaniases. A hallmark of Leishmania biology is its capacity to adapt to a variety of unpredictable fluctuations inside its human host, notably pharmacological interventions, thus, causing drug resistance. Here we investigated mechanisms of environmental adaptation using a comparative genomics approach by sequencing 10 new clinical isolates of the L. donovani, L. major, and L. tropica complexes that were sampled across eight distinct geographical regions. Our data provide new evidence that parasites adapt to environmental change in the field and in culture through a combination of chromosome and gene amplification that likely causes phenotypic variation and drives parasite fitness gains in response to environmental constraints. This novel form of gene expression regulation through genomic change compensates for the absence of classical transcriptional control in these early-branching eukaryotes and opens new venues for biomarker discovery.

Role of Human Macrophage Polarization in Inflammation during Infectious Diseases.

Experimental models have often been at the origin of immunological paradigms such as the M1/M2 dichotomy following macrophage polarization. However, this clear dichotomy in animal models is not as obvious in humans, and the separating line between M1-like and M2-like macrophages is rather represented by a continuum, where boundaries are still unclear. Indeed, human infectious diseases, are characterized by either a back and forth or often a mixed profile between the pro-inflammatory microenvironment (dominated by interleukin (IL)-1ß, IL-6, IL-12, IL-23 and Tumor Necrosis Factor (TNF)-α cytokines) and tissue injury driven by classically activated macrophages (M1-like) and wound healing driven by alternatively activated macrophages (M2-like) in an anti-inflammatory environment (dominated by IL-10, Transforming growth factor (TGF)-ß, chemokine ligand (CCL)1, CCL2, CCL17, CCL18, and CCL22). This review brews the complexity of the situation during infectious diseases by stressing on this continuum between M1-like and M2-like extremes. We first discuss the basic biology of macrophage polarization, function, and role in the inflammatory process and its resolution. Secondly, we discuss the relevance of the macrophage polarization continuum during infectious and neglected diseases, and the possibility to interfere with such activation states as a promising therapeutic strategy in the treatment of such diseases.

Designing a course model for distance-based online bioinformatics training in Africa: The H3ABioNet experience.

Africa is not unique in its need for basic bioinformatics training for individuals from a diverse range of academic backgrounds. However, particular logistical challenges in Africa, most notably access to bioinformatics expertise and internet stability, must be addressed in order to meet this need on the continent. H3ABioNet (www.h3abionet.org), the Pan African Bioinformatics Network for H3Africa, has therefore developed an innovative, free-of-charge "Introduction to Bioinformatics" course, taking these challenges into account as part of its educational efforts to provide on-site training and develop local expertise inside its network. A multiple-delivery-mode learning model was selected for this 3-month course in order to increase access to (mostly) African, expert bioinformatics trainers. The content of the course was developed to include a range of fundamental bioinformatics topics at the introductory level. For the first iteration of the course (2016), classrooms with a total of 364 enrolled participants were hosted at 20 institutions across 10 African countries. To ensure that classroom success did not depend on stable internet, trainers pre-recorded their lectures, and classrooms downloaded and watched these locally during biweekly contact sessions. The trainers were available via video conferencing to take questions during contact sessions, as well as via online "question and discussion" forums outside of contact session time. This learning model, developed for a resource-limited setting, could easily be adapted to other settings.

Comparative genomics of Tunisian Leishmania major isolates causing human cutaneous leishmaniasis with contrasting clinical severity.

Zoonotic cutaneous leishmaniasis caused by Leishmania (L.) major parasites affects urban and suburban areas in the center and south of Tunisia where the disease is endemo-epidemic. Several cases were reported in human patients for which infection due to L. major induced lesions with a broad range of severity. However, very little is known about the mechanisms underlying this diversity. Our hypothesis is that parasite genomic variability could, in addition to the host immunological background, contribute to the intra-species clinical variability observed in patients and explain the lesion size differences observed in the experimental model. Based on several epidemiological, in vivo and in vitro experiments, we focused on two clinical isolates showing contrasted severity in patients and BALB/c experimental mice model. We used DNA-seq as a high-throughput technology to facilitate the identification of genetic variants with discriminating potential between both isolates. Our results demonstrate that various levels of heterogeneity could be found between both L. major isolates in terms of chromosome or gene copy number variation (CNV), and that the intra-species divergence could surprisingly be related to single nucleotide polymorphisms (SNPs) and Insertion/Deletion (InDels) events. Interestingly, we particularly focused here on genes affected by both types of variants and correlated them with the observed gene CNV. Whether these differences are sufficient to explain the severity in patients is obviously still open to debate, but we do believe that additional layers of -omic information is needed to complement the genomic screen in order to draw a more complete map of severity determinants.

Treatment with synthetic lipophilic tyrosyl ester controls Leishmania major infection by reducing parasite load in BALB/c mice.

Synthesized lipophilic tyrosyl ester derivatives with increasing lipophilicity were effective against Leishmania (L.) major and Leishmania infantum species in vitro. These findings prompted us to test in vivo leishmanicidal properties of these molecules and their potential effect on the modulation of immune responses. The experimental BALB/c model of cutaneous leishmaniasis was used in this study. Mice were infected with L. major parasites and treated with three in vitro active tyrosyl esters derivatives. Among these tested tyrosylcaprate (TyC) compounds, only TyC10 exhibited an in vivo anti-leishmanial activity, when injected sub-cutaneously (s.c.). TyC10 treatment of L. major-infected BALB/c mice resulted in a decrease of lesion development and parasite load. TyC10 s.c. treatment of non-infected mice induced an imbalance in interferon γ/interleukin 4 (IFN-γ/IL-4) ratio cytokines towards a Th1 response. Our results indicate that TyC10 s.c. treatment improves lesions' healing and parasite clearance and may act on the cytokine balance towards a Th1 protective response by decreasing IL-4 and increasing IFN-γ transcripts. TyC10 is worthy of further investigation to uncover its mechanism of action that could lead to consider this molecule as a potential drug candidate.

Genetic micro-heterogeneity of Leishmania major in emerging foci of zoonotic cutaneous leishmaniasis in Tunisia.

Tunisia is endemic for zoonotic cutaneous leishmaniasis (ZCL), a parasitic disease caused by Leishmania (L.) major. ZCL displays a wide clinical polymorphism, with severe forms present more frequently in emerging foci where naive populations are dominant. In this study, we applied the multi-locus microsatellite typing (MLMT) using ten highly informative and discriminative markers to investigate the genetic structure of 35 Tunisian Leishmania (L.) major isolates collected from patients living in five different foci of Central Tunisia (two old and three emerging foci). Phylogenetic reconstructions based on genetic distances showed that nine of the ten tested loci were homogeneous in all isolates with homozygous alleles, whereas one locus (71AT) had a 58/64-bp bi-allelic profile with an allele linked to emerging foci. Promastigote-stage parasites with the 58-bp allele tend to be more resistant to in vitro complement lysis. These results, which stress the geographical dependence of the genetic micro-heterogeneity, may improve our understanding of the ZCL epidemiology and clinical outcome.

MicroRNA expression profile in human macrophages in response to Leishmania major infection.

BACKGROUND: Leishmania (L.) are intracellular protozoan parasites able to survive and replicate in the hostile phagolysosomal environment of infected macrophages. They cause leishmaniasis, a heterogeneous group of worldwide-distributed affections, representing a paradigm of neglected diseases that are mainly embedded in impoverished populations. To establish successful infection and ensure their own survival, Leishmania have developed sophisticated strategies to subvert the host macrophage responses. Despite a wealth of gained crucial information, these strategies still remain poorly understood. MicroRNAs (miRNAs), an evolutionarily conserved class of endogenous 22-nucleotide non-coding RNAs, are described to participate in the regulation of almost every cellular process investigated so far. They regulate the expression of target genes both at the levels of mRNA stability and translation; changes in their expression have a profound effect on their target transcripts. METHODOLOGY/PRINCIPAL FINDINGS: We report in this study a comprehensive analysis of miRNA expression profiles in L. major-infected human primary macrophages of three healthy donors assessed at different time-points post-infection (three to 24 h). We show that expression of 64 out of 365 analyzed miRNAs was consistently deregulated upon infection with the same trends in all donors. Among these, several are known to be induced by TLR-dependent responses. GO enrichment analysis of experimentally validated miRNA-targeted genes revealed that several pathways and molecular functions were disturbed upon parasite infection. Finally, following parasite infection, miR-210 abundance was enhanced in HIF-1α-dependent manner, though it did not contribute to inhibiting anti-apoptotic pathways through pro-apoptotic caspase-3 regulation. CONCLUSIONS/SIGNIFICANCE: Our data suggest that alteration in miRNA levels likely plays an important role in regulating macrophage functions following L. major infection. These results could contribute to better understanding of the dynamics of gene expression in host cells during leishmaniasis.

Methodology optimizing SAGE library tag-to-gene mapping: application to Leishmania.

BACKGROUND: Leishmaniasis are widespread parasitic-diseases with an urgent need for more active and less toxic drugs and for effective vaccines. Understanding the biology of the parasite especially in the context of host parasite interaction is a crucial step towards such improvements in therapy and control. Several experimental approaches including SAGE (Serial analysis of gene expression) have been developed in order to investigate the parasite transcriptome organisation and plasticity. Usual SAGE tag-to-gene mapping techniques are inadequate because almost all tags are normally located in the 3'-UTR outside the CDS, whereas most information available for Leishmania transcripts is restricted to the CDS predictions. The aim of this work is to optimize a SAGE libraries tag-to-gene mapping technique and to show how this development improves the understanding of Leishmania transcriptome. FINDINGS: The in silico method implemented herein was based on mapping the tags to Leishmania genome using BLAST then mapping the tags to their gene using a data-driven probability distribution. This optimized tag-to-gene mappings improved the knowledge of Leishmania genome structure and transcription. It allowed analyzing the expression of a maximal number of Leishmania genes, the delimitation of the 3' UTR of 478 genes and the identification of biological processes that are differentially modulated during the promastigote to amastigote differentiation. CONCLUSION: The developed method optimizes the assignment of SAGE tags in trypanosomatidae genomes as well as in any genome having polycistronic transcription and small intergenic regions.

EuPathDomains: the divergent domain database for eukaryotic pathogens.

Eukaryotic pathogens (e.g. Plasmodium, Leishmania, Trypanosomes, etc.) are a major source of morbidity and mortality worldwide. In Africa, one of the most impacted continents, they cause millions of deaths and constitute an immense economic burden. While the genome sequence of several of these organisms is now available, the biological functions of more than half of their proteins are still unknown. This is a serious issue for bringing to the foreground the expected new therapeutic targets. In this context, the identification of protein domains is a key step to improve the functional annotation of the proteins. However, several domains are missed in eukaryotic pathogens because of the high phylogenetic distance of these organisms from the classical eukaryote models. We recently proposed a method, co-occurrence domain detection (CODD), that improves the sensitivity of Pfam domain detection by exploiting the tendency of domains to appear preferentially with a few other favorite domains in a protein. In this paper, we present EuPathDomains (http://www.atgc-montpellier.fr/EuPathDomains/), an extended database of protein domains belonging to ten major eukaryotic human pathogens. EuPathDomains gathers known and new domains detected by CODD, along with the associated confidence measurements and the GO annotations that can be deduced from the new domains. This database significantly extends the Pfam domain coverage of all selected genomes, by proposing new occurrences of domains as well as new domain families that have never been reported before. For example, with a false discovery rate lower than 20%, EuPathDomains increases the number of detected domains by 13% in Toxoplasma gondii genome and up to 28% in Cryptospordium parvum, and the total number of domain families by 10% in Plasmodium falciparum and up to 16% in C. parvum genome. The database can be queried by protein names, domain identifiers, Pfam or Interpro identifiers, or organisms, and should become a valuable resource to decipher the protein functions of eukaryotic pathogens.

Application of Multi-SOM clustering approach to macrophage gene expression analysis.

The production of increasingly reliable and accessible gene expression data has stimulated the development of computational tools to interpret such data and to organize them efficiently. The clustering techniques are largely recognized as useful exploratory tools for gene expression data analysis. Genes that show similar expression patterns over a wide range of experimental conditions can be clustered together. This relies on the hypothesis that genes that belong to the same cluster are coregulated and involved in related functions. Nevertheless, clustering algorithms still show limits, particularly for the estimation of the number of clusters and the interpretation of hierarchical dendrogram, which may significantly influence the outputs of the analysis process. We propose here a multi level SOM based clustering algorithm named Multi-SOM. Through the use of clustering validity indices, Multi-SOM overcomes the problem of the estimation of clusters number. To test the validity of the proposed clustering algorithm, we first tested it on supervised training data sets. Results were evaluated by computing the number of misclassified samples. We have then used Multi-SOM for the analysis of macrophage gene expression data generated in vitro from the same individual blood infected with 5 different pathogens. This analysis led to the identification of sets of tightly coregulated genes across different pathogens. Gene Ontology tools were then used to estimate the biological significance of the clustering, which showed that the obtained clusters are coherent and biologically significant.

Simultaneous gene expression profiling in human macrophages infected with Leishmania major parasites using SAGE.

BACKGROUND: Leishmania (L) are intracellular protozoan parasites that are able to survive and replicate within the harsh and potentially hostile phagolysosomal environment of mammalian mononuclear phagocytes. A complex interplay then takes place between the macrophage (MPhi) striving to eliminate the pathogen and the parasite struggling for its own survival. To investigate this host-parasite conflict at the transcriptional level, in the context of monocyte-derived human MPhis (MDM) infection by L. major metacyclic promastigotes, the quantitative technique of serial analysis of gene expression (SAGE) was used. RESULTS: After extracting mRNA from resting human MPhis, Leishmania-infected human MPhis and L. major parasites, three SAGE libraries were constructed and sequenced generating up to 28,173; 57,514 and 33,906 tags respectively (corresponding to 12,946; 23,442 and 9,530 unique tags). Using computational data analysis and direct comparison to 357,888 publicly available experimental human tags, the parasite and the host cell transcriptomes were then simultaneously characterized from the mixed cellular extract, confidently discriminating host from parasite transcripts. This procedure led us to reliably assign 3,814 tags to MPhis' and 3,666 tags to L. major parasites transcripts. We focused on these, showing significant changes in their expression that are likely to be relevant to the pathogenesis of parasite infection: (i) human MPhis genes, belonging to key immune response proteins (e.g., IFNgamma pathway, S100 and chemokine families) and (ii) a group of Leishmania genes showing a preferential expression at the parasite's intra-cellular developing stage. CONCLUSION: Dual SAGE transcriptome analysis provided a useful, powerful and accurate approach to discriminating genes of human or parasitic origin in Leishmania-infected human MPhis. The findings presented in this work suggest that the Leishmania parasite modulates key transcripts in human MPhis that may be beneficial for its establishment and survival. Furthermore, these results provide an overview of gene expression at two developmental stages of the parasite, namely metacyclic promastigotes and intracellular amastigotes and indicate a broad difference between their transcriptomic profiles. Finally, our reported set of expressed genes will be useful in future rounds of data mining and gene annotation.