Literature DB >> 21151774

Functional analysis of VopF activity required for colonization in Vibrio cholerae.

Vincent C Tam1, Masato Suzuki, Margaret Coughlin, David Saslowsky, Kuntal Biswas, Wayne I Lencer, Shah M Faruque, John J Mekalanos.   

Abstract

Vibrio cholerae, a Gram-negative facultative pathogen, is the etiologic agent for the diarrheal disease cholera. We previously characterized a clinical isolate, AM-19226, that translocates a type III secretion system (T3SS) effector protein with actin-nucleating activity, VopF, into the host cells. From comparative genomic studies, we identified a divergent T3SS island in additional isolates which possess a VopF homolog, VopN. Unlike the VopF-mediated protrusion formation, VopN localizes to stress fiber in host cells similarly to VopL, which is present in the pandemic strain of Vibrio parahaemolyticus. Chimera and yeast two-hybrid studies indicated that the amino-terminal regions of VopF and VopN proteins interact with distinct host cell factors. We determined that AM-19226-infected cells are arrested at S phase of the cell cycle and that VopF/VopN are antiapoptotic factors. To understand how VopF may contribute to the pathogenesis of AM-19226, we examined the effect of VopF in an in vitro polarized-epithelial model and an in vivo adult rabbit diarrheal model. Within the T3SS pathogenicity island is VopE, a homolog of YopE from Yersinia, which has been shown to loosen tight junctions. In polarized intestinal epithelia, VopF and VopE compromised the integrity of tight junctions by inducing cortical actin depolymerization and aberrant localization of the tight-junction protein ZO-1. An assay for pathogenicity in the adult rabbit diarrhea model suggested that these effectors are involved in eliciting the diarrheal response in infected rabbits.

Entities:  

Mesh:

Substances:

Year:  2010        PMID: 21151774      PMCID: PMC2999938          DOI: 10.1128/mBio.00289-10

Source DB:  PubMed          Journal:  mBio            Impact factor:   7.867


INTRODUCTION

Vibrio cholerae, a Gram-negative facultative pathogen, is the causative agent for the diarrheal disease cholera. With more than 200 serogroups, this diverse organism has been able to persist in the aquatic environment while a subpopulation has evolved to cause human disease. The O1 and O139 serogroups have caused epidemics and pandemics using toxin-coregulated pilus (TCP) and cholera toxin (CT). While non-O1, non-O139 strains have caused sporadic disease globally (1–9), the virulence mechanisms are not fully understood, since most of these strains lack TCP and CT. Previously, we identified a type III secretion system (T3SS) in a clinical non-O1, non-O139 strain, AM-19226 (10). Furthermore, a functional T3SS was required for colonization of the small intestine in the infant mouse model (11). A type III effector, vopF, contains one FH1-like domain and three WH2 domains and was characterized as an actin nucleator. Delivery of VopF by infection or transfection induces actin-rich-protrusion formation in host cells (11). Recently, T3SS+ non-O1, non-O139 strains have been isolated in clinical surveys and epidemiology studies of acute diarrheal diseases (12). Additionally, non-O1, non-O139 strains clearly contribute genetically to newly emerging variants of V. cholerae, as exemplified by a horizontal gene transfer event between non-O1 and O1 strains resulting in the evolution of epidemic O139 strains (13, 14). Therefore, it is important to monitor and understand the pathogenic mechanisms of non-O1, non-O139 strains. From the Vibrio genome project, two V. cholerae genomes belonging to the O1 serogroup (N16961 and O395) have been completely sequenced, and 14 additional strains have been sequenced to 5- to 8-fold coverage. We have identified a divergent T3SS island in two non-O1, non-O139 strains, 1587 and 623-39. 1587 is a clinical V. cholerae strain isolated from an outbreak of cholera-like illness in Lima, Peru (2); 623-39 is a diarrheal-disease clinical isolate from Dhaka, Bangladesh (S. M. Faruque, unpublished data). Within the divergent T3SS island, we have identified a vopF homolog, which we named vopN. Interestingly, VopN localizes to the host stress fibers similarly to VopL, a VopF ortholog in Vibrio parahaemolyticus (15). While VopF and VopN have comparable actin nucleation activities in vitro, we determined through the use of chimera studies that the region prior to the FH1-like domain contained the information necessary for distinct localization of the two homologous proteins. A yeast two-hybrid screen identified the distinct interacting host proteins of VopN and VopF. Furthermore, it led us to discover that both proteins have protective effects against apoptosis. Finally, using our in vitro and in vivo models, we identified that both VopE and VopF together are responsible for disrupting cell polarity and loosening tight junctions in the T84 tissue culture model as well as the diarrheal response in the removable intestinal tie-adult rabbit diarrhea (RITARD) model.

RESULTS

Identification of a divergent T3SS pathogenicity island.

Vibrio cholerae is a diverse organism that persists in the aquatic environment. While most non-O1, non-O139 strains thrive exclusively in the environment, some have evolved pathogenic mechanisms to cause human disease independent of TCP and CT. Examination of the sequencing data of eight non-O1, non-O139 V. cholerae strains led to the discovery of a divergent T3SS island in strains 1587 and 623-39. The sequence homologies of the structural components between the T3SS pathogenicity islands in AM-19226 and 1587 are between 63.5% and 75.1% at the DNA sequence level (data not shown). The proteins are much more homologous (between 67 and 89.6% identity and 83.8 and 96.2% similarity) (data not shown). Due to the divergence between the nucleic acid sequences, this T3SS island was previously undetected by PCR, Southern hybridization, or genomic microarray using sequences from the AM-19226 T3SS island. Comparing the organization of the T3SS pathogenicity islands between V. cholerae AM-19226 and 1587 and V. parahaemolyticus T3SS2 (see Fig. S1 and S2 and Table S1 in the supplemental material), we have identified a central conserved region in which most of the structural components reside and two peripheral variable regions where the most diversity occurs and, coincidentally, where the putative effectors reside. The T3SS pathogenicity island identified in 1587 and 623-39 contains a divergent vopF homolog, which we named vopN (37% identity and 55% similarity at the amino acid sequence level). Amino acid sequence alignment between VopF and VopN from V. cholerae and VopL from V. parahaemolyticus (Fig. S3) indicated that VopN shares characteristics with both VopF and VopL.

VopN is an actin nucleator that localizes to the stress fibers.

