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


     


Microbiology 154 (2008), 1152-1160; DOI  10.1099/mic.0.2007/013953-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 CrossRef
Right arrow Citing Articles via Google Scholar
Google Scholar
Right arrow Articles by Bauer, M. E.
Right arrow Articles by Spinola, S. M.
Right arrow Search for Related Content
PubMed
Right arrow PubMed Citation
Right arrow Articles by Bauer, M. E.
Right arrow Articles by Spinola, S. M.
Agricola
Right arrow Articles by Bauer, M. E.
Right arrow Articles by Spinola, S. M.
Microbiology 154 (2008), 1152-1160; DOI  10.1099/mic.0.2007/013953-0
© 2008 Society for General Microbiology

Identification of Haemophilus ducreyi genes expressed during human infection

Margaret E. Bauer1,{dagger}, Kate R. Fortney2,{dagger}, Alistair Harrison3, Diane M. Janowicz2, Robert S. Munson3,4 and Stanley M. Spinola1,2,5,6

1 Department of Microbiology and Immunology, Indiana University School of Medicine, 635 Barnhill Drive, Room MS420, Indianapolis, IN 46202, USA
2 Department of Medicine, Indiana University School of Medicine, 545 Barnhill Drive, Room EH 435, Indianapolis, IN 46202, USA
3 The Center for Microbial Pathogenesis in The Research Institute at Nationwide Children's Hospital, 700 Children's Drive, Columbus, OH 43205-2696, USA
4 The Center for Microbial Interface Biology and Department of Pediatrics, The Ohio State University, 700 Children's Drive, Columbus, OH 43205-2696, USA
5 Department of Pathology and Laboratory Medicine, Indiana University School of Medicine, Indianapolis, IN 46202, USA
6 Center for Immunobiology, Indiana University School of Medicine, Indianapolis, IN 46202, USA

Correspondence
Margaret E. Bauer
mebauer{at}iupui.edu


    ABSTRACT
 TOP
 ABSTRACT
 INTRODUCTION
 METHODS
 RESULTS AND DISCUSSION
 REFERENCES
 
To identify Haemophilus ducreyi transcripts that are expressed during human infection, we used selective capture of transcribed sequences (SCOTS) with RNA isolated from pustules obtained from three volunteers infected with H. ducreyi, and with RNA isolated from broth-grown bacteria used to infect volunteers. With SCOTS, competitive hybridization of tissue-derived and broth-derived sequences identifies genes that may be preferentially expressed in vivo. Among the three tissue specimens, we identified 531 genes expressed in vivo. Southern blot analysis of 60 genes from each tissue showed that 87 % of the identified genes hybridized better with cDNA derived from tissue specimens than with cDNA derived from broth-grown bacteria. RT-PCR on nine additional pustules confirmed in vivo expression of 10 of 11 selected genes in other volunteers. Of the 531 genes, 139 were identified in at least two volunteers. These 139 genes fell into several functional categories, including biosynthesis and metabolism, regulation, and cellular processes, such as transcription, translation, cell division, DNA replication and repair, and transport. Detection of genes involved in anaerobic and aerobic respiration indicated that H. ducreyi likely encounters both microenvironments within the pustule. Other genes detected suggest an increase in DNA damage and stress in vivo. Genes involved in virulence in other bacterial pathogens and 32 genes encoding hypothetical proteins were identified, and may represent novel virulence factors. We identified three genes, lspA1, lspA2 and tadA, known to be required for virulence in humans. This is the first study to broadly define transcripts expressed by H. ducreyi in humans.


Abbreviations: SCOTS, selective capture of transcribed sequences

{dagger}These authors contributed equally to this work.

A supplementary figure illustrating the SCOTS procedure, and three supplementary tables listing tissue samples used in this study, primers used in this study and genes expressed in samples from at least two volunteers, are available with the online version of this paper.


    INTRODUCTION
 TOP
 ABSTRACT
 INTRODUCTION
 METHODS
 RESULTS AND DISCUSSION
 REFERENCES
 
Haemophilus ducreyi is the causative agent of chancroid, a genital ulcer disease. H. ducreyi facilitates both the acquisition and transmission of human immunodeficiency virus (HIV)-1 and contributes to the HIV-1 pandemic in certain regions of Africa and Asia (Steen, 2001Down).

To study the pathogenesis of H. ducreyi infection, we developed an experimental model of infection in which healthy adult volunteers are infected on the upper arm with 101–102 c.f.u. H. ducreyi. The clinical course of disease in the model accurately mimics natural disease through the pustular stage (Al-Tawfiq et al., 1998Down; Palmer et al., 1998Down; Spinola et al., 2002Down). The histology of experimental pustules and naturally acquired ulcers are nearly identical, and H. ducreyi maintains the same general relationship with host cells during experimental and natural infection (Bauer et al., 2006Down; Bong et al., 2002aDown). Pustules (n=10) that have been quantitatively cultured contain 1.8x105±3.6x105 total c.f.u. (mean±SD) of H. ducreyi, indicating that the bacteria replicate in the model.

Studies of virulence factors in H. ducreyi have relied on identification of gene function in vitro, followed by mutant/parent comparison trials in human or other models of disease. Of 20 putative virulence factors (Janowicz et al., 2006bDown; Spinola et al., 2002Down), nine genes or gene clusters have been shown to be required for full virulence in humans (Al-Tawfiq et al., 2000Down; Bong et al., 2001Down; Fortney et al., 2000Down; Fulcher et al., 2006Down; Janowicz et al., 2004Down, 2006aDown; Spinola et al., 2003bDown; S. M. Spinola and others, unpublished observations). Although these trials yield important information about the roles of specific putative virulence determinants in disease, this approach is impractical for identifying broad categories of genes important to the pathogenesis or survival of the organism. In other bacterial systems, identification of in vivo-expressed genes has led to the discovery of virulence factors that might otherwise not have been identified, because the genes were not expressed in vitro, their annotation did not suggest a role in virulence, or the complexity of their in vivo niche could not be adequately reproduced in vitro (Rediers et al., 2005Down). Thus, we sought a broader approach to examining H. ducreyi gene expression in vivo.

