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.

Hydrophobicity-Driven Increases in Editing in Mitochondrial mRNAs during the Evolution of Kinetoplastids.

Kinetoplastids are a diverse group of flagellates which exhibit editing by insertion/deletion of Us in the mitochondrial mRNAs. Some mRNAs require editing to build most of their coding sequences, a process known as pan-editing. Evidence suggests that pan-editing is an ancestral feature in kinetoplastids. Here, we investigate how the transition from nonedited to pan-edited states occurred. The mitochondrial mRNAs and protein sequences from nine kinetoplastids and related groups (diplonemids, euglenids, and jakobids) were analyzed. RNA editing increased protein hydrophobicity to extreme values by introducing Us in the second codon position, despite the absence of editing preferences related to codon position. In addition, hydrophobicity was maintained by purifying selection in species that lost editing by retroposition of the fully edited mRNA. Only a few hydrophobic to hydrophilic amino acid changes were inferred for such species. In the protein secondary structure, these changes occurred spatially close to other hydrophilic residues. The analysis of coevolving sites showed that multiple changes are required together for hydrophobicity to be lost, which suggest the proteins are locked into extended hydrophobicity. Finally, an analysis of the NAD7 protein-protein interactions showed they can also influence hydrophobicity increase in the protein and where editing can occur in the mRNA. In conclusion, our results suggest that protein hydrophobicity has influenced editing site selection and how editing expanded in mRNAs. In effect, the hydrophobicity increase was entrenched by a neutral ratchet moved by a mutational pressure to introduce Us, thus helping to explain both RNA editing increase and, possibly, persistence.

Genome data vs MLST for exploring intraspecific evolutionary history in bacteria: Much is not always better.

Genome-based phylogeny has been proposed to be more accurate than phylogeny based in a few genes as MLST-based phylogeny. However, much is not always better. Here we analyzed 368 complete genomes corresponding to 9 bacterial species in order to address intraspecific phylogeny. The studied species were: Burkholderia pseudomallei, Campylobacter jejuni, Chlamydia trachomatis, Helicobacter pylori, Klebsiella pneumoniae, Listeria monocytogenes, Salmonella enterica, Staphylococcus aureus and Streptococcus pyogenes. The intra-specific phylogenies were inferred using the complete genome sequences of different strains of these species and their MLST schemes. A supermatrix approach was used to infer maximum likelihood phylogenies in both cases. The phylogenetic incongruence between the supermatrix-based genome or MLST tree and individual trees (constructed from genome fragments or MLST genes, respectively) was analyzed. In supermatrix-based trees for genomes, most branches showed a high branch support; however, a high number of branches also showed high percentage of topologically incongruent individual trees. Interestingly, genome and MLST trees showed similar levels of incongruence in the phylogeny for each bacteria specie. Both genome and MLST approaches showed that C. trachomatis and S. aureus have a tree-like evolutionary history (low levels of internal incongruence). Instead, B. pseudomallei and S. pyogenes show high levels of incongruence (network-like evolutionary story) probably caused by HGT (horizontal gene transfer). Concluding, our analysis showed that: high branch supports obtained in genome phylogenies could be an artifact probably caused by data size; MLST is valid to address intraspecific phylogenetic structure; and, each species has its own evolutionary history, which could be affected by HGT to different extents.

Guide RNA Repertoires in the Main Lineages of Trypanosoma cruzi: High Diversity and Variable Redundancy Among Strains.

Trypanosoma cruzi, as other kinetoplastids, has a complex mechanism of editing of mitochondrial mRNAs that requires guide RNAs (gRNAs) coded in DNA minicircles in the kinetoplast. There are many variations on this mechanism among species. mRNA editing and gRNA repertoires are almost unknown in T. cruzi. Here, gRNAs were inferred based on deep-sequenced minicircle hypervariable regions (mHVRs) and editing cascades were rebuilt in strains belonging to the six main T. cruzi lineages. Inferred gRNAs were clustered according to their sequence similarity to constitute gRNA classes. Extreme diversity of gRNA classes was observed, which implied highly divergent gRNA repertoires among different lineages, even within some lineages. In addition, a variable gRNA class redundancy (i.e., different gRNA classes editing the same mRNA region) was detected among strains. Some strains had upon four times more gRNA classes than others. Such variations in redundancy affected gRNA classes of all mRNAs in a concerted way, i.e., there are correlated variations in the number of gRNAs classes editing each mRNA. Interestingly, cascades were incomplete for components of the respiratory complex I in several strains. Finally, gRNA classes of different strains may potentially edit mitochondrial mRNAs from other lineages in the same way as they edit their own mitochondrial mRNAs, which is a prerequisite for biparental inheritance of minicircle in hybrids. We propose that genetic exchange and biparental inheritance of minicircles combined with minicircle drift due to (partial) random segregation of minicircles during kDNA replication is a suitable hypothesis to explain the divergences among strains and the high levels of gRNA redundancy in some strains. In addition, our results support that the complex I may not be required in some stages in the life cycle as previously shown and that linkage (in the same minicircle) of gRNAs that edit different mRNAs may prevent gRNA class lost in such stage.

