Probing the genetic population structure of Trypanosoma cruzi with polymorphic microsatellites

Molecular Markers for the Phylogenetic Reconstruction of Trypanosoma cruzi: A Quantitative Review

Trypanosoma cruzi is the parasite responsible for Chagas disease, which has a significant amount of genetic diversification among the species complex. Many efforts are routinely made to characterize the genetic lineages of T. cruzi circulating in a particular geographic area. However, the genetic loci used to typify the genetic lineages of T. cruzi have not been consistent between studies. We report a quantitative analysis of the phylogenetic power that is acquired from the commonly used genetic loci that are employed for the typification of T. cruzi into its current taxonomic nomenclature. Based on three quantitative criteria (the number of phylogenetic informative characters, number of available reference sequences in public repositories, and accessibility to DNA sequences for their use as outgroup sequences), we examine and discuss the most appropriate genetic loci for the genetic typification of T. cruzi. Although the mini-exon gene is by far the locus that has been most widely used, it is not the most appropriate marker for the typification of T. cruzi based on the construction of a resolved phylogenetic tree. Overall, the mitochondrial COII-NDI locus stands out as the best molecular marker for this purpose, followed by the Cytochrome b and the Lathosterol oxidase genes.

New insights into Trypanosoma cruzi genetic diversity, and its influence on parasite biology and clinical outcomes

Chagas disease, caused by Trypanosoma cruzi, remains a serious public health problem worldwide. The parasite was subdivided into six distinct genetic groups, called "discrete typing units" (DTUs), from TcI to TcVI. Several studies have indicated that the heterogeneity of T. cruzi species directly affects the diversity of clinical manifestations of Chagas disease, control, diagnosis performance, and susceptibility to treatment. Thus, this review aims to describe how T. cruzi genetic diversity influences the biology of the parasite and/or clinical parameters in humans. Regarding the geographic dispersion of T. cruzi, evident differences were observed in the distribution of DTUs in distinct areas. For example, TcII is the main DTU detected in Brazilian patients from the central and southeastern regions, where there are also registers of TcVI as a secondary T. cruzi DTU. An important aspect observed in previous studies is that the genetic variability of T. cruzi can impact parasite infectivity, reproduction, and differentiation in the vectors. It has been proposed that T. cruzi DTU influences the host immune response and affects disease progression. Genetic aspects of the parasite play an important role in determining which host tissues will be infected, thus heavily influencing Chagas disease's pathogenesis. Several teams have investigated the correlation between T. cruzi DTU and the reactivation of Chagas disease. In agreement with these data, it is reasonable to suppose that the immunological condition of the patient, whether or not associated with the reactivation of the T. cruzi infection and the parasite strain, may have an important role in the pathogenesis of Chagas disease. In this context, understanding the genetics of T. cruzi and its biological and clinical implications will provide new knowledge that may contribute to additional strategies in the diagnosis and clinical outcome follow-up of patients with Chagas disease, in addition to the reactivation of immunocompromised patients infected with T. cruzi.

Fifteen Years after the Definition of Trypanosoma cruzi DTUs: What Have We Learned?

Trypanosoma cruzi, the protozoan causative of Chagas disease (ChD), exhibits striking genetic and phenotypic intraspecific diversity, along with ecoepidemiological complexity. Human-pathogen interactions lead to distinct clinical presentations of ChD. In 2009, an international consensus classified T. cruzi strains into six discrete typing units (DTUs), TcI to TcVI, later including TcBat, and proposed reproducible genotyping schemes for DTU identification. This article aims to review the impact of classifying T. cruzi strains into DTUs on our understanding of biological, ecoepidemiological, and pathogenic aspects of T. cruzi. We will explore the likely origin of DTUs and the intrinsic characteristics of each group of strains concerning genome organization, genomics, and susceptibility to drugs used in ChD treatment. We will also provide an overview of the association of DTUs with mammalian reservoirs, and summarize the geographic distribution, and the clinical implications, of prevalent specific DTUs in ChD patients. Throughout this review, we will emphasize the crucial roles of both parasite and human genetics in defining ChD pathogenesis and chemotherapy outcome.

Wide reference databases for typing Trypanosoma cruzi based on amplicon sequencing of the minicircle hypervariable region

BACKGROUND

Trypanosoma cruzi, the etiological agent of Chagas Disease, exhibits remarkable genetic diversity and is classified into different Discrete Typing Units (DTUs). Strain typing techniques are crucial for studying T. cruzi, because their DTUs have significant biological differences from one another. However, there is currently no methodological strategy for the direct typing of biological materials that has sufficient sensitivity, specificity, and reproducibility. The high diversity and copy number of the minicircle hypervariable regions (mHVRs) makes it a viable target for typing.

METHODOLOGY/PRINCIPAL FINDINGS

Approximately 24 million reads obtained by amplicon sequencing of the mHVR were analyzed for 62 strains belonging to the six main T. cruzi DTUs. To build reference databases of mHVR diversity for each DTU and to evaluate this target as a typing tool. Strains of the same DTU shared more mHVR clusters than strains of different DTUs, and clustered together. Different identity thresholds were used to build the reference sets of the mHVR sequences (85% and 95%, respectively). The 95% set had a higher specificity and was more suited for detecting co-infections, whereas the 85% set was excellent for identifying the primary DTU of a sample. The workflow's capacity for typing samples obtained from cultures, a set of whole-genome data, under various simulated PCR settings, in the presence of co-infecting lineages and for blood samples was also assessed.

CONCLUSIONS/SIGNIFICANCE

