Roles of selection and recombination in the evolution of type I restriction-modification systems in enterobacteria

Integrative DNA methylome and transcriptome analysis reveals DNA adenine methylation is involved in Salmonella enterica Typhimurium response to oxidative stress

IMPORTANCE

The intracellular pathogen Salmonella enterica serovar Typhimurium (S. Typhimurium) comes across a wide variety of stresses from entry to dissemination, such as reactive oxygen species. To adapt itself to oxidative stress, Salmonella must adopt various and complex strategies. In this study, we revealed that DNA adenine methyltransferase was essential for S. Typhimurium to survive in hydrogen peroxide. We then screened out oxidative stress-responsive genes that were potentially regulated by DNA methylation in S. Typhimurium. Our results show that the DNA methylome is highly stable throughout the genome, and the coupled change of m6A GATC with gene expression is identified in only a few positions, which suggests the complexity of the DNA methylation and gene expression regulation networks. The results may shed light on our understanding of m6A-mediated gene expression regulation in bacteria.

Salmonella Genomics in Public Health and Food Safety

The species Salmonella enterica comprises over 2,600 serovars, many of which are known to be intracellular pathogens of mammals, birds, and reptiles. It is now apparent that Salmonella is a highly adapted environmental microbe and can readily persist in a number of environmental niches, including water, soil, and various plant (including produce) species. Much of what is known about the evolution and diversity of nontyphoidal Salmonella serovars (NTS) in the environment is the result of the rise of the genomics era in enteric microbiology. There are over 340,000 Salmonella genomes available in public databases. This extraordinary breadth of genomic diversity now available for the species, coupled with widespread availability and affordability of whole-genome sequencing (WGS) instrumentation, has transformed the way in which we detect, differentiate, and characterize Salmonella enterica strains in a timely way. Not only have WGS data afforded a detailed and global examination of the molecular epidemiological movement of Salmonella from diverse environmental reservoirs into human and animal hosts, but they have also allowed considerable consolidation of the diagnostic effort required to test for various phenotypes important to the characterization of Salmonella. For example, drug resistance, serovar, virulence determinants, and other genome-based attributes can all be discerned using a genome sequence. Finally, genomic analysis, in conjunction with functional and phenotypic approaches, is beginning to provide new insights into the precise adaptive changes that permit persistence of NTS in so many diverse and challenging environmental niches.

Comparative genomics of Salmonella enterica serovar Enteritidis ST-11 isolated in Uruguay reveals lineages associated with particular epidemiological traits

Salmonella enterica serovar Enteritidis is a major cause of foodborne disease in Uruguay since 1995. We used a genomic approach to study a set of isolates from different sources and years. Whole genome phylogeny showed that most of the strains are distributed in two major lineages (E1 and E2), both belonging to MLST sequence type 11 the major ST among serovar Enteritidis. Strikingly, E2 isolates are over-represented in periods of outbreak abundance in Uruguay, while E1 span all epidemic periods. Both lineages circulate in neighbor countries at the same timescale as in Uruguay, and are present in minor numbers in distant countries. We identified allelic variants associated with each lineage. Three genes, ycdX, pduD and hsdM, have distinctive variants in E1 that may result in defective products. Another four genes (ybiO, yiaN, aas, aceA) present variants specific for the E2 lineage. Overall this work shows that S. enterica serovar Enteritidis strains circulating in Uruguay have the same phylogenetic profile than strains circulating in the region, as well as in more distant countries. Based on these results we hypothesize that the E2 lineage, which is more prevalent during epidemics, exhibits a combination of allelic variants that could be associated with its epidemic ability.

A comparative analysis of methylome profiles of Campylobacter jejuni sheep abortion isolate and gastroenteric strains using PacBio data

Campylobacter jejuni is a leading cause of human gastrointestinal disease and small ruminant abortions in the United States. The recent emergence of a highly virulent, tetracycline-resistant C. jejuni subsp. jejuni sheep abortion clone (clone SA) in the United States, and that strain's association with human disease, has resulted in a heightened awareness of the zoonotic potential of this organism. Pacific Biosciences' Single Molecule, Real-Time sequencing technology was used to explore the variation in the genome-wide methylation patterns of the abortifacient clone SA (IA3902) and phenotypically distinct gastrointestinal-specific C. jejuni strains (NCTC 11168 and 81-176). Several notable differences were discovered that distinguished the methylome of IA3902 from that of 11168 and 81-176: identification of motifs novel to IA3902, genome-specific hypo- and hypermethylated regions, strain level variability in genes methylated, and differences in the types of methylation motifs present in each strain. These observations suggest a possible role of methylation in the contrasting disease presentations of these three C. jejuni strains. In addition, the methylation profiles between IA3902 and a luxS mutant were explored to determine if variations in methylation patterns could be identified that might explain the role of LuxS-dependent methyl recycling in IA3902 abortifacient potential.

The other face of restriction: modification-dependent enzymes

The 1952 observation of host-induced non-hereditary variation in bacteriophages by Salvador Luria and Mary Human led to the discovery in the 1960s of modifying enzymes that glucosylate hydroxymethylcytosine in T-even phages and of genes encoding corresponding host activities that restrict non-glucosylated phage DNA: rglA and rglB (restricts glucoseless phage). In the 1980’s, appreciation of the biological scope of these activities was dramatically expanded with the demonstration that plant and animal DNA was also sensitive to restriction in cloning experiments. The rgl genes were renamed mcrA and mcrBC (modified cytosine restriction). The new class of modification-dependent restriction enzymes was named Type IV, as distinct from the familiar modification-blocked Types I–III. A third Escherichia coli enzyme, mrr (modified DNA rejection and restriction) recognizes both methylcytosine and methyladenine. In recent years, the universe of modification-dependent enzymes has expanded greatly. Technical advances allow use of Type IV enzymes to study epigenetic mechanisms in mammals and plants. Type IV enzymes recognize modified DNA with low sequence selectivity and have emerged many times independently during evolution. Here, we review biochemical and structural data on these proteins, the resurgent interest in Type IV enzymes as tools for epigenetic research and the evolutionary pressures on these systems.

Transposition of IS elements induced by electroporation of suicide plasmid in Acidithiobacillus caldus