A Novel Genotype and First Record of Trypanosoma lainsoni in Argentina.

Trypanosomes are a group of parasitic flagellates with medical and veterinary importance. Despite many species having been described in this genus, little is known about many of them. Here, we report a genetic and morphological characterization of trypanosomatids isolated from wild mammals from the Argentine Chaco region. Parasites were morphologically and ultrastructurally characterized by light microscopy and transmission electron microscopy. Additionally, 18s rRNA and gGAPDH genes were sequenced and analyzed using maximum likelihood and Bayesian inference. Morphological characterization showed clear characteristics associated with the Trypanosoma genus. The genetic characterization demonstrates that the studied isolates have identical sequences and a pairwise identity of 99% with Trypanosoma lainsoni, which belongs to the clade of lizards and snakes/rodents and marsupials. To date, this species had only been found in the Amazon region. Our finding represents the second report of T. lainsoni and the first record for the Chaco region. Furthermore, we ultrastructurally described for the first time the species. Finally, the host range of T. lainsoni was expanded (Leopardus geoffroyi, Carenivora, Felidae; and Calomys sp., Rodentia, Cricetidae), showing a wide host range for this species.

Evidence of hybridization, mitochondrial introgression and biparental inheritance of the kDNA minicircles in Trypanosoma cruzi I.

BACKGROUND: Genetic exchange in Trypanosoma cruzi is controversial not only in relation to its frequency, but also to its mechanism. Parasexual genetic exchange has been proposed based on laboratory hybrids, but population genomics strongly suggests meiosis in T. cruzi. In addition, mitochondrial introgression has been reported several times in natural isolates although its mechanism is not fully understood yet. Moreover, hybrid T. cruzi DTUs (TcV and TcVI) have inherited at least part of the kinetoplastic DNA (kDNA = mitochondrial DNA) from both parents. METHODOLOGY/PRINCIPAL FINDINGS: In order to address such topics, we sequenced and analyzed fourteen nuclear DNA fragments and three kDNA maxicircle genes in three TcI stocks which are natural clones potentially involved in events of genetic exchange. We also deep-sequenced (a total of 6,146,686 paired-end reads) the minicircle hypervariable region (mHVR) of the kDNA in such three strains. In addition, we analyzed the DNA content by flow cytometry to address cell ploidy. We observed that most polymorphic sites in nuclear loci showed a hybrid pattern in one cloned strain and the other two cloned strains were compatible as parental strains (or nearly related to the true parents). The three clones had almost the same ploidy and the DNA content was similar to the reference strain Sylvio (a nearly diploid strain). Despite maxicircle genes evolve faster than nuclear housekeeping ones, we detected no polymorphisms in the sequence of three maxicircle genes showing mito-nuclear discordance. Lastly, the hybrid stock shared 66% of its mHVR clusters with one putative parent and 47% with the other one; in contrast, the putative parental stocks shared less than 30% of the mHVR clusters between them. CONCLUSIONS/SIGNIFICANCE: The results suggest a reductive division, a natural hybridization, biparental inheritance of the minicircles in the hybrid and maxicircle introgression. The models including such phenomena and explaining the relationships between these three clones are discussed.

TcTASV Antigens of Trypanosoma cruzi: Utility for Diagnosis and High Accuracy as Biomarkers of Treatment Efficacy in Pediatric Patients.

