Multilocus Sequence Typing Datasets Research Articles

Evaluation of various distance computation methods for construction of haplotype-based phylogenies from large MLST datasets

Multi-locus sequence typing (MLST) is widely used to investigate genetic relationships among eukaryotic taxa, including parasitic pathogens. MLST analysis workflows typically involve construction of alignment-based phylogenetic trees – i.e., where tree structures are computed from nucleotide differences observed in a multiple sequence alignment (MSA). Notably, alignment-based phylogenetic methods require that all isolates/taxa are represented by a single sequence. When multiple loci are sequenced these sequences may be concatenated to produce one tree that includes information from all loci. Alignment-based phylogenetic techniques are robust and widely used yet possess some shortcomings, including how heterozygous sites are handled, intolerance for missing data (i.e., partial genotypes), and differences in the way insertions-deletions (indels) are scored/treated during tree construction. In certain contexts, ‘haplotype-based’ methods may represent a viable alternative to alignment-based techniques, as they do not possess the aforementioned limitations. This is namely because haplotype-based methods assess genetic similarity based on numbers of shared (i.e., intersecting) haplotypes as opposed to similarities in nucleotide composition observed in an MSA. For haplotype-based comparisons, choosing an appropriate distance statistic is fundamental, and several statistics are available to choose from. However, a comprehensive assessment of various available statistics for their ability to produce a robust haplotype-based phylogenetic reconstruction has not yet been performed. We evaluated seven distance statistics by applying them to extant MLST datasets from the gastrointestinal parasite Cyclospora cayetanensis and two species of pathogenic nematode of the genus Strongyloides. We compare the genetic relationships identified using each statistic to epidemiologic, geographic, and host metadata. We show that Barratt’s heuristic definition of genetic distance was the most robust among the statistics evaluated. Consequently, it is proposed that Barratt’s heuristic represents a useful approach for use in the context of challenging MLST datasets possessing features (i.e., high heterozygosity, partial genotypes, and indel or repeat-based polymorphisms) that confound or preclude the use of alignment-based methods.

Machine learning-based analyses support the existence of species complexes for Strongyloides fuelleborni and Strongyloides stercoralis.

Human strongyloidiasis is a serious disease mostly attributable to Strongyloides stercoralis and to a lesser extent Strongyloides fuelleborni, a parasite mainly of non-human primates. The role of animals as reservoirs of human-infecting Strongyloides is ill-defined, and whether dogs are a source of human infection is debated. Published multi-locus sequence typing (MLST) studies attempt to elucidate relationships between Strongyloides genotypes, hosts, and distributions, but typically examine relatively few worms, making it difficult to identify population-level trends. Combining MLST data from multiple studies is often impractical because they examine different combinations of loci, eliminating phylogeny as a means of examining these data collectively unless hundreds of specimens are excluded. A recently-described machine learning approach that facilitates clustering of MLST data may offer a solution, even for datasets that include specimens sequenced at different combinations of loci. By clustering various MLST datasets as one using this procedure, we sought to uncover associations among genotype, geography, and hosts that remained elusive when examining datasets individually. Multiple datasets comprising hundreds of S. stercoralis and S. fuelleborni individuals were combined and clustered. Our results suggest that the commonly proposed 'two lineage' population structure of S. stercoralis (where lineage A infects humans and dogs, lineage B only dogs) is an over-simplification. Instead, S. stercoralis seemingly represents a species complex, including two distinct populations over-represented in dogs, and other populations vastly more common in humans. A distinction between African and Asian S. fuelleborni is also supported here, emphasizing the need for further resolving these taxonomic relationships through modern investigations.

Analysis and Comparison of the Phylogenetic Diversity Within Helicobacter pylori Isolates from Iranian and Global Populations by Multi-Locus Sequence Typing

