Far-East Asian Toxoplasma isolates share ancestry with North and South/Central American recombinant lineages.

Toxoplasma gondii is a global protozoan pathogen. Clonal lineages predominate in Europe, North America, Africa, and China, whereas highly recombinant parasites are endemic in South/Central America. Far East Asian T. gondii isolates are not included in current global population genetic structure analyses at WGS resolution. Here we report a genome-wide population study that compared eight Japanese and two Chinese isolates against representative worldwide T. gondii genomes using POPSICLE, a novel population structure analyzing software. Also included were 7 genomes resurrected from non-viable isolates by target enrichment sequencing. Visualization of the genome structure by POPSICLE shows a mixture of Chinese haplogroup (HG) 13 haploblocks introgressed within the genomes of Japanese HG2 and North American HG12. Furthermore, two ancestral lineages were identified in the Japanese strains; one lineage shares a common ancestor with HG11 found in both Japanese strains and North American HG12. The other ancestral lineage, found in T. gondii isolates from a small island in Japan, is admixed with genetically diversified South/Central American strains. Taken together, this study suggests multiple ancestral links between Far East Asian and American T. gondii strains and provides insight into the transmission history of this cosmopolitan organism.

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.

Protocol for design, construction, and selection of genome phage (gPhage) display libraries.

This protocol describes the genomic phage (gPhage) display platform, a large-scale antigen and epitope mapping technique. We constructed a gPhage display peptide library of a eukaryotic organism, Trypanosoma cruzi (causative agent of Chagas disease), to map the antibody response landscape against the parasite. Here, we used an organism with a relatively large but intronless genome, although future applications could include other prevalent or (re)emerging infectious organisms carrying genomes with a limited number of introns. For complete details on the use and execution of this protocol, please refer to Teixeira et al. (2021).

Genomics and functional genomics in Leishmania and Trypanosoma cruzi: statuses, challenges and perspectives.

The availability of Trypanosomatid genomic data in public databases has opened myriad experimental possibilities that have contributed to a more comprehensive understanding of the biology of these parasites and their interactions with hosts. In this review, after brief remarks on the history of the Trypanosoma cruzi and Leishmania genome initiatives, we present an overview of the relevant contributions of genomics, transcriptomics and functional genomics, discussing the primary obstacles, challenges, relevant achievements and future perspectives of these technologies.

Genomics and functional genomics in Leishmania and Trypanosoma cruzi: statuses, challenges and perspectives

The availability of Trypanosomatid genomic data in public databases has opened myriad experimental possibilities that have contributed to a more comprehensive understanding of the biology of these parasites and their interactions with hosts. In this review, after brief remarks on the history of the Trypanosoma cruzi and Leishmania genome initiatives, we present an overview of the relevant contributions of genomics, transcriptomics and functional genomics, discussing the primary obstacles, challenges, relevant achievements and future perspectives of these technologies.

Identity-by-descent with uncertainty characterises connectivity of Plasmodium falciparum populations on the Colombian-Pacific coast.

Characterising connectivity between geographically separated biological populations is a common goal in many fields. Recent approaches to understanding connectivity between malaria parasite populations, with implications for disease control efforts, have used estimates of relatedness based on identity-by-descent (IBD). However, uncertainty around estimated relatedness has not been accounted for. IBD-based relatedness estimates with uncertainty were computed for pairs of monoclonal Plasmodium falciparum samples collected from five cities on the Colombian-Pacific coast where long-term clonal propagation of P. falciparum is frequent. The cities include two official ports, Buenaventura and Tumaco, that are separated geographically but connected by frequent marine traffic. Fractions of highly-related sample pairs (whose classification using a threshold accounts for uncertainty) were greater within cities versus between. However, based on both highly-related fractions and on a threshold-free approach (Wasserstein distances between parasite populations) connectivity between Buenaventura and Tumaco was disproportionally high. Buenaventura-Tumaco connectivity was consistent with transmission events involving parasites from five clonal components (groups of statistically indistinguishable parasites identified under a graph theoretic framework). To conclude, P. falciparum population connectivity on the Colombian-Pacific coast abides by accessibility not isolation-by-distance, potentially implicating marine traffic in malaria transmission with opportunities for targeted intervention. Further investigations are required to test this hypothesis. For the first time in malaria epidemiology (and to our knowledge in ecological and epidemiological studies more generally), we account for uncertainty around estimated relatedness (an important consideration for studies that plan to use genotype versus whole genome sequence data to estimate IBD-based relatedness); we also use threshold-free methods to compare parasite populations and identify clonal components. Threshold-free methods are especially important in analyses of malaria parasites and other recombining organisms with mixed mating systems where thresholds do not have clear interpretation (e.g. due to clonal propagation) and thus undermine the cross-comparison of studies.

