Identification of DNA motifs implicated in maintenance of bacterial core genomes by predictive modeling

RecBCD enzyme and Chi recombination hotspots as determinants of self vs. non-self: Myths and mechanisms

Bacteria face a challenge when DNA enters their cells by transformation, mating, or phage infection. Should they treat this DNA as an invasive foreigner and destroy it, or consider it one of their own and potentially benefit from incorporating new genes or alleles to gain useful functions? It is frequently stated that the short nucleotide sequence Chi (5' GCTGGTGG 3'), a hotspot of homologous genetic recombination recognized by Escherichia coli's RecBCD helicase-nuclease, allows E. coli to distinguish its DNA (self) from any other DNA (non-self) and to destroy non-self DNA, and that Chi is "over-represented" in the E. coli genome. We show here that these latter statements (dogmas) are not supported by available evidence. We note Chi's wide-spread occurrence and activity in distantly related bacterial species and phages. We illustrate multiple, highly non-random features of the genomes of E. coli and coliphage P1 that account for Chi's high frequency and genomic position, leading us to propose that P1 selects for Chi's enhancement of recombination, whereas E. coli selects for the preferred codons in Chi. We discuss other, previously described mechanisms for self vs. non-self determination involving RecBCD and for RecBCD's destruction of DNA that cannot recombine, whether foreign or domestic, with or without Chi.

Cleavage of viral DNA by restriction endonucleases stimulates the type II CRISPR-Cas immune response

Prokaryotic organisms have developed multiple defense systems against phages; however, little is known about whether and how these interact with each other. Here, we studied the connection between two of the most prominent prokaryotic immune systems: restriction-modification and CRISPR. While both systems employ enzymes that cleave a specific DNA sequence of the invader, CRISPR nucleases are programmed with phage-derived spacer sequences, which are integrated into the CRISPR locus upon infection. We found that restriction endonucleases provide a short-term defense, which is rapidly overcome through methylation of the phage genome. In a small fraction of the cells, however, restriction results in the acquisition of spacer sequences from the cleavage site, which mediates a robust type II-A CRISPR-Cas immune response against the methylated phage. This mechanism is reminiscent of eukaryotic immunity in which the innate response offers a first temporary line of defense and also activates a second and more robust adaptive response.

DNA Repair in Staphylococcus aureus

Staphylococcus aureus is a common cause of both superficial and invasive infections of humans and animals. Despite a potent host response and apparently appropriate antibiotic therapy, staphylococcal infections frequently become chronic or recurrent, demonstrating a remarkable ability of S. aureus to withstand the hostile host environment. There is growing evidence that staphylococcal DNA repair makes important contributions to the survival of the pathogen in host tissues, as well as promoting the emergence of mutants that resist host defenses and antibiotics. While much of what we know about DNA repair in S. aureus is inferred from studies with model organisms, the roles of specific repair mechanisms in infection are becoming clear and differences with Bacillus subtilis and Escherichia coli have been identified. Furthermore, there is growing interest in staphylococcal DNA repair as a target for novel therapeutics that sensitize the pathogen to host defenses and antibiotics. In this review, we discuss what is known about staphylococcal DNA repair and its role in infection, examine how repair in S. aureus is similar to, or differs from, repair in well-characterized model organisms, and assess the potential of staphylococcal DNA repair as a novel therapeutic target.

Evolution of Chi motifs in Proteobacteria

Homologous recombination is a key pathway found in nearly all bacterial taxa. The recombination complex not only allows bacteria to repair DNA double-strand breaks but also promotes adaption through the exchange of DNA between cells. In Proteobacteria, this process is mediated by the RecBCD complex, which relies on the recognition of a DNA motif named Chi to initiate recombination. The Chi motif has been characterized in Escherichia coli and analogous sequences have been found in several other species from diverse families, suggesting that this mode of action is widespread across bacteria. However, the sequences of Chi-like motifs are known for only five bacterial species: E. coli, Haemophilus influenzae, Bacillus subtilis, Lactococcus lactis, and Staphylococcus aureus. In this study, we detected putative Chi motifs in a large dataset of Proteobacteria and identified four additional motifs sharing high sequence similarity and similar properties to the Chi motif of E. coli in 85 species of Proteobacteria. Most Chi motifs were detected in Enterobacteriaceae and this motif appears well conserved in this family. However, we did not detect Chi motifs for the majority of Proteobacteria, suggesting that different motifs are used in these species. Altogether these results substantially expand our knowledge on the evolution of Chi motifs and on the recombination process in bacteria.

Ribosomal protein L4 of Lactobacillus rhamnosus LRB alters resistance to macrolides and other antibiotics

Lactobacillus rhamnosus is an important lactic acid bacterium that is predominantly used as a probiotic supplement. This bacterium secretes immunomodulatory and antibacterial peptides that are necessary for the probiotic trait. This organism also occupies diverse ecological niches, such as gastrointestinal tracts and the oral cavity. Several studies have shown that L. rhamnosus is prone to spontaneous genome rearrangement irrespective of the ecological origins. We previously characterized an oral isolate of L. rhamnosus, LRB, which is genetically closely related to the widely used probiotic strain L. rhamnosus LGG. In this study, we isolated a nontargeted mutant that was particularly sensitive to acid stress. Using next generation sequencing, we further mapped the putative mutations in the genome and found that the mutant had acquired a deletion of 75 base pairs in the rplD gene that encodes the large ribosomal subunit L4. The mutant had a growth defect at 37°C and at ambient temperature. Further antibiotic sensitivity analyses indicated that the mutant is relatively more resistant to erythromycin and chloramphenicol; two antibiotics that target the 50S subunit. In contrast, the mutant was more sensitive to tetracycline, which targets the 30S subunit. Thus, it appears that nontargeted mutations could significantly alter the antibiotic resistance profile of L. rhamnosus. Our study raises concern that probiotic use of L. rhamnosus should be carefully monitored to avoid unintended consequences.

