Tracing genomic variations in two highly virulent Yersinia enterocolitica strains with unequal ability to compete for host colonization

ABSTRACT: BACKGROUND: Yersinia enterocolitica is a gastrointestinal foodborne pathogen found worldwide and which especially affects infants and young children. While different bioserotypes have been associated with varying pathogenicity, research on Y. enterocolitica is mainly conducted on the highly virulent mouse-lethal strains of biotype 1B and serotype O:8. We demonstrate here that two Y. enterocolitica bioserotype 1B/O:8 strains, 8081 and WA-314, display different virulence and fitness properties in a mouse model. In vivo co-infection experiments revealed that strain WA-314 overcomes strain 8081 in the colonization of spleen and liver. To trace the reasons of this incongruity, we present here the first high-quality sequence of the whole genome of strain WA-314 and compare it to the published genome of strain 8081. RESULTS: Regions previously accepted as unique to strain 8081, like the YAPI and YGI-3 genomic islands, are absent from strain WA-314, confirming their strain-specificity. On the other hand, some fitness- and bacterial competition-associated features, such as a putative colicin cluster and a xenobiotic-acyltransferase-encoding gene, are unique to strain WA-314. Additional acquisitions of strain WA-314 are seven prophage-like regions. One of these prophages, the 28-kb P4-like prophage YWA-4, encodes a PilV-like protein that may be used for adhesion to and invasion of the intestinal cells. Furthermore, a putative autotransporter and two type 1 fimbrial proteins of strain WA-314 show a sequence similarity <50% with the orthologous proteins in strain 8081. The dissimilar sequences of these proteins indicate possible different functions or interaction modes, reflecting the specific adhesion properties of Y. enterocolitica strains 8081 and WA-314 and thus the different efficiency of host colonization. Further important differences were found in two pYV plasmid-encoded virulence factors, YopM and YscP. The impact of these differences on virulence is discussed. CONCLUSIONS: Our study emphasizes that the virulence of pathogens can be increased, by acquiring new genes and/or improving the function of essential virulence proteins, resulting in permanently hyper-virulent strains. This work also highlights the importance of addressing genetic and phenotypic variations among closely related bacterial strains, even those belonging to the same bioserotype.     

Detection and characterisation of pr1 virulent gene deficiencies in the insect pathogenic fungus Metarhizium anisopliae

The method of suppressive subtractive hybridization was employed to map out genomic differences between the highly pathogenic Yersinia enterocolitica (Ye) biogroup 1B, serotype O:8 strain (WA-314) and the closely related apathogenic Y. enterocolitica biogroup 1A, serotype O:5 strain (NF-O). A novel IS10-like element, IS1330, uncovered by this technique was found to be uniquely present in high copy numbers among the highly pathogenic Y. enterocolitica 1B strains, while a single copy of the element was found in the low pathogenic Ye biogroup 4 serotype O:3 strain. The 1321-bp repetitive element has 19-bp imperfect inverted terminal repeats and is bracketed by a 10-bp duplication of the target sequence. The predicted transposase shares high homology with the IS10 open reading frame of the large virulence plasmid pWR501, of Shigella flexneri, with IS10 transposase of Salmonella typhi, and with IS1999 (tnpA) of Pseudomonas aeruginosa. The IS1330 tnp gene is transcribed in vitro and in vivo in HeLa cells. At least one copy of IS1330 flanks the recently described chromosomal type III secretion cluster in Y. enterocolitica WA-314, O:8, and future studies should shed light on whether this novel transposase mediates transposition events in highly pathogenic Y. enterocolitica strains, thus enhancing the genetic plasticity of this species.     

Characterization of insertion sequence IS892 and related elements from the cyanobacterium Anabaena sp. strain PCC 7120.

IS892, one of the several insertion sequence (IS) elements discovered in Anabaena sp. strain PCC 7120 (Y. Cai and C. P. Wolk, J. Bacteriol. 172:3138-3145, 1990), is 1,675 bp with 24-bp near-perfect inverted terminal repeats and has two open reading frames (ORFs) that could code for proteins of 233 and 137 amino acids. Upon insertion into target sites, this IS generates an 8-bp directly repeated target duplication. A 32-bp sequence in the region between ORF1 and ORF2 is similar to the sequence of the inverted termini. Similar inverted repeats are found within each of those three segments, and the sequences of these repeats bear some similarity to the 11-bp direct repeats flanking the 11-kb insertion interrupting the nifD gene of this strain (J. W. Golden, S. J. Robinson, and R. Haselkorn, Nature [London] 314:419-423, 1985). A sequence similar to that of a binding site for the Escherichia coli integration host factor is found about 120 bp from the left end of IS892. Partial nucleotide sequences of active IS elements IS892N and IS892T, members of the IS892 family from the same Anabaena strain, were shown to be very similar to the sequence of IS892.     

