Mitochondrial DNA variation in Aedes aegypti (Diptera: Culicidae) mosquitoes from Jeddah, Saudi Arabia

Wolbachia (Hertig 1936) (Rickettsiales: Ehrlichiaceae) has emerged as a valuable biocontrol tool in the fight against dengue by suppressing the transmission of the virus through mosquitoes. Monitoring the dynamics of Wolbachia is crucial for evaluating the effectiveness of release programs. Mitochondrial (mtDNA) markers serve as important tools for molecular tracking of infected mitochondrial backgrounds over time but require an understanding of the variation in release sites. In this study, we investigated the mitochondrial lineages of Aedes aegypti (Linnaeus 1762) in Jeddah, Saudi Arabia, which is a prospective release site for the "wAlbBQ" Wolbachia-infected strain of this mosquito species. We employed a combination of comprehensive mitogenomic analysis (including all protein-coding genes) and mtDNA marker analysis (cox1 and nad5) using data collected from Jeddah. We combined our mitogenome and mtDNA marker data with those from previous studies to place mitochondrial variation in Saudi Arabia into a broader global context. Our findings revealed the presence of 4 subclades that can be broadly categorized into 2 major mitochondrial lineages. Ae. aegypti mosquitoes from Jeddah belonged to both major lineages. Whilst mitogenomic data offered a higher resolution for distinguishing Jeddah mosquitoes from the wAlbBQ strain, the combination of cox1 and nad5 mtDNA markers alone proved to be sufficient. This study provides the first important characterization of Ae. aegypti mitochondrial lineages in Saudi Arabia and offers essential baseline information for planning future molecular monitoring efforts during the release of Wolbachia-infected mosquitoes.

Adaptive evolution to the natural and anthropogenic environment in a global invasive crop pest, the cotton bollworm

The cotton bollworm, Helicoverpa armigera, is set to become the most economically devastating crop pest in the world, threatening food security and biosafety as its range expands across the globe. Key to understanding the eco-evolutionary dynamics of H. armigera, and thus its management, is an understanding of population connectivity and the adaptations that allow the pest to establish in unique environments. We assembled a chromosome-scale reference genome and re-sequenced 503 individuals spanning the species range to delineate global patterns of connectivity, uncovering a previously cryptic population structure. Using a genome-wide association study (GWAS) and cell line expression of major effect loci, we show that adaptive changes in a temperature- and light-sensitive developmental pathway enable facultative diapause and that adaptation of trehalose synthesis and transport underlies cold tolerance in extreme environments. Incorporating extensive pesticide resistance monitoring, we also characterize a suite of novel pesticide and Bt resistance alleles under selection in East China. These findings offer avenues for more effective management strategies and provide insight into how insects adapt to variable climatic conditions and newly colonized environments.

Close Kin Dyads Indicate Intergenerational Dispersal and Barriers

AbstractThe movement of individuals through continuous space is typically constrained by dispersal ability and dispersal barriers. A range of approaches have been developed to investigate these. Kindisperse is a new approach that infers recent intergenerational dispersal (σ) from close kin dyads and appears particularly useful for investigating taxa that are difficult to observe individually. This study, focusing on the mosquito Aedes aegypti, shows how the same close kin data can also be used for barrier detection. We empirically demonstrate this new extension of the method using genome-wide sequence data from 266 Ae. aegypti. First, we use the spatial distribution of full-sib dyads collected within one generation to infer past movements of ovipositing female mosquitoes. These dyads indicated the relative barrier strengths of two roads and performed favorably against alternative genetic methods for detecting barriers. We then use Kindisperse to quantify recent intergenerational dispersal (σ=81.5-197.1 m generation^-1/2) from the difference in variance between the sib and the first cousin spatial distributions and, from this, estimate effective population density (ρ=833-4,864 km^-2). Dispersal estimates showed general agreement with those from mark-release-recapture studies. Barriers, σ, ρ, and neighborhood size (331-526) can inform forthcoming releases of dengue-suppressing Wolbachia bacteria into this mosquito population.

Spatio-Temporal Modelling Informing Wolbachia Replacement Releases in a Low Rainfall Climate

