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

Diversity and dissemination of viruses in pathogenic protozoa

Viruses are the most abundant biological entities on Earth and play a significant role in the evolution of many organisms and ecosystems. In pathogenic protozoa, the presence of viruses has been linked to an increased risk of treatment failure and severe clinical outcome. Here, we studied the molecular epidemiology of the zoonotic disease cutaneous leishmaniasis in Peru and Bolivia through a joint evolutionary analysis of Leishmania braziliensis and their dsRNA Leishmania virus 1. We show that parasite populations circulate in tropical rainforests and are associated with single viral lineages that appear in low prevalence. In contrast, groups of hybrid parasites are geographically and ecologically more dispersed and associated with an increased prevalence, diversity and spread of viruses. Our results suggest that parasite gene flow and hybridization increased the frequency of parasite-virus symbioses, a process that may change the epidemiology of leishmaniasis in the region.

Insights into the Genetic Diversity of Leishmania (Viannia) panamensis in Panama, Inferred via Multilocus Sequence Typing (MLST)

Leishmaniasis is a disease caused by parasites of the genus Leishmania and transmitted by sand fly vectors. Tegumentary leishmaniasis is the most prevalent clinical outcome in Latin America, afflicting people from 18 countries. In Panama, the annual incidence rate of leishmaniasis is as high as 3000 cases, representing a major public health problem. In endemic regions, L. panamensis is responsible for almost eighty percent of human cases that present different clinical outcomes. These differences in disease outcomes could be the result of the local interplay between L. panamensis variants and human hosts with different genetic backgrounds. The genetic diversity of L. panamensis in Panama has only been partially explored, and the variability reported for this species is based on few studies restricted to small populations and/or with poor resolutive markers at low taxonomic levels. Accordingly, in this study, we explored the genetic diversity of sixty-nine L. panamensis isolates from different endemic regions of Panama, using an MLST approach based on four housekeeping genes (Aconitase, ALAT, GPI and HSP70). Two to seven haplotypes per locus were identified, and regional differences in the genetic diversity of L. panamensis were observed. A genotype analysis evidenced the circulation of thirteen L. panamensis genotypes, a fact that might have important implications for the local control of the disease.

Parasite hybridization promotes spreading of endosymbiotic viruses

Viruses are the most abundant biological entities on Earth and play a significant role in the evolution of many organisms and ecosystems. In pathogenic protozoa, the presence of endosymbiotic viruses has been linked to an increased risk of treatment failure and severe clinical outcome. Here, we studied the molecular epidemiology of the zoonotic disease cutaneous leishmaniasis in Peru and Bolivia through a joint evolutionary analysis of Leishmania braziliensis parasites and their endosymbiotic Leishmania RNA virus. We show that parasite populations circulate in isolated pockets of suitable habitat and are associated with single viral lineages that appear in low prevalence. In contrast, groups of hybrid parasites were geographically and ecologically dispersed, and commonly infected from a pool of genetically diverse viruses. Our results suggest that parasite hybridization, likely due to increased human migration and ecological perturbations, increased the frequency of endosymbiotic interactions known to play a key role in disease severity.

Selective whole-genome amplification reveals population genetics of Leishmania braziliensis directly from patient skin biopsies

In Brazil, Leishmania braziliensis is the main causative agent of the neglected tropical disease, cutaneous leishmaniasis (CL). CL presents on a spectrum of disease severity with a high rate of treatment failure. Yet the parasite factors that contribute to disease presentation and treatment outcome are not well understood, in part because successfully isolating and culturing parasites from patient lesions remains a major technical challenge. Here we describe the development of selective whole genome amplification (SWGA) for Leishmania and show that this method enables culture-independent analysis of parasite genomes obtained directly from primary patient skin samples, allowing us to circumvent artifacts associated with adaptation to culture. We show that SWGA can be applied to multiple Leishmania species residing in different host species, suggesting that this method is broadly useful in both experimental infection models and clinical studies. SWGA carried out directly on skin biopsies collected from patients in Corte de Pedra, Bahia, Brazil, showed extensive genomic diversity. Finally, as a proof-of-concept, we demonstrated that SWGA data can be integrated with published whole genome data from cultured parasite isolates to identify variants unique to specific geographic regions in Brazil where treatment failure rates are known to be high. SWGA provides a relatively simple method to generate Leishmania genomes directly from patient samples, unlocking the potential to link parasite genetics with host clinical phenotypes.