Background: Helicobacter pylori is a spiral-shaped bacterium that chronically resides in the human gastric epithelium. Helicobacter pylori strains are genetically diverse organisms encoding informative genotypic polymorphisms that can reflect both ancestral and human migrations history. Objectives: In this study we aimed to investigate the phylogenetic diversity of Iranian H. pylori isolates by using multi-locus sequence typing (MLST). Methods: A total of 37 H. pylori isolates cultured from biopsy specimens from patients who referred to endoscopy unit at Taleghani Hospital in Tehran during 2015 to 2016 were included in this study. The MLST was performed for all the isolates using seven housekeeping genes, including atpA, efp, mutY, ppa, trpC, ureI and yphC . Phylogenetic analysis was performed using sequence datasets of our H. pylori strains intermingled with MLST datasets of neighboring countries and global populations. Results: Our results showed that approximately 38% (14/37) of the Iranian H. pylori strains were identical in at least three gene loci, with the ureI gene (48.6%, 18/37) found as the most identical allele. Additionally, the phylogenetic analyses showed that Iranian strains fall into distinct clusters, and were intermingled mostly amongst the hpEurope isolates, including isolates from Turkey, Palestine, Israel, Lebanon, Egypt, Greece, Germany, Netherlands, Spain, UK, Finland, and Italy. The present study revealed that Iranian H. pylori strains are drastically diverse, and are originally comparable to the ancestry of the hpEurope populations. Conclusions: In conclusion, although similarity of the Iranian H. pylori strains is shown in few gene loci comparable to those mostly find in the hpEurope population, diversity at nucleotide sequences and STs is common among them.

Wolbachia infection in Argentinean populations of Anastrepha fraterculus sp1: preliminary evidence of sex ratio distortion by one of two strains

BackgroundWolbachia, one of the most abundant taxa of intracellular Alphaproteobacteria, is widespread among arthropods and filarial nematodes. The presence of these maternally inherited bacteria is associated with modifications of host fitness, including a variety of reproductive abnormalities, such as cytoplasmic incompatibility, thelytokous parthenogenesis, host feminization and male-killing. Wolbachia has attracted much interest for its role in biological, ecological and evolutionary processes as well as for its potential use in novel and environmentally-friendly strategies for the control of insect pests and disease vectors including a major agricultural pest, the South American fruit fly, Anastrepha fraterculus Wiedemann (Diptera: Tephritidae).ResultsWe used wsp, 16S rRNA and a multilocus sequence typing (MLST) scheme including gatB, coxA, hcpA, fbpA, and ftsZ genes to detect and characterize the Wolbachia infection in laboratory strains and wild populations of A. fraterculus from Argentina. Wolbachia was found in all A. fraterculus individuals studied. Nucleotide sequences analysis of wsp gene allowed the identification of two Wolbachia nucleotide variants (named wAfraCast1_A and wAfraCast2_A). After the analysis of 76 individuals, a high prevalence of the wAfraCast2_A variant was found both, in laboratory (82%) and wild populations (95%). MLST analysis identified both Wolbachia genetic variants as sequence type 13. Phylogenetic analysis of concatenated MLST datasets clustered wAfraCast1/2_A in the supergroup A. Paired-crossing experiments among single infected laboratory strains showed a phenotype specifically associated to wAfraCast1_A that includes slight detrimental effects on larval survival, a female-biased sex ratio; suggesting the induction of male-killing phenomena, and a decreased proportion of females producing descendants that appears attributable to the lack of sperm in their spermathecae.ConclusionsWe detected and characterized at the molecular level two wsp gene sequence variants of Wolbachia both in laboratory and wild populations of A. fraterculus sp.1 from Argentina. Crossing experiments on singly-infected A. fraterculus strains showed evidence of a male killing-like mechanism potentially associated to the wAfraCast1_A - A. fraterculus interactions. Further mating experiments including antibiotic treatments and the analysis of early and late immature stages of descendants will contribute to our understanding of the phenotypes elicited by the Wolbachia variant wAfraCast1_A in A. fraterculus sp.1.

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.

MLSTest: Novel software for multi-locus sequence data analysis in eukaryotic organisms

