I-CeuI recognition sites in the rrn operons of the Bacillus subtilis 168 chromosome: inherent landmarks for genome analysis

The linear chromosome of the plant-pathogenic mycoplasma 'Candidatus Phytoplasma mali'

Background

Phytoplasmas are insect-transmitted, uncultivable bacterial plant pathogens that cause diseases in hundreds of economically important plants. They represent a monophyletic group within the class Mollicutes (trivial name mycoplasmas) and are characterized by a small genome with a low GC content, and the lack of a firm cell wall. All mycoplasmas, including strains of 'Candidatus (Ca.) Phytoplasma asteris' and 'Ca. P. australiense', examined so far have circular chromosomes, as is the case for almost all walled bacteria.

Results

Our work has shown that 'Ca. Phytoplasma mali', the causative agent of apple proliferation disease, has a linear chromosome. Linear chromosomes were also identified in the closely related provisional species 'Ca. P. pyri' and 'Ca. P. prunorum'. The chromosome of 'Ca. P. mali' strain AT is 601,943 bp in size and has a GC content of 21.4%. The chromosome is further characterized by large terminal inverted repeats and covalently closed hairpin ends. Analysis of the protein-coding genes revealed that glycolysis, the major energy-yielding pathway supposed for 'Ca. P. asteris', is incomplete in 'Ca. P. mali'. Due to the apparent lack of other metabolic pathways present in mycoplasmas, it is proposed that maltose and malate are utilized as carbon and energy sources. However, complete ATP-yielding pathways were not identified. 'Ca. P. mali' also differs from 'Ca. P. asteris' by a smaller genome, a lower GC content, a lower number of paralogous genes, fewer insertions of potential mobile DNA elements, and a strongly reduced number of ABC transporters for amino acids. In contrast, 'Ca. P. mali' has an extended set of genes for homologous recombination, excision repair and SOS response than 'Ca. P. asteris'.

Conclusion

The small linear chromosome with large terminal inverted repeats and covalently closed hairpin ends, the extremely low GC content and the limited metabolic capabilities reflect unique features of 'Ca. P. mali', not only within phytoplasmas, but all mycoplasmas. It is expected that the genome information obtained here will contribute to a better understanding of the reduced metabolism of phytoplasmas, their fastidious nutrition requirements that prevented axenic cultivation, and the mechanisms involved in pathogenicity.

Combining two genomes in one cell: stable cloning of the Synechocystis PCC6803 genome in the Bacillus subtilis 168 genome

Cloning the whole 3.5-megabase (Mb) genome of the photosynthetic bacterium Synechocystis PCC6803 into the 4.2-Mb genome of the mesophilic bacterium Bacillus subtilis 168 resulted in a 7.7-Mb composite genome. We succeeded in such unprecedented large-size cloning by progressively assembling and editing contiguous DNA regions that cover the entire Synechocystis genome. The strain containing the two sets of genome grew only in the B. subtilis culture medium where all of the cloning procedures were carried out. The high structural stability of the cloned Synechocystis genome was closely associated with the symmetry of the bacterial genome structure of the DNA replication origin (oriC) and its termination (terC) and the exclusivity of Synechocystis ribosomal RNA operon genes (rrnA and rrnB). Given the significant diversity in genome structure observed upon horizontal DNA transfer in nature, our stable laboratory-generated composite genome raised fundamental questions concerning two complete genomes in one cell. Our megasize DNA cloning method, designated megacloning, may be generally applicable to other genomes or genome loci of free-living organisms.

Comparative analysis of physical maps of four Bacillus subtilis (natto) genomes

The complete SfiI and I-CeuI physical maps of four Bacillus subtilis (natto) strains, which were previously isolated as natto (fermented soybean) starters, were constructed to elucidate the genome structure. Not only the similarity in genome size and organization but also the microheterogeneity of the gene context was revealed. No large-scale genome rearrangements among the four strains were indicated by mapping of the genes, including 10 rRNA operons (rrn) and relevant genes required for natto production, to the loci corresponding to those of the B. subtilis strain Marburg 168. However, restriction fragment length polymorphism and the presence or absence of strain-specific DNA sequences, such as the prophages SP beta, skin element, and PBSX, as well as the insertion element IS4Bsu1, could be used to identify one of these strains as a Marburg type and the other three strains as natto types. The genome structure and gene heterogeneity were also consistent with the type of indigenous plasmids harbored by the strains.