A conformational switch in response to Chi converts RecBCD from phage destruction to DNA repair

The RecBCD complex plays key roles in phage DNA degradation, CRISPR array acquisition (adaptation) and host DNA repair. The switch between these roles is regulated by a DNA sequence called Chi. We report cryo-EM structures of the Escherichia coli RecBCD complex bound to several different DNA forks containing a Chi sequence, including one in which Chi is recognised and others in which it is not. The Chi-recognised structure shows conformational changes in regions of the protein that contact Chi and reveals a tortuous path taken by the DNA. Sequence specificity arises from interactions with both the RecC subunit and the sequence itself. These structures provide molecular details for how Chi is recognised and insights into the changes that occur in response to Chi binding that switch RecBCD from bacteriophage destruction and CRISPR spacer acquisition, to constructive host DNA repair.

Genetics of Lactococci

Lactococcus lactis is the best characterized species among the lactococci, and among the most consumed food-fermenting bacteria worldwide. Thanks to their importance in industrialized food production, lactococci are among the lead bacteria understood for fundamental metabolic pathways that dictate growth and survival properties. Interestingly, lactococci belong to the Streptococcaceae family, which includes food, commensal and virulent species. As basic metabolic pathways (e.g., respiration, metal homeostasis, nucleotide metabolism) are now understood to underlie virulence, processes elucidated in lactococci could be important for understanding pathogen fitness and synergy between bacteria. This chapter highlights major findings in lactococci and related bacteria, and covers five themes: distinguishing features of lactococci, metabolic capacities including the less known respiration metabolism in Streptococcaceae, factors and pathways modulating stress response and fitness, interbacterial dialogue via metabolites, and novel applications in health and biotechnology.

Genome Variation in the Model Halophilic Bacterium Salinibacter ruber

The halophilic bacterium Salinibacter ruber is an abundant and ecologically important member of halophilic communities worldwide. Given its broad distribution and high intraspecific genetic diversity, S. ruber is considered one of the main models for ecological and evolutionary studies of bacterial adaptation to hypersaline environments. However, current insights on the genomic diversity of this species is limited to the comparison of the genomes of two co-isolated strains. Here, we present a comparative genomic analysis of eight S. ruber strains isolated at two different time points in each of two different Mediterranean solar salterns. Our results show an open pangenome with contrasting evolutionary patterns in the core and accessory genomes. We found that the core genome is shaped by extensive homologous recombination (HR), which results in limited sequence variation within population clusters. In contrast, the accessory genome is modulated by horizontal gene transfer (HGT), with genomic islands and plasmids acting as gateways to the rest of the genome. In addition, both types of genetic exchange are modulated by restriction and modification (RM) or CRISPR-Cas systems. Finally, genes differentially impacted by such processes reveal functional processes potentially relevant for environmental interactions and adaptation to extremophilic conditions. Altogether, our results support scenarios that conciliate “Neutral” and “Constant Diversity” models of bacterial evolution.

Biochemical characterization of RecBCD enzyme from an Antarctic Pseudomonas species and identification of its cognate Chi (χ) sequence

Pseudomonas syringae Lz4W RecBCD enzyme, RecBCD^Ps, is a trimeric protein complex comprised of RecC, RecB, and RecD subunits. RecBCD enzyme is essential for P. syringae growth at low temperature, and it protects cells from low temperature induced replication arrest. In this study, we show that the RecBCD^Ps enzyme displays distinct biochemical behaviors. Unlike E. coli RecBCD enzyme, the RecD subunit is indispensable for RecBCD^Ps function. The RecD motor activity is essential for the Chi-like fragments production in P. syringae, highlighting a distinct role for P. syringae RecD subunit in DNA repair and recombination process. Here, we demonstrate that the RecBCD^Ps enzyme recognizes a unique octameric DNA sequence, 5′-GCTGGCGC-3′ (Chi^Ps) that attenuates nuclease activity of the enzyme when it enters dsDNA from the 3′-end. We propose that the reduced translocation activities manifested by motor-defective mutants cause cold sensitivity in P. syrinage; emphasizing the importance of DNA processing and recombination functions in rescuing low temperature induced replication fork arrest.

CRISPR-Cas systems exploit viral DNA injection to establish and maintain adaptive immunity

CRISPR-Cas systems provide protection against viral¹ and plasmid² infection by capturing short DNA sequences from these invaders and integrating them into the CRISPR locus of the prokaryotic host¹. These sequences, known as spacers, are transcribed into short RNA guides^3–5 that specify the cleavage site of Cas nucleases in the genome of the invader^6–8. When spacer sequences are acquired during viral infection is not known. To investigate this, we followed spacer acquisition in Staphylococcus aureus cells harboring a type II CRISPR-Cas9 system after infection with the staphylococcal bacteriophage ϕ12. We found that new spacers are acquired immediately following infection preferentially from the cos site, the viral free DNA end that is first injected into the cell. Analysis of spacer acquisition after infection with mutant phages demonstrated that most spacers are acquired during DNA injection, but not during other stages of the viral cycle that produce free DNA ends, such as DNA replication or packaging. Finally, we showed that spacers acquired from early-injected genomic regions, which direct Cas9 cleavage of the viral DNA immediately after infection, provide better immunity than spacers acquired from late-injected regions. Our results reveal that CRISPR-Cas systems exploit the phage life cycle to generate a pattern of spacer acquisition that ensures the success of the CRISPR immune response.

Unexpected DNA context-dependence identifies a new determinant of Chi recombination hotspots