The VopF ortholog in V. parahaemolyticus, VopL, has been shown to nucleate actin in vitro while forming stress fibers in infected or transfected cells (15). To determine whether VopN induces protrusion formation like VopF or stress fiber formation like VopL, we transfected CHO cells with a vopN::gfp fusion construct. Actin was visualized in fixed cells by staining with Alexa Fluor 588-conjugated phalloidin, and we observed that VopN-green fluorescent protein (GFP) localized to actin stress fibers (Fig. 1A).
FIG 1

VopN is an actin nucleator that localizes to the actin stress fibers. (A) CHO cells were transfected with vopF::pAcGFP-N1 or vopN::pAcGFP-N1 to express VopF-GFP or VopN-GFP fusion proteins, respectively. Overlay images of Texas Red phalloidin (red) and VopF-GFP or VopN-GFP (green) are shown. Bar = 10 µm. (B) Immunogold electron micrographs of vopF::gfp- or vopN::gfp-transfected CHO cells. The two right images are higher-magnification images of the boxed areas. Red arrows point to gold labeling. Bar = 2 µm (left images) or 500 nm (right images). (C) Pyrene actin assay of VopF and VopN. VopF or VopN protein (100, 10, or 1 nM) was incubated with 2 µM monomeric pyrene-labeled actin. Actin polymerization was measured as arbitrary fluorescence intensity (in arbitrary units [a.u.]) over time (seconds), with excitation and emission wavelengths of 365 and 407 nm, respectively. A mixture of actin and pyrene actin, labeled as actin alone, served as a negative control.

VopN is an actin nucleator that localizes to the actin stress fibers. (A) CHO cells were transfected with vopF::pAcGFP-N1 or vopN::pAcGFP-N1 to express VopF-GFP or VopN-GFP fusion proteins, respectively. Overlay images of Texas Red phalloidin (red) and VopF-GFP or VopN-GFP (green) are shown. Bar = 10 µm. (B) Immunogold electron micrographs of vopF::gfp- or vopN::gfp-transfected CHO cells. The two right images are higher-magnification images of the boxed areas. Red arrows point to gold labeling. Bar = 2 µm (left images) or 500 nm (right images). (C) Pyrene actin assay of VopF and VopN. VopF or VopN protein (100, 10, or 1 nM) was incubated with 2 µM monomeric pyrene-labeled actin. Actin polymerization was measured as arbitrary fluorescence intensity (in arbitrary units [a.u.]) over time (seconds), with excitation and emission wavelengths of 365 and 407 nm, respectively. A mixture of actin and pyrene actin, labeled as actin alone, served as a negative control. To gain better insight into the nature of the cellular structures that VopN and VopF induce in host cells, we conducted immunogold electron microscopy. Anti-GFP antibodies and 10-nm gold particle-conjugated antibodies were used to label VopF-GFP or VopN-GFP. vopF::gfp-transfected cells showed gold labeling enrichment along unbranched actin-rich protrusions. vopN::gfp-transfected cells, on the other hand, showed gold labeling enrichment along bundled actin cables at the cell dorsal surface (Fig. 1B). To confirm that VopN can nucleate actin in vitro, we determined the activity of VopN in the pyrene actin assay. Fast-performance liquid chromatography (FPLC)-purified VopN-His6 protein has actin nucleation and polymerization activity in vitro similar to that of VopF (Fig. 1C).

Chimera studies for determining the domain required for inducing different cellular phenotypes.

The discovery of VopF and VopN presents an interesting cell-biological question of how two homologous proteins are localized to different cellular structures. We constructed VopN/VopF hybrids in which we fused the amino-terminal half of VopF or VopN to the carboxy-terminal half of VopN or VopF, respectively. Transfection experiments with the hybrid proteins showed that VopF/N are enriched at protrusions similarly to VopF-transfected cells while the VopN/F hybrid protein localizes to stress fiber similarly to VopN-transfected cells (Fig. 2A). These results suggest that the amino-terminal halves of these type III effectors play an important role in determining protein subcellular localization. To systematically narrow down the domain or region required for stress fiber or protrusion localization, we conducted chimera studies where we fused different fragments of VopN/F to the remaining corresponding fragments of VopF/N. Transfection of these gfp fusion constructs led us to conclude that none of the actin-modulating domains (the FH1-like domain and three WH2 domains) contain information to differentiate the properties of the two proteins. Rather, the regions prior to the FH1-like domains are sufficient to determine the distinct cellular structure localization in host cells (Fig. 2B). To confirm that amino-terminal amino acids (aa) 1 to 84 of VopN and 1 to 135 of VopF are sufficient to differentiate the proteins from each other, we fused those fragments directly to GFP and determined their localization in transfected cells. While VopF1-135-GFP showed uniformly distributed punctate localization across the cells, VopN1-84-GFP localized to the actin stress fibers (Fig. 2C).
FIG 2

Domain analysis of VopF and VopN. (A) Transfection of hybrid genes comprising vopF and vopN into CHO cells. Overlay images of GFP fusion proteins (green) and actin (red) are displayed. (B) Transfection of chimera genes comprising vopF and vopN into CHO cells. Overlay images of GFP fusion proteins (green) and actin (red) are displayed. Micrographs on the left are higher-magnification images of the boxed areas. Since all four chimera genes (vopF and -N and vopN and -F) showed the same phenotype, only the bottom set of chimera genes used to generate the fluorescent micrographs is shown (asterisks). (C) Transfection of VopF1-135-GFP and VopN1-85-GFP into CHO cells. Overlay images of nuclei (cyan), GFP fusion proteins (green), and actin (red) are displayed. Bar = 10 µm.

Domain analysis of VopF and VopN. (A) Transfection of hybrid genes comprising vopF and vopN into CHO cells. Overlay images of GFP fusion proteins (green) and actin (red) are displayed. (B) Transfection of chimera genes comprising vopF and vopN into CHO cells. Overlay images of GFP fusion proteins (green) and actin (red) are displayed. Micrographs on the left are higher-magnification images of the boxed areas. Since all four chimera genes (vopF and -N and vopN and -F) showed the same phenotype, only the bottom set of chimera genes used to generate the fluorescent micrographs is shown (asterisks). (C) Transfection of VopF1-135-GFP and VopN1-85-GFP into CHO cells. Overlay images of nuclei (cyan), GFP fusion proteins (green), and actin (red) are displayed. Bar = 10 µm.

VopF and VopN interact with different host proteins.

