Replacing and additive horizontal gene transfer in Streptococcus

Investigating Additive and Replacing Horizontal Gene Transfers Using Phylogenies and Whole Genomes

Horizontal gene transfer (HGT) is fundamental to microbial evolution and adaptation. When a gene is horizontally transferred, it may either add itself as a new gene to the recipient genome (possibly displacing nonhomologous genes) or replace an existing homologous gene. Currently, studies do not usually distinguish between "additive" and "replacing" HGTs, and their relative frequencies, integration mechanisms, and specific roles in microbial evolution are poorly understood. In this work, we develop a novel computational framework for large-scale classification of HGTs as either additive or replacing. Our framework leverages recently developed phylogenetic approaches for HGT detection and classifies HGTs inferred between terminal edges based on gene orderings along genomes and phylogenetic relationships between the microbial species under consideration. The resulting method, called DART, is highly customizable and scalable and can classify a large fraction of inferred HGTs with high confidence and statistical support. Our application of DART to a large dataset of thousands of gene families from 103 Aeromonas genomes provides insights into the relative frequencies, functional biases, and integration mechanisms of additive and replacing HGTs. Among other results, we find that (i) the relative frequency of additive HGT increases with increasing phylogenetic distance, (ii) replacing HGT dominates at shorter phylogenetic distances, (iii) additive and replacing HGTs have strikingly different functional profiles, (iv) homologous recombination in flanking regions of a novel gene may be a frequent integration mechanism for additive HGT, and (v) phages and mobile genetic elements likely play an important role in facilitating additive HGT.

Flexible genomic island conservation across freshwater and marine Methylophilaceae

The evolutionary trajectory of Methylophilaceae includes habitat transitions from freshwater sediments to freshwater and marine pelagial that resulted in genome reduction (genome-streamlining) of the pelagic taxa. However, the extent of genetic similarities in the genomic structure and microdiversity of the two genome-streamlined pelagic lineages (freshwater "Ca. Methylopumilus" and the marine OM43 lineage) has so far never been compared. Here, we analyzed complete genomes of 91 "Ca. Methylopumilus" strains isolated from 14 lakes in Central Europe and 12 coastal marine OM43 strains. The two lineages showed a remarkable niche differentiation with clear species-specific differences in habitat preference and seasonal distribution. On the other hand, we observed a synteny preservation in their genomes by having similar locations and types of flexible genomic islands (fGIs). Three main fGIs were identified: a replacement fGI acting as phage defense, an additive fGI harboring metabolic and resistance-related functions, and a tycheposon containing nitrogen-, thiamine-, and heme-related functions. The fGIs differed in relative abundances in metagenomic datasets suggesting different levels of variability ranging from strain-specific to population-level adaptations. Moreover, variations in one gene seemed to be responsible for different growth at low substrate concentrations and a potential biogeographic separation within one species. Our study provides a first insight into genomic microdiversity of closely related taxa within the family Methylophilaceae and revealed remarkably similar dynamics involving mobile genetic elements and recombination between freshwater and marine family members.

New variant strain of Streptococcus canis with Lancefield group C isolated from canine otitis externa

Every basic course in microbiology teaches us, Streptococcus canis always tests positive for Lancefield group G. Surprisingly, we identified a strain of S. canis with Lancefield group C, cultured from a dog with otitis externa after lateral ear canal resection. Whole genome sequencing data and analysis points towards a horizontal gene transfer event between S. canis and S. dysgalactiae. Although these species are closely related, gene transfer in this region of the genome of S. canis has not been described before. The value of technologies as MALDI-TOF MS and sequencing in microbiological diagnostics will grow as more diverse streptococci arise that do not always conform anymore to the classical Lancefield group typing.

Leaky barriers to gene sharing between locally co-existing coagulase-negative Staphylococcus species

Coagulase-negative Staphylococcus (CoNS) are opportunistic pathogens implicated in many human and animal infections. The evolutionary history of CoNS remains obscure because of the historical lack of recognition for their clinical importance and poor taxonomic sampling. Here, we sequenced the genomes of 191 CoNS isolates representing 15 species sampled from diseased animals diagnosed in a veterinary diagnostic laboratory. We found that CoNS are important reservoirs of diverse phages, plasmids and mobilizable genes encoding antimicrobial resistance, heavy metal resistance, and virulence. Frequent exchange of DNA between certain donor-recipient partners suggests that specific lineages act as hubs of gene sharing. We also detected frequent recombination between CoNS regardless of their animal host species, indicating that ecological barriers to horizontal gene transfer can be surmounted in co-circulating lineages. Our findings reveal frequent but structured patterns of transfer that exist within and between CoNS species, which are driven by their overlapping ecology and geographical proximity.