The discovery and characterization of novel parasite antigens to improve the diagnosis of Trypanosoma cruzi by serological methods and for accurate and rapid follow-up of treatment efficiency are still needed. TcTASV is a T. cruzi-specific multigene family, whose products are expressed on the parasite stages present in the vertebrate host. In a previous work, a mix of antigens from subfamilies TcTASV-A and TcTASV-C (Mix A + C) was sensitive and specific to identify dogs with active infection of high epidemiological relevance. Here, TcTASV-A and TcTASV-C were assayed separately as well as together (Mix A + C) in an ELISA format on human samples. The Mix A + C presented moderate sensitivity (78%) but high diagnostic accuracy with a 100% of specificity, evaluated on healthy, leishmaniasic, and Strongyloides stercoralis infected patients. Moreover, antibody levels of pediatric patients showed-2 years posttreatment-diminished reactivity against the Mix A + C (P < 0.0001), pointing TcTASV antigens as promising tools for treatment follow-up.

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 105 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.

MLST Reveals a Separate and Novel Clonal Group for Acidovorax avenae Strains Causing Red Stripe in Sugarcane from Argentina.

Acidovorax spp. cause a wide range of economically important diseases in monocotyledonous and dicotyledonous plants, including sugarcane, corn, rice, oat, millet, foxtail watermelon, and orchid. In Argentina, the red stripe disease of sugarcane caused by Acidovorax avenae affects 30% of the milling stems with important economic losses. To explore the genetic diversity of this bacterium associated with red stripe in Argentina, multilocus sequence typing (MLST) was applied. This study included 15 local strains isolated from four different sugarcane planting regions and selected after random amplified polymorphic DNA analysis and reference strains of A. citrulli, A. avenae, and A. oryzae to investigate their phylogenetic relationships. MLST analysis resulted in five sequence types among the sugarcane A. avenae strains which constitute a clonal complex, meaning a common and close origin. Sugarcane strains were related to A. avenae from other hosts and distant to A. citrulli. Signals of frequent recombination in several lineages of A. avenae was detected and we observed that A. oryzae is closely related to A. avenae strains. This study provides valuable data in the field of epidemiological and evolutionary investigations of novel clone of A. avenae strains causing sugarcane red stripe. The knowledge of the genetic diversity and strain-host specificity are important to select the genotypes with the best response to the red stripe disease.

Introgression of the Kinetoplast DNA: An Unusual Evolutionary Journey in Trypanosoma cruzi.

INTRODUCTION: Phylogenetic relationships between different lineages of Trypanosoma cruzi, the agent of Chagas disease, have been controversial for several years. However, recent phylogenetic and phylogenomic analyses clarified the nuclear relationships among such lineages. However, incongruence between nuclear and kinetoplast DNA phylogenies has emerged as a new challenge. This incongruence implies several events of mitochondrial introgression at evolutionary level. However, the mechanism that gave origin to introgressed lineages is unknown. Here, I will review and discuss how maxicircles of the kinetoplast were horizontally and vertically transferred between different lineages of T. cruzi. CONCLUSION: Finally, I will discuss what we know - and what we don't - about the kDNA transference and inheritance in the context of sexual reproduction in this parasite.

Epidemiological modeling of Trypanosoma cruzi: Low stercorarian transmission and failure of host adaptive immunity explain the frequency of mixed infections in humans.

People living in areas with active vector-borne transmission of Chagas disease have multiple contacts with its causative agent, Trypanosoma cruzi. Reinfections by T. cruzi are possible at least in animal models leading to lower or even hardly detectable parasitaemia. In humans, although reinfections are thought to have major public health implications by increasing the risk of chronic manifestations of the disease, there is little quantitative knowledge about their frequency and the timing of parasite re-inoculation in the course of the disease. Here, we implemented stochastic agent-based models i) to estimate the rate of re-inoculation in humans and ii) to assess how frequent are reinfections during the acute and chronic stages of the disease according to alternative hypotheses on the adaptive immune response following a primary infection. By using a hybrid genetic algorithm, the models were fitted to epidemiological data of Argentinean rural villages where mixed infections by different genotypes of T. cruzi reach 56% in humans. To explain this percentage, the best model predicted 0.032 (0.008-0.042) annual reinfections per individual with 98.4% of them occurring in the chronic phase. In addition, the parasite escapes to the adaptive immune response mounted after the primary infection in at least 20% of the events of re-inoculation. With these low annual rates, the risks of reinfection during the typically long chronic stage of the disease stand around 14% (4%-18%) and 60% (21%-70%) after 5 and 30 years, with most individuals being re-infected 1-3 times overall. These low rates are better explained by the weak efficiency of the stercorarian mode of transmission than a highly efficient adaptive immune response. Those estimates are of particular interest for vaccine development and for our understanding of the higher risk of chronic disease manifestations suffered by infected people living in endemic areas.