Homologous recombination occurs especially frequently near special chromosomal sites called hotspots. In Escherichia coli, Chi hotspots control RecBCD enzyme, a protein machine essential for the major pathway of DNA break-repair and recombination. RecBCD generates recombinogenic single-stranded DNA ends by unwinding DNA and cutting it a few nucleotides to the 3′ side of 5′ GCTGGTGG 3′, the sequence historically equated with Chi. To test if sequence context affects Chi activity, we deep-sequenced the products of a DNA library containing 10 random base-pairs on each side of the Chi sequence and cut by purified RecBCD. We found strongly enhanced cutting at Chi with certain preferred sequences, such as A or G at nucleotides 4–7, on the 3′ flank of the Chi octamer. These sequences also strongly increased Chi hotspot activity in E. coli cells. Our combined enzymatic and genetic results redefine the Chi hotspot sequence, implicate the nuclease domain in Chi recognition, indicate that nicking of one strand at Chi is RecBCD's biologically important reaction in living cells, and enable more precise analysis of Chi's role in recombination and genome evolution.

Ebolavirus Database: Gene and Protein Information Resource for Ebolaviruses

Ebola Virus Disease (EVD) is a life-threatening haemorrhagic fever in humans. Even though there are many reports on EVD, the protein precursor functions and virulent factors of ebolaviruses remain poorly understood. Comparative analyses of Ebolavirus genomes will help in the identification of these important features. This prompted us to develop the Ebolavirus Database (EDB) and we have provided links to various tools that will aid researchers to locate important regions in both the genomes and proteomes of Ebolavirus. The genomic analyses of ebolaviruses will provide important clues for locating the essential and core functional genes. The aim of EDB is to act as an integrated resource for ebolaviruses and we strongly believe that the database will be a useful tool for clinicians, microbiologists, health care workers, and bioscience researchers.

Coevolution of the Organization and Structure of Prokaryotic Genomes

The cytoplasm of prokaryotes contains many molecular machines interacting directly with the chromosome. These vital interactions depend on the chromosome structure, as a molecule, and on the genome organization, as a unit of genetic information. Strong selection for the organization of the genetic elements implicated in these interactions drives replicon ploidy, gene distribution, operon conservation, and the formation of replication-associated traits. The genomes of prokaryotes are also very plastic with high rates of horizontal gene transfer and gene loss. The evolutionary conflicts between plasticity and organization lead to the formation of regions with high genetic diversity whose impact on chromosome structure is poorly understood. Prokaryotic genomes are remarkable documents of natural history because they carry the imprint of all of these selective and mutational forces. Their study allows a better understanding of molecular mechanisms, their impact on microbial evolution, and how they can be tinkered in synthetic biology.

YOC, A new strategy for pairwise alignment of collinear genomes

Background

Comparing and aligning genomes is a key step in analyzing closely related genomes. Despite the development of many genome aligners in the last 15 years, the problem is not yet fully resolved, even when aligning closely related bacterial genomes of the same species. In addition, no procedures are available to assess the quality of genome alignments or to compare genome aligners.

Results

We designed an original method for pairwise genome alignment, named YOC, which employs a highly sensitive similarity detection method together with a recent collinear chaining strategy that allows overlaps. YOC improves the reliability of collinear genome alignments, while preserving or even improving sensitivity. We also propose an original qualitative evaluation criterion for measuring the relevance of genome alignments. We used this criterion to compare and benchmark YOC with five recent genome aligners on large bacterial genome datasets, and showed it is suitable for identifying the specificities and the potential flaws of their underlying strategies.

Conclusions

The YOC prototype is available at https://github.com/ruricaru/YOC. It has several advantages over existing genome aligners: (1) it is based on a simplified two phase alignment strategy, (2) it is easy to parameterize, (3) it produces reliable genome alignments, which are easier to analyze and to use.

Electronic supplementary material

The online version of this article (doi:10.1186/s12859-015-0530-3) contains supplementary material, which is available to authorized users.

Haemophilus influenzae Genome Database (HIGDB): a single point web resource for Haemophilus influenzae

BACKGROUND

Haemophilus influenzae (H. Influenzae) is the causative agent of pneumonia, bacteraemia and meningitis. The organism is responsible for large number of deaths in both developed and developing countries. Even-though the first bacterial genome to be sequenced was that of H. Influenzae, there is no exclusive database dedicated for H. Influenzae. This prompted us to develop the Haemophilus influenzae Genome Database (HIGDB).

METHODS

All data of HIGDB are stored and managed in MySQL database. The HIGDB is hosted on Solaris server and developed using PERL modules. Ajax and JavaScript are used for the interface development.

RESULTS

The HIGDB contains detailed information on 42,741 proteins, 18,077 genes including 10 whole genome sequences and also 284 three dimensional structures of proteins of H. influenzae. In addition, the database provides "Motif search" and "GBrowse". The HIGDB is freely accessible through the URL: http://bioserver1.physics.iisc.ernet.in/HIGDB/.

DISCUSSION

The HIGDB will be a single point access for bacteriological, clinical, genomic and proteomic information of H. influenzae. The database can also be used to identify DNA motifs within H. influenzae genomes and to compare gene or protein sequences of a particular strain with other strains of H. influenzae.

Streptococcus pneumoniae Genome Database (SPGDB): a database for strain specific comparative analysis of Streptococcus pneumoniae genes and proteins

Streptococcus pneumoniae causes pneumonia, septicemia and meningitis. S. pneumoniae is responsible for significant mortality both in children and in the elderly. In recent years, the whole genome sequencing of various S. pneumoniae strains have increased manifold and there is an urgent need to provide organism specific annotations to the scientific community. This prompted us to develop the Streptococcus pneumoniae Genome Database (SPGDB) to integrate and analyze the completely sequenced and available S. pneumoniae genome sequences. Further, links to several tools are provided to compare the pool of gene and protein sequences, and proteins structure across different strains of S. pneumoniae. SPGDB aids in the analysis of phenotypic variations as well as to perform extensive genomics and evolutionary studies with reference to S. pneumoniae. The database will be updated at regular intervals and is freely accessible through the URL: http://pranag.physics.iisc.ernet.in/SPGDB/.

