Environmental Conditions May Shape the Patterns of Genomic Variations in Leishmania panamensis

Genomic diversity and genetic variation of Leishmania panamensis within its endemic range

Species belonging to the Leishmania (Viannia) subgenus are important causative agents of cutaneous and mucocutaneous leishmaniasis in Central and South America. These parasites possess several distinctive biological features that are influenced by their genetics, population structure, and genome instability. To date, several studies have revealed varying degrees of genetic diversity within Leishmania species. Particularly, in species of the L. (Viannia) subgenus, a generalized high intraspecific genetic diversity has been reported, although, conflicting conclusions have been drawn using different molecular techniques. Despite being the most common Leishmania species circulating in Panama and Colombia, few studies have analyzed clinical samples of Leishmania panamensis using whole-genome sequencing, and their restricted number of samples has limited the information they can provide to understand the population structure of L. panamensis. Here, we used next generation sequencing (NGS) to explore the genetic diversity of L. panamensis within its endemic range, analyzing data from 43 isolates of Colombian and Panamanian origin. Our results show the occurrence of three well-defined geographically correlated groups, and suggests the possible occurrence of additional phylogeographic groups. Furthermore, these results support the existence of a mixed mode of reproduction in L. panamensis, with varying frequencies of events of genetic recombination occurring primarily within subpopulations of closely related strains. This study offers important insights into the population genetics and reproduction mode of L. panamensis, paving the way to better understand their population structure and the emergence and maintenance of key eco-epidemiological traits.

Impact of Genetic Diversity and Genome Plasticity of Leishmania spp. in Treatment and the Search for Novel Chemotherapeutic Targets

Leishmaniasis is one of the major public health concerns in Latin America, Africa, Asia, and Europe. The absence of vaccines for human use and the lack of effective vector control programs make chemotherapy the main strategy to control all forms of the disease. However, the high toxicity of available drugs, limited choice of therapeutic agents, and occurrence of drug-resistant parasite strains are the main challenges related to chemotherapy. Currently, only a small number of drugs are available for leishmaniasis treatment, including pentavalent antimonials (Sb^V), amphotericin B and its formulations, miltefosine, paromomycin sulphate, and pentamidine isethionate. In addition to drug toxicity, therapeutic failure of leishmaniasis is a serious concern. The occurrence of drug-resistant parasites is one of the causes of therapeutic failure and is closely related to the diversity of parasites in this genus. Owing to the enormous plasticity of the genome, resistance can occur by altering different metabolic pathways, demonstrating that resistance mechanisms are multifactorial and extremely complex. Genetic variability and genome plasticity cause not only the available drugs to have limitations, but also make the search for new drugs challenging. Here, we examined the biological characteristics of parasites that hinder drug discovery.

Unpicking the Roles of DNA Damage Protein Kinases in Trypanosomatids

To preserve genome integrity when faced with DNA lesions, cells activate and coordinate a multitude of DNA repair pathways to ensure timely error correction or tolerance, collectively called the DNA damage response (DDR). These interconnecting damage response pathways are molecular signal relays, with protein kinases (PKs) at the pinnacle. Focused efforts in model eukaryotes have revealed intricate aspects of DNA repair PK function, including how they direct DDR pathways and how repair reactions connect to wider cellular processes, including DNA replication and transcription. The Kinetoplastidae, including many parasites like Trypanosoma spp. and Leishmania spp. (causative agents of debilitating, neglected tropical infections), exhibit peculiarities in several core biological processes, including the predominance of multigenic transcription and the streamlining or repurposing of DNA repair pathways, such as the loss of non-homologous end joining and novel operation of nucleotide excision repair (NER). Very recent studies have implicated ATR and ATM kinases in the DDR of kinetoplastid parasites, whereas DNA-dependent protein kinase (DNA-PKcs) displays uncertain conservation, questioning what functions it fulfills. The wide range of genetic manipulation approaches in these organisms presents an opportunity to investigate DNA repair kinase roles in kinetoplastids and to ask if further kinases are involved. Furthermore, the availability of kinase inhibitory compounds, targeting numerous eukaryotic PKs, could allow us to test the suitability of DNA repair PKs as novel chemotherapeutic targets. Here, we will review recent advances in the study of trypanosomatid DNA repair kinases.

Comparative analysis of the transcriptional responses of five Leishmania species to trivalent antimony

BACKGROUND

Leishmaniasis is a neglected tropical disease caused by several species of Leishmania. The resistance phenotype of these parasites depends on the characteristics of each species, which contributes to increased therapeutic failures. Understanding the mechanism used by the parasite to survive under treatment pressure in order to identify potential common and specific therapeutic targets is essential for the control of leishmaniasis. The aim of this study was to investigate the expression profiles and potential shared and specific resistance markers of the main Leishmania species of medical importance [subgenus L. (Leishmania): L. donovani, L. infantum and L. amazonensis; subgenus L. (Viannia): L. panamensis and L. braziliensis)] resistant and sensitive to trivalent stibogluconate (SbIII).

METHODS

We conducted comparative analysis of the transcriptomic profiles (only coding sequences) of lines with experimentally induced resistance to SbIII from biological replicates of five Leishmania species available in the databases of four articles based on ortholog attribution. Simultaneously, we carried out functional analysis of ontology and reconstruction of metabolic pathways of the resulting differentially expressed genes (DEGs).

RESULTS