Transposition insertional mutagenesis of the insertion sequences (IS elements) was discovered for the first time in Acidithiobacillus caldus (A. caldus), when A. caldus MTH-04 hsdM (type I restriction-modification system M-subunit) mutant was constructed by electroporation of a suicide plasmid. The IS element, specifically inserting into hsdM gene, was analyzed, identified, and named ISAtc2. The transposition frequency of ISAtc2 was ranged from 4% to 7%, and no reverse mutation occurred in the mutants after 50 generations of proliferation without selective pressure. These results revealed that transposition of IS elements on A. caldus chromosome could regulate the gene expression and metabolic pathways by gene inactivation, gene loss and gene acquisition. Therefore, the transposition of IS elements in A. caldus may be an important and unique regulation mechanism for adaptation to the living condition.

Movement of DNA sequence recognition domains between non-orthologous proteins

Comparisons of proteins show that they evolve through the movement of domains. However, in many cases, the underlying mechanisms remain unclear. Here, we observed the movements of DNA recognition domains between non-orthologous proteins within a prokaryote genome. Restriction-modification (RM) systems, consisting of a sequence-specific DNA methyltransferase and a restriction enzyme, contribute to maintenance/evolution of genomes/epigenomes. RM systems limit horizontal gene transfer but are themselves mobile. We compared Type III RM systems in Helicobacter pylori genomes and found that target recognition domain (TRD) sequences are mobile, moving between different orthologous groups that occupy unique chromosomal locations. Sequence comparisons suggested that a likely underlying mechanism is movement through homologous recombination of similar DNA sequences that encode amino acid sequence motifs that are conserved among Type III DNA methyltransferases. Consistent with this movement, incongruence was observed between the phylogenetic trees of TRD regions and other regions in proteins. Horizontal acquisition of diverse TRD sequences was suggested by detection of homologs in other Helicobacter species and distantly related bacterial species. One of these RM systems in H. pylori was inactivated by insertion of another RM system that likely transferred from an oral bacterium. TRD movement represents a novel route for diversification of DNA-interacting proteins.

The phage-host arms race: shaping the evolution of microbes

Bacteria, the most abundant organisms on the planet, are outnumbered by a factor of 10 to 1 by phages that infect them. Faced with the rapid evolution and turnover of phage particles, bacteria have evolved various mechanisms to evade phage infection and killing, leading to an evolutionary arms race. The extensive co-evolution of both phage and host has resulted in considerable diversity on the part of both bacterial and phage defensive and offensive strategies. Here, we discuss the unique and common features of phage resistance mechanisms and their role in global biodiversity. The commonalities between defense mechanisms suggest avenues for the discovery of novel forms of these mechanisms based on their evolutionary traits.

Characterization of a restriction modification system from the commensal Escherichia coli strain A0 34/86 (O83:K24:H31)

Background

Type I restriction-modification (R-M) systems are the most complex restriction enzymes discovered to date. Recent years have witnessed a renaissance of interest in R-M enzymes Type I. The massive ongoing sequencing programmes leading to discovery of, so far, more than 1 000 putative enzymes in a broad range of microorganisms including pathogenic bacteria, revealed that these enzymes are widely represented in nature. The aim of this study was characterisation of a putative R-M system EcoA0ORF42P identified in the commensal Escherichia coli A0 34/86 (O83: K24: H31) strain, which is efficiently used at Czech paediatric clinics for prophylaxis and treatment of nosocomial infections and diarrhoea of preterm and newborn infants.

Results

We have characterised a restriction-modification system EcoA0ORF42P of the commensal Escherichia coli strain A0 34/86 (O83: K24: H31). This system, designated as EcoAO83I, is a new functional member of the Type IB family, whose specificity differs from those of known Type IB enzymes, as was demonstrated by an immunological cross-reactivity and a complementation assay. Using the plasmid transformation method and the RM search computer program, we identified the DNA recognition sequence of the EcoAO83I as GGA(8N)ATGC. In consistence with the amino acids alignment data, the 3' TRD component of the recognition sequence is identical to the sequence recognized by the EcoEI enzyme. The A-T (modified adenine) distance is identical to that in the EcoAI and EcoEI recognition sites, which also indicates that this system is a Type IB member. Interestingly, the recognition sequence we determined here is identical to the previously reported prototype sequence for Eco377I and its isoschizomers.

Conclusion

Putative restriction-modification system EcoA0ORF42P in the commensal Escherichia coli strain A0 34/86 (O83: K24: H31) was found to be a member of the Type IB family and was designated as EcoAO83I. Combination of the classical biochemical and bacterial genetics approaches with comparative genomics might contribute effectively to further classification of many other putative Type-I enzymes, especially in clinical samples.

High allelic diversity in the methyltransferase gene of a phase variable type III restriction-modification system has implications for the fitness of Haemophilus influenzae

Phase variable restriction-modification (R-M) systems are widespread in Eubacteria. Haemophilus influenzae encodes a phase variable homolog of Type III R-M systems. Sequence analysis of this system in 22 non-typeable H.influenzae isolates revealed a hypervariable region in the central portion of the mod gene whereas the res gene was conserved. Maximum likelihood (ML) analysis indicated that most sites outside this hypervariable region experienced strong negative selection but evidence of positive selection for a few sites in adjacent regions. A phylogenetic analysis of 61 Type III mod genes revealed clustering of these H.influenzae mod alleles with mod genes from pathogenic Neisseriae and, based on sequence analysis, horizontal transfer of the mod-res complex between these species. Neisserial mod alleles also contained a hypervariable region and all mod alleles exhibited variability in the repeat tract. We propose that this hypervariable region encodes the target recognition domain (TRD) of the Mod protein and that variability results in alterations to the recognition sequence of this R-M system. We argue that the high allelic diversity and phase variable nature of this R-M system have arisen due to selective pressures exerted by diversity in bacteriophage populations but also have implications for other fitness attributes of these bacterial species.

Crystal structure of a putative type I restriction-modification S subunit from Mycoplasma genitalium

