The massive accumulation of genome-sequences in public databases promoted the proliferation of genome-level phylogenetic analyses in many areas of biological research. However, due to diverse evolutionary and genetic processes, many loci have undesirable properties for phylogenetic reconstruction. These, if undetected, can result in erroneous or biased estimates, particularly when estimating species trees from concatenated datasets. To deal with these problems, we developed GET_PHYLOMARKERS, a pipeline designed to identify high-quality markers to estimate robust genome phylogenies from the orthologous clusters, or the pan-genome matrix (PGM), computed by GET_HOMOLOGUES. In the first context, a set of sequential filters are applied to exclude recombinant alignments and those producing anomalous or poorly resolved trees. Multiple sequence alignments and maximum likelihood (ML) phylogenies are computed in parallel on multi-core computers. A ML species tree is estimated from the concatenated set of top-ranking alignments at the DNA or protein levels, using either FastTree or IQ-TREE (IQT). The latter is used by default due to its superior performance revealed in an extensive benchmark analysis. In addition, parsimony and ML phylogenies can be estimated from the PGM. We demonstrate the practical utility of the software by analyzing 170 Stenotrophomonas genome sequences available in RefSeq and 10 new complete genomes of Mexican environmental S. maltophilia complex (Smc) isolates reported herein. A combination of core-genome and PGM analyses was used to revise the molecular systematics of the genus. An unsupervised learning approach that uses a goodness of clustering statistic identified 20 groups within the Smc at a core-genome average nucleotide identity (cgANIb) of 95.9% that are perfectly consistent with strongly supported clades on the core- and pan-genome trees. In addition, we identified 16 misclassified RefSeq genome sequences, 14 of them labeled as S. maltophilia, demonstrating the broad utility of the software for phylogenomics and geno-taxonomic studies. The code, a detailed manual and tutorials are freely available for Linux/UNIX servers under the GNU GPLv3 license at https://github.com/vinuesa/get_phylomarkers. A docker image bundling GET_PHYLOMARKERS with GET_HOMOLOGUES is available at https://hub.docker.com/r/csicunam/get_homologues/, which can be easily run on any platform.

The genus Stenotrophomonas (Gammaproteobacteria) has a broad environmental distribution. Stenotrophomonas maltophilia 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 S. maltophilia 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: S. maltophilia and Smc1 to Smc4. Only S. maltophilia 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 S. acidaminiphila, S. humi and S. terrae. They were significantly more susceptible to antibiotics than S. maltophilia. 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.

Antibiotic resistance is a major human health problem. While health care facilities are main contributors to the emergence, evolution and spread of antibiotic resistance, other ecosystems are involved in such dissemination. Wastewater, farm animals and pets have been considered important contributors to the development of antibiotic resistance. Herein, we review the impact of wildlife in such problem. Current evidence supports that the presence of antibiotic resistance genes and/or antibiotic resistant bacteria in wild animals is a sign of anthropic pollution more than of selection of resistance. However, once antibiotic resistance is present in the wild, wildlife can contribute to its transmission across different ecosystems. Further, the finding that antibiotic resistance genes, currently causing problems at hospitals, might spread through horizontal gene transfer among the bacteria present in the microbiomes of ubiquitous animals as cockroaches, fleas or rats, supports the possibility that these organisms might be bioreactors for the horizontal transfer of antibiotic resistance genes among human pathogens. The contribution of wildlife in the spread of antibiotic resistance among different hosts and ecosystems occurs at two levels. Firstly, in the case of non-migrating animals, the transfer will take place locally; a One Health problem. Paradigmatic examples are the above mentioned animals that cohabit with humans and can be reservoirs and vehicles for antibiotic resistance dissemination. Secondly, migrating animals, such as gulls, fishes or turtles may participate in the dissemination of antibiotic resistance across different geographic areas, even between different continents, which constitutes a Global Health issue.

The use and misuse of antibiotics have made antibiotic-resistant bacteria widespread nowadays, constituting one of the most relevant challenges for human health at present. Among these bacteria, opportunistic pathogens with an environmental, non-clinical, primary habitat stand as an increasing matter of concern at hospitals. These organisms usually present low susceptibility to antibiotics currently used for therapy. They are also proficient in acquiring increased resistance levels, a situation that limits the therapeutic options for treating the infections they cause. In this article, we analyse the most predominant opportunistic pathogens with an environmental origin, focusing on the mechanisms of antibiotic resistance they present. Further, we discuss the functions, beyond antibiotic resistance, that these determinants may have in the natural ecosystems that these bacteria usually colonize. Given the capacity of these organisms for colonizing different habitats, from clinical settings to natural environments, and for infecting different hosts, from plants to humans, deciphering their population structure, their mechanisms of resistance and the role that these mechanisms may play in natural ecosystems is of relevance for understanding the dissemination of antibiotic resistance under a One-Health point of view.