Inadequacies of minimum spanning trees in molecular epidemiology

Outbreak of Salmonella Typhimurium linked to Swedish pre-washed rocket salad, Sweden, September to November 2022

In September 2022, the Public Health Agency of Sweden observed an increase in domestic Salmonella Typhimurium cases through the Swedish electronic notification system, and an outbreak strain was identified with whole genome sequencing. Overall, 109 cases with symptom onset between 17 September and 24 November 2022 were reported from 20 of 21 Swedish regions. The median age of cases was 52 years (range 4-87 years) and 62% were female. A case-control study found cases to be associated with consumption of rocket salad (adjusted odds ratio (aOR) = 4.9; 95% confidence interval (CI): 2.4-10, p value < 0.001) and bagged mixed salad (aOR = 4.0; 95% CI: 1.9-8.1, p value < 0.001). Trace-back, supported by Finnish authorities who identified the Swedish outbreak strain in a Finnish cluster during the same time period, identified rocket salad, cultivated, pre-washed and pre-packed in Sweden as the likely source of the outbreak. No microbiological analyses of rocket salad were performed. Our investigation indicates that bagged leafy greens such as rocket salad, regardless of pre-washing procedures in the production chain, may contain Salmonella and cause outbreaks, posing a health risk to consumers. We emphasise the need for primary producers of leafy greens to identify possible contamination points to prevent outbreaks.

Investigation of an international water polo tournament in Czechia as a potential source for early introduction of the SARS-CoV-2 Omicron variant into Belgium, Switzerland and Germany, November 2021

BackgroundThe earliest recognised infections by the SARS-CoV-2 Omicron variant (Pango lineage B.1.1.529) in Belgium and Switzerland suggested a connection to an international water polo tournament, held 12-14 November 2021 in Brno, Czechia.AimTo study the arrival and subsequent spread of the Omicron variant in Belgium and Switzerland, and understand the overall importance of this international sporting event on the number of infections in the two countries.MethodsWe performed intensive forward and backward contact tracing in both countries, supplemented by phylogenetic investigations using virus sequences of the suspected infection chain archived in public databases.ResultsThrough contact tracing, we identified two and one infected athletes of the Belgian and Swiss water polo teams, respectively, and subsequently also three athletes from Germany. In Belgium and Switzerland, four and three secondary infections, and three and one confirmed tertiary infections were identified. Phylogenetic investigation demonstrated that this sporting event played a role as the source of infection, but without a direct link with infections from South Africa and not as a superspreading event; the virus was found to already be circulating at that time in the countries involved.ConclusionThe SARS-CoV-2 Omicron variant started to circulate in Europe several weeks before its identification in South Africa on 24 November 2021. Accordingly, it can be assumed that travel restrictions are usually implemented too late to prevent the spread of newly detected SARS-CoV-2 variants to other regions. Phylogenetic analysis may modify the perception of an apparently clear result of intensive contact tracing.

Graph-Based Reconstruction and Analysis of Disease Transmission Networks Using Viral Genomic Data

Understanding the patterns of viral disease transmissions helps establish public health policies and aids in controlling and ending a disease outbreak. Classical methods for studying disease transmission dynamics that rely on epidemiological data, such as times of sample collection and duration of exposure intervals, struggle to provide desired insight due to limited informativeness of such data. A more precise characterization of disease transmissions may be acquired from sequencing data that reveal genetic distance between viral genomes in patient samples. Indeed, genetic distance between viral strains present in hosts contains valuable information about transmission history, thus motivating the design of methods that rely on genomic data to reconstruct a directed disease transmission network, detect transmission clusters, and identify significant network nodes (e.g., super-spreaders). In this article, we present a novel end-to-end framework for the analysis of viral transmissions utilizing viral genomic (sequencing) data. The proposed framework groups infected hosts into transmission clusters based on the reconstructed viral strains infecting them; the genetic distance between a pair of hosts is calculated using Earth Mover's Distance, and further used to infer transmission direction between the hosts. To quantify the significance of a host in the transmission network, the importance score is calculated by a graph convolutional autoencoder. The viral transmission network is represented by a directed minimum spanning tree utilizing the Edmond's algorithm modified to incorporate constraints on the importance scores of the hosts. The proposed framework outperforms state-of-the-art techniques for the analysis of viral transmission dynamics in several experiments on semiexperimental as well as experimental data.

Genotyping and biofilm formation of Mycoplasma hyopneumoniae and their association with virulence

Mycoplasma hyopneumoniae, the causative agent of swine respiratory disease, demonstrates differences in virulence. However, factors associated with this variation remain unknown. We herein evaluated the association between differences in virulence and genotypes as well as phenotype (i.e., biofilm formation ability). Strains 168 L, RM48, XLW-2, and J show low virulence and strains 232, 7448, 7422, 168, NJ, and LH show high virulence, as determined through animal challenge experiments, complemented with in vitro tracheal mucosa infection tests. These 10 strains with known virulence were then subjected to classification via multilocus sequence typing (MLST) with three housekeeping genes, P146-based genotyping, and multilocus variable-number tandem-repeat analysis (MLVA) of 13 loci. MLST and P146-based genotyping identified 168, 168 L, NJ, and RM48 as the same type and clustered them in a single branch. MLVA assigned a different sequence type to each strain. Simpson’s index of diversity indicates a higher discriminatory ability for MLVA. However, no statistically significant correlation was found between genotypes and virulence. Furthermore, we investigated the correlation between virulence and biofilm formation ability. The strains showing high virulence demonstrate strong biofilm formation ability, while attenuated strains show low biofilm formation ability. Pearson correlation analysis revealed a significant positive correlation between biofilm formation ability and virulence. To conclude, there was no association between virulence and our genotyping data, but virulence was found to be significantly associated with the biofilm formation ability of M. hyopneumoniae.

Molecular epidemiology of community- and hospital-associated Clostridioides difficile infections in Jönköping, Sweden, October 2017 - March 2018

Clostridioides difficile infections (CDIs) in Sweden are mostly hospital-associated (HA) with limited knowledge regarding community-associated (CA) infections. Here, we investigated the molecular epidemiology of clinical isolates of CA-CDI and HA-CDI in a Swedish county. Data and isolates (n = 156) of CDI patients (n = 122) from Jönköping county, October 2017-March 2018, were collected and classified as CA (without previous hospital care or onset ≤2 days after admission or >12 weeks after discharge from hospital) or HA (onset >3 days after hospital admission or within 4 weeks after discharge). Molecular characterization of isolates included PCR ribotyping (n = 156 isolates) and whole genome sequencing with single nucleotide polymorphisms (SNP) analysis (n = 53 isolates). We classified 47 patients (39%) as CA-CDI and 75 (61%) as HA-CDI. Between CA-CDI and HA-CDI patients, we observed no statistically significant differences regarding gender, age, 30-day mortality or recurrence. Ribotype 005 (RR 3.1; 95% CI: 1.79-5.24) and 020 (RR 2.5; 95% CI: 1.31-4.63) were significantly associated with CA-CDI. SNP analysis identified seven clusters (0-2 SNP difference) involving 17/53 isolates of both CA-CDI and HA-CDI. Molecular epidemiology differed between CA-CDI and HA-CDI and WGS analysis suggests transmission of CDI within and between hospitals and communities.

