Rational Design of DNA Sequence-Based Strategies for Subtyping <i>Listeria monocytogenes</i>

<jats:title>ABSTRACT</jats:title> <jats:p> The ability to differentiate bacteria beyond the species level is essential for identifying and tracking infectious disease outbreaks and to improve our knowledge of the population genetics, epidemiology, and ecology of bacterial pathogens. Commonly used subtyping methods, such as serotyping, phage typing, ribotyping, and pulsed-field gel electrophoresis, can yield ambiguous results that are difficult to standardize and share among laboratories. DNA sequence-based subtyping strategies can reduce interpretation ambiguity. We report the development of a rational approach for designing sequence-based subtyping methods. <jats:italic>Listeria monocytogenes</jats:italic> was selected as the model organism for testing the efficacy of this approach. Two housekeeping genes ( <jats:italic>recA</jats:italic> and <jats:italic>prs</jats:italic> ), one stress response gene ( <jats:italic>sigB</jats:italic> ), two virulence genes ( <jats:italic>actA</jats:italic> and <jats:italic>inlA</jats:italic> ), and two intergenic regions ( <jats:italic>hly</jats:italic> - <jats:italic>mpl</jats:italic> and <jats:italic>plcA</jats:italic> - <jats:italic>hly</jats:italic> ) were sequenced for 15 <jats:italic>L. monocytogenes</jats:italic> isolates. Isolates were chosen from a representative collection of more than 1,000 <jats:italic>L. monocytogenes</jats:italic> isolates to reflect the genetic diversity of this species. DNA sequences were aligned, and sliding window analyses were performed for each gene to define 600-bp-long regions that were (i) most polymorphic (using ProSeq) or (ii) most discriminatory (using a new algorithm implemented in WINDOWMIN). Complete gene sequences for <jats:italic>actA</jats:italic> (1,929 bp) and <jats:italic>inlA</jats:italic> (2,235 bp) provided the highest discrimination (identifying 15 and 14 allelic types, respectively). WINDOWMIN allowed identification of 600-bp regions within these genes that provided similar discriminatory power (yielding 15 and 13 allelic types, respectively). The most discriminatory 600-bp fragments identified in the housekeeping and stress response genes differentiated the isolates into 8 to 10 subtypes; intergenic region sequences yielded 8 and 12 allelic types based on 335- and 242-bp sequences for <jats:italic>hly</jats:italic> - <jats:italic>mpl</jats:italic> and <jats:italic>plcA</jats:italic> - <jats:italic>hly</jats:italic> , respectively. Regions identified as most polymorphic were not necessarily most discriminatory; therefore, application of the WINDOWMIN algorithm provided a powerful tool for determining the best target regions for DNA sequence-based subtyping. Our specific results also show that inclusion of virulence gene target sequences in a DNA sequence-based subtyping scheme for <jats:italic>L. monocytogenes</jats:italic> is necessary to achieve maximum subtype differentiation. </jats:p>