Genometrics as an essential tool for the assembly of whole genome sequences: the example of the chromosome of Bifidobacterium longum NCC2705

Integration of the pSLT Plasmid into the Salmonella Chromosome Results in a Temperature-Sensitive Growth Defect Due to Aberrant DNA Replication

A mutant of Salmonella enterica serovar Typhimurium was isolated that simultaneously affected two metabolic pathways as follows: NAD metabolism and DNA repair. The mutant was isolated as resistant to a nicotinamide analog and as temperature-sensitive for growth on minimal glucose medium. In this mutant, Salmonella's 94-kb virulence plasmid pSLT had recombined into the chromosome upstream of the NAD salvage pathway gene pncA This insertion blocked most transcription of pncA, which reduced uptake of the nicotinamide analog. The pSLT insertion mutant also exhibited phenotypes associated with induction of the SOS DNA repair system, including an increase in filamentous cells, higher exonuclease III and catalase activities, and derepression of SOS gene expression. Genome sequencing revealed increased read coverage extending out from the site of pSLT insertion. The two pSLT replication origins are likely initiating replication of the chromosome near the normal replication terminus. Too much replication initiation at the wrong site is probably causing the observed growth defects. Accordingly, deletion of both pSLT replication origins restored growth at higher temperatures.IMPORTANCE In studies that insert a second replication origin into the chromosome, both origins are typically active at the same time. In contrast, the integrated pSLT plasmid initiated replication in stationary phase after normal chromosomal replication had finished. The gradient in read coverage extending out from a single site could be a simple but powerful tool for studying replication and detecting chromosomal rearrangements. This technique may be of particular value when a genome has been sequenced for the first time to verify correct assembly.

Applying Shannon's information theory to bacterial and phage genomes and metagenomes

All sequence data contain inherent information that can be measured by Shannon's uncertainty theory. Such measurement is valuable in evaluating large data sets, such as metagenomic libraries, to prioritize their analysis and annotation, thus saving computational resources. Here, Shannon's index of complete phage and bacterial genomes was examined. The information content of a genome was found to be highly dependent on the genome length, GC content, and sequence word size. In metagenomic sequences, the amount of information correlated with the number of matches found by comparison to sequence databases. A sequence with more information (higher uncertainty) has a higher probability of being significantly similar to other sequences in the database. Measuring uncertainty may be used for rapid screening for sequences with matches in available database, prioritizing computational resources, and indicating which sequences with no known similarities are likely to be important for more detailed analysis.

Genomic insights into bifidobacteria

Since the discovery in 1899 of bifidobacteria as numerically dominant microbes in the feces of breast-fed infants, there have been numerous studies addressing their role in modulating gut microflora as well as their other potential health benefits. Because of this, they are frequently incorporated into foods as probiotic cultures. An understanding of their full interactions with intestinal microbes and the host is needed to scientifically validate any health benefits they may afford. Recently, the genome sequences of nine strains representing four species of Bifidobacterium became available. A comparative genome analysis of these genomes reveals a likely efficient capacity to adapt to their habitats, with B. longum subsp. infantis exhibiting more genomic potential to utilize human milk oligosaccharides, consistent with its habitat in the infant gut. Conversely, B. longum subsp. longum exhibits a higher genomic potential for utilization of plant-derived complex carbohydrates and polyols, consistent with its habitat in an adult gut. An intriguing observation is the loss of much of this genome potential when strains are adapted to pure culture environments, as highlighted by the genomes of B. animalis subsp. lactis strains, which exhibit the least potential for a gut habitat and are believed to have evolved from the B. animalis species during adaptation to dairy fermentation environments.

M-protein and other intrinsic virulence factors of Streptococcus pyogenes are encoded on an ancient pathogenicity island

BACKGROUND

The increasing number of completely sequenced bacterial genomes allows comparing their architecture and genetic makeup. Such new information highlights the crucial role of lateral genetic exchanges in bacterial evolution and speciation.

RESULTS

Here we analyzed the twelve sequenced genomes of Streptococcus pyogenes by a naïve approach that examines the preferential nucleotide usage along the chromosome, namely the usage of G versus C (GC-skew) and T versus A (TA-skew). The cumulative GC-skew plot presented an inverted V-shape composed of two symmetrical linear segments, where the minimum and maximum corresponded to the origin and terminus of DNA replication. In contrast, the cumulative TA-skew presented a V-shape, which segments were interrupted by several steep slopes regions (SSRs), indicative of a different nucleotide composition bias. Each S. pyogenes genome contained up to nine individual SSRs, encompassing all described strain-specific prophages. In addition, each genome contained a similar unique non-phage SSR, the core of which consisted of 31 highly homologous genes. This core includes the M-protein, other mga-related factors and other virulence genes, totaling ten intrinsic virulence genes. In addition to a high content in virulence-related genes and to a peculiar nucleotide bias, this SSR, which is 47 kb-long in a M1GAS strain, harbors direct repeats and a tRNA gene, suggesting a mobile element. Moreover, its complete absence in a M-protein negative group A Streptococcus natural isolate demonstrates that it could be spontaneously lost, but in vitro deletion experiments indicates that its excision occurred at very low rate. The stability of this SSR, combined to its presence in all sequenced S. pyogenes sequenced genome, suggests that it results from an ancient acquisition.

