Evolutionary Processes in the Emergence and Recent Spread of the Syphilis Agent, Treponema pallidum.

The incidence of syphilis has risen worldwide in the last decade in spite of being an easily treated infection. The causative agent of this sexually transmitted disease is the bacterium Treponema pallidum subspecies pallidum (TPA), very closely related to subsp. pertenue (TPE) and endemicum (TEN), responsible for the human treponematoses yaws and bejel, respectively. Although much focus has been placed on the question of the spatial and temporary origins of TPA, the processes driving the evolution and epidemiological spread of TPA since its divergence from TPE and TEN are not well understood. Here, we investigate the effects of recombination and selection as forces of genetic diversity and differentiation acting during the evolution of T. pallidum subspecies. Using a custom-tailored procedure, named phylogenetic incongruence method, with 75 complete genome sequences, we found strong evidence for recombination among the T. pallidum subspecies, involving 12 genes and 21 events. In most cases, only one recombination event per gene was detected and all but one event corresponded to intersubspecies transfers, from TPE/TEN to TPA. We found a clear signal of natural selection acting on the recombinant genes, which is more intense in their recombinant regions. The phylogenetic location of the recombination events detected and the functional role of the genes with signals of positive selection suggest that these evolutionary processes had a key role in the evolution and recent expansion of the syphilis bacteria and significant implications for the selection of vaccine candidates and the design of a broadly protective syphilis vaccine.

Molecular characterization of Treponema pallidum subsp. pallidum in Switzerland and France with a new multilocus sequence typing scheme.

Syphilis is an important public health problem and an increasing incidence has been noted in recent years. Characterization of strain diversity through molecular data plays a critical role in the epidemiological understanding of this re-emergence. We here propose a new high-resolution multilocus sequence typing (MLST) scheme for Treponema pallidum subsp. pallidum (TPA). We analyzed 30 complete and draft TPA genomes obtained directly from clinical samples or from rabbit propagated strains to identify suitable typing loci and tested the new scheme on 120 clinical samples collected in Switzerland and France. Our analyses yielded three loci with high discriminatory power: TP0136, TP0548, and TP0705. Together with analysis of the 23S rRNA gene mutations for macrolide resistance, we propose these loci as MLST for TPA. Among clinical samples, 23 allelic profiles as well as a high percentage (80% samples) of macrolide resistance were revealed. The new MLST has higher discriminatory power compared to previous typing schemes, enabling distinction of TPA from other treponemal bacteria, distinction between the two main TPA clades (Nichols and SS14), and differentiation of strains within these clades.

Origin of modern syphilis and emergence of a pandemic Treponema pallidum cluster.

The abrupt onslaught of the syphilis pandemic that started in the late fifteenth century established this devastating infectious disease as one of the most feared in human history1. Surprisingly, despite the availability of effective antibiotic treatment since the mid-twentieth century, this bacterial infection, which is caused by Treponema pallidum subsp. pallidum (TPA), has been re-emerging globally in the last few decades with an estimated 10.6 million cases in 2008 (ref. 2). Although resistance to penicillin has not yet been identified, an increasing number of strains fail to respond to the second-line antibiotic azithromycin3. Little is known about the genetic patterns in current infections or the evolutionary origins of the disease due to the low quantities of treponemal DNA in clinical samples and difficulties in cultivating the pathogen4. Here, we used DNA capture and whole-genome sequencing to successfully interrogate genome-wide variation from syphilis patient specimens, combined with laboratory samples of TPA and two other subspecies. Phylogenetic comparisons based on the sequenced genomes indicate that the TPA strains examined share a common ancestor after the fifteenth century, within the early modern era. Moreover, most contemporary strains are azithromycin-resistant and are members of a globally dominant cluster, named here as SS14-Ω. The cluster diversified from a common ancestor in the mid-twentieth century subsequent to the discovery of antibiotics. Its recent phylogenetic divergence and global presence point to the emergence of a pandemic strain cluster.

Transition from one- to two-dimensional development facilitates maintenance of multicellularity.

Filamentous organisms represent an example where incomplete separation after cell division underlies the development of multicellular formations. With a view to understanding the evolution of more complex multicellular structures, we explore the transition of multicellular growth from one to two dimensions. We develop a computational model to simulate multicellular development in populations where cells exhibit density-dependent division and death rates. In both the one- and two-dimensional contexts, multicellular formations go through a developmental cycle of growth and subsequent decay. However, the model shows that a transition to a higher dimension increases the size of multicellular formations and facilitates the maintenance of large cell clusters for significantly longer periods of time. We further show that the turnover rate for cell division and death scales with the number of iterations required to reach the stationary multicellular size at equilibrium. Although size and life cycles of multicellular organisms are affected by other environmental and genetic factors, the model presented here evaluates the extent to which the transition of multicellular growth from one to two dimensions contributes to the maintenance of multicellular structures during development.

