Microbiology
HOME HELP FEEDBACK SUBSCRIPTIONS ARCHIVE SEARCH TABLE OF CONTENTS
 QUICK SEARCH:   [advanced]


     


Microbiology 153 (2007), 1897-1906; DOI  10.1099/mic.0.2006/003020-0
This Article
Right arrow Abstract Freely available
Right arrow Full Text (PDF)
Right arrow Supplementary data
Right arrow Alert me when this article is cited
Right arrow Alert me if a correction is posted
Right arrow Citation Map
Services
Right arrow Email this article to a friend
Right arrow Similar articles in this journal
Right arrow Similar articles in PubMed
Right arrow Alert me to new issues of the journal
Right arrow Download to citation manager
Right arrow reprints & permissions
Citing Articles
Right arrow Citing Articles via HighWire
Right arrow Citing Articles via CrossRef
Right arrow Citing Articles via Google Scholar
Google Scholar
Right arrow Articles by Rio, R. V. M.
Right arrow Articles by Graf, J.
Right arrow Search for Related Content
PubMed
Right arrow PubMed Citation
Right arrow Articles by Rio, R. V. M.
Right arrow Articles by Graf, J.
Agricola
Right arrow Articles by Rio, R. V. M.
Right arrow Articles by Graf, J.
Microbiology 153 (2007), 1897-1906; DOI  10.1099/mic.0.2006/003020-0
© 2007 Society for General Microbiology

Characterization of a catalase gene from Aeromonas veronii, the digestive-tract symbiont of the medicinal leech

Rita V. M. Rio1, Matthias Anderegg2 and Joerg Graf1,2

1 Department of Molecular and Cell Biology, University of Connecticut, 91 N. Eagleville Road, Unit-3125, Storrs, CT 06269-3125, USA
2 Institute for Infectious Diseases, University of Berne, Friedbühlstr. 51, CH-3010 Berne, Switzerland

Correspondence
Joerg Graf
joerg.graf{at}uconn.edu


    ABSTRACT
 TOP
 ABSTRACT
 INTRODUCTION
 METHODS
 RESULTS
 DISCUSSION
 REFERENCES
 
The catalase gene katA of the medicinal leech symbiont Aeromonas veronii bv. sobria was cloned, sequenced, and functionally characterized. Southern hybridization, using an A. veronii katA-specific hybridization probe, suggested the presence of a single gene copy in many Aeromonas species. A. veronii katA consisted of 1446 nt encoding a protein with a high degree of similarity to the small-subunit group III bacterial catalases. A catalase-null mutant (JG186) was constructed through gene-replacement mutagenesis. In the parent strain (HM21R), catalase activity was only detected in extracts of cells grown to early exponential phase following H2O2 induction, in which the ability to induce activity was inversely related to optical density. In contrast, induced JG186 cells were very sensitive to oxidative stress, with survival being affected even at low H2O2 concentrations. In contrast to the findings of previous reports of other symbiotic systems, the catalase mutant was not defective in its ability to competitively colonize or persist within its host, in both co-inoculation and sole-colonization assays. This body of evidence suggests either that oxidative stress, in the form of H2O2 exposure, is not encountered by the microbial partner under the examined symbiotic conditions or that compensatory mechanisms exist. The data suggest that although many colonization factors reoccur, each symbiotic system has also evolved specific mechanisms that affect symbiont–host dynamics.


Abbreviations: CI, competitive index; ROS, reactive oxygen species

The GenBank/EMBL/DDBJ accession no. for the A. veronii katA sequence is EF028076.

A figure showing the results of transformation of E. coli ZK918 with purified plasmid DNA from the A. veronii genomic library, after plating on MacConkey agar, is available as supplementary data with the online version of this paper.


    INTRODUCTION
 TOP
 ABSTRACT
 INTRODUCTION
 METHODS
 RESULTS
 DISCUSSION
 REFERENCES
 
Physiological adaptations typically reflect the selection pressures encountered within ecological niches. For example, catalases are responsible for the decomposition of hydrogen peroxide to water and oxygen (2H2O2->2H2O+O2), and have presumably evolved to protect against cellular damage arising from exposure to H2O2. The two main factors that influence the expression of catalase gene(s) in bacteria are the exposure to sublethal H2O2 levels during exponential growth (Barnes et al., 1999Down) and entry into stationary phase (reviewed by Loewen, 1997Down). Reactive oxygen species (ROS) may also be powerful microbicidal weapons implemented by eukaryotic cells during infection. Highly diffusible oxidants, such as H2O2, enter membranes and can damage a variety of cell targets (Franzon et al., 1990Down; Hassett & Cohen, 1989Down), including DNA, RNA, proteins and lipids, leading to the destruction of microbes. In addition, ROS are generated continuously during aerobic growth by the electron transport chain, leading to the production of H2O2 (Cabiscol et al., 2000Down; Gonzalez-Flecha & Demple, 1995Down). The toxicity of H2O2 necessitates that aerobic microbes or those whose life cycle relies on host infection encode a functional catalase within their genomes.

Aeromonas veronii is a facultative anaerobe which has a propensity to colonize the digestive tracts of a variety of hosts, including humans, leeches and mosquitoes, with manifestations of infection ranging from pathogenesis to mutualism (Graf et al., 2006Down; Janda & Abbott, 1998Down). Recently culture-independent analysis of the digestive-tract microbiota of the medicinal leech (Hirudo verbana, one of three species within the Hirudo medicinalis complex; Hirudinea: Arhynchobdellida: Hirudinidae) (Apakupakul et al., 1999Down; Siddall & Burreson, 1998Down; Siddall et al., 2001Down; Trontelj & Utevsky, 2005Down) has revealed the presence of a currently uncultured Rikenella-like species (Worthen et al., 2006Down; Kikuchi & Graf, 2007Down). These two symbionts are the dominant members of the leech digestive-tract microbiota. Putative functions for the Aeromonas symbiont within the medicinal leech include: (i) aiding in the digestion of the blood meal; (ii) providing essential nutrients lacking in the exclusive blood diet, such as B-complex vitamins; (iii) providing ‘colonization resistance’, in which A. veronii prevents colonization by other potentially harmful micro-organisms, thus preventing the putrefaction of blood and permitting long-term storage; and (iv) priming the microenvironment for the obligate anaerobic Rikenella-like symbiont (Worthen et al., 2006Down; reviewed by Graf, 2002Down).