Several molecular techniques have been developed to identify in vivo-expressed bacterial genes (Chiang et al., 1999Down; Mahan et al., 2000Down; Rediers et al., 2005Down; Shelburne & Musser, 2004Down). In selecting a method for our model, we were limited by both practical and regulatory constraints. Biosafety and Food and Drug Administration regulations preclude us from using promoter-trap or mutagenesis strategies that utilize plasmids, mobilizable cassettes or reporter genes in humans. The amount of bacteria in a pustule yields too little RNA for microarray analysis, employing methodologies currently in use. However, H. ducreyi transcripts have been amplified from pustules by RT-PCR (Throm & Spinola, 2001Down). Selective capture of transcribed sequences (SCOTS) was designed to identify scarce bacterial mRNA in the presence of large amounts of eukaryotic RNA (recently reviewed by Daigle et al., 2002Down). With SCOTS, in vivo-derived cDNAs are captured onto bacterial chromosomal DNA by hybridization. Importantly, prehybridization of chromosomal DNA with in vitro-derived cDNAs reduces capture of in vivo-derived sequences that are also transcribed in vitro. SCOTS has been employed successfully to examine gene expression in multiple organisms and conditions: in Salmonella typhi, Mycobacterium tuberculosis and Mycobacterium avium after infection of macrophages (Daigle et al., 2001Down; Graham & Clark-Curtiss, 1999Down; Hou et al., 2002Down); in animal models of Actinobacillus pleuropneumoniae and Escherichia coli (Baltes & Gerlach, 2004Down; Dozois et al., 2003Down); and in Helicobacter pylori in human gastric biopsies (Graham et al., 2002Down). The technique has also been used to examine differential gene expression in Listeria monocytogenes in response to temperature changes (Liu et al., 2002Down) and to compare gene expression between related bacterial species or strains (Dozois et al., 2003Down; Morrow et al., 1999Down).

In the present study, we employed the SCOTS procedure to capture H. ducreyi transcripts that are expressed during experimental human infection. This study has led to the identification of several genes that are required for infection in humans.


    METHODS
 TOP
 ABSTRACT
 INTRODUCTION
 METHODS
 RESULTS AND DISCUSSION
 REFERENCES
 
Tissues.
Pustules for SCOTS were obtained by biopsy from three women who had participated in human-challenge trials (see Supplementary Table S1). Nine additional pustules were obtained from seven men and two women for RT-PCR (Supplementary Table S1). All pustules were obtained 6–9 days after inoculation when the subjects reported pain. Informed consent for participation and for HIV serology was obtained from the volunteers in accordance with the human experimentation guidelines of the US Department of Health and Human Services and the Institutional Review Board of Indiana University-Purdue University at Indianapolis.

RNA preparation.
Whole biopsies were immediately placed in RNAlater (Qiagen). After 30 min, tissues were homogenized in buffer RLT provided in the RNeasy Fibrous Tissue kit (Qiagen) using either a tissue homogenizer or a bead beater (Biospec Products) with 2.4 mm Zirconia beads (Biospec). RNA was extracted with the RNeasy Fibrous Tissue kit following the manufacturer's directions, except that lysozyme (400 µg ml–1) was added to the proteinase K digestion step to lyse the bacteria. In vitro-derived RNA was obtained from a culture of H. ducreyi 35 000 HP (HP, human passaged) (Al-Tawfiq et al., 1998Down) grown in broth and used to inoculate human volunteers, as described previously (Al-Tawfiq et al., 1998Down; Spinola et al., 1994Down, 1996Down). At mid-exponential phase (OD660=0.2), 4 ml bacteria were collected by centrifugation and suspended in 1 ml Ultraspec RNA (Biotecx Laboratories). Total RNA was extracted following the manufacturer's directions. After isolation, all RNA samples were treated with DNase (Ambion) following the manufacturer's directions. Integrity of the RNA samples was determined with an Agilent Bioanalyser (Agilent Technologies).

cDNA synthesis.
RNA was converted to double-stranded cDNA by random priming with the Advantage RT-for-PCR kit (Clontech Laboratories), following the manufacturer's instructions, and second-strand synthesis as described elsewhere (Froussard, 1992Down; Hou et al., 2002Down). The primers had defined 5' sequences followed by nine random 3' bases. Primers with different 5' sequences were used to differentially tag the 5' ends of cDNA derived from tissue (5'-GACACTCTCGAGACATGAGCGGTACC-3') and from broth-grown bacteria (5'-CCTCTGAAGGTTCCAGAATCGATAG-3'). After synthesis, each cDNA pool was amplified by PCR using primers consisting of the appropriate 5' tag without the random 3' sequences. The amplimers were precipitated and suspended in 10 mM HEPES–1 mM EDTA to a final concentration of 375 ng µl–1.

rDNA plasmid construction.
A 5.8 kb locus containing 5S, 16S and 23S rRNA genes was amplified by PCR from 35 000 HP using the following primers: (forward) 5'-TTGCAAAGAGGGAAAGATATCCTATAATGG-3' and (reverse) 5'-TAGTACCCAATTACCTTATTTGTTGTTTGG-3'. The amplicon was cloned into the TA cloning vector pCR-XL-TOPO (Invitrogen) to make pSCOTS1 and propagated in Escherichia coli TOP10F' (Invitrogen).

SCOTS.
The SCOTS procedure (illustrated in Supplementary Fig. S1) (Daigle et al., 2002Down; Faucher et al., 2006Down; Graham & Clark-Curtiss, 1999Down; Graham et al., 2002Down; Liu et al., 2002Down; Morrow et al., 1999Down) was followed as described by Hou et al. (2002)Down. Graham & Clark-Curtiss (1999)Down showed that three rounds of SCOTS was sufficient to effectively eliminate eukaryotic sequences and normalize the relative levels of transcripts in each cDNA pool. Because the M. tuberculosis chromosome contains one set of rRNA genes, while the H. ducreyi chromosome contains six rRNA gene clusters, we increased the molar ratio of blocking rDNA 10-fold from the original protocol, to ensure an excess of blocking rDNA (Graham & Clark-Curtiss, 1999Down). With this modification, a pilot study confirmed that three rounds of SCOTS eliminated detectable rRNA-derived sequences from our cDNA pools (data not shown).

