Phylogenomic, structural, and cell biological analyses reveal that Stenotrophomonas maltophilia replicates in acidified Rab7A-positive vacuoles of Acanthamoeba castellanii.

Acanthamoeba species are clinically relevant free-living amoebae (FLA) ubiquitously found in soil and water bodies. Metabolically active trophozoites graze on diverse microbes via phagocytosis. However, functional studies on Rab GTPases (Rabs), which are critical for controlling vesicle trafficking and maturation, are scarce for this FLA. This knowledge gap can be partly explained by the limited genetic tools available for Acanthamoeba cell biology. Here, we developed plasmids to generate fusions of A. castellanii strain Neff proteins to the N- or C-termini of mEGFP and mCherry2. Phylogenomic and structural analyses of the 11 Neff Rab7 paralogs found in the RefSeq assembly revealed that eight of them had non-canonical sequences. After correcting the gene annotation for the Rab7A ortholog, we generated a line stably expressing an mEGFP-Rab7A fusion, demonstrating its correct localization to acidified macropinocytic and phagocytic vacuoles using fluorescence microscopy live cell imaging (LCI). Direct labeling of live Stenotrophomonas maltophilia ESTM1D_MKCAZ16_6a (Sm18) cells with pHrodo Red, a pH-sensitive dye, demonstrated that they reside within acidified, Rab7A-positive vacuoles. We constructed new mini-Tn7 delivery plasmids and tagged Sm18 with constitutively expressed mScarlet-I. Co-culture experiments of Neff trophozoites with Sm18::mTn7TC1_Pc_mScarlet-I, coupled with LCI and microplate reader assays, demonstrated that Sm18 underwent multiple replication rounds before reaching the extracellular medium via non-lytic exocytosis. We conclude that S. maltophilia belongs to the class of bacteria that can use amoeba as an intracellular replication niche within a Stenotrophomonas-containing vacuole that interacts extensively with the endocytic pathway.IMPORTANCEDiverse Acanthamoeba lineages (genotypes) are of increasing clinical concern, mainly causing amoebic keratitis and granulomatous amebic encephalitis among other infections. S. maltophilia ranks among the top 10 most prevalent multidrug-resistant opportunistic nosocomial pathogens and is a recurrent member of the microbiome hosted by Acanthamoeba and other free-living amoebae. However, little is known about the molecular strategies deployed by Stenotrophomonas for an intracellular lifestyle in amoebae and other professional phagocytes such as macrophages, which allow the bacterium to evade the immune system and the action of antibiotics. Our plasmids and easy-to-use microtiter plate co-culture assays should facilitate investigations into the cellular microbiology of Acanthamoeba interactions with Stenotrophomonas and other opportunistic pathogens, which may ultimately lead to the discovery of new molecular targets and antimicrobial therapies to combat difficult-to-treat infections caused by these ubiquitous microbes.

Extensive Cryptic Diversity and Ecological Associations Uncovered among Mexican and Global Collections of Naegleria and Vermamoeba Species by 18S Ribosomal DNA, Internal Transcribed Spacer, and Cytochrome Oxidase Subunit I Sequence Analysis.

Free-living amoebae (FLA) are phagocytic protists that play crucial roles in microbial communities as significant microbial grazers. However, our current knowledge of their diversity, ecology, and population genetic structures is marginal due to the shallow and biased sampling of ecosystems and the use of few, poorly resolving molecular markers. Thirty-two FLA were isolated from soil and water samples collected across representative ecosystems of the State of Morelos in Central Mexico, including the drinking water distribution system (DWDS) from the state capital. We classified our isolates as members of Acanthamoeba, Vermamoeba, Naegleria, and Tetramitus by 18S ribosomal DNA (rDNA) sequencing. Vermamoeba isolates were recovered exclusively from the DWDS samples. In contrast, Naegleria strains displayed a broad distribution in soil and water samples across the natural ecosystems. We used a combination of phylogenetic and population genetic analyses of internal transcribed spacer (ITS) and cytochrome oxidase subunit I (COI) sequences from our isolates and a comprehensive set of reference sequences to analyze the currently known diversity of Naegleria spp. Significant associations were uncovered between the most prevalent lineages of Naegleria and Vermamoeba and broad ecological and geographical variables at regional and global levels. The population structure and cryptic diversity within the Naegleria galeacystis-Naegleria americana and Vermamoeba vermiformis species complexes were thoroughly analyzed. Our results prove that the genus Vermamoeba, which was previously thought to consist of only one species, actually encompasses at least seven widely distributed species, as indicated by consistent evidence from Bayesian phylogenetics, two species-delimitation programs, and population genetics analyses. IMPORTANCE Our study sheds new light on the population genetic structure of V. vermiformis and diverse Naegleria species. Using improved molecular markers and advanced analytical approaches, we discovered that N. americana, previously considered a single species, actually contains multiple distinct lineages, as revealed by COI sequencing. These lineages are highly differentiated, with little gene flow between them. Our findings demonstrate that the genus Vermamoeba holds multiple cryptic species, requiring a significant taxonomic revision in light of multilocus sequence analyses. These results advance our understanding of the ecology, molecular systematics, and biogeography of these genera and species complexes at both regional and global scales. This study has significant implications for diagnosing amoebal infections and evaluating health risks associated with FLA in domestic and recreational waters.