The crystal structure of the eubacteria Mycoplasma genitalium ORF MG438 polypeptide, determined by multiple anomalous dispersion and refined at 2.3 A resolution, reveals the organization of S subunits from the Type I restriction and modification system. The structure consists of two globular domains, with about 150 residues each, separated by a pair of 40 residue long antiparallel alpha-helices. The globular domains correspond to the variable target recognition domains (TRDs), as previously defined for S subunits on sequence analysis, while the two helices correspond to the central (CR1) and C-terminal (CR2) conserved regions, respectively. The structure of the MG438 subunit presents an overall cyclic topology with an intramolecular 2-fold axis that superimposes the N and the C-half parts, each half containing a globular domain and a conserved helix. TRDs are found to be structurally related with the small domain of the Type II N6-adenine DNA MTase TaqI. These relationships together with the structural peculiarities of MG438, in particular the presence of the intramolecular quasi-symmetry, allow the proposal of a model for S subunits recognition of their DNA targets in agreement with previous experimental results. In the crystal, two subunits of MG438 related by a crystallographic 2-fold axis present a large contact area mainly involving the symmetric interactions of a cluster of exposed hydrophobic residues. Comparison with the recently reported structure of an S subunit from the archaea Methanococcus jannaschii highlights the structural features preserved despite a sequence identity below 20%, but also reveals important differences in the globular domains and in their disposition with respect to the conserved regions.

Comparative analyses of Xanthomonas and Xylella complete genomes

Computational analyses of four bacterial genomes of the Xanthomonadaceae family reveal new unique genes that may be involved in adaptation, pathogenicity, and host specificity. The Xanthomonas genus presents 3636 unique genes distributed in 1470 families, while Xylella genus presents 1026 unique genes distributed in 375 families. Among Xanthomonas-specific genes, we highlight a large number of cell wall degrading enzymes, proteases, and iron receptors, a set of energy metabolism genes, second copy of the type II secretion system, type III secretion system, flagella and chemotactic machinery, and the xanthomonadin synthesis gene cluster. Important genes unique to the Xylella genus are an additional copy of a type IV pili gene cluster and the complete machinery of colicin V synthesis and secretion. Intersections of gene sets from both genera reveal a cluster of genes homologous to Salmonella's SPI-7 island in Xanthomonas axonopodis pv citri and Xylella fastidiosa 9a5c, which might be involved in host specificity. Each genome also presents important unique genes, such as an HMS cluster, the kdgT gene, and O-antigen in Xanthomonas axonopodis pv citri; a number of avrBS genes and a distinct O-antigen in Xanthomonas campestris pv campestris, a type I restriction-modification system and a nickase gene in Xylella fastidiosa 9a5c, and a type II restriction-modification system and two genes related to peptidoglycan biosynthesis in Xylella fastidiosa temecula 1. All these differences imply a considerable number of gene gains and losses during the divergence of the four lineages, and are associated with structural genome modifications that may have a direct relation with the mode of transmission, adaptation to specific environments and pathogenicity of each organism.

Diversity within the Campylobacter jejuni type I restriction-modification loci

The type I restriction-modification (hsd) systems of 73 Campylobacter jejuni strains were characterized according to their DNA and amino acid sequences, and/or gene organization. A number of new genes were identified which are not present in the sequenced strain NCTC 11168. The closely related organism Helicobacter pylori has three type I systems; however, no evidence was found that C. jejuni strains contain multiple type I systems, although hsd loci are present in at least two different chromosomal locations. Also, unlike H. pylori, intervening ORFs are present, in some strains, between hsdR and hsdS and between hsdS and hsdM. No definitive function can be ascribed to these ORFs, designated here as rloA-H (R-linked ORF) and mloA-B (M-linked ORF). Based on parsimony analysis of amino acid sequences to assess character relatedness, the C. jejuni type I R-M systems are assigned to one of three families: 'IAB', 'IC' or 'IF'. This study confirms that HsdM proteins within a family are highly conserved but share little homology with HsdM proteins from other families. The 'IC' hsd loci are >99 % identical at the nucleotide level, as are the 'IF' hsd loci. Additionally, whereas the nucleotide sequences of the 'IAB' hsdR and hsdM genes show a high degree of similarity, the nucleotide sequences of the 'IAB' hsdS and rlo genes vary considerably. This diversity suggests that recombination between 'IAB' hsd loci would lead not only to new hsdS alleles but also to the exchange of rlo genes; five C. jejuni hsd loci are presumably the result of such recombination. The importance of these findings with regard to the evolution of C. jejuni type I R-M systems is discussed.

Functional characterization of a catabolic plasmid from polychlorinated- biphenyl-degrading Rhodococcus sp. strain RHA1

Rhodococcus sp. strain RHA1, a potent polychlorinated-biphenyl (PCB)-degrading strain, contains three linear plasmids ranging in size from 330 to 1,100 kb. As part of a genome sequencing project, we report here the complete sequence and characterization of the smallest and least-well-characterized of the RHA1 plasmids, pRHL3. The plasmid is an actinomycete invertron, containing large terminal inverted repeats with a tightly associated protein and a predicted open reading frame (ORF) that is similar to that of a mycobacterial rep gene. The pRHL3 plasmid has 300 putative genes, almost 21% of which are predicted to have a catabolic function. Most of these are organized into three clusters. One of the catabolic clusters was predicted to include limonene degradation genes. Consistent with this prediction, RHA1 grew on limonene, carveol, or carvone as the sole carbon source. The plasmid carries three cytochrome P450-encoding (CYP) genes, a finding consistent with the high number of CYP genes found in other actinomycetes. Two of the CYP genes appear to belong to novel families; the third belongs to CYP family 116 but appears to belong to a novel class based on the predicted domain structure of its reductase. Analyses indicate that pRHL3 also contains four putative "genomic islands" (likely to have been acquired by horizontal transfer), insertion sequence elements, 19 transposase genes, and a duplication that spans two ORFs. One of the genomic islands appears to encode resistance to heavy metals. The plasmid does not appear to contain any housekeeping genes. However, each of the three catabolic clusters contains related genes that appear to be involved in glucose metabolism.

Cellular localization of Type I restriction-modification enzymes is family dependent

Cellular localization of Type I restriction-modification enzymes EcoKI, EcoAI, and EcoR124I-the most frequently studied representatives of IA, IB, and IC families-was analyzed by immunoblotting of subcellular fractions isolated from Escherichia coli strains harboring the corresponding hsd genes. EcoR124I shows characteristics similar to those of EcoKI. The complex enzymes are associated with the cytoplasmic membrane via DNA interaction as documented by the release of the Hsd subunits from the membrane into the soluble fraction following benzonase treatment. HsdR subunits of the membrane-bound enzymes EcoKI and EcoR124I are accessible, though to a different extent, at the external surface of cytoplasmic membrane as shown by trypsinization of intact spheroplasts. EcoAI strongly differs from EcoKI and EcoR124I, since neither benzonase nor trypsin affects its association with the cytoplasmic membrane. Possible reasons for such a different organization are discussed in relation of the control of the restriction-modification activities in vivo.