Single molecule approaches to monitor the recognition and resection of double-stranded DNA breaks during homologous recombination

The fate of a cell depends on its ability to repair the many double-stranded DNA breaks (DSBs) that occur during normal metabolism. Improper DSB repair may result in genomic instability, cancer, or other genetic diseases. The repair of a DSB can be initiated by the recognition and resection of a duplex DNA end to form a 3'-terminated single-stranded DNA overhang. This task is carried out by different single-strand exonucleases, endonucleases, and helicases that work in a coordinated manner. This manuscript reviews the different single-molecule approaches that have been employed to characterize the structural features of these molecular machines, as well as the intermediates and products formed during the process of DSB repair. Imaging techniques have unveiled the structural organization of complexes involved in the tethering and recognition of DSBs. In addition to that static picture, single molecule studies on the dynamics of helicase-nuclease complexes responsible for the processive resection of DSBs have provided detailed mechanistic insights into their function.

Manipulating or superseding host recombination functions: a dilemma that shapes phage evolvability

Phages, like many parasites, tend to have small genomes and may encode autonomous functions or manipulate those of their hosts'. Recombination functions are essential for phage replication and diversification. They are also nearly ubiquitous in bacteria. The E. coli genome encodes many copies of an octamer (Chi) motif that upon recognition by RecBCD favors repair of double strand breaks by homologous recombination. This might allow self from non-self discrimination because RecBCD degrades DNA lacking Chi. Bacteriophage Lambda, an E. coli parasite, lacks Chi motifs, but escapes degradation by inhibiting RecBCD and encoding its own autonomous recombination machinery. We found that only half of 275 lambdoid genomes encode recombinases, the remaining relying on the host's machinery. Unexpectedly, we found that some lambdoid phages contain extremely high numbers of Chi motifs concentrated between the phage origin of replication and the packaging site. This suggests a tight association between replication, packaging and RecBCD-mediated recombination in these phages. Indeed, phages lacking recombinases strongly over-represent Chi motifs. Conversely, phages encoding recombinases and inhibiting host recombination machinery select for the absence of Chi motifs. Host and phage recombinases use different mechanisms and the latter are more tolerant to sequence divergence. Accordingly, we show that phages encoding their own recombination machinery have more mosaic genomes resulting from recent recombination events and have more diverse gene repertoires, i.e. larger pan genomes. We discuss the costs and benefits of superseding or manipulating host recombination functions and how this decision shapes phage genome structure and evolvability.

Bacterial viruses, called bacteriophages, are extremely abundant in the biosphere. They have key roles in the regulation of bacterial populations and in the diversification of bacterial genomes. Among these viruses, lambdoid phages are very abundant in enterobacteria and exchange genetic material very frequently. This latter process is thought to increase phage diversity and therefore facilitate adaptation to hosts. Recombination is also essential for the replication of many lambdoid phages. Lambdoids have been described to encode their own recombination genes and inhibit their hosts'. In this study, we show that lambdoids are split regarding their capacity to encode autonomous recombination functions and that this affects the abundance of recombination-related sequence motifs. Half of the phages encode an autonomous system and inhibit their hosts'. The trade-off between superseding and manipulating the hosts' recombination functions has important consequences. The phages encoding autonomous recombination functions have more diverse gene repertoires and recombine more frequently. Viruses, as many other parasites, have small genomes and depend on their hosts for several housekeeping functions. Hence, they often face trade-offs between supersession and manipulation of molecular machineries. Our results suggest these trade-offs may shape viral gene repertoires, their sequence composition and even influence their evolvability.

On the mechanism of recombination hotspot scanning during double-stranded DNA break resection

Double-stranded DNA break repair by homologous recombination is initiated by resection of free DNA ends to produce a 3'-ssDNA overhang. In bacteria, this reaction is catalyzed by helicase-nuclease complexes such as AddAB in a manner regulated by specific recombination hotspot sequences called Crossover hotspot instigator (Chi). We have used magnetic tweezers to investigate the dynamics of AddAB translocation and hotspot scanning during double-stranded DNA break resection. AddAB was prone to stochastic pausing due to transient recognition of Chi-like sequences, unveiling an antagonistic relationship between DNA translocation and sequence-specific DNA recognition. Pauses at bona fide Chi sequences were longer, were nonexponentially distributed, and resulted in an altered velocity upon restart of translocation downstream of Chi. We propose a model for the recognition of Chi sequences to explain the origin of pausing during failed and successful hotspot recognition.

Google matrix analysis of DNA sequences

For DNA sequences of various species we construct the Google matrix [Formula: see text] of Markov transitions between nearby words composed of several letters. The statistical distribution of matrix elements of this matrix is shown to be described by a power law with the exponent being close to those of outgoing links in such scale-free networks as the World Wide Web (WWW). At the same time the sum of ingoing matrix elements is characterized by the exponent being significantly larger than those typical for WWW networks. This results in a slow algebraic decay of the PageRank probability determined by the distribution of ingoing elements. The spectrum of [Formula: see text] is characterized by a large gap leading to a rapid relaxation process on the DNA sequence networks. We introduce the PageRank proximity correlator between different species which determines their statistical similarity from the view point of Markov chains. The properties of other eigenstates of the Google matrix are also discussed. Our results establish scale-free features of DNA sequence networks showing their similarities and distinctions with the WWW and linguistic networks.

Dialects of the DNA uptake sequence in Neisseriaceae