Outbreak of Salmonella Newport linked to imported frozen cooked crayfish in dill brine, Sweden, July to November 2019

In autumn 2019, the Public Health Agency of Sweden identified a cluster of Salmonella Newport cases by whole genome sequencing (WGS). Cases' distribution in place and time indicated a nation-wide ongoing outbreak. An investigation was initiated to identify the source and prevent further cases. We conducted a case-case study based on notified salmonellosis cases and a Salmonella trawling questionnaire, comparing 20 outbreak cases and 139 control cases. Food exposures were compared by adjusted odds ratios (aOR) with 95% confidence interval (CI) using logistic regression. Implicated foods were sampled. Outbreak cases were more likely to have consumed crayfish (aOR = 26; 95% CI: 6.3-105). One specific brand of imported frozen, pre-cooked whole crayfish in dill brine was identified as the source. Salmonella Newport was later detected in different batches from retail and in one sample from border control. Isolates from food samples clustered with the human outbreak strain by WGS. Although the retailer made a complete recall, two more cases were identified long afterwards. This investigation demonstrated the successful use of a case-case study and targeted microbiological testing to identify the source. The immediate action taken by the retailer was important to confirm the source and stop the outbreak.

Identification of a canine coronavirus in Australian racing Greyhounds

Coronavirus infection can cause a range of syndromes, which in dogs can include mild-to-severe enteritis that generally resolves rapidly. Fatalities can occur from coinfection with other pathogens, including canine parvovirus. Between late December 2019 and April 2020, canine coronavirus (CCoV) was detected in Australian racing Greyhounds that displayed signs of gastrointestinal disease. The CCoV was genotyped using high-throughput sequencing, recovering 98.3% of a type IIb CCoV, generally thought to cause a mild but highly contagious enteric disease. The Australian CCoV was almost identical (99.9%, whole-genome sequence) to another CCoV associated with an outbreak of severe vomiting in dogs in the United Kingdom at the same time (December 2019-March 2020).

Trends in the molecular epidemiology and population genetics of emerging Sporothrix species