Pangenome Analysis of Plant Transcripts and Coding Sequences.

The pangenome of a species is the sum of the genomes of its individuals. As coding sequences often represent only a small fraction of each genome, analyzing the pangene set can be a cost-effective strategy for plants with large genomes or highly heterozygous species. Here, we describe a step-by-step protocol to analyze plant pangene sets with the software GET_HOMOLOGUES-EST . After a short introduction, where the main concepts are illustrated, the remaining sections cover the installation and typical operations required to analyze and annotate pantranscriptomes and gene sets of plants. The recipes include instructions on how to call core and accessory genes, how to compute a presence-absence pangenome matrix, and how to identify and analyze private genes, present only in some genotypes. Downstream phylogenetic analyses are also discussed.

An Epistatic Network Describes oppA and glgB as Relevant Genes for Mycobacterium tuberculosis.

Mycobacterium tuberculosis is an acid-fast bacterium that causes tuberculosis worldwide. The role of epistatic interactions among different loci of the M. tuberculosis genome under selective pressure may be crucial for understanding the disease and the molecular basis of antibiotic resistance acquisition. Here, we analyzed polymorphic loci interactions by applying a model-free method for epistasis detection, SpydrPick, on a pan-genome-wide alignment created from a set of 254 complete reference genomes. By means of the analysis of an epistatic network created with the detected epistatic interactions, we found that glgB (α-1,4-glucan branching enzyme) and oppA (oligopeptide-binding protein) are putative targets of co-selection in M. tuberculosis as they were associated in the network with M. tuberculosis genes related to virulence, pathogenesis, transport system modulators of the immune response, and antibiotic resistance. In addition, our work unveiled potential pharmacological applications for genotypic antibiotic resistance inherent to the mutations of glgB and oppA as they epistatically interact with fprA and embC, two genes recently included as antibiotic-resistant genes in the catalog of the World Health Organization. Our findings showed that this approach allows the identification of relevant epistatic interactions that may lead to a better understanding of M. tuberculosis by deciphering the complex interactions of molecules involved in its metabolism, virulence, and pathogenesis and that may be applied to different bacterial populations.

A Novel Group of Promiscuous Podophages Infecting Diverse Gammaproteobacteria from River Communities Exhibits Dynamic Intergenus Host Adaptation.

Phages are generally described as species specific or even strain specific, implying an inherent limitation for some to be maintained and spread in diverse bacterial communities. Moreover, phage isolation and host range determination rarely consider the phage ecological context, likely biasing our notion on phage specificity. Here we isolated and characterized a novel group of six promiscuous phages, named Atoyac, existing in rivers and sewage by using a diverse collection of over 600 bacteria retrieved from the same environments as potential hosts. These podophages isolated from different regions in Mexico display a remarkably broad host range, infecting bacteria from six genera: Aeromonas, Pseudomonas, Yersinia, Hafnia, Escherichia, and Serratia Atoyac phage genomes are ∼42 kb long and highly similar to each other, but not to those currently available in genome and metagenome public databases. Detailed comparison of the phages' efficiency of plating (EOP) revealed variation among bacterial genera, implying a cost associated with infection of distant hosts, and between phages, despite their sequence similarity. We show, through experimental evolution in single or alternate hosts of different genera, that efficiency of plaque production is highly dynamic and tends toward optimization in hosts rendering low plaque formation. However, adaptation to distinct hosts differed between similar phages; whereas one phage optimized its EOP in all tested hosts, the other reduced plaque production in one host, suggesting that propagation in multiple bacteria may be key to maintain promiscuity in some viruses. Our study expands our knowledge of the virosphere and uncovers bacterium-phage interactions overlooked in natural systems.IMPORTANCE In natural environments, phages coexist and interact with a broad variety of bacteria, posing a conundrum for narrow-host-range phage maintenance in diverse communities. This context is rarely considered in the study of host-phage interactions, typically focused on narrow-host-range viruses and their infectivity in target bacteria isolated from sources distinct to where the phages were retrieved from. By studying phage-host interactions in bacteria and viruses isolated from river microbial communities, we show that novel phages with promiscuous host range encompassing multiple bacterial genera can be found in the environment. Assessment of hundreds of interactions in diverse hosts revealed that similar phages exhibit different infection efficiency and adaptation patterns. Understanding host range is fundamental in our knowledge of bacterium-phage interactions and their impact on microbial communities. The dynamic nature of phage promiscuity revealed in our study has implications in different aspects of phage research such as horizontal gene transfer or phage therapy.