Design of a AFLP-PCR and PCR-RFLP test that identify the majority of discrete typing units of Trypanosoma cruzi.

BACKGROUND: Chagas disease, caused by the intracellular parasite Trypanosoma cruzi, is one of the most important parasitological infections in the Americas. It is estimated to infect approximately 6 million people from mostly low income countries in Latin America, although recent infections have been reported in southern US states. Several studies have described an extensive genetic diversity among T. cruzi isolates throughout its geographic distribution in the American continent. This diversity has been correlated with the pathology developed during an infection. However, due to a lack of a single reliable test, current diagnosis practices of the disease are not straightforward since several different tests are applied. The use of current genomic sequence data allows for the selection of molecular markers (MM) that have the ability to identify the Discrete Typing Unit (DTU) of T. cruzi in a given infection, without the need of any sequencing reaction. METHODOLOGY/PRINCIPAL FINDINGS: Applying three criteria on the genomic sequencing data of four different phylogenetic lineages of T. cruzi, we designed several molecular tests that can be used for the molecular typing of the parasite. The criteria used were: (1) single-copy orthologs of T. cruzi, (2) T. cruzi unique loci, and (3) T. cruzi polymorphic loci. All criteria combined allowed for the selection of 15 MM, 12 of which were confirmed to be functional and replicable in the laboratory with sylvatic samples. Furthermore, one MM produced distinct polymerase chain reaction (PCR) amplicon sizes among distinct T. cruzi DTUs, allowing the use of a AFLP-PCR test to distinguish DTUs I, II/IV, V and VI. Whereas two MM can differentiate DTUs I, II, IV and V/VI out of the six current DTUs with a PCR-RFLP test. CONCLUSIONS/SIGNIFICANCE: The designed molecular tests provide a practical and inexpensive molecular typing test for the majority of DTUs of T. cruzi, excluding the need to perform any sequencing reaction. This provides the scientific community with an additional specific, quick and inexpensive test that can enhance the understanding of the correlation between the DTU of T. cruzi and the pathology developed during the infection.

Fatty acid-binding proteins in Echinococcus spp.: the family has grown.

Fatty acid-binding proteins (FABPs) are small intracellular proteins that reversibly bind fatty acids and other hydrophobic ligands. In cestodes, due to their inability to synthesise fatty acids de novo, FABPs have been proposed as essential proteins, and thus, as possible drug targets and/or carriers against these parasites. We performed data mining in Echinococcus multilocularis and Echinococcus granulosus genomes in order to test whether this family of proteins is more complex than previously reported. By exploring the genomes of E. multilocularis and E. granulosus, six genes coding for FABPs were found in each organism. In the case of E. granulosus, all of them have different coding sequences, whereas in E. multilocularis, two of the genes code for the same protein. Remarkably, one of the genes (in both cestodes) encodes a FABP with a C-terminal extension unusual for this family of proteins. The newly described genes present variations in their structure in comparison with previously described FABP genes in Echinococcus spp. The coding sequences for E. multilocularis were validated by cloning and sequencing. Moreover, differential expression patterns of FABPs were observed at different stages of the life cycle of E. multilocularis by exploring transcriptomic data from several sources. In summary, FABP family in cestodes is far more complex than previously thought and includes new members that seem to be only present in flatworms.

A molecular barcode to inform the geographical origin and transmission dynamics of Plasmodium vivax malaria.

Although Plasmodium vivax parasites are the predominant cause of malaria outside of sub-Saharan Africa, they not always prioritised by elimination programmes. P. vivax is resilient and poses challenges through its ability to re-emerge from dormancy in the human liver. With observed growing drug-resistance and the increasing reports of life-threatening infections, new tools to inform elimination efforts are needed. In order to halt transmission, we need to better understand the dynamics of transmission, the movement of parasites, and the reservoirs of infection in order to design targeted interventions. The use of molecular genetics and epidemiology for tracking and studying malaria parasite populations has been applied successfully in P. falciparum species and here we sought to develop a molecular genetic tool for P. vivax. By assembling the largest set of P. vivax whole genome sequences (n = 433) spanning 17 countries, and applying a machine learning approach, we created a 71 SNP barcode with high predictive ability to identify geographic origin (91.4%). Further, due to the inclusion of markers for within population variability, the barcode may also distinguish local transmission networks. By using P. vivax data from a low-transmission setting in Malaysia, we demonstrate the potential ability to infer outbreak events. By characterising the barcoding SNP genotypes in P. vivax DNA sourced from UK travellers (n = 132) to ten malaria endemic countries predominantly not used in the barcode construction, we correctly predicted the geographic region of infection origin. Overall, the 71 SNP barcode outperforms previously published genotyping methods and when rolled-out within new portable platforms, is likely to be an invaluable tool for informing targeted interventions towards elimination of this resilient human malaria.