In all sexual organisms, adaptations exist that secure the safe reassortment of homologous alleles and prevent the intrusion of potentially hazardous alien DNA. Some bacteria engage in a simple form of sex known as transformation. In the human pathogen Neisseria meningitidis and in related bacterial species, transformation by exogenous DNA is regulated by the presence of a specific DNA Uptake Sequence (DUS), which is present in thousands of copies in the respective genomes. DUS affects transformation by limiting DNA uptake and recombination in favour of homologous DNA. The specific mechanisms of DUS-dependent genetic transformation have remained elusive. Bioinformatic analyses of family Neisseriaceae genomes reveal eight distinct variants of DUS. These variants are here termed DUS dialects, and their effect on interspecies commutation is demonstrated. Each of the DUS dialects is remarkably conserved within each species and is distributed consistent with a robust Neisseriaceae phylogeny based on core genome sequences. The impact of individual single nucleotide transversions in DUS on meningococcal transformation and on DNA binding and uptake is analysed. The results show that a DUS core 5'-CTG-3' is required for transformation and that transversions in this core reduce DNA uptake more than two orders of magnitude although the level of DNA binding remains less affected. Distinct DUS dialects are efficient barriers to interspecies recombination in N. meningitidis, N. elongata, Kingella denitrificans, and Eikenella corrodens, despite the presence of the core sequence. The degree of similarity between the DUS dialect of the recipient species and the donor DNA directly correlates with the level of transformation and DNA binding and uptake. Finally, DUS-dependent transformation is documented in the genera Eikenella and Kingella for the first time. The results presented here advance our understanding of the function and evolution of DUS and genetic transformation in bacteria, and define the phylogenetic relationships within the Neisseriaceae family.

Diverse functions of restriction-modification systems in addition to cellular defense

Restriction-modification (R-M) systems are ubiquitous and are often considered primitive immune systems in bacteria. Their diversity and prevalence across the prokaryotic kingdom are an indication of their success as a defense mechanism against invading genomes. However, their cellular defense function does not adequately explain the basis for their immaculate specificity in sequence recognition and nonuniform distribution, ranging from none to too many, in diverse species. The present review deals with new developments which provide insights into the roles of these enzymes in other aspects of cellular function. In this review, emphasis is placed on novel hypotheses and various findings that have not yet been dealt with in a critical review. Emerging studies indicate their role in various cellular processes other than host defense, virulence, and even controlling the rate of evolution of the organism. We also discuss how R-M systems could have successfully evolved and be involved in additional cellular portfolios, thereby increasing the relative fitness of their hosts in the population.

Detecting correlations among functional-sequence motifs

Sequence motifs are words of nucleotides in DNA with biological functions, e.g., gene regulation. Identification of such words proceeds through rejection of Markov models on the expected motif frequency along the genome. Additional biological information can be extracted from the correlation structure among patterns of motif occurrences. In this paper a log-linear multivariate intensity Poisson model is estimated via expectation maximization on a set of motifs along the genome of E. coli K12. The proposed approach allows for excitatory as well as inhibitory interactions among motifs and between motifs and other genomic features like gene occurrences. Our findings confirm previous stylized facts about such types of interactions and shed new light on genome-maintenance functions of some particular motifs. We expect these methods to be applicable to a wider set of genomic features.

Molecular determinants responsible for recognition of the single-stranded DNA regulatory sequence, χ, by RecBCD enzyme

The RecBCD enzyme is important for both restriction of foreign DNA and recombinational DNA repair. Switching enzyme function from the destructive antiviral state to the productive recombinational state is regulated by the recombination hotspot, χ (5'-GCTGGTGG-3'). Recognition of χ is unique in that it is recognized as a specific sequence within single-stranded DNA (ssDNA) during DNA translocation and unwinding by RecBCD. The molecular determinants of χ recognition and the subsequent alteration in function are unknown. Consequently, we mutated residues within the RecC subunit that comprise a channel where ssDNA is thought to be scanned for a χ sequence. These mutants were characterized in vivo with regard to χ recognition, UV-sensitivity, phage degradation, and recombination proficiency. Of 38 residues mutated, 11 were previously undescribed mutations that altered χ recognition. The mutants fell into two classes: five that failed to respond to χ, and six that suggested a relaxed specificity for χ recognition. The location of the first set of mutations defines a recognition structure responsible for sequence-specific binding of ssDNA. The second set defines a highly conserved structure, linked to the recognition structure, which we hypothesize regulates conversion of RecBCD from a molecular machine that destroys DNA to one that repairs it. These findings offer insight into the evolution of enzymes with alternate χ recognition specificities.

Co-evolution of segregation guide DNA motifs and the FtsK translocase in bacteria: identification of the atypical Lactococcus lactis KOPS motif

Bacteria use the global bipolarization of their chromosomes into replichores to control the dynamics and segregation of their genome during the cell cycle. This involves the control of protein activities by recognition of specific short DNA motifs whose orientation along the chromosome is highly skewed. The KOPS motifs act in chromosome segregation by orienting the activity of the FtsK DNA translocase towards the terminal replichore junction. KOPS motifs have been identified in γ-Proteobacteria and in Bacillus subtilis as closely related G-rich octamers. We have identified the KOPS motif of Lactococcus lactis, a model bacteria of the Streptococcaceae family harbouring a compact and low GC% genome. This motif, 5′-GAAGAAG-3, was predicted in silico using the occurrence and skew characteristics of known KOPS motifs. We show that it is specifically recognized by L. lactis FtsK in vitro and controls its activity in vivo. L. lactis KOPS is thus an A-rich heptamer motif. Our results show that KOPS-controlled chromosome segregation is conserved in Streptococcaceae but that KOPS may show important variation in sequence and length between bacterial families. This suggests that FtsK adapts to its host genome by selecting motifs with convenient occurrence frequencies and orientation skews to orient its activity.

Insights into Chi recognition from the structure of an AddAB-type helicase-nuclease complex

Homologous recombination DNA repair requires double-strand break resection by helicase–nuclease enzymes. The crystal structure of bacterial AddAB in complex with DNA substrates shows that it employs an inactive helicase site to recognize ‘Chi' recombination hotspot sequences that regulate resection.

