Evolutionary Genetic Analysis Uncovers Multiple Species with Distinct Habitat Preferences and Antibiotic Resistance Phenotypes in the Complex

The genus () has a broad environmental distribution. is its best known species because it is a globally emerging, multidrug-resistant (MDR), opportunistic pathogen. Members of this species are known to display high genetic, ecological and phenotypic diversity, forming the so-called complex (Smc). Heterogeneous resistance and virulence phenotypes have been reported for environmental Smc isolates of diverse ecological origin. We hypothesized that this heterogeneity could be in part due to the potential lumping of several cryptic species in the Smc. Here we used state-of-the-art phylogenetic and population genetics methods to test this hypothesis based on the multilocus dataset available for the genus at pubmlst.org. It was extended with sequences from complete and draft genome sequences to assemble a comprehensive set of reference sequences. This framework was used to analyze 108 environmental isolates obtained in this study from the sediment and water column of four rivers and streams in Central Mexico, affected by contrasting levels of anthropogenic pollution. The aim of the study was to identify species in this collection, defined as genetically cohesive sequence clusters, and to determine the extent of their genetic, ecological and phenotypic differentiation. The multispecies coalescent, coupled with Bayes factor analysis was used to delimit species borders, together with population genetic structure analyses, recombination and gene flow estimates between sequence clusters. These analyses consistently revealed that the Smc contains at least 5 significantly differentiated lineages: and Smc1 to Smc4. Only was found to be intrinsically MDR, all its members expressing metallo-β-lactamases (MBLs). The other Smc lineages were not MDR and did not express MBLs. We also obtained isolates related to , and . They were significantly more susceptible to antibiotics than . We demonstrate that the sympatric lineages recovered display significantly differentiated habitat preferences, antibiotic resistance profiles and β-lactamase expression phenotypes, as shown by diverse multivariate analyses and robust univariate statistical tests. We discuss our data in light of current models of bacterial speciation, which fit these data well, stressing the implications of species delimitation in ecological, evolutionary and clinical research.