Analogous to vertebrate innate immune responses, the introduction of Gram-negative bacteria and their by-products into leech wounds activates an inflammatory response involving the infiltration of macrophage-like cells (de Eguileor et al., 1999Down, 2000aDown, bDown). During phagocytosis or following stimulation with a wide variety of agents, macrophages undergo respiratory bursts that are characterized by the production and release of ROS into the extracellular milieu (Forman & Torres, 2002Down; Park, 2003Down). ROS production is essential for the increased bactericidal capability of stimulated macrophages (Johnston & Kitagawa, 1985Down). Bacteria can protect themselves against host-produced ROS by upregulating the genes encoding protective enzymes, such as superoxide dismutase, peroxidase and catalase. While the importance of these protective enzymes is well established for pathogens (Franzon et al., 1990Down; Mandell, 1975Down; Zheng et al., 1992Down), they have also been shown to be crucial for the successful colonization of the light organ of the Hawaiian bobtail squid Euprymna scolopes by its extracellular, mutualistic symbiont Vibrio fischeri (Visick & Ruby, 1998Down), and for the regulation of the infection of root nodules by the nitrogen-fixing Sinorhizobium meliloti (Santos et al., 2001Down).

We report the cloning and functional characterization of the A. veronii gene encoding the antioxidant enzyme catalase katA. We further examine the presence of katA homologues in other Aeromonas species. Implementing a targeted mutant analysis approach, the importance of katA for A. veronii survival following exposure to damaging oxidative-stress conditions, and for proliferation and persistence within the medicinal leech, was evaluated.


    METHODS
 TOP
 ABSTRACT
 INTRODUCTION
 METHODS
 RESULTS
 DISCUSSION
 REFERENCES
 
Bacterial strains, plasmids and growth conditions.
The sources and properties of the bacterial strains and plasmids used in this study are listed in Table 1Down. The bacteria were grown at 200 r.p.m. in Luria–Bertani (LB) medium at 37 and 30 °C for Escherichia coli and Aeromonas species, respectively (Sambrook & Russell, 2001Down). The growth rates were determined in LB containing rifampicin (Rf) at 30 °C.


View this table:
[in this window]
[in a new window]

 
Table 1. Strains and plasmids used in this study

 
The A. veronii strain HM21R is a spontaneous Rfr mutant derived from HM21, an isolate from the leech digestive tract (Graf, 1999Down). HM21RS was derived from HM21R by plating 108 cells from an overnight culture on LB [Rf, streptomycin (Sm)] plates. Where appropriate, antibiotics were added at the following concentrations: ampicillin (Ap), 100 µg ml–1; chloramphenicol (Cm), 1 µg ml–1 for A. veronii and 10 µg ml–1 for E. coli; kanamycin (Km), 100 µg ml–1; Sm, 100 µg ml–1; and Rf, 100 µg ml–1 for Aeromonas selection and 10 µg ml–1 for maintenance.

Animals.
The medicinal leeches used in this study were obtained from LeechesUSA and Zaug GmbH (Biebertal, Germany). The animals were starved for at least 3 months prior to delivery and maintained without feeding in leech tanks at 25±1 °C (Graf, 1999Down).

Isolation of genomic DNA.
Genomic DNA was isolated using a modification of the cetyltrimethylammonium bromide (CTAB) method (Nelson & Selander, 1994Down) by preheating the samples to 65 °C for 5 min prior to lysis with SDS. For large-scale DNA isolation, the volume was increased 30-fold, samples were preheated to 65 °C for 10 min, and the precipitated DNA was isolated using a glass Pasteur pipette.

Construction of A. veronii genomic library.
Genomic DNA was partially digested with Sau3A and separated by electrophoresis in low-melting-point agarose. DNA fragments (~6.5–7.6 kb) were purified by phenol extraction and cloned into the dephosphorylated BamHI site of pBSIIKS+. The ligations were transformed into E. coli XL-1 Blue MRF' super competent cells (Stratagene) through heat shock at 42 °C.

Complementation of the E. coli rpoS mutant.
Purified plasmid DNA from the A. veronii genomic library was transformed into calcium-competent E. coli rpoS mutant ZK918 (Sambrook & Russell, 2001Down). Transformants were plated on MacConkey agar (Ap, Km) and screened for red colonies. Strains exhibiting this phenotype were picked and patched onto MacConkey agar plates, where the colour phenotype was reassessed. The colonies were additionally examined for the production of catalase by observing the release of O2 bubbles following the addition of 3 % H2O2 (Supplementary Fig. S1).

Characterization of the complementing plasmids.
The complementing plasmids were isolated using Plasmid Mini or Plasmid Midi kits (Qiagen) and sequenced using a combination of subcloning and primer walking (Central DNA Sequencing Facility, Department of Clinical Research, University of Berne, and the Biotechnology Center, University of Connecticut). The DNA sequences were aligned using Vector NTI (Invitrogen) and compared to the databases using the BLAST 2.2.13 algorithm (Altschul et al., 1997Down).

Southern analysis.
DNA (Aeromonas hydrophila ATCC 14715, Aeromonas salmonicida CDC 0434-84, Aeromonas caviae ATCC 15468, Aeromonas media CDC 0862-83, Aeromonas eucrenophila ATCC 23309, Aeromonas sobria CIP 7433, A. veronii bv. sobria CDC 0437-84, Aeromonas veronii bv. sobria HM21, A. veronii bv. veronii ATCC 35624, Aeromonas schubertii ATCC 43700 and Aeromonas allosaccharophila LMG 140549) was digested with PstI for 6 h and separated by agarose gel electrophoresis. The loading of an equivalent amount of DNA in the various lanes was confirmed by ethidium bromide staining. The DNA was transferred onto a nylon membrane through overnight capillary transfer and fixed to the membrane with 0.4 M NaOH (Sambrook & Russell, 2001Down).

A 743 bp hybridization probe was PCR-amplified from HM21 DNA using PkatF, 5'-TCG ACA ACA ACA ACA GCC TCA C-3', and PkatR1, 5'-CAC CTG CAC ATA GAC TTT CCA GC-3', as described below, and subsequently labelled using the Amersham ECL Direct System (Amersham Life Science). The membrane was prehybridized at 42 °C in the Amersham ECL direct hybridization buffer. After 1 h, the probe was added and allowed to hybridize overnight. The membrane was washed twice for 20 min at 42 °C in the primary wash buffer and twice for 5 min in the secondary wash buffer. The membrane was incubated with detection reagents 1 and 2 for 1 min and exposed to X-ray film.

