This Article
Right arrow Abstract Freely available
Right arrow Full Text (PDF)
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 E-mail 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 arrowReprints and Permissions
Right arrow Copyright Information
Right arrow Books from ASM Press
Right arrow MicrobeWorld
Citing Articles
Right arrow Citing Articles via Google Scholar
Google Scholar
Right arrow Articles by Qin, J.-H.
Right arrow Articles by Guo, X.-K.
Right arrow Search for Related Content
PubMed
Right arrow PubMed Citation
Right arrow Articles by Qin, J.-H.
Right arrow Articles by Guo, X.-K.

 Previous Article  |  Next Article 

Infection and Immunity, June 2008, p. 2411-2419, Vol. 76, No. 6
0019-9567/08/$08.00+0     doi:10.1128/IAI.01730-07
Copyright © 2008, American Society for Microbiology. All Rights Reserved.

Identification of a Novel Prophage-Like Gene Cluster Actively Expressed in Both Virulent and Avirulent Strains of Leptospira interrogans Serovar Lai {triangledown}

Jin-Hong Qin,1,{dagger} Qing Zhang,2,{dagger} Zhi-Ming Zhang,3 Yi Zhong,1 Yang Yang,1 Bao-Yu Hu,1 Guo-Ping Zhao,2,3* and Xiao-Kui Guo1*

Department of Microbiology and Parasitology, Institutes of Medical Sciences, Shanghai Jiao Tong University, School of Medicine, Shanghai 200025, China,1 Laboratory of Molecular Microbiology, Institute of Plant Physiology and Ecology, Shanghai Institute for Biological Sciences, Chinese Academy of Sciences, Shanghai 200032, China,2 Shanghai-MOST Key Laboratory of Health and Disease Genomics, Chinese National Human Genome Center at Shanghai and National Engineering Center for BioChip at Shanghai, Shanghai 201203, China3

Received 23 December 2007/ Returned for modification 7 January 2008/ Accepted 13 March 2008


arrow
ABSTRACT
 
DNA microarray analysis was used to compare the differential gene expression profiles between Leptospira interrogans serovar Lai type strain 56601 and its corresponding attenuated strain IPAV. A 22-kb genomic island covering a cluster of 34 genes (i.e., genes LA0186 to LA0219) was actively expressed in both strains but concomitantly upregulated in strain 56601 in contrast to that of IPAV. Reverse transcription-PCR assays proved that the gene cluster comprised five transcripts. Gene annotation of this cluster revealed characteristics of a putative prophage-like remnant with at least 8 of 34 sequences encoding prophage-like proteins, of which the LA0195 protein is probably a putative prophage CI-like regulator. The transcription initiation activities of putative promoter-regulatory sequences of transcripts I, II, and III, all proximal to the LA0195 gene, were further analyzed in the Escherichia coli promoter probe vector pKK232-8 by assaying the reporter chloramphenicol acetyltransferase (CAT) activities. The strong promoter activities of both transcripts I and II indicated by the E. coli CAT assay were well correlated with the in vitro sequence-specific binding of the recombinant LA0195 protein to the corresponding promoter probes detected by the electrophoresis mobility shift assay. On the other hand, the promoter activity of transcript III was very low in E. coli and failed to show active binding to the LA0195 protein in vitro. These results suggested that the LA0195 protein is likely involved in the transcription of transcripts I and II. However, the identical complete DNA sequences of this prophage remnant from these two strains strongly suggests that possible regulatory factors or signal transduction systems residing outside of this region within the genome may be responsible for the differential expression profiling in these two strains.


arrow
INTRODUCTION
 
Leptospirosis is a globally important zoonosis caused by infection of pathogenic Leptospira species, including L. alexanderi, L. borgpetersenii, L. interrogans sensu stricto, L. kirschneri, L. noguchii, L. santarosai, L. weilii, L. fainei, L. inadai, and L. meyeri (1, 24). It affects a wide range of mammalian hosts, and human beings and animals become infected through contact with urine-contaminated soil and water. Because of the large spectra of animal species that serve as reservoirs, leptospirosis is considered to be the most widely spread zoonosis.

Prophages are commonly found in sequenced bacterial genomes, many of which appeared to be defective and in a state of mutational decay (15, 18-20). However, both intact and defective prophages may play a special physiological role in the host bacteria and, thus far, they have been implicated in serotype conversion, pathogenesis, and phage immunity (13, 15, 16, 18, 20, 22, 29, 38, 53). Genomic sequences of L. interrogans have been determined for two strains of different serotypes (40, 43), but no intact phage sequence was identified in any one of them.

In the present study, we compared the gene expression profiles between the virulent L. interrogans serogroup Icterohaemorrhagiae serovar Lai type strain 56601 and its corresponding high-passage avirulent strain IPAV under the same in vitro growth conditions and noticed that a fragment of the genomes, ranging from genes LA0186 to LA0219 was differentially expressed. Further analysis showed that this fragment was a defective prophage remnant encoding many phage-related genes. We further identified the five transcripts within this prophage region and characterized the possible regulatory protein encoded by one of the prophage genes (gene LA0195).


arrow
MATERIALS AND METHODS
 
Bacterial strains, plasmids, and culture conditions. The bacterial strains and plasmids used in the present study are listed in Table 1. The virulent strain 56601 and the avirulent variant strain IPAV of L. interrogans serogroup Icterohaemorrhagiae serovar Lai were grown in EMJH liquid medium (33) at 28°C (used for DNA extraction), the optimum growth temperature of pathogenic Leptospira in vitro or 37°C (used for RNA extraction), a temperature similar to that encountered during invasion of mammalian hosts under aerobic conditions. Only mid-log-phase cultures of 100 ml each at a density of ~108 ml–1 were used for gene expression profiling. The cells were harvested by centrifugation at 10,000 x g for 10 min at 4°C.


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

 
TABLE 1. Bacterial strains and plasmids used in this study

Escherichia coli was routinely grown in LB medium at 37°C (46), except where indicated otherwise. Antibiotics, when required, were used at the following concentrations: ampicillin at 100 µg ml–1, kanamycin at 50 µg ml–1, and chloramphenicol at 34 µg ml–1.