Releases of Aedes aegypti carrying Wolbachia bacteria are known to suppress arbovirus transmission and reduce the incidence of vector-borne diseases. In planning for Wolbachia releases in the arid environment of Jeddah, Saudi Arabia, we collected entomological data with ovitraps across a 7-month period in four locations. Herein, we show that mosquito presence in basements does not differ from that of non-basement areas of buildings. In modelling mosquito presence across the study sites, we found the spatial structure to be statistically significant in one of the four sites, while a significant spatial structure was found for egg production data across three of the four sites. The length scales of the spatial covariance functions fitted to the egg production data ranged from 143 m to 574 m, indicating that high productivity regions can be extensive in size. Rank-correlation analyses indicated that mosquito presence tended to persist from the dry to wet season, but that egg production ranks at locations could reverse. The data suggest that, in Jeddah, the quality of the local environment for breeding can vary over time. The data support the feasibility of dry season releases but with release numbers needing to be flexible depending on local rates of invasion.

Unveiling biogeographic patterns in the worldwide distributed Ceratitis capitata (medfly) using populations genomics and microbiome composition

Invasive species are among the most important, growing threats to food security and agricultural systems. The Mediterranean medfly, Ceratitis capitata, is one of the most damaging representatives of a group of rapidly expanding species in the Tephritidae family, due to their wide host range and high invasiveness potential. Here, we used restriction site-associated DNA sequencing (RADseq) to investigate the population genomic structure and phylogeographic history of medflies collected from six sampling sites, including Africa (South Africa), the Mediterranean (Spain, Greece), Latin America (Guatemala, Brazil) and Australia. A total of 1,907 single nucleotide polymorphisms (SNPs) were used to identify two genetic clusters separating native and introduced ranges, consistent with previous findings. In the introduced range, all individuals were assigned to one genetic cluster except for those in Brazil, which showed introgression of an additional genetic cluster that also appeared in South Africa, and which could not be previously identified using microsatellite markers. Moreover, we assessed the microbial composition variations in medfly populations from selected sampling sites using amplicon sequencing of the 16S ribosomal RNA (V4 region). Microbiome composition and structure were highly similar across geographic regions and host plants, and only the Brazilian specimens showed increased diversity levels and a unique composition of its microbiome compared to other sampling sites. The unique SNP patterns and microbiome features in the Brazilian specimens could point to a direct migration route from Africa with subsequent adaptation of the microbiota to the specific conditions present in Brazil. These findings significantly improve our understanding of the evolutionary history of the global medfly invasions and their adaptation to newly colonised environments.

Voltage-sensitive sodium channel (Vssc) mutations associated with pyrethroid insecticide resistance in Aedes aegypti (L.) from two districts of Jeddah, Kingdom of Saudi Arabia: baseline information for a Wolbachia release program

Background

Dengue suppression often relies on control of the mosquito vector, Aedes aegypti, through applications of insecticides of which the pyrethroid group has played a dominant role. Insecticide resistance is prevalent in Ae. aegypti around the world, and the resulting reduction of insecticide efficacy is likely to exacerbate the impact of dengue. Dengue has been a public health problem in Saudi Arabia, particularly in Jeddah, since its discovery there in the 1990s, and insecticide use for vector control is widespread throughout the city. An alternative approach to insecticide use, based on blocking dengue transmission in mosquitoes by the endosymbiont Wolbachia, is being trialed in Jeddah following the success of this approach in Australia and Malaysia. Knowledge of insecticide resistance status of mosquito populations in Jeddah is a prerequisite for establishing a Wolbachia-based dengue control program as releases of Wolbachia mosquitoes succeed when resistance status of the release population is similar to that of the wild population.

Methods

WHO resistance bioassays of mosquitoes with deltamethrin, permethrin and DDT were used in conjunction with TaqMan^® SNP Genotyping Assays to characterize mutation profiles of Ae. aegypti.

Results

Screening of the voltage-sensitive sodium channel (Vssc), the pyrethroid target site, revealed mutations at codons 989, 1016 and 1534 in Ae. aegypti from two districts of Jeddah. The triple mutant homozygote (1016G/1534C/989P) was confirmed from Al Safa and Al Rawabi. Bioassays with pyrethroids (Type I and II) and DDT showed that mosquitoes were resistant to each of these compounds based on WHO definitions. An association between Vssc mutations and resistance was established for the Type II pyrethroid, deltamethrin, with one genotype (989P/1016G/1534F) conferring a survival advantage over two others (989S/1016V/1534C and the triple heterozygote). An indication of synergism of Type I pyrethroid activity with piperonyl butoxide suggests that detoxification by cytochrome P450s accounts for some of the pyrethroid resistance response in Ae. aegypti populations from Jeddah.