Fibrillarin evolution through the Tree of Life: Comparative genomics and microsynteny network analyses provide new insights into the evolutionary history of Fibrillarin.

Fibrillarin (FIB), a methyltransferase essential for life in the vast majority of eukaryotes, is involved in methylation of rRNA required for proper ribosome assembly, as well as methylation of histone H2A of promoter regions of rRNA genes. RNA viral progression that affects both plants and animals requires FIB proteins. Despite the importance and high conservation of fibrillarins, there little is known about the evolutionary dynamics of this small gene family. We applied a phylogenomic microsynteny-network approach to elucidate the evolutionary history of FIB proteins across the Tree of Life. We identified 1063 non-redundant FIB sequences across 1049 completely sequenced genomes from Viruses, Bacteria, Archaea, and Eukarya. FIB is a highly conserved single-copy gene through Archaea and Eukarya lineages, except for plants, which have a gene family expansion due to paleopolyploidy and tandem duplications. We found a high conservation of the FIB genomic context during plant evolution. Surprisingly, FIB in mammals duplicated after the Eutheria split (e.g., ruminants, felines, primates) from therian mammals (e.g., marsupials) to form two main groups of sequences, the FIB and FIB-like groups. The FIB-like group transposed to another genomic context and remained syntenic in all the eutherian mammals. This transposition correlates with differences in the expression patterns of FIB-like proteins and with elevated Ks values potentially due to reduced evolutionary constraints of the duplicated copy. Our results point to a unique evolutionary event in mammals, between FIB and FIB-like genes, that led to non-redundant roles of the vital processes in which this protein is involved.

Correction to: Genome analysis of Salmonella enterica subsp. diarizonae isolates from invasive human infections reveals enrichment of virulence-related functions in lineage ST1256.

An amendment to this paper has been published and can be accessed via the original article.

The Concerted Action of Two B3-Like Prophage Genes Excludes Superinfecting Bacteriophages by Blocking DNA Entry into Pseudomonas aeruginosa.

In this study, we describe seven vegetative phage genomes homologous to the historic phage B3 that infect Pseudomonas aeruginosa Like other phage groups, the B3-like group contains conserved (core) and variable (accessory) open reading frames (ORFs) grouped at fixed regions in their genomes; however, in either case, many ORFs remain without assigned functions. We constructed lysogens of the seven B3-like phages in strain Ps33 of P. aeruginosa, a novel clinical isolate, and assayed the exclusion phenotype against a variety of temperate and virulent superinfecting phages. In addition to the classic exclusion conferred by the phage immunity repressor, the phenotype observed in B3-like lysogens suggested the presence of other exclusion genes. We set out to identify the genes responsible for this exclusion phenotype. Phage Ps56 was chosen as the study subject since it excluded numerous temperate and virulent phages. Restriction of the Ps56 genome, cloning of several fragments, and resection of the fragments that retained the exclusion phenotype allowed us to identify two core ORFs, so far without any assigned function, as responsible for a type of exclusion. Neither gene expressed separately from plasmids showed activity, but the concurrent expression of both ORFs is needed for exclusion. Our data suggest that phage adsorption occurs but that phage genome translocation to the host's cytoplasm is defective. To our knowledge, this is the first report on this type of exclusion mediated by a prophage in P. aeruginosaIMPORTANCEPseudomonas aeruginosa is a Gram-negative bacterium frequently isolated from infected immunocompromised patients, and the strains are resistant to a broad spectrum of antibiotics. Recently, the use of phages has been proposed as an alternative therapy against multidrug-resistant bacteria. However, this approach may present various hurdles. This work addresses the problem that pathogenic bacteria may be lysogenized by phages carrying genes encoding resistance against secondary infections, such as those used in phage therapy. Discovering phage genes that exclude superinfecting phages not only assigns novel functions to orphan genes in databases but also provides insight into selection of the proper phages for use in phage therapy.