DNA manipulation. Genomic DNA of L. interrogans was extracted by using a bacterial DNA minikit (Watsonbiot, China). PCR fragments were purified from agarose gel with a QIAquick gel extraction kit (Qiagen). Restriction endonucleases were purchased from various suppliers and were used according to their specifications. Plasmids DNA from E. coli were prepared by using a plasmid minikit (Tiangen, China).

The genomic draft sequence of the L. interrogans serotype Lai avirulent variant strain IPAV was determined by 454 sequencing (unpublished data). The finished sequence corresponding to the gene cluster analyzed here was generated by direct Sanger sequencing of the PCR-amplified gap-filling fragments.

Bioinformatics analysis. The complete genomic sequence of virulent L. interrogans serovar Lai type strain 56601 is available (http://www.chge.sh.cn/lep/). The phage-like sequences were identified in silico via protein-protein BLAST search (http://www.ncbi.nlm.nih.gov/BLAST/) and were compared to the annotated DNA sequence of the corresponding genomic region of the strain IPVA. Conserved protein structural domains and motifs were searched online (http://www.ncbi.nlm.nih.gov/Structure/cdd/cdd.shtml), and the Pfam domains were also searched online (http://pfam.sanger.ac.uk/). Transcription operons and the corresponding promoters were predicted using Softberry software (Softberry, Mount Kisco, NY). Possible symmetrical repeats located within the putative promoter/operator regions of the perspective transcripts were detected by using an in silico analysis tool (http://www.mgs.bionet.nsc.ru/mgs/programs/OligoRep/InpForm.htm).

RNA isolation. L. interrogans strains were cultured in EMJH medium at 37°C to mid-log phase, and total RNA was extracted from triplicate independent samples. RNA was extracted by using TRIzol reagent (Invitrogen) according to the manufacturer's protocol. Contaminating DNA was digested with RQ1 RNase-free DNase (Promega). The treated RNA was purified by using a Qiagen RNeasy kit, while the quality was verified by agarose gel electrophoresis, and the quantity was determined spectrophotometrically (Eppendorf).

Microarray analysis. The microarray chip was prepared as previously described (42). In brief, the 3,528 protein coding sequences (CDSs) successfully amplified by PCR were printed in triplicate on the slides as probes. Three expression profiling experiments were performed via reverse transcription (RT) using Superscript (Invitrogen); in one experiment the cDNA of IPAV was labeled by Cy3 and the cDNA of 56601 was labeled by Cy5, while in the other two experiments the cDNA of 56601 was labeled by Cy3 and that of IPAV was labeled with Cy5. The unincorporated dye was removed by using a QIAquick nucleotide removal kit (Qiagen) as specified by the manufacturer's protocol. Samples were hybridized at 42°C for 16 h and then washed as described in the manual. Hybridization experiments were performed in triplicate by using cDNA derived from three different cultures of virulent strain 56601 and avirulent strain IPAV grown at 37°C. The hybridized slides were scanned and analyzed by Tiffsplit (Agilent) to calculate the signal intensities and to determine the presence or absence of each CDS. The data were then normalized, and their backgrounds were defined by using GeneSpring 4.0 (Silicon Genetics). The GeneSpring software was used to further analyze the transcription patterns. To identify genes with significantly altered expression levels, cutoff values for expression level ratios 1.5 were used to filter genes with changes (n-fold) greater than ±1.5 in three independent biological samples (27, 48, 50). A Student t test analysis of variance was used to compare the mean expression levels of the test and reference samples. Genes with significant differential expression levels (P < 0.05) were selected.

Real-time PCR. The same RNA samples subjected to microarray analysis were reverse transcribed to cDNA by using Superscript II (Invitrogen) to confirm the changes of expression for selected genes by quantitative PCR (Roche) as described previously (47) with Sybr green dye (Invitrogen). For each amplification run, the CT (threshold cycle) of 16S rRNA gene amplified from the corresponding sample was used as the internal control to normalize the tested gene amplicon for CT calculation (37, 47), and the relative fold changes were further calculated as described previously (37).

Operon transcript analysis. The AMV Reverse Transcription system (Promega) was used for RT experiments according to the manufacturer's instructions. PCR primers designed to detect transcripts are listed in Table 2. To confirm that the RNA sample was free of contaminated DNA, control PCR amplifications with RNA sample plus the primers described above were performed. Control PCR amplifications that produced a PCR fragment indicating the contamination of genomic DNA were excluded from further analysis.


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

 
TABLE 2. Primers used for RT-PCR analysis for transcript identification

Promoter analysis. Several domains of the promoter/operator region for a prospective transcript were cloned into promoter probe vector pKK232-8 (Amersham Pharmacia Biotech, Sweden), upstream of the cat gene. PCR primers for the putative promoters of the perspective transcripts were as follows with BamHI and HindIII restriction sites introduced at the 5' ends of the primers (underlined): transcript I, LA0186-LA0194p-forward (5'-CGGGATCCTCCAACTCAGAATCAGAAAC-3') and LA0186-LA0194p-reverse (5'-CCCAAGCTTCTCAAAAAACTCATCGAAAT-3'); transcript II, LA0195p-forward (5'-CGGGATCCGGTATATGAATGTTGTCCGT-3') and LA0195p-reverse (5'-CCCAAGCTTCATTTTTAACGTATCGTTGT-3'); and transcript III, LA0198-LA0202p-forward (5'-CGGGATCCGATCCATTCATTCGTACCTT-3') and LA0198-LA0202p-reverse (5'-CCCAAGCTTCCTGTTTGGATCGAATACTA-3').

The PCR products and pKK232-8 were digested with BamHI plus HindIII and purified with agarose gel by using a QIAquick gel extraction kit (Qiagen). The PCR products were then ligated into promoter probe vector pKK232-8 individually creating three recombinant putative promoter-cat fusion plasmids designated pKP0194, pKP0195, and pKP0198 (Table 1).

To assay the promoter activities, E. coli XL1-Blue harboring putative promoter-cat fusion plasmid was grown at 37°C to mid-log phase in LB supplemented with ampicillin. The culture was then subcultured at a dilution of 1:100 in 5 ml of LB containing 100 µg of ampicillin ml–1. Chloramphenicol acetyltransferase (CAT) activities were determined by using an enzyme-linked immunosorbent assay test (Roche, Germany) according to the manufacturer's instructions on crude bacterial extracts prepared by sonication. E. coli. XL1-Blue transformed with pKK232-8 without insert was used as a control.

Cloning and overexpression of the LA0195 protein. Based on the DNA sequences of the LA0195 gene, the primers for cloning were designed with BamHI and HindIII restriction sites incorporated (underlined) as follows: LA0195-forward (5'-CGGGATCCGATGTCAAGAATTAATACGA-3) and LA0195-reverse (5'-CCCAAGCTTGGGTTATTTTTTCTTTCTAAA-3').

The PCR products were ligated to the His6 expression vector pET28b (Novagen) under the T7 promoter after digestion with BamHI and HindIII, creating the plasmid pET28b-0195. pET28b-0195 was introduced into E. coli BL21(DE3) competent cells as described previously (21). Induced by IPTG (isopropyl-β-D-thiogalactopyranoside), cloned LA0195 gene was overexpressed in E. coli, and the protein was purified on Ni-NTA agarose (Qiagen) in Tris-HCl buffer eluted by imidazole gradients. The protein was detected by sodium dodecyl sulfate-polyacrylamide gel electrophoresis.

Electrophoretic mobility shift assay (EMSA). Three probes were generated by PCR using 5'-biotin-labeled and unlabeled primers: probe 1, a 323-bp fragment containing a putative promoter of the transcript I sequence (primers 5'-CGGGATCCTCCAACTCAGAATCAGAAAC-3' and 5'-CCCAAGCTTCTCAAAAAACTCATCGAAAT-3'); probe 2, a 279-bp fragment containing a putative promoter of the transcript II sequence (primers 5'-CGGGATCCGGTATATGAATGTTGTCCGT-3' and 5'-CCCAAGCTTCATTTTTAACGTATCGTTGT-3'); and probe 3, a 371-bp fragment containing a putative promoter of the transcript III sequence (primers 5'-CGGGATCCGATCCATTCATTCGTACCTT-3' and 5'-CCCAAGCTTCCTGTTTGGATCGAATACTA-3').

Probes were purified with QIAquick gel extraction kit. An EMSA was performed by using a LightShift Chemiluminescent EMSA kit (Pierce). The purified protein was incubated with probe at 30°C for 20 min in a 20-µl binding reaction containing 1x binding buffer, 5 mM MgCl2, 2.5% glycerol, 0.05% NP-40, 1 µg of poly(dI-dC), and 10 fmol of biotin-labeled probe. A range of concentrations from 0 to 0.2 µg of the LA0195 protein were tested (details are provided in the legend to Fig. 4). Competitor experiments with 50- and 100-fold excesses of unlabeled probe as a specific competitor or poly(dI-dC) as a nonspecific competitor were used to demonstrate the specificity of LA0195 protein binding. Samples were subjected to electrophoresis at 120 V in 1% agarose gel in 0.5x Tris-borate-EDTA for 1.5 h; the gel was then electrophoretically transferred to nylon membrane at 380 mA for 60 min, and cross-linked DNA was transferred to the membrane by using a UV-light cross-linker (UVP, Upland, CA). After the membrane was cross-linked, biotin-labeled DNA was detected by chemiluminescence and exposed to X-ray film for 5 to 10 min.


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

 
FIG. 4. EMSAs for the binding of recombinant protein LA0195 to putative promoter probes representing the upstream regions of transcripts I and II. Each reaction (20 µl) containing 10 fmol of biotin-labeled promoter probe is presented in one lane of the gels. (A) Probe for transcript I with increasing concentrations of protein LA0195, with lanes 1 to 6 corresponding to 0, 0.01, 0.02, 0.05, 0.1, and 0.2 µg, respectively. (B) Probe for transcript I. Lane 1, biotin-labeled probe alone; lanes 2 and 3, biotin-labeled probe with protein LA0195 at 0.1 and 0.2 µg, respectively; lanes 4 and 5, biotin-labeled probe with 0.2 µg of protein LA0195 in the presence of 50- and 100-fold unlabeled probes, respectively; lanes 6 and 7, biotin-labeled probe with 0.2 µg of protein LA0195 in the presence of 50- and 100-fold nonspecific competitors, respectively. (C) Probe for transcript II with increasing concentrations of protein LA0195, with lanes 1 to 6 corresponding to 0, 0.01, 0.02, 0.05, 0.1, and 0.2 µg, respectively. (D) Probe for transcript II. Lanes 1 and 2, biotin-labeled probe with protein LA0195 at 0.1 and 0.2 µg, respectively; lanes 3 and 4, biotin-labeled probe with 0.2 µg of protein LA0195 in the presence of 50- and 100-fold unlabeled probes, respectively; lanes 5 and 6, biotin-labeled probe with 0.2 µg of protein LA0195 in the presence of 50- and 100-fold nonspecific competitor, respectively.


arrow
RESULTS
 
The LA0186 to LA0219 gene cluster of L. interrogans serovar Lai virulent strain 56601 was actively transcribed in its counterpart avirulent strain IPAV. Microarray-based comparative expression profiling of L. interrogans strain 56601 against its avirulent counterpart strain IPAV was performed with samples cultured in EMJH medium at 37°C. These microarray expression profiles have been deposited at the National Center for Biotechnology Information Gene Expression Omnibus database (http://www.ncbi.nlm.nih.gov/geo/info/linking.html), and the corresponding accession number is GSE10788. The results from all of the biological independent triplicates indicated that a gene cluster ranging from genes LA0186 to LA0219 was concomitantly expressed and upregulated in the type strain 56601 compared to its avirulent counterpart strain IPAV (Table 3). However, although the trends of the expressions were the same in three experiments, the fold difference in the third sample was lower than that of the first two samples. Real-time RT-PCR was performed on the third sample, which confirmed the expression profiling data and thus indicated that the discrepancy was not due to the operational error (data not shown).


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

 
TABLE 3. Functional annotation and expression profiling of the prophage-like genomic island (genes LA0186 to LA0219)

The complete DNA sequence of this gene cluster from strain IPAV was determined, which is identical to that from strain 56601 (data not shown). Therefore, the genetic determinants that probably accounted for the expression differences of this gene cluster in these two strains must reside in the genome rather than in this region. Further efforts were thus focused on identification of the gene cluster rather than a mechanistic study of the differential expression.

Gene annotation of the LA0186 to LA0219 gene cluster revealed characteristics of a putative prophage-like genomic island. The active concomitant expression of this gene cluster in both virulent and avirulent strains of L. interrogans serovar Lai led to a comprehensive reannotation of the gene cluster. The average (G+C) content of this gene cluster is 35.1% similar to that of the whole genome (36.0%). BLAST analysis of 34 CDSs within this cluster versus the nonredundant protein database and the conserved domain from the GenBank and Pfam showed that, although the majority of these CDSs (n = 21) are novel genes that failed to exhibit any similarities to proteins of known function in other organisms as previously annotated, in addition to the 5 previously annotated functional CDSs (Table 3), another 8 CDSs also exhibited homology with genes annotated by the National Center for Biotechnology Information.

Of the 13 annotated CDSs, 8 encoded proteins with various degrees of amino acid sequence homology with the ones likely related to phage and prophages (Table 3 and Fig. 1). The LA0205-, LA0209-, LA0211-, and LA0215-encoded proteins exhibit individually different degrees of homology with phage head or tail morphogenesis of different phages. The LA0189 gene was most similar to the gam gene of the Mu phage. The LA0193-encoded protein exhibits low homology with the putative phage integrase family of Pseudoalteromonas haloplanktis TAC125. The putative protein encoded by the LA0213 gene exhibits low homology with transposase, while the LA0195 gene encodes a transcriptional regulator exhibiting some homology with the prophage repressor CI. Based on the presence of phage functional and regulatory homology genes but the lack of plasmid replication or partitioning genes, we hypothesize that this fragment may not encode functions for autonomously extrachromosomal replication ability. Remarkably, at one end of this prophage-like cluster, the LA0219 protein was 70% similar to the LA0218 protein in amino acid sequence, which might be due to rearrangement.


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

 
FIG. 1. Linear presentation of the prophage-like genome island, with its five transcripts tested by RT-PCR. Upper rectangles indicate the transcripts transcribed from right to left, and lower rectangles are transcribed from left to right. Transcripts I to V were illustrated with all of the composed CDSs annotated indicated with three kinds of arrows corresponding to three categories of functional annotation. The eight genes putatively related to the prophages are illustrated (refer to the text for details). The color arrows represent the primers used for RT-PCR identification of the transcripts (see Fig. 2). Forward and reverse primers for a single PCR are labeled with the same color. Refer to the GenBank annotation for transcription directions.

The LA0186-LA0219 gene cluster comprised five transcripts. The genes within this gene cluster have both transcriptional directions shown by genomic annotation. Cotranscription of these genes was analyzed by RT-PCR using L. interrogans RNA as a template (Fig. 2). Primers for this analysis annealed to adjacent or nearby genes, so that amplification in RT-specific reactions would occur only if the genes were cotranscribed (55). The result of RT-PCR assay revealed five transcripts (Fig. 1 and 2). The gene cluster LA0186-LA0194 composed of one transcript was transcribed in a counterclockwise direction and was designated transcript I. The LA0195 gene was transcribed singly in a clockwise direction and was designated transcript II. Transcript III comprised genes LA0196 and LA0197 and was transcribed in an counterclockwise direction. The other two gene transcripts, LA0198-LA0202 and LA0203-LA0219, were transcribed in a clockwise direction and were designated transcripts IV and V, respectively.


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

 
FIG. 2. RT-PCR identification of cotranscribed putative operons within the LA0186-LA0219 gene cluster. RNA was isolated, and RT-PCR was performed as described in the text. Primer pairs spanning two or more of the genes (indicated above each lane) were used to determine whether adjacent or nearby genes are cotranscribed. Molecular size markers in kilobases are indicated. –RT, a negative control, in which the RT step was omitted.

The putative promoters of transcripts I and II are actively expressed in E. coli, whereas very low expression activity was detected for the putative promoter of transcript III. The putative promoter or regulatory regions of transcripts I, II, and III, all proximal to gene LA0195 but with different transcription directions, were identified by using the Softberry prokaryotic promoter prediction tools. The analysis results indicated that the specific promoter-regulatory sequences (–10 and –35 box) were located upstream of transcripts I, II, and III at distances of –146, –144, and –123 bp, respectively, from the translational start site. At the same time, we examined the sequence between transcripts I and II for possible symmetrical repeats by using an in silico analysis tool and found one direct repeat and one inverted repeat sequence in the corresponding regions as shown in Fig. 3.


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

 
FIG. 3. Putative –35 and –10 promoter sequences of transcripts I and II and the putative repeat sequences. Bases underlined in blue were putative –35 and –10 boxes of transcripts I and II. The translational start codon of transcripts I and II are underlined in red, and the translational directions are indicated. The putative repeat sequences are indicated in colors.

Due to the lack of effective genetic tools (25), L. interrogans promoter activity identification was previously tested in E. coli (36). The promoters described above were subcloned into the multicopy vector pKK232-8 (14) in order to generate promoter fusion plasmids pKP0194, pKP0195, and pKP0198 with the promoterless cat reporter gene in a positive orientation. The relative strengths of the promoters were measured by determining the CAT activities in E. coli. The putative promoters of transcripts I and II were detected with high levels of CAT enzyme activities, but very low level of CAT activity was detected for the putative promoter of transcript III (Table 4).


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

 
TABLE 4. E. coli CAT assay to determine the transcription activities of three putative L. interrogans prophage promoters determined in this study

EMSA showed that the LA0195 protein could bind to the putative promoters of transcripts I and II but not to the putative promoter of transcript III. The amino acid sequence of protein LA0195 is similar to that of the prophage repressor CI protein, and it contains a conserved DNA-binding domain. It is known that CI repressor can bind the specific DNA sequence upstream of and proximal to its own structure gene (7, 9, 34). Protein LA0195 was overexpressed and purified. By sodium dodecyl sulfate-polyacrylamide gel electrophoresis, we showed it to have a molecular mass of 21.6 kDa (data not shown).

Transcript I was transcribed head to head with the LA0195 gene. EMSA was performed with 10 fmol of the 5'-biotin-labeled probe 1 corresponding to the 323-bp putative promoter regulatory region of the transcript I by incubation with purified recombinant protein LA0195. As shown in Fig. 4A, a shifted band was observed when 0.02 µg of the LA0195 protein was added to the reaction mixture (Fig. 4A, lane 3), and the intensity of the shifted band became stronger along with the increase of LA0195 protein added up to 0.05 and 0.1 µg (Fig. 4A, lanes 4 and 5). This binding was effectively competed for by 50- and 100-fold excesses of unlabeled probe 1 (Fig. 4B, lanes 4 and 5), whereas the addition of 50- and 100-fold amounts of nonspecific competitor probe of poly(dI-dC) under the same experimental conditions could not prevent LA0195 protein binding (Fig. 4B, lanes 6 and 7). These two control experiments indicated that the binding is sequence specific.

The transcript II contains the LA0195 gene only. As shown in Fig. 4C, a shifted band was observed when 0.05 µg of the LA0195 protein was added to the reaction mixture with probe 2, corresponding to the 279-bp putative promoter of transcript II (Fig. 4C, lane 4). Competition analyses with the unlabeled probe and the nonspecific competitor probe (Fig. 4D, lanes 3 to 6) also indicated that this binding was sequence specific.

EMSA was also performed with the purified recombinant protein LA0195 with probe 3, corresponding to the 371-bp putative promoter of transcript III. No band shift was observed under the same experimental conditions described above (data not shown).


arrow
DISCUSSION
 
Prophages, including defective ones, may play an important role in the evolution of the host bacteria by contributing extra genetic components encoding certain biological properties, including virulence. Defective prophages have frequently been discovered in complete bacterial genomic sequences (35, 39, 51, 52), and many of them harbored some functional genes (6, 38, 49, 51). In the present study, a prophage-like gene cluster in L. interrogans was detected by whole-genome expression profiling and further characterized based on the identification of both cotranscription patterns and a CI-like phage repressor protein.

Active expression of genes covering the region from genes LA0186 to LA0219 was previously observed in whole-genome expression profiling experiments comparing in vitro growth conditions at 28°C versus those at 37°C (42). However, no concomitant up- or downregulations were constantly observed under these conditions. However, a concomitant differential expression in virulent strain 56601 versus the avirulent strains IPAV of serovar Lai was observed, albeit the level of variation was not quantitatively consistent. Because the DNA sequence of this studied region was identical between the two strains, we focused our analysis on the identification of the biological characters of this gene cluster rather than the mechanism that might account for the uncertain differential expression.

Transcription analysis showed that the gene cluster was expressed in five transcripts, i.e., organized into five operons, designated transcripts I to V (Fig. 2). Reannotation of this gene cluster assigned 13 out of the total 34 CDSs to some previously defined proteins. Among them, five CDSs encoding proteins apparently unrelated to prophages, such as the LA0202 gene, which encodes a pore-forming toxin protein that may permeabilize cell membranes by forming small pores (31). It has been reported these toxins are the most common class of bacterial protein toxins and are important for bacterial pathogenesis (2, 5).

Among the eight phage-related genes, LA0195, the only CDS of transcript II, is the only one encoding a putative regulatory protein based on BLAST analysis. Further, Pfam and Motif analysis indicates that LA0195 protein belongs to the HTH_3 family, which contains a helix-turn-helix domain ranging from the 20th to the 75th amino acid residues at its N terminus. It belongs to a large family of DNA-binding proteins that contain 71 members, including a bacterial plasmid copy number control protein, bacterial methylases, various bacteriophage transcription control proteins, and a vegetative specific protein from Dictyostelium discoideum (http://pfam.sanger.ac.uk/family?acc=PF01381). The molecular function of this large family is sequence-specific DNA binding. Remarkably, with respect to the amino acid sequence, the LA0195 protein is similar to the prophage repressor-like CI protein, which is included in this family. For most prophages, the main regulatory switch protein repressor CI usually coordinates the regulation with another protein, Cro. The cro gene is usually closely linked with cI and is transcribed in head-to-head directions to the cI gene, sharing a 14-bp dyad symmetry operator site for CI protein binding (3, 9, 10, 30, 51). However, neither cro nor the 14-bp dyad symmetry operator was found in this gene cluster. In contrast, the LA0195 gene is a singly transcribed gene, indicating that some substantial rearrangement might have taken place in this region since the integration of the prophage.

The putative promoter/operator sequence proximal to the LA0195 gene was further analyzed to test the probable regulatory function related with the putative trans-factor LA0195 gene via in vitro EMSA and in vivo CAT report activity assays in E. coli. The EMSA results showed that protein LA0195 may bind to the strong promoter (tested in vivo) of transcripts I and II, although both of them lack the 14-bp dyad symmetry operator of CI protein binding sites (7, 9, 10, 32, 34, 51). On the other hand, protein LA0195 failed to bind to the putative promoter region of transcript III, and this fact echoed the weak transcription initiation activity of this region tested in E. coli.

The predicted –35 and –10 consensus sequences of the putative promoters for transcripts I and II are both located in the coding regions (Fig. 3). Therefore, the binding of RNA polymerase to the promoter of transcript I located in the CDS of the LA0195 gene might hinder the transcription of transcript II, whose promoter was located in the CDS of the LA0194 gene, the first gene of transcript I. Therefore, one may speculate that the function of regulatory protein LA0195 is to control the transcription switch between transcripts I and II in response to environmental conditions.

Prophage is apparently acquired as easily as they are lost from the chromosome. Prophage genes without selective value to the host are therefore likely to be deleted (18). Bacterial phage particles have been observed in many different spirochetes (4, 8, 17, 23, 26, 28, 29, 41, 44, 45), and most of them can be induced by mitomycin C, ciprofloxacin, or N-methyl-N'-nitro-N-nitrosoguanidine (MNNG) treatment (54). Some of the bacteriophages in spirochetes that have been studied focus their function in the host, such as LE1 of Leptospira biflexa (11, 25); cp32s, a circular plasmid (23); phiBB-1, a bacteriophage of Borrelia burgdorferi (49); and VSH-1, prophage-like gene transfer agent of Brachyspira hyodysenteriae (38). However, no extra DNA content has been found to be part of the virulent L. interrogans genome except for a previously reported genomic island in L. interrogans serovar Lai, which could be excised from a hypothetical prophage (12). Therefore, our report here confirming the gene cluster LA0186-LA0219 as a prophage remnant in the L. interrogans chromosome is the first of its kind.

It was interesting to learn that the differential expression of this prophage-like remnant is unrelated to its DNA sequence. This suggests that complex mechanisms, such as signal transduction or the presence of trans-active coregulators for transcription initiation, may regulate the expression of this remnant corresponding to the environment changes. Further studies on the regulation of leptophage may also contribute to the evolution of the pathogenic mechanism of L. interrogans.


arrow
ACKNOWLEDGMENTS
 
This study was supported in part by the National Natural Science Foundation of China (grants 30370071 and 30670102), the National High Technology Research and Development Program of China (grant 2006AA02Z176), and the Shanghai Leading Academic Discipline Project (T0206).

We are grateful to M. Picardeau (Unite de Bacteriologie Moleculaire et Medicale, Institut Pasteur, Paris, France) and Yu-Feng Yao (Department of Microbiology and Parasitology, Shanghai Jiao Tong University School of Medicine) for thoughtful comments on the manuscript.


arrow
FOOTNOTES
 
* Corresponding author. Mailing address for G.-P. Zhao: Shanghai Institute for Biological Sciences, Chinese Academy of Sciences, Shanghai 200032, China. Phone: 86-21-54924002. Fax: 86-21-50801922. E-mail: gpzhao{at}sibs.ac.cn. Mailing address for X.-K. Guo: Department of Microbiology and Parasitology, Shanghai Jiao Tong University, School of Medicine, Shanghai 200025, China. Phone and fax: (86) 21-64453285. E-mail: microbiology{at}sjtu.edu.cn Back

{triangledown} Published ahead of print on 24 March 2008. Back

Editor: A. J. Bäumler

{dagger} J.-H.Q. and Q.Z. contributed equally to this study. Back


arrow
REFERENCES
 
    1
  1. Anonymous. 2003. Human leptospirosis: guidance for diagnosis, surveillance, and control. World Health Organization, Geneva, Switzerland.
  2. 2
  3. Aroian, R., and F. G. van der Goot. 2007. Pore-forming toxins and cellular non-immune defenses (CNIDs). Curr. Opin. Microbiol. 10:57-61.[CrossRef][Medline]
  4. 3
  5. Aziz, R. K., R. A. Edwards, W. W. Taylor, D. E. Low, A. McGeer, and M. Kotb. 2005. Mosaic prophages with horizontally acquired genes account for the emergence and diversification of the globally disseminated M1T1 clone of Streptococcus pyogenes. J. Bacteriol. 187:3311-3318.[Abstract/Free Full Text]
  6. 4
  7. Barbour, A. G., and S. F. Hayes. 1986. Biology of Borrelia species. Microbiol. Rev. 50:381-400.[Free Full Text]
  8. 5
  9. Basler, M., O. Knapp, J. Masin, R. Fiser, E. Maier, R. Benz, P. Sebo, and R. Osicka. 2007. Segments crucial for membrane translocation and pore-forming activity of Bordetella adenylate cyclase toxin. J. Biol. Chem. 282:12419-12429.[Abstract/Free Full Text]
  10. 6
  11. Beaber, J. W., B. Hochhut, and M. K. Waldor. 2002. Genomic and functional analyses of SXT, an integrating antibiotic resistance gene transfer element derived from Vibrio cholerae. J. Bacteriol. 184:4259-4269.[Abstract/Free Full Text]
  12. 7
  13. Beaber, J. W., and M. K. Waldor. 2004. Identification of operators and promoters that control SXT conjugative transfer. J. Bacteriol. 186:5945-5949.[Abstract/Free Full Text]
  14. 8
  15. Berthiaume, L., Y. Elazhary, R. Alain, and H. W. Ackermann. 1979. Bacteriophage-like particles associated with a spirochete. Can. J. Microbiol. 25:114-116.[Medline]
  16. 9
  17. Blatny, J. M., P. A. Risoen, D. Lillehaug, M. Lunde, and I. F. Nes. 2001. Analysis of a regulator involved in the genetic switch between lysis and lysogeny of the temperate Lactococcus lactis phage phi LC3. Mol. Genet. Genomics 265:189-197.[CrossRef][Medline]
  18. 10
  19. Blatny, J. M., M. Ventura, E. M. Rosenhaven, P. A. Risoen, M. Lunde, H. Brussow, and I. F. Nes. 2003. Transcriptional analysis of the genetic elements involved in the lysogeny/lysis switch in the temperate lactococcal bacteriophage phiLC3, and identification of the Cro-like protein ORF76. Mol. Genet. Genomics 269:487-498.[CrossRef][Medline]
  20. 11
  21. Bourhy, P., L. Frangeul, E. Couve, P. Glaser, I. Saint Girons, and M. Picardeau. 2005. Complete nucleotide sequence of the LE1 prophage from the spirochete Leptospira biflexa and characterization of its replication and partition functions. J. Bacteriol. 187:3931-3940.[Abstract/Free Full Text]
  22. 12
  23. Bourhy, P., L. Salaun, A. Lajus, C. Medigue, C. Boursaux-Eude, and M. Picardeau. 2007. A genomic island of the pathogen Leptospira interrogans serovar Lai can excise from its chromosome. Infect. Immun. 75:677-683.[Abstract/Free Full Text]
  24. 13
  25. Boyd, E. F., and H. Brussow. 2002. Common themes among bacteriophage-encoded virulence factors and diversity among the bacteriophages involved. Trends Microbiol. 10:521-529.[CrossRef][Medline]
  26. 14
  27. Brosius, J. 1984. Plasmid vectors for the selection of promoters. Gene 27:151-160.[CrossRef][Medline]
  28. 15
  29. Brussow, H., C. Canchaya, and W. D. Hardt. 2004. Phages and the evolution of bacterial pathogens: from genomic rearrangements to lysogenic conversion. Microbiol. Mol. Biol. Rev. 68:560-602.[Abstract/Free Full Text]
  30. 16
  31. Brussow, H., and R. W. Hendrix. 2002. Phage genomics: small is beautiful. Cell 108:13-16.[CrossRef][Medline]
  32. 17
  33. Calderaro, A., G. Dettori, L. Collini, P. Ragni, R. Grillo, P. Cattani, G. Fadda, and C. Chezzi. 1998. Bacteriophages induced from weakly beta-hemolytic human intestinal spirochaetes by mitomycin C. J. Basic Microbiol. 38:323-335.[CrossRef][Medline]
  34. 18
  35. Canchaya, C., G. Fournous, and H. Brussow. 2004. The impact of prophages on bacterial chromosomes. Mol. Microbiol. 53:9-18.[CrossRef][Medline]
  36. 19
  37. Canchaya, C., C. Proux, G. Fournous, A. Bruttin, and H. Brussow. 2003. Prophage genomics. Microbiol. Mol. Biol. Rev. 67:238-276.[Abstract/Free Full Text]
  38. 20
  39. Casjens, S. 2003. Prophages and bacterial genomics: what have we learned thus far? Mol. Microbiol. 49:277-300.[CrossRef][Medline]
  40. 21
  41. Chung, C. T., S. L. Niemela, and R. H. Miller. 1989. One-step preparation of competent Escherichia coli: transformation and storage of bacterial cells in the same solution. Proc. Natl. Acad. Sci. USA 86:2172-2175.[Abstract/Free Full Text]
  42. 22
  43. Eggers, C. H., B. J. Kimmel, J. L. Bono, A. F. Elias, P. Rosa, and D. S. Samuels. 2001. Transduction by phiBB-1, a bacteriophage of Borrelia burgdorferi. J. Bacteriol. 183:4771-4778.[Abstract/Free Full Text]
  44. 23
  45. Eggers, C. H., and D. S. Samuels. 1999. Molecular evidence for a new bacteriophage of Borrelia burgdorferi. J. Bacteriol. 181:7308-7313.[Abstract/Free Full Text]
  46. 24
  47. Faine, S., B. Adler, C. Bolin, and P. Perolat. 1999. Leptospira and leptospirosis, 2nd ed. Medisci, Melbourne, Australia.
  48. 25
  49. Girons, I. S., P. Bourhy, C. Ottone, M. Picardeau, D. Yelton, R. W. Hendrix, P. Glaser, and N. Charon. 2000. The LE1 bacteriophage replicates as a plasmid within Leptospira biflexa: construction of an L. biflexa-Escherichia coli shuttle vector. J. Bacteriol. 182:5700-5705.[Abstract/Free Full Text]
  50. 26
  51. Hayes, S. F., W. Burgdorfer, and A. G. Barbour. 1983. Bacteriophage in the Ixodes dammini spirochete, etiological agent of Lyme disease. J. Bacteriol. 154:1436-1439.[Abstract/Free Full Text]
  52. 27
  53. Hughes, T. R., M. J. Marton, A. R. Jones, C. J. Roberts, R. Stoughton, C. D. Armour, H. A. Bennett, E. Coffey, H. Dai, Y. D. He, M. J. Kidd, A. M. King, M. R. Meyer, D. Slade, P. Y. Lum, S. B. Stepaniants, D. D. Shoemaker, D. Gachotte, K. Chakraburtty, J. Simon, M. Bard, and S. H. Friend. 2000. Functional discovery via a compendium of expression profiles. Cell 102:109-126.[CrossRef][Medline]
  54. 28
  55. Humphrey, S. B., T. B. Stanton, and N. S. Jensen. 1995. Mitomycin C induction of bacteriophages from Serpulina hyodysenteriae and Serpulina innocens. FEMS Microbiol. Lett. 134:97-101.[CrossRef][Medline]
  56. 29
  57. Humphrey, S. B., T. B. Stanton, N. S. Jensen, and R. L. Zuerner. 1997. Purification and characterization of VSH-1, a generalized transducing bacteriophage of Serpulina hyodysenteriae. J. Bacteriol. 179:323-329.[Abstract/Free Full Text]
  58. 30
  59. Husson-Kao, C., J. Mengaud, B. Cesselin, D. van Sinderen, L. Benbadis, and M. P. Chapot-Chartier. 2000. The Streptococcus thermophilus autolytic phenotype results from a leaky prophage. Appl. Environ. Microbiol. 66:558-565.[Abstract/Free Full Text]
  60. 31
  61. Jimenez-Juarez, N., C. Munoz-Garay, I. Gomez, G. Saab-Rincon, J. Y. Damian-Almazo, S. S. Gill, M. Soberon, and A. Bravo. 2007. Bacillus thuringiensis Cry1Ab mutants affecting oligomer formation are non-toxic to Manduca sexta larvae. J. Biol. Chem. 282:21222-21229.[Abstract/Free Full Text]
  62. 32
  63. Johansen, A. H., L. Brondsted, and K. Hammer. 2003. Identification of operator sites of the CI repressor of phage TP901-1: evolutionary link to other phages. Virology 311:144-156.[CrossRef][Medline]
  64. 33
  65. Johnson, R. C., and V. G. Harris. 1967. Differentiation of pathogenic and saprophytic leptospires. I. Growth at low temperatures. J. Bacteriol. 94:27-31.[Abstract/Free Full Text]
  66. 34
  67. Kakikawa, M., S. Ohkubo, M. Syama, A. Taketo, and K. I. Kodaira. 2000. The genetic switch for the regulatory pathway of Lactobacillus plantarum phage (phi)g1e: characterization of the promoter P(L), the repressor gene cpg, and the cpg-encoded protein Cpg in Escherichia coli. Gene 242:155-166.[CrossRef][Medline]
  68. 35
  69. Kita, K., H. Kawakami, and H. Tanaka. 2003. Evidence for horizontal transfer of the EcoT38I restriction-modification gene to chromosomal DNA by the P2 phage and diversity of defective P2 prophages in Escherichia coli TH38 strains. J. Bacteriol. 185:2296-2305.[Abstract/Free Full Text]
  70. 36
  71. Lin, M., H. Dan, and Y. Li. 2004. Identification of a second flagellin gene and functional characterization of a {sigma}70-like promoter upstream of a Leptospira borgpetersenii flaB gene. Curr. Microbiol. 48:145-152.[CrossRef][Medline]
  72. 37
  73. Livak, K. J., and T. D. Schmittgen. 2001. Analysis of relative gene expression data using real-time quantitative PCR and the 2{Delta}{Delta}CT method. Methods 25:402-408.[CrossRef][Medline]
  74. 38
  75. Matson, E. G., M. G. Thompson, S. B. Humphrey, R. L. Zuerner, and T. B. Stanton. 2005. Identification of genes of VSH-1, a prophage-like gene transfer agent of Brachyspira hyodysenteriae. J. Bacteriol. 187:5885-5892.[Abstract/Free Full Text]
  76. 39
  77. Mehta, P., S. Casjens, and S. Krishnaswamy. 2004. Analysis of the lambdoid prophage element e14 in the Escherichia coli K-12 genome. BMC Microbiol. 4:4.[CrossRef][Medline]
  78. 40
  79. Nascimento, A. L., S. Verjovski-Almeida, M. A. Van Sluys, C. B. Monteiro-Vitorello, L. E. Camargo, L. A. Digiampietri, R. A. Harstkeerl, P. L. Ho, M. V. Marques, M. C. Oliveira, J. C. Setubal, D. A. Haake, and E. A. Martins. 2004. Genome features of Leptospira interrogans serovar Copenhageni. Braz. J. Med. Biol. Res. 37:459-477.[Medline]
  80. 41
  81. Neubert, U., M. Schaller, E. Januschke, W. Stolz, and H. Schmieger. 1993. Bacteriophages induced by ciprofloxacin in a Borrelia burgdorferi skin isolate. Zentralbl. Bakteriol. 279:307-315.[Medline]
  82. 42
  83. Qin, J. H., Y. Y. Sheng, Z. M. Zhang, P. He, H. B. Yu, Y. Yang, S. G. Liu, G. P. Zhao, and X. K. Guo. 2006. Genome-wide transcriptional analysis of temperature shift in Leptospira interrogans serovar Lai strain 56601. BMC Microbiol. 6:51.[CrossRef][Medline]
  84. 43
  85. Ren, S. X., G. Fu, X. G. Jiang, R. Zeng, Y. G. Miao, H. Xu, Y. X. Zhang, H. Xiong, G. Lu, L. F. Lu, H. Q. Jiang, J. Jia, Y. F. Tu, J. X. Jiang, W. Y. Gu, Y. Q. Zhang, Z. Cai, H. H. Sheng, H. F. Yin, Y. Zhang, G. F. Zhu, M. Wan, H. L. Huang, Z. Qian, S. Y. Wang, W. Ma, Z. J. Yao, Y. Shen, B. Q. Qiang, Q. C. Xia, X. K. Guo, A. Danchin, I. Saint Girons, R. L. Somerville, Y. M. Wen, M. H. Shi, Z. Chen, J. G. Xu, and G. P. Zhao. 2003. Unique physiological and pathogenic features of Leptospira interrogans revealed by whole-genome sequencing. Nature 422:888-893.[CrossRef][Medline]
  86. 44
  87. Ritchie, A. E., I. M. Robinson, L. A. Joens, and J. M. Kinyon. 1978. A bacteriophage for Treponema hyodysenteriae. Vet. Rec. 103:34-35.[Medline]
  88. 45
  89. Saint Girons, I., D. Margarita, P. Amouriaux, and G. Baranton. 1990. First isolation of bacteriophages for a spirochaete: potential genetic tools for Leptospira. Res. Microbiol. 141:1131-1138.[Medline]
  90. 46
  91. Sambrook, J., E. F. Fritsch, and T. Maniatis. 1989. Molecular cloning: a laboratory manual, 2nd ed. Cold Spring Harbor Laboratory Press, Cold Spring Harbor, NY.
  92. 47
  93. Schmittgen, T. D., B. A. Zakrajsek, A. G. Mills, V. Gorn, M. J. Singer, and M. W. Reed. 2000. Quantitative reverse transcription-polymerase chain reaction to study mRNA decay: comparison of endpoint and real-time methods. Anal. Biochem. 285:194-204.[CrossRef][Medline]
  94. 48
  95. Smoot, L. M., J. C. Smoot, M. R. Graham, G. A. Somerville, D. E. Sturdevant, C. A. Migliaccio, G. L. Sylva, and J. M. Musser. 2001. Global differential gene expression in response to growth temperature alteration in group A Streptococcus. Proc. Natl. Acad. Sci. USA 98:10416-10421.[Abstract/Free Full Text]
  96. 49
  97. Stanton, T. B. 2007. Prophage-like gene transfer agents-novel mechanisms of gene exchange for Methanococcus, Desulfovibrio, Brachyspira, and Rhodobacter species. Anaerobe 13:43-49.[CrossRef][Medline]
  98. 50
  99. Stintzi, A. 2003. Gene expression profile of Campylobacter jejuni in response to growth temperature variation. J. Bacteriol. 185:2009-2016.[Abstract/Free Full Text]
  100. 51
  101. Sumby, P., and M. K. Waldor. 2003. Transcription of the toxin genes present within the staphylococcal phage phiSa3ms is intimately linked with the phage's life cycle. J. Bacteriol. 185:6841-6851.[Abstract/Free Full Text]
  102. 52
  103. Tucker, C. P., and M. W. Heuzenroeder. 2004. ST64B is a defective bacteriophage in Salmonella enterica serovar Typhimurium DT64 that encodes a functional immunity region capable of mediating phage-type conversion. Int. J. Med. Microbiol. 294:59-63.[CrossRef][Medline]
  104. 53
  105. Wagner, P. L., and M. K. Waldor. 2002. Bacteriophage control of bacterial virulence. Infect. Immun. 70:3985-3993.[Free Full Text]
  106. 54
  107. Zhang, H., and R. T. Marconi. 2005. Demonstration of cotranscription and 1-methyl-3-nitroso-nitroguanidine induction of a 30-gene operon of Borrelia burgdorferi: evidence that the 32-kilobase circular plasmids are prophages. J. Bacteriol. 187:7985-7995.[Abstract/Free Full Text]
  108. 55
  109. Zuerner, R. L., R. A. Hartskeerl, H. van de Kemp, and A. E. Bal. 2000. Characterization of the Leptospira interrogans S10-spc-alpha operon. FEMS Microbiol. Lett. 182:303-308.[Medline]


Infection and Immunity, June 2008, p. 2411-2419, Vol. 76, No. 6
0019-9567/08/$08.00+0     doi:10.1128/IAI.01730-07
Copyright © 2008, American Society for Microbiology. All Rights Reserved.





This Article
Right arrow Abstract Freely available
Right arrow Full Text (PDF)
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 E-mail 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 arrowReprints and Permissions
Right arrow Copyright Information
Right arrow Books from ASM Press
Right arrow MicrobeWorld
Citing Articles
Right arrow Citing Articles via Google Scholar
Google Scholar
Right arrow Articles by Qin, J.-H.
Right arrow Articles by Guo, X.-K.
Right arrow Search for Related Content
PubMed
Right arrow PubMed Citation
Right arrow Articles by Qin, J.-H.
Right arrow Articles by Guo, X.-K.