Incidence, seasonal pattern, and clinical manifestations of Streptococcus dysgalactiae subspecies equisimilis bacteremia; a population-based study

Streptococcus dysgalactiae subspecies equisimilis (SDSE) is a human pathogen causing severe invasive infections. Population-based studies on SDSE bacteremia are limited. The purpose of this study was to investigate the incidence, seasonal pattern, clinical manifestations, and recurrence of SDSE bacteraemia. Records regarding patients aged ≥ 18 years with SDSE bacteremia in the Pirkanmaa health district in August 2015 to July 2018 were retrospectively reviewed. A total of 230 SDSE bacteremia episodes were identified, with 217 episodes (involving 211 patients) available for analysis. The mean annual incidence rate of SDSE bacteremia was 16.9/100 000 inhabitants. Most episodes (33%) were detected in the summer (June to August) (p = 0.058). Episodes with bacteremic cellulitis were statistically significantly more common during the summer compared with other seasons (p = 0.008). Cellulitis was the most common presenting clinical manifestation of SDSE bacteremia (68% of all episodes). Risk factors of recurring bacteremia were chronic eczema and/or skin erosion (OR 3.96 [95% CI 1.11-14.1]), heart disease (OR 3.56 [95% CI 1.22-10.4]), diabetes (OR 3.77 [95% CI 1.35-10.5]) and a history of cellulitis. We found a remarkably high incidence of SDSE bacteraemia in the Pirkanmaa health district. Bacteraemic cellulitis, which was the predominant clinical manifestation is more often occurred in the summer. Risk factors of recurring SDSE bacteremia were a history of cellulitis, chronic eczema or skin erosion, diabetes, and heart disease.

Major risk factors for Streptococcus dysgalactiae subsp. equisimilis bacteremia: a population-based study

BACKGROUND

Streptococcus dysgalactiae subspecies equisimilis is a human pathogen causing severe invasive infections. Detailed information on S. dysgalactiae subsp. equisimilis bacteremia and especially of predisposing factors are lacking. The purpose of the study is to investigate the risk factors of S. dysgalactiae subsp. equisimilis bacteremia compared to the general population in Finland.

METHODS

We retrospectively reviewed all patients older than 18 years with S. dysgalactiae subsp. equisimilis bacteremia in the Pirkanmaa health district from August 2015 to July 2018. The risk factors for S. dysgalactiae subsp. equisimilis bacteremia were investigated with respect to the normal population in Finland using the Finhealth study data provided by the Finnish institute for health and welfare. The study group was matched with the Finhealth study by age and sex.

RESULTS

Altogether 230 cases of S. dysgalactiae subsp. equisimilis bacteremia were detected. The medical records of 217 episodes of S. dysgalactiae subsp. equisimilis bacteremia (involving 211 patients) were available for analysis. Obesity was a statistically significant risk factor for S. dysgalactiae subsp. equisimilis bacteremia (Odds Ratio 2.96 [95% CI 2.22-3.96]). Diabetes and coronary artery disease were also associated with an increased risk of S. dysgalactiae subsp. equisimilis bacteremia (OR 4.82 [95% CI 3.62-6.42]) and (OR 3.03 [95% CI 2.18-4.19]).

CONCLUSIONS

We found obesity, diabetes, and coronary artery disease to be associated with an increased risk for S. dysgalactiae subsp. equisimilis bacteremia. These results provide an increased understanding of risk factors for S. dysgalactiae subsp. equisimilis bacteremia.

Assessment of the Relationship between Clinical Manifestation and Pathogenic Potential of Streptococcus pyogenes Strains-Distribution of Genes and Genotypes of Toxins