Quarantine as a public health measure against an emerging infectious disease: syphilis in Zurich at the dawn of the modern era (1496-1585).

Syphilis is considered as one of the most devastating sexually transmitted diseases in human history. Based on historical records, the "Böse Blattern" (German for "Evil Pocks") spread through Europe after 1495 and shared symptoms with what we know today as syphilis. Many cities took measures to protect their population. Here, transliterations of archival documents from the 15(th) and 16(th) century (provided in the appendix) are used to trace the steps taken by the governing authorities in Zurich to deal with this emerging infectious disease. One of the central measures taken by the city was to establish a quarantine facility referred to as the "Blatternhaus". The city doctors, including the well-known physician and naturalist Conrad Gessner, oversaw the obligatory quarantine and treatment of patients with symptoms. Treatment could range from better nutrition, herbal remedies and skin ointments to aggressive heat therapy and "smoking". Furthermore, the affliction was suspected as a sexually acquired disease, hence prostitutes and infected foreigners were extradited from the city. Meanwhile, the church used its social influence to promote a more "Christian" behavior. In several respects, the public health measures taken in the 16(th) century against a new and menacing epidemic do not diverge in their basic rationale from approaches used during the 20(th) century and today.

Dynamics of a Tularemia Outbreak in a Closely Monitored Free-Roaming Population of Wild House Mice.

Infectious disease outbreaks can be devastating because of their sudden occurrence, as well as the complexity of monitoring and controlling them. Outbreaks in wildlife are even more challenging to observe and describe, especially when small animals or secretive species are involved. Modeling such infectious disease events is relevant to investigating their dynamics and is critical for decision makers to accomplish outbreak management. Tularemia, caused by the bacterium Francisella tularensis, is a potentially lethal zoonosis. Of the few animal outbreaks that have been reported in the literature, only those affecting zoo animals have been closely monitored. Here, we report the first estimation of the basic reproduction number R0 of an outbreak in wildlife caused by F. tularensis using quantitative modeling based on a susceptible-infected-recovered framework. We applied that model to data collected during an extensive investigation of an outbreak of tularemia caused by F. tularensis subsp. holarctica (also designated as type B) in a closely monitored, free-roaming house mouse (Mus musculus domesticus) population in Switzerland. Based on our model and assumptions, the best estimated basic reproduction number R0 of the current outbreak is 1.33. Our results suggest that tularemia can cause severe outbreaks in small rodents. We also concluded that the outbreak self-exhausted in approximately three months without administrating antibiotics.

Chromothripsis-like patterns are recurring but heterogeneously distributed features in a survey of 22,347 cancer genome screens.

BACKGROUND: Chromothripsis is a recently discovered phenomenon of genomic rearrangement, possibly arising during a single genome-shattering event. This could provide an alternative paradigm in cancer development, replacing the gradual accumulation of genomic changes with a "one-off" catastrophic event. However, the term has been used with varying operational definitions, with the minimal consensus being a large number of locally clustered copy number aberrations. The mechanisms underlying these chromothripsis-like patterns (CTLP) and their specific impact on tumorigenesis are still poorly understood. RESULTS: Here, we identified CTLP in 918 cancer samples, from a dataset of more than 22,000 oncogenomic arrays covering 132 cancer types. Fragmentation hotspots were found to be located on chromosome 8, 11, 12 and 17. Among the various cancer types, soft-tissue tumors exhibited particularly high CTLP frequencies. Genomic context analysis revealed that CTLP rearrangements frequently occurred in genomes that additionally harbored multiple copy number aberrations (CNAs). An investigation into the affected chromosomal regions showed a large proportion of arm-level pulverization and telomere related events, which would be compatible to a number of underlying mechanisms. We also report evidence that these genomic events may be correlated with patient age, stage and survival rate. CONCLUSIONS: Through a large-scale analysis of oncogenomic array data sets, this study characterized features associated with genomic aberrations patterns, compatible to the spectrum of "chromothripsis"-definitions as previously used. While quantifying clustered genomic copy number aberrations in cancer samples, our data indicates an underlying biological heterogeneity behind these chromothripsis-like patterns, beyond a well defined "chromthripsis" phenomenon.

Phenotypic diversity of multicellular filamentation in oral Streptococci.