Resistant lines for each species had differential responses in metabolic processes, compound binding, and membrane components concerning their sensitive counterpart. One hundred and thirty-nine metabolic pathways were found, with the three main pathways comprising cysteine and methionine metabolism, glycolysis, and the ribosome. Differentially expressed orthologous genes assigned to species-specific responses predominated, with 899 self-genes. No differentially expressed genes were found in common among the five species. Two common upregulated orthologous genes were found among four species (L. donovani, L. braziliensis, L. amazonensis, and L. panamensis) related to an RNA-binding protein and the NAD(P)H cytochrome-B5-oxidoreductase complex, associated with transcriptional control and de novo synthesis of linoleic acid, critical mechanisms in resistance to antimonials.

CONCLUSION

Herein, we identified potential species-specific genes related to resistance to SbIII. Therefore, we suggest that future studies consider a treatment scheme that is species-specific. Despite the limitations of our study, this is the first approach toward unraveling the pan-genus genetic mechanisms of resistance in leishmaniasis.

Genomic Diversification, Structural Plasticity, and Hybridization in Leishmania (Viannia) braziliensis

Leishmania (Viannia) braziliensis is an important Leishmania species circulating in several Central and South American countries. Among Leishmania species circulating in Brazil, Argentina and Colombia, L. braziliensis has the highest genomic variability. However, genomic variability at the whole genome level has been only studied in Brazilian and Peruvian isolates; to date, no Colombian isolates have been studied. Considering that in Colombia, L. braziliensis is a species with great clinical and therapeutic relevance, as well as the role of genetic variability in the epidemiology of leishmaniasis, we analyzed and evaluated intraspecific genomic variability of L. braziliensis from Colombian and Bolivian isolates and compared them with Brazilian isolates. Twenty-one genomes were analyzed, six from Colombian patients, one from a Bolivian patient, and 14 Brazilian isolates downloaded from public databases. The results obtained of Phylogenomic analysis showed the existence of four well-supported clades, which evidenced intraspecific variability. The whole-genome analysis revealed structural variations in the somy, mainly in the Brazilian genomes (clade 1 and clade 3), low copy number variations, and a moderate number of single-nucleotide polymorphisms (SNPs) in all genomes analyzed. Interestingly, the genomes belonging to clades 2 and 3 from Colombia and Brazil, respectively, were characterized by low heterozygosity (~90% of SNP loci were homozygous) and regions suggestive of loss of heterozygosity (LOH). Additionally, we observed the drastic whole genome loss of heterozygosity and possible hybridization events in one genome belonging to clade 4. Unique/shared SNPs between and within the four clades were identified, revealing the importance of some of them in biological processes of L. braziliensis. Our analyses demonstrate high genomic variability of L. braziliensis in different regions of South America, mainly in Colombia and suggest that this species exhibits striking genomic diversity and a capacity of genomic hybridization; additionally, this is the first study to report whole-genome sequences of Colombian L. braziliensis isolates.

In-depth quantitative proteomics uncovers specie-specific metabolic programs in Leishmania (Viannia) species

Leishmania species are responsible for a broad spectrum of diseases, denominated Leishmaniasis, affecting over 12 million people worldwide. During the last decade, there have been impressive efforts for sequencing the genome of most of the pathogenic Leishmania spp. as well as hundreds of strains, but large-scale proteomics analyses did not follow these achievements and the Leishmania proteome remained mostly uncharacterized. Here, we report a comprehensive comparative study of the proteomes of strains representing L. braziliensis, L. panamensis and L. guyanensis species. Proteins extracted by SDS-mediated lysis were processed following the multi-enzyme digestion-filter aided sample preparation (FASP) procedure and analysed by high accuracy mass spectrometry. "Total Protein Approach" and "Proteomic Ruler" were applied for absolute quantification of proteins. Principal component analysis demonstrated very high reproducibility among biological replicates and a very clear differentiation of the three species. Our dataset comprises near 7000 proteins, representing the most complete Leishmania proteome yet known, and provides a comprehensive quantitative picture of the proteomes of the three species in terms of protein concentration and copy numbers. Analysis of the abundance of proteins from the major energy metabolic processes allow us to highlight remarkably differences among the species and suggest that these parasites depend on distinct energy substrates to obtain ATP. Whereas L. braziliensis relies the more on glycolysis, L. panamensis and L. guyanensis seem to depend mainly on mitochondrial respiration. These results were confirmed by biochemical assays showing opposite profiles for glucose uptake and O2 consumption in these species. In addition, we provide quantitative data about different membrane proteins, transporters, and lipids, all of which contribute for significant species-specific differences and provide rich substrate for explore new molecules for diagnosing purposes. Data are available via ProteomeXchange with identifier PXD017696.

Intraspecific Genomic Divergence and Minor Structural Variations in Leishmania (Viannia) panamensis

Leishmania (Viannia) panamensis is one of the most important Leishmania species associated with cutaneous leishmaniasis (CL) in Latin America. Despite its wide geographic distribution and pathogenic potential in humans and animals, the genomic variability of this species is low compared with other Leishmania species circulating in the same geographical area. No studies have reported a detailed analysis of the whole genome of L. panamensis from clinical isolates using DNA high-throughput sequencing to clarify its intraspecific genomic variability or plausible divergence. Therefore, this study aimed to evaluate the intraspecific genomic variability of L. panamensis from Colombia and Panama. A total of 22 genomes were analyzed, 19 from Colombian patients with CL and three genomes from Panama obtained from public databases. The phylogenomic analysis revealed the potential existence of three well-supported clades as evidence of intraspecific divergence. Additionally, the whole-genome analysis showed low structural variations in terms of ploidy, copy number variations, and single-nucleotide polymorphisms (SNPs). SNPs shared among all clades were identified, revealing their importance in different biological processes of L. panamensis. The findings not only expand our knowledge of intraspecific genomic variability of one of the most important Leishmania species in South America but also highlights the possible existence of different clades/lineages/subpopulations across a geographic scale.