Each round of SCOTS used 1.2 µg biotinylated chromosomal DNA, 66 µg rDNA and 3 µg of the appropriate cDNA pool. In order to obtain sufficient material and to minimize bias introduced by individual PCRs, the first round of SCOTS for each cDNA pool consisted of 10 independent reactions that were pooled after PCR amplification. Two reactions were carried out and pooled in each of the subsequent rounds.

The tissue-derived and broth-derived cDNAs were each subjected separately to three rounds of SCOTS to generate four cDNA pools. The three tissue-derived cDNA pools were individually subjected to three rounds of competitive SCOTS in which the chromosomal DNA was prehybridized with the broth-derived cDNA pool to block capture of genes expressed both in vitro and in vivo. The final pool of selectively captured cDNAs was cloned, and the clones were screened for inserts on Luria–Bertani (LB) agar plates supplemented with ampicillin (50 µg ml–1), X-Gal (40 µg ml–1) and IPTG (20 mM). Individual clones containing inserts were picked into LB broth in 96-well plates. Following overnight culture at 37 °C, the plates were subjected to plasmid isolation and automated sequencing in a 96-well format, using universal primers (M13 forward and reverse), exactly as described previously (Munson et al., 2004Down). Three or four 96-well plates of clones were sequenced from each of the three cDNA libraries. Sequences of cloned inserts were compared with the H. ducreyi genome using the BLASTN algorithm.

Southern analysis.
Inserts from cDNA clones were amplified by PCR using M13 primers, separated electrophoretically and transferred to nitrocellulose by capillary action. Normalized, tissue-derived and broth-derived cDNA pools were labelled using the PCR DIG Probe Synthesis kit (Roche) for use as probes. Blots were hybridized with the labelled probes and developed using the DIG Easy Hyb system (Roche).

RT-PCR.
RNA was extracted from pustules as described above or as described previously (Throm & Spinola, 2001Down). RT-PCR was performed to amplify H. ducreyi genes shown previously to be expressed in vivo (Throm & Spinola, 2001Down) or on the genes listed in Table 1Down. Each PCR reaction was performed alongside controls that received the same amount of RNA template but no reverse transcriptase; all such controls were negative. Each reaction was performed once on each tissue indicated in Table 1Down. Primers are listed in Supplementary Table S2 available with the online version of this paper.


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

 
Table 1. Verification of in vivo gene transcription

 

    RESULTS AND DISCUSSION
 TOP
 ABSTRACT
 INTRODUCTION
 METHODS
 RESULTS AND DISCUSSION
 REFERENCES
 
RNA isolation
RNA was obtained from pustules of three volunteers (subjects 238, 240 and 250) who were experimentally infected with H. ducreyi (Supplementary Table S1). Pustules contain necrotic cells that may release RNases; RNAlater was required to stabilize the specimens. The integrity of RNA in the samples was determined by analysis with an Agilent Bioanalyser, which provides a spectrophotometric tracing and a virtual gel image of the RNA preparation (Fig. 1Down). The eukaryotic RNA was intact (Fig. 1Down), but the prokaryotic RNA levels were too low to visualize in these samples. Amplification of mRNA from several H. ducreyi genes, including cdtB, losB, dsrA and lspA1, by RT-PCR was successful (data not shown), suggesting that the prokaryotic RNA was also intact. Intact RNA was also derived from a broth culture of 35 000 HP that was used to infect volunteers (Fig. 1Down). The four samples were independently subjected to SCOTS to produce four cDNA pools. The in vitro-derived cDNA pool was then used to block capture of each in vivo-derived cDNA pool in competitive SCOTS. Competitive SCOTS theoretically identifies transcripts that may be exclusively expressed or more abundant in one condition than in another. Increased transcript abundance may be due to higher rates of transcription, variations in mRNA degradation, and quality of the RNA samples. Thus, our results were not likely to be biased by differences in the quality of RNA among the samples.


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

 
Fig. 1. Composite Agilent Bioanalyser virtual gel image of RNA samples used in the study. Lanes: 1, RNA ladder with band sizes in bp; 2, representative tissue-derived RNA sample; 3, broth-derived RNA sample. Major eukaryotic and prokaryotic rRNA species are indicated.

 
Sequence analysis of clones obtained after competitive SCOTS
We sequenced inserts of 924 clones obtained after competitive SCOTS. All inserts contained H. ducreyi DNA. Only one insert was derived from H. ducreyi rRNA, indicating a blocking efficiency of >99 %. DNA in eight inserts was derived from tRNAs. A few clones contained sequences that mapped to intergenic regions or pseudogenes; the remaining clones had inserts containing DNA homologous to known or putative ORFs in the H. ducreyi genome.

The 924 clones contained sequences from 531 genes, representing nearly one-third of the genome. We were somewhat surprised by the large number of captured transcripts, as H. ducreyi only infects humans and has no other known environmental niche, and presumably would not require the complex regulation observed in bacteria that infect humans and other species or occupy environmental niches. However, within the human host, H. ducreyi infects a variety of tissues, including stratified squamous or mucosal epithelium and lymph nodes, external sites such as the foreskin and labia, and internal sites such as the vagina and cervix. The cutaneous immune response to the organism is complex and consists of serum components, polymorphonuclear neutrophils, macrophages, T cells and dendritic cells. In addition, H. ducreyi infects men and women, who differ in their susceptibility to infection (Bong et al., 2002bDown; Spinola et al., 2003aDown). The H. ducreyi genome contains at least two intact systems that could sense extracytoplasmic stress; specifically, the Cpx two-component system and the sigma E (RpoE) system. Thus, H. ducreyi likely regulates its gene expression in response to different host pressures and microenvironments, as has been shown for other strict human pathogens such as Bordetella pertussis, which utilizes the bvg two-component system to regulate virulence genes.

