High-resolution melting genotyping of Enterococcus faecium based on multilocus sequence typing derived single nucleotide polymorphisms

We have developed a single nucleotide polymorphism (SNP) nucleated high-resolution melting (HRM) technique to genotype Enterococcus faecium. Eight SNPs were derived from the E. faecium multilocus sequence typing (MLST) database and amplified fragments containing these SNPs were interrogated by HRM. We tested the HRM genotyping scheme on 85 E. faecium bloodstream isolates and compared the results with MLST, pulsed-field gel electrophoresis (PFGE) and an allele specific real-time PCR (AS kinetic PCR) SNP typing method. In silico analysis based on predicted HRM curves according to the G+C content of each fragment for all 567 sequence types (STs) in the MLST database together with empiric data from the 85 isolates demonstrated that HRM analysis resolves E. faecium into 231 "melting types" (MelTs) and provides a Simpson's Index of Diversity (D) of 0.991 with respect to MLST. This is a significant improvement on the AS kinetic PCR SNP typing scheme that resolves 61 SNP types with D of 0.95. The MelTs were concordant with the known ST of the isolates. For the 85 isolates, there were 13 PFGE patterns, 17 STs, 14 MelTs and eight SNP types. There was excellent concordance between PFGE, MLST and MelTs with Adjusted Rand Indices of PFGE to MelT 0.936 and ST to MelT 0.973. In conclusion, this HRM based method appears rapid and reproducible. The results are concordant with MLST and the MLST based population structure.     

Molecular characterization of Enterococcus faecium isolated from hospitalized patients in Korea

AIMS: Multilocus sequence typing (MLST) was performed for vancomycin-resistant Enterococcus faecium (VREF) from diverse geographical areas in Korea to obtain insights into the genetic relationships with other molecular profiles. To understand the diversity of lineages, vancomycin-susceptible E. faecium (VSEF) were included. METHODS AND RESULTS: A total of 60 E. faecium isolates were analysed by MLST and esp profile. Molecular typing of Tn1546 of 30 VREF strains was evaluated by overlapping PCR of Tn1546 and DNA sequencing. Seven sequence types (ST) were found among 30 VSEF isolates, and four STs were found among 30 VREF isolates. The types most frequently encountered were ST 78 (26 isolates) and ST 203 (16 isolates). Of the 60 E. faecium isolates, 35 isolates were positive for the esp gene. On molecular typing of Tn1546, all VREF isolates were divided into four main types. Strains with the same ST showed divergence in Tn1546 types and strains with the same Tn1546 type represented different STs. CONCLUSIONS: An association between Tn1546 typing and MLST was not found. SIGNIFICANCE AND IMPACT OF THE STUDY: These results suggest that the horizontal spread of Tn1546 between strains plays a major role in the dissemination of vancomycin resistance in Korea.     

Multilocus Sequence Typing (MLST) for Characterization of Enterobacter cloacae

Enterobacter cloacae is an important emerging pathogen, which sometime causes respiratory infection, surgical site infection, urinary infection, sepsis, and outbreaks at neonatal units. We have developed a multilocus sequence typing (MLST) scheme utilizing seven housekeeping genes and evaluated the performance in 101 clinical isolates. The MLST scheme yielded 83 sequence types (ST) including 78 novel STs found in the clinical isolates. These findings supported the robustness of the MLST scheme developed in this study.     

Multilocus sequence typing reveals high genetic diversity and epidemic population structure for the fish pathogen Yersinia ruckeri