Genome plasticity among related ++Lactococcus strains: identification of genetic events associated with macrorestriction polymorphisms

The genomic diversity of nine strains of the Lactococcus lactis subsp. cremoris (NCDO712, NCDO505, NCDO2031, NCDO763, MMS36, C2, LM0230, LM2301, and MG1363) was studied by macrorestriction enzyme analysis using pulsed-field gel electrophoresis. These strains were considered adequate for the investigation of genomic plasticity because they have been described as belonging to the same genetic lineage. Comparison of ApaI and SmaI genome fingerprints of each strain revealed the presence of several macrorestriction fragment length polymorphisms (RFLPs), despite a high degree of similarity of the generated restriction patterns. The physical map of the MG1363 chromosome was used to establish a genome map of the other strains and allocate the RFLPs to five regions. Southern hybridization analysis correlated the polymorphic regions with genetic events such as chromosomal inversion, integration of prophage DNA, and location of the transposon-like structures carrying conjugative factor or oligopeptide transport system.

Chromosome mapping of the sweet potato little leaf phytoplasma reveals genome heterogeneity within the phytoplasmas

To further understand the genomic diversity and genetic architecture of phytoplasmas, a physical and genetic map of the sweet potato little leaf (SPLL) strain V4 phytoplasma chromosome was determined. PFGE was used to determine the size of the SPLL-V4 genome, which was estimated to be 622 kb. A physical map was prepared by two-dimensional reciprocal digestions using the restriction endonucleases BssHII, Smal, Eagl and I-Ceul. Sixteen cleavage sites were located on the map. Southern hybridizations of digested SPLL-V4 chromosomal DNA were done using random clones and PCR-amplified genes as probes. This confirmed fragment positions and located the two rRNA operons and the linked fus/tuf genes encoding elongation factors G and Tu, respectively, on the physical map. An inversion of one of the rRNA operons was observed from hybridization data. Sequence analysis of one of the random clones identified a gid gene encoding a glucose-inhibited division protein. Digestions of the tomato big bud (TBB) phytoplasma chromosome with the same four enzymes revealed genome heterogeneity when compared to the closely related SPLL-V4, and a preliminary chromosome size for the TBB phytoplasma of 662 kb was estimated. This mapping information has revealed that significant genome diversity exists within the phytoplasmas.

The molecular biology of pasteurella multocida

Pasteurella multocida is an important veterinary and opportunistic human pathogen. The species is diverse and complex with respect to antigenic variation, host predeliction and pathogenesis. Certain serological types are the aetiologic agents of severe pasteurellosis, such as fowl cholera in domestic and wild birds, bovine haemorrhagic septicaemia and porcine atrophic rhinitis. The recent application of molecular methods such as the polymerase chain reaction, restriction endonuclease analysis, ribotyping, pulsed-field gel electrophoresis, gene cloning, characterisation and recombinant protein expression, mutagenesis, plasmid and bacteriophage analysis and genomic mapping, have greatly increased our understanding of P. multocida and has provided researchers with a number of molecular tools to study pathogenesis and epidemiology at a molecular level.

Fate of unstable Bacillus subtilis subgenome: re-integration and amplification in the main genome

The plastic Bacillus subtilis genome was dissected into two physically separate genomes, the 3.9 Mb main genome and the 0.3 Mb subgenome. DNA replication of the main genome was initiated from the normal replication origin (oriC) and that of the subgenome was from a 7.2 kb oriN-containing fragment artificially inserted. When the 7.2 kb fragment was shortened to a 1.5 kb fragment that contains oriN but lacks the segregational function, the subgenome became unstable and was rapidly lost from the cell, producing inviable cells due to the loss of essential genes carried by the subgenome. Stable survivors were isolated in which the subgenome had re-integrated and multiplied in the main genome. These results suggest that a reduced genetic stability of the subgenome induces size variation of the B. subtilis genome.

Construction of a physical and preliminary genetic map of Aeromonas hydrophila JMP636