Confirmation that the genes identified by SCOTS are expressed in vivo
To confirm that the identified genes were expressed in vivo, 60 clones from each tissue were analysed by Southern blotting, which has been used as a confirmatory test in most other SCOTS-based studies (Baltes & Gerlach, 2004Down; Daigle et al., 2001Down; Graham & Clark-Curtiss, 1999Down; Haydel et al., 2002Down; Hou et al., 2002Down; Morrow et al., 1999Down). Duplicate blots of cloned inserts were probed with normalized cDNA pools from tissues or broth cultures after three rounds of SCOTS, and the signals from the two probes were compared. Representative results are shown in Fig. 2Down. Of 180 probed clones, 87 % hybridized with the normalized tissue-derived H. ducreyi cDNA but either did not hybridize with the normalized broth-derived H. ducreyi cDNA (Fig. 2Down, lanes 1, 2 and 4–6) or showed a much stronger signal with the tissue-derived cDNA than with the broth-derived cDNA (Fig. 2Down, lane 7). These data show that the competitive SCOTS procedure preferentially captured sequences that were present in the in vivo-derived cDNA pool. Of the clones that did not show a stronger signal with the tissue probe, only one (Fig. 2Down, lane 8) showed a similar signal with both probes, while the remainder showed no signal with either probe and were therefore inconclusive (Fig. 2Down, lane 3).


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

 
Fig. 2. Composite Southern blot representative of analysis of sequences identified by SCOTS. Cloned inserts were amplified by PCR, electrophoretically transferred, and probed with broth-derived cDNA (a) or tissue-derived cDNA (b), as described in Methods. Lanes: 1–7, amplified inserts of SCOTS-derived clones; 8, cloned insert expressed in broth and in vivo; 9, unlabelled probe as positive control.

 
To further confirm the SCOTS data, we also performed RT-PCR on RNA isolated from H. ducreyi-infected pustules obtained from nine additional volunteers (Supplementary Table S1). Ten of 11 (91 %) of the genes captured by SCOTS were amplified from one or more additional tissues (Table 1Up). The one gene not detected by RT-PCR was identified in two of the subjects who provided tissue for SCOTS. The data suggest that genes identified by SCOTS are transcribed in multiple volunteers.

We considered but could not use quantitative real-time RT-PCR (qRT-PCR) to confirm our findings. qRT-PCR requires the ability to normalize results from each sample to an objective reference reflecting total bacterial genomes or RNA levels, or to a gene constitutively expressed under the conditions being compared (Shelburne & Musser, 2004Down). Our in vivo RNA is almost entirely eukaryotic, while our in vitro RNA is entirely prokaryotic RNA. RNAlater was required to stabilize the tissue samples, which were processed entirely for RNA. Thus, we have no information about the number of bacteria or bacterial genomes in the samples. Unfortunately, multiple attempts to process biopsies simultaneously for DNA and RNA so that we could quantify bacterial genomes yielded degraded RNA. The identity of genes constitutively expressed by H. ducreyi in vitro and in vivo is unknown. In the absence of this information, 16S rRNA is commonly used to normalize bacterial RNA samples. However, 16S rRNA levels are highly dependent on growth conditions and nutrient availability (Gourse et al., 1996Down). The doubling time of 35 000 HP in broth is approximately 90 min, while the estimated doubling time of 35 000 HP in lesions is 16.5±3.8 h (Throm & Spinola, 2001Down). Thus, we cannot assume that the 16S rRNA content of the bacteria grown in vitro is similar to that of the bacteria in a pustule. Although the chief limitation of SCOTS is that it is qualitative and we did not prove that the captured transcripts were upregulated in vivo, the captured genes must be expressed in vivo.

Redundancy of gene expression among volunteers
Fig. 3Down shows the distribution of the 531 identified genes among the three volunteers. Although the majority of transcripts were identified in cDNA from only one volunteer, 139 transcripts were identified in at least two volunteers, and 28 transcripts were identified as being expressed in all three volunteers (Fig. 3Down). Analysis of the functional categories of the identified genes demonstrated similar patterns of gene expression in all three volunteers (Table 2Down).


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

 
Fig. 3. Venn diagram of the distribution of genes identified by SCOTS in tissues from three volunteers. The number of genes identified in each tissue is indicated; shaded areas show the overlap of genes identified in more than one tissue.

 

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

 
Table 2. Functional categories of in vivo expressed genes

 
With the exception of microarray analysis, none of the published in vivo expression procedures, including SCOTS as well as promoter-trap and mutagenic strategies, is comprehensive for the entire genome (Rediers et al., 2005Down; Shelburne & Musser, 2004Down). Data from SCOTS are limited to the number of clones sequenced. Thus, we do not know whether 139 genes represents the true level of overlap of gene expression from host to host, or whether analysis of additional clones would have yielded greater overlap. RT-PCR of additional tissues showed that many transcripts captured in only one tissue in SCOTS were present in additional volunteers, suggesting that we would have found more overlap had we sequenced more clones.

Analysis of genes expressed in at least two volunteers
Of the 139 genes identified in at least two tissues (see Supplementary Table S3), 68 of 74 (92 %) tested hybridized to a greater extent with tissue-derived than with broth-derived cDNA. Although the RT-PCR data suggested that there may be greater overlap among the three volunteers, observed host effects on H. ducreyi infection indicate that some of the detected genes could be expressed only in a subset of volunteers (Bong et al., 2002bDown; Spinola et al., 2003aDown). We therefore conservatively confined our examination of the potential functional significance of in vivo expressed genes to those from more than one volunteer.

The 139 genes identified in at least two volunteers fell into a variety of functional categories, based on homology with previously characterized bacterial genes (Table 2Up, final column). The largest group was the broad category of biosynthetic and metabolic genes. There were also genes involved in a number of cellular processes, including cell division, DNA replication and repair, transcription, and translation. Some identified genes, including those encoding heat-shock proteins or transporters, can play important roles in responses to stress or nutrient deprivation. Additionally, two regulatory genes were expressed in vivo, including the sensor kinase of a two-component regulator. We also identified 32 genes encoding hypothetical proteins. Homologues of at least 21 of the genes (shown in bold type in Supplementary Table S3) have been identified by promoter-trap strategies as in vivo-expressed during infections caused by other pathogens (Rediers et al., 2005Down). Overall, these data are consistent with in vivo bacterial gene induction analyses, in which many categories of genes are affected in addition to traditional virulence factors (Rediers et al., 2005Down; Shelburne & Musser, 2004Down), and suggest that H. ducreyi responds to the host in a global fashion for in vivo survival.