Tetra-amino-acid tandem repeats are involved in HsdS complementation in type IC restriction-modification systems

All known type I restriction and modification (R-M) systems of Escherichia coli and Salmonella enterica belong to one of four discrete families: type IA, IB, IC or ID. The classification of type I systems from a wide range of other genera is mainly based on complementation and molecular evidence derived from the comparison of the amino acid similarity of the corresponding subunits. This affiliation was seldom based on the strictest requirement for membership of a family, which depends on relatedness as demonstrated by complementation tests. This paper presents data indicating that the type I NgoAV R-M system from Neisseria gonorrhoeae, despite the very high identity of HsdM and HsdR subunits with members of the type IC family, does not show complementation with E. coli type IC R-M systems. Sequence analysis of the HsdS subunit of several different potential type IC R-M systems shows that the presence of different tetra-amino-acid sequence repeats, e.g. TAEL, LEAT, SEAL, TSEL, is characteristic for type IC R-M systems encoded by distantly related bacteria. The other regions of the HsdS subunits potentially responsible for subunit interaction are also different between a group of distantly related bacteria, but show high similarity within these bacteria. Complementation between the NgoAV R-M system and members of the EcoR124 R-M family can be restored by changing the tetra-amino-acid repeat within the HsdS subunit. The authors propose that the type IC family of R-M systems could consist of several complementation subgroups whose specificity would depend on differences in the conserved regions of the HsdS polypeptide.

Cassette-like variation of restriction enzyme genes in Escherichia coli C and relatives

A surprising result of comparative bacterial genomics has been the large amount of DNA found to be present in one strain but not in another of the same species. We examine in detail one location where gene content varies extensively, the restriction cluster in Escherichia coli. This region is designated the Immigration Control Region (ICR) for the density and variability of restriction functions found there. To better define the boundaries of this variable locus, we determined the sequence of the region from a restrictionless strain, E.coli C. Here we compare the 13.7 kb E.coli C sequence spanning the site of the ICR with corresponding sequences from five E.coli strains and Salmonella typhimurium LT2. To discuss this variation, we adopt the term 'framework' to refer to genes that are stable components of genomes within related lineages, while 'migratory' genes are transient inhabitants of the genome. Strikingly, seven different migratory DNA segments, encoding different sets of genes and gene fragments, alternatively occupy a single well-defined location in the seven strains examined. The flanking framework genes, yjiS and yjiA, display approximately normal patterns of conservation. The patterns observed are consistent with the action of a site-specific recombinase. Since no nearby gene codes for a likely recombinase of known families, such a recombinase must be of a new family or unlinked.

Tracking EcoKI and DNA fifty years on: a golden story full of surprises

1953 was a historical year for biology, as it marked the birth of the DNA helix, but also a report by Bertani and Weigle on 'a barrier to infection' of bacteriophage lambda in its natural host, Escherichia coli K-12, that could be lifted by 'host-controlled variation' of the virus. This paper lay dormant till Nobel laureate Arber and PhD student Dussoix showed that the lambda DNA was rejected and degraded upon infection of different bacterial hosts, unless it carried host-specific modification of that DNA, thus laying the foundations for the phenomenon of restriction and modification (R-M). The restriction enzyme of E.coli K-12, EcoKI, was purified in 1968 and required S-adenosylmethionine (AdoMet) and ATP as cofactors. By the end of the decade there was substantial evidence for a chromosomal locus hsdK with three genes encoding restriction (R), modification (M) and specificity (S) subunits that assembled into a large complex of >400 kDa. The 1970s brought the message that EcoKI cut away from its DNA recognition target, to which site the enzyme remained bound while translocating the DNA past itself, with concomitant ATP hydrolysis and subsequent double-strand nicks. This translocation event created clearly visible DNA loops in the electron microscope. EcoKI became the archetypal Type I R-M enzyme with curious DNA translocating properties reminiscent of helicases, recognizing the bipartite asymmetric site AAC(N6)GTGC. Cloning of the hsdK locus in 1976 facilitated molecular understanding of this sophisticated R-M complex and in an elegant 'pas de deux' Murray and Dryden constructed the present model based on a large body of experimental data plus bioinformatics. This review celebrates the golden anniversary of EcoKI and ends with the exciting progress on the vital issue of restriction alleviation after DNA damage, also first reported in 1953, which involves intricate control of R subunit activity by the bacterial proteasome ClpXP, important results that will keep scientists on the EcoKI track for another 50 years to come.

Complex restriction enzymes: NTP-driven molecular motors

Survival is assuredly the prime directive for all living organisms either as individuals or as a species. One of the main challenges encountered by bacterial populations is the danger of bacteriophage attacks, since infection of a single bacterium may rapidly propagate, decimating the entire population. In order to protect themselves against this acute threat, bacteria have developed an array of defence mechanisms, which range from preventing the infection itself via interference with bacteriophage adsorption to the cell surface and prevention of phage DNA injection, to degradation of the injected phage DNA. This last defence mechanism is catalysed by the bacterial restriction-modification (R-M) systems, and in particular, by nucleoside 5'-triphosphate (NTP)-dependent restriction enzymes, e.g. type I and type III R-M systems or the modification-dependent endonucleases. Type I and type III restriction systems have dual properties. They may either act as methylases and protect the host's own DNA against restriction by methylating specific residues, or they catalyse ATP-dependent endonuclease activity so that invading foreign DNA lacking the host-specific methylation is degraded. These defence mechanism systems are further complemented by the presence of methylation-dependent, GTP-dependent endonucleases, that restricts specifically methylated DNA. Although all three types of endonucleases are structurally very different, they share a common functional mechanism. They recognise and bind to specific DNA sequences but do not cleave DNA within those target sites. They belong to the general class of DNA motor proteins, which use the free energy associated with nucleoside 5'-triphosphate hydrolysis to translocate DNA so that the subsequent DNA cleavage event occurs at a distance from the endonuclease recognition site. Moreover, DNA cleavage appears to be a random process triggered upon stalling of the DNA translocation process and requiring dimerisation of the bound endonucleases for a concerted break of both DNA strands. In this review, we present a detailed description and analysis of the functional mechanism of the three known NTP-dependent restriction systems: type I and type III restriction-modification enzymes, as well as the methylation-dependent McrBC endonuclease.