Construction of the catalase mutants.
For the construction of the catalase mutants, two internal fragments of katA were amplified by PCR using PkatF and either PkatR1 (amplifying a 743 bp product and resulting in {Delta}katAS) or PkatR2, 5'-AGG AGA AGA GTC GCC CTT G-3' (amplifying a 966 bp product and resulting in {Delta}katAL). Single amplicons were obtained for each of the two primer sets. The PCR products were blunt-ended using T4 DNA polymerase and ligated into SmaI-digested pBSII, yielding pJG54 (containing the shorter PCR product, {Delta}katAS) and pJG55 (containing the longer PCR product, {Delta}katAL). The cloned fragments were introduced into the {pi}-dependent R6K origin of replication suicide vector pKAS32, using XbaI–EcoRI ends, resulting in pJG57 ({Delta}katAS) and pJG58 ({Delta}katAL). The suicide plasmids were transformed into S17-1 {lambda}pir and introduced into HM21R by conjugation with shaking (Simon et al., 1983Down). Donor (5x107 c.f.u. between 0.4 and 0.8 OD600) and recipient (2x108 c.f.u. between 0.25 and 0.4 OD600) cells were harvested and spotted on an LB agar plate. Following an overnight incubation at 30 °C, transconjugants were selected on LB agar (Rf, Km).

Zymography.
Whole-cell lysates were prepared from the cells as described in Methods, Sensitivity to H2O2, below. The cells were lysed by sonication (three bursts, 15 s) and centrifuged at 4 °C for 30 min (12 000 g). Supernatants were loaded onto a 10 %, w/v, acrylamide gel and separated at 150 V for 45 min. The gel was washed extensively in double-distilled (dd) H2O and for 10 min in 0.0155 M H2O2. Subsequently, the gel was stained in a solution containing 1 %, w/v, ferric chloride and 1 %, w/v, potassium ferric cyanide for the visualization of catalase activity (Barnes et al., 1999Down).

Analysis of catalase activity.
Paired 5 ml aliquots of cultures were removed at various optical densities and either induced with 0.05 mM H2O2 or not induced. The aliquots were incubated for an additional 30 min before being placed on ice. A crude cell extract was prepared by sonication. Specific catalase activity was determined by measuring the removal of H2O2 from the crude cell extract for 60 s and dividing by the amount of protein, as determined by a Lowry protein assay (Lowry et al., 1951Down; Beers & Sizer, 1952Down; Visick & Ruby, 1998Down). Briefly, the cell pellet was resuspended in phosphate buffer, pH 7.0, containing 5 mM EDTA, 10 %, v/v, glycerol and 25 µM PMSF. The change in A240 was monitored four consecutive times (15 s per period) following the addition of 0.5 volumes of 59 mM H2O2 to 1 volume of crude extract. The mean rates were converted to U mg–1 min–1.

Reverse transcription of katA.
Total RNA was isolated from 2.5x108 A. veronii (JG186 and HM21R) cells grown to stationary or early mid-exponential phase, employing the Qiagen RNeasy protocol for total RNA isolation from bacteria. An aliquot of the mid-exponential-phase cells was exposed to a sublethal dose of H2O2 (0.05 mM) for 30 min at 30 °C. The optional on-column RNase-free DNase I (Qiagen) was used to remove contaminating DNA. After RNA isolation, traces of contaminating DNA were further eliminated with an RNase-free DNase I treatment. Random hexamer primers, iScript reverse transcriptase (Bio-Rad) and 0.5 µg total RNA were utilized for first-strand cDNA synthesis. Subsequent PCR was performed with katAF', 5'-GAC AAG ACC CTG CAC AGC-3', and katAR', 5'-CGC TCA TTG GCG TCG TTG-3'. The absence of DNA contamination was verified by PCR using RNA template lacking a reverse-transcription step.

Sensitivity to H2O2.
The sensitivity of the katA mutant JG186 and its parent strain HM21R to H2O2 was determined for cells prepared in three different ways. The cells were grown either to stationary or to early exponential phase. A portion of the early exponential-phase cells was exposed to a sublethal dose of H2O2, as described above. Cells from both strains and at each growth condition were diluted to the same density (5x106 c.f.u. ml–1), challenged with various H2O2 concentrations for 30 min at 30 °C, serially diluted and plated. The percentage survival was determined by dividing the count obtained following the exposure to H2O2 of each strain by the count obtained prior to challenge.

Complementation of JG186.
A 1.8 kb amplicon that contained katA and its promoter and termination regions was PCR-amplified using katAFullF', 5'-GAC CAC ATC ACC GTT CTC CA-3', and katAFullR', 5'-TGC TCG CAA TAG AAA ACG GG-3'. The amplicon was cloned into pCR 2.1 (Invitrogen), resulting in pRRJG5. The cloned 1.8 kb fragment was moved into the broad-host-range vector pMMB207 (Morales et al., 1991Down) from pRRJG5 by utilizing the EcoRI sites. The resulting plasmid, pRRSR1, was electroporated into E. coli DH5{alpha} {lambda} pir cells. The presence of the insert was verified through DNA sequencing and restriction digests. pRRSR1 and pMMB207 were each introduced into HM21R and JG186 through conjugation. The resulting strains were induced and challenged with 0.8 and 1.0 mM H2O2, respectively, as described above.

Colonization assay.
The competition assay used in this study, which compared the ability of JG186 and a competitor strain to colonize the medicinal leech, was similar to the assay we described previously (Graf, 1999Down), except that the two strains were simultaneously fed to the leech in a single blood meal. Briefly, leeches were fed heparinized sheep blood (25 U ml–1, Sigma) that had been stored for 18–24 h at 23 °C. Cells (500 c.f.u. ml –1 of JG186 and HM21RS) were added to blood that had been preheated to 56 °C, which inactivates the antimicrobial properties of the complement system, which can kill sensitive bacteria (Indergand & Graf, 2000Down). The animals were fed 5 ml blood and individuals were subsequently incubated at 23 °C. At 6 and 18 h post-blood meal, animals were sacrificed and the intraluminal fluid was collected from the crop of the digestive tract. Serially diluted samples were plated on LB plates containing the appropriate antibiotics and incubated overnight at 30 °C, allowing for comparative counts of the introduced microbiota. The limit of detection was 10 c.f.u. ml–1. Similarly, sole-colonization assays were performed to rule out the possibility that the katA mutation was circumvented by the endogenous catalase of HM21RS. For the sole-colonization assay, the inoculum was between 200 and 500 c.f.u. ml–1 of either JG186 or HM21R, and animals were sacrificed 18 h post-blood meal to determine colonization levels. For each time point, at least three animals were used.