Streptococcus pyogenes is one of the most important species among beta-haemolytic streptococci, causing human infections of different localization. It is isolated from clinical specimens relatively frequently. In this study, the frequency and co-occurrence of toxin genes (speA, speB, speC, speH, speJ, speK) among 147 S. pyogenes strains were evaluated, using real-time PCR. In addition, the relationship between the occurrence of these genes and the origin of S. pyogenes strains from selected clinical material was assessed. The speB gene was present with the highest incidence (98.6%), while the speK gene was the least frequent (8.2%) among the tested strains. Based on the presence of the detected genes, the distribution of 17 genotypes was determined. The most common (21.8%), was speA (−) speB (+) speC (−) speH (−) speJ (−) speK (−) genotype. Furthermore, significant variation in the presence of some genes and genotypes of toxins in S. pyogenes strains isolated from different types of clinical material was found. There is a considerable variety and disproportion between the frequency of individual genes and genotypes of toxins in S. pyogenes strains. The relationship between the origin of S. pyogenes isolates and the presence of toxins genes indicates their pathogenic potential in the development of infections of selected localization.

Transposable elements promote the evolution of genome streamlining

Eukaryotes and prokaryotes have distinct genome architectures, with marked differences in genome size, the ratio of coding/non-coding DNA, and the abundance of transposable elements (TEs). As TEs replicate independently of their hosts, the proliferation of TEs is thought to have driven genome expansion in eukaryotes. However, prokaryotes also have TEs in intergenic spaces, so why do prokaryotes have small, streamlined genomes? Using an in silico model describing the genomes of single-celled asexual organisms that coevolve with TEs, we show that TEs acquired from the environment by horizontal gene transfer can promote the evolution of genome streamlining. The process depends on local interactions and is underpinned by rock-paper-scissors dynamics in which populations of cells with streamlined genomes beat TEs, which beat non-streamlined genomes, which beat streamlined genomes, in continuous and repeating cycles. Streamlining is maladaptive to individual cells, but improves lineage viability by hindering the proliferation of TEs. Streamlining does not evolve in sexually reproducing populations because recombination partially frees TEs from the deleterious effects they cause. This article is part of the theme issue 'The secret lives of microbial mobile genetic elements'.

Indirect identification of horizontal gene transfer

Several implicit methods to infer horizontal gene transfer (HGT) focus on pairs of genes that have diverged only after the divergence of the two species in which the genes reside. This situation defines the edge set of a graph, the later-divergence-time (LDT) graph, whose vertices correspond to genes colored by their species. We investigate these graphs in the setting of relaxed scenarios, i.e., evolutionary scenarios that encompass all commonly used variants of duplication-transfer-loss scenarios in the literature. We characterize LDT graphs as a subclass of properly vertex-colored cographs, and provide a polynomial-time recognition algorithm as well as an algorithm to construct a relaxed scenario that explains a given LDT. An edge in an LDT graph implies that the two corresponding genes are separated by at least one HGT event. The converse is not true, however. We show that the complete xenology relation is described by an rs-Fitch graph, i.e., a complete multipartite graph satisfying constraints on the vertex coloring. This class of vertex-colored graphs is also recognizable in polynomial time. We finally address the question “how much information about all HGT events is contained in LDT graphs” with the help of simulations of evolutionary scenarios with a wide range of duplication, loss, and HGT events. In particular, we show that a simple greedy graph editing scheme can be used to efficiently detect HGT events that are implicitly contained in LDT graphs.

Can mobile genetic elements rescue genes from extinction?

Bacteria and other prokaryotes evolve primarily through rapid changes in their gene content by quickly losing and gaining genes whenever an ecological opportunity emerges. As gene loss and horizontal gene transfer (HGT) appear to be the most common events across the prokaryotic tree of life, we need to think beyond gradual sequence evolution if we wish to understand the microbial world. Especially genes that reside on mobile genetic elements (MGEs) may spread much more rapidly through a microbial population than genes that reside on the bacterial chromosome. This raises the question: why are some genes associated with MGEs, while others are not? Here, I briefly review a recently proposed class of genes for which we have coined the term "rescuable genes". The fitness effect of carrying these genes is so small, either constantly or on average, that they are prone to be lost from a microbial population. I argue that HGT, even when costly to the individual cells, may play an important role in maintaining these rescuable genes in microbial communities.

Microbial genomic island discovery, visualization and analysis

