The Landscape of Realized Homologous Recombination in Pathogenic Bacteria

Recombination enhances the adaptive potential of organisms by allowing genetic variants to be tested on multiple genomic backgrounds. Its distribution in the genome can provide insight into the evolutionary forces that underlie traits, such as the emergence of pathogenicity. Here, we examined landscapes of realized homologous recombination of 500 genomes from ten bacterial species and found all species have “hot” regions with elevated rates relative to the genome average. We examined the size, gene content, and chromosomal features associated with these regions and the correlations between closely related species. The recombination landscape is variable and evolves rapidly. For example in Salmonella, only short regions of around 1 kb in length are hot whereas in the closely related species Escherichia coli, some hot regions exceed 100 kb, spanning many genes. Only Streptococcus pyogenes shows evidence for the positive correlation between GC content and recombination that has been reported for several eukaryotes. Genes with function related to the cell surface/membrane are often found in recombination hot regions but E. coli is the only species where genes annotated as “virulence associated” are consistently hotter. There is also evidence that some genes with “housekeeping” functions tend to be overrepresented in cold regions. For example, ribosomal proteins showed low recombination in all of the species. Among specific genes, transferrin-binding proteins are recombination hot in all three of the species in which they were found, and are subject to interspecies recombination

Evolution of DNA Double-Strand Break Repair by Gene Conversion: Coevolution Between a Phage and Restriction-Modification System

The necessity to repair genome damage has been considered to be an immediate factor responsible for the origin of sex. Indeed, attack by a cellular restriction enzyme of invading DNA from several bacteriophages initiates recombinational repair by gene conversion if there is homologous DNA. In the present work, we modeled the interaction between a bacteriophage and a bacterium carrying a restriction enzyme as antagonistic coevolution. We assume a locus on the bacteriophage genome has either a restriction-sensitive or a -resistant allele, and another locus determines whether it is recombination/repair-proficient or -defective. A restriction break can be repaired by a co-infecting phage genome if one of them is recombination/repair-proficient. We define the fitness of phage (resistant/sensitive and repair-positive/-negative) genotypes and bacterial (restriction-positive/-negative) genotypes by assuming random encounter of the genotypes, with a given probabilities of single and double infections, and the costs of resistance, repair and restriction. Our results show the evolution of the repair allele depends on b1 / b0 , the ratio of the burst size b1 under damage to host cell physiology induced by an unrepaired double-strand break to the default burst size b0 . It was not until this effect was taken into account that the evolutionary advantage of DNA repair became apparent

Genome-wide survey of codons under diversifying selection in a highly recombining bacterial species, Helicobacter pylori

Selection has been a central issue in biology in eukaryotes as well as prokaryotes. Inference of selection in recombining bacterial species, compared with clonal ones, has been a challenge. It is not known how codons under diversifying selection are distributed along the chromosome or among functional categories or how frequently such codons are subject to mutual homologous recombination. Here, we explored these questions by analysing genes present in >90% among 29 genomes of Helicobacter pylori, one of the bacterial species with the highest mutation and recombination rates. By a method for recombining sequences, we identified codons under diversifying selection (dN/dS > 1), which were widely distributed and accounted for ∼0.2% of all the codons of the genome. The codons were enriched in genes of host interaction/cell surface and genome maintenance (DNA replication, recombination, repair, and restriction modification system). The encoded amino acid residues were sometimes found adjacent to critical catalytic/binding residues in protein structures. Furthermore, by estimating the intensity of homologous recombination at a single nucleotide level, we found that these codons appear to be more frequently subject to recombination. We expect that the present study provides a new approach to population genomics of selection in recombining prokaryotes

Genomic structure and insertion sites of Helicobacter pylori prophages from various geographical origins

