Hugo O. Valdivia

Natural Leishmania Infection of Lutzomyia auraensis in Madre de Dios, Peru, Detected by a Fluorescence Resonance Energy Transfer-Based Real-Time Polymerase Chain Reaction

Leishmania species of the Viannia subgenus are responsible for most cases of New World tegumentary leishmaniasis. However, little is known about the vectors involved in disease transmission in the Amazon regions of Peru. We used a novel real-time polymerase chain reaction (PCR) to assess Leishmania infections in phlebotomines collected in rural areas of Madre de Dios, Peru. A total of 1,299 non-blood fed female sand flies from 33 species were captured by using miniature CDC light traps. Lutzomyia auraensis was the most abundant species (63%) in this area. Seven of 164 pools were positive by PCR for Leishmania by kinetoplast DNA. The real-time PCR identified four Lu. auraensis pools as positive for L. (Viannia) lainsoni and L. (V.) braziliensis. The minimum infection prevalence for Lu. auraensis was estimated to be 0.6% (95% confidence interval = 0.20-1.42%). Further studies are needed to assess the importance of Lu. auraensis in the transmission of New World tegumentary leishmaniasis in hyperendemic areas of Peru.

Chromosomal copy number variation reveals differential levels of genomic plasticity in distinct Trypanosoma cruzi strains

BackgroundTrypanosoma cruzi, the etiologic agent of Chagas disease, is currently divided into six discrete typing units (DTUs), named TcI–TcVI. CL Brener, the reference strain of the T. cruzi genome project, is a hybrid with a genome assembled into 41 putative chromosomes. Gene copy number variation (CNV) is well documented as an important mechanism to enhance gene expression and variability in T. cruzi. Chromosomal CNV (CCNV) is another level of gene CNV in which whole blocks of genes are expanded simultaneously. Although the T. cruzi karyotype is not well defined, several studies have demonstrated a significant variation in the size and content of chromosomes between different T. cruzi strains. Despite these studies, the extent of diversity in CCNV among T. cruzi strains based on a read depth coverage analysis has not been determined.ResultsWe identify the CCNV in T. cruzi strains from the TcI, TcII and TcIII DTUs, by analyzing the depth coverage of short reads from these strains using the 41 CL Brener chromosomes as reference. This study led to the identification of a broader extent of CCNV in T. cruzi than was previously speculated. The TcI DTU strains have very few aneuploidies, while the strains from TcII and TcIII DTUs present a high degree of chromosomal expansions. Chromosome 31, which is the only chromosome that is supernumerary in all six T. cruzi samples evaluated in this study, is enriched with genes related to glycosylation pathways, highlighting the importance of glycosylation to parasite survival.ConclusionsIncreased gene copy number due to chromosome amplification may contribute to alterations in gene expression, which represents a strategy that may be crucial for parasites that mainly depend on post-transcriptional mechanisms to control gene expression.

Comparative genomic analysis of Leishmania (Viannia) peruviana and Leishmania (Viannia) braziliensis.

BackgroundThe Leishmania (Viannia) braziliensis complex is responsible for most cases of New World tegumentary leishmaniasis. This complex includes two closely related species but with different geographic distribution and disease phenotypes, L. (V.) peruviana and L. (V.) braziliensis. However, the genetic basis of these differences is not well understood and the status of L. (V.) peruviana as distinct species has been questioned by some.Here we sequenced the genomes of two L. (V.) peruviana isolates (LEM1537 and PAB-4377) using Illumina high throughput sequencing and performed comparative analyses against the L. (V.) braziliensis M2904 reference genome. Comparisons were focused on the detection of Single Nucleotide Polymorphisms (SNPs), insertions and deletions (INDELs), aneuploidy and gene copy number variations.ResultsWe found 94,070 variants shared by both L. (V.) peruviana isolates (144,079 in PAB-4377 and 136,946 in LEM1537) against the L. (V.) braziliensis M2904 reference genome while only 26,853 variants separated both L. (V.) peruviana genomes.Analysis in coding sequences detected 26,750 SNPs and 1,513 indels shared by both L. (V.) peruviana isolates against L. (V.) braziliensis M2904 and revealed two L. (V.) braziliensis pseudogenes that are likely to have coding potential in L. (V.) peruviana. Chromosomal read density and allele frequency profiling showed a heterogeneous pattern of aneuploidy with an overall disomic tendency in both L. (V.) peruviana isolates, in contrast with a trisomic pattern in the L. (V.) braziliensis M2904 reference.Read depth analysis allowed us to detect more than 368 gene expansions and 14 expanded gene arrays in L. (V.) peruviana, and the likely absence of expanded amastin gene arrays.ConclusionsThe greater numbers of interspecific SNP/indel differences between L. (V.) peruviana and L. (V.) braziliensis and the presence of different gene and chromosome copy number variations support the classification of both organisms as closely related but distinct species.The extensive nucleotide polymorphisms and differences in gene and chromosome copy numbers in L. (V.) peruviana suggests the possibility that these may contribute to some of the unique features of its biology, including a lower pathology and lack of mucosal development.