Conclusions

The results provide a baseline for monitoring and management of resistance as well as knowledge of Vssc genotype frequencies required in Wolbachia release populations to ensure homogeneity with the target field population. Vssc mutation haplotypes observed show some similarity with those from Ae. aegypti in southeast Asia and the Indo-Pacific, but the presence of the triple mutant haplotype in three genotypes indicates that the species in this region may have a unique population history.

Graphical Abstract

Supplementary Information

The online version contains supplementary material available at 10.1186/s13071-021-04867-3.

On species delimitation, hybridization and population structure of cassava whitefly in Africa

The Bemisia cassava whitefly complex includes species that cause severe crop damage through vectoring cassava viruses in eastern Africa. Currently, this whitefly complex is divided into species and subgroups (SG) based on very limited molecular markers that do not allow clear definition of species and population structure. Based on 14,358 genome-wide SNPs from 62 Bemisia cassava whitefly individuals belonging to sub-Saharan African species (SSA1, SSA2 and SSA4), and using a well-curated mtCOI gene database, we show clear incongruities in previous taxonomic approaches underpinned by effects from pseudogenes. We show that the SSA4 species is nested within SSA2, and that populations of the SSA1 species comprise well-defined south-eastern (Madagascar, Tanzania) and north-western (Nigeria, Democratic Republic of Congo, Burundi) putative sub-species. Signatures of allopatric incipient speciation, and the presence of a 'hybrid zone' separating the two putative sub-species were also detected. These findings provide insights into the evolution and molecular ecology of a highly cryptic hemipteran insect complex in African, and allow the systematic use of genomic data to be incorporated in the development of management strategies for this cassava pest.

Genomic insights into a population of introduced European rabbits Oryctolagus cuniculus in Australia and the development of genetic resistance to rabbit hemorrhagic disease virus

The European rabbit (Oryctolagus cuniculus) is one of the most devastating invasive species in Australia. Since the 1950s, myxoma virus (MYXV) and rabbit haemorrhagic disease virus (RHDV) have been used to manage overabundant rabbit populations. Resistance to MYXV was observed within a few years of the release. More recently, resistance to lethal RHDV infection has also been reported, undermining the efficiency of landscape-scale rabbit control. Previous studies suggest that genetic resistance to lethal RHDV infection may differ locally between populations, yet the mechanisms of genetic resistance remain poorly understood. Here, we used genotyping by sequencing (GBS) data representing a reduced representation of the genome, to investigate Australian rabbit populations. Our aims were to understand the relationship between populations and identify possible genomic signatures of selection for RHDV resistance. One population we investigated had previously been reported to show levels of resistance to lethal RHDV infection. This population was compared to three other populations with lower or no previously reported RHDV resistance. We identified a set of novel candidate genes that could be involved in host-pathogen interactions such as virus binding and infection processes. These genes did not overlap with previous studies on RHDV resistance carried out in different rabbit populations, suggesting that multiple mechanisms are feasible. These findings provide useful insights into the different potential mechanisms of genetic resistance to RHDV virus which will inform future functional studies in this area.

Adaptive Introgression across Semipermeable Species Boundaries between Local Helicoverpa zea and Invasive Helicoverpa armigera Moths

Hybridization between invasive and native species has raised global concern, given the dramatic increase in species range shifts and pest outbreaks due to anthropogenic dispersal. Nevertheless, secondary contact between sister lineages of local and invasive species provides a natural laboratory to understand the factors that determine introgression and the maintenance or loss of species barriers. Here, we characterize the early evolutionary outcomes following secondary contact between invasive Helicoverpa armigera and native H. zea in Brazil. We carried out whole-genome resequencing of Helicoverpa moths from Brazil in two temporal samples: during the outbreak of H. armigera in 2013 and 2017. There is evidence for a burst of hybridization and widespread introgression from local H. zea into invasive H. armigera coinciding with H. armigera expansion in 2013. However, in H. armigera, the admixture proportion and the length of introgressed blocks were significantly reduced between 2013 and 2017, suggesting selection against admixture. In contrast to the genome-wide pattern, there was striking evidence for adaptive introgression of a single region from the invasive H. armigera into local H. zea, including an insecticide resistance allele that increased in frequency over time. In summary, despite extensive gene flow after secondary contact, the species boundaries are largely maintained except for the single introgressed region containing the insecticide-resistant locus. We document the worst-case scenario for an invasive species, in which there are now two pest species instead of one, and the native species has acquired resistance to pyrethroid insecticides through introgression.