Families of restriction enzymes: an analysis prompted by molecular and genetic data for type ID restriction and modification systems

Current genetic and molecular evidence places all the known type I restriction and modification systems of Escherichia coli and Salmonella enterica into one of four discrete families: type IA, IB, IC or ID. StySBLI is the founder member of the ID family. Similarities of coding sequences have identified restriction systems in E.coli and Klebsiella pneumoniae as probable members of the type ID family. We present complementation tests that confirm the allocation of EcoR9I and KpnAI to the ID family. An alignment of the amino acid sequences of the HsdS subunits of StySBLI and EcoR9I identify two variable regions, each predicted to be a target recognition domain (TRD). Consistent with two TRDs, StySBLI was shown to recognise a bipartite target sequence, but one in which the adenine residues that are the substrates for methylation are separated by only 6 bp. Implications of family relationships are discussed and evidence is presented that extends the family affiliations identified in enteric bacteria to a wide range of other genera.

Nucleoside triphosphate-dependent restriction enzymes

The known nucleoside triphosphate-dependent restriction enzymes are hetero-oligomeric proteins that behave as molecular machines in response to their target sequences. They translocate DNA in a process dependent on the hydrolysis of a nucleoside triphosphate. For the ATP-dependent type I and type III restriction and modification systems, the collision of translocating complexes triggers hydrolysis of phosphodiester bonds in unmodified DNA to generate double-strand breaks. Type I endonucleases break the DNA at unspecified sequences remote from the target sequence, type III endonucleases at a fixed position close to the target sequence. Type I and type III restriction and modification (R-M) systems are notable for effective post-translational control of their endonuclease activity. For some type I enzymes, this control is mediated by proteolytic degradation of that subunit of the complex which is essential for DNA translocation and breakage. This control, lacking in the well-studied type II R-M systems, provides extraordinarily effective protection of resident DNA should it acquire unmodified target sequences. The only well-documented GTP-dependent restriction enzyme, McrBC, requires methylated target sequences for the initiation of phosphodiester bond cleavage.

Genome sequence of enterohaemorrhagic Escherichia coli O157:H7

The bacterium Escherichia coli O157:H7 is a worldwide threat to public health and has been implicated in many outbreaks of haemorrhagic colitis, some of which included fatalities caused by haemolytic uraemic syndrome. Close to 75,000 cases of O157:H7 infection are now estimated to occur annually in the United States. The severity of disease, the lack of effective treatment and the potential for large-scale outbreaks from contaminated food supplies have propelled intensive research on the pathogenesis and detection of E. coli O157:H7 (ref. 4). Here we have sequenced the genome of E. coli O157:H7 to identify candidate genes responsible for pathogenesis, to develop better methods of strain detection and to advance our understanding of the evolution of E. coli, through comparison with the genome of the non-pathogenic laboratory strain E. coli K-12 (ref. 5). We find that lateral gene transfer is far more extensive than previously anticipated. In fact, 1,387 new genes encoded in strain-specific clusters of diverse sizes were found in O157:H7. These include candidate virulence factors, alternative metabolic capacities, several prophages and other new functions--all of which could be targets for surveillance.

Type I restriction systems: sophisticated molecular machines (a legacy of Bertani and Weigle)

Restriction enzymes are well known as reagents widely used by molecular biologists for genetic manipulation and analysis, but these reagents represent only one class (type II) of a wider range of enzymes that recognize specific nucleotide sequences in DNA molecules and detect the provenance of the DNA on the basis of specific modifications to their target sequence. Type I restriction and modification (R-M) systems are complex; a single multifunctional enzyme can respond to the modification state of its target sequence with the alternative activities of modification or restriction. In the absence of DNA modification, a type I R-M enzyme behaves like a molecular motor, translocating vast stretches of DNA towards itself before eventually breaking the DNA molecule. These sophisticated enzymes are the focus of this review, which will emphasize those aspects that give insights into more general problems of molecular and microbial biology. Current molecular experiments explore target recognition, intramolecular communication, and enzyme activities, including DNA translocation. Type I R-M systems are notable for their ability to evolve new specificities, even in laboratory cultures. This observation raises the important question of how bacteria protect their chromosomes from destruction by newly acquired restriction specifities. Recent experiments demonstrate proteolytic mechanisms by which cells avoid DNA breakage by a type I R-M system whenever their chromosomal DNA acquires unmodified target sequences. Finally, the review will reflect the present impact of genomic sequences on a field that has previously derived information almost exclusively from the analysis of bacteria commonly studied in the laboratory.

Molecular evolution of the Escherichia coli chromosome. V. Recombination patterns among strains of diverse origin

Incorporation patterns of donor DNA into recipient chromosomes following transduction or conjugation have been studied in the progeny of a variety of Escherichia coli crosses in which donor and recipient nucleotide sequences differ by 1-3%. Series of contiguous or variously spaced PCR fragments have been amplified from each recombinant chromosome and digested with a commercial restriction endonuclease previously shown to distinguish the respective parents in a given fragment. We conclude that entering donor DNA fragments are frequently abridged (cut and shortened) before incorporation, the cutting being due to restriction systems, and the shortening presumably due to exonuclease activity. Analysis of several backcrosses confirms, and extends to conjugation, the importance of restriction in E. coli recombination in nature. The transmission patterns in conjugation are similar to those of transduction, but (as expected) on a much larger scale. Asymmetric results of reciprocal crosses imply that mismatch frequency is not a major factor. Marked differences among the results of simple crosses according to parental strain combinations are consistent with observations that E. coli strains in nature vary dramatically in their restriction-modification systems.

Combinational variation of restriction modification specificities in Lactococcus lactis

Three genes coding for a type I R-M system related to the class C enzymes have been identified on the chromosome of Lactococcus lactis strain IL1403. In addition, plasmids were found that encode only the HsdS subunit that directs R-M specificity. The presence of these plasmids in IL1403 conferred a new R-M phenotype on the host, indicating that the plasmid-encoded HsdS is able to interact with the chromosomally encoded HsdR and HsdM subunits. Such combinational variation of type I R-M systems may facilitate the evolution of their specificity and thus reinforce bacterial resistance against invasive foreign unmethylated DNA.

A type IC restriction-modification system in Lactococcus lactis

Three genes coding for the endonuclease, methylase, and specificity subunits of a type I restriction-modification (R-M) system in the Lactococcus lactis plasmid pIL2614 have been characterized. Plasmid location, sequence homologies, and inactivation studies indicated that this R-M system is most probably of type IC.