Filamentous multicellular bacteria are among the most ancient multicellular organisms. They inhabit a great variety of environments and are present in the human body, including the oral cavity. Beside the selective advantages related to the larger size achieved through filamentation, the development of multicellular bacteria can be also driven by simple ecological factors such as birth and death rates at the cellular level. In order to extend earlier results obtained in aquatic species, we investigate the filamentation process of four different strains of oral streptococci, namely S. mutans, S. salivarius, S. oralis and S. anginosus. The results indicate differences in the capacities of different streptococcus species to form filaments, manifested in terms of length and the time-scale of filament elongation. The filamentation pattern of these oral streptococci resembles that of aquatic bacteria, whereby filaments reach a peak length during exponential growth and become short when the population reaches a steady state. Hence, this study validates that multicellularity can be an emergent property of filamentous bacteria of different ecological niches, and that phenotypic differences in filamentation can occur within species of the same genus, in this case oral streptococci. Moreover, given the role that specific oral streptococci can play in the etiology of oral diseases, these results can possibly open new perspectives in the study of the virulence properties of these species.

Correction: Use of a Four-Tiered Graph to Parse the Factors Leading to Phenotypic Clustering in Bacteria: A Case Study Based on Samples from the Aletsch Glacier.

[This corrects the article on p. e65059 in vol. 8.].

Use of a four-tiered graph to parse the factors leading to phenotypic clustering in bacteria: a case study based on samples from the Aletsch Glacier.

An understanding of bacterial diversity and evolution in any environment requires knowledge of phenotypic diversity. In this study, the underlying factors leading to phenotypic clustering were analyzed and interpreted using a novel approach based on a four-tiered graph. Bacterial isolates were organized into equivalence classes based on their phenotypic profile. Likewise, phenotypes were organized in equivalence classes based on the bacteria that manifest them. The linking of these equivalence classes in a four-tiered graph allowed for a quick visual identification of the phenotypic measurements leading to the clustering patterns deduced from principal component analyses. For evaluation of the method, we investigated phenotypic variation in enzyme production and carbon assimilation of members of the genera Pseudomonas and Serratia, isolated from the Aletsch Glacier in Switzerland. The analysis indicates that the genera isolated produce at least six common enzymes and can exploit a wide range of carbon resources, though some specialist species within the pseudomonads were also observed. We further found that pairwise distances between enzyme profiles strongly correlate with distances based on carbon profiles. However, phenotypic distances weakly correlate with phylogenetic distances. The method developed in this study facilitates a more comprehensive understanding of phenotypic clustering than what would be deduced from principal component analysis alone.

Evolution of multicellularity coincided with increased diversification of cyanobacteria and the Great Oxidation Event.

Cyanobacteria are among the most diverse prokaryotic phyla, with morphotypes ranging from unicellular to multicellular filamentous forms, including those able to terminally (i.e., irreversibly) differentiate in form and function. It has been suggested that cyanobacteria raised oxygen levels in the atmosphere around 2.45-2.32 billion y ago during the Great Oxidation Event (GOE), hence dramatically changing life on the planet. However, little is known about the temporal evolution of cyanobacterial lineages, and possible interplay between the origin of multicellularity, diversification of cyanobacteria, and the rise of atmospheric oxygen. We estimated divergence times of extant cyanobacterial lineages under Bayesian relaxed clocks for a dataset of 16S rRNA sequences representing the entire known diversity of this phylum. We tested whether the evolution of multicellularity overlaps with the GOE, and whether multicellularity is associated with significant shifts in diversification rates in cyanobacteria. Our results indicate an origin of cyanobacteria before the rise of atmospheric oxygen. The evolution of multicellular forms coincides with the onset of the GOE and an increase in diversification rates. These results suggest that multicellularity could have played a key role in triggering cyanobacterial evolution around the GOE.

Renaissance model of an epidemic with quarantine.

Quarantine is one possible solution to limit the propagation of an emerging infectious disease. Typically, infected individuals are removed from the population by avoiding physical contact with healthy individuals. A key factor for the success of a quarantine strategy is the carrying capacity of the facility. This is often a known parameter, while other parameters such as those defining the population structure are more difficult to assess. Here we develop a model where we explicitly introduce the carrying capacity of the quarantine facility into a susceptible-infected-recovered (SIR) framework. We show how the model can address the propagation and control of contact and sexually transmitted infections. We illustrate this by a case study of the city of Zurich during the 16th century, when it had to face an epidemic of syphilis. After Swiss mercenaries came back from a war in Naples in 1495, the authorities of the city addressed subsequent epidemics by, among others, placing infected members of the population in quarantine. Our results suggest that a modestly sized quarantine facility can successfully prevent or reduce an epidemic. However, false detection can present a real impediment for this solution. Indiscriminate quarantine of individuals can lead to the overfilling of the facility, and prevent the intake of infected individuals. This results in the failure of the quarantine policy. Hence, improving the rate of true over false detection becomes the key factor for quarantine strategies. Moreover, in the case of sexually transmitted infections, asymmetries in the male to female ratio, and the force of infection pertaining to each sex and class of sexual encounter can alter the effectiveness of quarantine measures. For example, a heterosexually transmitted disease that mainly affects one sex is harder to control in a population with more individuals of the opposite sex. Hence an imbalance in the sex ratios as seen in situations such as mining colonies, or populations at war, can present impediments for the success of quarantine policies.