Sporothrix (Ophiostomatales) comprises species that are pathogenic to humans and other mammals as well as environmental fungi. Developments in molecular phylogeny have changed our perceptions about the epidemiology, host-association, and virulence of Sporothrix. The classical agent of sporotrichosis, Sporothrix schenckii, now comprises several species nested in a clinical clade with S. brasiliensis, S. globosa, and S. luriei. To gain a more precise view of outbreaks dynamics, structure, and origin of genetic variation within and among populations of Sporothrix, we applied three sets of discriminatory AFLP markers (#3 EcoRI-GA/MseI-TT, #5 EcoRI-GA/MseI-AG, and #6 EcoRI-TA/MseI-AA) and mating-type analysis to a large collection of human, animal and environmental isolates spanning the major endemic areas. A total of 451 polymorphic loci were amplified in vitro from 188 samples, and revealed high polymorphism information content (PIC = 0.1765-0.2253), marker index (MI = 0.0001-0.0002), effective multiplex ratio (E = 15.1720-23.5591), resolving power (Rp = 26.1075-40.2795), discriminating power (D = 0.9766-0.9879), expected heterozygosity (H = 0.1957-0.2588), and mean heterozygosity (Havp  = 0.000007-0.000009), demonstrating the effectiveness of AFLP markers to speciate Sporothrix. Analysis using the program structure indicated three genetic clusters matching S. brasiliensis (population 1), S. schenckii (population 2), and S. globosa (population 3), with the presence of patterns of admixture amongst all populations. AMOVA revealed highly structured clusters (PhiPT = 0.458-0.484, P < 0.0001), with roughly equivalent genetic variability within (46-48 %) and between (52-54 %) populations. Heterothallism was the exclusive mating strategy, and the distributions of MAT1-1 or MAT1-2 idiomorphs were not significantly skewed (1:1 ratio) for S. schenckii (χ2 = 2.522; P = 0.1122), supporting random mating. In contrast, skewed distributions were found for S. globosa (χ2 = 9.529; P = 0.0020) with a predominance of MAT1-1 isolates, and regional differences were highlighted for S. brasiliensis with the overwhelming occurrence of MAT1-2 in Rio de Janeiro (χ2 = 14.222; P = 0.0002) and Pernambuco (χ2 = 7.364; P = 0.0067), in comparison to a higher prevalence of MAT1-1 in the Rio Grande do Sul (χ2 = 7.364; P = 0.0067). Epidemiological trends reveal the geographic expansion of cat-transmitted sporotrichosis due to S. brasiliensis via founder effect. These data support Rio de Janeiro as the centre of origin that has led to the spread of this disease to other regions in Brazil. Our ability to reconstruct the source, spread, and evolution of the ongoing outbreaks from molecular data provides high-quality information for decision-making aimed at mitigating the progression of the disease. Other uses include surveillance, rapid diagnosis, case connectivity, and guiding access to appropriate antifungal treatment.

Exploring genetic diversity, population structure, and phylogeography in Paracoccidioides species using AFLP markers

Paracoccidioidomycosis (PCM) is a life-threatening systemic fungal infection acquired after inhalation of Paracoccidioides propagules from the environment. The main agents include members of the P. brasiliensis complex (phylogenetically-defined species S1, PS2, PS3, and PS4) and P. lutzii. DNA-sequencing of protein-coding loci (e.g., GP43, ARF, and TUB1) is the reference method for recognizing Paracoccidioides species due to a lack of robust phenotypic markers. Thus, developing new molecular markers that are informative and cost-effective is key to providing quality information to explore genetic diversity within Paracoccidioides. We report using new amplified fragment length polymorphism (AFLP) markers and mating-type analysis for genotyping Paracoccidioides species. The bioinformatic analysis generated 144 in silico AFLP profiles, highlighting two discriminatory primer pairs combinations (#1 EcoRI-AC/MseI-CT and #2 EcoRI-AT/MseI-CT). The combinations #1 and #2 were used in vitro to genotype 165 Paracoccidioides isolates recovered from across a vast area of South America. Considering the overall scored AFLP markers in vitro (67-87 fragments), the values of polymorphism information content (PIC = 0.3345-0.3456), marker index (MI = 0.0018), effective multiplex ratio (E = 44.6788-60.3818), resolving power (Rp = 22.3152-34.3152), discriminating power (D = 0.5183-0.5553), expected heterozygosity (H = 0.4247-0.4443), and mean heterozygosity (H avp  = 0.00002-0.00004), demonstrated the utility of AFLP markers to speciate Paracoccidioides and to dissect both deep and fine-scale genetic structures. Analysis of molecular variance (AMOVA) revealed that the total genetic variance (65-66 %) was due to variability among P. brasiliensis complex and P. lutzii (PhiPT = 0.651-0.658, P < 0.0001), supporting a highly structured population. Heterothallism was the exclusive mating strategy, and the distributions of MAT1-1 or MAT1-2 idiomorphs were not significantly skewed (1:1 ratio) for P. brasiliensis s. str. (χ2 = 1.025; P = 0.3113), P. venezuelensis (χ2 = 0.692; P = 0.4054), and P. lutzii (χ2 = 0.027; P = 0.8694), supporting random mating within each species. In contrast, skewed distributions were found for P. americana (χ2 = 8.909; P = 0.0028) and P. restrepiensis (χ2 = 4.571; P = 0.0325) with a preponderance of MAT1-1. Geographical distributions confirmed that P. americana, P. restrepiensis, and P. lutzii are more widespread than previously thought. P. brasiliensis s. str. is by far the most widely occurring lineage in Latin America countries, occurring in all regions of Brazil. Our new DNA fingerprint assay proved to be rapid, reproducible, and highly discriminatory, to give insights into the taxonomy, ecology, and epidemiology of Paracoccidioides species, guiding disease-control strategies to mitigate PCM.

Investigations on spreading of PRRSV among swine herds by improved minimum spanning network analysis

In Hungary, the economic losses caused by porcine reproductive and respiratory syndrome virus (PRRSV) led to the launching of a national PRRSV Eradication Program. An important element of the program was investigating the spread of PRRSV among swine herds and the possible ways of introduction by sequencing of the open reading frame 5 (ORF5) gene. However, the classical phylogenetic tree presentation cannot explain several genetic relationships clearly, while more precise visualization can be represented by network tree diagram. In this paper, we describe a practical and easy-to-follow enriched minimum spanning similarity network application for improved representation of phylogenetic relations among viral strains. This method eliminated the necessity of applying a predefined, arbitrary cut-off or computationally extensive algorithms. The network-based visualization allowed processing and visualizing large amount of data equally for the laboratory, private and official veterinarians, and helped identify the potential connections between different viral sequences that support data-driven decisions in the eradication program. By applying network analysis, previously unknown epidemiological connections between infected herds were identified, and virus spreading was analyzed within short period of time. In our study, we successfully built and applied network analysis tools in the course of the Hungarian PRRSV Eradication Program.

Outbreak of gastroenteritis highlighting the diagnostic and epidemiological challenges of enteroinvasive Escherichia coli, County of Halland, Sweden, November 2017

An outbreak of gastroenteritis with 83 cases occurred at a conference venue in November 2017 in Halland County, Sweden. Stool samples from two venue visitors and a symptomatic secondary case attributed to household transmission were PCR-positive for the ipaH gene, a target found in both Shigella spp. and enteroinvasive Escherichia coli (EIEC). EIEC was isolated from stool samples and whole genome sequencing analysis confirmed EIEC O96:H19 to be the aetiological agent. A cohort study was conducted among venue attendees and employees and the findings implicated contaminated leafy greens as the vehicle of infection, however, no microbiological evidence could support the study results. Here, we report the investigation into the first recorded EIEC outbreak in Sweden and illustrate the challenges associated with the differential laboratory diagnostics of Shigella/EIEC in an outbreak setting.

Outbreak of unusual H2S-negative monophasic Salmonella Typhimurium strain likely associated with small tomatoes, Sweden, August to October 2019

Sweden is investigating an outbreak of monophasic Salmonella Typhimurium. Eighty-two nationally-distributed cases have been confirmed, with date of symptom onset between 28 August and 29 October. Cases were 51 years of age on average (range: 0–89) and the majority of cases were female (62%). A case–control study was conducted and suggested small tomatoes as source of the outbreak (adjusted odds ratio (OR): 10.8, 95% confidence interval (CI): 4.15-112.68, p value < 0.001), and a trace-back investigation led to a single, non-Swedish producer in Europe. Both the Salmonella strain and the source of the outbreak are rarely encountered in Europe. Results from investigation at the producer are pending.

Degree and centrality-based approaches in network-based variable selection: Insights from the Singapore Longitudinal Aging Study

We describe a network-based method to obtain a subset of representative variables from clinical data of subjects of the second Singapore Longitudinal Aging Study (SLAS-2), while preserving to a good extent the predictive performance of the full set with regards to a multi-faceted index of successful aging, SAGE. To examine differences in predictive performance of high-degree nodes (“hubs”) and high-centrality ones (“cores”), we implement four subsetting strategies (two degree-based, two centrality-based) and obtain four surrogate sets of variables, which we use as input features for machine learning models to predict the SAGE index of subjects. All four models have variables belonging to the physical, cardiovascular, cognitive and immunological domains among their fifteen most important predictors. A fifth domain (leisure-time activities, LTA) is also present in some form. From a comparison of the surrogate sets’ size and predictive performance, a centrality-based approach (selection of the most central variable-nodes within each cluster) yielded the smallest-sized surrogate set, while having high prediction accuracy (measured by its model’s area-under-curve, AUC) in comparison to its analogous degree-based strategy (selection of the highest-degree nodes per cluster). Inclusion of the next most-central variables yielded negligible changes in predictive performance while more than doubling the surrogate set size. The centrality-based approach thus yields a surrogate set which offers a good balance between number of variables and prediction performance, and can act as a representative subset of the SLAS-2 clinical dataset.

Clostridium difficile in patients attending tuberculosis hospitals in Cape Town, South Africa, 2014-2015

Background

Diarrhoea due to Clostridium difficile infection (CDI) poses a significant burden on healthcare systems around the world. However, there are few reports on the current status of the disease in sub-Saharan Africa.

Objectives

This study examined the occurrence of CDI in a South African population of tuberculosis patients, as well as the molecular epidemiology and antibiotic susceptibility profiles of C. difficile strains responsible for disease.

Methods

Toxigenic C. difficile in patients with suspected CDI attending two specialist tuberculosis hospitals in the Cape Town area were detected using a PCR-based diagnostic assay (Xpert^® C. difficile). C. difficile strains isolated from PCR-positive specimens were characterised by ribotyping, multilocus variable-number tandem-repeat analysis and antibiotic susceptibility testing.

Results

The period prevalence of CDI was approximately 70.07 cases per 1000 patient admissions. Strains belonging to ribotype 017 (RT017) made up over 95% of the patient isolates and all of them were multi-drug resistant. Multilocus variable-number tandem-repeat analysis revealed several clusters of highly related C. difficile RT017 strains present in tuberculosis patients in several wards at each hospital.

Conclusion

Tuberculosis patients represent a population that may be at an increased risk of developing CDI and, in addition, may constitute a multi-drug resistant reservoir of this bacterium. This warrants further investigation and surveillance of the disease in this patient group and other high-risk patient groups in sub-Saharan Africa.

Antimicrobial Resistance, Virulence Factors, and Genetic Profiles of Vibrio parahaemolyticus from Seafood

Vibrio parahaemolyticus is a widespread bacterium in the marine environment and is responsible for gastroenteritis in humans. Foodborne infections are mainly associated with the consumption of contaminated raw or undercooked fish and shellfish. The aim of this study was to determine the antimicrobial resistance, virulence factors, and genetic profiles of V. parahaemolyticus isolates from seafood originating from different countries. A total of 104 (17.5%) isolates were recovered from 595 analyzed samples. The isolates were tested for the presence of the tdh and trh genes, involved in the pathogenesis of V. parahaemolyticus infections in humans, and these genes were detected in 3 (2.9%) and 11 (10.6%) isolates, respectively. The trh-positive isolates also possessed the ure gene, which is responsible for urease production. Moreover, the activity of protease A was identified in all V. parahaemolyticus strains. Antimicrobial resistance revealed that most isolates were resistant to ampicillin (75.0%) and streptomycin (68.3%), whereas all strains were sensitive to chloramphenicol and tetracyclines. Most of the isolates (55.8%) showed resistance against two classes of antimicrobials, mainly to ampicillin and streptomycin (46.2%). Only one isolate displayed a multiresistant pattern. Genotypic analysis of V. parahaemolyticus revealed a high degree of diversity among the isolates tested. The pulsed-field gel electrophoresis (PFGE) method distinguished 73 clonal groups, and the most numerous group consisted of 7 strains. Sequencing by the multilocus sequence typing (MLST) method showed 76 sequence types (STs), of which ST481 and ST1361 were most frequently identified. In addition, 51 (67.1%) new sequence types were discovered and added to the PubMLST international database.IMPORTANCE The presence of V. parahaemolyticus in seafood may pose a risk for consumers, especially in countries where shellfish are eaten raw. In recent years, a significant increase of food poisoning caused by these bacteria has been also observed in Europe. Our results highlight the high level of V. parahaemolyticus contamination of seafood, along with the isolates being potentially pathogenic for humans. However, the first-line antimicrobials, such as tetracyclines and fluoroquinolones, remained highly effective against V. parahaemolyticus The monitoring of antimicrobial resistance of isolates is important to ensure the high efficacy in the treatment of human infections. Most of V. parahaemolyticus strains possessed new sequence types (STs), which showed the high genetic diversity of the isolates tested.

Microevolution of Streptococcus agalactiae ST-261 from Australia Indicates Dissemination via Imported Tilapia and Ongoing Adaptation to Marine Hosts or Environment

Streptococcus agalactiae (group B Streptococcus [GBS]) causes disease in a wide range of animals. The serotype Ib lineage is highly adapted to aquatic hosts, exhibiting substantial genome reduction compared with terrestrial conspecifics. Here, we sequence genomes from 40 GBS isolates, including 25 isolates from wild fish and captive stingrays in Australia, six local veterinary or human clinical isolates, and nine isolates from farmed tilapia in Honduras, and compared them with 42 genomes from public databases. Phylogenetic analysis based on nonrecombinant core-genome single nucleotide polymorphisms (SNPs) indicated that aquatic serotype Ib isolates from Queensland were distantly related to local veterinary and human clinical isolates. In contrast, Australian aquatic isolates are most closely related to a tilapia isolate from Israel, differing by only 63 core-genome SNPs. A consensus minimum spanning tree based on core-genome SNPs indicates the dissemination of sequence type 261 (ST-261) from an ancestral tilapia strain, which is congruent with several introductions of tilapia into Australia from Israel during the 1970s and 1980s. Pangenome analysis identified 1,440 genes as core, with the majority being dispensable or strain specific, with non-protein-coding intergenic regions (IGRs) divided among core and strain-specific genes. Aquatic serotype Ib strains have lost many virulence factors during adaptation, but six adhesins were well conserved across the aquatic isolates and might be critical for virulence in fish and for targets in vaccine development. The close relationship among recent ST-261 isolates from Ghana, the United States, and China with the Israeli tilapia isolate from 1988 implicates the global trade in tilapia seed for aquaculture in the widespread dissemination of serotype Ib fish-adapted GBS.IMPORTANCEStreptococcus agalactiae (GBS) is a significant pathogen of humans and animals. Some lineages have become adapted to particular hosts, and serotype Ib is highly specialized to fish. Here, we show that this lineage is likely to have been distributed widely by the global trade in tilapia for aquaculture, with probable introduction into Australia in the 1970s and subsequent dissemination in wild fish populations. We report here the variability in the polysaccharide capsule among this lineage but identify a cohort of common surface proteins that may be a focus of future vaccine development to reduce the biosecurity risk in international fish trade.

Cryptosporidium genotyping in Europe: The current status and processes for a harmonised multi-locus genotyping scheme

Due to the occurrence of genetic recombination, a reliable and discriminatory method to genotype Cryptosporidium isolates at the intra-species level requires the analysis of multiple loci, but a standardised scheme is not currently available. A workshop was held at the Robert Koch Institute, Berlin in 2016 that gathered 23 scientists with appropriate expertise (in either Cryptosporidium genotyping and/or surveillance, epidemiology or outbreaks) to discuss the processes for the development of a robust, standardised, multi-locus genotyping (MLG) scheme and propose an approach. The background evidence and main conclusions were outlined in a previously published report; the objectives of this further report are to describe 1) the current use of Cryptosporidium genotyping, 2) the elicitation and synthesis of the participants' opinions, and 3) the agreed processes and criteria for the development, evaluation and validation of a standardised MLG scheme for Cryptosporidium surveillance and outbreak investigations. Cryptosporidium was characterised to the species level in 7/12 (58%) participating European countries, mostly for human outbreak investigations. Further genotyping was mostly by sequencing the gp60 gene. A ranking exercise of performance and convenience criteria found that portability, biological robustness, typeability, and discriminatory power were considered by participants as the most important attributes in developing a multilocus scheme. The major barrier to implementation was lack of funding. A structured process for marker identification, evaluation, validation, implementation, and maintenance was proposed and outlined for application to Cryptosporidium, with prioritisation of Cryptosporidium parvum to support investigation of transmission in Europe.