The distinct localization patterns of VopF and VopN led us to hypothesize that the two proteins interact with different host proteins. We conducted a yeast two-hybrid screen using a HeLa cell cDNA library and determined that VopN interacts with filamin A and filamin B (Fig. 3A). Filamin is a family of actin binding proteins that cross-link actin and localize to the stress fiber in tissue culture cells. Additionally, filamin acts as a scaffolding protein for signaling proteins and bridges transmembrane proteins to the actin cytoskeleton (16). Examination of the VopN truncations indicated that the portion comprising aa 1 to 84 of the VopN protein could still interact with filamin while the portion comprising aa 1 to 50 of the VopN protein could not (Fig. 3A and data not shown). These data were further corroborated by transfection and colocalization experiments with filamin (Fig. 3D). Conducting the yeast two-hybrid screen using full-length vopF as the bait resulted in no significant hits other than actin. As expected, full-length vopF and vopN interacted with actin while the truncated fusions that did not contain any actin-modulating domains did not. When we conducted the yeast two-hybrid screen using the truncated vopF gene encoding the first 135 aa (without the actin-modulating domains), we discovered interactions with the cyclophilin G (PPIG), breast cancer-associated gene 3 (BCA3/AKIP), and copper metabolism (Murr1) domain-containing 1 (COMMD1) proteins (Fig. 3A). Interactions between the bacterial effectors and the host proteins were confirmed using coimmunoprecipitation (Fig. 3B). PPIG is a nuclear protein that is cell cycle regulated (17). It demonstrates distinct nuclear/cytoplasmic distribution patterns during the cell cycle (17). BCA3 and COMMD1 have been shown to repress NF-κB-dependent transcription (18, 19), and COMMD1 further accelerates the ubiquitination and degradation of NF-κB subunits (20). Indeed, VopF and aa 1 to 135 of VopF, but not VopN, had the activity to inhibit tumor necrosis factor alpha (TNF-α)-induced NF-κB activation (Fig. 3C), suggesting that the binding of VopF to these nuclear proteins in the cytoplasm could exclude them from the nucleus and deregulate NF-κB activity and NF-κB-regulated normal cell functions, such as cell cycle and cell death (21, 22).
FIG 3

VopF and VopN interact with distinct host proteins. (A) A yeast two-hybrid screen identified filamin as the cellular interactor for VopN and PPIG, COMMD1, and BCA3 as the cellular interactors for VopF. The Y184 Saccharomyces cerevisiae host strain carrying various preys fused to the Gal4 activation domain was mated with the AH109 strain carrying full-length vopF, the region encoding aa 1 to 135 of vopF, full-length vopN, or the region encoding aa 1 to 84 of vopN fused to the Gal4 DNA binding domain. The resulting diploids were plated and restreaked onto SD–Ade–His–Leu–Trp–X-α-Gal to select for colonies strongly expressing the Ade2, His3, and Mel1 reporters. (B) Immunoprecipitation results between VopF1-135-GFP and VopN1-84-GFP for various preys fused to the pProLabel tag. The presence of the prey proteins was assayed by chemiluminescence in the lysates (red) and anti-GFP-immunoprecipitated fractions (blue). (C) Luciferase assays of 293T cells cotransfected with the NF-κB reporter and the indicated GFP vectors (250 ng) and stimulated with or without TNF-α (100 ng/ml). (D) Colocalization of filamin and various GFP fusion proteins. Various GFP fusion constructs were transfected into CHO cells. Filamin was stained using antifilamin antibodies. Overlay images of nuclei (cyan), GFP fusion proteins (green), and filamin (red) are displayed. Bar = 10 µm.

VopF and VopN interact with distinct host proteins. (A) A yeast two-hybrid screen identified filamin as the cellular interactor for VopN and PPIG, COMMD1, and BCA3 as the cellular interactors for VopF. The Y184 Saccharomyces cerevisiae host strain carrying various preys fused to the Gal4 activation domain was mated with the AH109 strain carrying full-length vopF, the region encoding aa 1 to 135 of vopF, full-length vopN, or the region encoding aa 1 to 84 of vopN fused to the Gal4 DNA binding domain. The resulting diploids were plated and restreaked onto SD–AdeHisLeuTrp–X-α-Gal to select for colonies strongly expressing the Ade2, His3, and Mel1 reporters. (B) Immunoprecipitation results between VopF1-135-GFP and VopN1-84-GFP for various preys fused to the pProLabel tag. The presence of the prey proteins was assayed by chemiluminescence in the lysates (red) and anti-GFP-immunoprecipitated fractions (blue). (C) Luciferase assays of 293T cells cotransfected with the NF-κB reporter and the indicated GFP vectors (250 ng) and stimulated with or without TNF-α (100 ng/ml). (D) Colocalization of filamin and various GFP fusion proteins. Various GFP fusion constructs were transfected into CHO cells. Filamin was stained using antifilamin antibodies. Overlay images of nuclei (cyan), GFP fusion proteins (green), and filamin (red) are displayed. Bar = 10 µm.

AM-19226 affects cell cycle progression in vitro.

To investigate whether AM-19226 affects host cell responses and whether VopF plays a role in this context, we infected HeLa cells with V. cholerae for 4 h and then incubated the cells in medium with gentamicin for 72 h to assess the contribution of translocated effectors alone. We stained the fixed cells with propidium iodide (PI) to determine the cell cycle distribution of these cells (Fig. 4). While uninfected cells and type III mutant-infected cells have a normal cell cycle distribution, with the majority of cells in G1 phase, cells infected with VTV18 (AM-19226Δendo Δrtx Δhly Δhap) (11) have drastically shifted from G1 to S phase. An increase in the population of apoptotic cells (17%) was also observed in these cells. vopF mutant-infected cells contained an even higher proportion of apoptotic cells (31%) and almost evenly distributed populations of G1, S, and G2/M cells. Apoptotic cells were also visualized using terminal deoxynucleotidyltransferase-mediated dUTP-biotin nick end labeling (TUNEL) staining (see Fig. S4 in the supplemental material). The strain harboring a deletion of all actin-modulating domains of vopF (the VTV18 ΔvopF Dm2 mutant) showed a distribution similar to that of VTV18, suggesting that the actin nucleation domain of VopF is not involved in the antiapoptotic activity of the protein. Interestingly, replacing vopF with vopN resulted in a distribution similar to that of the type III mutant. These data suggest that VopF and VopN are both capable of reducing the type III-dependent apoptosis. Visualization of the infected cells suggested that the cells arrested at S phase have an enlarged nucleus and cell size indicative of cell cycle arrest. The accumulation of apoptotic cells instead of the restoration of normal cell cycle distribution that resulted from infection with the vopF deletion mutant suggests that there may be an unidentified cyclomodulin responsible for the cell cycle arrest.
FIG 4