Population genetics and migration pathways of the Mediterranean fruit fly Ceratitis capitata inferred with coalescent methods

Background

Invasive species are a growing threat to food biosecurity and cause significant economic losses in agricultural systems. Despite their damaging effect, they are attractive models for the study of evolution and adaptation in newly colonised environments. The Mediterranean fruit fly, Ceratitis capitata, as a member of the family Tephritidae, is one of the most studied invasive species feeding on many fruit crops in the tropics and subtropics worldwide. This study aims to determine the global macrogeographic population structure of Ceratitis capitata and reconstruct its potential migration routes.

Method

A partial mitochondrial cytochrome oxidase I gene from >400 individual medflies and 14 populations from four continents was sequenced and subjected to Bayesian demographic modelling.

Results

The Afrotropical populations (Kenya, South Africa and Ghana) harbour the majority of haplotypes detected, which also are highly divergent, in accordance with the presumed ancestral range of medflies in Sub-Saharan Africa. All other populations in the presumed non-native areas were dominated by a single haplotype also present in South Africa, in addition to a few, closely related haplotypes unique to a single local population or regional set, but missing from Africa. Bayesian coalescence methods revealed recent migration pathways from Africa to all continents, in addition to limited bidirectional migration among many local and intercontinental routes.

Conclusion

The detailed investigation of the recent migration history highlights the interconnectedness of affected crop production regions worldwide and pinpoints the routes and potential source areas requiring more specific quarantine measures.

The Trouble with MEAM2: Implications of Pseudogenes on Species Delimitation in the Globally Invasive Bemisia tabaci (Hemiptera: Aleyrodidae) Cryptic Species Complex

Molecular species identification using suboptimal PCR primers can over-estimate species diversity due to coamplification of nuclear mitochondrial (NUMT) DNA/pseudogenes. For the agriculturally important whitefly Bemisia tabaci cryptic pest species complex, species identification depends primarily on characterization of the mitochondrial DNA cytochrome oxidase I (mtDNA COI) gene. The lack of robust PCR primers for the mtDNA COI gene can undermine correct species identification which in turn compromises management strategies. This problem is identified in the B. tabaci Africa/Middle East/Asia Minor clade which comprises the globally invasive Mediterranean (MED) and Middle East Asia Minor I (MEAM1) species, Middle East Asia Minor 2 (MEAM2), and the Indian Ocean (IO) species. Initially identified from the Indian Ocean island of Réunion, MEAM2 has since been reported from Japan, Peru, Turkey and Iraq. We identified MEAM2 individuals from a Peruvian population via Sanger sequencing of the mtDNA COI gene. In attempting to characterize the MEAM2 mitogenome, we instead characterized mitogenomes of MEAM1. We also report on the mitogenomes of MED, AUS, and IO thereby increasing genomic resources for members of this complex. Gene synteny (i.e., same gene composition and orientation) was observed with published B. tabaci cryptic species mitogenomes. Pseudogene fragments matching MEAM2 partial mtDNA COI gene exhibited low frequency single nucleotide polymorphisms that matched low copy number DNA fragments (<3%) of MEAM1 genomes, whereas presence of internal stop codons, loss of expected stop codons and poor primer annealing sites, all suggested MEAM2 as a pseudogene artifact and so not a real species.

Standardized molecular diagnostic tool for the identification of cryptic species within the Bemisia tabaci complex

BACKGROUND

The whitefly Bemisia tabaci complex harbours over 40 cryptic species that have been placed in 11 phylogenetically distinct clades based on the molecular characterization of partial mitochondrial DNA COI (mtCOI) gene region. Four cryptic species are currently within the invasive clade, i.e. MED, MEAM1, MEAM2 and IO. Correct identification of these species is a critical step towards implementing reliable measures for plant biosecurity and border protection; however, no standardized B. tabaci-specific primers are currently available which has caused inconsistencies in the species identification processes.

RESULTS