The high-pathogenicity island of Yersinia pseudotuberculosis can be inserted into any of the three chromosomal asn tRNA genes

Pathogenicity islands (PAIs) have been identified in several bacterial species. A PAI called high-pathogenicity island (HPI) and carrying genes involved in iron acquisition (yersiniabactin system) has been previously identified in Yersinia enterocolitica and Yersinia pestis . In this study, the HPI of the third species of Yersinia pathogenic for humans, Y. pseudotuberculosis , has been characterized. We demonstrate that the HPI of strain IP32637 has a physical and genetic map identical to that of Y. pestis . A gene homologous to the bacteriophage P4 integrase gene is located downstream of the asn tRNA locus that borders the HPI of strain IP32637. This int gene is at the same position on the HPI of all three pathogenic Yersinia species. However, in contrast to Y. pestis 6/69, the HPI of Y. pseudotuberculosis IP32637 is not invariably adjacent to the pigmentation segment and can be inserted at a distance ≥ 190 kb from this segment. Also, in contrast to Y. pestis and Y. enterocolitica , the HPI of Y. pseudotuberculosis IP32637 can precisely excise from the chromosome, and, strikingly, it can be found inserted in any of the three asn tRNA loci present on the chromosome of this species, one of which is adjacent to the pigmentation segment. The pigmentation segment, which is present in Y. pestis but not in Y. enterocolitica , is also present and well conserved in all strains of Y. pseudotuberculosis studied. In contrast, the presence and size of the HPIs vary depending on the serotype of the strain: an entire HPI is found in strains of serotypes I only, a HPI with a 9 kb truncation in its left-hand part that carries the IS 100 sequence and the psn and ybtE genes characterizes the strains of serotype III, and no HPI is found in strains of serotypes II, IV and V.     

Cloning and characterisation of the serC and aroA genes of Yersinia enterocolitica, and construction of an aroA mutant

A gene library of Yersinia enterocolitica 8081 was constructed in the cosmid vector pHC79. Recombinants containing the aroA gene, encoding 5-enolpyruvylshikimate 3-phosphate synthase, were identified by complementation of the aroA mutation in Escherichia coli K-12 strain AB2829. All six recombinant plasmids which complemented aroA also complemented the serC mutation in E. coli K-12 strain KL282. Tn5 mutagenesis suggested serC encoding 3-phosphoserine aminotransferase was the proximal gene in an operon with aroA. The nucleotide sequence of a 3-kb HindII-EcoRV fragment encoding the two genes was determined. The serC and aroA open reading frames contain 362 and 428 codons, respectively, and the deduced amino acid sequences share 78% and 81% homology, respectively, with the corresponding E. coli genes. Sequence inspection revealed no obvious terminators or promoters in the intergenic region. The cloned Y. enterocolitica aroA gene was inactivated in vitro and reintroduced into the parental Y. enterocolitica 8081 strain using the suicide vector pJM703.1. Stable aroA insertion mutants of Y. enterocolitica were isolated.     

Genetic analysis of the yopE region of Yersinia spp.: identification of a novel conserved locus, yerA, regulating yopE expression.

The yopE gene of Yersinia pseudotuberculosis was recently sequenced, and YopE was identified as an indispensable virulence determinant when tested in a mouse model (A. Forsberg and H. Wolf-Watz, Mol. Microbiol. 2:121-133, 1988). In the study described here, the DNA sequences of the yopE genes of Yersinia pestis EV76 and Yersinia enterocolitica 8081 were determined and compared with that of the Y. pseudotuberculosis gene. Only two codons were found to differ, both leading to amino acid replacements, when the gene from Y. pestis was compared. These two replacements were also present in the gene from Y. enterocolitica; in addition, 18 other codons were found to differ. Thirteen of these substitutions led to amino acid replacements. Downstream of the yopE gene, the plasmid partition locus par was found to be conserved in all three species. In Y. enterocolitica 8081, the sequence homology was interrupted by a putative insertion sequence element inserted between the yopE gene and the par region at a position only 5 base pairs downstream of the yopE stop codon. Upstream of the yopE gene, 620 base pairs were conserved in the three species. This region contained a 130-amino-acid-long open reading frame reading in the opposite direction to the yopE gene and expressed a 14-kilodalton protein in minicells. An insertion mutation in this region constructed in Y. pseudotuberculosis expressed significantly lower amounts of YopE protein in vitro than did the corresponding wild type. The expression level could be restored by transcomplementation. This new locus was designated yerA, for yopE-regulating gene A. The yerA mutant was avirulent when mice were challenged by oral infection.     