Cell cycle analysis. HeLa cells were infected with V. cholerae at an MOI of 100 for 4 h. Cells were then incubated in gentamicin for 72 h before being fixed in 75% ethanol and stained with propidium iodide. Cells were analyzed by flow cytometry. Cells were also stained for actin (red) and nuclei (cyan) for microscopy analysis. Bar = 50 µm.

Cell cycle analysis. HeLa cells were infected with V. cholerae at an MOI of 100 for 4 h. Cells were then incubated in gentamicin for 72 h before being fixed in 75% ethanol and stained with propidium iodide. Cells were analyzed by flow cytometry. Cells were also stained for actin (red) and nuclei (cyan) for microscopy analysis. Bar = 50 µm.

Role of VopE and VopF in epithelial integrity and diarrheal response.

As induction of cell cycle arrest by bacterial pathogens has been hypothesized to alter epithelial integrity due to the lack of cell renewal and to prolong bacterial colonization through inhibition of cell shedding (23), we investigated the role of VopF and VopN in the context of our animal models. While we previously found that deletion of vopF results in a 1.5- to 2-log defect in a suckling mouse competition assay, vopN deletions in its native strains, 1587 and 1582 (another isolate from the same outbreak), surprisingly showed no defect in colonization when these strains competed with wild-type strains (Fig. 5A). To understand whether actin nucleation can complement a vopF deletion, we tested the ability of vopN to complement the vopF deletion in AM-19226. Addition of vopN into a vopF-deleted strain partially (1-log) complemented vopF, suggesting that while the actin nucleation activity of VopN and VopF is important, it may not be sufficient to promote efficient colonization in vivo (Fig. 5A). While both vopN and vopF are complemented by homologous recombination into the vopF locus, the level of expression and, more importantly, the translocation of VopF and VopN proteins may differ, which provides an alternate explanation for the phenotype. To understand the role of VopF and VopN in pathogenesis, they must be studied in the context of the native strains because of their unique combination of translocated effectors. Since a VopN deletion mutant of 1587 has no phenotypic defect in colonization, probably due to redundancy of the function of additional effectors, we focused on the role of VopF in its native strain, AM-19226.
FIG 5

Assay for pathogenicity in rabbits. (A) Competitive index results from the infant mouse colonization model. AMD, 1587, and 1582 show the AM-19226, 1587, and 1582 strains of V. cholerae, respectively. The competitive index was calculated as the ratio of the output of the test strain to the output of the reference strain over the ratio of the input of the test strain to the input of the reference strain. An asterisk denotes that the differences between the populations are statistically significant as determined by the Mann-Whitney U test (P < 0.05). (B) Percent decrease in TER after infection with V. cholerae. VTV18, a Δendo Δrtx Δhly Δhap mutant of AM-19226, was used as the wild-type strain. Values were derived by dividing the difference between the recorded resistance observed at the beginning and that observed at the end of the 180-min period by the initial resistance value. An asterisk denotes that the differences between the populations are statistically significant as determined by the Mann-Whitney U test (P < 0.05). (C) Confocal micrographs of T84 cells infected with V. cholerae strains. Overlay images of bacteria (cyan), actin (green), and ZO-1 (red) are shown. The top panel of each micrograph displays a cross section of the x-z plane (apical toward the top and basolateral toward the bottom), and the left panel displays a cross section of the y-z plane, while the main panel displays the x-y view onward to the apical side. (D) Diarrheal response of adult rabbits in the RITARD model.