We report three sets of polymerase chain reaction (PCR) primers developed to amplify the mtCOI region which can be used for genotyping MED, MEAM1 and IO species, and tested these primers on 91 MED, 35 MEAM1 and five IO individuals. PCR and sequencing of amplicons identified a total of 21, six and one haplotypes in MED, MEAM1 and IO respectively, of which six haplotypes were new to the B. tabaci database.

CONCLUSION

These primer pairs enabled standardization and robust molecular species identification via mtCOI screening of the targeted invasive cryptic species and will improve quarantine decisions. Use of this diagnostic tool could be extended to other species within the complex. © 2017 Society of Chemical Industry.

Genomic innovations, transcriptional plasticity and gene loss underlying the evolution and divergence of two highly polyphagous and invasive Helicoverpa pest species

BACKGROUND

Helicoverpa armigera and Helicoverpa zea are major caterpillar pests of Old and New World agriculture, respectively. Both, particularly H. armigera, are extremely polyphagous, and H. armigera has developed resistance to many insecticides. Here we use comparative genomics, transcriptomics and resequencing to elucidate the genetic basis for their properties as pests.

RESULTS

We find that, prior to their divergence about 1.5 Mya, the H. armigera/H. zea lineage had accumulated up to more than 100 more members of specific detoxification and digestion gene families and more than 100 extra gustatory receptor genes, compared to other lepidopterans with narrower host ranges. The two genomes remain very similar in gene content and order, but H. armigera is more polymorphic overall, and H. zea has lost several detoxification genes, as well as about 50 gustatory receptor genes. It also lacks certain genes and alleles conferring insecticide resistance found in H. armigera. Non-synonymous sites in the expanded gene families above are rapidly diverging, both between paralogues and between orthologues in the two species. Whole genome transcriptomic analyses of H. armigera larvae show widely divergent responses to different host plants, including responses among many of the duplicated detoxification and digestion genes.

CONCLUSIONS

The extreme polyphagy of the two heliothines is associated with extensive amplification and neofunctionalisation of genes involved in host finding and use, coupled with versatile transcriptional responses on different hosts. H. armigera's invasion of the Americas in recent years means that hybridisation could generate populations that are both locally adapted and insecticide resistant.

Identification and preliminary characterization of chemosensory perception-associated proteins in the melon fly Bactrocera cucurbitae using RNA-seq

An investigation into proteins involved in chemosensory perception in the melon fly, Bactrocera cucurbitae (Diptera: Tephritidae) is described here using a newly generated transcriptome dataset. The melon fly is a major agricultural pest, widely distributed in the Asia-Pacific region and some parts of Africa. For this study, a transcriptome dataset was generated using RNA extracted from 4-day-old adult specimens of the melon fly. The dataset was assembled and annotated via Gene Ontology (GO) analysis. Based on this and similarity searches to data from other species, a number of protein sequences putatively involved in chemosensory reception were identified and characterized in the melon fly. This included the highly conserved "Orco" along with a number of other less conserved odorant binding protein sequences. In addition, several sequences representing putative ionotropic and gustatory receptors were also identified. This study provides a foundation for future functional studies of chemosensory proteins in the melon fly and for making more detailed comparisons to other species. In the long term, this will ultimately help in the development of improved tools for pest management.

The complete mitochondrial DNA genome of a Chloridea (Heliothis) subflexa (Lepidoptera: Noctuidae) morpho-species

We investigate the complete mitogenome of a pheromone-trapped morpho-species of Chloridea subflexa from Brazil (initially identified by the Sanger sequencing of partial mtCOI gene) as 15 323 bp (KT598688) via next generation sequencing platform. The mitogenome has an A/T rich base composition (A: 40.4%; T: 40.3%; C: 11.5%; G: 7.8%), and included 13 protein-coding genes (PCGs), 22 tRNAs, 2 ribosomal RNAs and a putative replication region (ca. 323 bp). All PCGs start with a methionine (M) amino acid except the COI gene which has an arginine (R). The trnL2 and trn-Lys genes were partially embedded within the COII gene, while the trn-His gene was completely embedded within the ND4 gene. All PCGs ends with the "TAA" stop codon except ND3 which has a "TAG" stop codon.

The Importance of Multilocus Sequence Typing: Cautionary Tales from the Bacterium Xylella fastidiosa