Horizontal gene transfer (also called lateral gene transfer) is a major mechanism for microbial genome evolution, enabling rapid adaptation and survival in specific niches. Genomic islands (GIs), commonly defined as clusters of bacterial or archaeal genes of probable horizontal origin, are of particular medical, environmental and/or industrial interest, as they disproportionately encode virulence factors and some antimicrobial resistance genes and may harbor entire metabolic pathways that confer a specific adaptation (solvent resistance, symbiosis properties, etc). As large-scale analyses of microbial genomes increases, such as for genomic epidemiology investigations of infectious disease outbreaks in public health, there is increased appreciation of the need to accurately predict and track GIs. Over the past decade, numerous computational tools have been developed to tackle the challenges inherent in accurate GI prediction. We review here the main types of GI prediction methods and discuss their advantages and limitations for a routine analysis of microbial genomes in this era of rapid whole-genome sequencing. An assessment is provided of 20 GI prediction software methods that use sequence-composition bias to identify the GIs, using a reference GI data set from 104 genomes obtained using an independent comparative genomics approach. Finally, we present guidelines to assist researchers in effectively identifying these key genomic regions.

Slightly beneficial genes are retained by bacteria evolving DNA uptake despite selfish elements

Horizontal gene transfer (HGT) and gene loss result in rapid changes in the gene content of bacteria. While HGT aids bacteria to adapt to new environments, it also carries risks such as selfish genetic elements (SGEs). Here, we use modelling to study how HGT of slightly beneficial genes impacts growth rates of bacterial populations, and if bacterial collectives can evolve to take up DNA despite selfish elements. We find four classes of slightly beneficial genes: indispensable, enrichable, rescuable, and unrescuable genes. Rescuable genes - genes with small fitness benefits that are lost from the population without HGT - can be collectively retained by a community that engages in costly HGT. While this 'gene-sharing' cannot evolve in well-mixed cultures, it does evolve in a spatial population like a biofilm. Despite enabling infection by harmful SGEs, the uptake of foreign DNA is evolutionarily maintained by the hosts, explaining the coexistence of bacteria and SGEs.

SaGePhy: an improved phylogenetic simulation framework for gene and subgene evolution

SUMMARY

SaGePhy is a software package for improved phylogenetic simulation of gene and subgene evolution. SaGePhy can be used to generate species trees, gene trees and subgene or (protein) domain trees using a probabilistic birth-death process that allows for gene and subgene duplication, horizontal gene and subgene transfer and gene and subgene loss. SaGePhy implements a range of important features not found in other phylogenetic simulation frameworks/software. These include (i) simulation of subgene or domain level evolution inside one or more gene trees, (ii) simultaneous simulation of both additive and replacing horizontal gene/subgene transfers and (iii) probabilistic sampling of species tree and gene tree nodes, respectively, for gene- and domain-family birth. SaGePhy is open-source, platform independent and written in Java and Python.

AVAILABILITY AND IMPLEMENTATION

Executables, source code (open-source under the revised BSD license) and a detailed manual are freely available from http://compbio.engr.uconn.edu/software/sagephy/.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Molecular Epidemiology, Ecology, and Evolution of Group A Streptococci

The clinico-epidemiological features of diseases caused by group A streptococci (GAS) is presented through the lens of the ecology, population genetics, and evolution of the organism. The serological targets of three typing schemes (M, T, SOF) are themselves GAS cell surface proteins that have a myriad of virulence functions and a diverse array of structural forms. Horizontal gene transfer expands the GAS antigenic cell surface repertoire by generating numerous combinations of M, T, and SOF antigens. However, horizontal gene transfer of the serotype determinant genes is not unconstrained, and therein lies a genetic organization that may signify adaptations to a narrow ecological niche, such as the primary tissue reservoirs of the human host. Adaptations may be further shaped by selection pressures such as herd immunity. Understanding the molecular evolution of GAS on multiple levels-short, intermediate, and long term-sheds insight on mechanisms of host-pathogen interactions, the emergence and spread of new clones, rational vaccine design, and public health interventions.

Analysis of gene gain and loss in the evolution of predatory bacteria