The restriction-modification genes of Escherichia coli K-12 may not be selfish: they do not resist loss and are readily replaced by alleles conferring different specificities

Type II restriction and modification (R-M) genes have been described as selfish because they have been shown to impose selection for the maintenance of the plasmid that encodes them. In our experiments, the type I R-M system EcoKI does not behave in the same way. The genes specifying EcoKI are, however, normally residents of the chromosome and therefore our analyses were extended to monitor the deletion of chromosomal genes rather than loss of plasmid vector. If EcoKI were to behave in the same way as the plasmid-encoded type II R-M systems, the loss of the relevant chromosomal genes by mutation or recombination should lead to cell death because the cell would become deficient in modification enzyme and the bacterial chromosome would be vulnerable to the restriction endonuclease. Our data contradict this prediction; they reveal that functional type I R-M genes in the chromosome are readily replaced by mutant alleles and by alleles encoding a type I R-M system of different specificity. The acquisition of allelic genes conferring a new sequence specificity, but not the loss of the resident genes, is dependent on the product of an unlinked gene, one predicted [Prakash-Cheng, A., Chung, S. S. & Ryu, J. (1993) Mol. Gen. Genet. 241, 491-496] to be relevant to control of expression of the genes that encode EcoKI. Our evidence suggests that not all R-M systems are evolving as "selfish" units; rather, the diversity and distribution of the family of type I enzymes we have investigated require an alternative selective pressure.

A prediction of the amino acids and structures involved in DNA recognition by type I DNA restriction and modification enzymes

The S subunits of type I DNA restriction/modification enzymes are responsible for recognising the DNA target sequence for the enzyme. They contain two domains of approximately 150 amino acids, each of which is responsible for recognising one half of the bipartite asymmetric target. In the absence of any known tertiary structure for type I enzymes or recognisable DNA recognition motifs in the highly variable amino acid sequences of the S subunits, it has previously not been possible to predict which amino acids are responsible for sequence recognition. Using a combination of sequence alignment and secondary structure prediction methods to analyse the sequences of S subunits, we predict that all of the 51 known target recognition domains (TRDs) have the same tertiary structure. Furthermore, this structure is similar to the structure of the TRD of the C5-cytosine methyltransferase, Hha I, which recognises its DNA target via interactions with two short polypeptide loops and a beta strand. Our results predict the location of these sequence recognition structures within the TRDs of all type I S subunits.

KpnAI, a new type I restriction-modification system in Klebsiella pneumoniae

The KpnAI restriction-modification (R-M) system has been identified in Klebsiella pneumoniae strain M5a1. The restriction gene of KpnAI was first cloned into pBR322 using an r-m+ M5a1 derivative and phage SBS for screening. Subsequently, an adjacent DNA fragment showing modification activity was cloned into pUC19. A total of 7.2 kb DNA sequencing data revealed three open reading frames, corresponding to hsdR, hsdM and hsdS genes of type I R-M systems. The predicted hsdR, hsdM and hsdS-coded peptides shared 95%, 98% and 44% identity, respectively, with the corresponding peptides of the recently identified StySBLI system, a prototype of the type ID family. This high homology suggests that KpnAI is also a member of the type ID family. The KpnAI system seems to be the first type I system identified in Klebsiella species.

Avoidance of palindromic words in bacterial and archaeal genomes: a close connection with restriction enzymes

Short palindromic sequences (4, 5 and 6 bp palindromes) are avoided at a statistically significant level in the genomes of several bacteria, including the completely sequenced Haemophilus influenzae and Synechocystis sp. genomes and in the complete genome of the archaeon Methanococcus jannaschii. In contrast, there is only moderate avoidance of palindromes in the small genome of the bacterium Mycoplasma genitalium and no detectable avoidance in the genomes of chloroplasts and mitochondria. The sites for type II restriction-modification enzymes detected in the given species tend to be among the most avoided palindromes in a particular genome, indicating a direct connection between the avoidance of short oligonucleotide words and restriction-modification systems with the respective specificity. Palindromes corresponding to sites for restriction enzymes from other species are also avoided, albeit less significantly, suggesting that in the course of evolution bacterial DNA has been exposed to a wide spectrum of restriction enzymes, probably as the result of lateral transfer mediated by mobile genetic elements, such as plasmids and prophages. Palindromic words appear to accumulate in DNA once it becomes isolated from restriction-modification systems, as demonstrated by the case of organellar genomes. By combining these observations with protein sequence analysis, we show that the most avoided 4-palindrome and the most avoided 6-palindrome in the archaeon M.jannaschii are likely to be recognition sites for two novel restriction-modification systems.

The specificity of sty SKI, a type I restriction enzyme, implies a structure with rotational symmetry

The type I restriction and modification (R-M) enzyme from Salmonella enterica serovar kaduna ( Sty SKI) recognises the DNA sequence 5'-CGAT(N)7GTTA, an unusual target for a type I R-M system in that it comprises two tetranucleotide components. The amino target recognition domain (TRD) of Sty SKI recognises 5'-CGAT and shows 36% amino acid identity with the carboxy TRD of Eco R124I which recognises the complementary, but degenerate, sequence 5'-RTCG. Current models predict that the amino and carboxy TRDs of the specificity subunit are in inverted orientations within a structure with 2-fold rotational symmetry. The complementary target sequences recognised by the amino TRD of Sty SKI and the carboxy TRD of Eco R124I are consistent with the predicted inverted positions of the TRDs. Amino TRDs of similar amino acid sequence have been shown to recognise the same nucleotide sequence. The similarity reported here, the first example of one between amino and carboxy TRDs, while consistent with a conserved mechanism of target recognition, offers additional flexibility in the evolution of sequence specificity by increasing the potential diversity of DNA targets for a given number of TRDs. Sty SKI identifies the first member of the IB family in Salmonella species.

Regulation of the activity of the type IC EcoR124I restriction enzyme