Statistical analysis.
The data were analysed using Graph Pad Prism 2.01. A competitive index (CI) was calculated in order to compare test strain concentrations in the intraluminal fluid following their introduction. The absolute and relative levels of each strain were monitored over time by sacrificing animals at predetermined time points and plating serial dilutions on antibiotic-containing plates. CI was calculated as follows:

Formula
A CI of 1 indicates that the mutant colonizes to the same level as the competitor strain, while CI <1 indicates that the mutant has a colonization defect. A two-tailed, one-sided t test was used to test whether the CI differed from 1. Student’s t test was performed to determine whether colonization levels significantly differed for mutants introduced in sole-colonization assays versus competition assays. Significant differences (P≤0.05) are reported.


    RESULTS
 TOP
 ABSTRACT
 INTRODUCTION
 METHODS
 RESULTS
 DISCUSSION
 REFERENCES
 
Isolation of plasmids that complement the rpoS mutant ZK918
The initial goal of this study was to isolate a functional rpoS homologue from a genomic library of A. veronii bv. sobria by complementing the E. coli rpoS mutant ZK918 (Bohannon et al., 1991Down). This mutant has two phenotypes that can be easily screened for and that have been complemented in other studies by rpoS homologues from Vibrio cholerae and Pseudomonas putida (Ramos-Gonzalez & Molin, 1998Down; Yildiz & Schoolnik, 1998Down). One of these phenotypes is the red colouration of colonies on MacConkey agar due to lactose fermentation that is dependent on the transcription of the {sigma}S-dependent promoter fusion bolAp1lacZ. The second phenotype is the release of O2 bubbles by colonies after the addition of H2O2 due to the expression of the {sigma}S-dependent HPII catalase katE.

Competent ZK918 were transformed with a genomic library from A. veronii. Of the 9000 transformants screened, three exhibited the red colouration and release of O2 after the addition of H2O2. The complementing plasmids pMA2, pMA5 and pEP2 were isolated and reintroduced into the rpoS mutant ZK918. All three plasmids again complemented ZK918, indicating that complementation was due to the plasmids and not to a secondary chromosomal mutation in the transformed strain. Analysis of these plasmids by restriction endonuclease digestion was consistent with all three plasmids carrying similar chromosomal regions and was confirmed by sequencing the ends of the inserted genomic DNA. Because all three plasmids complemented ZK918 and contained essentially the same chromosomal region, we focused on pMA5 in our subsequent analysis.

The region of pMA5 that was responsible for complementing the E. coli rpoS mutant was identified by deleting regions of the inserted DNA, transforming the subcloned plasmids into ZK918, and testing for red colony colour as well as the release of O2 bubbles following H2O2 exposure (Fig. 1aDown). Interestingly, transformation with pMA10 resulted in a red colour but failed to restore the ability of ZK918 to produce O2 bubbles (Supplementary Fig. S1), suggesting that two different loci may actually be responsible for the two phenotypes.


Figure 1
View larger version (38K):
[in this window]
[in a new window]

 
Fig. 1. (a) The region of pMA5 responsible for complementing the E. coli rpoS mutant was identified by deleting regions of the inserted DNA, transforming the subcloned plasmids into ZK918 and testing for red colony colour as well as the release of O2 bubbles following exposure to H2O2. The smallest plasmid that restored both phenotypes of ZK918 was pMA7. (b) Genetic map of the katA region of pMA7. Sequence analysis revealed one large ORF (ORF1) encoding katA, and two smaller ORFs (ORF2 and ORF3). Large arrows indicate the locations and directions of transcription of the identified genes. (c) Deduced amino acid sequence of A. veronii katA. The boxed region denotes catalase region 1, and the conserved tyrosine residue is in white on black type. The underlined region denotes catalase region 2, part of the catalytic site that contains a conserved histidine residue (asterisked).

 
DNA sequence analysis
The complementing region was sequenced in both directions using a combination of primer walking and subcloning. Database searches with the sequence of the cloned Aeromonas DNA using BLASTN and BLASTX showed no similarity to any sigma factors, suggesting that the inserted DNA does not encode a sigma factor. Further sequence analysis revealed one large ORF (ORF1) and two smaller ORFs (ORF2 and ORF3) that are transcribed in the opposite direction (Fig. 1bUp).

ORF1 consisted of 1446 nt encoding 482 amino acids with a predicted molecular mass of 54 kDa and isolectric point (pI) of 5.98. Putative –10 and –35 regions for the E. coli housekeeping sigma factor were found 27 and 52 bp, respectively, upstream of the presumptive start codon. A potential Shine–Dalgarno sequence (AGGAGA) was also detected 7 bp upstream of the presumptive start codon. Downstream of the stop codon a GC-rich interrupted dyad symmetry followed by a run of Ts was found, consistent with the stem–loop structure of a rho-independent terminator.

Comparison of the nucleotide sequence from ORF1 with the public databases using BLASTX revealed sequence similarity to the group III small-subunit bacterial catalases (Klotz et al., 1997Down), including HktE from Pasteurella multocida (70 % identity and 81 % similarity), a putative catalase from Photobacterium profundum (71 % identity and 79 % similarity) and the catalases of Bordetella pertussis, Bordetella parapertussis and Bordetella bronchiseptica (69 % identity and 83 % similarity). Furthermore, the deduced amino acid sequence was conserved along the entire protein length with a slightly lower degree of sequence identity at the carboxyl terminus. In addition, a PROSITE search revealed that the two regions conserved in functional catalases are also present in A. veronii protein. Catalase region 1 is responsible for haem binding and contains a conserved tyrosine residue (RLFSY.DTQ), while catalase region 2, a section of the catalytic site, contains a conserved histidine residue (F.R....ER..H..GSG) (Fig. 1cUp). In accordance with the high sequence identity, the retention of conserved residues and the release of bubbles following the addition of H2O2 in the complemented rpoS mutant ZK918, we designated the gene as katA (accession no. EF028076).

ORF2 had sequence similarity to a putative acetyltransferase of Bacillus subtilis, whereas ORF3 was homologous to a DNA gyrase inhibitor. A comparison of the A. veronii katA locus with the A. hydrophila genome revealed conservation in genetic organization (Seshadri et al., 2006Down).