Assay for pathogenicity in rabbits. (A) Competitive index results from the infant mouse colonization model. AMD, 1587, and 1582 show the AM-19226, 1587, and 1582 strains of V. cholerae, respectively. The competitive index was calculated as the ratio of the output of the test strain to the output of the reference strain over the ratio of the input of the test strain to the input of the reference strain. An asterisk denotes that the differences between the populations are statistically significant as determined by the Mann-Whitney U test (P < 0.05). (B) Percent decrease in TER after infection with V. cholerae. VTV18, a Δendo Δrtx Δhly Δhap mutant of AM-19226, was used as the wild-type strain. Values were derived by dividing the difference between the recorded resistance observed at the beginning and that observed at the end of the 180-min period by the initial resistance value. An asterisk denotes that the differences between the populations are statistically significant as determined by the Mann-Whitney U test (P < 0.05). (C) Confocal micrographs of T84 cells infected with V. cholerae strains. Overlay images of bacteria (cyan), actin (green), and ZO-1 (red) are shown. The top panel of each micrograph displays a cross section of the x-z plane (apical toward the top and basolateral toward the bottom), and the left panel displays a cross section of the y-z plane, while the main panel displays the x-y view onward to the apical side. (D) Diarrheal response of adult rabbits in the RITARD model. Since V. cholerae colonizes the single-cell layer of the polarized small-intestinal epithelium, we wanted to study the interaction of AM-19226 with host cells in the more biologically relevant polarized-T84-tissue-culture-cell model. After allowing the T84 cells to differentiate in Transwell inserts for 7 days, we infected the cells with various bacterial strains and assessed the changes in transepithelial resting potential (Eo) and transepithelial resistance (TER) to determine the epithelial integrity. In the polarized-T84-tissue-culture-cell model, a decrease of TER in the absence of increased Eo is indicative of the loosening of tight junctions. T84 cells infected with the wild-type strain VTV18 showed a 49% decrease in resistance (Fig. 5B). Cells infected with the type III structural mutant (the ΔvcsN2 mutant) were not affected (3.4% increase), showing a level comparable to that of the medium-alone control (0.7% increase). These data suggested that the decrease in TER is type III dependent. For O1 V. cholerae-infected cells, a similar decrease in TER is observed, which is dependent on cholera toxin (data not shown) (24). The T3SS island in AM-19226 carries another putative effector, which we named vopE because of its homology to the T3SS effector yopE from Yersinia species. VopE, like YopE, has a highly conserved GTPase-activating protein (GAP) domain. The GAP activity of YopE has been shown to inactivate RhoA, Rac1, and Cdc42 to depolymerize actin (25). VopE and VopF seem to have opposing functions, where one depolymerizes actin and the other nucleates actin. However, a pathogen possessing effectors with seemingly opposite functions has been observed previously in the Salmonella T3SS, where SptP is a GAP and SopE is a GTPase-exchanging factor (GEF) (26). The combination of SptP and SopE helps Salmonella to invade the host cells and subsequently recover the actin cytoskeleton. Deleting vopE has a minimal effect on colonization in the infant mouse model (data not shown). To compare the opposing functions of VopE and VopF in relation to epithelial integrity, we tested single and double mutants in the polarized-T84-cell system. Compared to the strongly decreased TER caused by VTV18, the vopE and vopF deletion mutants showed intermediate phenotypes (34.6% and 30.2%, respectively). The vopF vopE double mutant (24.6%) showed a more dramatic difference from VTV18. These data indicate that while vopE and vopF are involved in loosening of tight junctions, there are likely to be additional effectors involved. To study the actin cytoskeletal organization and disruption in cell polarity in infected cells, we fixed and stained the Transwell inserts for actin, ZO-1, and AM-19226 (11) (Fig. 2C). The medium-alone control and ΔvcsN2 mutant-infected cells showed a normal localization of ZO-1, a tight-junction protein. The x-z and y-z plane micrographs depicted that ZO-1 was localized to the apical surfaces of the polarized cells. A normal cortical actin network at the apical surface and the perijunctional actin ring was also observed. Cells infected with VTV18 showed massive depolymerization of actin, where minimal cortical actin was observed along with redistribution of ZO-1 toward the basolateral surface, consistent with loss of cell polarity and intercellular tight junctions. ΔvopF and ΔvopE mutant-infected cells also showed a certain degree of depolarization based on ZO-1 mislocalization but noticeably more-abundant cortical actin than wild-type-infected cells. In contrast, cells infected with the ΔvopEF mutant maintained the tight junction and cell polarity relatively well, showing somewhat normal distribution of ZO-1 and actin in large areas throughout the monolayer. Overall, these results are fully consistent with our functional measurement of epithelial barrier function as assessed by TER. Disruption of tight-junction integrity has been identified as the common mechanism for diarrhea in gastrointestinal diseases such as ulcerative colitis (27), Crohn’s disease (28), and gastritis caused by bacterial pathogens (29, 30). To correlate the results from the electrophysiology experiment with those from an in vivo model, we employed the removable intestinal tie-adult rabbit diarrhea (RITARD) assay to study the diarrheal response (Fig. 5D). Out of 10 animals infected with wild-type AM-19226, 9 had fatal diarrhea and 1 had nonfatal diarrhea. ΔvcsN2 mutant-infected animals had no symptoms (10/10). The ΔvopF mutant, the ΔvopE mutant, and the double mutants all had intermediate phenotypes, ranging from fatal and nonfatal diarrhea to no symptoms. These data concur with the in vitro electrophysiology data indicating that vopE and vopF together strongly contribute to the disruption of tight junctions and diarrhea in infected animals.

DISCUSSION

V.  cholerae is a Gram-negative bacterium that colonizes the small intestine and causes a massive diarrheal response. In the pandemic O1 and O139 strains, the toxin-coregulated pilus (TCP) is the major virulence factor that contributes to the colonization of the intestinal epithelium. After colonizing and replicating in the small intestine, the pathogen secretes cholera toxin, which is translocated into the host cells in an AB5 toxin-dependent manner. The O1 and O139 virulence mechanisms seem modular, requiring two separate proteins; however, the results from this study suggest that colonization and diarrheal induction are interconnected in T3SS+ non-O1, non-O139 strains, with a single protein contributing to both processes. Previously, we described a type III effector, VopF, which is an actin nucleator that induces protrusion formation in transfected or infected cells and is required for efficient colonization in the infant mouse model. In this study, we have identified a divergent T3SS pathogenicity island in additional non-O1, non-O139 strains, 1587 and 623-39, which were sequenced as part of the NIAID microbial sequencing project. This T3SS island encodes a VopF homolog, VopN. We have shown that VopN is an actin nucleator like VopF, but unlike VopF, it localizes to stress fiber, similarly to VopL from V. parahaemolyticus. Chimera and transfection studies suggested that the amino-terminal portions of the proteins, preceding the FH1-like domain, are responsible for the differing localizations of the proteins in host cells. A yeast two-hybrid study identified filamin as the interacting partner for VopN. Filamin is a high-molecular-weight actin-cross-linking protein that localizes to stress fibers in cultured cells. In addition to organizing actin filaments into bundled networks, filamins may also be required for formation of filopodia. Additionally, filamin acts as a scaffolding protein for signaling proteins and bridges transmembrane proteins to the actin cytoskeleton. Interestingly, the Salmonella type III effectors Ssph2 and Ssel have been shown to bind filamin and inhibit actin polymerization to avoid vacuole-associated actin polymerization (31). As bacterial pathogens have evolved numerous mechanisms to manipulate the actin cytoskeleton, filamin may be a suitable target for localizing effector proteins and modulating actin dynamics. By studying the interactions of host proteins and VopF, we discovered that AM-19226 affects the cell cycle progression of infected cells and that the presence of vopF prevents apoptosis. While the mechanism remains to be elucidated, actin nucleation and polymerization activity of VopF are not required for the antiapoptotic activity. Numerous bacterial effectors, such as cytolethal distending toxin (CDT) and cycle-inhibiting factor (Cif) (23), have been shown to affect the eukaryotic cell cycle. Arresting the cell cycle has been proposed to benefit pathogens infecting the intestinal epithelium (23), which normally undergoes a high rate of cell turnover, with a cellular life span of less than 5 days in humans. The newly identified VopN protein provided a tool for dissecting the actin nucleation function away from the cellular phenotype induced by the bacterial effectors. By complementing a vopF deletion with vopN, we determined that vopN can only partially complement the defect. Surprisingly, deletion of vopN in 1587 or 1582 had a minimal defect in colonization in the infant mouse model. The lack of a defect in colonization may be explained by the presence of a gene encoding a putative cytotoxic necrotizing factor (CNF). CNF belongs to a family of Rho GTPase-activating proteins, which includes dermonecrotic toxin (DNT) from Bordetella, SopE/SopE2 from Salmonella, IpgB1/2 from Shigella, and Map from enteropathogenic Escherichia coli (EPEC)/enterohemorrhagic E. coli (EHEC) (32, 33). CNF from E. coli has been shown to affect tight junctions (34) and forms filopodia and stress fibers (35). The presence of CNF may be functionally redundant to VopN or VopF even though CNF is found only in the 1587 T3SS and V. parahaemolyticus T3SS2 pathogenicity islands. Interestingly, CNF is also a cyclomodulin that induces multinucleation (23) and prevents apoptosis by remodeling the mitochondrial network (36). Although CNF is not present in the AM-19226 T3SS pathogenicity island, a type III effector named VopE was identified. VopE is a homolog of YopE from Yersinia. Similar to YopE, VopE has a putative GTPase-activating protein domain. By inactivating RhoA, Rac1, and Cdc42, YopE depolymerizes actin and inhibits phagocytosis by macrophages and other immune cells (25, 37). YopE has also been shown to depolarize the epithelial layer (38). Since V. cholerae colonizes the small intestine, which is lined with a single layer of polarized cells, we investigated how AM-19226 would interact with a polarized-T84-cell model. As the TER and confocal micrographs showed, AM-19226 translocates VopE and VopF to reorganize the actin cytoskeleton and depolarize the epithelium. Loosening the tight junctions has been a common strategy by which bacterial pathogens cause diarrhea in the human host. We investigated the diarrheal response to AM-19226 in the adult rabbit model and determined that the T3SS is required for causation of diarrhea; vopE and vopF both contribute to this phenotype. With the discovery of the T3SS in V. cholerae and its role in the virulence of non-O1, non-O139 strains, it is astonishing to observe how this bacterium has evolved two independent pathogenic mechanisms to cause similar diseases. The presence of a homologous T3SS pathogenicity island in other non-O1, non-O139 V. cholerae strains and an orthologous T3SS2 island in the sequenced pandemic V. parahaemolyticus strain confirms that the T3SS has emerged to be a viable conserved pathogenic mechanism among Vibrio species. The genetic events and mechanisms in which virulent strains gain the ability to cause epidemics and pandemics have yet to be elucidated, and further investigation of T3SS-dependent cholera-like gastroenteritis will help us combat the disease caused by a rapidly evolving and diverse pathogen.