A complex scenario of tuberculosis transmission is revealed through genetic and epidemiological surveys in Porto

Background

Tuberculosis (TB) incidence is decreasing worldwide and eradication is becoming plausible. In low-incidence countries, intervention on migrant populations is considered one of the most important strategies for elimination. However, such measures are inappropriate in European areas where TB is largely endemic, such as Porto in Portugal. We aim to understand transmission chains in Porto through a genetic characterization of Mycobacterium tuberculosis strains and through a detailed epidemiological evaluation of cases.

Methods

We genotyped the M. tuberculosis strains using the MIRU-VNTR system. We performed an evolutionary reconstruction of the genotypes with median networks, used in this context for the first time. TB cases from a period of two years were evaluated combining genetic, epidemiological and georeferencing information.

Results

The data reveal a unique complex scenario in Porto where the autochthonous population acts as a genetic reservoir of M. tuberculosis diversity with discreet episodes of transmission, mostly undetected using classical epidemiology alone.

Conclusions

Although control policies have been successful in decreasing incidence in Porto, the discerned complexity suggests that, for elimination to be a realistic goal, strategies need to be adjusted and coupled with a continuous genetic characterization of strains and detailed epidemiological evaluation, in order to successfully identify and interrupt transmission chains.

Shiga Toxin-Producing Escherichia coli O157 Shedding Dynamics in an Australian Beef Herd