Presence of catalase genes in Aeromonas species
We were interested in determining whether other Aeromonas species possess a similar catalase gene. An internal fragment of katA was amplified and used as a probe for Southern analysis of nine Aeromonas species (Fig. 2Down). With the exception of A. schubertii, all of the Aeromonas species examined revealed one band that hybridized to the katA probe, suggesting that these eight species possess a similar catalase gene. Interestingly, the Aeromonas leech isolate HM21 katA band co-migrated with that of A. veronii bv. veronii rather than that of A. veronii bv. sobria.


Figure 2
View larger version (83K):
[in this window]
[in a new window]

 
Fig. 2. Southern blot analysis of Aeromonas species DNA hybridized with a 743 bp A. veronii katA fragment. Ah, A. hydrophila ATCC 14715; As, A. salmonicida CDC 0434-84; Ac, A. caviae ATCC 15468; Am, A. media CDC 08262-83; Ae, A. eucrenophila ATCC 23309; Aso, A. sobria CIP 7433; Avs, A. veronii bv. sobria CDC 04374-84; Avv, A. veronii bv. veronii ATCC 35624; Asc; A. schubertii ATCC 43700; Aa, A. allosaccharophila LMG 140549; Hm, A. veronii bv. sobria HM21.

 
Construction of katA mutants
Deletion mutants were constructed by introducing internal fragments of katA, located on a suicide vector, through homologous recombination into the A. veronii chromosome. This results in two copies of katA with each having a deletion in either the 5' or the 3' region. The presumptive mutants were screened for inability to release air bubbles after the addition of H2O2, and the disruption of the katA gene in four of the Kmr and catalase-minus strains, JG185 and JG186 ({Delta}katAL derived) and JG183 and JG184 ({Delta}katAS derived), was confirmed by Southern blotting using a katA-specific probe. A shift in the molecular mass of the bands that hybridized with the katA probe was observed, indicating that the plasmid had incorporated into the katA locus (data not shown). We chose JG186 for further functional characterization. The growth rate of JG186 in LB did not differ from that of HM21R and HM21RS under similar conditions.

Detection of catalase activity
Facultative or obligate aerobic bacteria can possess multiple catalases. Our Southern analysis suggested the presence of only one catalase, but we wanted to verify catalase activity independently of sequence similarity through zymography. Whole-cell lysates of the parent strain and JG186 were obtained from cells grown to early exponential phase and stationary phase, and from early exponential-phase cells that were exposed to a sublethal concentration of H2O2 for 30 min (induced). Catalase activity was detected only from lysates of induced early exponential-phase HM21R cells (Fig. 3Down). No other samples exhibited catalase activity, further supporting the presence of only one A. veronii catalase that is inducible by oxidative stress. Furthermore, these results also demonstrated that JG186, the katA mutant, does not produce detectable levels of catalase.


Figure 3
View larger version (44K):
[in this window]
[in a new window]

 
Fig. 3. Zymogram analysis of catalase activity of HM21R and JG186 through cell growth. Catalase activity was specifically detected from HM21R early exponential-phase lysates that had been induced and subsequently challenged with H2O2. No other samples demonstrated catalase activity. Stat., stationary; EL, early exponential; Ind. EL, induced early exponential.

 
The specific catalase activity (U mg–1 min–1) of the parent strain and JG186, with and without exposure to sublethal levels of H2O2, was monitored over a range of cell densities (Fig. 4Down and data not shown). No catalase activity above background was ever detected for JG186. Interestingly, the ability of the wild-type strain to induce the catalase activity was inversely related to the optical density of the culture. It is possible that the specific catalase activity from cells at low optical density can be elevated by an underestimation of total protein level. Reverse-transcriptional analysis further confirmed katA expression by H2O2-induced HM21R (data not shown).


Figure 4
View larger version (14K):
[in this window]
[in a new window]

 
Fig. 4. Catalase activity (U mg–1 min–1) of HM21R and JG186 over a range of cell densities (OD600 0.1–2.7). Catalase activity was determined following the addition of H2O2, with and without prior exposure to sublethal levels. Error bars show SD, unless smaller than the symbols. Ind., induced; Non-ind, non-induced.

 
Sensitivity to H2O2
The sensitivity and protective response of the wild-type strain (HM21R) to oxidative stress was determined by exposing the cells to increasing concentrations of H2O2 (0.01–10 mM) and monitoring their survival. The cells were grown to early exponential phase or stationary phase, and subsequently a portion of the early exponential-phase cells was exposed to a sublethal concentration of H2O2 (induced cells). Non-induced HM21R early exponential-phase cells (Fig. 5aDown) were less resistant to increasing H2O2 concentrations than induced (Fig. 5bDown) and stationary-phase (Fig. 5cDown) cells. Both induced and stationary-phase HM21R were adversely affected commencing at 1 mM H2O2, with survival rates of 88 and 70 %, respectively (Fig. 5b, cDown).


Figure 5
View larger version (12K):
[in this window]
[in a new window]

 
Fig. 5. Representative percentage survival of (a) early exponential-phase HM21R and JG186, (b) induced early exponential-phase HM21R and JG186, and (c) stationary-phase HM21R and JG186, following exposure to various H2O2 concentrations (0.01–10 mM), based on three independent trials. EL, early exponential; Ind. EL, induced early exponential; Stat., stationary.

 
We then determined the importance of katA in providing oxidative stress protection by exposing JG186 to similar concentrations of H2O2. Interestingly, no survival differences between JG186 and HM21R were detected in stationary-phase cells (Fig. 5cUp). In contrast to induced HM21R, induced JG186 cells were much more sensitive to H2O2, with survival negatively affected at concentrations as low as 0.1 mM. These results demonstrated that katA is critical in providing an inducible protection against exogenous H2O2 in vitro. However, in a similar manner to HM21R, non-induced early exponential-phase JG186 cells were less resistant than their induced counterparts, with survival dropping to 2 % (Fig. 5aUp) in comparison to 38 % (Fig. 5bUp) at 0.5 mM H2O2. These results suggest that A. veronii bv. sobria retains another inducible mechanism in addition to katA to safeguard against oxidative stress, and that during stationary phase, protection is provided through catalase-independent means.