Microbial identification methods have evolved rapidly over the last few decades. One such method is multilocus sequence typing (MLST). MLST is a powerful tool for understanding the evolutionary dynamics of pathogens and to gain insight into their genetic diversity. We illustrate the importance of accurate typing by reporting on three problems that have arisen in the study of a single bacterial species, the plant pathogen Xylella fastidiosa. Two of these were particularly serious since they concerned contamination of important research material that has had detrimental consequences for Xylella research: the contamination of DNA used in the sequencing of an X. fastidiosa genome (Ann-1) with DNA from another X. fastidiosa strain, and the unrecognized mislabeling of a strain (Temecula1) distributed from a culture collection (ATCC). We advocate the routine use of MLST to define strains maintained in culture collections and emphasize the importance of confirming the purity of DNA submitted for sequencing. We also present a third example that illustrates the value of MLST in guiding the choice of taxonomic types. Beyond these situations, there is a strong case for MLST whenever an isolate is used experimentally, especially where genotypic differences are suspected to influence the outcome.

Detection of the Acetylcholinesterase Insecticide Resistance Mutation (G328A) in Natural Populations of Ceratitis capitata

Wild Mediterranean fruit fly specimens collected from various regions worldwide were screened for the glycine to alanine (Gly->Ala) point mutation (G328A) in the acetylcholinesterase enzyme, presumably causing resistance to organophosphates. We found that the single nucleotide polymorphism (SNP) responsible for this amino acid change is located at the beginning of exon 6 of the Ccace2 gene. The identification of the exact location of the SNP permitted PCR primer design around this site and direct sequencing of the corresponding genomic region. We detected the resistance allele in natural Mediterranean fruit fly populations from Brazil and Spain, but not from other sites in four continents. The known treatment history of sites suggests that the resistance buildup is linked to organophosphate application in the field. The PCR-based detection provides a screening method useful for monitoring Mediterranean fruit fly insecticide resistance in local populations and improving pest management strategies accordingly.

Complete mitochondrial DNA genome of Bemisia tabaci cryptic pest species complex Asia I (Hemiptera: Aleyrodidae)

The complete length of the Asia I member of the Bemisia tabaci species complex mitochondrial DNA genome (mitogenome) is 15,210 bp (GenBank accession no. KJ778614) with an A-T biased nucleotide composition (A: 32.7%; T: 42.4%; G: 14.0%; C: 10.8%). The mitogenome consists of 13 protein-coding genes (PCGs), 22 transfer RNAs (tRNAs), 2 ribosomal RNA (rRNAs) and a 467 bp putative control region which also includes the A+T rich repeat region. All PCGs have an ATA (n = 8) or ATG (n = 5) start codon. Gene synteny of Asia I is overall similar to B. afer and two other members of the B. tabaci species complex Mediterranean and New World 1, and contains the tRNA-Ser2 located between the Cytb and ND1 genes found in Mediterranean and New World 1, but which is absent in B. afer. The orientation of the tRNA-Arg in Asia I is on the "plus" strand and differed from Mediterranean which is found on the "minus" strand. The Asia I mitogenome size is currently ranked the second smallest after B. afer (14,968 bp) followed by New World 1 (15,322 bp) and Mediterranean (15,632 bp).

The importance of multilocus sequence typing: cautionary tales from the bacterium Xylella fastidiosa

Microbial identification methods have evolved rapidly over the last few decades. One such method is multilocus sequence typing (MLST). MLST is a powerful tool for understanding the evolutionary dynamics of pathogens and to gain insight into their genetic diversity. We illustrate the importance of accurate typing by reporting on three problems that have arisen in the study of a single bacterial species, the plant pathogen Xylella fastidiosa. Two of these were particularly serious since they concerned contamination of important research material that has had detrimental consequences for Xylella research: the contamination of DNA used in the sequencing of an X. fastidiosa genome (Ann-1) with DNA from another X. fastidiosa strain, and the unrecognized mislabeling of a strain (Temecula1) distributed from a culture collection (ATCC). We advocate the routine use of MLST to define strains maintained in culture collections and emphasize the importance of confirming the purity of DNA submitted for sequencing. We also present a third example that illustrates the value of MLST in guiding the choice of taxonomic types. Beyond these situations, there is a strong case for MLST whenever an isolate is used experimentally, especially where genotypic differences are suspected to influence the outcome.