Shiga toxin-producing Escherichia coli (STEC) O157 is an important foodborne pathogen that can be transmitted to humans both directly and indirectly from the feces of beef cattle, its primary reservoir. Numerous studies have investigated the shedding dynamics of E. coli O157 by beef cattle; however, the spatiotemporal trends of shedding are still not well understood. Molecular tools can increase the resolution through the use of strain typing to explore transmission dynamics within and between herds and identify strain-specific characteristics that may influence pathogenicity and spread. Previously, the shedding dynamics and molecular diversity, through the use of multilocus variable number of tandem repeat analysis (MLVA) of STEC O157, were separately investigated in an Australian beef herd over a 9-month study period. Variation in shedding was observed over time, and 33 MLVA types were identified. The study presented here combines the two datasets previously published with an aim to clarify the relationship between epidemiological variables and strain types. Three major genetic clusters (GCs) were identified that were significantly associated with the location of the cattle in different paddocks. No significant association between GCs and individual cow was observed. Results from this molecular epidemiological study provide evidence for herd-level clonal replacement over time that may have been triggered by movement to a new paddock. In conclusion, this study has provided further insight into STEC O157 shedding dynamics and pathogen transmission. Knowledge gaps remain regarding the relationship of strain types and the shedding dynamics of STEC O157 by beef cattle that could be further clarified through the use of whole-genome sequencing.

Improving the genotyping resolution of Cryptosporidium hominis subtype IbA10G2 using one step PCR-based amplicon sequencing

Cryptosporidium hominis gp60 subtype IbA10G2 is a common cause of cryptosporidiosis. This subtype is responsible for many waterborne outbreaks as well as sporadic cases and is considered virulent and highly important in the epidemiology of cryptosporidiosis. Due to low heterogeneity within the genome of C. hominis it has been difficult to identify epidemiological markers with higher resolution than gp60. However, new markers are required in order to improve outbreak investigations and studies of the transmission dynamics of this clinically important subtype. Based on the whole genome sequences of 17 C. hominis isolates, we have identified several differential loci and developed a new sequence based typing panel with higher resolution than gp60. An amplicon sequencing method was also developed which is based on a one-step PCR which can be sequenced using a Next Generation Sequencing (NGS) platform. Such a system provides a rapid and high-throughput workflow. A panel of nine loci with 10 single nucleotide variants (SNV) was selected and evaluated using clinical IbA10G2 isolates from sporadic, cluster and outbreak associated cases. The specimens were separated into 10 different genetic profiles named sequence types (STs). All isolates within an outbreak or cluster belonged to the same ST, including several samples from the two large waterborne outbreaks which occurred in Sweden between 2010 and 2011 indicating that these outbreaks might be linked. The results demonstrate the methods suitability for improved genotyping of C. hominis IbA10G2.

A Characterization of Minimum Spanning Tree-Like Metric Spaces

Recent years have witnessed a surge of biological interest in the minimum spanning tree (MST) problem for its relevance to automatic model construction using the distances between data points. Despite the increasing use of MST algorithms for this purpose, the goodness-of-fit of an MST to the data is often elusive because no quantitative criteria have been developed to measure it. Motivated by this, we provide a necessary and sufficient condition to ensure that a metric space on n points can be represented by a fully labeled tree on n vertices, and thereby determine when an MST preserves all pairwise distances between points in a finite metric space.

Molecular epidemiology of Corynebacterium pseudotuberculosis isolated from horses in California

Corynebacterium pseudotuberculosis biovar Equi is an important pathogen of horses. It is increasing in frequency in the United States, and is responsible for various clinical forms of infection, including external abscesses, internal abscesses of the abdominal or thoracic cavities, and ulcerative lymphangitis. The host/pathogen factors dictating the form or severity of infection are currently unknown. Our recent investigations have shown that genotyping C. pseudotuberculosis isolates using enterobacterial repetitive intergenic consensus (ERIC)-PCR is useful for understanding the evolutionary genetics of the species as well for molecular epidemiology studies. The aims of the present study were to assess (i) the genetic diversity of C. pseudotuberculosis strains isolated from horses in California, United States and (ii) the epidemiologic relationships among isolates. One hundred and seven C. pseudotuberculosis biovar Equi isolates from ninety-five horses, and two C. pseudotuberculosis biovar Ovis strains, C. pseudotuberculosis ATCC 19410^T type strain and C. pseudotuberculosis 1002 vaccine strain, were fingerprinted using the ERIC 1+2-PCR. C. pseudotuberculosis isolated from horses showed a high genetic diversity, clustering in twenty-seven genotypes with a diversity index of 0.91. Minimal spanning tree showed four major clonal complexes with a pattern of temporal clustering. Strains isolated from the same horse showed identical ERIC 1+2-PCR genotype, with the exception of two strains isolated from the same animal that showed distinct genotypes, suggesting a co-infection. We found no strong genetic signals related to clinical form (including internal versus external infections). However, temporal clustering of genotypes was observed.

Genetic Analysis of Phytophthora nicotianae Populations from Different Hosts Using Microsatellite Markers