MATERIALS AND METHODS

Cell lines and bacterial strains.

Chinese hamster ovary (CHO) cells (ATCC ccl-61) were cultivated in Ham’s F-12K medium supplemented with 10% fetal calf serum (FCS) and l-glutamine (2 mM) adjusted to contain 1.5 g/liter sodium bicarbonate. HeLa, HEK293, and 293T cells were cultivated in Dulbecco’s modified Eagle’s medium (DMEM) supplemented with 10% fetal bovine serum (FBS) and l-glutamine (2 mM). V. cholerae AM-19226, 1587, and 1582, non-O1, non-O139 strains of serogroup O39, O12 strains, and N16961, an O1 El Tor biotype, were used in this study. E. coli and V. cholerae strains were propagated in LB containing appropriate antibiotics at 37°C.

Strain construction.

Deletion mutations were constructed by in-frame deletions of the entire reading frame using a technique described previously (39). vopF and vopN genes were complemented into the genome at the vopF locus using homologous recombination.

Protein expression and purification.

vopF and vopN genes were cloned into pET28a (EMD Biosciences), which introduces an N-terminal fusion of the protein to a His6 tag. BL21(DE3) strains carrying expression plasmids were subcultured 1:1,000 into 1 liter of 2× yeast extract-tryptone (YT) medium growing at 37°C with aeration until the optical density at 600 nm reached 0.5 and were induced using 0.8 mM IPTG (isopropyl-β-d-thiogalactopyranoside) and incubated at 18°C overnight. Cells were harvested and suspended in 500 mM NaCl, 50 mM Tris (pH 7.5), 10% glycerol, 15 mM imidazole, 2 mM β-mercaptoethanol and lysed by lysozyme and sonication. The fusion protein was purified by affinity chromatography using His-Trap HP (GE Health Systems).

Pyrene actin assay.

The pyrene actin assay was described previously (40). Rabbit skeletal muscle actin was purified as described previously (41). Protein samples were dialyzed using Slide-a-lyzer (Pierce) in assay buffer for 2 h, and 100 nM of VopF or VopN protein was used in the assays unless specified.

Eukaryotic cell transfection.

vopF and vopN and various chimera genes were cloned into pAcGFP-N1 (Clontech) to make transfection vectors expressing GFP fusion proteins. pAcGFP-N1, which encodes a GFP from Aequorea coerulescens, was also used as a positive control for transfection. Plasmid DNA was extracted using an endotoxin-free maxiprep kit (Qiagen). CHO cells were seeded onto 18-mm glass coverslips (no. 1.5). jetPEI (polyPlus) transfection reagent was used as recommended by the manufacturer. At 24 h posttransfection, cells were fixed with 3.7% formaldehyde and permeabilized with 0.1% Triton X-100. Actin filaments were stained with Alexa Fluor 568-conjugated phalloidin (Invitrogen). The slides were then covered with coverslips and observed with a Nikon inverted fluorescence microscope.

Immunogold electron microscopy.

CHO cells were seeded and transfected as described above. At 24 h posttransfection, cells were permeabilized with 0.1% Triton X-100 and fixed in 4% paraformaldehyde-0.025% gluteraldehyde. After being blocked with 1% fish gelatin, cells were stained with mouse monoclonal anti-GFP antibodies (1:50 dilution) (JL-8; Clontech) for 1 h and then incubated with a bridging rabbit anti-mouse antibody (1:65 dilution) for 30 min. Lastly, colloidal gold (10 nm)-conjugated secondary antibodies were incubated with the samples for 1 h. A secondary fixation was performed with 3% gluteraldehyde. Samples were postfixed with 1% osmium and stained with 1% aqueous uranyl acetate. After dehydration, a propylene oxide-epon araldite mixture was used to infiltrate the samples, which were mounted and sectioned. Microscopy was performed using a Tecnai G2 Spirit BioTWIN, and images were captured with an AMT 2k charge-coupled-device (CCD) camera.

Yeast two-hybrid screen.

vopF and vopN genes were cloned into pGBT9 (Clontech) and transformed into AH109 strains. The bait strains were screened against a pretransformed human HeLa matchmaker cDNA library (Clontech) on synthetic defined (SD)/−Ade/−His/−Trp/−Leu/X-α-Gal (5-bromo-4-chloro-3-indolyl alpha-d-galactopyranoside) agar plates. Blue colonies were restreaked multiple times to confirm the phenotype. Plasmids were rescued from Saccharomyces cerevisiae and transformed into E. coli KC8 cells. Purified DNA was reintroduced into yeast cells to confirm interactions.