Six identified electron transport genes are involved in anaerobic respiration in other bacterial systems, including nitrate reductases napA and nrfA, L-lactate dehydrogenase lldD, anaerobic glycerol-3-phosphate dehydrogenase subunits glpB and glpC, and torY, a c-type cytochrome involved in trimethylamine N-oxide reduction. Similarly, two metabolic genes, pflB encoding pyruvate formate-lyase and ansB encoding L-asparaginase II, are homologous to genes that are induced under anaerobic conditions in E. coli (Jennings & Beacham, 1990Down; Sawers & Bock, 1988Down). In contrast, we also identified atpA, an aerobic electron transfer chain subunit, suggesting that some aerobic respiration may also occur. Typically, the environment within pustules is anaerobic (Hays & Mandell, 1974Down). However, H. ducreyi cells are found both deep within the pustule and at the pustule surface, where the environment may be more aerobic (Bauer & Spinola, 2000Down; Bauer et al., 2001Down). Expression of both aerobic and anaerobic respiration genes could reflect differential responses within the bacterial population to microenvironments with different oxygenation levels in vivo.

Six genes identified in this study are involved in amino acid biosynthesis, including argA, argC and argE from the arginine biosynthetic pathway, and carB, which is involved in synthesis of arginine and pyrimidines. Arginine synthesis genes have been identified during in vivo expression studies of several other bacterial pathogens (Rediers et al., 2005Down). Notably, homologues of argA and carAB were identified during in vivo expression studies in Vibrio cholerae and Salmonella enterica serovar Typhimurium, respectively (Camilli & Mekalanos, 1995Down; Mahan et al., 1993Down). Mutational analysis in these pathogens has demonstrated that argA and carAB are required for virulence of V. cholerae and S. enterica, respectively, in mouse models of infection (Camilli & Mekalanos, 1995Down; Mahan et al., 1993Down). These data suggest that in vivo, H. ducreyi, among other pathogens, may require increased arginine, which is a semi-essential amino acid for mammals and thus likely to be in limited supply in tissues. In addition to its utility as an amino acid, arginine is a substrate in other synthesis pathways, such as polyamine production (Cunin et al., 1986Down). Possibly, one or more arginine-requiring biosynthetic pathway(s) is induced in vivo, leading to an increased arginine requirement.

Eight genes identified in vivo are involved in DNA repair, including mutL, mutS, radA, recB, recD, recG, ung and uvrA, suggesting that H. ducreyi encounters DNA-damaging stresses during human infection. Six of the identified DNA repair genes were also identified during in vivo expression studies with other pathogens (Rediers et al., 2005Down). Additionally, the recBCD-encoded exonuclease V is required for repair of double-stranded DNA breaks and has been shown to be necessary in S. enterica for both survival in macrophages and virulence in mice (Cano et al., 2002Down). Thus, DNA repair could be an important mechanism for H. ducreyi survival in vivo.

In addition to genes involved in various cellular processes or regulation, we identified several genes that encode bacterial components that have been characterized by our laboratory and others as potential virulence determinants, including the large supernatant proteins LspA1 and LspA2, the cytolethal distending toxin subunit CdtB, tight adherence protein TadA, and outer-membrane proteins, including MOMP, predicted outer-membrane protein HD1655, predicted outer-membrane lipoprotein HD1094, and a homologue of Haemophilus influenzae outer-membrane protein OMP P1 (Table 2Up and Supplementary Table S3). The LspA proteins LspA1 and LspA2 share 86 % sequence identity and confer antiphagocytic activity to the organism (Vakevainen et al., 2003Down). We identified clones that contained inserts corresponding to identical regions shared by LspA1 and LspA2. Although expression of MOMP and CdtC is not required for pustule formation in humans (Throm et al., 2000Down; Young et al., 2001Down), expression of LspA1 and LspA2 and the TadA-containing operon is required for full virulence in the human model of H. ducreyi infection (Janowicz et al., 2004Down; Spinola et al., 2003bDown). Prompted by the findings from this study, we constructed mutants in HD1844, which encodes WecA, the first enzyme in the enterobacterial common antigen pathway, and in HD0192 (Table 1Up), which encodes a hypothetical lipoprotein. Both mutants were partially attenuated for pustule formation in humans (K. E. Banks and others, unpublished observations; D. M. Janowicz and others, unpublished observations). Thus, several captured transcripts were relevant to pathogenesis. Identification of momp, lspA1 or lspA2, cdtB, and tadA also showed that SCOTS detected some genes expressed by H. ducreyi in vitro (Nika et al., 2002Down; Stevens et al., 1999Down; Throm & Spinola, 2001Down; Ward et al., 1998Down). Southern blot analysis confirmed that momp, lspA2 and tadA were more abundant in the in vivo-derived cDNA pool (data not shown), suggesting that the expression of these genes is upregulated in vivo.

HD1895, a gene encoding a HmwC-like protein, was identified in this screen. In H. influenzae, HmwC is involved in glycosylation of the high-molecular-weight (HMW) adhesins, and a copy of this gene is present in both HMW gene clusters in H. influenzae (Grass et al., 2003Down; St. Geme & Grass, 1998Down). The N-terminal portions of LspA1 and LspA2 share homology with HMW1A and HMW2A (Ward et al., 1998Down), but there is no evidence of glycosylation of the LspA proteins, and the genes up- and downstream of HD1895 are not likely to encode surface-associated proteins. Thus, the substrate for the putative glycosyl transferase HmwC is unknown.

Several genes encoding enzymes involved in lipooligosaccharide (LOS) biosynthesis were identified, including galE, kdkA, rfaE and lpxK. Sialic acid is a component of the LOS of H. ducreyi. HD1669 and HD1670 were identified in the screen and have recently been shown to encode components of an ABC transporter that transports sialic acid (Post et al., 2005Down). HD1842 and HD1844 also encode proteins that are involved in complex carbohydrate biosynthesis. These genes are part of a gene cluster that has strong homology to genes involved in enterobacterial common antigen synthesis.

H. ducreyi has three Mu-like bacteriophage gene clusters. Genes in these regions have not been characterized, and no bacteriophage production or activity has been observed in H. ducreyi. However, seven genes from two of these clusters were identified in our study (Supplementary Table S3).