Characterization of the IS895 family of insertion sequences from the cyanobacterium Anabaena sp. strain PCC 7120.

A family of repetitive elements from the cyanobacterium Anabaena sp. strain PCC 7120 was identified through the proximity of one element to the psbAI gene. Four members of this seven-member family were isolated and shown to have structures characteristic of bacterial insertion sequences. Each element is approximately 1,200 bp in length, is delimited by a 30-bp inverted repeat, and contains two open reading frames in tandem on the same DNA strand. The four copies differ from each other by small insertions or deletions, some of which alter the open reading frames. By using a system designed to trap insertion elements, one of the elements, denoted IS895, was shown to be mobile. The target site was not duplicated upon insertion of the element. Two other filamentous cyanobacterial strains were also found to contain sequences homologous to IS895.     

Application of comparative phylogenomics to study the evolution of Yersinia enterocolitica and to identify genetic differences relating to pathogenicity

Yersinia enterocolitica, an important cause of human gastroenteritis generally caused by the consumption of livestock, has traditionally been categorized into three groups with respect to pathogenicity, i.e., nonpathogenic (biotype 1A), low pathogenicity (biotypes 2 to 5), and highly pathogenic (biotype 1B). However, genetic differences that explain variation in pathogenesis and whether different biotypes are associated with specific nonhuman hosts are largely unknown. In this study, we applied comparative phylogenomics (whole-genome comparisons of microbes with DNA microarrays combined with Bayesian phylogenies) to investigate a diverse collection of 94 strains of Y. enterocolitica consisting of 35 human, 35 pig, 15 sheep, and 9 cattle isolates from nonpathogenic, low-pathogenicity, and highly pathogenic biotypes. Analysis confirmed three distinct statistically supported clusters composed of a nonpathogenic clade, a low-pathogenicity clade, and a highly pathogenic clade. Genetic differences revealed 125 predicted coding sequences (CDSs) present in all highly pathogenic strains but absent from the other clades. These included several previously uncharacterized CDSs that may encode novel virulence determinants including a hemolysin, a metalloprotease, and a type III secretion effector protein. Additionally, 27 CDSs were identified which were present in all 47 low-pathogenicity strains and Y. enterocolitica 8081 but absent from all nonpathogenic 1A isolates. Analysis of the core gene set for Y. enterocolitica revealed that 20.8% of the genes were shared by all of the strains, confirming this species as highly heterogeneous, adding to the case for the existence of three subspecies of Y. enterocolitica. Further analysis revealed that Y. enterocolitica does not cluster according to source (host).     

Characterization of the ysa Pathogenicity Locus in the Chromosome of Yersinia enterocolitica and Phylogeny Analysis of Type III Secretion Systems

Several Gram negative bacteria use a complex system called "type III secretion system" (TTSS) to engage their host. The archetype of TTSS is the plasmid-encoded "Yop virulon" shared by the three species of pathogenic Yersinia (Y. pestis, Y. pseudotuberculosis, and Y. enterocolitica). A second TTSS, called Ysa (for Yersinia secretion apparatus) was recently described in Y. enterocolitica 8081, a strain from serotype O:8. In this study, we describe the ysa locus from A127/90, another strain of serotype O:8, and we extend the sequence to several new genes encoding Ysp proteins which are the substrates of this secretion system, and a putative chaperone SycB. According to the deduced protein sequences, the ysa system from A127/90 is identical to that of 8081. It is different from the chromosome-encoded TTSS of Y. pestis but is instead closely related to the Mxi-Spa TTSS of Shigella and to the SPI-1 encoded TTSS of Salmonella enterica. We further demonstrated that the ysa locus is only present in biotype IB strains of Y. enterocolitica. Including this new Ysa system, a phylogenetic analysis of the 26 known TTSSs was carried out, based on the sequence analysis of three conserved proteins. All the TTSSs fall into five different clusters. The phylogenetic tree of these TTSSs is completely different from the evolutionary tree based on 16S RNA, indicating that TTSSs have been distributed by horizontal transfer.     