Gene copy number variation and its significance in cyanobacterial phylogeny.

BACKGROUND: In eukaryotes, variation in gene copy numbers is often associated with deleterious effects, but may also have positive effects. For prokaryotes, studies on gene copy number variation are rare. Previous studies have suggested that high numbers of rRNA gene copies can be advantageous in environments with changing resource availability, but further association of gene copies and phenotypic traits are not documented. We used one of the morphologically most diverse prokaryotic phyla to test whether numbers of gene copies are associated with levels of cell differentiation. RESULTS: We implemented a search algorithm that identified 44 genes with highly conserved copies across 22 fully sequenced cyanobacterial taxa. For two very basal cyanobacterial species, Gloeobacter violaceus and a thermophilic Synechococcus species, distinct phylogenetic positions previously found were supported by identical protein coding gene copy numbers. Furthermore, we found that increased ribosomal gene copy numbers showed a strong correlation to cyanobacteria capable of terminal cell differentiation. Additionally, we detected extremely low variation of 16S rRNA sequence copies within the cyanobacteria. We compared our results for 16S rRNA to three other eubacterial phyla (Chroroflexi, Spirochaetes and Bacteroidetes). Based on Bayesian phylogenetic inference and the comparisons of genetic distances, we could confirm that cyanobacterial 16S rRNA paralogs and orthologs show significantly stronger conservation than found in other eubacterial phyla. CONCLUSIONS: A higher number of ribosomal operons could potentially provide an advantage to terminally differentiated cyanobacteria. Furthermore, we suggest that 16S rRNA gene copies in cyanobacteria are homogenized by both concerted evolution and purifying selection. In addition, the small ribosomal subunit in cyanobacteria appears to evolve at extraordinary slow evolutionary rates, an observation that has been made previously for morphological characteristics of cyanobacteria.

Genome sequence of the filamentous bacterium Fibrisoma limi BUZ 3T.

Fibrisoma limi strain BUZ 3(T), a Gram-negative bacterium, was isolated from coastal mud from the North Sea (Fedderwardersiel, Germany) and characterized using a polyphasic approach in 2011. The genome consists of a chromosome of about 7.5 Mb and three plasmids.

Genome sequence of Fibrella aestuarina BUZ 2(T), a filamentous marine bacterium.

Fibrella aestuarina BUZ 2(T) is the type strain of the recently characterized genus Fibrella. Here we report the draft genome sequence of this strain, which consists of a single scaffold representing the chromosome (with 11 gaps) and a 161-kb circular plasmid.

Advantages of the division of labour for the long-term population dynamics of cyanobacteria at different latitudes.

A fundamental advancement in the evolution of complexity is division of labour. This implies a partition of tasks among cells, either spatially through cellular differentiation, or temporally via a circadian rhythm. Cyanobacteria often employ either spatial differentiation or a circadian rhythm in order to separate the chemically incompatible processes of nitrogen fixation and photosynthesis. We present a theoretical framework to assess the advantages in terms of biomass production and population size for three species types: terminally differentiated (heterocystous), circadian, and an idealized species in which nitrogen and carbon fixation occur without biochemical constraints. On the basis of real solar irradiance data at different latitudes, we simulate population dynamics in isolation and in competition for light over a period of 40 years. Our results show that in isolation and regardless of latitude, the biomass of heterocystous cyanobacteria that optimally invest resources is comparable to that of the idealized unconstrained species. Hence, spatial division of labour overcomes biochemical constraints and enhances biomass production. In the circadian case, the strict temporal task separation modelled here hinders high biomass production in comparison with the heterocystous species. However, circadian species are found to be successful in competition for light whenever their resource investment prevents a waste of fixed nitrogen more effectively than do heterocystous species. In addition, we show the existence of a trade-off between population size and biomass accumulation, whereby each species can optimally invest resources to be proficient in biomass production or population growth, but not necessarily both. Finally, the model produces chaotic dynamics for population size, which is relevant to the study of cyanobacterial blooms.

