16S rDNA pyrosequencing analysis of bacterial community in heavy metals polluted soils.

Soil contamination with heavy metals is a widespread problem, especially prominent on grounds lying in the vicinity of mines, smelters, and other industrial facilities. Many such areas are located in Southern Poland; they are polluted mainly with Pb, Zn, Cd, or Cu, and locally also with Cr. As for now, little is known about most bacterial species thriving in such soils and even less about a core bacterial community--a set of taxa common to polluted soils. Therefore, we wanted to answer the question if such a set could be found in samples differing physicochemically and phytosociologically. To answer the question, we analyzed bacterial communities in three soil samples contaminated with Pb and Zn and two contaminated with Cr and lower levels of Pb and Zn. The communities were assessed with 16S rRNA gene fragments pyrosequencing. It was found that the samples differed significantly and Zn decreased both diversity and species richness at species and family levels, while plant species richness did not correlate with bacterial diversity. In spite of the differences between the samples, they shared many operational taxonomic units (OTUs) and it was possible to delineate the core microbiome of our sample set. The core set of OTUs comprised members of such taxa as Sphingomonas, Candidatus Solibacter, or Flexibacter showing that particular genera might be shared among sites ~40 km distant.

Molecular evolution in court: analysis of a large hepatitis C virus outbreak from an evolving source.

BACKGROUND: Molecular phylogenetic analyses are used increasingly in the epidemiological investigation of outbreaks and transmission cases involving rapidly evolving RNA viruses. Here, we present the results of such an analysis that contributed to the conviction of an anesthetist as being responsible for the infection of 275 of his patients with hepatitis C virus. RESULTS: We obtained sequences of the NS5B and E1-E2 regions in the viral genome for 322 patients suspected to have been infected by the doctor, and for 44 local, unrelated controls. The analysis of 4,184 cloned sequences of the E1-E2 region allowed us to exclude 47 patients from the outbreak. A subset of patients had known dates of infection. We used these data to calibrate a relaxed molecular clock and to determine a rough estimate of the time of infection for each patient. A similar analysis led to an estimate for the time of infection of the source. The date turned out to be 10 years before the detection of the outbreak. The number of patients infected was small at first, but it increased substantially in the months before the detection of the outbreak. CONCLUSIONS: We have developed a procedure to integrate molecular phylogenetic reconstructions of rapidly evolving viral populations into a forensic setting adequate for molecular epidemiological analysis of outbreaks and transmission events. We applied this procedure to a large outbreak of hepatitis C virus caused by a single source and the results obtained played a key role in the trial that led to the conviction of the suspected source.

Influence of the Escherichia coli oxyR gene function on lambda prophage maintenance.

In Escherichia coli hosts, hydrogen peroxide is one of the factors that may cause induction of lambda prophage. Here, we demonstrate that H2O2-mediated lambda prophage induction is significantly enhanced in the oxyR mutant host. The mRNA levels for cI gene expression were increased in a lambda lysogen in the presence of H2O2. On the other hand, stimulation of the p(M) promoter by cI857 overproduced from a multicopy plasmid was decreased in the DeltaoxyR mutant in the presence of H2O2 but not under normal growth conditions. The purified OxyR protein did bind specifically to the p(M) promoter region. This binding impaired efficiency of interaction of the cI protein with the OR3 site, while stimulating such a binding to OR2 and OR1 sites, in the regulatory region of the p(M) promoter. We propose that changes in cI gene expression, perhaps in combination with moderately induced SOS response, may be responsible for enhanced lambda prophage induction by hydrogen peroxide in the oxyR mutant. Therefore, OxyR seems to be a factor stimulating lambda prophage maintenance under conditions of oxidative stress. This proposal is discussed in the light of efficiency of induction of lambdoid prophages bearing genes coding for Shiga toxins.

Using evolutionary tools to refine the new hypervariable region 3 within the envelope 2 protein of hepatitis C virus.