Multi-locus sequence typing (MLST) is a frequently used genotyping method whose goal is the unambiguous assignment of microorganisms to genetic clusters. MLST typically involves analysis of DNA sequence results generated from several house-keeping gene loci. MLST remains the gold standard for molecular typing of many bacterial pathogens. Eukaryotic pathogens have also been the subject of MLST, however, few tools are available to deal with diploid sequence data. Here we present novel software for MLST data analysis tailored towards diploid Eukaryotes: MLSTest. This software meets various methods used in MLST and introduces some novel methodologies for the evaluation of the data set. In addition to construction of allelic profiles and basic clustering analysis, the MLSTest looks for network structures that suggest genetic exchange in BURST graphs. Additionally, it uses several simple methods for tree construction with the advantage of managing heterozygous or three-state sites. Additionally, the software analyses whether concatenation of fragments from different genes is suitable for the data set using different tests (bionj-incongruence length difference test, Templeton test). It evaluates how the incongruence is distributed across the tree using a variation of the localized incongruence length difference test based on a modified neighbour joining algorithm. We tested the last method in simulated datasets. We showed that is conservative (adequate type I error rate) and moderately to highly powerful as well as useful to localize incongruences in two bacterial and two eukaryotic MLST datasets. MLSTest was also designed for developing MLST schemes. It thus has tools to optimize locus combinations and to reduce the number of targets required for typing. MLSTest also analyses whether the discriminatory power of the typing scheme is increased by including more loci. We evaluated the software over simulated and real datasets from bacterial and eukaryotic microorganisms. The software is freely available at http://www.ipe.unsa.edu.ar/software.

Geographically Structured Populations of Cryptococcus neoformans Variety grubii in Asia Correlate with HIV Status and Show a Clonal Population Structure

Cryptococcosis is an important fungal disease in Asia with an estimated 140,000 new infections annually the majority of which occurs in patients suffering from HIV/AIDS. Cryptococcus neoformans variety grubii (serotype A) is the major causative agent of this disease. In the present study, multilocus sequence typing (MLST) using the ISHAM MLST consensus scheme for the C. neoformans/C. gattii species complex was used to analyse nucleotide polymorphisms among 476 isolates of this pathogen obtained from 8 Asian countries. Population genetic analysis showed that the Asian C. neoformans var. grubii population shows limited genetic diversity and demonstrates a largely clonal mode of reproduction when compared with the global MLST dataset. HIV-status, sequence types and geography were found to be confounded. However, a correlation between sequence types and isolates from HIV-negative patients was observed among the Asian isolates. Observations of high gene flow between the Middle Eastern and the Southeastern Asian populations suggest that immigrant workers in the Middle East were originally infected in Southeastern Asia.

Practicalities of Using Non-Local or Non-Recent Multilocus Sequence Typing Data for Source Attribution in Space and Time of Human Campylobacteriosis

In this study, 1208 Campylobacter jejuni and C. coli isolates from humans and 400 isolates from chicken, collected in two separate periods over 12 years in The Netherlands, were typed using multilocus sequence typing (MLST). Statistical evidence was found for a shift of ST frequencies in human isolates over time. The human MLST data were also compared to published data from other countries to determine geographical variation. Because only MLST typed data from chicken, taken from the same time point and spatial location, were available in addition to the human data, MLST datasets for other Campylobacter reservoirs from selected countries were used. The selection was based on the degree of similarity of the human isolates between countries. The main aim of this study was to better understand the consequences of using non-local or non-recent MLST data for attributing domestically acquired human Campylobacter infections to specific sources of origin when applying the asymmetric island model for source attribution. In addition, a power-analysis was done to find the minimum number of source isolates needed to perform source attribution using an asymmetric island model. This study showed that using source data from other countries can have a significant biasing effect on the attribution results so it is important to carefully select data if the available local data lack in quality and/or quantity. Methods aimed at reducing this bias were proposed.

Molecular epidemiology of Italian clinicalCryptococcus neoformansvar.grubiiisolates