Selective Whole-Genome Amplification Is a Robust Method That Enables Scalable Whole-Genome Sequencing of Plasmodium vivax from Unprocessed Clinical Samples.

ABSTRACT Whole-genome sequencing (WGS) of microbial pathogens from clinical samples is a highly sensitive tool used to gain a deeper understanding of the biology, epidemiology, and drug resistance mechanisms of many infections. However, WGS of organisms which exhibit low densities in their hosts is challenging due to high levels of host genomic DNA (gDNA), which leads to very low coverage of the microbial genome. WGS of Plasmodium vivax, the most widely distributed form of malaria, is especially difficult because of low parasite densities and the lack of an ex vivo culture system. Current techniques used to enrich P. vivax DNA from clinical samples require significant resources or are not consistently effective. Here, we demonstrate that selective whole-genome amplification (SWGA) can enrich P. vivax gDNA from unprocessed human blood samples and dried blood spots for high-quality WGS, allowing genetic characterization of isolates that would otherwise have been prohibitively expensive or impossible to sequence. We achieved an average genome coverage of 24×, with up to 95% of the P. vivax core genome covered by ≥5 reads. The single-nucleotide polymorphism (SNP) characteristics and drug resistance mutations seen were consistent with those of other P. vivax sequences from a similar region in Peru, demonstrating that SWGA produces high-quality sequences for downstream analysis. SWGA is a robust tool that will enable efficient, cost-effective WGS of P. vivax isolates from clinical samples that can be applied to other neglected microbial pathogens. IMPORTANCE Malaria is a disease caused by Plasmodium parasites that caused 214 million symptomatic cases and 438,000 deaths in 2015. Plasmodium vivax is the most widely distributed species, causing the majority of malaria infections outside sub-Saharan Africa. Whole-genome sequencing (WGS) of Plasmodium parasites from clinical samples has revealed important insights into the epidemiology and mechanisms of drug resistance of malaria. However, WGS of P. vivax is challenging due to low parasite levels in humans and the lack of a routine system to culture the parasites. Selective whole-genome amplification (SWGA) preferentially amplifies the genomes of pathogens from mixtures of target and host gDNA. Here, we demonstrate that SWGA is a simple, robust method that can be used to enrich P. vivax genomic DNA (gDNA) from unprocessed human blood samples and dried blood spots for cost-effective, high-quality WGS. Malaria is a disease caused by Plasmodium parasites that caused 214 million symptomatic cases and 438,000 deaths in 2015. Plasmodium vivax is the most widely distributed species, causing the majority of malaria infections outside sub-Saharan Africa. Whole-genome sequencing (WGS) of Plasmodium parasites from clinical samples has revealed important insights into the epidemiology and mechanisms of drug resistance of malaria. However, WGS of P. vivax is challenging due to low parasite levels in humans and the lack of a routine system to culture the parasites. Selective whole-genome amplification (SWGA) preferentially amplifies the genomes of pathogens from mixtures of target and host gDNA. Here, we demonstrate that SWGA is a simple, robust method that can be used to enrich P. vivax genomic DNA (gDNA) from unprocessed human blood samples and dried blood spots for cost-effective, high-quality WGS.

Comparative genomics of canine-isolated Leishmania (Leishmania) amazonensis from an endemic focus of visceral leishmaniasis in Governador Valadares, southeastern Brazil.

Leishmaniasis is a highly diverse group of diseases caused by kinetoplastid of the genus Leishmania. These parasites are taxonomically diverse, with human pathogenic species separated into two subgenera according to their development site inside the alimentary tract of the sand fly insect vector. The disease encompasses a variable spectrum of clinical manifestations with tegumentary or visceral symptoms. Among the causative species in Brazil, Leishmania (Leishmania) amazonensis is an important etiological agent of human cutaneous leishmaniasis that accounts for more than 8% of all cases in endemic regions. L. (L.) amazonensis is generally found in the north and northeast regions of Brazil. Here, we report the first isolation of L. (L.) amazonensis from dogs with clinical manifestations of visceral leishmaniasis in Governador Valadares, an endemic focus in the southeastern Brazilian State of Minas Gerais where L. (L.) infantum is also endemic. These isolates were characterized in terms of SNPs, chromosome and gene copy number variations, confirming that they are closely related to a previously sequenced isolate obtained in 1973 from the typical Northern range of this species. The results presented in this article will increase our knowledge of L. (L.) amazonensis-specific adaptations to infection, parasite survival and the transmission of this Amazonian species in a new endemic area of Brazil.