Draft genomes of Blastocystis subtypes from human samples of Colombia

BACKGROUND

Blastocystis is one of the most common eukaryotic microorganisms colonizing the intestines of both humans and animals, but the conditions under which it may be a pathogen are unclear.

METHODS

To study the genomic characteristics of circulating subtypes (ST) in Colombia, we established nine xenic cultures from Blastocystis isolated from human fecal samples, we identified 10 different subtypes, since one sample had a mixed infection. Thus, the genomes of the subtypes ST1 (n = 3), ST2 (n = 1), ST3 (n = 2), ST6 (n = 1), ST7 (n = 1), and ST8 (n = 2) were sequenced using Illumina and Oxford Nanopore Technologies (ONT).

RESULTS

Analyses of these draft nuclear genomes indicated remarkable diversity in terms of genome size and guanine-cytosine (GC) content among the compared STs. Illumina sequencing-only draft genomes contained 824 to 2077 scaffolds, with total genome size ranging from 12 to 13.2 Mb and N₅₀ values ranging from 10,585 to 29,404 base pairs (bp). The genome of one ST1 isolate was sequenced using ONT. This assembly was more contiguous, with a size of 20 million base pairs (Mb) spread over 116 scaffolds, and an N₅₀ of 248,997 bp.

CONCLUSION

This work represents one of the few large-scale comparative genomic analyses of Blastocystis isolates, providing an additional glimpse into its genomic diversity.

Impact of collaborative actions of Leishmania (Viannia) braziliensis subpopulations on the infection profile

This study focuses on the role of the population structure of Leishmania spp. on the adaptive capacity of the parasite. Herein, we investigate the contribution of subpopulations of the L. (V.) braziliensis Thor strain (Thor03, Thor10 and Thor22) in the profile of murine macrophages infection. Infection assays were performed with binary combinations of these subpopulations at stationary phases. The initial interaction time showed major effects on the combination assays, as demonstrated by the significant increase in the infection rate at 5 h. Based on the endocytic index (EI), Thor10 (EI = 563.6) and Thor03 (EI = 497) showed a higher infection load compared to Thor22 (EI = 227.3). However, the EI decreased in Thor03 after 48 h (EI = 447) and 72 h (EI = 388.3) of infection, and showed changes in the infection level in all Thor10/Thor22 combinations. Assays with CellTrace CFSE-labelled Thor22 promastigotes indicated an increase (~1.5 fold) in infection by this subpopulation in the presence of Thor10 when compared to the infection profile of Thor03/Thor22 combinations in the same proportions. In addition, the potential of these subpopulations, alone or in binary combinations, to modulate the expression of cytokines and nitric oxide (NO) in vitro was investigated. Lower NO and tumour necrosis factor-α production levels were observed for all Thor10/Thor22 combinations at 24 h compared to these subpopulations alone. In contrast, Thor03/Thor22 combination assays increased IL-10 production at this time. Collectively, these results provide in vitro evidence on the potential of L. (V.) braziliensis population structure to play a relevant role in a host infection by this parasite.

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.

Genome plasticity driven by aneuploidy and loss of heterozygosity in Trypanosoma cruzi