Restriction-modification (R-M) systems must regulate the expression of their genes so that the chromosomal genome is modified at all times by the methyltransferase to protect the host cell from the potential lethal action of the cognate restriction endonuclease. Since type I R-M systems can be transferred to non-modified Escherichia coli cells by conjugation or transformation without killing the recipient, they must have some means to regulate their restriction activity upon entering a new host cell to avoid restriction of unprotected host DNA and cell death. This is especially true for EcoR124I, a type IC family member, which is coded for by a conjugative plasmid. Control of EcoR124I restriction activity is most likely at the post-translational level as the transfer of the EcoR124I system into a recipient cell that already expressed the HsdR subunit of this system was not a lethal event. Additionally, the kinetics of restriction activity upon transfer of the genes coding for the EcoR124I RM system to a recipient cell are the same, irrespective of the modification state of the recipient cell or the presence or absence of the EcoR124I HsdR subunit in the new host cells. The mechanism controlling the restriction activity of a type IC R-M system upon transfer to a new host cell is different from that controlling the chromosomally coded type IA and IB R-M systems. The previously discovered hsdC mutant, which affects the establishment of the type IA system EcoKI, was shown to affect the establishment of the type IB system EcoAI, but to have no influence on EcoR124I.

Restriction by EcoKI is enhanced by co-operative interactions between target sequences and is dependent on DEAD box motifs

One subunit of both type I and type III restriction and modification enzymes contains motifs characteristic of DEAD box proteins, which implies that these enzymes may be DNA helicases. This subunit is essential for restriction, but not modification. The current model for restriction by both types of enzyme postulates that DNA cutting is stimulated when two enzyme complexes bound to neighbouring target sequences meet as the consequence of ATP-dependent DNA translocation. For type I enzymes, this model is supported by in vitro experiments, but the predicted co-operative interactions between targets have not been detected by assays that monitor restriction in vivo. The experiments reported here clearly establish the required synergistic effect but, in contrast to earlier experiments, they use Escherichia coli K-12 strains deficient in the restriction alleviation function associated with the Rac prophage. In bacteria with elevated levels of EcoKI the co-operative interactions are obscured, consistent with co-operation between free enzyme and that bound at target sites. We have made changes in three of the motifs characteristic of DEAD box proteins, including motif III, which in RecG is implicated in the migration of Holliday junctions. Conservative changes in each of the three motifs impair restriction.

The diversity of alleles at the hsd locus in natural populations of Escherichia coli

In enteric bacteria three discrete families of type I restriction and modification systems (IA, IB and ID) are encoded by alleles of the serB-linked hsd locus. Probes specific for each of the three families were used to monitor the distribution of related systems in 37 of the 72 wild-type Escherichia coli strains comprising the ECOR collection. All 25 members of group A in this collection were screened; 12 were probe-positive, nine have hsd genes in the IA family, two in the IB and one in the ID. Twelve strains, representing all groups other than A, were screened; five were probe-positive, one has hsd genes in the IA family, one in the IB and three in the ID. The type ID genes are the first representatives of this family in E. coli, the probe-negative strains could have alternative families of hsd genes. The type IA and IB systems added at least five new specificities to the five already identified in natural isolates of E. coli. The distribution of alleles is inconsistent with the dendrogram of the bacterial strains derived from other criteria. This discrepancy and the dissimilar coding sequences of allelic hsd genes both imply lateral transfer of hsd genes.

Genetic map of Salmonella typhimurium, edition VIII

We present edition VIII of the genetic map of Salmonella typhimurium LT2. We list a total of 1,159 genes, 1,080 of which have been located on the circular chromosome and 29 of which are on pSLT, the 90-kb plasmid usually found in LT2 lines. The remaining 50 genes are not yet mapped. The coordinate system used in this edition is neither minutes of transfer time in conjugation crosses nor units representing "phage lengths" of DNA of the transducing phage P22, as used in earlier editions, but centisomes and kilobases based on physical analysis of the lengths of DNA segments between genes. Some of these lengths have been determined by digestion of DNA by rare-cutting endonucleases and separation of fragments by pulsed-field gel electrophoresis. Other lengths have been determined by analysis of DNA sequences in GenBank. We have constructed StySeq1, which incorporates all Salmonella DNA sequence data known to us. StySeq1 comprises over 548 kb of nonredundant chromosomal genomic sequences, representing 11.4% of the chromosome, which is estimated to be just over 4,800 kb in length. Most of these sequences were assigned locations on the chromosome, in some cases by analogy with mapped Escherichia coli sequences.

In vivo restriction by LlaI is encoded by three genes, arranged in an operon with llaIM, on the conjugative Lactococcus plasmid pTR2030

The LlaI restriction and modification (R/M) system is encoded on pTR2030, a 46.2-kb conjugative plasmid from Lactococcus lactis. The llaI methylase gene, sequenced previously, encodes a functional type IIS methylase and is located approximately 5 kb upstream from the abiA gene, encoding abortive phage resistance. In this study, the sequence of the region between llaIM and abiA was determined and revealed four consecutive open reading frames (ORFs). Northern (RNA) analysis showed that the four ORFs were part of a 7-kb operon with llaIM and the downstream abiA gene on a separate transcriptional unit. The deduced protein sequence of ORF2 revealed a P-loop consensus motif for ATP/GTP-binding sites and a three-part consensus motif for GTP-binding proteins. Data bank searches with the deduced protein sequences for all four ORFs revealed no homology except for ORF2 with MerB, in three regions that coincided with the GTP-binding motifs in both proteins. To phenotypically analyze the llaI operon, a 9.0-kb fragment was cloned into a high-copy-number lactococcal shuttle vector, pTRKH2. The resulting construct, pTRK370, exhibited a significantly higher level of in vivo restriction and modification in L. lactis NCK203 than the low-copy-number parental plasmid, pTR2030. A combination of deletion constructions and frameshift mutations indicated that the first three ORFs were involved in LlaI restriction, and they were therefore designated llaI.1, llaI.2, and llaI.3. Mutating llaI.1 completely abolished restriction, while disrupting llaI.2 or llaI.3 allowed an inefficient restriction of phage DNA to occur, manifested primarily by a variable plaque phenotype. ORF4 had no discernible effect on in vivo restriction. A frameshift mutation in llaIM proved lethal to L. lactis NCK203, implying that the restriction component was active without the modification subunit. These results suggested that the LlaI R/M system is unlike any other R/M system studied to date and has diverged from the type IIS class of restriction enzymes by acquiring some characteristics reminiscent of type I enzymes.

Transduction, restriction and recombination patterns in Escherichia coli