Calpains of Leishmania braziliensis: genome analysis, differential expression, and functional analysis.

BACKGROUND: Calpains are proteins belonging to the multi-gene family of calcium-dependent cysteine peptidases that undergo tight on/off regulation, and uncontrolled proteolysis of calpains is associated with severe human pathologies. Calpain orthologues are expanded and diversified in the trypanosomatids genome. OBJECTIVES: Here, we characterised calpains in Leishmania braziliensis, the main causative agent of cutaneous leishmaniasis in Brazil. METHODS/FINDINGS: In total, 34 predicted calpain-like genes were identified. After domain structure evaluation, reverse transcription-quantitative polymerase chain reaction (RT-qPCR) during in vitro metacyclogenesis revealed (i) five genes with enhanced expression in the procyclic stage, (ii) one augmented gene in the metacyclic stage, and (iii) one procyclic-exclusive transcript. Western blot analysis revealed that an antibody against a consensus-conserved peptide reacted with multiple calpain-like proteins, which is consistent with the multi-gene family characteristic. Flow cytometry and immunocytochemistry analyses revealed the presence of calpain-like molecules mainly in the cytoplasm, to a lesser extent in the plasma membrane, and negligible levels in the nucleus, which are all consistent with calpain localisation. Eventually, the calpain inhibitor MDL28170 was used for functional studies revealing (i) a leishmaniostatic effect, (ii) a reduction in the association index in mouse macrophages, (iii) ultra-structural alterations conceivable with autophagy, and (iv) an enhanced expression of the virulence factor GP63. CONCLUSION: This report adds novel insights into the domain structure, expression, and localisation of L. braziliensis calpain-like molecules.

Genome Sequence Data of Six Isolates of Phytophthora capsici from Mexico.

Phytophthora capsici is an oomycete plant pathogen with a wide host range. Worldwide, P. capsici is known for causing the principal disease of chili pepper crops. Our goal was to expand the available genome resources for this diverse pathogen by generating whole-genome sequences for six isolates of P. capsici from Mexico.

Nanopore Sequencing Significantly Improves Genome Assembly of the Protozoan Parasite Trypanosoma cruzi.

Chagas disease was described by Carlos Chagas, who first identified the parasite Trypanosoma cruzi from a 2-year-old girl called Berenice. Many T. cruzi sequencing projects based on short reads have demonstrated that genome assembly and downstream comparative analyses are extremely challenging in this species, given that half of its genome is composed of repetitive sequences. Here, we report de novo assemblies, annotation, and comparative analyses of the Berenice strain using a combination of Illumina short reads and MinION long reads. Our work demonstrates that Nanopore sequencing improves T. cruzi assembly contiguity and increases the assembly size in ∼16 Mb. Specifically, we found that assembly improvement also refines the completeness of coding regions for both single-copy genes and repetitive transposable elements. Beyond its historical and epidemiological importance, Berenice constitutes a fundamental resource because it now constitutes a high-quality assembly available for TcII (clade C), a prevalent lineage causing human infections in South America. The availability of Berenice genome expands the known genetic diversity of these parasites and reinforces the idea that T. cruzi is intraspecifically divided in three main clades. Finally, this work represents the introduction of Nanopore technology to resolve complex protozoan genomes, supporting its subsequent application for improving trypanosomatid and other highly repetitive genomes.

Technical and economic efficiency of methods for extracting genomic DNA from Meloidogyne javanica.