Trypanosoma cruzi the causative agent of Chagas disease shows a marked genetic diversity and divided into at least six Discrete Typing Units (DTUs). High intra genetic variability has been observed in the TcI DTU, the most widely distributed DTU, where patterns of genomic diversity can provide information on ecological and evolutionary processes driving parasite population structure and genome organization. Chromosomal aneuploidies and rearrangements across multigene families represent an evidence of T. cruzi genome plasticity. We explored genomic diversity among 18 Colombian T. cruzi I clones and 15 T. cruzi I South American strains. Our results confirm high genomic variability, heterozygosity and presence of a clade compatible with the TcI_dom genotype, described for strains from humans in Colombia and Venezuela. TcI showed high structural plasticity across the geographical region studied. Differential events of whole and segmental aneuploidy (SA) along chromosomes even between clones from the same strain were found and corroborated by the depth and allelic frequency. We detected loss of heterozygosity (LOH) events in different chromosomes, however, the size and location of segments under LOH varied between clones. Genes adjacent to breakpoints were evaluated, and retrotransposon hot spot genes flanked the beginning of segmental aneuploidies. Our results suggest that T. cruzi genomes, like those of Leishmania, may have a highly unstable structure and there is now an urgent need to design experiments to explore any potential adaptive role for the plasticity observed.

Novel Insights Into Leishmania (Viannia) braziliensis In Vitro Fitness Guided by Temperature Changes Along With Its Subtilisins and Oligopeptidase B

Proteases are virulence factors with a recognized impact on the Leishmania spp. life cycle. This study considers a set of analyses measuring phenotypic factors of L. (V.) braziliensis clinical isolates as promastigotes growth curves, murine peritoneal macrophages infection, inflammatory mediators production, and serine proteases gene expression (subtilisin 13: S13, subtilisin 28: S28, oligopeptidase B: OPB) assessing these isolates’ fitness on in vitro conditions. Parasites had different behavior during the early growth phase from day zero to day three, and all isolates reached the stationary growth phase between days four and seven. Macrophages infection showed two tendencies, one of decreased infection rate and number of parasites per macrophage (Infection Index <1000) and another with a constant infection index (≥1400). TNF-α (≥10 pg/mL) detected in infections by 75% of isolates, IL-6 (≥80 pg/mL) by 30% of isolates and low levels of NO (≥0.01µM) in almost all infections. Gene expression showed higher values of S13 (≥2RQ) in the intracellular amastigotes of all the isolates evaluated. On the contrary, S28 expression was low (≤1RQ) in all isolates. OPB expression was different between promastigotes and intracellular amastigotes, being significantly higher (≥2RQ) in the latter form of 58% of the isolates. Predictive structural assays of S13 and OPB were performed to explore temperature influence on gene expression and the encoded proteases. Gene expression data is discussed based on in silico predictions of regulatory regions that show plasticity in the linearity index of secondary structures of S13 and OPB 3’-untranslated regions of mRNA, dependent on temperature changes. While hairpin structures suggest an active region of mRNA for both genes above 26°C, pseudoknot structure found in S13 is an indication of a particular profile of this gene at mammalian host temperatures (37°C). Furthermore, the predicted 3D structures are in accordance with the influence of these temperatures on the catalytic site stability of both enzymes, favoring their action over peptide substrates. Data gathered here suggest that L. (V.) braziliensis serine proteases can be influenced by the temperature conditions affecting parasite fitness throughout its life cycle.

Sex under pressure: stress facilitates Leishmania in vitro hybridization

The selection of Leishmania hybrids in axenic culture was considered rare until recently, when Louradour and Ferreira et al., demonstrated that induced DNA damage facilitates genetic exchange, resulting in full genome tetraploid progenies in vitro. Meiosis-related gene homologues HAP2, GEX1, and RAD51 were found to be involved, opening new avenues for functional genomic studies.

Filling the gaps in Leishmania naiffi and Leishmania guyanensis genome plasticity