Predatory bacteria are ubiquitously distributed in nature in including in aquatic environments, sewage, intestinal tracts of animals and humans, rhizophere and, soils. However, our understanding of their evolutionary history is limited. Results of recent studies have shown that acquiring novel genes is a major force driving bacterial evolution. Therefore, to gain a better understanding of the impact of gene gain and loss in the evolution of bacterial predators, this study employed comparative genomic approaches to identify core-set gene families and species-specific gene families, and model gene gain and loss events among 11 genomes that represented diverse lineages. In total, 1977 gene families were classified. Of these 509 (pattern 11111111111) were present all of the 11 species. Among the non-core set gene families, 52 were present only in saltwater bacteria predators and had no ortholog in the other genomes. Similarly 109 and 44 were present only in the genomes of Micavibrio spp. and Bdellovibrio spp., respectively. In this study, the gain loss mapping engine GLOOME was selected to analyze and estimate the expectations and probabilities of both gain and loss events in the predatory bacteria. In total, 354 gene families were involved in significant gene gain events, and 407 gene families were classified into gene loss events with high supported value. Moreover, 18 families from the core set gene family were identified as putative genes under positive selection. The results of this study suggest that acquisition of particular genes that encode functional proteins in metabolism and cellular processes and signaling, especially ABC systems, may help bacterial predators adapt to surrounding environmental changes and present different predation strategies for survival in their habitats.

Temporal trends of β-haemolytic streptococcal osteoarticular infections in western Norway

BACKGROUND

Beta-haemolytic streptococci are important contributors to the global burden of osteoarticular infections (OAI). Knowledge on the disease traits specific for streptococcal OAI, however, remains scarce. We wished to explore temporal trends of OAI caused by Group A Streptococci (GAS), Group B Streptococci (GBS) and Group C and G Streptococci (GCGS), and furthermore, to describe the associated host and pathogen characteristics.

METHODS

All cases of microbiologically verified β-haemolytic streptococcal OAI in Health Region Bergen, Norway, in the period 1999-2013 were retrospectively identified. Clinical data were extracted from medical records. Microbial isolates were submitted to antibiotic susceptibility testing and molecular typing.

RESULTS

A total of 24 GAS, 45 GBS and 42 GCGS acute OAI were identified. The cumulative incidence of GCGS OAI, but not GAS or GBS OAI, increased significantly from the first to the last 5-year period (IRR 5.7, p = 0.0003), with the annual incidence peaking at 1.9/100 000 in 2013. GAS OAI generally produced the most acute and severe clinical presentation, whereas GBS and GCGS predominantly affected the elderly, and were significantly associated with the presence of host risk factors of systemic and focal origin, respectively.

CONCLUSIONS

We found a significantly increasing incidence of GCGS OAI, likely related to the presence of host susceptibility factors, including prosthetic material and pre-existing joint disease. With an increasing application of therapeutic and diagnostic bone and joint procedures, the rising trend of OAI caused by GCGS is likely to continue. Sustained epidemiological attentiveness to GCGS seems warranted.

Population genomics reveals additive and replacing horizontal gene transfers in the emerging pathogen Dickeya solani

BACKGROUND

Dickeya solani is an emerging pathogen that causes soft rot and blackleg diseases in several crops including Solanum tuberosum, but little is known about its genomic diversity and evolution.

RESULTS

We combined Illumina and PacBio technologies to complete the genome sequence of D. solani strain 3337 that was used as a reference to compare with 19 other genomes (including that of the type strain IPO2222(T)) which were generated by Illumina technology. This population genomic analysis highlighted an unexpected variability among D. solani isolates since it led to the characterization of two distinct sub-groups within the D. solani species. This approach also revealed different types of variations such as scattered SNP/InDel variations as well as replacing and additive horizontal gene transfers (HGT). Infra-species (between the two D. solani sub-groups) and inter-species (between D. solani and D. dianthicola) replacing HGTs were observed. Finally, this work pointed that genetic and functional variation in the motility trait could contribute to aggressiveness variability in D. solani.

CONCLUSIONS

This work revealed that D. solani genomic variability may be caused by SNPs/InDels as well as replacing and additive HGT events, including plasmid acquisition; hence the D. solani genomes are more dynamic than that were previously proposed. This work alerts on precautions in molecular diagnosis of this emerging pathogen.

Draft Genome Sequences of Four Genetically Distinct Human Isolates of Streptococcus dysgalactiae subsp. equisimilis