Production of distinct labdane-type diterpenoids using a novel cryptic labdane-like cluster from Streptomyces thermocarboxydus K155.

Bioinformatic mining of the Streptomyces thermocarboxydus K155 genome predicted the presence of four synthases for the production of geosmin, hopene, albaflavenone, and a type B-type A diterpenoid system like that described for labdane-related diterpenoids (LRD). The lrd cluster was comprised by an operon of four genes (lrdABDC). This cluster seemed to be silent in the wild-type strain, as neither labdane nor terpene-like compounds were detected by UPLC-TOF-MS and GC-MS analyses in both culture supernatants and mycelial extracts. Heterologous expression of the lrdABDC cluster in a defective cyslabdan producer (Streptomyces cyslabdanicus K04-0144Δcld) generated 8,17-epoxy-7-hydroxy labda-12,14-diene and cyslabdan. The same cluster expressed in the strains Streptomyces coelicolor M1152, Streptomyces peucetius var. caesius, and Streptomyces avermitilis SUKA22 produced the general intermediary labda-8(17), 12(E),14-triene [(E)-biformene]. Besides (E)-biformene, S. coelicolor M1152 and S. avermitilis SUKA22 produced two and three different labdane-type diterpenoids, underlying the relevance of the genetic background of the Streptomyces host in product formation.

Population analysis of D6-like plasmid prophage variants associated with specific IncC plasmid types in the emerging Salmonella Typhimurium ST213 genotype.

The Salmonella enterica serovar Typhimurium sequence type 213 (ST213) emerged as a predominant genotype in Mexico. It is characterized by harboring multidrug resistance (MDR) IncC plasmids (previously IncA/C) and the lack of the Salmonella virulence plasmid (pSTV). Here we show that the D6-like plasmid prophage is present in most of the ST213 strains. We used the reported nucleotide sequence of YU39 plasmid (pYU39_89) to design a PCR typing scheme for the D6-like plasmid prophages, and determined the complete nucleotide sequences for the D6-like prophages of three additional ST213 strains (YU07-18, SL26 and SO21). Two prophage variants were described: i) a complete prophage, containing homologous sequences for most of the genetic modules described in P1 and D6 phages, which most likely allow for the lytic and lysogenic lifestyles; and ii) an incomplete prophage, lacking a 15 kb region containing morphogenesis genes, suggesting that it is defective. The tail fiber gene inversion region was the most divergent one between D6 and pYU39_89 genomes, suggesting the production of a distinct set of tail fibers, which could be involved in host range preferences. A glutaminyl-tRNA synthetase gene (glnS), which could be involved in providing host cell increased fitness or plasmid maintenance functions, was found in all D6-like genomes. Population level analysis revealed a biogeographic pattern of distribution of these plasmid-phages and specific associations with variants of MDR IncC plasmids. Statistically significant associations were found between the two prophage variants (p75 or p89), the type of IncC plasmids (I or II) and geographic isolation regions (Sonora, San Luis Potosí, Michoacán and Yucatán). This work integrates results from molecular typing, genomics and epidemiology to provide a broad overview for the evolution of an emergent Salmonella genotype.

Global phylogeography and ancient evolution of the widespread human gut virus crAssphage.

Microbiomes are vast communities of microorganisms and viruses that populate all natural ecosystems. Viruses have been considered to be the most variable component of microbiomes, as supported by virome surveys and examples of high genomic mosaicism. However, recent evidence suggests that the human gut virome is remarkably stable compared with that of other environments. Here, we investigate the origin, evolution and epidemiology of crAssphage, a widespread human gut virus. Through a global collaboration, we obtained DNA sequences of crAssphage from more than one-third of the world's countries and showed that the phylogeography of crAssphage is locally clustered within countries, cities and individuals. We also found fully colinear crAssphage-like genomes in both Old-World and New-World primates, suggesting that the association of crAssphage with primates may be millions of years old. Finally, by exploiting a large cohort of more than 1,000 individuals, we tested whether crAssphage is associated with bacterial taxonomic groups of the gut microbiome, diverse human health parameters and a wide range of dietary factors. We identified strong correlations with different clades of bacteria that are related to Bacteroidetes and weak associations with several diet categories, but no significant association with health or disease. We conclude that crAssphage is a benign cosmopolitan virus that may have coevolved with the human lineage and is an integral part of the normal human gut virome.