Differences in cell division rates drive the evolution of terminal differentiation in microbes.

Multicellular differentiated organisms are composed of cells that begin by developing from a single pluripotent germ cell. In many organisms, a proportion of cells differentiate into specialized somatic cells. Whether these cells lose their pluripotency or are able to reverse their differentiated state has important consequences. Reversibly differentiated cells can potentially regenerate parts of an organism and allow reproduction through fragmentation. In many organisms, however, somatic differentiation is terminal, thereby restricting the developmental paths to reproduction. The reason why terminal differentiation is a common developmental strategy remains unexplored. To understand the conditions that affect the evolution of terminal versus reversible differentiation, we developed a computational model inspired by differentiating cyanobacteria. We simulated the evolution of a population of two cell types -nitrogen fixing or photosynthetic- that exchange resources. The traits that control differentiation rates between cell types are allowed to evolve in the model. Although the topology of cell interactions and differentiation costs play a role in the evolution of terminal and reversible differentiation, the most important factor is the difference in division rates between cell types. Faster dividing cells always evolve to become the germ line. Our results explain why most multicellular differentiated cyanobacteria have terminally differentiated cells, while some have reversibly differentiated cells. We further observed that symbioses involving two cooperating lineages can evolve under conditions where aggregate size, connectivity, and differentiation costs are high. This may explain why plants engage in symbiotic interactions with diazotrophic bacteria.

The effect of polyandry on a distorter system with differential viabilities in the sexes.

The presence of selfish genetic elements can have fatal consequences for populations that harbor them. In the well known t haplotype in wild house mice, large proportions of the population die from t/t recessive lethal effects. Due to strong advantages at the gamete level (drive), t haplotypes nevertheless occur at substantial frequencies. The stable presence of a lethal is not the only effect of the t. It also distorts the fate of mutations that differentially affect male and female survival and reproduction (such as in sexual conflict), by giving male selective effects a strong advantage over female selective effects. In a recent study, we proposed polyandry as a potential counterstrategy against t deleterious effects. Here, we show that (1) the efficiency of polyandry in reducing the t frequency strongly depends on the selective context and (2) polyandry helps to reduce male-biased leverage in sex dependent selection.

Transgenerational propagation and quantitative maintenance of paternal centromeres depends on Cid/Cenp-A presence in Drosophila sperm.

In Drosophila melanogaster, as in many animal and plant species, centromere identity is specified epigenetically. In proliferating cells, a centromere-specific histone H3 variant (CenH3), named Cid in Drosophila and Cenp-A in humans, is a crucial component of the epigenetic centromere mark. Hence, maintenance of the amount and chromosomal location of CenH3 during mitotic proliferation is important. Interestingly, CenH3 may have different roles during meiosis and the onset of embryogenesis. In gametes of Caenorhabditis elegans, and possibly in plants, centromere marking is independent of CenH3. Moreover, male gamete differentiation in animals often includes global nucleosome for protamine exchange that potentially could remove CenH3 nucleosomes. Here we demonstrate that the control of Cid loading during male meiosis is distinct from the regulation observed during the mitotic cycles of early embryogenesis. But Cid is present in mature sperm. After strong Cid depletion in sperm, paternal centromeres fail to integrate into the gonomeric spindle of the first mitosis, resulting in gynogenetic haploid embryos. Furthermore, after moderate depletion, paternal centromeres are unable to re-acquire normal Cid levels in the next generation. We conclude that Cid in sperm is an essential component of the epigenetic centromere mark on paternal chromosomes and it exerts quantitative control over centromeric Cid levels throughout development. Hence, the amount of Cid that is loaded during each cell cycle appears to be determined primarily by the preexisting centromeric Cid, with little flexibility for compensation of accidental losses.

Evolution of cyanobacterial morphotypes: Taxa required for improved phylogenomic approaches.

Within prokaryotes cyanobacteria represent one of the oldest and morphologically most diverse phyla on Earth. The rise of oxygen levels in the atmosphere 2.32-2.45 billion years ago is assigned to the photosynthetic activity of ancestors from this phylum. Subsequently cyanobacteria were able to adapt to various habitats evolving a comprehensive set of different morphotypes. In a recent study we showed that this evolution is not a gradual transition from simple unicellular to more complex multicellular forms as often assumed. Instead complexity was lost several times and regained at least once. An understanding of the genetic basis of these transitions would be further strengthened by phylogenomic approaches. However, considering that new methods for phylogenomic analyses are emerging, it is unfortunate that genomes available today are comprised of an unbalanced sampling of taxa. We propose avenues to remedy this by identifying taxa that would improve the representation of phylogenetic diversity in this phylum.