Insufficient and irregular data reports on Leishmaniasis, issuing from the developing world, have left much to be desired in terms of understanding the molecular signatures producing distinct infectious phenotypes of the disease. Herein, we report on the complete genome sequencing of Leishmania naiffi and Leishmania guyanensis, sampled from patients in regions of Colombia and Venezuela. In this study, the isolates of cutaneous lesions from both species presented limited structural variation at the chromosomal level, low gene copy number variation, and high genetic heterogeneity. We compared these sequences to the reference genomes hitherto related from Brazil and French Guyana. Although of the same species, we note a consequential genomic disparity between the Venezuelan and French Guyanese isolates of L. guyanensis. Although less significant on the global schema of cutaneous and mucosal disease, such genomic studies of L. naiffi and L. guyanensis substantiate the gaps in understanding of the molecular architecture and multivariate clinical pictures of Leishmaniasis, on an international scale.

Development of an Amplicon-Based Next-Generation Sequencing Protocol to Identify Leishmania Species and Other Trypanosomatids in Leishmaniasis Endemic Areas

Trypanosomatid infections are an important public health threat affecting many low-income countries across the tropics, particularly in the Americas. Trypanosomatids can infect many vertebrate, invertebrate, and plant species and play an important role as human pathogens. Among these clinically relevant pathogens are species from the genera Leishmania and Trypanosoma. Mixed trypanosomatid infections remain a largely unexplored phenomenon. Herein, we describe the application of an amplicon-based next-generation sequencing (NGS) assay to detect and identify trypanosomatid species in mammalian reservoirs, human patients, and sand fly vectors throughout regions of Leishmania endemicity. Sixty-five samples from different departments of Colombia, including two samples from Venezuela, were analyzed: 49 samples from cutaneous leishmaniasis (CL) patients, 8 from sand flies, 2 from domestic reservoirs (Canis familiaris), and 6 from wild reservoirs (Phyllostomus hastatus). DNA from each sample served to identify the presence of trypanosomatids through conventional PCR using heat shock protein 70 (HSP70) gene as the target. PCR products underwent sequencing by Sanger sequencing and NGS, and trypanosomatid species were identified by using BLASTn against a reference database built from trypanosomatid-derived HSP70 sequences. The alpha and beta diversity indexes of amplicon sequence variants were calculated for each group. The results revealed the presence of mixed infections with more than two Leishmania species in 34% of CL samples analyzed. Trypanosoma cruzi was identified in samples from wild reservoirs, as well as in sand fly vectors. Coinfection events with three different Leishmania species were identified in domestic reservoirs. These findings depose the traditional paradigm of leishmaniasis as being a single-species-driven infection and redraw the choreography of host-pathogen interaction in the context of multiparasitism. Further research is needed to decipher how coinfections may influence disease progression. This knowledge is key to developing an integrated approach for diagnosis and treatment.

IMPORTANCE Traditionally, there has been a frequent, yet incorrect assumption that phlebotomine vectors, animal reservoirs, and human hosts are susceptible to Leishmania infection by a single parasite species. However, current evidence supports that these new vector-parasite-reservoir associations lend vectors and reservoirs greater permissiveness to certain Leishmania species, thus promoting the appearance of coinfection events, particularly in disease-endemic regions. The present study describes the application of an amplicon-based next-generation sequencing (NGS) assay to detect and identify trypanosomatid species in mammalian reservoirs, human patients, and sand fly vectors from regions of endemicity for leishmaniasis. This changes our understanding of the clinical course of leishmaniasis in areas of endemicity.

Variation in Leishmania chemokine suppression driven by diversification of the GP63 virulence factor