In bacterial cells, processing of double-stranded DNA breaks for repair by homologous recombination is dependent upon the recombination hotspot sequence Chi and is catalysed by either an AddAB- or RecBCD-type helicase–nuclease. Here, we report the crystal structure of AddAB bound to DNA. The structure allows identification of a putative Chi-recognition site in an inactivated helicase domain of the AddB subunit. By generating mutant protein complexes that do not respond to Chi, we show that residues responsible for Chi recognition are located in positions equivalent to the signature motifs of a conventional helicase. Comparison with the related RecBCD complex, which recognizes a different Chi sequence, provides further insight into the structural basis for sequence-specific ssDNA recognition. The structure suggests a simple mechanism for DNA break processing, explains how AddAB and RecBCD can accomplish the same overall reaction with different sets of functional modules and reveals details of the role of an Fe–S cluster in protein stability and DNA binding.

Robustness assessment of whole bacterial genome segmentations

Comparison of closely related bacterial genomes has revealed the presence of highly conserved sequences forming a "backbone" that is interrupted by numerous, less conserved, DNA fragments. Segmentation of bacterial genomes into backbone and variable regions is particularly useful to investigate, among other things, bacterial genome evolution. Several software tools have been designed to compare complete bacterial chromosomes and a few online databases store pre-computed genome comparisons. However, very few statistical methods are available to evaluate the reliability of these software tools and to compare the results obtained with them. To fill this gap, we have developed two local scores to measure the robustness of bacterial genome segmentations. Our method uses a simulation procedure based on random perturbations of the compared genomes. The two scores described in this article provide useful information and are easy to implement, and their interpretation is intuitive. We show that they are suited to discriminate between robust and non-robust segmentations when genome aligners such as MAUVE and MGA are used.

An alternative approach to multiple genome comparison

Genome comparison is now a crucial step for genome annotation and identification of regulatory motifs. Genome comparison aims for instance at finding genomic regions either specific to or in one-to-one correspondance between individuals/strains/species. It serves e.g. to pre-annotate a new genome by automatically transfering annotations from a known one. However, efficiency, flexibility and objectives of current methods do not suit the whole spectrum of applications, genome sizes and organizations. Innovative approaches are still needed. Hence, we propose an alternative way of comparing multiple genomes based on segmentation by similarity. In this framework, rather than being formulated as a complex optimization problem, genome comparison is seen as a segmentation question for which a single optimal solution can be found in almost linear time. We apply our method to analyse three strains of a virulent pathogenic bacteria, Ehrlichia ruminantium, and identify 92 new genes. We also find out that a substantial number of genes thought to be strain specific have potential orthologs in the other strains. Our solution is implemented in an efficient program, qod, equipped with a user-friendly interface, and enables the automatic transfer of annotations betwen compared genomes or contigs (Video in Supplementary Data). Because it somehow disregards the relative order of genomic blocks, qod can handle unfinished genomes, which due to the difficulty of sequencing completion may become an interesting characteristic for the future. Availabilty: http://www.atgc-montpellier.fr/qod.

Comparative supragenomic analyses among the pathogens Staphylococcus aureus, Streptococcus pneumoniae, and Haemophilus influenzae using a modification of the finite supragenome model

BACKGROUND

Staphylococcus aureus is associated with a spectrum of symbiotic relationships with its human host from carriage to sepsis and is frequently associated with nosocomial and community-acquired infections, thus the differential gene content among strains is of interest.

RESULTS

We sequenced three clinical strains and combined these data with 13 publically available human isolates and one bovine strain for comparative genomic analyses. All genomes were annotated using RAST, and then their gene similarities and differences were delineated. Gene clustering yielded 3,155 orthologous gene clusters, of which 2,266 were core, 755 were distributed, and 134 were unique. Individual genomes contained between 2,524 and 2,648 genes. Gene-content comparisons among all possible S. aureus strain pairs (n = 136) revealed a mean difference of 296 genes and a maximum difference of 476 genes. We developed a revised version of our finite supragenome model to estimate the size of the S. aureus supragenome (3,221 genes, with 2,245 core genes), and compared it with those of Haemophilus influenzae and Streptococcus pneumoniae. There was excellent agreement between RAST's annotations and our CDS clustering procedure providing for high fidelity metabolomic subsystem analyses to extend our comparative genomic characterization of these strains.

CONCLUSIONS

Using a multi-species comparative supragenomic analysis enabled by an improved version of our finite supragenome model we provide data and an interpretation explaining the relatively larger core genome of S. aureus compared to other opportunistic nasopharyngeal pathogens. In addition, we provide independent validation for the efficiency and effectiveness of our orthologous gene clustering algorithm.

Impact of ciprofloxacin exposure on Staphylococcus aureus genomic alterations linked with emergence of rifampin resistance

Intensive use of antimicrobial agents in health care settings not only leads to the selection of multiresistant nosocomial isolates of Staphylococcus aureus but may also promote endogenous, resistance-conferring mutations in bacterial genes that encode drug targets. We evaluated the spectrum of rifampin resistance-conferring mutations in cultures of methicillin-susceptible S. aureus (MSSA) or methicillin-resistant S. aureus (MRSA) strains exposed in vitro to sub-MICs of ciprofloxacin. Growth of ciprofloxacin-susceptible MRSA strain MRGR3 and ciprofloxacin-resistant MSSA strain RA1 (a NCTC 8325 derivative) in the presence of 1/2× or 1/4× MIC of ciprofloxacin led to higher frequencies of rifampin-resistant mutants on agar supplemented with rifampin (0.25 mg/liter) than under ciprofloxacin-free conditions. While rifampin-resistant mutants from ciprofloxacin-free cultures essentially showed single-amino-acid substitutions, a significant proportion of rifampin-resistant mutants from ciprofloxacin-exposed cultures displayed in-frame deletions or insertions in the rpoB gene at several positions of the rifampin resistance cluster I. In-frame deletions or insertions were also recorded in rpoB cluster I of rifampin-resistant mutants from ciprofloxacin-exposed cultures of mutS and mutL DNA repair mutants of ciprofloxacin-resistant S. aureus strain RA1. Frequencies of rifampin-resistant mutants grown under ciprofloxacin-free medium were higher for mutant strains RA1 mutS2 and RA1 mutL, but not RA1 recA, than for their parent RA1. In conclusion, ciprofloxacin-mediated DNA damage in S. aureus, as exemplified by the wide diversity of deletions or insertions in rpoB, suggests the occurrence of major, quinolone-mediated disturbances in DNA fork progression and replication repair. Besides promoting antibiotic resistance, accumulation of unrepaired DNA replication errors, including insertions and deletions, may also contribute to potentially lethal mutations.