Coimmunoprecipitation.

Genes encoding the VopF/VopN-interacting proteins were cloned into the pProLabel-C vector using InFusion cloning (Clontech). HEK293 cells were cotransfected with vopF1-135::pAcGFP or vopN1-84::pAcGFP constructs and various prey fusion constructs. At 36 h posttransfection, cells were lysed using M-PER mammalian protein extraction reagent (Thermo Scientific) in the presence of protease inhibitors (Roche). GFP fusion proteins were immunoprecipitated using polyclonal anti-GFP antibodies (Clontech) with protein A Dynabeads (Invitrogen). ProLabel fusion protein was detected using a chemiluminescence ProLabel assay (Clontech). The fluorescence and chemiluminescence of the cell lysates and immunoprecipitated fractions were determined using an Envision3 plate reader (PerkinElmer).

Luciferase assay.

Luciferase assays were performed using the dual-luciferase reporter assay system (Promega) according to the manufacturer’s instructions. The levels of the firefly luciferase activities were measured and normalized to the activity of phRL-TK-derived Renilla luciferase. The results are expressed as the means ± standard errors (SE) from triplicate experiments.

Cell cycle analysis and TUNEL staining.

The cell cycle distribution and apoptotic population of HeLa cells were determined as previously described (42, 43). Briefly, ~5 × 108 HeLa cells were infected with V. cholerae strains at a multiplicity of infection (MOI) of 100 for 4 h. Cells were washed with phosphate-buffered saline (PBS) and incubated in medium with 100 µg/ml gentamicin for 72 h. Cells were then washed, trypsinized, fixed in 75% ethanol, and stained in a propidium iodide-RNase A solution (Invitrogen). Flow cytometric analysis was performed with a FACSCalibur flow cytometer (Becton Dickinson) using the PI emission (630 nm) for DNA quantification. Data from >104 cells were collected and analyzed using CELLQUEST software (Becton Dickinson). Apoptotic cells were stained using an APO-BrdU TUNEL assay kit (Invitrogen).

Polarized-T84-epithelial-cell infection.

Electrophysiology experiments were carried out as previously described (24). Briefly, T84 cells were cultured in 0.33-cm[2] Transwell inserts in T84 medium. Experiments were performed 10 to 12 days after the plating. Cells were rinsed in Hanks’ balanced salt solution and were transferred into CO2-independent medium with 1% serum. Transwell inserts were kept at 37°C without CO2 on a plate warmer or in an incubator for the course of the experiments. Isc and TER were measured at 10-min intervals using a dual voltage clamp device. V. cholerae strains were subcultured to log phase in LB at 37°C, washed, and diluted to 109 CFU per ml. In all experiments, 10 µl of bacteria was pipetted into the apical chamber of the Transwell insert. After the electrophysiology experiment, cells were fixed using 3.7% formaldehyde and permeabilized using 0.1% Triton X-100. Rabbit anti-AM-19226 antibodies, mouse anti-ZO-1 antibodies (Zymogen), and Alexa Fluor 568-conjugated goat anti-rabbit antibodies (Invitrogen) were used to visualize the bacteria or VopF protein. Alexa Fluor 633-conjugated phalloidin (Invitrogen) was used to stain the actin. The slides were then covered with coverslips and observed on a Nikon spinning-disk confocal microscope.

Infant mouse colonization.

The infant mouse colonization assay was described previously (44). V. cholerae bacteria (1 × 106) were inoculated intragastrically into 5-day-old CD-1 mice (Charles River Laboratories). Infected pups were killed after 20 h, and intestinal colonized bacteria were quantified as described previously (44). A competition assay was performed by coinoculation with the ΔlacZ wild-type AM-19226, 1587, or 1582 strain. Differentiation of wild-type or mutant strains was based on blue/white colonies on LB plates containing X-Gal (5-bromo-4-chloro-3-indolyl-β-d-galactopyranoside). The competitive index was calculated as the ratio of the number of CFU of the test strain to the number of CFU of the reference strain divided by the ratio of the size of the inoculum of the test strain to the size of the inoculum of the reference strain. If no CFU were recovered from the intestine, a colony value of 1 was assigned to the lowest dilution.

Rabbit diarrheal assay.

Diarrheal response of rabbits to V. cholerae was assayed by the removable intestinal tie-adult rabbit diarrhea (RITARD) model, using adult New Zealand White rabbits weighing 1.5 to 2.7 kg (45). After inoculation, rabbits were observed at 6-h intervals for 7 days for death, diarrheal severity, and shedding of challenge organisms. All rabbits that died were subjected to postmortem examination to assess the amount of fluid in the intestine.

Ethics statement.

The animal experiments were performed with protocols approved by the Harvard Medical School Office for Research Protection Standing Committee on Animals. This committee approved the animal protocol used in this study. The Harvard Medical School animal management program is accredited by the Association for the Assessment and Accreditation of Laboratory Animal Care, International (AAALAC), and meets National Institutes of Health standards as set forth in the Guide for the Care and Use of Laboratory Animals (DHHS [NIH] publication no. 85-23, revised 1996). The institution also accepts as mandatory the PHS Policy on Humane Care and Use of Laboratory Animals by Awardee Institutions and NIH Principles for the Utilization and Care of Vertebrate Animals Used in Testing, Research, and Training. There is on file with the Office of Laboratory Animal Welfare (OLAW) an approved Assurance of Compliance (A3431-01).

Sequence and comparative genomic analysis.

The BLAST-Extend-Repraze (BER) algorithm and BLAST searches were used to compare nucleic acid and amino acid sequence homologies (46, 47). Artemis version 8 from the Sanger Institute was used to visualize and annotate sequence data (48). ACT from the Sanger Institute was used to visualize BLASTn and tBLASTx comparisons (49).

Statistical analysis.

The paired differences of experimental groups were compared using the nonparametric Mann-Whitney U test. A P value of <0.05 was considered statistically significant.
  49 in total

1.  Positive selection vectors for allelic exchange.

Authors:  K Skorupski; R K Taylor
Journal:  Gene       Date:  1996-02-22       Impact factor: 3.688