The envelope 2 protein of hepatitis C virus (HCV) presents three hypervariable regions, named HVR1, HVR2 and HVR3, in which the presence of antigenic sites has been described. Genetic variability in these regions may reflect the generation of escape mutants as a consequence of the immune response. Therefore, these regions would tend to accumulate amino acid changes along the infection process, an effect that could be accelerated by antiviral treatments. In this study, we have analyzed the E1-E2 region of 23 HCV patients non-responders to antiviral treatment, 7 of which were infected with subtype 1a, 15 with subtype 1b, and 1 with a new HCV-1 subtype, before and after 6 and/or 12 months of peg-interferon+ribavirin treatment. We have sequenced about 100 clones from each sample, analyzing a total of 4906 sequences. A detailed analysis of the evolutionary forces acting along the genome region studied confirmed the existence of the three hypervariable regions, characterized by significant changes in amino acid composition between samples taken at different times from the same patient and a high number of sites evolving under positive selection. Moreover, for the recently described HVR3, our results suggest that its location could be restricted to residues 434-450, instead of the originally postulated 431-466.

Statistical measures of uncertainty for branches in phylogenetic trees inferred from molecular sequences by using model-based methods.

In recent years, the emphasis of theoretical work on phylogenetic inference has shifted from the development of new tree inference methods to the development of methods to measure the statistical support for the topologies. This paper reviews 3 approaches to assign support values to branches in trees obtained in the analysis of molecular sequences: the bootstrap, the Bayesian posterior probabilities for clades, and the interior branch tests. In some circumstances, these methods give different answers. It should not be surprising: their assumptions are different. Thus the interior branch tests assume that a given topology is true and only consider if a particular branch length is longer than zero. If a tree is incorrect, a wrong branch (a low bootstrap or Bayesian support may be an indication) may have a non-zero length. If the substitution model is oversimplified, the length of a branch may be overestimated, and the Bayesian support for the branch may be inflated. The bootstrap, on the other hand, approximates the variance of the data under the real model of sequence evolution, because it involves direct resampling from this data. Thus the discrepancy between the Bayesian support and the bootstrap support may signal model inaccuracy. In practical application, use of all 3 methods is recommended, and if discrepancies are observed, then a careful analysis of their potential origins should be made.

Topology testing of phylogenies using least squares methods.

BACKGROUND: The least squares (LS) method for constructing confidence sets of trees is closely related to LS tree building methods, in which the goodness of fit of the distances measured on the tree (patristic distances) to the observed distances between taxa is the criterion used for selecting the best topology. The generalized LS (GLS) method for topology testing is often frustrated by the computational difficulties in calculating the covariance matrix and its inverse, which in practice requires approximations. The weighted LS (WLS) allows for a more efficient albeit approximate calculation of the test statistic by ignoring the covariances between the distances. RESULTS: The goal of this paper is to assess the applicability of the LS approach for constructing confidence sets of trees. We show that the approximations inherent to the WLS method did not affect negatively the accuracy and reliability of the test both in the analysis of biological sequences and DNA-DNA hybridization data (for which character-based testing methods cannot be used). On the other hand, we report several problems for the GLS method, at least for the available implementation. For many data sets of biological sequences, the GLS statistic could not be calculated. For some data sets for which it could, the GLS method included all the possible trees in the confidence set despite a strong phylogenetic signal in the data. Finally, contrary to WLS, for simulated sequences GLS showed undercoverage (frequent non-inclusion of the true tree in the confidence set). CONCLUSION: The WLS method provides a computationally efficient approximation to the GLS useful especially in exploratory analyses of confidence sets of trees, when assessing the phylogenetic signal in the data, and when other methods are not available.

Biodiversity of bacteriophages: morphological and biological properties of a large group of phages isolated from urban sewage.