Cryptococcus neoformans variety grubii is the major etiological agent of cryptococcal meningitis in both immunocompromised and immunocompetent patients. The current PCR-based molecular methods are not sufficient to discriminate among the different populations of this yeast. Therefore, the aim of the present study was to investigate the genotypes of the Italian clinical C. neoformans var. grubii isolates by multilocus sequence typing (MLST). A total of 53 isolates, each representative of a single case, were studied. Genotyping was performed using the ISHAM Cryptococcus MLST consensus scheme and the results were compared to the publically available global C. neoformans var. grubii MLST dataset. A total of 16 genotypes were identified; 14 were new genotypes, one was identical to sequence type (ST) ST81, which had been previously reported from Thailand, and one to ST23 already identified in Uganda, the USA and Korea. Sequence type ST61 was the most numerous, including 16 isolates. Network phylogenetic analysis showed that the Italian isolates could be divided into at least three clusters with similarities with those recovered in Africa, Asia and Americas. Distribution of the STs among the isolates could not be correlated to the hospital in which they were recovered or to the HIV status of the patients. The majority of the isolates belonged to the molecular type VNI; three belonged to the rare molecular type VNII and one to the VNB group, which until now had not been described in Europe. The results reveal that the Italian C. neoformans var. grubii population presents a distinct variability, displaying a high number of new genotypes, and probably recombines sexually.

Phylogenetic evaluation of the ‘Typhimurium’ complex of Salmonella strains using a seven-gene multi-locus sequence analysis

Salmonella enterica comprises over 2500 serovars, many of which are significant foodborne pathogens in humans. The ability to subtype these microbes is difficult due to the highly clonal nature of many Salmonella strains and a lack of congruence among traditional typing approaches. This work examines the phylogenetic utility of a multi-locus sequence typing (MLST) approach to discriminate between members of a closely related collection of salmonellae, the Salmonella reference collection A (SARA). This 72 strain collection, referred to as the ‘Typhimurium’ complex, consists of S. Typhimurium and its four closest serological relatives. In this analysis, nucleotide sequences from seven housekeeping genes ( aroC, dnaN, hemD, hisD, purE, sucA, and thrA) were PCR amplified, sequenced, and combined into a single concatenated character matrix providing 3360 bp for cladistic analysis. The resultant most parsimonious tree yielded seven clades of Salmonella strains that partitioned largely along serovar divisions within the collection except for five ‘Paratyphi B’ strains, two ‘Saintpaul’ strains, and two ‘Typhimurium’ strains. Convergence in the SARA tree was approximately 20% indicating that the vast majority of sequence changes were phylogenetically informative. Despite a high consistency among nucleotide substitutions, analysis of congruence identified several SARA strains with recombinant alleles in the concatenated matrix. These findings point to important differences among phylogenetic contributions made by the individual genes comprising this MLST dataset.

Mating is rare within as well as between clades of the human pathogen Candida albicans

Candida albicans is a diploid yeast that can undergo mating and a parasexual cycle, but is apparently unable to undergo meiosis. Characterization of the population structure of C. albicans has shown that reproduction is largely clonal and that mating, if it occurs, is rare or limited to genetically related isolates. Because molecular typing has delineated distinct clades in C. albicans, we have tested whether recombination was common within clades, but rare between clades. Two hundred and three C. albicans isolates have been subjected to multilocus sequence typing (MLST) and the haplotypes at heterozygous MLST genotypes characterized. The C. albicans isolates were distributed among nine clades, of which five corresponded to those previously identified by Ca3 fingerprinting. In each of these clades with more than 10 isolates, polymorphic nucleotide positions located on between 3 and 4 of the six loci were in Hardy–Weinberg disequilibrium. Moreover, each of these polymorphic sites contained excess heterozygotes. This was confirmed by an expanded analysis performed on a recently published MLST dataset for 1044 isolates. On average, 66% of polymorphic positions in the individual clades were in significant excess of heterozygotes over the five clades. These data indicate that mating within clades as well as self-fertilization are both limited and that C. albicans clades do not represent a collection of cryptic species. The study of haplotypes at heterozygous loci performed on our dataset indicates that loss of heterozygosity events due to mitotic recombination is moderately common in natural populations of C. albicans. The maintenance of substantial heterozygosity despite relatively frequent loss of heterozygosity could result from a selective advantage conferred by heterozygosity.

Molecular Evolution Perspectives on Intraspecific Lateral DNA Transfer of Topoisomerase and Gyrase Loci in Streptococcus pneumoniae , with Implications for Fluoroquinolone Resistance Development and Spread