We present reference databases of mHVR sequences and an optimized typing workflow for T. cruzi including a simple online tool for deep amplicon sequencing analysis (https://ntomasini.github.io/cruzityping/). The results show that the workflow displays an equivalent resolution to that of the other typing methods. Owing to its specificity, sensitivity, relatively low cost, and simplicity, the proposed workflow could be an alternative for screening different types of samples.

Unveiling challenges in real-time PCR strategies for detecting treatment failure: observations from clinical trials on chronic Chagas disease

Chagas disease (CD) caused by Trypanosoma cruzi remains a Neglected Tropical Disease with limited access to diagnosis and treatment, particularly for chronically infected patients. Clinical trials are underway to improve treatment using new drugs or different regimens, and Real-Time PCR is used to assess the parasitological response as a surrogate biomarker. However, PCR-based strategies have limitations due to the complex nature of T. cruzi infection. The parasite exhibits asynchronous replication, different strains and clones, and diverse tissue tropism, making it challenging to determine optimal timeline points for monitoring treatment response. This mini-review explores factors that affect PCR-based monitoring and summarizes the endpoints used in clinical trials for detecting treatment failure. Serial sampling and cumulative PCR results may improve sensitivity in detecting parasitemia and treatment failure in these trials.

Genetic variation and clonal diversity in floating aquatic plants: Comparative genomic analysis of water hyacinth species in their native range

Many eukaryotic organisms reproduce by sexual and asexual reproduction. Genetic diversity in populations can be strongly dependent on the relative importance of these two reproductive modes. Here, we compare the amounts and patterns of genetic diversity in related water hyacinths that differ in their propensity for clonal propagation - highly clonal Eichhornia crassipes and moderately clonal E. azurea (Pontederiaceae). Our comparisons involved genotype-by-sequencing (GBS) of 137 E. crassipes ramets from 60 locations (193,495 nucleotide sites) and 118 E. azurea ramets from 53 locations (198,343 nucleotide sites) among six hydrological basins in central South America, the native range of both species. We predicted that because of more prolific clonal propagation, E. crassipes would exhibit lower clonal diversity than E. azurea. This prediction was supported by all measures of clonal diversity that we examined. Eichhornia crassipes also had a larger excess of heterozygotes at variant sites, another signature of clonality. However, genome-wide heterozygosity was not significantly different between the species. Eichhornia crassipes had weaker spatial genetic structure and lower levels of differentiation among hydrological basins than E. azurea, probably because of higher clonality and more extensive dispersal of its free-floating life form. Our findings for E. crassipes contrast with earlier studies from the invasive range which have reported very low levels of clonal diversity and extensive geographic areas of genetic uniformity.

Microsatellite-Based Genotyping, Analysis of Population Structure, Presence of Trichomonas vaginalis Virus (TVV) and Mycoplasma hominis in T. vaginalis Isolates from Southwest of Turkey

Background:

The present study aimed to determine genetic diversity of Trichomonas vaginalis (T. vaginalis) isolates with microsatellite markers in Turkey (Nov 2015 to 2016) and to create a web-based microsatellite typing (MT) approach for the global interpretation of the data. In addition, the endosymbiosis of Mycoplasma hominis (M. hominis) and T. vaginalis virus (TVV) in the isolates was also examined.

Methods:

The allele sizes for each locus were calculated and microsatellite types were determined according to the allele profiles. The population structure was examined with Bayesian clustering method. A website (http://mttype.adu.edu.tr) was created for collection and sharing of microsatellite data. Presence of TVV and M. hominis in T. vaginalis isolates were investigated with electrophoresis and PCR.

Results:

Of 630 vaginal samples T. vaginalis was detected in 30 (4.7%) and those were used for further analysis. The structure produced by a clustering algorithm revealed eight genetic groups. The typing of isolates according to microsatellites revealed 23 different microsatellite types. Three clones were determined among isolates (MT10 16.7%; MT18 10% and MT3 6.7%). The frequency of TVV and M. hominis was 16.6% (n=5) and 20% (n=6), respectively.

Conclusion:

Presence of three clones among 30 T. vaginalis isolates indicated that microsatellite-based genotyping was efficient to determine the clonal distribution of T. vaginalis isolates. Therefore, a promising tool might be developed further and adapted to the studies dealing with molecular epidemiology of T. vaginalis. Microsatellite data from forthcoming studies will be deposited and presented on the website. In addition, we also presented the frequency of two endosymbionts in T. vaginalis isolates for the first time in Turkey.

Genomic Organization and Generation of Genetic Variability in the RHS (Retrotransposon Hot Spot) Protein Multigene Family in Trypanosoma cruzi

Retrotransposon Hot Spot (RHS) is the most abundant gene family in Trypanosoma cruzi, with unknown function in this parasite. The aim of this work was to shed light on the organization and expression of RHS in T. cruzi. The diversity of the RHS protein family in T. cruzi was demonstrated by phylogenetic and recombination analyses. Transcribed sequences carrying the RHS domain were classified into ten distinct groups of monophyletic origin. We identified numerous recombination events among the RHS and traced the origins of the donors and target sequences. The transcribed RHS genes have a mosaic structure that may contain fragments of different RHS inserted in the target sequence. About 30% of RHS sequences are located in the subtelomere, a region very susceptible to recombination. The evolution of the RHS family has been marked by many events, including gene duplication by unequal mitotic crossing-over, homologous, as well as ectopic recombination, and gene conversion. The expression of RHS was analyzed by immunofluorescence and immunoblotting using anti-RHS antibodies. RHS proteins are evenly distributed in the nuclear region of T. cruzi replicative forms (amastigote and epimastigote), suggesting that they could be involved in the control of the chromatin structure and gene expression, as has been proposed for T. brucei.

Differential Modulation of Mouse Heart Gene Expression by Infection With Two Trypanosoma cruzi Strains: A Transcriptome Analysis

The protozoan Trypanosoma cruzi (T. cruzi) is a well-adapted parasite to mammalian hosts and the pathogen of Chagas disease in humans. As both host and T. cruzi are highly genetically diverse, many variables come into play during infection, making disease outcomes difficult to predict. One important challenge in the field of Chagas disease research is determining the main factors leading to parasite establishment in the chronic stage in some organs, mainly the heart and/or digestive system. Our group previously showed that distinct strains of T. cruzi (JG and Col1.7G2) acquired differential tissue distribution in the chronic stage in dually infected BALB/c mice. To investigate changes in the host triggered by the two distinct T. cruzi strains, we assessed the gene expression profiles of BALB/c mouse hearts infected with either JG, Col1.7G2 or an equivalent mixture of both parasites during the initial phase of infection. This study demonstrates the clear differences in modulation of host gene expression by both parasites. Col1.7G2 strongly activated Th1-polarized immune signature genes, whereas JG caused only minor activation of the host immune response. Moreover, JG strongly reduced the expression of genes encoding ribosomal proteins and mitochondrial proteins related to the electron transport chain. Interestingly, the evaluation of gene expression in mice inoculated with a mixture of the parasites produced expression profiles with both up- and downregulated genes, indicating the coexistence of both parasite strains in the heart during the acute phase. This study suggests that different strains of T. cruzi may be distinguished by their efficiency in activating the immune system, modulating host energy metabolism and reactive oxygen species production and decreasing protein synthesis during early infection, which may be crucial for parasite persistence in specific organs.

Immigration and establishment of Trypanosoma cruzi in Arequipa, Peru

Changing environmental conditions, including those caused by human activities, reshape biological communities through both loss of native species and establishment of non-native species in the altered habitats. Dynamic interactions with the abiotic environment impact both immigration and initial establishment of non-native species into these altered habitats. The repeated emergence of disease systems in urban areas worldwide highlights the importance of understanding how dynamic migratory processes affect the current and future distribution and abundance of pathogens in urban environments. In this study, we examine the pattern of invasion of Trypanosoma cruzi-the causative agent of human Chagas disease-in the city of Arequipa, Peru. Phylogenetic analyses of 136 T. cruzi isolates from Arequipa and other South American locations suggest that only one T. cruzi lineage established a population in Arequipa as all T. cruzi isolated from vectors in Arequipa form a recent monophyletic group within the broader South American phylogeny. We discuss several hypotheses that may explain the limited number of established T. cruzi lineages despite multiple introductions of the parasite.

Elucidating diversity in the class composition of the minicircle hypervariable region of Trypanosoma cruzi: New perspectives on typing and kDNA inheritance

Background

Trypanosoma cruzi, the protozoan causative of Chagas disease, is classified into six main Discrete Typing Units (DTUs): TcI-TcVI. This parasite has around 10⁵ copies of the minicircle hypervariable region (mHVR) in their kinetoplastic DNA (kDNA). The genetic diversity of the mHVR is virtually unknown. However, cross-hybridization assays using mHVRs showed hybridization only between isolates belonging to the same genetic group. Nowadays there is no methodologic approach with a good sensibility, specificity and reproducibility for direct typing on biological samples. Due to its high copy number and apparently high diversity, mHVR becomes a good target for typing.

Methodology/Principal findings

Around 22 million reads, obtained by amplicon sequencing of the mHVR, were analyzed for nine strains belonging to six T. cruzi DTUs. The number and diversity of mHVR clusters was variable among DTUs and even within a DTU. However, strains of the same DTU shared more mHVR clusters than strains of different DTUs and clustered together. In addition, hybrid DTUs (TcV and TcVI) shared similar percentages (1.9–3.4%) of mHVR clusters with their parentals (TcII and TcIII). Conversely, just 0.2% of clusters were shared between TcII and TcIII suggesting biparental inheritance of the kDNA in hybrids. Sequencing at low depth (20,000–40,000 reads) also revealed 95% of the mHVR clusters for each of the analyzed strains. Finally, the method revealed good correlation in cluster identity and abundance between different replications of the experiment (r = 0.999).

Conclusions/Significance

Our work sheds light on the sequence diversity of mHVRs at intra and inter-DTU level. The mHVR amplicon sequencing workflow described here is a reproducible technique, that allows multiplexed analysis of hundreds of strains and results promissory for direct typing on biological samples in a future. In addition, such approach may help to gain knowledge on the mechanisms of the minicircle evolution and phylogenetic relationships among strains.

Chagas disease is an important public health problem in Latin America showing a wide diversity of clinical manifestations and epidemiological patterns. It is caused by the parasite Trypanosoma cruzi. This parasite is genetically diverse and classified into six main lineages. However, the relationship between intra-specific genetic diversity and clinical or epidemiological features is not clear, mainly because low sensitivity for direct typing on biological samples. For this reason, genetic markers with high copy number are required to achieve sensitivity. Here, we deep sequenced and analyzed a DNA region present in the large mitochondria of the parasite (named as mHVR, 10⁵ copies per parasite) from strains belonging to the six main lineages in order to analyze mHVR diversity and to evaluate its usefulness for typing. Despite the high sequence diversity, strains of the same lineage shared more sequences than strains of different lineages. Curiously, hybrid lineages shared mHVR sequences with both parents suggesting that mHVR (and DNA minicircles from the mitochondria) are inherited from both parentals. The mHVR amplicon sequencing workflow proposed here is reproducible and, potentially, it would be useful for typing hundreds of biological samples at time. It also provides a valuable approach to perform evolutionary and functional studies.

Sexual reproduction in a natural Trypanosoma cruzi population

Background

Sexual reproduction provides an evolutionary advantageous mechanism that combines favorable mutations that have arisen in separate lineages into the same individual. This advantage is especially pronounced in microparasites as allelic reassortment among individuals caused by sexual reproduction promotes allelic diversity at immune evasion genes within individuals which is often essential to evade host immune systems. Despite these advantages, many eukaryotic microparasites exhibit highly-clonal population structures suggesting that genetic exchange through sexual reproduction is rare. Evidence supporting clonality is particularly convincing in the causative agent of Chagas disease, Trypanosoma cruzi, despite equally convincing evidence of the capacity to engage in sexual reproduction.

Methodology/ Principle Findings

In the present study, we investigated two hypotheses that can reconcile the apparent contradiction between the observed clonal population structure and the capacity to engage in sexual reproduction by analyzing the genome sequences of 123 T. cruzi isolates from a natural population in Arequipa, Peru. The distribution of polymorphic markers within and among isolates provides clear evidence of the occurrence of sexual reproduction. Large genetic segments are rearranged among chromosomes due to crossing over during meiosis leading to a decay in the genetic linkage among polymorphic markers compared to the expectations from a purely asexually-reproducing population. Nevertheless, the population structure appears clonal due to a high level of inbreeding during sexual reproduction which increases homozygosity, and thus reduces diversity, within each inbreeding lineage.

Conclusions/ Significance

These results effectively reconcile the apparent contradiction by demonstrating that the clonal population structure is derived not from infrequent sex in natural populations but from high levels of inbreeding. We discuss epidemiological consequences of this reproductive strategy on genome evolution, population structure, and phenotypic diversity of this medically important parasite.

The rearrangement of alleles among individuals in a population during sexual reproduction maintains high allelic diversity within individuals in a population at polymorphic genes. Allelic diversity within individuals can be particularly important for parasites as it enhances their ability to evade host immune systems. Despite the potential benefits of sexual reproduction for parasites, natural populations of the protozoan parasite—and causative agent of human Chagas disease—Trypanosoma cruzi, exhibit clonal population structures indicative of asexual reproduction. This is particularly surprising as T. cruzi has the capacity for sexual reproduction. Here, we resolve this apparent contradiction by sequencing whole genomes of 123 T. cruzi isolates from a natural population in Arequipa, Peru. Evidence of past sexual reproduction and allelic rearrangements are common in this T. cruzi population. However, the majority of sexual reproduction events occur between close relatives resulting in an apparent clonal population structure. Sexual reproduction with distant relatives in areas with greater strain diversity has the potential to affect public health by increasing diversity in immune evasion genes within individuals and enhancing within-host survival, rapidly diversifying antigens that could affect the sensitivity of serological diagnostics, and by generating diversity in pathogenicity or drug resistance.

Molecular genotyping, diversity studies and high-resolution molecular markers unveiled by microsatellites in Giardia duodenalis

BACKGROUND

Giardia duodenalis (synonyms G. lamblia and G. intestinalis) is an enteric protozoan parasite of a wide range of mammalian hosts, including humans and various domestic and wild animals. There is considerable genetic variability in G. duodenalis and isolates of this parasite have been divided into eight genetic assemblages. Microsatellites markers can be used to discriminate isolates with a high level of sensitivity. This study was conducted to identify and characterize genomic microsatellites (simple sequence repeats-SSRs), sequences of one- to six-nucleotide motifs repeated in tandem, present in the available genomes of G. duodenalis and to develop new markers that can serve as a tool for detection and for characterizing the genetic diversity of this parasite.

METHODOLOGY/ PRINCIPAL FINDINGS

For each genetic assemblage, polymorphism levels for the microsatellite markers were evaluated. After performing the analysis using the MISA and SciRoKo software, 1,853 simple sequence repeats (SSRs) were identified. In all the genomes, trinucleotide repeats were the most common class followed by tetranucleotide. Many of the SSR loci are assemblage-specific, and 36 SSR loci shared among all the genomes were identified. Together with hypothetical proteins, variant-specific surface proteins represented nearly half of the annotated SSR loci. The results regarding the most common repeat among the SSRs led us to infer that positive selection occurred to avoid frameshift mutations. Additionally, based on inter- and intra-genetic assemblages polymorphism analyses, we unveiled previously undetected genetic variation, indicating that the microsatellite markers we developed are useful molecular tools for epidemiological inferences based on population genetics patterns and processes.

CONCLUSIONS

There is increasing demand for the development of new molecular markers and for the characterization of pathogens at a higher resolution level. In this study, we present 60 G. duodenalis microsatellites markers that exhibited high polymerase chain reaction (PCR) amplification efficiency among the different genetic assemblages. Twenty of these markers presented nucleotide sequence polymorphisms and may be used as a genotyping tool. The monomorphic markers can be used for the detection of the parasite at the species and genetic assemblage level. These polymorphic markers revealed a genetic diversity that was previously undetectable, thus they can be considered valuable molecular tools for high resolution markers in future studies investigating Giardia and may also be used for epidemiological inferences based on populations genetics patterns and processes.

Recent advances in trypanosomatid research: genome organization, expression, metabolism, taxonomy and evolution

Unicellular flagellates of the family Trypanosomatidae are obligatory parasites of invertebrates, vertebrates and plants. Dixenous species are aetiological agents of a number of diseases in humans, domestic animals and plants. Their monoxenous relatives are restricted to insects. Because of the high biological diversity, adaptability to dramatically different environmental conditions, and omnipresence, these protists have major impact on all biotic communities that still needs to be fully elucidated. In addition, as these organisms represent a highly divergent evolutionary lineage, they are strikingly different from the common 'model system' eukaryotes, such as some mammals, plants or fungi. A number of excellent reviews, published over the past decade, were dedicated to specialized topics from the areas of trypanosomatid molecular and cell biology, biochemistry, host-parasite relationships or other aspects of these fascinating organisms. However, there is a need for a more comprehensive review that summarizing recent advances in the studies of trypanosomatids in the last 30 years, a task, which we tried to accomplish with the current paper.

Cardiomyocyte oxidants production may signal to T. cruzi intracellular development

Chagas disease, caused by the protozoan parasite Trypanosoma cruzi, presents a variable clinical course, varying from asymptomatic to serious debilitating pathologies with cardiac, digestive or cardio-digestive impairment. Previous studies using two clonal T. cruzi populations, Col1.7G2 (T. cruzi I) and JG (T. cruzi II) demonstrated that there was a differential tissue distribution of these parasites during infection in BALB/c mice, with predominance of JG in the heart. To date little is known about the mechanisms that determine this tissue selection. Upon infection, host cells respond producing several factors, such as reactive oxygen species (ROS), cytokines, among others. Herein and in agreement with previous data from the literature we show that JG presents a higher intracellular multiplication rate when compared to Col1.7G2. We also showed that upon infection cardiomyocytes in culture may increase the production of oxidative species and its levels are higher in cultures infected with JG, which expresses lower levels of antioxidant enzymes. Interestingly, inhibition of oxidative stress severely interferes with the intracellular multiplication rate of JG. Additionally, upon H2O2-treatment increase in intracellular Ca2+ and oxidants were observed only in JG epimastigotes. Data presented herein suggests that JG and Col1.7G2 may sense extracellular oxidants in a distinct manner, which would then interfere differently with their intracellular development in cardiomyocytes.

Geographical, landscape and host associations of Trypanosoma cruzi DTUs and lineages

Background

The evolutionary history and ecological associations of Trypanosoma cruzi, the need to identify genetic markers that can distinguish parasite subpopulations, and understanding the parasite’s evolutionary and selective processes have been the subject of a significant number of publications since 1998, the year when the first DNA sequence analysis for the species was published.

Methods

The current analysis systematizes and re-analyzes this original research, focusing on critical methodological and analytical variables and results that have given rise to interpretations of putative patterns of genetic diversity and diversification of T. cruzi lineages, discrete typing units (DTUs), and populations, and their associations with hosts, vectors, and geographical distribution that have been interpreted as evidence for parasite subpopulation specificities.

Results

Few studies use hypothesis-driven or quantitative analysis for T. cruzi phylogeny (16/58 studies) or phylogeography (10/13). Among these, only one phylogenetic and five phylogeographic studies analyzed molecular markers directly from tissues (i.e. not from isolates). Analysis of T. cruzi DTU or lineage niche and its geographical projection demonstrate extensive sympatry among all clades across the continent and no significant niche differences among DTUs. DTU beta-diversity was high, indicating diverse host assemblages across regions, while host dissimilarity was principally due to host species turnover and to a much lesser degree to nestedness. DTU-host order specificities appear related to trophic or microenvironmental interactions.

Conclusions

More rigorous study designs and analyses will be required to discern evolutionary processes and the impact of landscape modification on population dynamics and risk for T. cruzi transmission to humans.

Electronic supplementary material

The online version of this article (doi:10.1186/s13071-016-1918-2) contains supplementary material, which is available to authorized users.

Mining SNPs in extracellular vesicular transcriptome of Trypanosoma cruzi: a step closer to early diagnosis of neglected Chagas disease

One of the newest and strongest members of intercellular communicators, the Extracellular vesicles (EVs) and their enclosed RNAs; Extracellular RNAs (exRNAs) have been acknowledged as putative biomarkers and therapeutic targets for various diseases. Although a very deep insight has not been possible into the physiology of these vesicles, they are believed to be involved in cell-to-cell communication and host-pathogen interactions. EVs might be significantly helpful in discovering biomarkers for possible target identification as well as prognostics, diagnostics and developing vaccines. In recent studies, highly bioactive EVs have drawn attention of parasitologists for being able to communicate between different cells and having likeliness of reflecting both source and target environments. Next-generation sequencing (NGS) has eased the way to have a deeper insight into these vesicles and their roles in various diseases. This article arises from bioinformatics-based analysis and predictive data mining of transcriptomic (RNA-Seq) data of EVs, derived from different life stages of Trypanosoma cruzi; a causing agent of neglected Chagas disease. Variants (Single Nucleotide Polymorphisms (SNPs)) were mined from Extracellular vesicular transcriptomic data and functionally analyzed using different bioinformatics based approaches. Functional analysis showed the association of these variants with various important factors like Trans-Sialidase (TS), Alpha Tubulin, P-Type H+-ATPase, etc. which, in turn, are associated with disease in different ways. Some of the ‘candidate SNPs’ were found to be stage-specific, which strengthens the probability of finding stage-specific biomarkers. These results may lead to a better understanding of Chagas disease, and improved knowledge may provide further development of the biomarkers for prognosis, diagnosis and drug development for treating Chagas disease.

Microsatellites for microbiologists

Microsatellites are repeating sequences of 2-6base pairs of DNA. Currently, they are used as molecular markers in many organisms, specifically in genetic studies analyzing kinship and population structure. In addition, they can be used to study gene duplication and/or deletion. Although they are used in studies on microorganisms including fungi, bacteria, protists, and archaea, it appears that these genetic markers are not being utilized to their full microbiological potential. Microsatellites have many advantages over other genetic markers currently in use as they are in general species specific, and therefore, cross-contamination by nontarget organisms is rare. Furthermore, microsatellites are suitable for use with fast and cheap DNA extraction methods, with ancient DNA or DNA from hair and fecal samples used in noninvasive sampling, making them widely available as a genetic marker. Microsatellites have already proven to be a useful tool for evolutionary studies of pathogenic microorganisms such as Candida albicans and Helicobacter pylori, and the onset of new sequencing techniques (such as 454, PACBIO, and mini-ion sequencing) means the ability to detect such markers will become less time consuming and cheaper, thus further expanding their potential to answer important microbial ecology questions.

Population Structure in Naegleria fowleri as Revealed by Microsatellite Markers

Naegleria sp. is a free living amoeba belonging to the Heterolobosea class. Over 40 species of Naegleria were identified and recovered worldwide in different habitats such as swimming pools, freshwater lakes, soil or dust. Among them, N. fowleri, is a human pathogen responsible for primary amoeboic meningoencephalitis (PAM). Around 300 cases were reported in 40 years worldwide but PAM is a fatal disease of the central nervous system with only 5% survival of infected patients. Since both pathogenic and non pathogenic species were encountered in the environment, detection and dispersal mode are crucial points in the fight against this pathogenic agent. Previous studies on identification and genotyping of N. fowleri strains were focused on RAPD analysis and on ITS sequencing and identified 5 variants: euro-american, south pacific, widespread, cattenom and chooz. Microsatellites are powerful markers in population genetics with broad spectrum of applications (such as paternity test, fingerprinting, genetic mapping or genetic structure analysis). They are characterized by a high degree of length polymorphism. The aim of this study was to genotype N. fowleri strains using microsatellites markers in order to track this population and to better understand its evolution. Six microsatellite loci and 47 strains from different geographical origins were used for this analysis. The microsatellite markers revealed a level of discrimination higher than any other marker used until now, enabling the identification of seven genetic groups, included in the five main genetic groups based on the previous RAPD and ITS analyses. This analysis also allowed us to go further in identifying private alleles highlighting intra-group variability. A better identification of the N. fowleri isolates could be done with this type of analysis and could allow a better tracking of the clinical and environmental N. fowleri strains.

Current drug therapy and pharmaceutical challenges for Chagas disease

One of the most significant health problems in the American continent in terms of human health, and socioeconomic impact is Chagas disease, caused by the protozoan parasite Trypanosoma cruzi. Infection was originally transmitted by reduviid insects, congenitally from mother to fetus, and by oral ingestion in sylvatic/rural environments, but blood transfusions, organ transplants, laboratory accidents, and sharing of contaminated syringes also contribute to modern day transmission. Likewise, Chagas disease used to be endemic from Northern Mexico to Argentina, but migrations have earned it global. The parasite has a complex life cycle, infecting different species, and invading a variety of cells - including muscle and nerve cells of the heart and gastrointestinal tract - in the mammalian host. Human infection outcome is a potentially fatal cardiomyopathy, and gastrointestinal tract lesions. In absence of a vaccine, vector control and treatment of patients are the only tools to control the disease. Unfortunately, the only drugs now available for Chagas' disease, Nifurtimox and Benznidazole, are relatively toxic for adult patients, and require prolonged administration. Benznidazole is the first choice for Chagas disease treatment due to its lower side effects than Nifurtimox. However, different strategies are being sought to overcome Benznidazole's toxicity including shorter or intermittent administration schedules-either alone or in combination with other drugs. In addition, a long list of compounds has shown trypanocidal activity, ranging from natural products to specially designed molecules, re-purposing drugs commercialized to treat other maladies, and homeopathy. In the present review, we will briefly summarize the upturns of current treatment of Chagas disease, discuss the increment on research and scientific publications about this topic, and give an overview of the state-of-the-art research aiming to produce an alternative medication to treat T. cruzi infection.

Characterizing the genetic diversity of the monkey malaria parasite Plasmodium cynomolgi

Plasmodium cynomolgi is a malaria parasite that typically infects Asian macaque monkeys, and humans on rare occasions. P. cynomolgi serves as a model system for the human malaria parasite Plasmodium vivax, with which it shares such important biological characteristics as formation of a dormant liver stage and a preference to invade reticulocytes. While genomes of three P. cynomolgi strains have been sequenced, genetic diversity of P. cynomolgi has not been widely investigated. To address this we developed the first panel of P. cynomolgi microsatellite markers to genotype eleven P. cynomolgi laboratory strains and 18 field isolates from Sarawak, Malaysian Borneo. We found diverse genotypes among most of the laboratory strains, though two nominally different strains were found to be genetically identical. We also investigated sequence polymorphism in two erythrocyte invasion gene families, the reticulocyte binding protein and Duffy binding protein genes, in these strains. We also observed copy number variation in rbp genes.

Identification and Characterization of Microsatellite Markers Derived from the Whole Genome Analysis of Taenia solium

Background

Infections with Taenia solium are the most common cause of adult acquired seizures worldwide, and are the leading cause of epilepsy in developing countries. A better understanding of the genetic diversity of T. solium will improve parasite diagnostics and transmission pathways in endemic areas thereby facilitating the design of future control measures and interventions. Microsatellite markers are useful genome features, which enable strain typing and identification in complex pathogen genomes. Here we describe microsatellite identification and characterization in T. solium, providing information that will assist in global efforts to control this important pathogen.

Methods

For genome sequencing, T. solium cysts and proglottids were collected from Huancayo and Puno in Peru, respectively. Using next generation sequencing (NGS) and de novo assembly, we assembled two draft genomes and one hybrid genome. Microsatellite sequences were identified and 36 of them were selected for further analysis. Twenty T. solium isolates were collected from Tumbes in the northern region, and twenty from Puno in the southern region of Peru. The size-polymorphism of the selected microsatellites was determined with multi-capillary electrophoresis. We analyzed the association between microsatellite polymorphism and the geographic origin of the samples.

Results

The predicted size of the hybrid (proglottid genome combined with cyst genome) T. solium genome was 111 MB with a GC content of 42.54%. A total of 7,979 contigs (>1,000 nt) were obtained. We identified 9,129 microsatellites in the Puno-proglottid genome and 9,936 in the Huancayo-cyst genome, with 5 or more repeats, ranging from mono- to hexa-nucleotide. Seven microsatellites were polymorphic and 29 were monomorphic within the analyzed isolates. T. solium tapeworms were classified into two genetic groups that correlated with the North/South geographic origin of the parasites.

Conclusions/Significance

The availability of draft genomes for T. solium represents a significant step towards the understanding the biology of the parasite. We report here a set of T. solium polymorphic microsatellite markers that appear promising for genetic epidemiology studies.

Taenia solium, the pork tapeworm, is an important pathogen as it is a major cause of acquired epilepsy in developing countries. The parasite was eliminated from most developed countries decades ago due to improvement in sanitary conditions but it remains a common infection across Asia, Africa and Latin America. Identification of genetic variants within T. solium will enable to study the genetic epidemiology, distribution and movement of this parasite within endemic communities, which will ultimately facilitate the design of control strategies to reduce the health and economic burden of disease.

Microsatellites have been used in other parasites to identify genetic variants. In this study, we partially sequenced the genome of T. solium and identified microsatellites widely distributed in the genome using bioinformatics tools. We evaluated the distribution of these microsatellites collected from 20 tapeworms from the north and 20 tapeworms from the south of Peru. We identified seven polymorphic microsatellites, and evaluated their capacity to differentiate genetic variants of T. solium. Interestingly, tapeworms from the North and South of Peru showed different genotypes, suggesting its use as a potential marker to differentiate geographic origin.

The population genetics of Trypanosoma cruzi revisited in the light of the predominant clonal evolution model

Comparing the population structure of Trypanosoma cruzi with that of other pathogens, including parasitic protozoa, fungi, bacteria and viruses, shows that the agent of Chagas disease shares typical traits with many other species, related to a predominant clonal evolution (PCE) pattern: statistically significant linkage disequilibrium, overrepresented multilocus genotypes, near-clades (genetic subdivisions somewhat blurred by occasional genetic exchange/hybridization) and "Russian doll" patterns (PCE is observed, not only at the level of the whole species, but also, within the near-clades). Moreover, T. cruzi population structure exhibits linkage with the diversity of several strongly selected genes, with gene expression profiles, and with some major phenotypic traits. We discuss the evolutionary significance of these results, and their implications in terms of applied research (molecular epidemiology/strain typing, analysis of genes of interest, vaccine and drug design, immunological diagnosis) and of experimental evolution. Lastly, we revisit the long-term debate of describing new species within the T. cruzi taxon.

Assessment of population genetic structure in the arbovirus vector midge, Culicoides brevitarsis (Diptera: Ceratopogonidae), using multi-locus DNA microsatellites

Bluetongue virus (BTV) is a major pathogen of ruminants that is transmitted by biting midges (Culicoides spp.). Australian BTV serotypes have origins in Asia and are distributed across the continent into two distinct episystems, one in the north and another in the east. Culicoides brevitarsis is the major vector of BTV in Australia and is distributed across the entire geographic range of the virus. Here, we describe the isolation and use of DNA microsatellites and gauge their ability to determine population genetic connectivity of C. brevitarsis within Australia and with countries to the north. Eleven DNA microsatellite markers were isolated using a novel genomic enrichment method and identified as useful for genetic analyses of sampled populations in Australia, northern Papua New Guinea (PNG) and Timor-Leste. Significant (P < 0.05) population genetic subdivision was observed between all paired regions, though the highest levels of genetic sub-division involved pair-wise tests with PNG (PNG vs. Australia (FST = 0.120) and PNG vs. Timor-Leste (FST = 0.095)). Analysis of multi-locus allelic distributions using STRUCTURE identified a most probable two-cluster population model, which separated PNG specimens from a cluster containing specimens from Timor-Leste and Australia. The source of incursions of this species in Australia is more likely to be Timor-Leste than PNG. Future incursions of BTV positive C. brevitarsis into Australia may be genetically identified to their source populations using these microsatellite loci. The vector's panmictic genetic structure within Australia cannot explain the differential geographic distribution of BTV serotypes.

Assessment of population genetic structure in the arbovirus vector midge, Culicoides brevitarsis (Diptera: Ceratopogonidae), using multi-locus DNA microsatellites

Bluetongue virus (BTV) is a major pathogen of ruminants that is transmitted by biting midges (Culicoides spp.). Australian BTV serotypes have origins in Asia and are distributed across the continent into two distinct episystems, one in the north and another in the east. Culicoides brevitarsis is the major vector of BTV in Australia and is distributed across the entire geographic range of the virus. Here, we describe the isolation and use of DNA microsatellites and gauge their ability to determine population genetic connectivity of C. brevitarsis within Australia and with countries to the north. Eleven DNA microsatellite markers were isolated using a novel genomic enrichment method and identified as useful for genetic analyses of sampled populations in Australia, northern Papua New Guinea (PNG) and Timor-Leste. Significant (P < 0.05) population genetic subdivision was observed between all paired regions, though the highest levels of genetic sub-division involved pair-wise tests with PNG (PNG vs. Australia (F_ST = 0.120) and PNG vs. Timor-Leste (F_ST = 0.095)). Analysis of multi-locus allelic distributions using STRUCTURE identified a most probable two-cluster population model, which separated PNG specimens from a cluster containing specimens from Timor-Leste and Australia. The source of incursions of this species in Australia is more likely to be Timor-Leste than PNG. Future incursions of BTV positive C. brevitarsis into Australia may be genetically identified to their source populations using these microsatellite loci. The vector’s panmictic genetic structure within Australia cannot explain the differential geographic distribution of BTV serotypes.

Between a bug and a hard place: Trypanosoma cruzi genetic diversity and the clinical outcomes of Chagas disease

Over the last 30 years, concomitant with successful transnational disease control programs across Latin America, Chagas disease has expanded from a neglected, endemic parasitic infection of the rural poor to an urbanized chronic disease, and now a potentially emergent global health problem. Trypanosoma cruzi infection has a highly variable clinical course, ranging from complete absence of symptoms to severe and often fatal cardiovascular and/or gastrointestinal manifestations. To date, few correlates of clinical disease progression have been identified. Elucidating a putative role for T. cruzi strain diversity in Chagas disease pathogenesis is complicated by the scarcity of parasites in clinical specimens and the limitations of our contemporary genotyping techniques. This article systematically reviews the historical literature, given our current understanding of parasite genetic diversity, to evaluate the evidence for any association between T. cruzi genotype and chronic clinical outcome, risk of congenital transmission or reactivation and orally transmitted outbreaks.

Trypanosoma rangeli displays a clonal population structure, revealing a subdivision of KP1(-) strains and the ancestry of the Amazonian group

Assessment of the genetic variability and population structure of Trypanosoma rangeli, a non-pathogenic American trypanosome, was carried out through microsatellite and single-nucleotide polymorphism analyses. Two approaches were used for microsatellite typing: data mining in expressed sequence tag /open reading frame expressed sequence tags libraries and PCR-based Isolation of Microsatellite Arrays from genomic libraries. All microsatellites found were evaluated for their abundance, frequency and usefulness as markers. Genotyping of T. rangeli strains and clones was performed for 18 loci amplified by PCR from expressed sequence tag/open reading frame expressed sequence tags libraries. The presence of single-nucleotide polymorphisms in the nuclear, multi-copy, spliced leader gene was assessed in 18 T. rangeli strains, and the results show that T. rangeli has a predominantly clonal population structure, allowing a robust phylogenetic analysis. Microsatellite typing revealed a subdivision of the KP1(-) genetic group, which may be influenced by geographical location and/or by the co-evolution of parasite and vectors occurring within the same geographical areas. The hypothesis of parasite-vector co-evolution was corroborated by single-nucleotide polymorphism analysis of the spliced leader gene. Taken together, the results suggest three T. rangeli groups: (i) the T. rangeli Amazonian group; (ii) the T. rangeli KP1(-) group; and (iii) the T. rangeli KP1(+) group. The latter two groups possibly evolved from the Amazonian group to produce KP1(+) and KP1(-) strains.

Molecular genotyping of Trypanosoma cruzi for lineage assignment and population genetics

Trypanosoma cruzi, the etiological agent of Chagas disease, remains a major public health problem in Latin America. Infection with T. cruzi is lifelong and can lead to a spectrum of pathological sequelae ranging from subclinical to lethal cardiac and/or gastrointestinal complications. Isolates of T. cruzi can be assigned to six genetic lineages or discrete typing units (DTUs), which are broadly associated with disparate ecologies, transmission cycles, and geographical distributions. This extensive genetic diversity is also believed to contribute to the clinical variation observed among chagasic patients. Unravelling the population structure of T. cruzi is fundamental to understanding Chagas disease epidemiology, developing control strategies, and resolving the relationship between parasite genotype and clinical prognosis. To date, no single, widely validated, genetic target allows unequivocal resolution to DTU-level. In this chapter we present standardized methods for strain DTU assignment using PCR-restriction fragment length polymorphism analysis (PCR-RFLP) and nuclear multilocus sequence typing (MLST). PCR-RFLPs have the advantages of simplicity and reproducibility, requiring limited expertise and few laboratory consumables. MLST data are more laborious to generate but more informative; DNA sequences are readily transferable between research groups and amenable to recombination detection and intra-lineage analyses. We also recommend a mitochondrial (maxicircle) MLST scheme and a panel of 28 microsatellite loci for higher resolution population genetics studies. Due to the scarcity of T. cruzi in blood and tissue, all of these genotyping techniques have limited sensitivity when applied directly to clinical or biological specimens, particularly when targets are single (MLST) or low copy number (PCR-RFLPs). We therefore describe essential protocols to isolate parasites, derive biological clones, and extract T. cruzi genomic DNA from field and clinical samples.

Species-specific markers for the differential diagnosis of Trypanosoma cruzi and Trypanosoma rangeli and polymorphisms detection in Trypanosoma rangeli

Trypanosoma cruzi and Trypanosoma rangeli are kinetoplastid parasites which are able to infect humans in Central and South America. Misdiagnosis between these trypanosomes can be avoided by targeting barcoding sequences or genes of each organism. This work aims to analyze the feasibility of using species-specific markers for identification of intraspecific polymorphisms and as target for diagnostic methods by PCR. Accordingly, primers which are able to specifically detect T. cruzi or T. rangeli genomic DNA were characterized. The use of intergenic regions, generally divergent in the trypanosomatids, and the serine carboxypeptidase gene were successful. Using T. rangeli genomic sequences for the identification of group-specific polymorphisms and a polymorphic AT(n) dinucleotide repeat permitted the classification of the strains into two groups, which are entirely coincident with T. rangeli main lineages, KP1 (+) and KP1 (-), previously determined by kinetoplast DNA (kDNA) characterization. The sequences analyzed totalize 622 bp (382 bp represent a hypothetical protein sequence, and 240 bp represent an anonymous sequence), and of these, 581 (93.3%) are conserved sites and 41 bp (6.7%) are polymorphic, with 9 transitions (21.9%), 2 transversions (4.9%), and 30 (73.2%) insertion/deletion events. Taken together, the species-specific markers analyzed may be useful for the development of new strategies for the accurate diagnosis of infections. Furthermore, the identification of T. rangeli polymorphisms has a direct impact in the understanding of the population structure of this parasite.

Genomic and proteomic approaches for Chagas’ disease: critical analysis of diagnostic methods

Trypanosoma cruzi is the etiologic agent of Chagas' disease, a chronic inflammatory condition that results in heart and digestive complications. The first draft of the parasite genome is now complete and it is expected that, along with the published genomic and proteomic analyses discussed herein, it will lead to the identification of crucial genes and proteins directly associated with disease. This article reviews the current research trends addressing host-parasite interaction, parasite genetic variability and diagnosis. These advances will certainly bring about major developments not only in our understanding of Trypanosoma cruzi biology, but also in the application of new technologies to disease prevention and control.

A call to arms: on refining Plasmodium vivax microsatellite marker panels for comparing global diversity

BACKGROUND

Microsatellite (MS) markers have become an important tool for studying the population diversity, evolutionary history and multiplicity of infection (MOI) of malaria parasite infections. MS are typically selected on the basis of being highly polymorphic. However, it is known that the polymorphic potential (mutability) of each marker can vary as much as two orders of magnitude, which radically changes how diversity is represented in the genome from one marker to the next. Over the past decade, approximately 240 Plasmodium vivax MS have been published, comprising nine major panels of markers. Inconsistent usage of each panel has resulted in a surfeit of descriptive genetic diversity data that are largely incomparable between populations. The objective of this study was to statistically evaluate the quality of individual MS markers in order to validate a refined panel of markers that will provide a balanced picture of P. vivax population diversity.

METHODS

All previously published data, including genetic diversity indices, MS parameters, and population parameters, were assembled from 18 different global studies into a flat file to facilitate statistical analysis and modelling using JMP® Genomics 6.0 (SAS Institute Inc, Cary, NC, USA). Statistical modeling was employed to down-select markers with extreme variation among the mean number of alleles, expected heterozygosity, maximum repeat length and/or chromosomal location of the repeat. Individual MS were analysed by step-down whole model linear regression and standard least squares fit models, both stratified by annual parasite incidence to identify MS markers with values significantly different from the mean.

RESULTS

Of the 42 MS under evaluation in this study, 18 (nine high priority) were identified as ideal candidates for measuring population diversity between global regions, while five (two high priority) additional markers were identified as candidates for MOI studies.

CONCLUSIONS

MS diversity was found to be a function of endemicity and motif structure. Evaluation of individual MS permitted the assembly of a refined panel of markers that can be reliably utilized in the field to compare population structures between global regions.

Evidence of substantial recombination among Trypanosoma cruzi II strains from Minas Gerais

Due to the scarcity of evidence of sexuality in Trypanosoma cruzi, the causative agent of Chagas disease, it has been general accepted that the parasite reproduction is essentially clonal with infrequent genetic recombination. This assumption is mainly supported by indirect evidence, such as Hardy-Weinberg imbalances, linkage disequilibrium and a strong correlation between independent sets of genetic markers of T. cruzi populations. However, because the analyzed populations are usually isolated from different geographic regions, the possibility of population substructuring as generating these genetic marker imbalances cannot be eliminated. To investigate this possibility, we firstly compared the allele frequencies and haplotype networks using seven different polymorphic loci (two from mitochondrial and five from different nuclear chromosomes) in two groups of TcII strains: one including isolates obtained from different regions in Latin America and the other including isolates obtained only from patients of the Minas Gerais State in Brazil. Our hypothesis was that if the population structure is essentially clonal, Hardy-Weinberg disequilibrium and a sharp association between the clusters generated by analyzing independent markers should be observed in both strain groups, independent of the geographic origin of the samples. The results demonstrated that the number of microsatellite loci in linkage disequilibrium decreased from 4 to 1 when only strains from Minas Gerais were analyzed. Moreover, we did not observed any correlation between the clusters when analyzing the nuclear and mitochondrial loci, suggesting independent inheritance of these markers among the Minas Gerais strains. Besides, using a second subset of five physically linked microsatellite loci and the Minas Gerais strains, we could also demonstrate evidence of homologous recombination roughly proportional to the relative distance among them. Taken together, our results do not support a clonal population structure for T. cruzi, particularly in TcII, which coexists in the same geographical area, suggesting that genetic exchanges among these strains may occur more frequently than initially expected.

Putative panmixia in restricted populations of Trypanosoma cruzi isolated from wild Triatoma infestans in Bolivia

Trypanosoma cruzi, the causative agent of Chagas disease, is subdivided into six discrete typing units (DTUs; TcI-TcVI) of which TcI is ubiquitous and genetically highly variable. While clonality is the dominant mode of propagation, recombinant events play a significant evolutive role. Recently, foci of wild Triatoma infestans have been described in Bolivia, mainly infected by TcI. Hence, for the first time, we evaluated the level of genetic exchange within TcI natural potentially panmictic populations (single DTU, host, area and sampling time). Seventy-nine TcI stocks from wild T. infestans, belonging to six populations were characterized at eight microsatellite loci. For each population, Hardy-Weinberg equilibrium (HWE), linkage disequilibrium (LD), and presence of repeated multilocus genotypes (MLG) were analyzed by using a total of seven statistics, to test the null hypothesis of panmixia (H0). For three populations, none of the seven statistics allowed to rejecting H0; for another one the low size did not allow us to conclude, and for the two others the tests have given contradictory results. Interestingly, apparent panmixia was only observed in very restricted areas, and was not observed when grouping populations distant of only two kilometers or more. Nevertheless it is worth stressing that for the statistic tests of "HWE", in order to minimize the type I error (i. e. incorrect rejection of a true H0), we used the Bonferroni correction (BC) known to considerably increase the type II error ( i. e. failure to reject a false H0). For the other tests (LD and MLG), we did not use BC and the risk of type II error in these cases was acceptable. Thus, these results should be considered as a good indicator of the existence of panmixia in wild environment but this must be confirmed on larger samples to reduce the risk of type II error.

Microsatellite loci-based distribution of Trypanosoma cruzi genotypes from Chilean chronic Chagas disease patients and Triatoma infestans is concordant with a specific host-parasite association hypothesis

The objective of this study was to investigate if there is specific host-parasite association in Chilean populations of Trypanosoma cruzi. For this purpose, two groups of parasites were analyzed, one from chronic chagasic patients, and the other from Triatoma infestans triatomines in three regions of the country. The first group consisted of four types of samples: parasites from peripheral blood of non-cardiopathic T. cruzi infected patients (NB); parasites from their corresponding xenodiagnosis (NX); parasites from peripheral blood of T. cruzi infected cardiopathic patients (CB) and parasites from their xenodiagnostics (CX). The T. infestans sample in turn was from three regions: III, V and M (Metropolitan). The genetic differentiation by the Fisher exact method, the lineage distribution of the samples, the molecular phylogeny and the frequency of multiclonality were analysed. The results show that not only are the groups of T. cruzi clones from Chagas disease patients and vectors genetically differentiated, but also all the sub-groups (NB, NX, CB and CX) from the III, V and M regions. The analysis of lineage distribution was concordant with the above results, because significant differences among the percentages of TcI, TcIII and hybrids (TcV or TcVI) were observed. The phylogenetic reconstruction with these Chilean T. cruzi samples was coherent with the above results because the four chagasic samples clustered together in a node with high bootstrap support, whereas the three triatomine samples (III, V and M) were located apart from that node. The topology of the tree including published T. cruzi clones and isolates was concordant with the known topology, which confirmed that the results presented here are correct and are not biased by experimental error. Taken together the results presented here are concordant with a specific host-parasite association between some Chilean T. cruzi populations.

Interclonal variations in the molecular karyotype of Trypanosoma cruzi: chromosome rearrangements in a single cell-derived clone of the G strain

Trypanosoma cruzi comprises a pool of populations which are genetically diverse in terms of DNA content, growth and infectivity. Inter- and intra-strain karyotype heterogeneities have been reported, suggesting that chromosomal rearrangements occurred during the evolution of this parasite. Clone D11 is a single-cell-derived clone of the T. cruzi G strain selected by the minimal dilution method and by infecting Vero cells with metacyclic trypomastigotes. Here we report that the karyotype of clone D11 differs from that of the G strain in both number and size of chromosomal bands. Large chromosomal rearrangement was observed in the chromosomes carrying the tubulin loci. However, most of the chromosome length polymorphisms were of small amplitude, and the absence of one band in clone D11 in relation to its reference position in the G strain could be correlated to the presence of a novel band migrating above or below this position. Despite the presence of chromosomal polymorphism, large syntenic groups were conserved between the isolates. The appearance of new chromosomal bands in clone D11 could be explained by chromosome fusion followed by a chromosome break or interchromosomal exchange of large DNA segments. Our results also suggest that telomeric regions are involved in this process. The variant represented by clone D11 could have been induced by the stress of the cloning procedure or could, as has been suggested for Leishmania infantum, have emerged from a multiclonal, mosaic parasite population submitted to frequent DNA amplification/deletion events, leading to a 'mosaic' structure with different individuals having differently sized versions of the same chromosomes. If this is the case, the variant represented by clone D11 would be better adapted to survive the stress induced by cloning, which includes intracellular development in the mammalian cell. Karyotype polymorphism could be part of the T. cruzi arsenal for responding to environmental pressure.

LSSP-PCR of Trypanosoma cruzi: how the single primer sequence affects the kDNA signature

Background

Low-stringency single specific primer PCR (LSSP-PCR) is a highly sensitive and discriminating technique that has been extensively used to genetically characterize Trypanosoma cruzi populations in the presence of large amounts of host DNA. To ensure high sensitivity, in most T. cruzi studies, the variable regions of the naturally amplified kinetoplast DNA (kDNA) minicircles were targeted, and this method translated the intraspecific polymorphisms of these molecules into specific and reproducible kDNA signatures. Although the LSSP-PCR technique is reproducible under strict assay conditions, the complex banding pattern generated can be significantly altered by even a single-base change in the target DNA. Our survey of the literature identified eight different primers with similar, if not identical, names that have been used for kDNA amplification and LSSP-PCR of T. cruzi. Although different primer sequences were used in these studies, many of the authors cited the same reference report to justify their primer choice. We wondered whether these changes in the primer sequence could affect also the parasite LSSP-PCR profiles.

Findings

To answer this question we compared the kDNA signatures obtained from three different and extensively studied T. cruzi populations with the eight primers found in the literature. Our results clearly demonstrate that even minimal modifications in the oligonucleotide sequences, especially in the 3′ or 5′ end, can significantly change the kDNA signature of a T. cruzi strain.

Conclusions

These results highlight the necessity of careful preservation of primer nomenclature and sequence when reproducing an LSSP-PCR work to avoid confusion and allow comparison of results among different laboratories.

Microsatellite and mini-exon analysis of Mexican human DTU I Trypanosoma cruzi strains and their susceptibility to nifurtimox and benznidazole

Chagas disease is caused by the protozoan parasite Trypanosoma cruzi, and it affects as many as 10 million people in North and South America, where it represents a major public health problem. T. cruzi is a parasite with high genetic diversity, and it has been grouped into 6 discrete typing units (DTUs), designated as T. cruzi I (TcI) to T. cruzi VI (TcVI). Mexican isolates from humans and from vector insects have been primarily found to be TcI, and these isolates are likely to be the strains that cause the clinical manifestations observed in Mexico. However, genetic characterization and drug susceptibility assays are limited in Mexican TcI strains. In this work, 24 Mexican T. cruzi strains, obtained primarily from humans, were studied with 7 locus microsatellites and mini-exon gene by PCR. Also, drug susceptibility was evaluated by growth and mobility assays. All of the human strains belonged to TcI, and they could be further grouped through microsatellite analysis into 2 subgroups (microsatellite genotypes 1 and 2), which were not related to the host clinical status or biological origin of the strain. Two strains, both from wild mammals, belonged to the TcII-TcVI groups; these strains and the CL Brener strain constituted microsatellite genotype 3. The number of alleles in each locus was lower than reported for South American strains, and a departure from the Hardy-Weinberg equilibrium was observed. The susceptibility of these strains to nifurtimox and benznidazole was heterogeneous. T. cruzi strains characterized as microsatellite genotypes 2 and 3 were significantly more susceptible to benznidazole than strains of microsatellite genotype 1. Only 1 Mexican strain resistant to both drugs was found in this study.

Unequivocal identification of subpopulations in putative multiclonal Trypanosoma cruzi strains by FACs single cell sorting and genotyping

Trypanosoma cruzi, the etiological agent of Chagas disease, is a polymorphic species. Evidence suggests that the majority of the T. cruzi populations isolated from afflicted humans, reservoir animals, or vectors are multiclonal. However, the extent and the complexity of multiclonality remain to be established, since aneuploidy cannot be excluded and current conventional cloning methods cannot identify all the representative clones in an infection. To answer this question, we adapted a methodology originally described for analyzing single spermatozoids, to isolate and study single T. cruzi parasites. Accordingly, the cloning apparatus of a Fluorescence-Activated Cell Sorter (FACS) was used to sort single T. cruzi cells directly into 96-wells microplates. Cells were then genotyped using two polymorphic genomic markers and four microsatellite loci. We validated this methodology by testing four T. cruzi populations: one control artificial mixture composed of two monoclonal populations--Silvio X10 cl1 (TcI) and Esmeraldo cl3 (TcII)--and three naturally occurring strains, one isolated from a vector (A316A R7) and two others derived from the first reported human case of Chagas disease. Using this innovative approach, we were able to successfully describe the whole complexity of these natural strains, revealing their multiclonal status. In addition, our results demonstrate that these T. cruzi populations are formed of more clones than originally expected. The method also permitted estimating of the proportion of each subpopulation of the tested strains. The single-cell genotyping approach allowed analysis of intrapopulation diversity at a level of detail not achieved previously, and may thus improve our comprehension of population structure and dynamics of T. cruzi. Finally, this methodology is capable to settle once and for all controversies on the issue of multiclonality.

Genetic modulation in Be-78 and Y Trypanosoma cruzi strains after long-term infection in Beagle dogs revealed by molecular markers

The genetic profile of Trypanosoma cruzi was evaluated in parasite populations isolated from Beagle dogs experimentally infected with Be-78 and Y strains that present distinct biological and genetic characteristics. Molecular characterization of the isolates obtained 30days and 2years after infection was carried out. For typing MLEE, sequence polymorphisms of the mitochondrial cytochrome oxidase subunit II gene (COII) and RAPD profiles were used. The profiles of MLEE were the same for the parental Be-78 strains as their respective isolates. However, changes of MLEE profile were observed in two T. cruzi isolates from dogs inoculated with Y strain. Changes in the mitochondrial DNA (COII) and RAPD profiles of the Y strain were also observed. The dendogram constructed by UPGMA with RAPD results indicated two major branches. Global data show that the genetic modulation in polyclonal strains during the long-term infection occurred and was strain-dependent. This study still suggests that each host (here each dog) harbors a determinate T. cruzi population that may change or be modulated throughout long-term infection. This might to hinder the observation of correlation between the genetics of T. cruzi and their biological properties and behavior in different host species due to the complexity of the parasite-host interaction in which probably the genetic background of both should be considered.

The revised Trypanosoma cruzi subspecific nomenclature: rationale, epidemiological relevance and research applications

The protozoan Trypanosoma cruzi, its mammalian reservoirs, and vectors have existed in nature for millions of years. The human infection, named Chagas disease, is a major public health problem for Latin America. T. cruzi is genetically highly diverse and the understanding of the population structure of this parasite is critical because of the links to transmission cycles and disease. At present, T. cruzi is partitioned into six discrete typing units (DTUs), TcI-TcVI. Here we focus on the current status of taxonomy-related areas such as population structure, phylogeographical and eco-epidemiological features, and the correlation of DTU with natural and experimental infection. We also summarize methods for DTU genotyping, available for widespread use in endemic areas. For the immediate future multilocus sequence typing is likely to be the gold standard for population studies. We conclude that greater advances in our knowledge on pathogenic and epidemiological features of these parasites are expected in the coming decade through the comparative analysis of the genomes from isolates of various DTUs.

Geographical structuring of Trypanosoma cruzi populations from Chilean Triatoma infestans triatomines and their genetic relationship with other Latino American counterparts

In order to obtain more information about the population structure of Chilean Trypanosoma cruzi, and their genetic relationship with other Latino American counterparts, we performed the study of T. cruzi samples detected in the midgut content of Triatoma infestans insects from three endemic regions of Chile. The genetic characteristics of these samples were analysed using microsatellite markers and PCR conditions that allow the detection of predominant T. cruzi clones directly in triatomine midgut content. Population genetic analyses using the Fisher's exact method, analysis of molecular variance (AMOVA) and the determination of F(ST) showed that the northern T. cruzi population sample was genetically differentiated from the two southern population counterparts. Further analysis showed that the cause of this genetic differentiation was the asymmetrical distribution of TcIII T. cruzi predominant clones. Considering all triatomines from the three regions, the most frequent predominant lineages were TcIII (38%), followed by TcI (34%) and hybrid (8%). No TcII lineage was observed along the predominant T. cruzi clones. The best phylogenetic reconstruction using the shared allelic genetic distance was concordant with the population genetic analysis and tree topology previously described studying foreign samples. The correlation studies showed that the lineage TcIII from the III region was genetically differentiated from the other two, and this differentiation was correlated with geographical distance including Chilean and mainly Brazilian samples. It will be interesting to investigate whether this geographical structure may be related with different clinical manifestation of Chagas disease.

Sexual recombination is a signature of a persisting malaria epidemic in Peru

Background

The aim of this study was to consider the impact that multi-clone, complex infections have on a parasite population structure in a low transmission setting. In general, complexity of infection (minimum number of clones within an infection) and the overall population level diversity is expected to be minimal in low transmission settings. Additionally, the parasite population structure is predicted to be clonal, rather than sexual due to infrequent parasite inoculation and lack of recombination between genetically distinct clones. However, in this low transmission of the Peruvian Amazon, complex infections are becoming more frequent, in spite of decreasing infection prevalence. In this study, it was hypothesized that sexual recombination between distinct clonal lineages of Plasmodium falciparum parasites were altering the subpopulation structure and effectively maintaining the population-level diversity.

Methods

Fourteen microsatellite markers were chosen to describe the genetic diversity in 313 naturally occurring P. falciparum infections from Peruvian Amazon. The population and subpopulation structure was characterized by measuring: clusteredness, expected heterozygosity (H_e), allelic richness, private allelic richness, and linkage disequilibrium. Next, microsatellite haplotypes and alleles were correlated with P. falciparum merozoite surface protein 1 Block 2 (Pfmsp1-B2) to examine the presence of recombinant microsatellite haplotypes.

Results

The parasite population structure consists of six genetically diverse subpopulations of clones, called "clusters". Clusters 1, 3, 4, and 6 have unique haplotypes that exceed 70% of the total number of clones within each cluster, while Clusters 2 and 5 have a lower proportion of unique haplotypes, but still exceed 46%. By measuring the H_e, allelic richness, and private allelic richness within each of the six subpopulations, relatively low levels of genetic diversity within each subpopulation (except Cluster 4) are observed. This indicated that the number of alleles, and not the combination of alleles, are limited. Next, the standard index of association (I_A^S) was measured, which revealed a significant decay in linkage disequilibrium (LD) associated with Cluster 6, which is indicative of independent assortment of alleles. This decay in LD is a signature of this subpopulation approaching linkage equilibrium by undergoing sexual recombination. To trace possible recombination events, the two most frequent microsatellite haplotypes observed over time (defined by either a K1 or Mad20) were selected as the progenitors and then potential recombinants were identified in within the natural population.

Conclusions

Contrary to conventional low transmission models, this study provides evidence of a parasite population structure that is superficially defined by a clonal backbone. Sexual recombination does occur and even arguably is responsible for maintaining the substructure of this population.