A large scale analysis presented in this article focuses on biological and physiological variety of bacteriophages. A collection of 83 bacteriophages, isolated from urban sewage and able to propagate in cells of different bacterial hosts, has been obtained (60 infecting Escherichia coli, 10 infecting Pseudomonas aeruginosa, 4 infecting Salmonella enterica, 3 infecting Staphylococcus sciuri, and 6 infecting Enterococcus faecalis). High biological diversity of the collection is indicated by its characteristics, both morphological (electron microscopic analyses) and biological (host range, plaque size and morphology, growth at various temperatures, thermal inactivation, sensitivity to low and high pH, sensitivity to osmotic stress, survivability upon treatment with organic solvents and detergents), and further supported by hierarchical cluster analysis. By the end of the research no larger collection of phages from a single environmental source investigated by these means had been found. The finding was confirmed by whole genome analysis of 7 selected bacteriophages. Moreover, particular bacteriophages revealed unusual biological features, like the ability to form plaques at low temperature (4 °C), resist high temperature (62 °C or 95 °C) or survive in the presence of an organic solvents (ethanol, acetone, DMSO, chloroform) or detergent (SDS, CTAB, sarkosyl) making them potentially interesting in the context of biotechnological applications.

Pseudorandomness of gene expression: a new evo-devo theory of ageing.

In contrast to the first part of life (development), ageing appears to be under less strict genetic control. The precise timing of events so characteristic of development seems to loosen its grasp, while stochastic and environmental factors seem to become the dominant force. Evolutionary theories put forward a decreasing evolutionary pressure over the course of life as the reason behind this pattern, yet dissenting views on ageing as a genetically programmed process linger. In this paper we address this dissent by presenting insights from an artificial evolutionary-developmental system, ET, and propose a new evo-devo theory of ageing-a theory that sees ageing as a continuation of development in the postreproductive period. In this theory both development and ageing are under genetic control. Nonetheless, while gene expression patterns that drive development are optimised by evolution, patterns that drive ageing are not optimised, because evolutionary pressure decreases with age. For these reasons, during ageing the changes orchestrated by genes are "pseudorandom"- deterministic but erratic-and their effects on an individual's health are more likely to be detrimental than beneficial. As such, they contribute to the continuous deterioration of bodily functions that characterise ageing.

Embryogenesis, morphogens and cancer stem cells: putting the puzzle together.

This paper describes a model which puts together three key elements of cancer theory: the analogies between embryogenesis and carcinogenesis, the role played in both processes by morphogens and related pathways, and the recently emerged paradigm of cancer stem cells. The model is called Epigenetic Tracking. Originally conceived as a model of embryonic development, it was later extended to interpret other aspects of biology, such as the presence of junk DNA, the phenomenon of ageing and the process of cancer formation. In this work we deepen our vision of carcinogenesis, and propose a novel hypothesis on the role of morphogen-processing pathways. According to the hypothesis, the interplay of these pathways leads in stem cells to the production of new transcription factors, which act as drivers of cellular differentiation. The disruption of these pathways, caused by mutations in specific genes, would represent the first and most distinctive event in the carcinogenic process. Our hypothesis allows us to make testable predictions on patterns of gene mutations involved in carcinogenesis. Our hypothesis also suggests that cancer stem cells can stay dormant until they are activated in a process that resembles activation of stem cells during tissue repair or at a specific time during development.

Low virus to prokaryote ratios in the cold: benthic viruses and prokaryotes in a subpolar marine ecosystem (Hornsund, Svalbard).

The density and spatial distribution of benthic viruses and prokaryotes in relation to biotic and abiotic factors were investigated in sediment cores collected in Hornsund, a permanently cold fjord on the West coast of Svalbard, Norway. The cores were obtained from the mouth of the fjord to the central basin, along a longitudinal transect. The results of our analyses showed lower densities of viruses (0.2 x 10(8) to 5.4 x 10(8) virus-like particles/g) and lower virus-to-prokaryote ratios (0.2-0.6, with the exception of the uppermost layer in the central basin, where the ratio was about 1.2) at the study site than generally found in the temperate areas, despite the relatively high organic matter content in subpolar sediments. Variations in benthic viral and prokaryote abundances along gradients of particle sedimentation rates, phytopigment concentrations, and macrobenthic species composition together suggested the influence of particle sedimentation and macrobenthic bioturbation on the abundance and spatial distribution ofprokaryotes and viruses in cold habitats.