We present the full genomic sequences, insertion sites and phylogenetic analysis of 28 prophages found in H. pylori isolates from patients of distinct disease types, ranging from gastritis to gastric cancer, and geographic origins, covering most continents. The gentic diversity of H pylori is known to be influenced by these genomic elements including prophages who’s geneomes range from 22.6 to 33.0 Kbp. There was a high conservation of integration site shared in over 50% of cases with greater than 40% or prophage genomes harbouring insertion sequences (IS). Furthermore prophage genomes present a robust phylogeographic pattern, revealing four distinct clusters: one African, one Asian and two European prophage populations. There was evidence of recombination within the genome of some prophages, which resulted in genome mosaics composed by different populations, which may yield additional H. pylori phenotypes

Identification of biomolecule mass transport and binding rate parameters in living cells by inverse modeling

BACKGROUND: Quantification of in-vivo biomolecule mass transport and reaction rate parameters from experimental data obtained by Fluorescence Recovery after Photobleaching (FRAP) is becoming more important. METHODS AND RESULTS: The Osborne-Moré extended version of the Levenberg-Marquardt optimization algorithm was coupled with the experimental data obtained by the Fluorescence Recovery after Photobleaching (FRAP) protocol, and the numerical solution of a set of two partial differential equations governing macromolecule mass transport and reaction in living cells, to inversely estimate optimized values of the molecular diffusion coefficient and binding rate parameters of GFP-tagged glucocorticoid receptor. The results indicate that the FRAP protocol provides enough information to estimate one parameter uniquely using a nonlinear optimization technique. Coupling FRAP experimental data with the inverse modeling strategy, one can also uniquely estimate the individual values of the binding rate coefficients if the molecular diffusion coefficient is known. One can also simultaneously estimate the dissociation rate parameter and molecular diffusion coefficient given the pseudo-association rate parameter is known. However, the protocol provides insufficient information for unique simultaneous estimation of three parameters (diffusion coefficient and binding rate parameters) owing to the high intercorrelation between the molecular diffusion coefficient and pseudo-association rate parameter. Attempts to estimate macromolecule mass transport and binding rate parameters simultaneously from FRAP data result in misleading conclusions regarding concentrations of free macromolecule and bound complex inside the cell, average binding time per vacant site, average time for diffusion of macromolecules from one site to the next, and slow or rapid mobility of biomolecules in cells. CONCLUSION: To obtain unique values for molecular diffusion coefficient and binding rate parameters from FRAP data, we propose conducting two FRAP experiments on the same class of macromolecule and cell. One experiment should be used to measure the molecular diffusion coefficient independently of binding in an effective diffusion regime and the other should be conducted in a reaction dominant or reaction-diffusion regime to quantify binding rate parameters. The method described in this paper is likely to be widely used to estimate in-vivo biomolecule mass transport and binding rate parameters

Environmental Stress-Dependent Effects of Deletions Encompassing Hsp70Ba on Canalization and Quantitative Trait Asymmetry in Drosophila melanogaster

Hsp70 genes may influence the expression of wing abnormalities in Drosophila melanogaster but their effects on variability in quantitative characters and developmental instability are unclear. In this study, we focused on one of the six Hsp70 genes, Hsp70Ba, and investigated its effects on within-and among-individual variability in orbital bristle number, sternopleural bristle number, wing size and wing shape under different environmental conditions. To do this, we studied a newly constructed deletion, Df(3R)ED5579, which encompasses Hsp70Ba and nine non-Hsp genes, in the heterozygous condition and another, Hsp70Ba304, which deletes only Hsp70Ba, in the homozygous condition. We found no significant effect of both deletions on within-individual variation quantified by fluctuating asymmetry (FA) of morphological traits. On the other hand, the Hsp70Ba304/Hsp70Ba304 genotype significantly increased among-individual variation quantified by coefficient of variation (CV) of bristle number and wing size in female, while the Df(3R)ED5579 heterozygote showed no significant effect. The expression level of Hsp70Ba in the deletion heterozygote was 6 to 20 times higher than in control homozygotes, suggesting that the overexpression of Hsp70Ba did not influence developmental stability or canalization significantly. These findings suggest that the absence of expression of Hsp70Ba increases CV of some morphological traits and that HSP70Ba may buffer against environmental perturbations on some quantitative traits