Phylogenomic Rhizobium Species Are Structured by a Continuum of Diversity and Genomic Clusters.

The bacterial genus Rhizobium comprises diverse symbiotic nitrogen-fixing species associated with the roots of plants in the Leguminosae family. Multiple genomic clusters defined by whole genome comparisons occur within Rhizobium, but their equivalence to species is controversial. In this study we investigated such genomic clusters to ascertain their significance in a species phylogeny context. Phylogenomic inferences based on complete sets of ribosomal proteins and stringent core genome markers revealed the main lineages of Rhizobium. The clades corresponding to R. etli and R. leguminosarum species show several genomic clusters with average genomic nucleotide identities (ANI > 95%), and a continuum of divergent strains, respectively. They were found to be inversely correlated with the genetic distance estimated from concatenated ribosomal proteins. We uncovered evidence of a Rhizobium pangenome that was greatly expanded, both in its chromosomes and plasmids. Despite the variability of extra-chromosomal elements, our genomic comparisons revealed only a few chromid and plasmid families. The presence/absence profile of genes in the complete Rhizobium genomes agreed with the phylogenomic pattern of species divergence. Symbiotic genes were distributed according to the principal phylogenomic Rhizobium clades but did not resolve genome clusters within the clades. We distinguished some types of symbiotic plasmids within Rhizobium that displayed different rates of synonymous nucleotide substitutions in comparison to chromosomal genes. Symbiotic plasmids may have been repeatedly transferred horizontally between strains and species, in the process displacing and substituting pre-existing symbiotic plasmids. In summary, the results indicate that Rhizobium genomic clusters, as defined by whole genomic identities, might be part of a continuous process of evolutionary divergence that includes the core and the extrachromosomal elements leading to species formation.

Minimal standards for the description of new genera and species of rhizobia and agrobacteria.

Herein the members of the Subcommittee on Taxonomy of Rhizobia and Agrobacteria of the International Committee on Systematics of Prokaryotes review recent developments in rhizobial and agrobacterial taxonomy and propose updated minimal standards for the description of new species (and genera) in these groups. The essential requirements (minimal standards) for description of a new species are (1) a genome sequence of at least the proposed type strain and (2) evidence for differentiation from other species based on genome sequence comparisons. It is also recommended that (3) genetic variation within the species is documented with sequence data from several clearly different strains and (4) phenotypic features are described, and their variation documented with data from a relevant set of representative strains. Furthermore, it is encouraged that information is provided on (5) nodulation or pathogenicity phenotypes, as appropriate, with relevant gene sequences. These guidelines supplement the current rules of general bacterial taxonomy, which require (6) a name that conforms to the International Code of Nomenclature of Prokaryotes, (7) validation of the name by publication either directly in the International Journal of Systematic and Evolutionary Microbiology or in a validation list when published elsewhere, and (8) deposition of the type strain in two international culture collections in separate countries.

Genome analysis of Salmonella enterica subsp. diarizonae isolates from invasive human infections reveals enrichment of virulence-related functions in lineage ST1256.

BACKGROUND: Salmonella enterica subsp. diarizonae (IIIb) is frequently isolated from the environment, cold-blooded reptiles, sheep and humans; however only a few studies describe the isolation of this subspecies from invasive human infections. The factors contributing to this unusual behavior are currently unknown. RESULTS: We report here the genome features of two diarizonae strains, SBO13 and SBO27, isolated from endocervical tissue collected post-abortion and from cerebrospinal fluid of a newborn child, respectively, in the city of Santa Cruz, Bolivia. Although isolated six years apart, SBO27 in 2008 and SBO13 in 2014, both strains belong to the same sequence type 1256 (ST1256) and show a high degree of genome conservation sharing more than 99% of their genes, including the conservation of a ~ 10 kb plasmid. A prominent feature of the two genomes is the presence of 24 genomic islands (GIs), in addition to 10 complete Salmonella pathogenicity islands (SPI) and fragments of SPI-7, a pathogenicity island first reported in the human-adapted serovar Typhi. Some of the GIs identified in SBO13 and SBO27 harbor genes putatively encoding auto-transporters involved in adhesion, lipopolysaccharide modifying enzymes, putative toxins, pili-related proteins, efflux pumps, and several putative membrane cation transport related-genes, among others. These two Bolivian isolates also share genes encoding the type-III secretion system effector proteins SseK2, SseK3 and SlrP with other diarizonae sequence types (ST) mainly-associated with infections in humans. The sseK2, sseK3 and slrP genes were either absent or showing frameshift mutations in a significant proportion of genomes from environmental diarizonae isolates. CONCLUSIONS: The comparative genomic study of two diarizonae strains isolated in Bolivia from human patients uncovered the presence of many genes putatively related to virulence. The statistically-significant acquisition of a unique combination of these functions by diarizonae strains isolated from humans may have impacted the ability of these isolates to successfully infect the human host.