Human migration and the spread of malaria parasites to the New World

We examined the mitogenomes of a large global collection of human malaria parasites to explore how and when Plasmodium falciparum and P. vivax entered the Americas. We found evidence of a significant contribution of African and South Asian lineages to present-day New World malaria parasites with additional P. vivax lineages appearing to originate from Melanesia that were putatively carried by the Australasian peoples who contributed genes to Native Americans. Importantly, mitochondrial lineages of the P. vivax-like species P. simium are shared by platyrrhine monkeys and humans in the Atlantic Forest ecosystem, but not across the Amazon, which most likely resulted from one or a few recent human-to-monkey transfers. While enslaved Africans were likely the main carriers of P. falciparum mitochondrial lineages into the Americas after the conquest, additional parasites carried by Australasian peoples in pre-Columbian times may have contributed to the extensive diversity of extant local populations of P. vivax.

Gene and Chromosomal Copy Number Variations as an Adaptive Mechanism Towards a Parasitic Lifestyle in Trypanosomatids

Trypanosomatids are a group of kinetoplastid parasites including some of great public health importance, causing debilitating and life-long lasting diseases that affect more than 24 million people worldwide. Among the trypanosomatids, Trypanosoma cruzi, Trypanosoma brucei and species from the Leishmania genus are the most well studied parasites, due to their high prevalence in human infections. These parasites have an extreme genomic and phenotypic variability, with a massive expansion in the copy number of species-specific multigene families enrolled in host-parasite interactions that mediate cellular invasion and immune evasion processes. As most trypanosomatids are heteroxenous, and therefore their lifecycles involve the transition between different hosts, these parasites have developed several strategies to ensure a rapid adaptation to changing environments. Among these strategies, a rapid shift in the repertoire of expressed genes, genetic variability and genome plasticity are key mechanisms. Trypanosomatid genomes are organized into large directional gene clusters that are transcribed polycistronically, where genes derived from the same polycistron may have very distinct mRNA levels. This particular mode of transcription implies that the control of gene expression operates mainly at post-transcriptional level. In this sense, gene duplications/losses were already associated with changes in mRNA levels in these parasites. Gene duplications also allow the generation of sequence variability, as the newly formed copy can diverge without loss of function of the original copy. Recently, aneuploidies have been shown to occur in several Leishmania species and T. cruzi strains. Although aneuploidies are usually associated with debilitating phenotypes in superior eukaryotes, recent data shows that it could also provide increased fitness in stress conditions and generate drug resistance in unicellular eukaryotes. In this review, we will focus on gene and chromosomal copy number variations and their relevance to the evolution of trypanosomatid parasites.

Exploration of Plasmodium vivax transmission dynamics and recurrent infections in the Peruvian Amazon using whole genome sequencing

BackgroundPlasmodium vivax poses a significant challenge to malaria elimination due to its ability to cause relapsed infections from reactivation of dormant liver parasites called hypnozoites. We analyzed 69 P. vivax whole genome sequences obtained from subjects residing in three different villages along the Peruvian Amazon. This included 23 paired P. vivax samples from subjects who experienced recurrent P. vivax parasitemia following observed treatment with chloroquine and primaquine.MethodsGenomic DNA was extracted from whole blood samples collected from subjects. P. vivax DNA was enriched using selective whole genome amplification and whole genome sequencing. We used single nucleotide polymorphisms (SNPs) from the core P. vivax genome to determine characteristics of the parasite population using discriminant analysis of principal components, maximum likelihood estimation of individual ancestries, and phylogenetic analysis. We estimated the relatedness of the paired samples by calculating the number of segregating sites and using a hidden Markov model approach to estimate identity by descent.ResultsWe present a comprehensive dataset of population genetics of Plasmodium vivax in the Peruvian Amazonian. We define the parasite population structure in this region and demonstrate a novel method for distinguishing homologous relapses from reinfections or heterologous relapses with improved accuracy. The parasite population in this area was quite diverse with an estimated five subpopulations and evidence of a highly heterogeneous ancestry of some of the isolates, similar to previous analyses of P. vivax in this region. Pairwise comparison of recurrent infections determined that there were 12 homologous relapses and 3 likely heterologous relapses with highly related parasites. To the best of our knowledge, this is the first large-scale study to evaluate recurrent P. vivax infections using whole genome sequencing.ConclusionsWhole genome sequencing is a high-resolution tool that can identify P. vivax homologous relapses with increased sensitivity, while also providing data about drug resistance and parasite population genetics. This information is important for evaluating the efficacy of known and novel antirelapse medications in endemic areas and thus advancing the campaign to eliminate malaria.