Plant parasitic nematodes reduce the production of agricultural crops. Species diagnosis is essential to predict losses, determine economic damage levels and develop integrated pest management programs. DNA extraction techniques need to be improved for precise and rapid molecular diagnosis of nematodes. The objective of the present study was to evaluate the efficiency of DNA extraction and amplification by PCR, cost and execution time by Chelex, Worm Lysis Buffer Method (WLB), Holterman Lysis Buffer Method (HLB) and FastDNA methods for nematodes of the Meloidogyne genus. The qualitative and quantitative efficiency of DNA extraction varied between methods. The band size of the amplified PCR product with WLB, Chelex and HLB methods was 590 bp. Extraction with the FastDNA is not recommended for DNA extraction from nematodes because it results in a low DNA concentration without bands in PCR amplification, besides presenting high cost. The efficiency of the WLB method to extracting DNA from Meloidogyne javanica was greater, ensuring a higher concentration and purity of the extracted material and guaranteeing lower costs and greater ease of PCR amplification.

Calpains of Leishmania braziliensis: genome analysis, differential expression, and functional analysis

BACKGROUND Calpains are proteins belonging to the multi-gene family of calcium-dependent cysteine peptidases that undergo tight on/off regulation, and uncontrolled proteolysis of calpains is associated with severe human pathologies. Calpain orthologues are expanded and diversified in the trypanosomatids genome. OBJECTIVES Here, we characterised calpains in Leishmania braziliensis, the main causative agent of cutaneous leishmaniasis in Brazil. METHODS/FINDINGS In total, 34 predicted calpain-like genes were identified. After domain structure evaluation, reverse transcription-quantitative polymerase chain reaction (RT-qPCR) during in vitro metacyclogenesis revealed (i) five genes with enhanced expression in the procyclic stage, (ii) one augmented gene in the metacyclic stage, and (iii) one procyclic-exclusive transcript. Western blot analysis revealed that an antibody against a consensus-conserved peptide reacted with multiple calpain-like proteins, which is consistent with the multi-gene family characteristic. Flow cytometry and immunocytochemistry analyses revealed the presence of calpain-like molecules mainly in the cytoplasm, to a lesser extent in the plasma membrane, and negligible levels in the nucleus, which are all consistent with calpain localisation. Eventually, the calpain inhibitor MDL28170 was used for functional studies revealing (i) a leishmaniostatic effect, (ii) a reduction in the association index in mouse macrophages, (iii) ultra-structural alterations conceivable with autophagy, and (iv) an enhanced expression of the virulence factor GP63. CONCLUSION This report adds novel insights into the domain structure, expression, and localisation of L. braziliensis calpain-like molecules.

Independent regulation of Plasmodium falciparum rif gene promoters.

All Plasmodium species express variant antigens which may mediate immune escape in the vertebrate host. In Plasmodium falciparum, the rif gene family encodes variant antigens which are partly exposed on the infected red blood cell surface and may function as virulence factors. Not all rif genes are expressed at the same time and it is unclear what controls rif gene expression. In this work, we addressed global rif transcription using plasmid vectors with two drug resistance markers, one controlled by a rif 5' upstream region and the second by a constitutively active promoter. After spontaneous integration into the genome of one construct, we observed that the resistance marker controlled by the rif 5' upstream region was expressed dependent on the applied drug pressure. Then, the global transcription of rif genes in these transfectants was compared in the presence or absence of drugs. The relative transcript quantities of all rif loci did not change profoundly between strains grown with or without drug. We conclude that either there is no crosstalk between rif loci or that the elusive system of allelic exclusion of rif gene transcription is not controlled by their 5' upstream region alone.

Identification of Leishmania (Viannia) species and clinical isolates of Leishmania (Leishmania) amazonensis from Brazil using PCR-RFLP of the heat-shock protein 70 gene reveals some unexpected observations.

Hsp70 is a cytoplasmic heat-shock protein, encoded by a multicopy tandemly repeated gene that has recently been gaining popularity as a valuable marker for typing Leishmania species. In this study, we used a previously described hsp70 PCR-RFLP method for identifying Brazilian Leishmania isolates. We identified two distinct L. (L.) amazonensis hsp70 alleles that resulted in two different RFLP patterns. Also, we found RFLP polymorphisms amongst L. (Viannia) naiffi strains. The profiles of both L. (V.) shawi and L. (V.) lindenbergi were very similar to those of other L. (Viannia) species. The observations described herein reflect the polymorphism found within species of Leishmania and indicate that results from this hsp70 PCR-RFLP method should be used with caution when typing isolates from clinical cases of leishmaniasis and Leishmania species from Brazil.

In silico analysis of the EF-hand proteins in the genome of Giardia intestinalis assembly A.