Our results were similar to those reported in a recent SCOTS-based study of in vivo gene expression by Actinobacillus pleuropneumoniae in a porcine model of infection (Baltes & Gerlach, 2004Down). A total of 46 A. pleuropneumoniae genes expressed in lung specimens pooled from four pigs were identified. The H. ducreyi genome contains homologues of 28 of the 46 A. pleuropneumoniae genes identified in the porcine study, and 12 of the 28 homologues were also expressed in our three human volunteers in vivo. These homologues included a gene encoding an outer-membrane protein, six biosynthetic genes, two regulatory genes, a stress response gene and genes encoding two transporters. While H. ducreyi and A. pleuropneumoniae have different host ranges and cause different diseases, they are members of the Pasteurellaceae, and both are extracellular pathogens that rely on similar strategies such as resistance to phagocytosis for in vivo survival.

In summary, our data suggest that H. ducreyi gene expression undergoes multiple changes within the human host compared with in vitro broth culture. These data also indicate that one or more regulatory mechanisms of H. ducreyi functions in vivo, and define a number of candidates for genes important to in vivo survival or virulence of the organism, including several hypothetical genes as well as genes required for virulence in other bacteria. Future studies will include defining regulatory mechanisms, identifying host factors to which the organism responds and defining the roles in virulence of genes expressed in vivo.


    ACKNOWLEDGEMENTS
 
This work was supported by Public Health Service grants AI27863 and AI31494 from the National Institute of Allergy and Infectious Diseases (NIAID). The human challenge trials were supported by Public Health Service grant MO1RR00750 to the General Clinical Research Center at Indiana University. D. M. J. was supported by NIH grants T32 AI007637 and K08 AI074657 from NIAID. M. E. B. was partially supported by the Indiana Genomics Initiative (INGEN) of Indiana University, which is supported in part by Lilly Endowment. We thank Josephine Clark-Curtiss for training on the SCOTS procedure, many helpful discussions, and critical review of the manuscript. We thank Joan Hou for training and trouble-shooting on the SCOTS procedure, Carisa Townsend for technical assistance, Genshi Zhao, James Graham and France Daigle for helpful discussions, and Tricia Humphreys, Barbara Van Der Pol and Byron Batteiger for critical review of the manuscript.

Edited by: M. P. Stevens


    REFERENCES
 TOP
 ABSTRACT
 INTRODUCTION
 METHODS
 RESULTS AND DISCUSSION
 REFERENCES
 
Al-Tawfiq, J. A., Thornton, A. C., Katz, B. P., Fortney, K. R., Todd, K. D., Hood, A. F. & Spinola, S. M. (1998). Standardization of the experimental model of Haemophilus ducreyi infection in human subjects. J Infect Dis 178, 1684–1687.[CrossRef][Medline]

Al-Tawfiq, J. A., Fortney, K. R., Katz, B. P., Elkins, C. & Spinola, S. M. (2000). An isogenic hemoglobin receptor-deficient mutant of Haemophilus ducreyi is attenuated in the human model of experimental infection. J Infect Dis 181, 1049–1054.[CrossRef][Medline]

Baltes, N. & Gerlach, G. F. (2004). Identification of genes transcribed by Actinobacillus pleuropneumoniae in necrotic porcine lung tissue by using selective capture of transcribed sequences. Infect Immun 72, 6711–6716.[Abstract/Free Full Text]

Bauer, M. E. & Spinola, S. M. (2000). Localization of Haemophilus ducreyi at the pustular stage of disease in the human model of infection. Infect Immun 68, 2309–2314.[Abstract/Free Full Text]

Bauer, M. E., Goheen, M. P., Townsend, C. A. & Spinola, S. M. (2001). Haemophilus ducreyi associates with phagocytes, collagen, and fibrin and remains extracellular throughout infection of human volunteers. Infect Immun 69, 2549–2557.[Abstract/Free Full Text]

Bauer, M. E., Townsend, C. A., Ronald, A. R. & Spinola, S. M. (2006). Localization of Haemophilus ducreyi in naturally acquired chancroidal ulcers. Microbes Infect 8, 2465–2468.[CrossRef][Medline]

Bong, C. T. H., Throm, R. E., Fortney, K. R., Katz, B. P., Hood, A. F., Elkins, C. & Spinola, S. M. (2001). DsrA-deficient mutant of Haemophilus ducreyi is impaired in its ability to infect human volunteers. Infect Immun 69, 1488–1491.[Abstract/Free Full Text]

Bong, C. T. H., Bauer, M. E. & Spinola, S. M. (2002a). Haemophilus ducreyi: clinical features, epidemiology, and prospects for disease control. Microbes Infect 4, 1141–1148.[CrossRef][Medline]

Bong, C. T. H., Harezlak, J., Katz, B. P. & Spinola, S. M. (2002b). Men are more susceptible to pustule formation than women in the experimental model of Haemophilus ducreyi infection. Sex Transm Dis 29, 114–118.[Medline]

Camilli, A. & Mekalanos, J. J. (1995). Use of recombinase gene fusions to identify Vibrio cholerae genes induced during infection. Mol Microbiol 18, 671–683.[CrossRef][Medline]

Cano, D. A., Pucciarelli, M. G., García-del Portillo, F. & Casadesús, J. (2002). Role of the RecBCD recombination pathway in Salmonella virulence. J Bacteriol 184, 592–595.[Abstract/Free Full Text]

Chiang, S. L., Mekalanos, J. J. & Holden, D. W. (1999). In vivo genetic analysis of bacterial virulence. Annu Rev Microbiol 53, 129–154.[CrossRef][Medline]

Cunin, R., Glansdorff, N., Piérard, A. & Stalon, V. (1986). Biosynthesis and metabolism of arginine in bacteria. Microbiol Rev 50, 314–352.[Free Full Text]

Daigle, F., Graham, J. E. & Curtiss, R., III (2001). Identification of Salmonella typhi genes expressed within macrophages by selective capture of transcribed sequences (SCOTS). Mol Microbiol 41, 1211–1222.[CrossRef][Medline]

Daigle, F., Hou, J. Y. & Clark-Curtiss, J. E. (2002). Microbial gene expression elucidated by selective capture of transcribed sequences (SCOTS). Methods Enzymol 358, 108–122.[CrossRef][Medline]