Isolation of a new insertion element of Yersinia intermedia closely related to remnants of mobile genetic elements present on Yersinia plasmids harboring the Yop virulon

A new insertion element present in two alleles, designated IS1635.1 and IS1635.2, was identified on a plasmid of a Yersinia intermedia strain by hybridization with the Yersinia enterocolitica pYV virulence plasmid. IS1635.1 and IS1635.2 are 861 bp long, carry imperfect inverted terminal repeats and possess a single open reading frame encoding a putative transposase of the IS6 family. A truncated IS1635 element is present immediately downstream of element IS1635.2. The capacity of the IS1635 elements to mediate transposition in Yersinia was demonstrated with a R6K-derived suicide vector, where a kanamycin resistance gene had been inserted between IS1635.1 and IS1635.2. Hybridization and sequence alignments showed that remnants of IS1635-like insertion elements harboring large deletions and point mutations are present on the Yop virulon harboring plasmids of pathogenic Yersinia strains. In a few cases, the IS1635 element has also been found on plasmids of apathogenic Yersinia strains.     

Modulation of complement activityin Vitro andin Vivo byYersinia wild and mutant strains

The ability of released proteins (Yops) and surface lipopolysaccharides (LPS) from the wild-type strain Yersinia enterocolitica 8081-L2, serotype 0:8 to influence the complement activity was determined. Yops and LPS from wild-type and mutant strains showed different ability to affect the classical pathway (CP) functional complement activity in vitro. The serum CP activity was inhibited during the infection induced with six Y. enterocolitica and three Y. pseudotuberculosis strains in rabbits. The changed complement activity might be of importance for the course of Yersinia infections.     

Molecular characterization of beta-lactamase genes blaA and blaB of Yersinia enterocolitica biovar 1A

The beta-lactamase genes blaA and blaB were detected by PCR amplification in strains of Yersinia enterocolitica biovar 1A isolated from India, Germany, France and the USA. Both genes were detected in all strains. Polymerase chain reaction-restriction fragment length polymorphism revealed genetic heterogeneity in blaA but not in blaB. Cluster analysis of blaA restriction profiles grouped the strains into three groups. The blaA gene of Y. enterocolitica biovar 1A showed a high degree of sequence homology to that of Y. enterocolitica 8081 (biovar 1B) and Y. enterocolitica Y-56 (biovar 4), whereas homology was low with class A beta-lactamase genes of other members of the family Enterobacteriaceae. The pI 8.7 of enzyme Bla-A of Y. enterocolitica biovar 1A was similar to that of biovars 2, 3 and 4. The enzyme Bla-B focused at 6.8 and 7.1, indicating that biovar 1A strains produced a 'B-like' enzyme. This is the first study to have investigated the genetic heterogeneity of the beta-lactamase genes of Y. enterocolitica.     

Characterization of IS1541-like elements in Yersinia enterocolitica and Yersinia pseudotuberculosis

We characterized Yersinia enterocolitica and Yersinia pseudotuberculosis insertion sequences related to insertion sequence 1541, recently identified in Yersinia pestis. For each of the two species, two insertion sequence copies were cloned and sequenced. Genetic elements from Y. pseudotuberculosis were almost identical to insertion sequence 1541, whereas these from Y. enterocolitica were less related. Phylogenetic analysis of the putative transposases encoded by insertion sequences from the three pathogenic members of the genus Yersinia showed that they clustered with those encoded by Escherichia coli and Salmonella enterica elements belonging to the insertion sequence 200/insertion sequence 605 group. Insertion sequences originating from Y. pestis and Y. pseudotuberculosis constitute a monophyletic lineage distinct from that of Y. enterocolitica.     

Multilocus variable number tandem repeat analysis as a tool to discern genetic relationships among strains of Yersinia enterocolitica biovar 1A

Aims:  To identify variable number tandem repeat (VNTR)-containing loci, and to use multilocus VNTR (MLVA) to discern genetic relationships among strains of Yersinia enterocolitica biovar 1A isolated from diverse sources.
Methods and Results:  The whole genome sequence of Y. enterocolitica 8081 was analysed and eight VNTR loci with repeat sizes between 4 and 9 bp, and each containing more than four repeat copies were selected for MLVA typing of 88 strains of Y. enterocolitica . Of these, four loci were polymorphic and generated 26 MLVA genotypes among 81 strains of Y. enterocolitica biovar 1A. MLVA was found to be quite discriminatory (DI = 0·87). Cluster analysis and population modelling using minimum spanning tree (MST) clearly clustered Y. enterocolitica biovar 1A into two major groups.
Conclusions:  The MLVA is easy to perform and can be used to discern clonal relationship among strains of Y. enterocolitica . Also the phylogenetic relationships obtained with MLVA genotypes were in good agreement with those established by other typing methods.
Significance and Impact of the Study:  The MLVA method reported is relatively more discriminatory than the other genotyping methods and has the potential to be used as an epidemiological tool for the study of strains of Y. enterocolitica biovar 1A.     