TipMT: Identification of PCR-based taxon-specific markers

BackgroundMolecular genetic markers are one of the most informative and widely used genome features in clinical and environmental diagnostic studies. A polymerase chain reaction (PCR)-based molecular marker is very attractive because it is suitable to high throughput automation and confers high specificity. However, the design of taxon-specific primers may be difficult and time consuming due to the need to identify appropriate genomic regions for annealing primers and to evaluate primer specificity.ResultsHere, we report the development of a Tool for Identification of Primers for Multiple Taxa (TipMT), which is a web application to search and design primers for genotyping based on genomic data. The tool identifies and targets single sequence repeats (SSR) or orthologous/taxa-specific genes for genotyping using Multiplex PCR. This pipeline was applied to the genomes of four species of Leishmania (L. amazonensis, L. braziliensis, L. infantum and L. major) and validated by PCR using artificial genomic DNA mixtures of the Leishmania species as templates. This experimental validation demonstrates the reliability of TipMT because amplification profiles showed discrimination of genomic DNA samples from Leishmania species.ConclusionsThe TipMT web tool allows for large-scale identification and design of taxon-specific primers and is freely available to the scientific community at http://200.131.37.155/tipMT/.

Detection of multiple circulating Leishmania species in Lutzomyia longipalpis in the city of Governador Valadares, southeastern Brazil

Background Leishmaniasis encompasses a group of diverse clinical diseases caused by protozoan parasites of the Leishmania genus. This disease is a major public health problem in the New World affecting people exposed in endemic regions. The city of Governador Valadares (Minas Gerais/Brazil) is a re-emerging area for visceral leishmaniasis, with 191 human cases reported from 2008 to 2017 and a lethality rate of 14.7%. The transmission of the parasite occurs intensely in this region with up to 22% of domestic dogs with positive serology for the visceral form. Lu. longipalpis is one of the most abundant sand fly species in this area. Despite this scenario, so far there is no information regarding the circulating Leishmania species in the insect vector Lutzomyia longipalpis in this focus. Methodology/Principal Findings We collected 616 female Lutzomyia longipalpis sand flies between January and September 2015 in the Vila Parque Ibituruna neighborhood (Governador Valadares/MG), which is located on a transitional area between the sylvatic and urban environments with residences built near a preserved area. After DNA extraction of individual sand flies, the natural Leishmania infections in Lu. longipalpis were detected by end-point PCR, using primers derived from kDNA sequences, specific for L. (Leishmania) or L. (Viannia) subgenus. The sensitivity of these PCR reactions was 0.1 pg of DNA for each Leishmania subgenus and the total infection rate of 16.2% (100 positive specimens). Species-specific PCR detected the presence of multiple Leishmania species in infected Lu. longipalpis specimens in Governador Valadares, including L. amazonensis (n=3), L. infantum (n=28), L. (Viannia) spp. (n=20), coinfections with L. infantum and L. (Viannia) spp. (n=5), and L. (Leishmania) spp (n=44). Conclusions Our results demonstrate that multiple Leishmania species circulate in Lu. longipalpis in Governador Valadares and reveal a potential increasing risk of transmission of the different circulating parasite species. This information is a key factor for planning surveillance and effective control strategies against leishmaniasis in this endemic focus. Author summary Leishmaniasis is a neglected tropical disease transmitted to mammals by the bite of sand flies infected with parasites of the Leishmania genus. This disease affects millions of people in various regions of the world, including Brazil. The municipality of Governador Valadares (Minas Gerais/Brazil) is a re-emergent focus of intense transmission of leishmaniasis, with a high number of human cases and a high prevalence of infected domestic dogs. To develop better leishmaniasis control strategies for the region, we performed a surveillance study of Lu. longipalpis, the main vector of visceral leishmaniasis in Brazil, and identified circulating species of Leishmania in this insect vector. We estimate that the natural infection rate of Lu. longipalpis for these parasites was of 16.2% in the study area. We also detected the presence of multiple circulating Leishmania species (L. amazonensis, L. infantum and Viannia subgenus) in Lu. longipalpis in Governador Valadares city, including 5 sand flies coinfected with L. infantum and L. (Viannia). Thus, our results reinforce the need for a rigid and systematic control of the sand flies monitoring in this area, due to the potential risk of transmission of different species of the Leishmania parasites.