Giardia intestinalis is a parasite that inhabits the small intestine of humans and other mammals, causing a disease that can manifest itself with acute diarrhea. This parasite is an early divergent eukaryote with a compact genome and a life cycle composed of two distinct cell types: the trophozoite, the replicative form, and the cyst, the infectious form. Signal transduction pathways implicated in differentiation processes of G. intestinalis are largely unknown. Calcium, considered an essential messenger in cell signaling, has been shown to regulate a myriad of key cell processes including metabolism, motility, and exocytosis, among other important functions, through calcium-binding proteins (CaBPs). The most important and largest family of CaBPs is the EF-hand protein family. To investigate the nature of calcium signaling pathways present in this protozoan, an in silico analysis of the genome to identify genes encoding EF-hand proteins was undertaken. Twenty-eight sequences containing EF-hand domains were found; most of which have only a pair of domains, and half of the sequences were divergent or unique to Giardia. In addition, the transcription pattern for eight genes encoding EF-hand proteins was assessed during encystation. It was found that all the genes were differentially transcribed suggesting a different function in this process. The in silico results suggest that in G. intestinalis, calcium is involved in the regulation of protein phosphorylation through kinases and phosphatases.

Population genetic analysis of Giardia duodenalis: genetic diversity and haplotype sharing between clinical and environmental sources.

Giardia duodenalis is a flagellated intestinal protozoan responsible for infections in various hosts including humans and several wild and domestic animals. Few studies have correlated environmental contamination and clinical infections in the same region. The aim of this study was to compare groups of Giardia duodenalis from clinical and environmental sources through population genetic analyses to verify haplotype sharing and the degree of genetic similarity among populations from clinical and environmental sources in the metropolitan region of Campinas. The results showed high diversity of haplotypes and substantial genetic similarity between clinical and environmental groups of G. duodenalis. We demonstrated sharing of Giardia genotypes among the different populations studied. The comparison between veterinary and human sequences led us to identify new zoonotic genotypes, including human isolates from genetic assemblage C. The application of a population genetic analysis in epidemiological studies allows quantification of the degree of genetic similarity among populations of Giardia duodenalis from different sources of contamination. The genetic similarity of Giardia isolates among human, veterinary, and environmental groups reinforced the correlation between clinical and environmental isolates in this region, which is of great importance for public health.

Whole Genome Sequencing of Field Isolates Reveals Extensive Genetic Diversity in Plasmodium vivax from Colombia.

Plasmodium vivax is the most prevalent malarial species in South America and exerts a substantial burden on the populations it affects. The control and eventual elimination of P. vivax are global health priorities. Genomic research contributes to this objective by improving our understanding of the biology of P. vivax and through the development of new genetic markers that can be used to monitor efforts to reduce malaria transmission. Here we analyze whole-genome data from eight field samples from a region in Cordóba, Colombia where malaria is endemic. We find considerable genetic diversity within this population, a result that contrasts with earlier studies suggesting that P. vivax had limited diversity in the Americas. We also identify a selective sweep around a substitution known to confer resistance to sulphadoxine-pyrimethamine (SP). This is the first observation of a selective sweep for SP resistance in this species. These results indicate that P. vivax has been exposed to SP pressure even when the drug is not in use as a first line treatment for patients afflicted by this parasite. We identify multiple non-synonymous substitutions in three other genes known to be involved with drug resistance in Plasmodium species. Finally, we found extensive microsatellite polymorphisms. Using this information we developed 18 polymorphic and easy to score microsatellite loci that can be used in epidemiological investigations in South America.

Expanding the tool box for genetic manipulation of Trypanosoma cruzi.

Trypanosoma cruzi is a protozoan parasite that causes Chagas disease, an illness that affects 6-7 million people and for which there is no effective drug therapy or vaccine. The publication of its complete genome sequence allowed a rapid advance in molecular studies including in silico screening of genes involved with pathogenicity as well as molecular targets for the development of new diagnostic methods, drug therapies and prophylactic vaccines. Alongside with in silico genomic analyses, methods to study gene function in this parasite such as gene deletion, overexpression, mutant complementation and reporter gene expression have been largely explored. More recently, the use of genome-wide strategies is producing a shift towards a global perspective on gene function studies, with the examination of the expression and biological roles of gene networks in different stages of the parasite life cycle and under different contexts of host parasite interactions. Here we describe the molecular tools and protocols currently available to perform genetic manipulation of the T. cruzi genome, with emphasis on recently described strategies of gene editing that will facilitate large-scale functional genomic analyses. These new methodologies are long overdue, since more efficient protocols for genetic manipulation in T. cruzi are urgently needed for a better understanding of the biology of this parasite and molecular processes involved with the complex and often harmful, interaction with its human host.