Leishmaniasis is a neglected tropical disease with diverse outcomes ranging from self-healing lesions, to progressive non-healing lesions, to metastatic spread and destruction of mucous membranes. Although resolution of cutaneous leishmaniasis is a classic example of type-1 immunity leading to self-healing lesions, an excess of type-1 related inflammation can contribute to immunopathology and metastatic spread. Leishmania genetic diversity can contribute to variation in polarization and robustness of the immune response through differences in both pathogen sensing by the host and immune evasion by the parasite. In this study, we observed a difference in parasite chemokine suppression between the Leishmania (L.) subgenus and the Viannia (V.) subgenus, which is associated with severe immune-mediated pathology such as mucocutaneous leishmaniasis. While Leishmania (L.) subgenus parasites utilize the virulence factor and metalloprotease glycoprotein-63 (gp63) to suppress the type-1 associated host chemokine CXCL10, L. (V.) panamensis did not suppress CXCL10. To understand the molecular basis for the inter-species variation in chemokine suppression, we used in silico modeling to identify a putative CXCL10-binding site on GP63. The putative CXCL10 binding site is in a region of gp63 under significant positive selection, and it varies from the L. major wild-type sequence in all gp63 alleles identified in the L. (V.) panamensis reference genome. Mutating wild-type L. (L.) major gp63 to the L. (V.) panamensis sequence at the putative binding site impaired cleavage of CXCL10 but not a non-specific protease substrate. Notably, Viannia clinical isolates confirmed that L. (V.) panamensis primarily encodes non-CXCL10-cleaving gp63 alleles. In contrast, L. (V.) braziliensis has an intermediate level of activity, consistent with this species having more equal proportions of both alleles. Our results demonstrate how parasite genetic diversity can contribute to variation in immune responses to Leishmania spp. infection that may play critical roles in the outcome of infection.

Leishmaniasis is a neglected tropical disease caused by Leishmania parasites and spread by the bites of infected sand flies. Most cases of leishmaniasis present as self-healing sores that are resolved by a balanced immune response. Other cases of leishmaniasis involve spread to sites distant from the original bite, including damage of the inner surfaces of the mouth and nose. These cases of leishmaniasis involve an excessive immune response. Leishmania parasites produce virulence factor proteins, such as GP63, to trick the immune system into mounting a weaker response. GP63 specifically degrades signaling proteins that attract and activate certain immune cells. Here, we demonstrate that Leishmania parasite species have evolved to differ in their ability to degrade signaling proteins. In Leishmania species known to cause more immune-mediated tissue damage, the GP63 virulence factor has evolved to not degrade specific immune signaling proteins, thus attracting, and activating more immune cells. Our results demonstrate how diversity among Leishmania parasite species can contribute to variation in immune responses that may play critical roles in the outcome of infection.

Distinct microbial communities colonize tonsillar squamous cell carcinoma

Squamous cell carcinoma of the tonsil is one of the most frequent cancers of the oropharynx. The escalating rate of tonsil cancer during the last decades is associated with the increase of high risk-human papilloma virus (HR-HPV) infections. While the microbiome in oropharyngeal malignant diseases has been characterized to some extent, the microbial colonization of HR-HPV-associated tonsil cancer remains largely unknown. Using 16S rRNA gene amplicon sequencing, we have characterized the microbiome of human palatine tonsil crypts in patients suffering from HR-HPV-associated tonsil cancer in comparison to a control cohort of adult sleep apnea patients. We found an increased abundance of the phyla Firmicutes and Actinobacteria in tumor patients, whereas the abundance of Spirochetes and Synergistetes was significantly higher in the control cohort. Furthermore, the accumulation of several genera such as Veillonella, Streptococcus and Prevotella_7 in tonsillar crypts was associated with tonsil cancer. In contrast, Fusobacterium, Prevotella and Treponema_2 were enriched in sleep apnea patients. Machine learning-based bacterial species analysis indicated that a particular bacterial composition in tonsillar crypts is tumor-predictive. Species-specific PCR-based validation in extended patient cohorts confirmed that differential abundance of Filifactor alocis and Prevotella melaninogenica is a distinct trait of tonsil cancer. This study shows that tonsil cancer patients harbor a characteristic microbiome in the crypt environment that differs from the microbiome of sleep apnea patients on all phylogenetic levels. Moreover, our analysis indicates that profiling of microbial communities in distinct tonsillar niches provides microbiome-based avenues for the diagnosis of tonsil cancer.