In all, 231 isolates of Phytophthora nicotianae representing 14 populations from different host genera, including agricultural crops (Citrus, Nicotiana, and Lycopersicon), potted ornamental species in nurseries (Lavandula, Convolvulus, Myrtus, Correa, and Ruta), and other plant genera were characterized using simple-sequence repeat markers. In total, 99 multilocus genotypes (MLG) were identified, revealing a strong association between genetic grouping and host of recovery, with most MLG being associated with a single host genus. Significant differences in the structure of populations were revealed but clonality prevailed in all populations. Isolates from Citrus were found to be genetically related regardless of their geographic origin and were characterized by high genetic uniformity and high inbreeding coefficients. Higher variability was observed for other populations and a significant geographical structuring was determined for isolates from Nicotiana. Detected differences were related to the propagation and cultivation systems of different crops. Isolates obtained from Citrus spp. are more likely to be distributed worldwide with infected plant material whereas Nicotiana and Lycopersicon spp. are propagated by seed, which would not contribute to the spread of the pathogen and result in a greater chance for geographic isolation of lineages. With regard to ornamental species in nurseries, the high genetic variation is likely the result of the admixture of diverse pathogen genotypes through the trade of infected plant material from various geographic origins, the presence of several hosts in the same nursery, and genetic recombination through sexual reproduction of this heterothallic species.

Using genomics data to reconstruct transmission trees during disease outbreaks

Genetic sequence data from pathogens present a novel means to investigate the spread of infectious disease between infected hosts or infected premises, complementing traditional contact-tracing approaches, and much recent work has gone into developing methods for this purpose. The objective is to recover the epidemic transmission tree, which identifies who infected whom. This paper reviews the various approaches that have been taken. The first step is to define a measure of difference between sequences, which must be done while taking into account such factors as recombination and convergent evolution. Three broad categories of method exist, of increasing complexity: those that assume no withinhost genetic diversity or mutation, those that assume no within-host diversity but allow mutation, and those that allow both. Until recently, the assumption was usually made that every host in the epidemic could be identified, but this is now being relaxed, and some methods are intended for sparsely sampled data, concentrating on the identification of pairs of sequences that are likely to be the result of direct transmission rather than inferring the complete transmission tree. Many of the procedures described here are available to researchers as free software.

Corynebacterium pseudotuberculosis may be under anagenesis and biovar Equi forms biovar Ovis: a phylogenic inference from sequence and structural analysis

BACKGROUND

Corynebacterium pseudotuberculosis can be classified into two biovars or biovars based on their nitrate-reducing ability. Strains isolated from sheep and goats show negative nitrate reduction and are termed biovar Ovis, while strains from horse and cattle exhibit positive nitrate reduction and are called biovar Equi. However, molecular evidence has not been established so far to understand this difference, specifically if these C. pseudotuberculosis strains are under an evolutionary process.

RESULTS

The ERIC 1 + 2 Minimum-spanning tree from 367 strains of C. pseudotuberculosis showed that the great majority of biovar Ovis strains clustered together, but separately from biovar Equi strains that also clustered amongst themselves. Using evolutionarily conserved genes (rpoB, gapA, fusA, and rsmE) and their corresponding amino acid sequences, we analyzed the phylogenetic relationship among eighteen strains of C. pseudotuberculosis belonging to both biovars Ovis and Equi. Additionally, conserved point mutation based on structural variation analysis was also carried out to elucidate the genotype-phenotype correlations and speciation. We observed that the biovars are different at the molecular phylogenetic level and a probable anagenesis is occurring slowly within the species C. pseudotuberculosis.

CONCLUSIONS

Taken together the results suggest that biovar Equi is forming the biovar Ovis. However, additional analyses using other genes and other bacterial strains are required to further support our anagenesis hypothesis in C. pseudotuberculosis.

Limited Sexual Reproduction and Quick Turnover in the Population Genetic Structure of Phytophthora infestans in Fujian, China

The mating system plays an important role in the spatiotemporal dynamics of pathogen populations through both its direct and indirect impact on the generation and distribution of genetic variation. Here, we used a combination of microsatellite and phenotypic markers to investigate the spatiotemporal distribution of genetic variation in Phytophthora infestans isolates collected from Fujian, China and to determine the role of sexual reproduction in the dynamics. Although the pathogen populations in this region were dominated by self-fertile genotypes, sexual reproduction only occurred occasionally and its contributions to the population genetic structure of P. infestans and epidemics of late blight in the region were limited. Only 49 genotypes were detected among the 534 isolates assayed and the pathogen populations displayed significant heterozygosity excess. Hierarchical analysis revealed that 21.42% of genetic variation was attributed to the difference among sampling years while only 4.45% was attributed to the difference among locations, suggesting temporal factors play a more important role in the population genetic dynamics of P. infestans than spatial factors in this region. We propose that clonal reproduction, combined with founder effects and long distance dispersal of sporangia, is responsible for the observed pattern of spatiotemporal dynamics in P. infestans.

Not seeing the forest for the trees: size of the minimum spanning trees (MSTs) forest and branch significance in MST-based phylogenetic analysis

Trees, including minimum spanning trees (MSTs), are commonly used in phylogenetic studies. But, for the research community, it may be unclear that the presented tree is just a hypothesis, chosen from among many possible alternatives. In this scenario, it is important to quantify our confidence in both the trees and the branches/edges included in such trees. In this paper, we address this problem for MSTs by introducing a new edge betweenness metric for undirected and weighted graphs. This spanning edge betweenness metric is defined as the fraction of equivalent MSTs where a given edge is present. The metric provides a per edge statistic that is similar to that of the bootstrap approach frequently used in phylogenetics to support the grouping of taxa. We provide methods for the exact computation of this metric based on the well known Kirchhoff’s matrix tree theorem. Moreover, we implement and make available a module for the PHYLOViZ software and evaluate the proposed metric concerning both effectiveness and computational performance. Analysis of trees generated using multilocus sequence typing data (MLST) and the goeBURST algorithm revealed that the space of possible MSTs in real data sets is extremely large. Selection of the edge to be represented using bootstrap could lead to unreliable results since alternative edges are present in the same fraction of equivalent MSTs. The choice of the MST to be presented, results from criteria implemented in the algorithm that must be based in biologically plausible models.

Genetic diversity of Brucella ovis isolates from Rio Grande do Sul, Brazil, by MLVA16

Background

Ovine epididymitis is predominantly associated with Brucella ovis infection. Molecular characterization of Brucella spp. achieved by multi-locus variable number of tandem repeats (VNTR) analyses (MLVA) have proved to be a powerful tool for epidemiological trace-back studies. Thus, the aim of this study was to evaluate the genetic diversity of Brucella ovis isolates from Rio Grande do Sul State, Brazil, by MLVA16.

Findings