β-Hemolytic group C and group G streptococci (GCS-GGS; Streptococcus dysgalactiae subsp. equisimilis) emerged as human pathogens in the late 1970s. We report here the draft genome sequences of four genetically distinct human strains of GCS-GGS isolated between the 1960s and 1980s. Comparative analysis of these genomes may provide a deeper understanding of GCS-GGS genome and virulence evolution.

Molecular epidemiology and genomics of group A Streptococcus

Streptococcus pyogenes (group A Streptococcus; GAS) is a strict human pathogen with a very high prevalence worldwide. This review highlights the genetic organization of the species and the important ecological considerations that impact its evolution. Recent advances are presented on the topics of molecular epidemiology, population biology, molecular basis for genetic change, genome structure and genetic flux, phylogenomics and closely related streptococcal species, and the long- and short-term evolution of GAS. The application of whole genome sequence data to addressing key biological questions is discussed.

Improved gene tree error correction in the presence of horizontal gene transfer

MOTIVATION

The accurate inference of gene trees is a necessary step in many evolutionary studies. Although the problem of accurate gene tree inference has received considerable attention, most existing methods are only applicable to gene families unaffected by horizontal gene transfer. As a result, the accurate inference of gene trees affected by horizontal gene transfer remains a largely unaddressed problem.

RESULTS

In this study, we introduce a new and highly effective method for gene tree error correction in the presence of horizontal gene transfer. Our method efficiently models horizontal gene transfers, gene duplications and losses, and uses a statistical hypothesis testing framework [Shimodaira-Hasegawa (SH) test] to balance sequence likelihood with topological information from a known species tree. Using a thorough simulation study, we show that existing phylogenetic methods yield inaccurate gene trees when applied to horizontally transferred gene families and that our method dramatically improves gene tree accuracy. We apply our method to a dataset of 11 cyanobacterial species and demonstrate the large impact of gene tree accuracy on downstream evolutionary analyses.

AVAILABILITY AND IMPLEMENTATION

An implementation of our method is available at http://compbio.mit.edu/treefix-dtl/

CONTACT

: mukul@engr.uconn.edu or manoli@mit.edu

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Beta-haemolytic group A, C and G streptococcal infections in Western Norway: a 15-year retrospective survey

Pyogenic streptococci cause significant morbidity and mortality, and the incidence of invasive group C and G streptococcal disease appears to be increasing. In this retrospective study we describe the epidemiological characteristics of invasive group A, C and G, along with non-invasive group C and G streptococcal infections in Western Norway from 1999 to 2013. A total of 512 invasive streptococcal infections were identified, of these 297 (58%) were group A (GAS), 24 (5%) group C (GCS) and 188 (37%) group G streptococci (GGS). In the non-invasive group, 4935 GCS and GGS-infections were identified. GCS and GGS were treated as one group (GCGS) for statistical purposes. All microbial categories displayed increasing incidence with age, seasonal variation and a male predominance. The incidence of invasive GCGS infections increased significantly from 1.4/100,000 inhabitants in 1999 to 6.3/100,000 in 2013 (p <0.001). Conversely, the annual rates of invasive GAS infection exhibited marked fluctuations, ranging from 2.7/100,000 (2000) to 8.3/100,000 (1999), but no significant temporal trends were observed. The incidence of non-invasive GCGS infections decreased significantly during the study period (p <0.001). The most frequently encountered emm-types among the 209 iGAS-isolates analysed were emm1 (24%), emm3 (14%) and emm28 (14%); whereas stG643 (19%), stG485 (15%) and stG6 (13%) were most prevalent among the 122 iGCGS-isolates available for typing. The increasing burden of invasive β-haemolytic streptococcal disease in our community calls for sustained attentiveness to the clinical and molecular aspects of GAS, GCS and GGS infections.

Genomic heterogeneity and ecological speciation within one subspecies of Bacillus subtilis