Colonization of the medicinal leech
In symbiotic relationships, mechanisms that control the spatial (anatomical) localization and the proliferation of the symbiotic flora are critical for the maintenance of homeostasis (Rio et al., 2006Down). One common mechanism employed to control infections of pathogenic bacteria is to create an environment of oxidative stress that kills sensitive bacteria. Catalase has been shown to be required for the normal symbiotic competence of V. fischeri and Eu. scolopes, being induced by both oxidative stress and the approach to stationary phase (Visick & Ruby, 1998Down). Accordingly, we were interested to determine whether the loss of katA would lead to a reduced ability of A. veronii to colonize H. medicinalis. Using the competition assay, no significant differences in the c.f.u. ml–1 of JG186 and HM21RS at 6 h (a two-tailed, one-sided t test was used to test whether the CI differed from 1, with P=0.1827) following introduction were detected, suggesting that at least under these conditions, the loss of katA was not essential for symbiotic establishment (data not shown). Furthermore, no significant differences in the c.f.u. ml–1 of JG186 and HM21RS at 18 h following feeding (a two-tailed, one-sided t test was used to test whether the CI differed from 1, with P=0.3762) were detected, indicating that symbiotic persistence remains unaltered following the loss of catalase activity (data not shown). Additionally, a sole-colonization assay indicated no difference in JG186 colonization levels in comparison to HM21R at 18 h, suggesting that the mutant did not rely on the catalase activity of the competitor strain for colonization (P=0.27) (data not shown).


    DISCUSSION
 TOP
 ABSTRACT
 INTRODUCTION
 METHODS
 RESULTS
 DISCUSSION
 REFERENCES
 
The A. veronii katA belongs to the group III catalases, which are predominantly found in proteobacteria, many of which lead a parasitic lifestyle (Klotz et al., 1997Down). It is interesting to note that bacteria which harbour a group III catalase and inhabit a restricted environment typically contain only a single catalase isozyme (Klotz et al., 1997Down), supporting our finding of only one catalase gene, katA, within the A. veronii genome. Similarly, the A. hydrophila genome contains only one katA gene (Seshadri et al., 2006Down), and zymographic analysis of A. salmonicida also reveals only one catalase (Barnes et al., 1999Down). Although this study found no definitive role for A. veronii katA in the medicinal leech symbiosis, it is likely that katA has functional importance to counter ROS produced as by-products during aerobic growth (Loewen, 1997Down). The lack of a hybridization signal from A. schubertii strain ATCC 43700 suggests either that the katA homologue is undetectable with the probe utilized due to the evolutionary distance from A. veronii or possibly that a homologue is lacking (Yanez et al., 2003Down; Soler et al., 2004Down). The differences observed in katA hybridization patterns provide a promising tool for distinguishing the two A. veronii bvs (sobria versus veronii) that could previously only be differentiated using an array of biochemical assays (Graf, 1999Down) or by DNA sequencing of gyrB and rpoD (Yanez et al., 2003Down; Soler et al., 2004Down).

Significant protection against oxidant killing was provided by katA during H2O2 challenges, particularly during exponential growth. Correspondingly, catalase activity was detected only from early exponential-phase HM21R cell lysates that had been induced with a sublethal concentration of H2O2. These induced cells were more protected against oxidative stress when subsequently challenged with even greater H2O2 concentrations. Similar challenges proved detrimental to non-induced HM21R and the katA mutant JG186, supporting the importance of both catalase activity and an adaptive antioxidant response that facilitates protection against subsequent lethal doses. Interestingly, the induction of katA was detected only in cells during early exponential-phase growth, differing from that reported for A. salmonicida, in which induction during mid-exponential and early stationary phase provides the greatest protection against H2O2 challenge (Barnes et al., 1999Down). Also, A. salmonicida is protected against higher levels of H2O2 (Barnes et al., 1999Down), although this could be due to differences in experimental design. In addition to the protection provided by KatA, a catalase-independent defensive mechanism against exogenous H2O2 was observed in stationary-phase A. veronii cells. This additional protection may be due to the upregulation of other antioxidant enzymes, such as the glutathione and thiol peroxidases that are present in the A. hydrophila genome (Seshadri et al., 2006Down). If such peroxidases are also present in the A. veronii genome, they could contribute to the protection detected during the stationary phase. Additionally, the low metabolic activity during stationary phase reduces the generation of ROS during respiration, and thus even basal levels of antioxidant enzymes may provide sufficient protection against oxidative stress (Barnes et al., 1999Down). The expression of A. veronii katA from pRRSR1 was unable to restore the wild-type phenotype. However, a polar effect of the katA mutation is unlikely, given that following the rho-independent terminator, according to the A. hydrophila genome, the downstream gene, a conserved hypothetical protein, is transcribed in the opposite direction. Our A. veronii draft genome sequence supports conservation in the genetic organization of the katA locus.

Oxygen and its reactive derivatives have been demonstrated to play significant roles in promoting the efficiency and specificity of non-digestive-tract symbiotic systems, such as those between S. meliloti and leguminous plants, photosynthetic algae and marine invertebrates, and V. fischeri and Eu. scolopes (reviewed by Ruby & McFall-Ngai, 1999Down). Unlike catalase mutants in these systems, the A. veronii catalase mutant was not defective in its ability to colonize or persist within its host digestive tract, at least under the conditions tested. This suggests either that oxidative stress, in the form of H2O2 exposure, is not encountered by the microbial partner during the colonization of the leech digestive tract or that other compensatory mechanisms exist, such as the previously mentioned peroxidases. Examining whether catalase is crucial for A. veronii infection in other associations, such as the cooperative zebrafish digestive tract, and whether virulence is attenuated in the mouse model may provide further insights. Similarly, B. pertussis does not require catalase for persistence within human polymorphonuclear leukocytes (PMNs), although phagocytosis is associated with a respiratory burst that involves the generation of O2 and H2O2 (DeShazer et al., 1994Down). It is possible that O2 introduced with the fresh blood meal is removed rapidly in the leech digestive tract by the aerobic metabolism of Aeromonas, generating an environment that does not permit the generation of ROS and is suitable for the anaerobic Rikenella-like symbiont that likely does not possess antioxidant enzymes. Synergistic interactions between the two microbes resulting in enhanced mixed-species, polysaccharide-embedded microcolony formation have been suggested to occur within the leech digestive tract (Kikuchi & Graf, 2007). The apparent lack of oxidative stress suggests that other mechanisms must be responsible for ensuring the specificity, establishment and maintenance of this unusually simple digestive-tract community.


    ACKNOWLEDGEMENTS
 
We thank Eva Pitek and Cuong Phang for help with initial experiments and Adam Silver and Rolf Troller for excellent technical assistance. This research was supported by the National Science Foundation Career award MCB 0448052 to J. G.

Edited by: J. G. Shaw


    REFERENCES
 TOP
 ABSTRACT
 INTRODUCTION
 METHODS
 RESULTS
 DISCUSSION
 REFERENCES
 