Phylogenomics of Trypanosoma cruzi: Few evidence of TcI/TcII mosaicism in TcIII challenges the hypothesis of an ancient TcI/TcII hybridization.

Phylogenetic relationships among major lineages of Trypanosoma cruzi are still debatable. Particularly, it is controversial the origin of two main lineages: TcIII and TcIV. Some authors proposed that these lineages have been the result of an ancient hybridization between TcI and TcII, and this was one of the most accepted evolutionary models in the scientific community for several years. In the present paper we analyse several genomes of T. cruzi in order to examine if there is evidence supporting that TcIII is an ancient TcI/TcII hybrid. Our results show that TcIII is mainly related to TcI and not to TcII and there is few evidence of mosaicism for TcIII. Our results challenge the hypothesis of the ancient TcI/TcII hybridization.

Evolution of Trypanosoma cruzi: clarifying hybridisations, mitochondrial introgressions and phylogenetic relationships between major lineages

Several different models of Trypanosoma cruzi evolution have been proposed. These models suggest that scarce events of genetic exchange occurred during the evolutionary history of this parasite. In addition, the debate has focused on the existence of one or two hybridisation events during the evolution of T. cruzi lineages. Here, we reviewed the literature and analysed available sequence data to clarify the phylogenetic relationships among these different lineages. We observed that TcI, TcIII and TcIV form a monophyletic group and that TcIII and TcIV are not, as previously suggested, TcI-TcII hybrids. Particularly, TcI and TcIII are sister groups that diverged around the same time that a widely distributed TcIV split into two clades (TcIVS and TcIVN). In addition, we collected evidence that TcIII received TcIVS kDNA by introgression on several occasions. Different demographic hypotheses (surfing and asymmetrical introgression) may explain the origin and expansion of the TcIII group. Considering these hypotheses, genetic exchange should have been relatively frequent between TcIII and TcIVS in the geographic area in which their distributions overlapped. In addition, our results support the hypothesis that two independent hybridisation events gave rise to TcV and TcVI. Consequently, TcIVS kDNA was first transferred to TcIII and later to TcV and TcVI in TcII/TcIII hybridisation events.

Evolution of Trypanosoma cruzi: clarifying hybridisations, mitochondrial introgressions and phylogenetic relationships between major lineages.

Several different models of Trypanosoma cruzi evolution have been proposed. These models suggest that scarce events of genetic exchange occurred during the evolutionary history of this parasite. In addition, the debate has focused on the existence of one or two hybridisation events during the evolution of T. cruzi lineages. Here, we reviewed the literature and analysed available sequence data to clarify the phylogenetic relationships among these different lineages. We observed that TcI, TcIII and TcIV form a monophyletic group and that TcIII and TcIV are not, as previously suggested, TcI-TcII hybrids. Particularly, TcI and TcIII are sister groups that diverged around the same time that a widely distributed TcIV split into two clades (TcIVS and TcIVN). In addition, we collected evidence that TcIII received TcIVS kDNA by introgression on several occasions. Different demographic hypotheses (surfing and asymmetrical introgression) may explain the origin and expansion of the TcIII group. Considering these hypotheses, genetic exchange should have been relatively frequent between TcIII and TcIVS in the geographic area in which their distributions overlapped. In addition, our results support the hypothesis that two independent hybridisation events gave rise to TcV and TcVI. Consequently, TcIVS kDNA was first transferred to TcIII and later to TcV and TcVI in TcII/TcIII hybridisation events.

Experimental evidence of biological interactions among different isolates of Trypanosoma cruzi from the Chaco Region.