Fluoroquinolones are an important class of antibiotics for the treatment of infections arising from the gram-positive respiratory pathogen Streptococcus pneumoniae. Although there is evidence supporting interspecific lateral DNA transfer of fluoroquinolone target loci, no studies have specifically been designed to assess the role of intraspecific lateral transfer of these genes in the spread of fluoroquinolone resistance. This study involves a comparative evolutionary perspective, in which the evolutionary history of a diverse set of S. pneumoniae clinical isolates is reconstructed from an expanded multilocus sequence typing data set, with putative recombinants excluded. This control history is then assessed against networks of each of the four fluoroquinolone target loci from the same isolates. The results indicate that although the majority of fluoroquinolone target loci from this set of 60 isolates are consistent with a clonal dissemination hypothesis, 3 to 10% of the sequences are consistent with an intraspecific lateral transfer hypothesis. Also evident were examples of interspecific transfer, with two isolates possessing a parE-parC gene region arising from viridans group streptococci. The Spain 23F-1 clone is the most dominant fluoroquinolone-nonsusceptible clone in this set of isolates, and the analysis suggests that its members act as frequent donors of fluoroquinolone-nonsusceptible loci. Although the majority of fluoroquinolone target gene sequences in this set of isolates can be explained on the basis of clonal dissemination, a significant number are more parsimoniously explained by intraspecific lateral DNA transfer, and in situations of high S. pneumoniae population density, such events could be an important means of resistance spread.

EBURST: inferring patterns of evolutionary descent among clusters of related bacterial genotypes from multilocus sequence typing data.

The introduction of multilocus sequence typing (MLST) for the precise characterization of isolates of bacterial pathogens has had a marked impact on both routine epidemiological surveillance and microbial population biology. In both fields, a key prerequisite for exploiting this resource is the ability to discern the relatedness and patterns of evolutionary descent among isolates with similar genotypes. Traditional clustering techniques, such as dendrograms, provide a very poor representation of recent evolutionary events, as they attempt to reconstruct relationships in the absence of a realistic model of the way in which bacterial clones emerge and diversify to form clonal complexes. An increasingly popular approach, called BURST, has been used as an alternative, but present implementations are unable to cope with very large data sets and offer crude graphical outputs. Here we present a new implementation of this algorithm, eBURST, which divides an MLST data set of any size into groups of related isolates and clonal complexes, predicts the founding (ancestral) genotype of each clonal complex, and computes the bootstrap support for the assignment. The most parsimonious patterns of descent of all isolates in each clonal complex from the predicted founder(s) are then displayed. The advantages of eBURST for exploring patterns of evolutionary descent are demonstrated with a number of examples, including the simple Spain(23F)-1 clonal complex of Streptococcus pneumoniae, "population snapshots" of the entire S. pneumoniae and Staphylococcus aureus MLST databases, and the more complicated clonal complexes observed for Campylobacter jejuni and Neisseria meningitidis.

The relative contributions of recombination and mutation to the divergence of clones of Neisseria meningitidis.

Multilocus sequence typing (MLST) is a recently developed nucleotide sequence-based method for the definitive assignment of isolates within bacterial populations to specific clones. MLST uses the same principles as multilocus enzyme electrophoresis and provides data that can be used to investigate aspects of the population genetics and evolution of bacterial species. We used an MLST data set consisting of the sequences of approximately 450-bp fragments from seven housekeeping loci from a large strain collection of Neisseria meningitidis to estimate the relative impact of recombination compared with point mutation in the diversification of N. meningitidis clonal complexes. 126 meningococcal isolates were assigned to 10 clonal complexes, 9 of which contained minor clonal variants. The allelic variation within each complex was classified as a recombinational exchange or a putative point mutation through a comparison of the sequences of each variant allele with that of the allele typically found in the clonal complex. The nine clonal complexes contained a total of 23 allelic variants, and analysis of the sequences of these variant alleles revealed that a single nucleotide site in a meningococcal housekeeping gene is at least 80-fold more likely to change as a result of recombination than as a result of mutation. This value is estimated to be 10-50-fold for Escherichia coli and approximately 50-fold for Streptococcus pneumoniae.