Altschul, S. F., Madden, T. L., Schäffer, A. A., Zhang, J., Zhang, Z., Miller, W. & Lipman, D. J. (1997). Gapped BLAST and PSI-BLAST: a new generation of protein database search programs. Nucleic Acids Res 25, 3389–3402.[Abstract/Free Full Text]

Apakupakul, K., Siddall, M. E. & Burreson, E. M. (1999). Higher level relationships of leeches (Annelida: Clitellata: Euhirudinea) based on morphology and gene sequences. Mol Phylogenet Evol 12, 350–359.[CrossRef][Medline]

Barnes, A. C., Bowden, T. J., Horne, M. T. & Ellis, A. E. (1999). Peroxide-inducible catalase in Aeromonas salmonicida subsp. salmonicida protects against exogenous hydrogen peroxide and killing by activated rainbow trout, Oncorhynchus mykiss L., macrophages. Microb Pathog 26, 149–158.[CrossRef][Medline]

Beers, R. F. & Sizer, I. W. (1952). A spectrophotometric method for measuring the breakdown of hydrogen peroxide by catalase. J Biol Chem 195, 133–140.[Free Full Text]

Bohannon, D. E., Connell, N., Keener, J., Tormo, A., Espinosa-Urgel, M., Zambrano, M. M. & Kolter, R. (1991). Stationary-phase-inducible ‘gearbox’ promoters: differential effects of katF mutations and role of sigma 70. J Bacteriol 173, 4482–4492.[Abstract/Free Full Text]

Cabiscol, E., Tamarit, J. & Ros, J. (2000). Oxidative stress in bacteria and protein damage by reactive oxygen species. Int Microbiol 3, 3–8.[Medline]

de Eguileor, M., Tettamanti, G., Grimaldi, A., Boselli, A., Scari, G., Valvassori, R., Cooper, E. L. & Lanzavecchia, G. (1999). Histopathological changes after induced injury in leeches. J Invertebr Pathol 74, 14–28.[CrossRef][Medline]

de Eguileor, M., Grimaldi, A., Tettamanti, G., Valvassori, R., Cooper, E. L. & Lanzavecchia, G. (2000a). Lipopolysaccharide-dependent induction of leech leukocytes that cross-react with vertebrate cellular differentiation markers. Tissue Cell 32, 437–445.[CrossRef][Medline]

de Eguileor, M., Grimaldi, A., Tettamanti, G., Valvassori, R., Cooper, E. L. & Lanzavecchia, G. (2000b). Different types of response to foreign antigens by leech leukocytes. Tissue Cell 32, 40–48.[CrossRef][Medline]

DeShazer, D., Wood, G. E. & Friedman, R. L. (1994). Molecular characterization of catalase from Bordetella pertussis: identification of the katA promoter in an upstream insertion sequence. Mol Microbiol 14, 123–130.[CrossRef][Medline]

Forman, H. J. & Torres, M. (2002). Reactive oxygen species and cell signaling: respiratory burst in macrophage signaling. Am J Respir Crit Care Med 166, S4–S8.[Abstract/Free Full Text]

Franzon, V. L., Arondel, J. & Sansonetti, P. J. (1990). Contribution of superoxide dismutase and catalase activities to Shigella flexneri pathogenesis. Infect Immun 58, 529–535.[Abstract/Free Full Text]

Gonzalez-Flecha, B. & Demple, B. (1995). Metabolic sources of hydrogen peroxide in aerobically growing Escherichia coli. J Biol Chem 270, 13681–13687.[Abstract/Free Full Text]

Graf, J. (1999). Symbiosis of Aeromonas veronii biovar sobria and Hirudo medicinalis, the medicinal leech: a novel model for digestive tract associations. Infect Immun 67, 1–7.[Abstract/Free Full Text]

Graf, J. (2002). The effects of symbionts on the physiology of Hirudo medicinalis, the medicinal leech. Invertebr Reprod Dev 41, 269–275.

Graf, J., Kikuchi, Y. & Rio, R. V. (2006). Leeches and their microbiota: naturally simple symbiosis models. Trends Microbiol 14, 365–371.[CrossRef][Medline]

Hassett, D. J. & Cohen, M. S. (1989). Bacterial adaptation to oxidative stress: implications for pathogenesis and interaction with phagocytic cells. FASEB J 3, 2574–2582.[Abstract]

Indergand, S. & Graf, J. (2000). Ingested blood contributes to the specificity of the symbiosis of Aeromonas veronii biovar sobria and Hirudo medicinalis, the medicinal leech. Appl Environ Microbiol 66, 4735–4741.[Abstract/Free Full Text]

Janda, J. M. & Abbott, S. L. (1998). Evolving concepts regarding the genus Aeromonas: an expanding panorama of species, disease presentations, and unanswered questions. Clin Infect Dis 27, 332–344.[Medline]

Johnston, R. B., Jr & Kitagawa, S. (1985). Molecular basis for the enhanced respiratory burst of activated macrophages. Fed Proc 44, 2927–2932.[Medline]

Kikuchi, Y. & Graf, J. (2007). Spatial and temporal population dynamics of a naturally-occurring, two-species microbial community inside the digestive-tract of the medicinal leech. Appl Environ Microbiol 73, 1984–1991.[Abstract/Free Full Text]

Klotz, M. G., Klassen, G. R. & Loewen, P. C. (1997). Phylogenetic relationships among prokaryotic and eukaryotic catalases. Mol Biol Evol 14, 951–958.[Abstract]

Loewen, P. C. (1997). Bacterial catalases. In Oxidative Stress and the Molecular Biology of Antioxidant Defenses, pp. 273–308. Edited by J. G. Scandalios. Cold Spring Harbor, NY: Cold Spring Harbor Laboratory.

Lowry, O. H., Rosebrough, N. J., Farr, A. L. & Randall, R. J. (1951). Protein measurement with the folin phenol reagent. J Biol Chem 193, 265–275.[Free Full Text]

Mandell, G. L. (1975). Catalase, superoxide, dismutase, and virulence of Staphylococcus aureus. In vitro and in vivo studies with emphasis on staphylococcal-leukocyte interaction. J Clin Invest 55, 561–566.[Medline]

Morales, V. M., Backman, A. & Bagdasarian, M. (1991). A series of wide-host-range low-copy-number vectors that allow direct screening for recombinants. Gene 97, 39–47.[CrossRef][Medline]

Nelson, K. & Selander, R. K. (1994). Analysis of genetic variation by polymerase chain reaction-based nucleotide sequencing. In Bacterial Pathogenesis, Part A, Identification and Regulation of Virulence Factors, pp. 174–183. Edited by V. L. Clark & P. M. Bavoil. San Diego, CA: Academic Press.