Yersinia ruckeri is the causative agent of enteric redmouth in fish and one of the major bacterial pathogens causing losses in salmonid aquaculture. Previously typing methods, including restriction enzyme analysis, pulsed-field gel electrophoresis and multilocus enzyme electrophoresis (MLEE) have indicated a clonal population structure. In this work, we describe a multilocus sequence typing (MLST) scheme for Y.ruckeri based on the internal fragment sequence of six housekeeping genes. This MLST scheme was applied to 103 Y.ruckeri strains from diverse geographic areas and hosts as well as environmental sources. Sequences obtained from this work were deposited and are available in a public database (http://publmst.org/yruckeri/). Thirty different sequence types (ST) were identified, 21 of which were represented by a single isolate, evidencing high genetic diversity. ST2 comprised more than one-third of the isolates and was most frequently observed among isolates from trout. Two major clonal complexes (CC) were identified by eBURST analysis showing a common evolutionary origin for 94 isolates forming 21 STs into CC1 and for 6 isolates of 6 STs in the CC2. It was also possible to associate some unique ST with isolates from recent outbreaks in vaccinated salmonid fish.     

MLST clustering of Campylobacter jejuni isolates from patients with gastroenteritis,reactive arthritis and Guillain–Barré syndrome

Aims: To determine the diversity and population structure of Campylobacter jejuni (C. jejuni) isolates from Danish patients and to examine the association between multilocus sequence typing types and different clinical symptoms including gastroenteritis (GI), Guillain–Barré syndrome (GBS) and reactive arthritis (RA). Methods and Results: Multilocus sequence typing (MLST) was used to characterize 122 isolates, including 18 from patients with RA and 8 from patients with GBS. The GI and RA isolates were collected in Denmark during 2002–2003 and the GBS isolates were obtained from other countries. In overall, 51 sequence types (STs) were identified within 18 clonal complexes (CCs). Of these three CCs, ST‐21, ST‐45 and ST‐22 clonal complexes accounted for 64 percent of all isolates. The GBS isolates in this study significantly grouped into the ST‐22 clonal complex, consistent with the PubMLST database isolates. There was no significant clustering of the RA isolates. Conclusions: Isolates from Denmark were found to be highly genetically diverse. GBS isolates grouped significantly with clonal complex ST‐22, but the absence of clustering of RA isolates indicated that the phylogenetic background for this sequela could not be reconstructed using variation in MLST loci. Possibly, putative RA‐associated genes may vary, by recombination or expression differences, independent of MLST loci. Significance and Impact of the Study: MLST typing of C. jejuni isolates from Danish patients with gastroenteritis confirmed that the diversity of clones in Denmark is comparable to that in other European countries. Furthermore, a verification of the grouping of GBS isolates compared to RA isolates provides information about evolution of the bacterial population resulting in this important sequela.     

Investigating Burkholderia cepacia complex populations recovered from Italian maize rhizosphere by multilocus sequence typing

The Burkholderia cepacia complex (BCC) comprises at least nine closely related species of abundant environmental microorganisms. Some of these species are highly spread in the rhizosphere of several crop plants, particularly of maize; additionally, as opportunistic pathogens, strains of the BCC are capable of colonizing humans. We have developed and validated a multilocus sequence typing (MLST) scheme for the BCC. Although widely applied to understand the epidemiology of bacterial pathogens, MLST has seen limited application to the population analysis of species residing in the natural environment; we describe its novel application to BCC populations within maize rhizospheres. 115 BCC isolates were recovered from the roots of different maize cultivars from three different Italian regions over a 9-year period (1994-2002). A total of 44 sequence types (STs) were found of which 41 were novel when compared with existing MLST data which encompassed a global database of 1000 clinical and environmental strains representing nearly 400 STs. In this study of rhizosphere isolates approximately 2.5 isolates per ST was found, comparable to that found for the whole BCC population. Multilocus sequence typing also resolved inaccuracies associated with previous identification of the maize isolates based on recA gene restriction fragment length polymorphims and species-specific polymerase chain reaction. The 115 maize isolates comprised the following BCC species groups, B. ambifaria (39%), BCC6 (29%), BCC5 (10%), B. pyrrocinia (8%), B. cenocepacia IIIB (7%) and B. cepacia (6%), with BCC5 and BCC6 potentially constituting novel species groups within the complex. Closely related clonal complexes of strains were identified within B. cepacia, B. cenocepacia IIIB, BCC5 and BCC6, with one of the BCC5 clonal complexes being distributed across all three sampling sites. Overall, our analysis demonstrates that the maize rhizosphere harbours a massive diversity of novel BCC STs, so that their addition to our global MLST database increased the ST diversity by 10%.     

Molecular Characterization of Streptococcus agalactiae Isolated from Bovine Mastitis in Eastern China

One hundred and two Streptococcus agalactiae (group B streptococcus [GBS]) isolates were collected from dairy cattle with subclinical mastitis in Eastern China during 2011. Clonal groups were established by multilocus sequence typing (MLST) and pulsed-field gel electrophoresis (PFGE), respectively. Capsular polysaccharides (CPS), pilus and alpha-like-protein (Alp) family genes were also characterized by molecular techniques. MLST analysis revealed that these isolates were limited to three clonal groups and were clustered in six different lineages, i.e. ST (sequence type) 103, ST568, ST67, ST301, ST313 and ST570, of which ST568 and ST570 were new genotypes. PFGE analysis revealed this isolates were clustered in 27 PFGE types, of which, types 7, 8, 14, 15, 16, 18, 23 and 25 were the eight major types, comprising close to 70% (71/102) of all the isolates. The most prevalent sequence types were ST103 (58% isolates) and ST568 (31% isolates), comprising capsular genotype Ia isolates without any of the detected Alp genes, suggesting the appearance of novel genomic backgrounds of prevalent strains of bovine S. agalactiae. All the strains possessed the pilus island 2b (PI-2b) gene and the prevalent capsular genotypes were types Ia (89% isolates) and II (11% isolates), the conserved pilus type providing suitable data for the development of vaccines against mastitis caused by S. agalactiae.     

Multilocus Sequence Typing of Uropathogenic ESBL-Producing <Emphasis Type="Italic">Escherichia coli</Emphasis> Isolated in a Brazilian Community

In the present study, multilocus sequence typing (MLST), pulsed-field gel electrophoresis (PFGE) and polymerase chain reaction detection of three resistance genes were combined to characterize seven uropathogenic E. coli isolated from outpatients. Selected portions of seven housekeeping and three antibiotic-resistance genes of the isolates were sequenced. The seven isolates were classified into four different sequence types (STs) by MLST and five PGFE types. Three isolates had a novel allelic profile representing a new ST designated as ST528 and showed the same PFGE and resistance genes. Two isolates, both characterized as ST359, were differentiated by PFGE and shared only one of the antibiotic-resistance genes studied. Comparison of MLST results with those of PFGE and resistance genes demonstrated that Escherichia coli had acquired different antibiotic-resistance genes and DNA rearrangements, causing alterations in PFGE patterns but maintaining the same ST. Furthermore, this article also reports the first detection of a CTX-M-2 ESBL E. coli and SHV-5 in a Brazilian community.     

Core genome conservation of Staphylococcus haemolyticus limits sequence based population structure analysis

The notoriously multi-resistant Staphylococcus haemolyticus is an emerging pathogen causing serious infections in immunocompromised patients. Defining the population structure is important to detect outbreaks and spread of antimicrobial resistant clones. Currently, the standard typing technique is pulsed-field gel electrophoresis (PFGE). In this study we describe novel molecular typing schemes for S. haemolyticus using multi locus sequence typing (MLST) and multi locus variable number of tandem repeats (VNTR) analysis. Seven housekeeping genes (MLST) and five VNTR loci (MLVF) were selected for the novel typing schemes. A panel of 45 human and veterinary S. haemolyticus isolates was investigated. The collection had diverse PFGE patterns (38 PFGE types) and was sampled over a 20 year-period from eight countries. MLST resolved 17 sequence types (Simpsons index of diversity [SID]=0.877) and MLVF resolved 14 repeat types (SID=0.831). We found a low sequence diversity. Phylogenetic analysis clustered the isolates in three (MLST) and one (MLVF) clonal complexes, respectively. Taken together, neither the MLST nor the MLVF scheme was suitable to resolve the population structure of this S. haemolyticus collection. Future MLVF and MLST schemes will benefit from addition of more variable core genome sequences identified by comparing different fully sequenced S. haemolyticus genomes.     

Multi-Locus Sequence Typing of a Geographically and Temporally Diverse Sample of the Highly Clonal Human Pathogen Bartonella quintana

Bartonella quintana is a re-emerging pathogen and the causative agent of a variety of disease manifestations in humans including trench fever. Various typing methods have been developed for B. quintana, but these tend to be limited by poor resolution and, in the case of gel-based methods, a lack of portability. Multilocus sequence typing (MLST) has been used to study the molecular epidemiology of a large number of pathogens, including B. henselae, a close relative of B. quintana. We developed a MLST scheme for B. quintana based on the 7 MLST loci employed for B. henselae with two additional loci to cover underrepresented regions of the B. quintana chromosome. A total of 16 B. quintana isolates spanning over 60 years and three continents were characterized. Allelic variation was detected in five of the nine loci. Although only 8/4270 (0.002%) of the nucleotide sites examined were variable over all loci, these polymorphisms resolved the 16 isolates into seven sequence types (STs). We also demonstrate that MLST can be applied on uncultured isolates by direct PCR from cardiac valve tissue, and suggest this method presents a promising approach for epidemiological studies in this highly clonal organism. Phylogenetic and clustering analyses suggest that two of the seven STs form a distinct lineage within the population.     

Polyclonal non multiresistant Staphylococcus aureus isolates from clinical cases of infection occurring in Palermo, Italy, during a one-year surveillance period

ABSTRACT: BACKGROUND: The evolving epidemiology of methicillin resistant Staphylococcus aureus (MRSA) is characterized by the emergence of infections caused by non multiresistant MRSA carrying staphylococcal chromosomal cassette (SCC)mec IV or V in the healthcare settings. A molecular epidemiological analysis of non multiresistant MRSA isolates from four acute general hospitals was performed in Palermo, Italy, during a one year period. METHODS: For the purpose of the study, MRSA isolates were defined as non multiresistant when they were susceptible to at least three classes of non beta-lactam antibiotics. Seventy-five isolates were submitted to antimicrobial susceptibility testing, multilocus sequence typing (MLST) and polymerase chain reaction (PCR) for SCCmec, accessory gene regulator (agr) groups, arginine catabolic mobile element (ACME) and Panton Valentine leukocidin (PVL) toxin genes. For epidemiological typing, Multiple-Locus Variable-Number Tandem Repeat Fingerprinting (MLVF) was performed on all isolates and pulsed field gel electrophoresis (PFGE) on ST8 isolates. RESULTS: Non multiresistant MRSA isolates were isolated from all hospitals. Resistances to ciprofloxacin, macrolides and tetracycline were the most prevalent. MLST attributed 46 isolates with ST22, 13 with ST8, eight with ST1, three with ST50 and three with ST398. SCCmec type IV was found in all isolates. PVL was detected in one ST22 isolate. All isolates tested negative for the ACME element. MLVF identified 31 different patterns, some subtype clusters ranging in size between two and 22 isolates. The closely related PFGE patterns of the ST8 isolates differed from USA300. CONCLUSIONS: A polyclonal circulation of non multiresistant MRSA along with blurring of boundaries between healthcare associated (HA)-MRSA and community associated (CA)-MRSA appear to be occurring in our epidemiological setting. A better understanding of spread of MRSA with the support of molecular typing can provide invaluable information in the epidemiological, microbiological and clinical fields.     

Multilocus sequence typing and antimicrobial resistance in Enterococcus faecium isolates from fresh produce

The purpose of the present study was to determine the relatedness of Enterococcus faecium isolates from fresh produce to E. faecium strains from other sources by using multi-locus sequence typing (MLST) and to determine the antimicrobial resistance of the isolates. MLST analysis of 22 E. faecium isolates from fresh produce revealed 7 different sequence types (ST 22, ST 26, ST 43, ST 46, ST 55, ST 94 and ST 296). Most isolates belonged to ST 296 (40.9 %), followed by ST 94 (27.3 %). All isolates were sensitive to vancomycin and to imipenem, and only one was resistant to ampicillin (MIC 32 mg/l). However, all were resistant to cefotaxime and ceftazidine. E. faecium isolates from fresh produce were inhibited by quaternary compounds (benzalkonium chloride, cetrimide, hexadecylpyridinium chloride, didecyldimethylammonium bromide), biguanides (chlorhexidine), polyguanides [poly-(hexamethylene guanidinium) hydrochloride], bisphenols (triclosan, hexachlorophene) and biocidal solutions of P3 oxonia and P3 topax 66. Didecyldimethylammonium bromide and triclosan were the least effective biocides in growth inhibition, while hexadecylpyridinium chloride was the most effective. Results from MLST typing and antibiotic resistance suggest that the studied E. faecium isolates from fresh produce are not related to the clinically-relevant clonal complex CC17.     

The 75-kilodalton antigen of Bartonella bacilliformis is a structural homolog of the cell division protein FtsZ.

A genomic library of Bartonella bacilliformis was constructed and screened with human anti-Bartonella serum from a patient with the chronic, verruga peruana phase of bartonellosis. An immunoreactive clone isolated from this library was found to code for a 591-amino-acid protein with a high degree of sequence similarity to the FtsZ family of proteins. The degree of amino acid identity between the B. bacilliformis protein (FtsZ[Bb]) and the other FtsZ proteins is especially pronounced over the N-terminal 321 amino acids (N-terminal domain) of the sequence, with values ranging from 45% identity for the homolog from Micrococcus luteus (FtsZ[Ml]) to 91% identity for the homolog from Rhizobium melliloti, (FtsZ[Rm1]). All of the functional domains required for FtsZ activity are conserved in FtsZ(Bb) and are located within the N-terminal domain of the protein. FtsZ(Bb) is approximately twice as large as most of the other FtsZ proteins previously reported, a property it shares with FtsZ(Rm1). Like the Rhizobium homolog, FtsZ(Bb) has a C-terminal region of approximately 256 amino acids that is absent in the other FtsZ proteins. Evidence is presented that implicates this region in the protein's antigenicity and suggests that, unlike most other FtsZ homologs, FtsZ(Bb) is at least partly exposed at the cell surface. PCR analysis revealed that an ftsZ gene similar in size to the B. bacilliformis gene is present in Bartonella henselae, a bacterium that is closely related to B. bacilliformis.     

Multi-locus sequence typing of Bartonella henselae isolates from three continents reveals hypervirulent and feline-associated clones

Bartonella henselae is a zoonotic pathogen and the causative agent of cat scratch disease and a variety of other disease manifestations in humans. Previous investigations have suggested that a limited subset of B. henselae isolates may be associated with human disease. In the present study, 182 human and feline B. henselae isolates from Europe, North America and Australia were analysed by multi-locus sequence typing (MLST) to detect any associations between sequence type (ST), host species and geographical distribution of the isolates. A total of 14 sequence types were detected, but over 66% (16/24) of the isolates recovered from human disease corresponded to a single genotype, ST1, and this type was detected in all three continents. In contrast, 27.2% (43/158) of the feline isolates corresponded to ST7, but this ST was not recovered from humans and was restricted to Europe. The difference in host association of STs 1 (human) and 7 (feline) was statistically significant (P< or =0.001). eBURST analysis assigned the 14 STs to three clonal lineages, which contained two or more STs, and a singleton comprising ST7. These groups were broadly consistent with a neighbour-joining tree, although splits decomposition analysis was indicative of a history of recombination. These data indicate that B. henselae lineages differ in their virulence properties for humans and contribute to a better understanding of the population structure of B. henselae.     

Nonencapsulated or nontypeable Haemophilus influenzae are more likely than their encapsulated or serotypeable counterparts to have mutations in their fucose operon

Population biology of Haemophilus influenzae can be studied by multilocus sequence typing (MLST), and isolates are assigned sequence types (STs) based on nucleotide sequence variations in seven housekeeping genes, including fucK. However, the ST cannot be assigned if one of the housekeeping genes is absent or cannot be detected by the current protocol. Occasionally, strains of H. influenzae have been reported to lack the fucK gene. In this study, we examined the prevalence of this mutation among our collection of H. influenzae isolates. Of the 704 isolates studied, including 282 encapsulated and 422 nonencapsulated isolates, nine were not typeable by MLST owing to failure to detect the fucK gene. All nine fucK-negative isolates were nonencapsulated and belonged to various biotypes. DNA sequencing of the fucose operon region confirmed complete deletion of genes in the operon in seven of the nine isolates, while in the remaining two isolates, some of the genes were found intact or in parts. The significance of these findings is discussed.     

Multilocus Sequence Typing and Virulence Profiles in Uropathogenic Escherichia coli Isolated from Cats in the United States

The population structure, virulence, and antimicrobial resistance of uropathogenic E. coli (UPEC) from cats are rarely characterized. The aim of this study was to compare and characterize the UPEC isolated from cats in four geographic regions of USA in terms of their multilocus sequence typing (MLST), virulence profiles, clinical signs, antimicrobial resistance and phylogenetic grouping. The results showed that a total of 74 E. coli isolates were typed to 40 sequence types with 10 being novel. The most frequent phylogenetic group was B2 (n = 57). The most frequent sequence types were ST73 (n = 12) and ST83 (n = 6), ST73 was represented by four multidrug resistant (MDR) and eight non-multidrug resistant (SDR) isolates, and ST83 were significantly more likely to exhibit no drug resistant (NDR) isolates carrying the highest number of virulence genes. Additionally, MDR isolates were more diverse, and followed by SDR and NDR isolates in regards to the distribution of the STs. afa/draBC was the most prevalent among the 29 virulence-associated genes. Linking virulence profile and antimicrobial resistance, the majority of virulence-associated genes tested were more prevalent in NDR isolates, and followed by SDR and MDR isolates. Twenty (50%) MLST types in this study have previously been associated with human isolates, suggesting that these STs are potentially zoonotic. Our data enhanced the understanding of E. coli population structure and virulence association from cats. The diverse and various combinations of virulence-associated genes implied that the infection control may be challenging.     

Comparison of genotypes and antibiotic resistance of Campylobacter jejuni isolated from humans and slaughtered chickens in Switzerland

Aims: To get an overview of genotypes and antibiotic resistances in Swiss Campylobacter jejuni implicated in human gastroenteritis and to examine the association with isolates from chickens. Methods and Results: Multilocus sequence typing (MLST) and flaB typing were applied to 136 human clinical isolates. Phenotypic resistance to 12 antimicrobials and genotypic resistance to macrolides and quinolones were determined. MLST resulted in 35 known and six new sequence types (ST). The flaB analysis revealed 35 different types, which – in combination with MLST – increased the resolution of the typing approach. Resistance to quinolones, tetracycline and ampicillin was found in 37·5, 33·1 and 8·1% of the isolates, respectively, whereas macrolide resistance was found only once. Genotypic and phenotypic resistance correlated in all cases. A comparison to Camp. jejuni isolated from slaughtered chickens was performed. While 86% of the quinolone‐sensitive human isolates showed overlapping MLST‐flaB types with those of chicken origin, resistant strains showed only 39% of matching types. Conclusion: Mainly quinolone‐sensitive Camp. jejuni strains implicated in human campylobacteriosis showed matching genotypes with isolates originating from chickens. Significance and Impact of the Study: A large proportion of human cases in Switzerland are likely to originate from domestic chickens, confirming that prevention measures in the poultry production are important.     

Multicolor Melting Curve Analysis-Based Multilocus Melt Typing of Vibrio parahaemolyticus

Vibrio parahaemolyticus is the leading cause of seafood-borne gastroenteritis outbreaks. To track the source of these diseases in a timely manner, a high throughput typing method is critical. We hereby describe a novel genotyping method for V. parahaemolyticus, termed multilocus melt typing (MLMT), based on multilocus sequence typing (MLST). MLMT utilizes melting curve analysis to interrogate the allelic types of a set of informative single nucleotide polymorphisms (SNPs) derived from the housekeeping genes used in MLST. For each SNP, one allelic type generates distinct T_m values, which are converted into a binary code. Multiple SNPs thus generate a series of binary codes, forming a melt type (MT) corresponding with a sequence type (ST) of MLST. Using a set of 12 SNPs, the MLMT scheme could resolve 218 V.parahaemolyticus isolates into 50 MTs corresponding with 56 STs. The discriminatory power of MLMT and MLST was similar with Simpson’s index of diversity of 0.638 and 0.646, respectively. The global (adjusted Rand index = 0.982) and directional congruence (adjusted Wallace coefficient, MT→ST = 0.965; ST→MT = 1.000) between the two typing approaches was high. The entire procedure of MLMT could be finished within 3 h with negligible hands on time in a real-time PCR machine. We conclude that MLMT provides a reliable and efficient approach for V. parahaemolyticus genotyping and might also find use in other pathogens.     

Detecting genetic introgression: high levels of intersubspecific recombination found in Xylella fastidiosa in Brazil

Documenting the role of novel mutation versus homologous recombination in bacterial evolution, and especially in the invasion of new hosts, is central to understanding the long-term dynamics of pathogenic bacteria. We used multilocus sequence typing (MLST) to study this issue in Xylella fastidiosa subsp. pauca from Brazil, a bacterium causing citrus variegated chlorosis (CVC) and coffee leaf scorch (CLS). All 55 citrus isolates typed (plus one coffee isolate) defined three similar sequence types (STs) dominated by ST11 (85%), while the remaining 22 coffee isolates defined two STs, mainly ST16 (74%). This low level of variation masked unusually large allelic differences (>1% divergence with no intermediates) at five loci (leuA, petC, malF, cysG, and holC). We developed an introgression test to detect whether these large differences were due to introgression via homologous recombination from another X. fastidiosa subspecies. Using additional sequencing around these loci, we established that the seven randomly chosen MLST targets contained seven regions of introgression totaling 2,172 bp of 4,161 bp (52%), only 409 bp (10%) of which were detected by other recombination tests. This high level of introgression suggests the hypothesis that X. fastidiosa subsp. pauca became pathogenic on citrus and coffee (crops cultivated in Brazil for several hundred years) only recently after it gained genetic variation via intersubspecific recombination, facilitating a switch from native hosts. A candidate donor is the subspecies infecting plum in the region since 1935 (possibly X. fastidiosa subsp. multiplex). This hypothesis predicts that nonrecombinant native X. fastidiosa subsp. pauca (not yet isolated) does not cause disease in citrus or coffee.