A physical and preliminary genetic map of the Aeromonas hydrophila JMP636 chromosome has been constructed. The topology of the genome was predicted to be circular as chromosomal DNA did not migrate from the origin during PFGE unless linearized by S1 nuclease. Cleavage of the chromosome with PacI and PmeI produced 23 and 14 fragments, respectively, and enabled calculation of the genome size at 4.5 Mb. Digestion of the chromosome with I-CeuI produced 10 fragments, indicating that 10 rrl (23S) genes were likely to be present. Hybridizations between DNA fragments generated with PacI, PmeI and I-CeuI were used to initially determine the relationship between these segments. To accurately map genes previously characterized from JMP636, the suicide vector pJP5603 was modified to introduce restriction sites for PacI and PmeI, producing pJP9540. Following cloning of genes into this vector and recombinational insertion into the JMP636 chromosome, PacI and PmeI cleavage determined the location of genes within macrorestriction fragments with the additional bands produced forming hybridization probes. From the data generated, it was possible to form a physical map comprising all the fragments produced by PacI and PmeI, and assign the contig of I-CeuI fragments on this map. The preliminary genetic map defines the location of six loci for degradative enzymes previously characterized from JMP636, while the locations of the 10 sets of ribosomal genes were assigned with less accuracy from hybridization data.

Physical and genetic map of the Pasteurella multocida A:1 chromosome

A physical and genetic map of the Pasteurella multocida A:1 genome was generated by using the restriction enzymes ApaI, CeuI, and NotI. The positions of 23 restriction sites and 32 genes, including 5 rrn operons, were localized on the 2.35-Mbp single circular chromosome. This report presents the first genetic and physical map for this genus.

A natural large chromosomal inversion in Lactococcus lactis is mediated by homologous recombination between two insertion sequences

Comparative analysis of chromosomal macrorestriction polymorphism of the two closely related Lactococcus lactis subsp. cremoris strains MG1363 and NCDO763 revealed the presence of a large inversion covering half of the genome. To determine what kind of genetic element could be implicated in this rearrangement, the two inversion junctions of MG1363 and NCDO763 chromosomes were cloned and characterized. Nucleotide sequence analysis showed the presence of one copy of the lactococcal IS905 element in each junction. Each copy of this element contained the same nucleotide mutation that inactivates the putative transposase. Comparison of the sequences surrounding the insertion sequence demonstrated that the large inversion arose from a single-step homologous recombination event between the two defective copies of the IS905 element. The large inversion presumably conferred no selective disadvantage on strain NCDO763 because this rearrangement did not alter the oriC-terC symmetry of the chromosome and the local genetic environment.

Physical map of the Bacillus subtilis 166 genome: evidence for the inversion of an approximately 1900 kb continuous DNA segment, the translocation of an approximately 100 kb segment and the duplication of a 5 kb segment

An I-CeuI-NotI-SfiI endonuclease map of the Bacillus subtilis 166 genome was constructed. It was almost identical to that of B. subtilis 168 except for the inversion of an approximately 1900 kb DNA segment, the translocation of an approximately 100 kb segment and the duplication of a 5 kb segment. Continuity of the inverted segment was investigated by direct measurement of the distances between the two genomic loci where I-SceI recognition sites were created in the 168 and the 166 genomes. Size difference of the I-SceI fragments between the two strains fully demonstrated the inversion of an approximately 1900 kb long 'continuous' DNA segment and the 'location' of the two inversion junctions in the genome. The 100 kb DNA segment including the lysogenic SP beta prophage was translocated close to one of the inversion junctions and was probably associated with the duplication of a 5 kb segment. These rearrangements are consistent with those indicated by genetic analyses.

Genomic fluidity of Bordetella pertussis assessed by a new method for chromosomal mapping