Isolation of highly heat-resistant Listeria monocytogenes variants by use of a kinetic modeling-based sampling scheme

Stable high-hydrostatic-pressure (HHP)-resistant Listeria monocytogenes LO28 variants were previously isolated and characterized. These HHP variants were also more resistant to heat. In addition, nonlinear heat inactivation kinetics pointed toward the existence of heat-resistant variants, although these could not be isolated so far. In this study, we used kinetic modeling of inactivation curves of two isolated HHP variants and their wild type, and this revealed that the probability of finding resistant variants should depend on the nature of the inactivation treatment and the time of exposure. At specific heat and HHP conditions, resistant LO28 and EGDe variants were indeed isolated. Resistant LO28 variants were even isolated after a heat inactivation at 72°C in milk, and these variants showed high resistance to standard pasteurization conditions. The increased resistance of part of the isolated LO28 and EGDe variants was due to mutations in their ctsR genes. For the variants whose ctsR genes and upstream regions were not altered, the mechanisms leading to increased resistance remain to be elucidated. This research showed the strength of kinetic modeling in unraveling the causes of nonlinear inactivation and facilitating the isolation of heat-resistant L. monocytogenes variants.

DNA motifs that sculpt the bacterial chromosome

During the bacterial cell cycle, the processes of chromosome replication, DNA segregation, DNA repair and cell division are coordinated by precisely defined events. Tremendous progress has been made in recent years in identifying the mechanisms that underlie these processes. A striking feature common to these processes is that non-coding DNA motifs play a central part, thus 'sculpting' the bacterial chromosome. Here, we review the roles of these motifs in the mechanisms that ensure faithful transmission of genetic information to daughter cells. We show how their chromosomal distribution is crucial for their function and how it can be analysed quantitatively. Finally, the potential roles of these motifs in bacterial chromosome evolution are discussed.

The processing of double-stranded DNA breaks for recombinational repair by helicase-nuclease complexes

Double-stranded DNA breaks are prepared for recombinational repair by nucleolytic digestion to form single-stranded DNA overhangs that are substrates for RecA/Rad51-mediated strand exchange. This processing can be achieved through the activities of multiple helicases and nucleases. In bacteria, the function is mainly provided by a stable multi-protein complex of which there are two structural classes; AddAB- and RecBCD-type enzymes. These helicase-nucleases are of special interest with respect to DNA helicase mechanism because they are exceptionally powerful DNA translocation motors, and because they serve as model systems for both single molecule studies and for understanding how DNA helicases can be coupled to other protein machinery. This review discusses recent developments in our understanding of the AddAB and RecBCD complexes, focussing on their distinctive strategies for processing DNA ends. We also discuss the extent to which bacterial DNA end resection mechanisms may parallel those used in eukaryotic cells.

Genome analysis and genome-wide proteomics of Thermococcus gammatolerans, the most radioresistant organism known amongst the Archaea

The genome sequence of Thermococcus gammatolerans, a radioresistant archaeon, is described; a proteomic analysis reveals that radioresistance may be due to unknown DNA repair enzymes.

Background

Thermococcus gammatolerans was isolated from samples collected from hydrothermal chimneys. It is one of the most radioresistant organisms known amongst the Archaea. We report the determination and annotation of its complete genome sequence, its comparison with other Thermococcales genomes, and a proteomic analysis.

Results

T. gammatolerans has a circular chromosome of 2.045 Mbp without any extra-chromosomal elements, coding for 2,157 proteins. A thorough comparative genomics analysis revealed important but unsuspected genome plasticity differences between sequenced Thermococcus and Pyrococcus species that could not be attributed to the presence of specific mobile elements. Two virus-related regions, tgv1 and tgv2, are the only mobile elements identified in this genome. A proteogenome analysis was performed by a shotgun liquid chromatography-tandem mass spectrometry approach, allowing the identification of 10,931 unique peptides corresponding to 951 proteins. This information concurrently validates the accuracy of the genome annotation. Semi-quantification of proteins by spectral count was done on exponential- and stationary-phase cells. Insights into general catabolism, hydrogenase complexes, detoxification systems, and the DNA repair toolbox of this archaeon are revealed through this genome and proteome analysis.

Conclusions

This work is the first archaeal proteome investigation done at the stage of primary genome annotation. This archaeon is shown to use a large variety of metabolic pathways even under a rich medium growth condition. This proteogenomic study also indicates that the high radiotolerance of T. gammatolerans is probably due to proteins that remain to be characterized rather than a larger arsenal of known DNA repair enzymes.

MOSAIC: an online database dedicated to the comparative genomics of bacterial strains at the intra-species level

BACKGROUND

The recent availability of complete sequences for numerous closely related bacterial genomes opens up new challenges in comparative genomics. Several methods have been developed to align complete genomes at the nucleotide level but their use and the biological interpretation of results are not straightforward. It is therefore necessary to develop new resources to access, analyze, and visualize genome comparisons.

DESCRIPTION