Evolution of robustness to damage in artificial 3-dimensional development.

GReaNs is an Artificial Life platform we have built to investigate the general principles that guide evolution of multicellular development and evolution of artificial gene regulatory networks. The embryos develop in GReaNs in a continuous 3-dimensional (3D) space with simple physics. The developmental trajectories are indirectly encoded in linear genomes. The genomes are not limited in size and determine the topology of gene regulatory networks that are not limited in the number of nodes. The expression of the genes is continuous and can be modified by adding environmental noise. In this paper we evolved development of structures with a specific shape (an ellipsoid) and asymmetrical pattering (a 3D pattern inspired by the French flag problem), and investigated emergence of the robustness to damage in development and the emergence of the robustness to noise. Our results indicate that both types of robustness are related, and that including noise during evolution promotes higher robustness to damage. Interestingly, we have observed that some evolved gene regulatory networks rely on noise for proper behaviour.

A model of evolution of development based on germline penetration of new "no-junk" DNA.

There is a mounting body of evidence that somatic transposition may be involved in normal development of multicellular organisms and in pathology, especially cancer. Epigenetic Tracking (ET) is an abstract model of multicellular development, able to generate complex 3-dimensional structures. Its aim is not to model the development of a particular organism nor to merely summarise mainstream knowledge on genetic regulation of development. Rather, the goal of ET is to provide a theoretical framework to test new postulated genetic mechanisms, not fully established yet in mainstream biology. The first proposal is that development is orchestrated through a subset of cells which we call driver cells. In these cells, the cellular state determines a specific pattern of gene activation which leads to the occurrence of developmental events. The second proposal is that evolution of development is affected by somatic transposition events. We postulate that when the genome of a driver cell does not specify what developmental event should be undertaken when the cell is in a particular cellular state, somatic transposition events can reshape the genome, build new regulatory regions, and lead to a new pattern of gene activation in the cell. Our third hypothesis, not supported yet by direct evidence, but consistent with some experimental observations, is that these new "no-junk" sequences-regulatory regions created by transposable elements at new positions in the genome-can exit the cell and enter the germline, to be incorporated in the genome of the progeny. We call this mechanism germline penetration. This process allows heritable incorporation of novel developmental events in the developmental trajectory. In this paper we will present the model and link these three postulated mechanisms to biological observations.

Diversity of tailed phages in Baltic Sea sediment: large number of siphoviruses with extremely long tails.

We present the first attempt at quantitative analysis of morphological diversity of tailed viruses obtained from marine sediments without ultracentrifugation or enrichment on specific host strains. Sandy mud samples were collected in the Gulf of Gdansk in the spring, autumn and winter. VLPs were analyzed by transmission electron microscopy. The distribution of three groups of tailed phages was similar in all seasons (Siphoviridae: 52% on average; Myoviridae: 42%; Podoviridae: 6%). 19% of siphoviruses had prolate heads. Interestingly, 11% of siphoviral particles had tails longer than 300 nm, and 6% longer than 600 nm.

Metal and antibiotic resistance of bacteria isolated from the Baltic Sea.

The resistance of 49 strains of bacteria isolated from surface Baltic Sea waters to 11 antibiotics was analyzed and the resistance of selected strains to three metal ions (Ni2+, Mn2+, Zn2+) was tested. Most isolates belonged to Gammaproteobacteria (78%), while Alphaproteobacteria (8%), Actinobacteria (10%), and Bacteroidetes (4%) were less abundant. Even though previous reports suggested relationships between resistance and the presence of plasmids or the ability to produce pigments, no compelling evidence for such relationships was obtained for the strains isolated in this work. In particular, strains resistant to multiple antibiotics did not carry plasmids more frequently than sensitive strains. A relation between resistance and the four aminoglycosides tested (gentamycin, kanamycin, neomycin, and streptomycin), but not to spectinomycin, was demonstrated. This observation is of interest given that spectinomycin is not always classified as an aminoglycoside because it lacks a traditional sugar moiety. Statistical analysis indicated relationships between resistance to some antibiotics (ampicillin and erythromycin, chloramphenicol and erythromycin, chloramphenicol and tetracycline, erythromycin and tetracycline), suggesting the linkage of resistance genes for antibiotics belonging to different classes. The effects of NiSO4, ZnCl2 and MnCl2 on various media suggested that the composition of Marine Broth might result in low concentrations of Mn2+ due to chemical interactions that potentially lead to precipitation.