The genomic organization of Bordetella pertussis strains has been examined by using a new method. This method does not depend on the prior determination of a restriction map of the bacterial chromosome but is based on the ability to measure directly the distance between two genes. This is accomplished through the integration at each gene of a suicide vector containing a cleavage site for the intron-encoded endonuclease I-SceI, which is not otherwise found in the chromosome. Integration is mediated by homologous recombination between the chromosomal and cloned plasmid copies of a gene of interest. Digestion with I-SceI gives rise to a fragment the size of which represents the distance between the two genes. Multiple pairwise determinations within a set of genes provide sufficient information to derive a map of the relative gene positions. Mapping a set of 11 to 13 genes for five strains of B. pertussis and one strain of B. parapertussis revealed extensive divergence of gene order between B. pertussis Tohama I, B. pertussis 18-323, and B. parapertussis ATCC 15311. Less extensive divergence of gene order was observed between B. pertussis Tohama I and B. pertussis Tohama III, BP165, and Wellcome 28, with most of the observed differences explainable by large inversions.

Homing endonucleases: keeping the house in order

Homing endonucleases are rare-cutting enzymes encoded by introns and inteins. They have striking structural and functional properties that distinguish them from restriction enzymes. Nomenclature conventions analogous to those for restriction enzymes have been developed for the homing endonucleases. Recent progress in understanding the structure and function of the four families of homing enzymes is reviewed. Of particular interest are the first reported structures of homing endonucleases of the LAGLIDADG family. The exploitation of the homing enzymes in genome analysis and recombination research is also summarized. Finally, the evolution of homing endonucleases is considered, both at the structure-function level and in terms of their persistence in widely divergent biological systems.

Evolutionarily conserved and functionally important residues in the I-CeuI homing endonuclease

Two approaches were used to discern critical amino acid residues for the function of the I- Ceu I homing endonuclease: sequence comparison of subfamilies of homologous proteins and genetic selection. The first approach revealed residues potentially involved in catalysis and DNA recognition. Because I- Ceu I is lethal in Escherichia coli , enzyme variants not perturbing cell viability were readily selected from an expression library. A collection of 49 variants with single amino acid substitutions at 37 positions was assembled. Most of these positions are clustered within or around the LAGLI-DADG dodecapeptide and the TQH sequence, two motifs found in all protein subfamilies examined. The Km and kcat values of the wild-type and nine variant enzymes synthesized in vitro were determined. Three variants, including one showing a substitution of the glutamine residue in the TQH motif, revealed no detectable endonuclease activity; five others showed reduced activity compared to the wild-type enzyme; whereas the remaining variant cleaved the top strand about three times more efficiently than the wild-type. Our results not only confirm recent reports indicating that amino acids in the LAGLI-DADG dodecapeptide are functionally critical, but they also suggest that some residues outside this motif directly participate in catalysis.

The Bacillus subtilis 168 chromosome from sspE to katA

We have cloned and sequenced a 24.5 kb region of the Bacillus subtilis 168 chromosome spanning the sspE and katA genes. The region contains a ribosomal RNA operon, rrnD, a tRNA gene set, trnD and 17 ORFs, 16 with putative ribosome-binding sites. Four of the ORFs (ORF2, ORF14, ORF16 and ORF17) match to known B. subtilis genes (sspE, thiA, senS and katA). Eight of the remaining ORF products show similarities with proteins present in the databases, including an ATP-binding transport protein, a glutamate-1-semialdehyde aminotransferase, a thiol-specific antioxidant protein, a mitomycin radical oxidase and a ferric uptake regulation protein.

Experimental surgery to create subgenomes of Bacillus subtilis 168

The 4,188-kb circular genome of Bacillus subtilis 168 was artificially dissected into two stable circular chromosomes in vivo, one being the 3,878-kb main genome and the other the 310-kb subgenome that was recovered as covalently closed circular DNA in CsCl-ethidium bromide ultracentrifugation. The minimal requirements to physically separate the 310-kb DNA segment out of the genome were two interrepeat homologous sequences and an origin of DNA replication between them. The subgenome originated from the 1,255-1, 551-kb region of the B. subtilis genome was essential for the cell to survive because the subgenome was not lost from the cell. The finding that the B. subtilis genome has a potential to be divided and the resulting two replicons stably maintained may shed light on origins and formation mechanisms of giant plasmids or second chromosomes present in many bacteria. Similar excision or its reversal process, i.e., integration of large sized covalently closed circular DNA pieces into the main genome, implies significant roles of subgenomes in the exchange of genetic information and size variation of bacterial genomes in bacterial evolution.