Virulence-associated plasmids from Yersinia enterocolitica and Yersinia pestis.

A 44-megadalton plasmid associated with virulence and Ca2+ dependence from Yersinia enterocolitica 8081 was compared at the molecular level with a 47-megadalton plasmid associated with Ca2+ dependence from Yersinia pestis EV76. The plasmids were found to share 55% deoxyribonucleic acid sequence homology distributed over approximately 80% of the plasmid genomes. One region in which the plasmids differed was found to contain sequences concerned with essential plasmid functions. Forty-five mutants of Y. pestis were isolated which had spontaneously acquired the ability to grow on calcium-free medium (Ca2+ independence). Of these mutants, 21 were cured of their 47-megadalton plasmid, whereas the remaining had either suffered a major deletion of the plasmid or had a 2.2-kilobase insertion located in one of two adjacent BamHI restriction fragments encompassing approximately 9 kilobases. The inserted sequence was found at numerous sites on the Y. pestis chromosome and on all three plasmids in the strain and may represent a Y. pestis insertion sequence element.     

Mechanisms of spontaneous mutagenesis: an analysis of the spectrum of spontaneous mutation in the Escherichia coli lacI gene

We have obtained via DNA sequence analysis a spectrum of 174 spontaneous mutations occurring in the lac I gene of Escherichia coli. The spectrum comprised base substitution, frameshift, deletion, duplication and insertion mutations, of which the relative contributions to spontaneous mutation could be estimated. Two thirds of all lacI mutations occurred in the frameshift hotspot site. An analysis of the local DNA sequence suggested that the intensity of this hotspot may depend on structural features of the DNA that extend beyond those permitted by the repeated tetramer at this site. Deletions comprised the largest non-hotspot class (37%). They could be divided into two subclasses, depending on whether they included the lac operator sequence; the latter was found to be a preferred site for deletion endpoints. Most of the deletions internal to the lacI gene were associated with the presence of directly or invertedly repeated sequences capable of accounting for their endpoints. Base substitutions comprised 34% of the non-hotspot events. Unlike the base substitution spectrum obtained via nonsense mutations, G . C----A . T transitions do not predominate. A new base substitution hotspot was discovered at position +6 in the lac operator; its intensity may reflect specific features of the operator DNA. IS1 insertion mutations contributed 12% of the non-hotspot mutations and occurred dispersed throughout the gene in both orientations. Since the lacI gene is not A + T-rich, the contribution of IS1 insertion to spontaneous mutation in general might be underestimated. Single-base frameshift mutations were found only infrequently. In general, they did not occur in runs of a common base. Instead, their occurrence seemed based on the "perfection" of direct or inverted repeats in the local DNA sequence. Three (tandem) duplication events were recovered. No repeated sequences were found that might have determined their endpoints.     

Horizontal transfer of the high-pathogenicity island of Yersinia pseudotuberculosis

The horizontal transfer of genetic elements plays a major role in bacterial evolution. The high-pathogenicity island (HPI), which codes for an iron uptake system, is present and highly conserved in various Enterobacteriaceae, suggesting its recent acquisition by lateral gene transfer. The aim of this work was to determine whether the HPI has kept its ability to be transmitted horizontally. We demonstrate here that the HPI is indeed transferable from a donor to a recipient Yersinia pseudotuberculosis strain. This transfer was observable only when the donor and recipient bacteria were cocultured at low temperatures in a liquid medium. When optimized conditions were used (bacteria actively growing in an iron-deprived medium at 4 degrees C), the frequency of HPI transfer reached approximately 10(-8). The island was transferable to various serotype I strains of Y. pseudotuberculosis and to Yersinia pestis, but not to Y. pseudotuberculosis strains of serotypes II and IV or to Yersinia enterocolitica. Upon transfer, the HPI was inserted almost systematically into the asn3 tRNA locus. Acquisition of the HPI resulted in the loss of the resident island, suggesting an incompatibility between two copies of the HPI within the same strain. Transfer of the island did not require a functional HPI-borne insertion-excision machinery and was RecA dependent in the recipient but not the donor strain, suggesting that integration of the island into the recipient chromosome occurs via a mechanism of homologous recombination. This lateral transfer also involved the HPI-adjacent sequences, leading to the mobilization of a chromosomal region at least 46 kb in size.     