Analyzing the functional divergence of Slo1 and Slo3 channel subfamilies.

Slo1 and Slo3 encode close paralogues of the Slo potassium (K+) channels family. Despite their evolutionary relatedness, Slo1 and Slo3 channels show marked functional differences and evolutionary dynamics. Whereas Slo1 is a highly conserved and widely expressed channel, Slo3 is a rapidly evolving channel restricted to sperm. However, the molecular mechanisms behind the structural-functional differences of Slo1 and Slo3 channels are unknown. In this study, we explored the functional divergence of Slo1 and Slo3 subfamilies in vertebrates and examined the structure-function relationships of our predictions using experimental data. We found that ∼25% of sites between Slo1 and Slo3 underwent altered functional constraints, affecting some residues with important roles in Slo1 channel gating. Because functional divergence was principally generated by accelerated evolution of Slo3 after gene duplication, we explored selective forces behind Slo3 diversification. We observed that Slo3 subjected was principally subjected to relaxation of purifying selection, but we also identified several sites evolving under positive selection in the cytosolic domain of this channel . Concerning Slo1, this channel presented strong purifying selection. Whether residues evolving under different selection in Slo1 and Slo3 are responsible for functional differences observed between these channels, as well as among Slo3 orthologs, remains to be established.

The Genomic Basis of Intrinsic and Acquired Antibiotic Resistance in the Genus Serratia.

Serratia marcescens, a member of the Enterobacteriaceae family, was long thought to be a non-pathogenic bacterium prevalent in environmental habitats. Together with other members of this genus, it has emerged in recent years as an opportunistic nosocomial pathogen causing various types of infections. One important feature of pathogens belonging to this genus is their intrinsic and acquired resistance to a variety of antibiotic families, including ß-lactam, aminoglycosides, quinolones and polypeptide antibiotics. The aim of this study was to elucidate which genes participate in the intrinsic and acquired antibiotic resistance of this genus in order to determine the Serratia genus resistome. We performed phylogenomic and comparative genomic analyses using 32 Serratia spp. genomes deposited in the NCBI GenBank from strains isolated from different ecological niches and different lifestyles. S. marcescens strain SmUNAM836, which was previously isolated from a Mexican adult with obstructive pulmonary disease, was included in this study. The results show that most of the antibiotic resistance genes (ARGs) were found on the chromosome, and to a lesser degree, on plasmids and transposons acquired through horizontal gene transfer. Four strains contained the gyrA point mutation in codon Ser83 that confers quinolone resistance. Pathogenic and environmental isolates presented a high number of ARGs, especially genes associated with efflux systems. Pathogenic strains, specifically nosocomial strains, presented more acquired resistance genes than environmental isolates. We may conclude that the environment provides a natural reservoir for antibiotic resistance, which has been underestimated in the medical field.

GET_PHYLOMARKERS, a Software Package to Select Optimal Orthologous Clusters for Phylogenomics and Inferring Pan-Genome Phylogenies, Used for a Critical Geno-Taxonomic Revision of the Genus Stenotrophomonas.

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.

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

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.

Complete Genome Sequences of Eight Rhizobium Symbionts Associated with Common Bean (Phaseolus vulgaris).

We present here the high-quality complete genome sequences of eight strains of Rhizobium-nodulating Phaseolus vulgaris Comparative analyses showed that some of them belonged to different genomic and evolutionary lineages with common symbiotic properties. Two novel symbiotic plasmids (pSyms) with P. vulgaris specificity are reported here.