A New Target Organ of Leishmania (Viannia) braziliensis Chronic Infection: The Intestine

Leishmania (Viannia) braziliensis is one of the main causes of cutaneous leishmaniasis in the Americas. This species presents genetic polymorphism that can cause destructive lesions in oral, nasal, and oropharyngeal tracts. In a previous study, the parasite caused several histopathological changes to hamster ileums. Our study evaluates immune response components, morphological changes, and effects on neurons in the ileums of hamsters infected by three different strains of L. (V.) braziliensis in two infection periods. For the experiment, we separated hamsters into four groups: a control group and three infected groups. Infected hamsters were euthanized 90- or 120-days post infection. We used three strains of L. (V.) braziliensis: the reference MHOM/BR/1975/M2903 and two strains isolated from patients who had different responses to Glucantime^® treatment (MHOM/BR/2003/2314 and MHOM/BR/2000/1655). After laparotomy, ileums were collected for histological processing, biochemical analysis, and evaluation of neurons in the myenteric and submucosal plexuses of the enteric nervous system (ENS). The results demonstrated the increase of blood leukocytes after the infection. Optical microscopy analysis showed histopathological changes with inflammatory infiltrates, edemas, ganglionitis, and Leishmania amastigotes in the ileums of infected hamsters. We observed changes in the organ histoarchitecture of infected hamsters when compared to control groups, such as thicker muscular and submucosa layers, deeper and wider crypts, and taller and broader villi. The number of intraepithelial lymphocytes and TGF-β-immunoreactive cells increased in all infected groups when compared to the control groups. Mast cells increased with longer infection periods. The infection also caused remodeling of intestinal collagen and morphometry of myenteric and submucosal plexus neurons; but this effect was dependent on infection duration. Our results show that L. (V.) braziliensis infection caused time-dependent alterations in hamster ileums. This was demonstrated by the reduction of inflammatory cells and the increase of tissue regeneration factors at 120 days of infection. The infected groups demonstrated different profiles in organ histoarchitecture, migration of immune cells, and morphometry of ENS neurons. These findings suggest that the small intestine (or at least the ileum) is a target organ for L. (V.) braziliensis infection, as the infection caused changes that were dependent on duration and strain.

Serine proteases profiles of Leishmania (Viannia) braziliensis clinical isolates with distinct susceptibilities to antimony

Glucantime (Sb^V) is the first-line treatment against American Tegumentary Leishmaniasis. Resistance cases to this drug have been reported and related to host characteristics and parasite phenotypes. In this study, 12 Leishmania (Viannia) braziliensis isolates from patients that presented clinical cure (Responders—R) and relapse or therapeutic failure (Non-responders—NR) after treatment with antimony, were analyzed. These parasites were assessed by in vitro susceptibility to Sb^III and Sb^V, serine proteases activity measured with substrate (z-FR-AMC) and specific inhibitors (TLCK, AEBSF and PMSF). In vitro susceptibility of axenic amastigotes to Sb^III showed a significant difference between R and NR groups. The protease assays showed that TLCK inhibited almost 100% of activity in both axenic amastigotes and promastigotes while AEBSF inhibited around 70%, and PMSF showed lower inhibition of some isolates. Principal component and clustering analysis performed with these data yielded one homogeneous cluster with only NR isolates and three heterogeneous clusters with R and NR isolates. Additionally, differential expression of subtilisins (LbrM.13.0860 and LbrM.28.2570) and TXNPx (LbrM.15.1080) was evaluated in promastigotes and axenic amastigotes from both groups. The results showed a higher expression of LbrM.13.0860 and LbrM.15.1080 genes in axenic amastigotes, while LbrM.28.2570 gene had the lowest expression in all isolates, regardless of the parasite form. The data presented here show a phenotypic heterogeneity among the parasites, suggesting that exploration of in vitro phenotypes based on Sb^III and serine proteases profiles can aid in the characterization of L. (V.) braziliensis clinical isolates.