Here we present recent developments on MOSAIC, a generalist comparative bacterial genome database. This database provides the bacteriologist community with easy access to comparisons of complete bacterial genomes at the intra-species level. The strategy we developed for comparison allows us to define two types of regions in bacterial genomes: backbone segments (i.e., regions conserved in all compared strains) and variable segments (i.e., regions that are either specific to or variable in one of the aligned genomes). Definition of these segments at the nucleotide level allows precise comparative and evolutionary analyses of both coding and non-coding regions of bacterial genomes. Such work is easily performed using the MOSAIC Web interface, which allows browsing and graphical visualization of genome comparisons.

CONCLUSION

The MOSAIC database now includes 493 pairwise comparisons and 35 multiple maximal comparisons representing 78 bacterial species. Genome conserved regions (backbones) and variable segments are presented in various formats for further analysis. A graphical interface allows visualization of aligned genomes and functional annotations. The MOSAIC database is available online at http://genome.jouy.inra.fr/mosaic.

The MatP/matS site-specific system organizes the terminus region of the E. coli chromosome into a macrodomain

The organization of the Escherichia coli chromosome into insulated macrodomains influences the segregation of sister chromatids and the mobility of chromosomal DNA. Here, we report that organization of the Terminus region (Ter) into a macrodomain relies on the presence of a 13 bp motif called matS repeated 23 times in the 800-kb-long domain. matS sites are the main targets in the E. coli chromosome of a newly identified protein designated MatP. MatP accumulates in the cell as a discrete focus that colocalizes with the Ter macrodomain. The effects of MatP inactivation reveal its role as main organizer of the Ter macrodomain: in the absence of MatP, DNA is less compacted, the mobility of markers is increased, and segregation of Ter macrodomain occurs early in the cell cycle. Our results indicate that a specific organizational system is required in the Terminus region for bacterial chromosome management during the cell cycle.

Faster exact Markovian probability functions for motif occurrences: a DFA-only approach

BACKGROUND

The computation of the statistical properties of motif occurrences has an obviously relevant application: patterns that are significantly over- or under-represented in genomes or proteins are interesting candidates for biological roles. However, the problem is computationally hard; as a result, virtually all the existing motif finders use fast but approximate scoring functions, in spite of the fact that they have been shown to produce systematically incorrect results. A few interesting exact approaches are known, but they are very slow and hence not practical in the case of realistic sequences.

RESULTS

We give an exact solution, solely based on deterministic finite-state automata (DFA), to the problem of finding the whole relevant part of the probability distribution function of a simple-word motif in a homogeneous (biological) sequence. Out of that, the z-value can always be computed, while the P-value can be obtained either when it is not too extreme with respect to the number of floating-point digits available in the implementation, or when the number of pattern occurrences is moderately low. In particular, the time complexity of the algorithms for Markov models of moderate order (0 < or = m < or = 2) is far better than that of Nuel, which was the fastest similar exact algorithm known to date; in many cases, even approximate methods are outperformed.

CONCLUSIONS

DFA are a standard tool of computer science for the study of patterns; previous works in biology propose algorithms involving automata, but there they are used, respectively, as a first step to write a generating function, or to build a finite Markov-chain imbedding (FMCI). In contrast, we directly rely on DFA to perform the calculations; thus we manage to obtain an algorithm which is both easily interpretable and efficient. This approach can be used for exact statistical studies of very long genomes and protein sequences, as we illustrate with some examples on the scale of the human genome.

FtsK-dependent dimer resolution on multiple chromosomes in the pathogen Vibrio cholerae

Unlike most bacteria, Vibrio cholerae harbors two distinct, nonhomologous circular chromosomes (chromosome I and II). Many features of chromosome II are plasmid-like, which raised questions concerning its chromosomal nature. Plasmid replication and segregation are generally not coordinated with the bacterial cell cycle, further calling into question the mechanisms ensuring the synchronous management of chromosome I and II. Maintenance of circular replicons requires the resolution of dimers created by homologous recombination events. In Escherichia coli, chromosome dimers are resolved by the addition of a crossover at a specific site, dif, by two tyrosine recombinases, XerC and XerD. The process is coordinated with cell division through the activity of a DNA translocase, FtsK. Many E. coli plasmids also use XerCD for dimer resolution. However, the process is FtsK-independent. The two chromosomes of the V. cholerae N16961 strain carry divergent dimer resolution sites, dif1 and dif2. Here, we show that V. cholerae FtsK controls the addition of a crossover at dif1 and dif2 by a common pair of Xer recombinases. In addition, we show that specific DNA motifs dictate its orientation of translocation, the distribution of these motifs on chromosome I and chromosome II supporting the idea that FtsK translocation serves to bring together the resolution sites carried by a dimer at the time of cell division. Taken together, these results suggest that the same FtsK-dependent mechanism coordinates dimer resolution with cell division for each of the two V. cholerae chromosomes. Chromosome II dimer resolution thus stands as a bona fide chromosomal process.

During proliferation, DNA synthesis, chromosome segregation, and cell division must be coordinated to ensure the stable inheritance of the genetic material. In eukaryotes, this is achieved by checkpoint mechanisms that delay certain steps until others are completed. No such temporal separation exists in bacteria, which can undergo overlapping replication cycles. The eukaryotic cell cycle is particularly well suited to the management of multiple chromosomes, with the same replication initiation and segregation machineries operating on all the chromosomes, while the bacterial cell cycle is linked to genomes of less complexity, most bacteria harboring a single chromosome. The discovery of bacteria harboring multiple circular chromosomes, such as V. cholerae, raised therefore a considerable interest for the mechanisms ensuring the synchronous management of different replicons. Here, we took advantage of our knowledge of chromosome dimer resolution, the only bacterial segregation process for which coordination with cell division is well understood, to investigate one of the mechanisms ensuring the synchronous management of the smaller, plasmid-like, and larger, chromosome-like, replicons of V. cholerae.