Park, J. B. (2003). Phagocytosis induces superoxide formation and apoptosis in macrophages. Exp Mol Med 35, 325–335.[Medline]

Ramos-Gonzalez, M. I. & Molin, S. (1998). Cloning, sequencing, and phenotypic characterization of the rpoS gene from Pseudomonas putida KT2440. J Bacteriol 180, 3421–3431.[Abstract/Free Full Text]

Rio, R. V., Wu, Y., Filardo, G. & Aksoy, S. (2006). Dynamics of multiple symbiont density regulation during host development: tsetse fly and its microbial flora. Proc Biol Sci 273, 805–814.[CrossRef][Medline]

Ruby, E. G. & McFall-Ngai, M. J. (1999). Oxygen-utilizing reactions and symbiotic colonization of the squid light organ by Vibrio fischeri. Trends Microbiol 7, 414–420.[CrossRef][Medline]

Sambrook, J. & Russell, D. W. (2001). Molecular Cloning: a Laboratory Manual, 3rd edn. Cold Spring Harbor, NY: Cold Spring Harbor Laboratory.

Santos, R., Herouart, D., Sigaud, S., Touati, D. & Puppo, A. (2001). Oxidative burst in alfalfa–Sinorhizobium meliloti symbiotic interaction. Mol Plant Microbe Interact 14, 86–89.[Medline]

Seshadri, R., Joseph, S. W., Chopra, A. K., Sha, J., Shaw, J., Graf, J., Haft, D., Wu, M., Ren, Q. & other authors (2006). Genome sequence of Aeromonas hydrophila ATCC 7966T: the jack of all trades. J Bacteriol 188, 8272–8282.[Abstract/Free Full Text]

Siddall, M. E. & Burreson, E. M. (1998). Phylogeny of leeches (Hirudinea) based on mitochondrial cytochrome c oxidase subunit I. Mol Phylogenet Evol 9, 156–162.[CrossRef][Medline]

Siddall, M. E., Apakupakul, K., Burreson, E. M., Coates, K. A., Erseus, C., Gelder, S. R., Kallersjo, M. & Trapido-Rosenthal, H. (2001). Validating Livanow: molecular data agree that leeches, Branchiobdellidans, and Acanthobdella peledina form a monophyletic group of oligochaetes. Mol Phylogenet Evol 21, 346–351.[CrossRef][Medline]

Simon, R., Priefer, U. & Puhler, A. (1983). A broad host range mobilization system for in vivo genetic engineering: transposon mutagenesis in Gram negative bacteria. Biotechnology 1, 784–791.[CrossRef]

Skorupski, K. & Taylor, R. K. (1996). Positive selection vectors for allelic exchange. Gene 169, 47–52.[CrossRef][Medline]

Soler, L., Yanez, M. A., Chacon, M. R., Aguilera-Arreola, M. G., Catalan, V., Figueras, M. J. & Martinez-Murcia, A. J. (2004). Phylogenetic analysis of the genus Aeromonas based on two housekeeping genes. Int J Syst Evol Microbiol 54, 1511–1519.[Abstract/Free Full Text]

Trontelj, P. & Utevsky, S. Y. (2005). Celebrity with a neglected taxonomy: molecular systematics of the medicinal leech (genus Hirudo). Mol Phylogenet Evol 34, 616–624.[CrossRef][Medline]

Visick, K. L. & Ruby, E. G. (1998). The periplasmic, group III catalase of Vibrio fischeri is required for normal symbiotic competence and is induced both by oxidative stress and by approach to stationary phase. J Bacteriol 180, 2087–2092.[Abstract/Free Full Text]

Worthen, P. L., Gode, C. J. & Graf, J. (2006). Culture-independent characterization of the digestive-tract microbiota of the medicinal leech reveals a tripartite symbiosis. Appl Environ Microbiol 72, 4775–4781.[Abstract/Free Full Text]

Yanez, M. A., Catalan, V., Apraiz, D., Figueras, M. J. & Martinez-Murcia, A. J. (2003). Phylogenetic analysis of members of the genus Aeromonas based on gyrB gene sequences. Int J Syst Evol Microbiol 53, 875–883.[Abstract/Free Full Text]

Yildiz, F. H. & Schoolnik, G. K. (1998). Role of rpoS in stress survival and virulence of Vibrio cholerae. J Bacteriol 180, 773–784.[Abstract/Free Full Text]

Zheng, H. Y., Hassett, D. J., Bean, K. & Cohen, M. S. (1992). Regulation of catalase in Neisseria gonorrhoeae. Effects of oxidant stress and exposure to human neutrophils. J Clin Invest 90, 1000–1006.[Medline]

Received 2 October 2006; revised 31 January 2007; accepted 4 February 2007.


This article has been cited by other articles:


Home page
Appl. Environ. Microbiol.Home page
A. S. Laufer, M. E. Siddall, and J. Graf
Characterization of the Digestive-Tract Microbiota of Hirudo orientalis, a European Medicinal Leech
Appl. Envir. Microbiol., October 1, 2008; 74(19): 6151 - 6154.
[Abstract] [Full Text] [PDF]


Home page
J. Bacteriol.Home page
A. C. Silver, N. M. Rabinowitz, S. Kuffer, and J. Graf
Identification of Aeromonas veronii Genes Required for Colonization of the Medicinal Leech, Hirudo verbana
J. Bacteriol., October 1, 2007; 189(19): 6763 - 6772.
[Abstract] [Full Text] [PDF]


This Article
Right arrow Abstract Freely available
Right arrow Full Text (PDF)
Right arrow Supplementary data
Right arrow Alert me when this article is cited
Right arrow Alert me if a correction is posted
Right arrow Citation Map
Services
Right arrow Email this article to a friend
Right arrow Similar articles in this journal
Right arrow Similar articles in PubMed
Right arrow Alert me to new issues of the journal
Right arrow Download to citation manager
Right arrow reprints & permissions
Citing Articles
Right arrow Citing Articles via HighWire
Right arrow Citing Articles via CrossRef
Right arrow Citing Articles via Google Scholar
Google Scholar
Right arrow Articles by Rio, R. V. M.
Right arrow Articles by Graf, J.
Right arrow Search for Related Content
PubMed
Right arrow PubMed Citation
Right arrow Articles by Rio, R. V. M.
Right arrow Articles by Graf, J.
Agricola
Right arrow Articles by Rio, R. V. M.
Right arrow Articles by Graf, J.