Plasmids derived from Gifsy-1/Gifsy-2, lambdoid prophages contributing to the virulence of Salmonella enterica serovar Typhimurium: implications for the evolution of replication initiation proteins of lambdoid phages and enterobacteria.

Gifsy-1 and Gifsy-2 are lambdoid prophages which contribute to the virulence of Salmonella enterica serovar Typhimurium. The nucleotide sequence of the replication region of both prophages is identical, and similar in organization to the replication region of bacteriophage lambda. To investigate the replication of the Gifsy phages and the relationship between Gifsy and host chromosome replication, a plasmid which contained all the genes and regulatory sequences required for autonomous replication in bacterial cells was constructed. This plasmid, pGifsy, was stably maintained in Escherichia coli cells. The helicase loader of the Gifsy phages is very similar to the DnaC protein of the host, a feature characteristic of a large group of prophages common in the sequenced genomes of pathogenic enterobacteria. This DnaC-like protein showed no similarity to the helicase loader of bacteriophage lambda and closely related phages. Interestingly, unlike plasmids derived from bacteriophage lambda (lambda plasmids), pGifsy did not require a gene encoding the putative helicase loader for replication, although deletion of this gene resulted in a decrease in plasmid copy number. Under these conditions, it was shown that the plasmid utilized the helicase loader coded by the host. On the other hand, the viral protein could not substitute for DnaC in bacterial chromosome replication. The results of the current study support the hypothesis that the enterobacterial helicase loader is of viral origin. This hypothesis explains why the gene for DnaC, the protein central to both replication initiation and replication restart in E. coli, is present in the genomes of Escherichia, Shigella, Salmonella and Buchnera, but not in the genomes of related enterobacteria.

Analysis of the overdispersed clock in the short-term evolution of hepatitis C virus: Using the E1/E2 gene sequences to infer infection dates in a single source outbreak.

The assumption of a molecular clock for dating events from sequence information is often frustrated by the presence of heterogeneity among evolutionary rates due, among other factors, to positively selected sites. In this work, our goal is to explore methods to estimate infection dates from sequence analysis. One such method, based on site stripping for clock detection, was proposed to unravel the clocklike molecular evolution in sequences showing high variability of evolutionary rates and in the presence of positive selection. Other alternatives imply accommodating heterogeneity in evolutionary rates at various levels, without eliminating any information from the data. Here we present the analysis of a data set of hepatitis C virus (HCV) sequences from 24 patients infected by a single individual with known dates of infection. We first used a simple criterion of relative substitution rate for site removal prior to a regression analysis. Time was regressed on maximum likelihood pairwise evolutionary distances between the sequences sampled from the source individual and infected patients. We show that it is indeed the fastest evolving sites that disturb the molecular clock and that these sites correspond to positively selected codons. The high computational efficiency of the regression analysis allowed us to compare the site-stripping scheme with random removal of sites. We demonstrate that removing the fast-evolving sites significantly increases the accuracy of estimation of infection times based on a single substitution rate. However, the time-of-infection estimations improved substantially when a more sophisticated and computationally demanding Bayesian method was used. This method was used with the same data set but keeping all the sequence positions in the analysis. Consequently, despite the distortion introduced by positive selection on evolutionary rates, it is possible to obtain quite accurate estimates of infection dates, a result of especial relevance for molecular epidemiology studies.

Weighted least-squares likelihood ratio test for branch testing in phylogenies reconstructed from distance measures.