Dozois, C. M., Daigle, F. & Curtiss, R., III (2003). Identification of pathogen-specific and conserved genes expressed in vivo by an avian pathogenic Escherichia coli strain. Proc Natl Acad Sci U S A 100, 247–252.[Abstract/Free Full Text]

Faucher, S. P., Porwollik, S., Dozois, C. M., McClelland, M. & Daigle, F. (2006). Transcriptome of Salmonella enterica serovar Typhi within macrophages revealed through the selective capture of transcribed sequences. Proc Natl Acad Sci U S A 103, 1906–1911.[Abstract/Free Full Text]

Fortney, K. R., Young, R. S., Bauer, M. E., Katz, B. P., Hood, A. F., Munson, R. S., Jr & Spinola, S. M. (2000). Expression of peptidoglycan-associated lipoprotein is required for virulence in the human model of Haemophilus ducreyi infection. Infect Immun 68, 6441–6448.[Abstract/Free Full Text]

Froussard, P. (1992). A random-PCR method (rPCR) to construct whole cDNA library from low amounts of RNA. Nucleic Acids Res 20, 2900[Free Full Text]

Fulcher, R. A., Cole, L. E., Janowicz, D. M., Toffer, K. L., Fortney, K. R., Katz, B. P., Orndorff, P. E., Spinola, S. M. & Kawula, T. H. (2006). Expression of Haemophilus ducreyi collagen binding outer membrane protein NcaA is required for virulence in swine and human challenge models of chancroid. Infect Immun 74, 2651–2658.[Abstract/Free Full Text]

Gourse, R. L., Gaal, T., Bartlett, M. S., Appleman, J. A. & Ross, W. (1996). rRNA transcription and growth rate-dependent regulation of ribosome synthesis in Escherichia coli. Annu Rev Microbiol 50, 645–677.[CrossRef][Medline]

Graham, J. E. & Clark-Curtiss, J. E. (1999). Identification of Mycobacterium tuberculosis RNAs synthesized in response to phagocytosis by human macrophages by selective capture of transcribed sequences (SCOTS). Proc Natl Acad Sci U S A 96, 11554–11559.[Abstract/Free Full Text]

Graham, J. E., Peek, R. M., Jr, Krishna, U. & Cover, T. L. (2002). Global analysis of Helicobacter pylori gene expression in human gastric mucosa. Gastroenterology 123, 1637–1648.[CrossRef]

Grass, S., Buscher, A. Z., Swords, W. E., Apicella, M. A., Barenkamp, S. J., Ozchlewski, N. & St. Geme, J. W., III (2003). The Haemophilus influenzae HMW1 adhesin is glycosylated in a process that requires HMW1C and phosphoglucomutase, an enzyme involved in lipooligosaccharide biosynthesis. Mol Microbiol 48, 737–751.[CrossRef][Medline]

Haydel, S. E., Benjamin, W. H., Jr, Dunlap, N. E. & Clark-Curtiss, J. E. (2002). Expression, autoregulation, and DNA binding properties of the Mycobacterium tuberculosis TrcR response regulator. J Bacteriol 184, 2192–2203.[Abstract/Free Full Text]

Hays, R. C. & Mandell, G. L. (1974). pO2, pH, and redox potential of experimental abscesses. Proc Soc Exp Biol Med 147, 29–30.[Medline]

Hou, J. Y., Graham, J. E. & Clark-Curtiss, J. E. (2002). Mycobacterium avium genes expressed during growth in human macrophages detected by selective capture of transcribed sequences (SCOTS). Infect Immun 70, 3714–3726.[Abstract/Free Full Text]

Janowicz, D. M., Fortney, K. R., Katz, B. P., Latimer, J. L., Deng, K., Hansen, E. J. & Spinola, S. M. (2004). Expression of the LspA1 and LspA2 proteins by Haemophilus ducreyi is required for virulence in human volunteers. Infect Immun 72, 4528–4533.[Abstract/Free Full Text]

Janowicz, D. M., Leduc, I., Fortney, K. R., Katz, B. P., Elkins, C. & Spinola, S. M. (2006a). A DltA mutant of Haemophilus ducreyi is partially attenuated in its ability to cause pustules in human volunteers. Infect Immun 74, 1394–1397.[Abstract/Free Full Text]

Janowicz, D. M., Luke, N. R., Fortney, K. R., Katz, B. P., Campagnari, A. A. & Spinola, S. M. (2006b). Expression of OmpP2A and OmpP2B is not required for pustule formation by Haemophilus ducreyi in human volunteers. Microb Pathog 40, 110–115.[CrossRef][Medline]

Jennings, M. P. & Beacham, I. R. (1990). Analysis of the Escherichia coli gene encoding L-asparaginase II, ansB, and its regulation by cyclic AMP receptor. J Bacteriol 172, 1491–1498.[Abstract/Free Full Text]

Liu, S., Graham, J. E., Bigelow, L., Morse, P. D., II & Wilkinson, B. J. (2002). Identification of Listeria monocytogenes genes expressed in response to growth at low temperature. Appl Environ Microbiol 68, 1697–1705.[Abstract/Free Full Text]

Mahan, M. J., Slauch, J. M. & Mekalanos, J. J. (1993). Selection of bacterial virulence genes that are specifically induced in host tissues. Science 259, 686–688.[Abstract/Free Full Text]

Mahan, M. J., Heithoff, D. M., Sinsheimer, R. L. & Low, D. A. (2000). Assessment of bacterial pathogenesis by analysis of gene expression in the host. Annu Rev Genet 34, 139–164.[CrossRef][Medline]

Morrow, B. J., Graham, J. E. & Curtiss, R., III (1999). Genomic subtractive hybridization and selective capture of transcribed sequences identify a novel Salmonella typhimurium fimbrial operon and putative transcriptional regulator that are absent from the Salmonella typhi genome. Infect Immun 67, 5106–5116.[Abstract/Free Full Text]

Munson, R. S., Jr, Harrison, A., Gillaspy, A., Ray, W. C., Carson, M., Armbruster, D., Gipson, J., Gipson, M., Johnson, L. & other authors (2004). Partial analysis of the genomes of two nontypeable Haemophilus influenzae otitis media isolates. Infect Immun 72, 3002–3010.[Abstract/Free Full Text]