CONCLUSION

Thus, this non-phagic SSR is compatible with a pathogenicity island, acquired before S. pyogenes speciation. Its potential excision might bear relevance for vaccine development, because vaccines targeting M-protein might select for M-protein-negative variants that still carry other virulence determinants.

Completely sequenced genomes of pathogenic bacteria: a review

Six out of ten completely sequenced bacterial genomes are pathogenic or opportunistic bacteria. The genome sequence of at least one strain of all the principal pathogenic bacteria will soon be available. This information should enable us to identify genes that encode virulence factors. As these genes are potential targets for drugs and vaccines, their identification should have considerable repercussions on prevention, diagnosis, and treatment of the main bacterial infectious diseases. Comparison of genome sequences of several strains of the same species should allow identification of the genetic clues responsible for the differing behavior of related bacterial pathogens. This article reviews the genomes from pathogenic bacteria that have been or are currently being sequenced, describes the main tasks to be accomplished after a genome sequence becomes available, and discusses the benefits of having the genome sequence of bacterial pathogens.

Separating the effects of mutation and selection in producing DNA skew in bacterial chromosomes

Background

Many bacterial chromosomes display nucleotide asymmetry, or skew, between the leading and lagging strands of replication. Mutational differences between these strands result in an overall pattern of skew that is centered about the origin of replication. Such a pattern could also arise from selection coupled with a bias for genes coded on the leading strand. The relative contributions of selection and mutation in producing compositional skew are largely unknown.

Results

We describe a model to quantify the contribution of mutational differences between the leading and lagging strands in producing replication-induced skew. When the origin and terminus of replication are known, the model can be used to estimate the relative accumulation of G over C and of A over T on the leading strand due to replication effects in a chromosome with bidirectional replication arms. The model may also be implemented in a maximum likelihood framework to estimate the locations of origin and terminus. We find that our estimations for the origin and terminus agree very well with the location of genes that are thought to be associated with the replication origin. This indicates that our model provides an accurate, objective method of determining the replication arms and also provides support for the hypothesis that these genes represent an ancestral cluster of origin-associated genes.

Conclusion

The model has several advantages over other methods of analyzing genome skew. First, it quantifies the role of mutation in generating skew so that its effect on composition, for example codon bias, can be assessed. Second, it provides an objective method for locating origin and terminus, one that is based on chromosome-wide accumulation of leading vs lagging strand nucleotide differences. Finally, the model has the potential to be utilized in a maximum likelihood framework in order to analyze the effect of chromosome rearrangements on nucleotide composition.

Solving ambiguities in contig assembly of Idiomarina loihiensis L2TR chromosome by in silico analyses

Nucleotide composition analyses of bacterial genomes such as cumulative GC skew highlight the atypical, strongly asymmetric architecture of the recently published chromosome of Idiomarina loihiensis L2TR, suggesting that an inversion of a 600-kb chromosomal segment occurred. The presence of 3.4-kb inverted repeated sequences at the borders of the putative rearrangement supports this hypothesis. Reverting in silico this segment restores (1) a symmetric chromosome architecture; (2) the co-orientation of transcription of all rRNA operons with DNA replication; and (3) a better conservation of gene order between this chromosome and other gamma-proteobacterial ones. Finally, long-range PCRs encompassing the ends of the 600-kb segment reveal the existence of the reverted configuration but not of the published one. This demonstrates how cumulative nucleotide-skew analyses can validate genome assemblies.

The clc element of Pseudomonas sp. strain B13, a genomic island with various catabolic properties

Pseudomonas sp. strain B13 is a bacterium known to degrade chloroaromatic compounds. The properties to use 3- and 4-chlorocatechol are determined by a self-transferable DNA element, the clc element, which normally resides at two locations in the cell's chromosome. Here we report the complete nucleotide sequence of the clc element, demonstrating the unique catabolic properties while showing its relatedness to genomic islands and integrative and conjugative elements rather than to other known catabolic plasmids. As far as catabolic functions, the clc element harbored, in addition to the genes for chlorocatechol degradation, a complete functional operon for 2-aminophenol degradation and genes for a putative aromatic compound transport protein and for a multicomponent aromatic ring dioxygenase similar to anthranilate hydroxylase. The genes for catabolic functions were inducible under various conditions, suggesting a network of catabolic pathway induction. For about half of the open reading frames (ORFs) on the clc element, no clear functional prediction could be given, although some indications were found for functions that were similar to plasmid conjugation. The region in which these ORFs were situated displayed a high overall conservation of nucleotide sequence and gene order to genomic regions in other recently completed bacterial genomes or to other genomic islands. Most notably, except for two discrete regions, the clc element was almost 100% identical over the whole length to a chromosomal region in Burkholderia xenovorans LB400. This indicates the dynamic evolution of this type of element and the continued transition between elements with a more pathogenic character and those with catabolic properties.