Closely related bacterial genomes usually differ in gene content, suggesting that nearly every strain in nature may be ecologically unique. We have tested this hypothesis by sequencing the genomes of extremely close relatives within a recognized taxon and analyzing the genomes for evidence of ecological distinctness. We compared the genomes of four Death Valley isolates plus the laboratory strain W23, all previously classified as Bacillus subtilis subsp. spizizenii and hypothesized through multilocus analysis to be members of the same ecotype (an ecologically homogeneous population), named putative ecotype 15 (PE15). These strains showed a history of positive selection on amino acid sequences in 38 genes. Each of the strains was under a different regimen of positive selection, suggesting that each strain is ecologically unique and represents a distinct ecological speciation event. The rate of speciation appears to be much faster than can be resolved with multilocus sequencing. Each PE15 strain contained unique genes known to confer a function for bacteria. Remarkably, no unique gene conferred a metabolic system or subsystem function that was not already present in all the PE15 strains sampled. Thus, the origin of ecotypes within this clade shows no evidence of qualitative divergence in the set of resources utilized. Ecotype formation within this clade is consistent with the nanoniche model of bacterial speciation, in which ecotypes use the same set of resources but in different proportions, and genetic cohesion extends beyond a single ecotype to the set of ecotypes utilizing the same resources.

Genetic diversity and virulence properties of Streptococcus dysgalactiae subsp. equisimilis from different sources

A recent increase in virulence of pathogenic Streptococcus dysgalactiae subsp. equisimilis (SDSE) has been widely proposed. Such an increase may be partly explained by the acquisition of new virulence traits by horizontal gene transfer from related streptococci such as Streptococcus pyogenes (GAS) and Streptococcus agalactiae (GBS). A collection of 54 SDSE strains isolated in Italy in the years 2000–2010 from different sources (paediatric throat carriage, invasive and non-invasive diseases) was characterized by emm typing and pulsed-field gel electrophoresis (PFGE) analysis. The virulence repertoire was evaluated by PCR for the presence of GAS superantigen (spe) genes, the streptolysin S (sagA) gene, the group G fibronectin-binding protein (gfbA) gene and GAS–GBS alpha-like protein family (alp) genes; moreover, the ability to invade human epithelial cells was investigated. Resistance to tetracycline, erythromycin and clindamycin was assessed. The combined use of emm typing and PFGE proved to be a reliable strategy for the epidemiological analysis of SDSE isolates. The most frequent emm types were the same as those more frequently reported in other studies, thus indicating the diffusion of a limited number of a few successful emm types fit to disseminate in humans. The speG gene was detected in SDSE strains of different genetic backgrounds. Erythromycin resistance determined by the erm(T) gene, and the unusual, foggy MLSB phenotype, observed in one and seven strains, respectively, have never previously, to our knowledge, been reported in SDSE. Moreover, a new member of the alp family was identified. The identification of new antibiotic and virulence determinants, despite the small size of the sample analysed, shows the importance of constant attention to monitoring the extent of lateral gene transfer in this emerging pathogen.

A molecular study on bacterial resistance to arsenic-toxicity in surface and underground waters of Latium (Italy)

Latium, a region in central Italy, is known for its extensive volcanic areas that make a significant contribution to the arsenic (As) contamination of freshwater environments, even though some degree of As water pollution may be caused by human activities. The information available on indigenous As-resistant prokaryotes in aquatic environments of Latium is, however, still limited. In this study, we describe new bacteria that are resistant to arsenic toxicity and were isolated from the surface waters of Lake Vico and the Sacco River, two groundwater systems in Latium, as well as from bottled natural mineral water from the same region. The 16S rRNA gene sequence analysis for the As-resistant strains in lake and river waters points to a prevalence of β- and γ-Proteobacteria, while α-Proteobacteria, Firmicutes and Bacteroidetes are represented to a lesser extent. By contrast, solely γ-Proteobacteria were isolated from groundwater samples. The presence of Actinobacteria was documented exclusively in bottled mineral water. In addition, we conducted a DNA sequence-based study on the gene codifying arsB, an As(III) efflux membrane protein pump related to arsenic resistance, for all the As-resistant bacterial isolates. A phylogenetic analysis was carried out on the newly sequenced 16S rRNA genes and arsB in the present study as well as on an additional 16S rRNA/arsB dataset we obtained previously from Lake Albano, from the Tiber and from a well in Bassano Romano located in Latium (Davolos and Pietrangeli, 2011). Overall, the phylogenetic diversity of As-resistant bacteria in underground water was very limited if compared with lentic and lotic waters. Lastly, our molecular data support the hypothesis that the horizontal gene transfer of ars in As-containing freshwater environments is not limited to closely-related genomes, but also occurs between bacteria that are distant from an evolutionary viewpoint, thereby indicating that such genetic events may be considered a source of microbial resistance to arsenic-toxicity.