Nika, J. R., Latimer, J. L., Ward, C. K., Blick, R. J., Wagner, N. J., Cope, L. D., Mahairas, G. G., Munson, R. S., Jr & Hansen, E. J. (2002). Haemophilus ducreyi requires the flp gene cluster for microcolony formation in vitro. Infect Immun 70, 2965–2975.[Abstract/Free Full Text]

Palmer, K. L., Schnizlein-Bick, C. T., Orazi, A., John, K., Chen, C.-Y., Hood, A. F. & Spinola, S. M. (1998). The immune response to Haemophilus ducreyi resembles a delayed-type hypersensitivity reaction throughout experimental infection of human subjects. J Infect Dis 178, 1688–1697.[CrossRef][Medline]

Post, D. M., Mungur, R., Gibson, B. W. & Munson, R. S., Jr (2005). Identification of a novel sialic acid transporter in Haemophilus ducreyi. Infect Immun 73, 6727–6735.[Abstract/Free Full Text]

Rediers, H., Rainey, P. B., Vanderleyden, J. & De Mot, R. (2005). Unraveling the secret lives of bacteria: use of in vivo expression technology and differential fluorescence induction promoter traps as tools for exploring niche-specific gene expression. Microbiol Mol Biol Rev 69, 217–261.[Abstract/Free Full Text]

Sawers, G. & Bock, A. (1988). Anaerobic regulation of pyruvate formate-lyase from Escherichia coli K-12. J Bacteriol 170, 5330[Abstract/Free Full Text]

Shelburne, S. A. & Musser, J. M. (2004). Virulence gene expression in vivo. Curr Opin Microbiol 7, 283–289.[CrossRef][Medline]

Spinola, S. M., Wild, L. M., Apicella, M. A., Gaspari, A. A. & Campagnari, A. A. (1994). Experimental human infection with Haemophilus ducreyi. J Infect Dis 169, 1146–1150.[Medline]

Spinola, S. M., Orazi, A., Arno, J. N., Fortney, K., Kotylo, P., Chen, C.-Y., Campagnari, A. A. & Hood, A. F. (1996). Haemophilus ducreyi elicits a cutaneous infiltrate of CD4 cells during experimental human infection. J Infect Dis 173, 394–402.[Medline]

Spinola, S. M., Bauer, M. E. & Munson, R. S., Jr (2002). Immunopathogenesis of Haemophilus ducreyi infection (chancroid). Infect Immun 70, 1667–1676.[Free Full Text]

Spinola, S. M., Bong, C. T. H., Faber, A. L., Fortney, K. R., Bennett, S. L., Townsend, C. A., Zwickl, B. E., Billings, S. D., Humphreys, T. L. & other authors (2003a). Differences in host susceptibility to disease progression in the human challenge model of Haemophilus ducreyi infection. Infect Immun 71, 6658–6663.[Abstract/Free Full Text]

Spinola, S. M., Fortney, K. R., Katz, B. P., Latimer, J. L., Mock, J. R., Vakevainen, M. & Hansen, E. J. (2003b). Haemophilus ducreyi requires an intact flp gene cluster for virulence in humans. Infect Immun 71, 7178–7182.[Abstract/Free Full Text]

St. Geme, J. W., III & Grass, S. (1998). Secretion of the Haemophilus influenzae HMW1 and HMW2 adhesins involves a periplasmic intermediate and requires the HMWB and HMWC proteins. Mol Microbiol 27, 617–630.[CrossRef][Medline]

Steen, R. (2001). On eradicating chancroid. Bull World Health Organ 79, 818–826.[Medline]

Stevens, M. K., Latimer, J. L., Lumbley, S. R., Ward, C. K., Cope, L. D., Lagergård, T. & Hansen, E. J. (1999). Characterization of a Haemophilus ducreyi mutant deficient in expression of cytolethal distending toxin. Infect Immun 67, 3900–3908.[Abstract/Free Full Text]

Throm, R. E. & Spinola, S. M. (2001). Transcription of candidate virulence genes of Haemophilus ducreyi during infection of human volunteers. Infect Immun 69, 1483–1487.[Abstract/Free Full Text]

Throm, R. E., Al-Tawfiq, J. A., Fortney, K. R., Katz, B. P., Hood, A. F., Hansen, E. J. & Spinola, S. M. (2000). Evaluation of an isogenic MOMP-deficient mutant in the human model of Haemophilus ducreyi infection. Infect Immun 68, 2602–2607.[Abstract/Free Full Text]

Vakevainen, M., Greenberg, S. & Hansen, E. J. (2003). Inhibition of phagocytosis by Haemophilus ducreyi requires expression of the LspA1 and LspA2 proteins. Infect Immun 71, 5994–6003.[Abstract/Free Full Text]

Ward, C. K., Lumbley, S. R., Latimer, J. L., Cope, L. D. & Hansen, E. J. (1998). Haemophilus ducreyi secretes a filamentous hemagglutinin-like protein. J Bacteriol 180, 6013–6022.[Abstract/Free Full Text]

Young, R. S., Fortney, K. R., Gelfanova, V., Phillips, C. L., Katz, B. P., Hood, A. F., Latimer, J. L., Munson, R. S., Jr, Hansen, E. J. & Spinola, S. M. (2001). Expression of cytolethal distending toxin and hemolysin are not required for pustule formation by Haemophilus ducreyi in human volunteers. Infect Immun 69, 1938–1942.[Abstract/Free Full Text]

Received 11 October 2007; revised 16 December 2007; accepted 18 December 2007.



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 CrossRef
Right arrow Citing Articles via Google Scholar
Google Scholar
Right arrow Articles by Bauer, M. E.
Right arrow Articles by Spinola, S. M.
Right arrow Search for Related Content
PubMed
Right arrow PubMed Citation
Right arrow Articles by Bauer, M. E.
Right arrow Articles by Spinola, S. M.
Agricola
Right arrow Articles by Bauer, M. E.
Right arrow Articles by Spinola, S. M.


HOME HELP FEEDBACK SUBSCRIPTIONS ARCHIVE SEARCH TABLE OF CONTENTS
INT J SYST EVOL MICROBIOL MICROBIOLOGY J GEN VIROL
J MED MICROBIOL ALL SGM JOURNALS
Copyright © 2008 Society for General Microbiology.