Detection and enumeration of virulent Yersinia enterocolitica in food by DNA colony hybridization.

A portion of a 44-megadalton plasmid found in Yersinia enterocolitica 8081 was used as a genetic probe to differentiate virulent and nonvirulent strains of the species. A DNA colony hybridization technique was employed. Three BamHI restriction endonuclease fragments labeled with 32P by nick translation were hybridized to lysed colonies of pure cultures, mixtures of virulent and nonvirulent cells, and portions of a food sample artificially contaminated with virulent Y. enterocolitica. The results of the colony hybridization test for virulence were the same as those obtained by the autoagglutination and suckling mouse tests. DNA colony hybridization detects pathogenic Y. enterocolitica in foods without the need for enrichment or pure cultures.     

Transposition of IS10 from the host Escherichia coli genome to a plasmid may lead to cloning artefacts

During recloning of Nicotiana tabacum L. repetitive sequence R8.3 in Escherichia coli, a modified clone that differed from the original by the insertion of an IS10 sequence was unintentionally produced. The insert was flanked by a 9-bp direct repeat derived from the R8.3 sequence, the 9-bp duplication of acceptor DNA in the site of insertion being a characteristic of IS10 transposition events. A database search using the FASTA program showed IS10 and other prokaryotic IS elements inserted into numerous eukaryotic clones. Unexpectedly, the IS10, which is not a natural component of the E. coli genome, appeared to be by far the most frequent contaminant of DNA databases among several IS sequences tested. In the GenEMBL database, the IS10 query sequence yielded positive scores with more than 500 eukaryotic clones. Insertions of shortened IS10 sequences having only one intact terminal inverted repeat were commonly found. Most full-length IS10 insertions (32 out of 40 analyzed) were flanked by 9-bp direct repeats having the consensus 5'-NPuCNN-NGPyN-3' with a strong preference for 5'-TGCTNA-GNN-3'. One insertion was flanked by an inverted repeat of more than 400 bp in length. PCR amplification and Southern analysis revealed the presence of IS10 sequences in E. coli strains commonly used for DNA cloning, including some reported to be Tn10-free. No IS10-specific PCR product was obtained with N. tabacum or human DNA. Our data suggest that transposition of IS10 elements may accompany cloning steps, particularly into large BAC vectors. This might lead to the relatively frequent contamination of DNA databases by this bacterial sequence. It is estimated that one in approximately every thousand eukaryotic clone in the databases is contaminated by IS-derived sequences. We recommend checking submitted sequences for the presence of IS10 and other IS elements. In addition, DNA databases should be corrected by removing contaminating IS sequences.     

IS1631 occurrence in Bradyrhizobium japonicum highly reiterated sequence-possessing strains with high copy numbers of repeated sequences RSalpha and RSbeta.

From Bradyrhizobium japonicum highly reiterated sequence-possessing (HRS) strains indigenous to Niigata and Tokachi in Japan with high copy numbers of the repeated sequences RSalpha and RSbeta (K. Minamisawa, T. Isawa, Y. Nakatsuka, and N. Ichikawa, Appl. Environ. Microbiol. 64:1845-1851, 1998), several insertion sequence (IS)-like elements were isolated by using the formation of DNA duplexes by denaturation and renaturation of total DNA, followed by treatment with S1 nuclease. Most of these sequences showed structural features of bacterial IS elements, terminal inverted repeats, and homology with known IS elements and transposase genes. HRS and non-HRS strains of B. japonicum differed markedly in the profiles obtained after hybridization with all the elements tested. In particular, HRS strains of B. japonicum contained many copies of IS1631, whereas non-HRS strains completely lacked this element. This association remained true even when many field isolates of B. japonicum were examined. Consequently, IS1631 occurrence was well correlated with B. japonicum HRS strains possessing high copy numbers of the repeated sequence RSalpha or RSbeta. DNA sequence analysis indicated that IS1631 is 2,712 bp long. In addition, IS1631 belongs to the IS21 family, as evidenced by its two open reading frames, which encode putative proteins homologous to IstA and IstB of IS21, and its terminal inverted repeat sequences with multiple short repeats.