A variety of analytical methods is available for branch testing in distance-based phylogenies. However, these methods are rarely used, possibly because the estimation of some of their statistics, especially the covariances, is not always feasible. We show that these difficulties can be overcome if some simplifying assumptions are made, namely distance independence. The weighted least-squares likelihood ratio test (WLS-LRT) we propose is easy to perform, using only the distances and some of their associated variances. If no variances are known, the use of the Felsenstein F-test, also based on weighted least squares, is discussed. Using simulated data and a data set of 43 mammalian mitochondrial sequences we demonstrate that the WLS-LRT performs as well as the generalized least-squares test, and indeed better for a large number of taxa data set. We thus show that the assumption of independence does not negatively affect the reliability or the accuracy of the least-squares approach. The results of the WLS-LRT are no worse than the results of the bootstrap methods, such as the Felsenstein bootstrap selection probability test and the Dopazo test. We also show that WLS-LRT can be applied in instances where other analytical methods are inappropriate. This point is illustrated by analyzing the relationships between human immunodeficiency virus type 1 (HIV-1) sequences isolated from various organs of different individuals.

Evolution of lambdoid replication modules.

Comparison of the putative iteron-binding proteins of lambdoid phages allows us to propose that in the case of lambdoid replication modules, the units on which natural selection acts do not coincide with the open reading frames. Rather, the first replication gene is split into two segments, and its 3' part (corresponding to the C-terminal domain of the iteron-binding protein) forms one unit with the second gene. We also propose from the phylogenetic analysis of phage-encoded homologs of E. coli DnaB and DnaC, that the recombination with the host sequences is not frequent. Accessory ATP-ases for helicase loading (E. coli DnaC homologs) may not be universal replication proteins. Our analysis may suggest that the bacterial helicase loaders might be of phage origin. The comparison of DnaC homologs of enterobacteria and enterobacterial phages supports the experimental data on residues important in interaction with DnaB. We propose that construction of plasmids carrying the replication origins of lambdoid prophages could be useful not only in further research on DNA replication but also on the role of these prophages in shuttling genes for bacterial virulence. The phage replication sequences could be also useful for identification of clinical enterobacterial isolates.

The optimal eukaryotic signal for translation initiation from non-AUG codons, present upstream of bacteriophage lambda P cistron, is inactive in Escherichia coli.

Expression of the replication genes of bacteriophage lambda, O and P, is believed to be translationally coupled. However, it was previously noted that, under conditions of amino acid starvation, when O is not synthesized, P continues to be expressed at a relatively high level. The results presented in this report, contrary to the previously presented hypothesis, suggest that an AGACUGGAU sequence (an optimal context for translation initiation from non-AUG codons in eukaryotes, and present upstream the P cistron) is inactive in Escherichia coli. Comparative sequence analysis confirms that such a signal is unlikely to be important for P synthesis. Instead, a weak Shine-Dalgarno sequence may be present upstream the P cistron, and be active in the absence of O gene expression.

Vibrio harveyi bioluminescence plays a role in stimulation of DNA repair.

Although the genetics and biochemistry of bacterial luminescence have been investigated extensively, the biological role of this phenomenon remains unclear. Here it is shown that luxA, luxB and luxD mutants (unable to emit light) of the marine bacterium Vibrio harveyi are significantly more sensitive to UV irradiation when cultivated in the dark after irradiation than when cultivated under a white fluorescent lamp. This difference was much less pronounced in the wild-type (luminescent) V. harveyi strain. Survival of UV-irradiated Escherichia coli wild-type cells depended on subsequent cultivation conditions (in the dark or in the presence of external light). However, after UV irradiation, the percentage of surviving E. coli cells that bear V. harveyi genes responsible for luminescence was significantly higher than that of non-luminescent E. coli, irrespective of the subsequent cultivation conditions. Moreover, it is demonstrated that luminescence of V. harveyi can be stimulated by UV irradiation even in diluted cultures, under conditions when light emission by these bacteria is normally impaired due to quorum sensing regulation. It is proposed that luminescent bacteria have an internal source of light which could be used in DNA repair by a photoreactivation process. Therefore, production of internal light ensuring effective DNA repair seems to be at least one of the biological functions of bacterial luminescence.