MLVA16 genotyping identified thirteen distinct genotypes and a Hunter-Gaston diversity index of 0.989 among the fourteen B. ovis genotyped strains. All B. ovis MLVA16 genotypes observed in the present study represented non-previously described profiles. Analyses of the eight conserved loci included in panel 1 (MLVA8) showed three different genotypes, two new and one already described for B. ovis isolates. Among ten B. ovis isolates from same herd only two strains had identical pattern, whereas the four isolates with no epidemiologic information exhibited a single MLVA16 pattern each. Analysis of minimal spanning tree, constructed using the fourteen B. ovis strains typed in this study together with all nineteen B. ovis MLVA16 genotypes available in the MLVAbank 2014, revealed the existence of two clearly distinct major clonal complexes.

Conclusions

In conclusion, the results of the present study showed a high genetic diversity among B. ovis field isolates from Rio Grande do Sul State, Brazil, by MLVA16.

Genetic relatedness of Brucella suis biovar 2 isolates from hares, wild boars and domestic pigs

Porcine brucellosis generally manifests as disorders in reproductive organs potentially leading to serious losses in the swine industry. Brucella suis biovar 2 is endemic in European wild boar (Sus scrofa) and hare (Lepus europeus, Lepus capensis) populations, thus these species may play a significant role in disease spread and serve as potential sources of infection for domestic pigs. The aim of this study was an epidemiologic analysis of porcine brucellosis in Hungary and a comparative analysis of B. suis bv. 2 strains from Europe using multiple-locus variable-number tandem repeat analysis (MLVA). MLVA-16 and its MLVA-11 subset were used to determine the genotypes of 68 B. suis bv. 2 isolates from Hungary and results were then compared to European MLVA genotypes. The analyses indicated relatively high genetic diversity of B. suis bv. 2 in Hungary. Strains isolated from hares and wild boars from Hungary showed substantial genetic divergence, suggesting separate lineages in each host and no instances of cross species infections. The closest relatives of strains from Hungarian wild boars and domestic pigs were mainly in the isolates from German and Croatian boars and pigs. The assessment of the European MLVA genotypes of wild boar isolates generally showed clustering based on geographic origin. The hare strains were relatively closely related to one another and did not cluster based on geographic origin. The limited relationships between geographic origin and genotype in isolates from hares might be the result of cross-border live animal translocation. The results could also suggest that certain B. suis strains are more adapted to hares. Across Europe, isolates from domestic pigs were closely related to isolates originating from both hares and wild boars, supporting the idea that wild animals are a source of brucellosis in domestic pigs.

Evaluation of ERIC-PCR as genotyping method for Corynebacterium pseudotuberculosis isolates

The aim of this study was to evaluate the Enterobacterial Repetitive Intergenic Consensus (ERIC-PCR) as a tool for molecular typing of C. pseudotuberculosis isolates from eight different hosts in twelve countries. Ninety-nine C. pseudotuberculosis field strains, one type strain (ATCC 19410^T) and one vaccine strain (1002) were fingerprinted using the ERIC-1R and ERIC-2 primers, and the ERIC-1R+ERIC-2 primer pair. Twenty-nine different genotypes were generated by ERIC 1-PCR, 28 by ERIC 2-PCR and 35 by ERIC 1+2-PCR. The discriminatory index calculated for ERIC 1, ERIC 2, and ERIC 1+2-PCR was 0.89, 0.86, and 0.92, respectively. Epidemiological concordance was established for all ERIC-PCR assays. ERIC 1+2-PCR was defined as the best method based on suitability of the amplification patterns and discriminatory index. Minimal spanning tree for ERIC 1+2-PCR revealed three major clonal complexes and clustering around nitrate-positive (biovar Equi) and nitrate-negative (biovar Ovis) strains. Therefore, ERIC 1+2-PCR proved to be the best technique evaluated in this study for genotyping C. pseudotuberculosis strains, due to its usefulness for molecular epidemiology investigations.

Population genetic structure of the Lyme disease vector Ixodes scapularis at an apparent spatial expansion front

Modeling and empirical evidence suggests that Lyme disease is undergoing geographic expansion from principal foci in the midwestern and northeastern United States. Virginia is at the southern edge of the current expansion zone and has seen dramatic rise in human Lyme disease cases since 2007, potentially owing to a recent increase in vector abundance. Ixodes scapularis is known throughout the eastern US but behavioral or physiological variation between northern and southern lineages might lead northern-variant ticks to more frequently parasitize humans. We hypothesized that recent spatial and numerical increase in Lyme disease cases is associated with demographic and/or spatial expansion of I. scapularis and that signals of these phenomena would be detectable and discernable in population genetic signals. In summer and fall 2011, we collected nymphal I. scapularis by drag sampling and adult I. scapularis from deer carcasses at hunting check stations at nine sites arranged along an east-west transect through central Virginia. We analyzed 16S mtDNA sequences data from up to 24 I. scapularis individuals collected from each site and detected a total of 24 haplotypes containing 29 segregating sites. We found no evidence for population genetic structure among these sites but we did find strong signals of both demographic and spatial expansion throughout our study system. We found two haplotypes (one individual each) representing a lineage of ticks that is only found in the southeastern United States, with the remaining individuals representing a less genetically diverse clade that is typical of the northern United States, but that has also been detected in the American South. Taken together, these results lead us to conclude that I. scapularis populations in Virginia are expanding and that this expansion may account for recent observed increases in Lyme disease.

Multilocus sequence typing (MLST) of Staphylococcus aureus

MLST is a widely accepted method of sequence-based typing that relies on analysis of relatively conserved genes that encode essential proteins. For Staphylococcus aureus the level of discrimination provided by MLST is sufficient to provide a relatively detailed picture of the global dissemination of the pathogen. The method is not restrictive in the precise methodology used to acquire the sequences, but the method of assigning types requires that the data be of high quality. Excellent web-based tools have been developed and are curated by the groups that launched MLST. These tools have allowed the scheme to be maintained as a coherent global asset and assist users in the analysis of their data.

Clustering acinetobacter strains by optical mapping

Optical mapping is a technique that produces an ordered restriction map of a bacterial or eukaryotic chromosome. We have developed a new method, the BOP method, to compare experimental optical maps with in silico optical maps of complete genomes to infer the presence/absence of short DNA sequences (bops) in each genome. The BOP method, as implemented by the Optical Mapping suite of four programs, circumvents the necessity of whole-genome multiple alignments and permits reliable strain typing and clustering on the basis of optical maps. We have applied the Optical Mapping Suite to 125 strains of Acinetobacter sp., including 11 completely sequenced genomes and 114 Acinetobacter complex from three US military hospitals. We found that optical mapping completely resolves all 125 strains. Signal to noise analysis showed that when the 125 strains were considered together almost 1/3 of the experimental fragments were misidentified. We found that the set of 125 genomes could be divided into three clusters, two of which included sequenced genomes. Signal to noise analysis after clustering showed that only 3.5% of the experimental restriction fragments were misidentified. Minimum spanning trees of the two clusters that included sequenced genomes are presented. The programs we have developed provide a more rigorous approach for analyzing optical map data than previously existed.