2.  Deamidation of Cdc42 and Rac by Escherichia coli cytotoxic necrotizing factor 1: activation of c-Jun N-terminal kinase in HeLa cells.

Authors:  M Lerm; J Selzer; A Hoffmeyer; U R Rapp; K Aktories; G Schmidt
Journal:  Infect Immun       Date:  1999-02       Impact factor: 3.441

3.  A high proportion of Vibrio cholerae strains isolated from children with diarrhoea in Bangkok, Thailand are multiple antibiotic resistant and belong to heterogenous non-O1, non-O139 O-serotypes.

Authors:  A Dalsgaard; A Forslund; L Bodhidatta; O Serichantalergs; C Pitarangsi; L Pang; T Shimada; P Echeverria
Journal:  Epidemiol Infect       Date:  1999-04       Impact factor: 2.451

4.  Distribution and virulence of Vibrio cholerae belonging to serogroups other than O1 and O139: a nationwide survey.

Authors:  A K Mukhopadhyay; P K Saha; S Garg; S K Bhattacharya; T Shimada; T Takeda; Y Takeda; G B Nair
Journal:  Epidemiol Infect       Date:  1995-02       Impact factor: 2.451

5.  NF-kappaB function in growth control: regulation of cyclin D1 expression and G0/G1-to-S-phase transition.

Authors:  M Hinz; D Krappmann; A Eichten; A Heder; C Scheidereit; M Strauss
Journal:  Mol Cell Biol       Date:  1999-04       Impact factor: 4.272

6.  Cluster of cases of clinical cholera due to Vibrio cholerae 010 in east Delhi.

Authors:  S Rudra; R Mahajan; M Mathur; K Kathuria; V Talwar
Journal:  Indian J Med Res       Date:  1996-02       Impact factor: 2.375

7.  Alterations in Vibrio cholerae motility phenotypes correlate with changes in virulence factor expression.

Authors:  C L Gardel; J J Mekalanos
Journal:  Infect Immun       Date:  1996-06       Impact factor: 3.441

8.  Molecular analysis of non-O1, non-O139 Vibrio cholerae associated with an unusual upsurge in the incidence of cholera-like disease in Calcutta, India.

Authors:  C Sharma; M Thungapathra; A Ghosh; A K Mukhopadhyay; A Basu; R Mitra; I Basu; S K Bhattacharya; T Shimada; T Ramamurthy; T Takeda; S Yamasaki; Y Takeda; G B Nair
Journal:  J Clin Microbiol       Date:  1998-03       Impact factor: 5.948

9.  The sixth and seventh cholera pandemics are due to independent clones separately derived from environmental, nontoxigenic, non-O1 Vibrio cholerae.

Authors:  D K Karaolis; R Lan; P R Reeves
Journal:  J Bacteriol       Date:  1995-06       Impact factor: 3.490

10.  Characterization of Vibrio cgolerae non-O1 serogroups obtained from an outbreak of diarrhea in Lima, Peru.

Authors:  A Dalsgaard; M J Albert; D N Taylor; T Shimada; R Meza; O Serichantalergs; P Echeverria
Journal:  J Clin Microbiol       Date:  1995-10       Impact factor: 5.948

View more
  27 in total

1.  Multiple forms of Spire-actin complexes and their functional consequences.

Authors:  Christine K Chen; Michael R Sawaya; Martin L Phillips; Emil Reisler; Margot E Quinlan
Journal:  J Biol Chem       Date:  2012-02-08       Impact factor: 5.157

Review 2.  Bacterial nucleators: actin' on actin.

Authors:  Joana N Bugalhão; Luís Jaime Mota; Irina S Franco
Journal:  Pathog Dis       Date:  2015-09-27       Impact factor: 3.166

3.  Dimeric WH2 domains in Vibrio VopF promote actin filament barbed-end uncapping and assisted elongation.

Authors:  Julien Pernier; Jozsef Orban; Balendu Sankara Avvaru; Antoine Jégou; Guillaume Romet-Lemonne; Bérengère Guichard; Marie-France Carlier
Journal:  Nat Struct Mol Biol       Date:  2013-08-04       Impact factor: 15.369

4.  Type 3 Secretion System Island Encoded Proteins Required for Colonization by Non-O1/non-O139 Serogroup Vibrio cholerae.

Authors:  Mudit Chaand; Kelly A Miller; Madeline K Sofia; Cory Schlesener; Jacob W A Weaver; Vibha Sood; Michelle Dziejman
Journal:  Infect Immun       Date:  2015-05-04       Impact factor: 3.441

5.  Vibrio cholerae T3SS effector VopE modulates mitochondrial dynamics and innate immune signaling by targeting Miro GTPases.

Authors:  Masato Suzuki; Olga Danilchanka; John J Mekalanos
Journal:  Cell Host Microbe       Date:  2014-10-23       Impact factor: 21.023

Review 6.  Vibrio variations on a type three theme.

Authors:  Kelly A Miller; Katharine F Tomberlin; Michelle Dziejman
Journal:  Curr Opin Microbiol       Date:  2019-01-31       Impact factor: 7.934

7.  Characterization of V. cholerae T3SS-dependent cytotoxicity in cultured intestinal epithelial cells.

Authors:  Kelly A Miller; Mudit Chaand; Stacy Gregoire; Takeshi Yoshida; Lisa A Beck; Andrei I Ivanov; Michelle Dziejman
Journal:  Cell Microbiol       Date:  2016-07-22       Impact factor: 3.715

8.  Filament assembly by Spire: key residues and concerted actin binding.

Authors:  Amy S Rasson; Justin S Bois; Duy Stephen L Pham; Haneul Yoo; Margot E Quinlan
Journal:  J Mol Biol       Date:  2014-09-16       Impact factor: 5.469

Review 9.  The WH2 Domain and Actin Nucleation: Necessary but Insufficient.

Authors:  Roberto Dominguez
Journal:  Trends Biochem Sci       Date:  2016-04-05       Impact factor: 13.807

10.  Regulation by ToxR-Like Proteins Converges on vttRB Expression To Control Type 3 Secretion System-Dependent Caco2-BBE Cytotoxicity in Vibrio cholerae.

Authors:  Kelly A Miller; Madeline K Sofia; Jacob W A Weaver; Christopher H Seward; Michelle Dziejman
Journal:  J Bacteriol       Date:  2016-05-13       Impact factor: 3.490

View more

北京卡尤迪生物科技股份有限公司 © 2022-2023.