Many infectious diseases arise from co-infections or re-infections with more than one genotype of the same pathogen. These mixed infections could alter host fitness, the severity of symptoms, success in pathogen transmission and the epidemiology of the disease. Trypanosoma cruzi, the etiological agent of Chagas disease, exhibits a high biological variability often correlated with its genetic diversity. Here, we developed an experimental approach in order to evaluate biological interaction between three T. cruzi isolates belonging to different Discrete Typing Units (DTUs TcIII, TcV and TcVI). These isolates were obtained from a restricted geographical area in the Chaco Region. Different mixed infections involving combinations of two isolates (TcIII + TcV, TcIII + TcVI and TcV + TcVI) were studied in a mouse model. The parameters evaluated were number of parasites circulating in peripheral blood, histopathology and genetic characterization of each DTU in different tissues by DNA hybridization probes. We found a predominance of TcVI isolate in blood and tissues respect to TcIII and TcV; and a decrease of the inflammatory response in heart when the damage of mice infected with TcVI and TcIII + TcVI mixture were compared. In addition, simultaneous presence of two isolates in the same tissue was not detected. Our results show that biological interactions between isolates with different biological behaviors lead to changes in their biological properties. The occurrence of interactions among different genotypes of T. cruzi observed in our mouse model suggests that these phenomena could also occur in natural cycles in the Chaco Region.

Multilocus sequence typing approach for a broader range of species of Leishmania genus: describing parasite diversity in Argentina.

Leishmaniasis is a vector-borne protozoan infection affecting over 350 million people around the world. In Argentina cutaneous leishmaniasis is endemic in nine provinces and visceral leishmaniasis is spreading from autochthonous transmission foci in seven provinces. However, there is limited information about the diversity of the parasite in this country. Implementation of molecular strategies for parasite typing, particularly multilocus sequence typing (MLST), represents an improved approach for genetic variability and population dynamics analyses. We selected six loci as candidates implemented in reference strains and Argentinean isolates. Phylogenetic analysis showed high correlation with taxonomic classification of the parasite. Autochthonous Leishmania (Viannia) braziliensis showed higher genetic diversity than L. (Leishmania) infantum but low support was obtained for intra-L. braziliensis complex variants suggesting the need of new loci that contribute to phylogenetic resolution for an improved MLST or nested-MLST scheme. This study represents the first characterization of genetic variability of Leishmania spp. in Argentina.

Trypanosoma cruzi diversity in the Gran Chaco: mixed infections and differential host distribution of TcV and TcVI.

The transmission cycles of Trypanosoma cruzi in the Gran Chaco are complex networks involving domestic and wild components, whose interrelationships are not well understood. Knowing the circuit of transmission of the different Discrete Typing Units (DTUs) of T. cruzi in the complex environment of the Chaco region is relevant to understanding how the different components (reservoirs, vectors, ecotopes) interact. In the present study we identified the DTUs infecting humans and dogs in two rural areas of the Gran Chaco in Argentina, using molecular methods which avoid parasite culture. Blood samples of humans and dogs were typified by PCR-DNA blotting and hybridization assays with five specific DNA probes (TcI, TcII, TcIII, TcV and TcVI). PCR analyses were performed on seropositive human and dog samples and showed the presence of T. cruzi DNA in 41.7% (98/235) and 53% (35/66) samples, respectively. The identification of infective DTUs was determined in 83.6% (82/98) and 91.4% (32/35) in human and dog samples, respectively. Single infections (36.7% - 36/98) and a previously not detected high proportion of mixed infections (47.9% - 47/98) were found. In a 15.3% (15/98) of samples the infecting DTU was not identified. Among the single infections TcV was the most prevalent DTU (30.6% - 30/98) in human samples; while TcVI (42.8% - 15/35) showed the highest prevalence in dog samples. TcV/TcVI was the most prevalent mixed infection in humans (32.6% - 32/98); and TcI/TcVI (14.3% - 5/35) in dogs. Significant associations between TcV with humans and TcVI with dogs were detected. For the first time, the presence of TcIII was detected in humans from this region. The occurrence of one human infected whit TcIII (a principally wild DTU) could be suggested the emergence of this, in domestic cycles in the Gran Chaco.

Preponderant clonal evolution of Trypanosoma cruzi I from Argentinean Chaco revealed by Multilocus Sequence Typing (MLST).