Chromosomal DNA from several Escherichia coli reference (ECOR) strains was transduced by bacteriophage P1 into E. coli strain K12 W3110 trpA33. Recombination patterns of the transductants were determined by restriction fragment length polymorphism over a 40-kb region centering on a single marker (trpA+) in the tryptophan operon. These experiments demonstrate that transduction between different strains of E. coli can result in recombinational replacements that are small in comparison to the entrant molecule (replacements average 8-14 kb, whereas P1 packages approximately 100 kb) often in a series of discrete segments. The transduction patterns generated resemble the natural mosaic sequence patterns of the ECOR strains described in previous work. Extensive polymorphisms in the restriction-modification systems of the ECOR strains are a possible explanation for the sequence patterns in nature. To test this possibility two transductants were back-transduced into strain K12 W3110 trpA33. The resulting patterns were strikingly different from the original transductions. The size of the replacements was greater, and no multiple replacements were observed, suggesting a role for restriction-modification systems in the transduction patterns and perhaps for the mosaic sequence patterns in nature.

Intergeneric transfer and recombination of the 6-phosphogluconate dehydrogenase gene (gnd) in enteric bacteria

The gnd gene, encoding 6-phosphogluconate dehydrogenase (EC 1.1.1.44), was sequenced in 87 strains of 15 species assigned to five nominal genera of the Enterobacteriaceae, including 36 isolates of Salmonella enterica and 32 strains of Escherichia coli. In S. enterica, the effective (realized) rate of recombination of horizontally transferred gnd sequences is only moderately higher than the rates for other chromosomal housekeeping genes. In contrast, recombination at gnd has occurred with such high frequency in Escherichia coli that the indicated evolutionary relationships among strains are not congruent with those estimated by sequence analysis of other genes and by multilocus enzyme electrophoresis. E. coli and S. enterica apparently have not exchanged gnd sequences, but those of several strains of E. coli have been imported from species of Citrobacter and Klebsiella. The relatively frequent exchange of gnd within and among taxonomic groups of the Enterobacteriaceae, compared with other housekeeping genes, apparently results from its close linkage with genes that are subject to diversifying selection, including those of the rfb region determining the structure of the O antigen polysaccharide.

Bacterial gene transfer by natural genetic transformation in the environment

Natural genetic transformation is the active uptake of free DNA by bacterial cells and the heritable incorporation of its genetic information. Since the famous discovery of transformation in Streptococcus pneumoniae by Griffith in 1928 and the demonstration of DNA as the transforming principle by Avery and coworkers in 1944, cellular processes involved in transformation have been studied extensively by in vitro experimentation with a few transformable species. Only more recently has it been considered that transformation may be a powerful mechanism of horizontal gene transfer in natural bacterial populations. In this review the current understanding of the biology of transformation is summarized to provide the platform on which aspects of bacterial transformation in water, soil, and sediments and the habitat of pathogens are discussed. Direct and indirect evidence for gene transfer routes by transformation within species and between different species will be presented, along with data suggesting that plasmids as well as chromosomal DNA are subject to genetic exchange via transformation. Experiments exploring the prerequisites for transformation in the environment, including the production and persistence of free DNA and factors important for the uptake of DNA by cells, will be compiled, as well as possible natural barriers to transformation. The efficiency of gene transfer by transformation in bacterial habitats is possibly genetically adjusted to submaximal levels. The fact that natural transformation has been detected among bacteria from all trophic and taxonomic groups including archaebacteria suggests that transformability evolved early in phylogeny. Probable functions of DNA uptake other than gene acquisition will be discussed. The body of information presently available suggests that transformation has a great impact on bacterial population dynamics as well as on bacterial evolution and speciation.

Sequence divergence of the murB and rrfB genes from Escherichia coli and Salmonella typhimurium

The murB gene of Salmonella typhimurium was cloned and found to be 75% and 82% identical to the DNA and protein sequences, respectively, of the same gene in Escherichia coli. These identities are among the lowest recorded between the two bacteria. Nevertheless, wild-type S. typhimurium murB complemented the known temperature-sensitive E. coli mutant, and wild-type E. coli murB complemented three temperature-sensitive mutants of S. typhimurium. The 5S rRNA gene, rrfB, and the region between murB and rrfB were also cloned and sequenced. The rrfB gene of S. typhimurium differs from rrfB of E. coli in only 2 of 120 nt, but the region between murB and rrfB has diverged greatly and includes a sequence that closely resembles a repetitive extragenic palindrome of the type normally associated with E. coli. Previous comparisons of gene divergence have suggested that the chromosomal mutation rate is lower in the vicinity of the origin of replication. However, the S. typhimurium murB gene, located 6 map minutes from the origin of replication, is highly substituted at synonymous sites and the sequence between murB and rrfB is significantly modified as well. Thus, murB is an exception to the general observation that genes near the origin of replication show less divergence than do genes elsewhere in the bacterial chromosome.

Evolution of the recA gene and the molecular phylogeny of bacteria

The DNA sequences of the recA gene from 25 strains of bacteria are known. The evolution of these recA gene sequences, and of the derived RecA protein sequences, is examined, with special reference to the effect of variations in genomic G + C content. From the aligned RecA protein sequences, phylogenetic trees have been drawn using both distance matrix and maximum parsimony methods. There is a broad concordance between these trees and those derived from other data (largely 16S ribosomal RNA sequences). There is a fair degree of certainty in the relationships among the "Purple" or Proteobacteria, but the branching pattern between higher taxa within the eubacteria cannot be reliably resolved with these data.

Conservation of motifs within the unusually variable polypeptide sequences of type I restriction and modification enzymes

Type I restriction enzymes comprise three subunits encoded by genes designated hsdR, hsdM, and hsdS; S confers sequence specificity. Three families of enzymes are known and within families, but not between, hsdM and hsdR are conserved. Consequently, interfamily comparisons of M and R sequences focus on regions of putative functional significance, while both inter- and intrafamily comparisons address the origin, nature and role of diversity of type I restriction systems. We have determined the sequence of the hsdR gene for EcoA, thus making available sequences of all three hsd genes of one representative from each family. The predicted R polypeptide sequences share conserved regions with one superfamily of putative helicases, so-called 'DEAD box' proteins; these conserved sequences may be associated with the ATP-dependent translocation of DNA that precedes restriction. We also present hsdM and hsdR sequences for EcoE, a member of the same family as EcoA. The sequences of the M and R genes of EcoA and EcoE are at least as divergent as typical genes from Escherichia coli and Salmonella, perhaps as the result of selection favouring diversity of restriction specificities combined with lateral transfer among different species.