When whole-genome alignments just won't work: kSNP v2 software for alignment-free SNP discovery and phylogenetics of hundreds of microbial genomes

Effective use of rapid and inexpensive whole genome sequencing for microbes requires fast, memory efficient bioinformatics tools for sequence comparison. The kSNP v2 software finds single nucleotide polymorphisms (SNPs) in whole genome data. kSNP v2 has numerous improvements over kSNP v1 including SNP gene annotation; better scaling for draft genomes available as assembled contigs or raw, unassembled reads; a tool to identify the optimal value of k; distribution of packages of executables for Linux and Mac OS X for ease of installation and user-friendly use; and a detailed User Guide. SNP discovery is based on k-mer analysis, and requires no multiple sequence alignment or the selection of a single reference genome. Most target sets with hundreds of genomes complete in minutes to hours. SNP phylogenies are built by maximum likelihood, parsimony, and distance, based on all SNPs, only core SNPs, or SNPs present in some intermediate user-specified fraction of targets. The SNP-based trees that result are consistent with known taxonomy. kSNP v2 can handle many gigabases of sequence in a single run, and if one or more annotated genomes are included in the target set, SNPs are annotated with protein coding and other information (UTRs, etc.) from Genbank file(s). We demonstrate application of kSNP v2 on sets of viral and bacterial genomes, and discuss in detail analysis of a set of 68 finished E. coli and Shigella genomes and a set of the same genomes to which have been added 47 assemblies and four “raw read” genomes of H104:H4 strains from the recent European E. coli outbreak that resulted in both bloody diarrhea and hemolytic uremic syndrome (HUS), and caused at least 50 deaths.

Pathotype and genetic diversity amongst Indian isolates of Xanthomonas oryzae pv. oryzae

A number of rice resistance genes, called Xa genes, have been identified that confer resistance against various strains of Xanthomonas oryzae pv. oryzae (Xoo), the causal agent of bacterial blight. An understanding of pathotype diversity within the target pathogen population is required for identifying the Xa genes that are to be deployed for development of resistant rice cultivars. Among 1024 isolates of Xoo collected from 20 different states of India, 11 major pathotypes were distinguished based on their reaction towards ten Xa genes (Xa1, Xa3, Xa4, xa5, Xa7, xa8, Xa10, Xa11, xa13, Xa21). Isolates belonging to pathotype III showing incompatible interaction towards xa8, xa13 and Xa21 and compatible interaction towards the rest of Xa genes formed the most frequent (41%) and widely distributed pathotype. The vast majority of the assayed Xoo isolates were incompatible with one or more Xa genes. Exceptionally, the isolates of pathotype XI were virulent on all Xa genes, but have restricted distribution. Considering the individual R-genes, Xa21 appeared as the most broadly effective, conferring resistance against 88 % of the isolates, followed in decreasing order by xa13 (84 %), xa8 (64 %), xa5 (30 %), Xa7 (17 %) and Xa4 (14 %). Fifty isolates representing all the eleven pathotypes were analyzed by southern hybridization to determine their genetic relatedness using the IS1112 repeat element of Xoo. Isolates belonging to pathotype XI were the most divergent. The results suggest that one RFLP haplotype that is widely distributed all over India and is represented in strains from five different pathotypes might be an ancestral haplotype. A rice line with xa5, xa13 and Xa21 resistance genes is resistant to all strains, including those belonging to pathotype XI. This three gene combination appears to be the most suitable Xa gene combination to be deployed in Indian rice cultivars.

Using complete genome comparisons to identify sequences whose presence accurately predicts clinically important phenotypes

In clinical settings it is often important to know not just the identity of a microorganism, but also the danger posed by that particular strain. For instance, Escherichia coli can range from being a harmless commensal to being a very dangerous enterohemorrhagic (EHEC) strain. Determining pathogenic phenotypes can be both time consuming and expensive. Here we propose a simple, rapid, and inexpensive method of predicting pathogenic phenotypes on the basis of the presence or absence of short homologous DNA segments in an isolate. Our method compares completely sequenced genomes without the necessity of genome alignments in order to identify the presence or absence of the segments to produce an automatic alignment of the binary string that describes each genome. Analysis of the segment alignment allows identification of those segments whose presence strongly predicts a phenotype. Clinical application of the method requires nothing more that PCR amplification of each of the set of predictive segments. Here we apply the method to identifying EHEC strains of E. coli and to distinguishing E. coli from Shigella. We show in silico that with as few as 8 predictive sequences, if even three of those predictive sequences are amplified the probability of being EHEC or Shigella is >0.99. The method is thus very robust to the occasional amplification failure for spurious reasons. Experimentally, we apply the method to screening a set of 98 isolates to distinguishing E. coli from Shigella, and EHEC from non-EHEC E. coli strains and show that all isolates are correctly identified.

Investigation of the population structure of the tick vector of Lyme disease Ixodes scapularis (Acari: Ixodidae) in Canada using mitochondrial cytochrome C oxidase subunit I gene sequences

Genotyping of Ixodes scapularis Say (Acari: Ixodidae) ticks could enhance understanding of the occurrence and genotypes of I. scapularis-borne pathogens. We investigated the utility of mitochondrial (mt) Cytochrome C Oxidase subunit I gene (cox1) sequences as a tool for understanding the population structure of I. scapularis collected in Canada, where we also investigated the geographic occurrence of different cox1 haplotypes. Sequences obtained from 414 ticks were one of 55 unique haplotypes, most of which grouped into one of six clades. Demographic analysis suggested that cox1 sequences have haplotype and nucleotide diversity comparable to other mt genes. All haplotypes were connected in a single minimum spanning network tree. Despite low fixation index values there were significant differences in the frequency of occurrence of haplotypes of different clades among four geographic regions: 1) Alberta to western Ontario, 2) eastern Ontario, 3) Quebec, and 4) Atlantic Provinces; suggesting that cox1 sequences could reveal population structure differences between I. scapularis in geographically separated populations of northeastern and midwestern North America. Spatial clusters of ticks of the same haplotype identified in regions of southern Quebec and southern Ontario where I. scapularis is invading were consistent with population bottlenecks associated with founder events. These findings suggest that cox1 sequences are useful for the study of I. scapularis population structure, are of sufficient diversity that spatial analyses of haplotypes can be used to identify where I. scapularis is emerging in southern Canada, and may be useful for exploring differences between northeastern and midwestern populations of I. scapularis.