Trypanosoma cruzi has been historically classified as a species with preponderant clonal evolution (PCE). However, with the advent of highly polymorphic markers and studies at geographically reduced scales, the PCE in T. cruzi was challenged. In fact, some studies have suggested that recombination in T. cruzi lineage I (TcI) is much more frequent than previously believed. Further analyses of TcI populations from different geographical regions of Latin America are needed to examine this hypothesis. In the present study, we contribute to this topic by analyzing the population structure of TcI from a restricted geographical area in the Chaco region, Argentina. We analyzed TcI isolates from different hosts and vectors using a Multilocus Sequence Typing (MLST) approach. These isolates were previously characterized by sequencing the spliced leader intergenic region (SL-IR). Low levels of incongruence and well-supported clusters for MLST dataset were obtained from the analyses. Moreover, high linkage disequilibrium was found and five repeated and overrepresented genotypes were detected. In addition, a good correspondence between SL-IR and MLST was observed which is expected under PCE. However, recombination is not ruled out because five out of 28 pairs of loci were incompatible with strict clonality and one possible genetic exchange event was detected. Overall, our results represent evidence of PCE in TcI from the study area. Finally, considering our findings we discuss the scenario for the genetic structure of TcI.

Optimized multilocus sequence typing (MLST) scheme for Trypanosoma cruzi.

Trypanosoma cruzi, the aetiological agent of Chagas disease possess extensive genetic diversity. This has led to the development of a plethora of molecular typing methods for the identification of both the known major genetic lineages and for more fine scale characterization of different multilocus genotypes within these major lineages. Whole genome sequencing applied to large sample sizes is not currently viable and multilocus enzyme electrophoresis, the previous gold standard for T. cruzi typing, is laborious and time consuming. In the present work, we present an optimized Multilocus Sequence Typing (MLST) scheme, based on the combined analysis of two recently proposed MLST approaches. Here, thirteen concatenated gene fragments were applied to a panel of T. cruzi reference strains encompassing all known genetic lineages. Concatenation of 13 fragments allowed assignment of all strains to the predicted Discrete Typing Units (DTUs), or near-clades, with the exception of one strain that was an outlier for TcV, due to apparent loss of heterozygosity in one fragment. Monophyly for all DTUs, along with robust bootstrap support, was restored when this fragment was subsequently excluded from the analysis. All possible combinations of loci were assessed against predefined criteria with the objective of selecting the most appropriate combination of between two and twelve fragments, for an optimized MLST scheme. The optimum combination consisted of 7 loci and discriminated between all reference strains in the panel, with the majority supported by robust bootstrap values. Additionally, a reduced panel of just 4 gene fragments displayed high bootstrap values for DTU assignment and discriminated 21 out of 25 genotypes. We propose that the seven-fragment MLST scheme could be used as a gold standard for T. cruzi typing, against which other typing approaches, particularly single locus approaches or systematic PCR assays based on amplicon size, could be compared.

How often do they have sex? A comparative analysis of the population structure of seven eukaryotic microbial pathogens.

The model of predominant clonal evolution (PCE) proposed for micropathogens does not state that genetic exchange is totally absent, but rather, that it is too rare to break the prevalent PCE pattern. However, the actual impact of this "residual" genetic exchange should be evaluated. Multilocus Sequence Typing (MLST) is an excellent tool to explore the problem. Here, we compared online available MLST datasets for seven eukaryotic microbial pathogens: Trypanosoma cruzi, the Fusarium solani complex, Aspergillus fumigatus, Blastocystis subtype 3, the Leishmania donovani complex, Candida albicans and Candida glabrata. We first analyzed phylogenetic relationships among genotypes within each dataset. Then, we examined different measures of branch support and incongruence among loci as signs of genetic structure and levels of past recombination. The analyses allow us to identify three types of genetic structure. The first was characterized by trees with well-supported branches and low levels of incongruence suggesting well-structured populations and PCE. This was the case for the T. cruzi and F. solani datasets. The second genetic structure, represented by Blastocystis spp., A. fumigatus and the L. donovani complex datasets, showed trees with weakly-supported branches but low levels of incongruence among loci, whereby genetic structuration was not clearly defined by MLST. Finally, trees showing weakly-supported branches and high levels of incongruence among loci were observed for Candida species, suggesting that genetic exchange has a higher evolutionary impact in these mainly clonal yeast species. Furthermore, simulations showed that MLST may fail to show right clustering in population datasets even in the absence of genetic exchange. In conclusion, these results make it possible to infer variable impacts of genetic exchange in populations of predominantly clonal micro-pathogens. Moreover, our results reveal different problems of MLST to determine the genetic structure in these organisms that should be considered.