Identification of type III secretion substrates of Chlamydia trachomatis using Yersinia enterocolitica as a heterologous system
- Maria da Cunha1, 3,
- Catarina Milho†1, 4,
- Filipe Almeida†1, 3,
- Sara V Pais1, 3,
- Vítor Borges2,
- Rui Maurício1, 5,
- Maria José Borrego2,
- João Paulo Gomes2 and
- Luís Jaime Mota1, 3Email author
© da Cunha et al.; licensee BioMed Central Ltd. 2014
Received: 6 October 2013
Accepted: 10 February 2014
Published: 17 February 2014
Chlamydia trachomatis is an obligate intracellular human pathogen causing ocular and urogenital infections that are a significant clinical and public health concern. This bacterium uses a type III secretion (T3S) system to manipulate host cells, through the delivery of effector proteins into their cytosol, membranes, and nucleus. In this work, we aimed to find previously unidentified C. trachomatis T3S substrates.
We first analyzed the genome of C. trachomatis L2/434 strain for genes encoding mostly uncharacterized proteins that did not appear to possess a signal of the general secretory pathway and which had not been previously experimentally shown to be T3S substrates. We selected several genes with these characteristics and analyzed T3S of the encoding proteins using Yersinia enterocolitica as a heterologous system. We identified 23 C. trachomatis proteins whose first 20 amino acids were sufficient to drive T3S of the mature form of β-lactamase TEM-1 by Y. enterocolitica. We found that 10 of these 23 proteins were also type III secreted in their full-length versions by Y. enterocolitica, providing additional support that they are T3S substrates. Seven of these 10 likely T3S substrates of C. trachomatis were delivered by Y. enterocolitica into host cells, further suggesting that they could be effectors. Finally, real-time quantitative PCR analysis of expression of genes encoding the 10 likely T3S substrates of C. trachomatis showed that 9 of them were clearly expressed during infection of host cells.
Using Y. enterocolitica as a heterologous system, we identified 10 likely T3S substrates of C. trachomatis (CT053, CT105, CT142, CT143, CT144, CT161, CT338, CT429, CT656, and CT849) and could detect translocation into host cells of CT053, CT105, CT142, CT143, CT161, CT338, and CT429. Therefore, we revealed several C. trachomatis proteins that could be effectors subverting host cell processes.
Chlamydiae are a large group of obligate intracellular bacteria that includes human pathogens (e.g. Chlamydia trachomatis or C. pneumoniae), animal pathogens (e.g. C. abortus, C. caviae, C. felis, or C. muridarum), or symbionts of free-living amoebae. Among Chlamydiae, C. trachomatis is a particular clinical and public health concern, being the leading cause of infectious blindness in developing countries  and the most prevalent sexually transmitted bacteria worldwide .
Like all Chlamydiae, C. trachomatis undergoes a developmental cycle involving the inter-conversion between two morphologically distinct forms: a non-replicative infectious form, the elementary body (EB), and a replicative non-infectious form, the reticulate body (RB) . Throughout its developmental cycle, C. trachomatis uses a type III secretion system (T3SS) to translocate several effector proteins across the host cell plasma membrane and the inclusion membrane [4, 5]. These T3S effectors are thought to play a central role in bacterial invasion [6, 7] and exit of host cells , and in the subversion of various host cell processes [9–16]. There are, however, chlamydial effectors, such as CPAF/CT858 or CT441, which are not T3S substrates .
Given their likely central role during infection, considerable efforts have been placed at identifying chlamydial effectors. This is not a trivial task because the amino acid sequence of most effectors does not display significant similarity to proteins of known function. Additionally, T3S substrates, which should comprise the bulk of Chlamydia effectors, contain no easily recognizable secretion signal. Moreover, in spite of the recent development of systems for transformation of Chlamydia[17, 18], for a long time no methods have been available for genetic manipulation of these bacteria. To overcome these obstacles, chlamydial effectors have been searched: i) by systematic phenotypic analyses of yeast Saccharomyces cerevisiae expressing individual chlamydial proteins ; ii) by using Salmonella, Shigella[15, 21–23], or Yersinia[13, 14, 24–27] as genetically tractable heterologous host bacteria carrying well characterized T3SSs; or iii) by complex computational predictions of T3S signals [28–30]. The subsequent use of specific antibodies enabled to detect translocation into host cells of some of the C. trachomatis proteins singled out in these searches, such as in the case of Tarp/CT456 , CT694 , CopN/CT089 , Cap1/CT529 , CT620 , CT621 [22, 32], CT711 , lipid-droplet associated (Lda) proteins Lda1/CT156, Lda2/CT163, and Lda3/CT473 , Nue/CT737 , or of a group of proteins containing a hydrophobic motif thought to mediate their insertion into the inclusion membrane (Inc proteins) [12, 34]. Moreover, the direct use of antibodies raised against particular C. trachomatis proteins (CT311, CT622, CT795, GlgA/CT798, HtrA/CT823, or Pgp3) revealed their presence in the host cell cytosol or nucleus of infected cells [35–40]. Finally, the in vitro deubiquitinase activity of Chla DUB1/CT868 and of Chla DUB2/CT867 , and the capacity of Chla DUB1/CT868 to suppress the NF-κB pathway in transfected cells , indicate that these two proteins should be effectors.
In this work, we have surveyed the genome of C. trachomatis mostly for genes encoding uncharacterized proteins that were not described before as T3S substrates. We then used Yersinia enterocolitica as a heterologous system to identify 10 novel likely T3S substrates of C. trachomatis and real-time quantitative PCR (RT-qPCR) to show that 9 of the genes encoding these proteins are clearly expressed during the bacterial developmental cycle. Furthermore, we showed that 7 of the 10 likely T3S substrates of C. trachomatis could be translocated into host cells by Y. enterocolitica. Therefore, we identified several novel putative effectors of C. trachomatis.
Cell culture, bacterial strains and growth conditions
HeLa 229 (ATCC) cells were maintained in Dulbecco’s modified Eagle’s medium (DMEM; Invitrogen) supplemented with 10% (v/v) foetal bovine serum (FBS; Invitrogen) at 37°C in a humidified atmosphere of 5% (v/v) CO2. C. trachomatis L2/434/Bu (from ATCC) was propagated in HeLa 229 cells using standard techniques . Escherichia coli TOP10 (Invitrogen) was used for construction and purification of the plasmids. Yersinia enterocolitica ΔHOPEMT (MRS40 pIML421 [yopH Δ1-352, yopO Δ65-558, yopP 23 , yopE 21 , yopM 23 , yopT 135 ]), deficient for the Yersinia T3S effectors YopH, O, P, E, M, and T, but T3S-proficient  and T3S-deficient Y. enterocolitica ΔHOPEMT ΔYscU (MRS40 pFA1001 [yopH Δ1-352, yopO Δ65-558, yopP 23 , yopE 21 , yopM 23 , yopT 135 , yscU Δ1-354 )  were used for T3S assays. The yscU gene encodes an essential component of the Y. enterocolitica T3S system, and the yscU Δ1-354 mutation is non-polar . E. coli or Y. enterocolitica were routinely grown in liquid or solid Luria-Bertani (LB) medium (NZYtech) with the appropriate antibiotics and supplements. Plasmids were introduced into E. coli or Y. enterocolitica by electroporation.
DNA manipulations, plasmids, and primers
The plasmids used in this work and their main characteristics are detailed in Additional file 1: Table S1. The DNA primers used in their construction are shown in Additional file 2: Table S2. Plasmids were constructed and purified with proof-reading Phusion DNA polymerase (Finnzymes), restriction enzymes (MBI Fermentas), T4 DNA Ligase (Invitrogen), DreamTaq DNA polymerase (MBI Fermentas), DNA clean & concentrator™-5 Kit and Zymoclean™ Gel DNA Recovery kit (Zymo Research), and purified with GeneElute Plasmid Miniprep kit (Sigma), according to the instructions of the manufacturers. In brief, to analyze T3S signals we constructed plasmids harboring hybrid genes encoding the first 10, 15, 20, or 40 amino acids of each protein (C. trachomatis proteins, SycT and YopE) and the mature form of TEM-1 β-lactamase (TEM-1) . These hybrids were made using as vector pLJM3, a low-copy plasmid which enables expression of the cloned genes driven by the promoter of the Y. enterocolitica yopE gene , either by overlapping PCR or by using a cloning strategy previously described for the construction of plasmids encoding Inc-TEM-1 hybrid proteins . To analyze secretion of full-length C. trachomatis proteins, we constructed plasmids expressing the proteins C-terminally tagged with a haemagglutinin (HA) epitope. For this, the genes were amplified by PCR from chromosomal DNA of strain L2/434/Bu using a reverse primer with a sequence complementary to the transcribed strand of the DNA encoding the HA-epitope. PCR products digested with the appropriate enzymes were ligated into pLJM3 . The accuracy of the nucleotide sequence of all the inserts in the constructed plasmids was checked by DNA sequencing.
Y. enterocolitica T3S assays
T3S assays were done as previously described . We used Y. enterocolitica ΔHOPEMT or ΔHOPEMT ΔYscU strains carrying the plasmids described in Additional file 1: Table S1. The proteins in bacterial pellets and culture supernatants were analyzed by immunoblotting, and the amounts of protein in bacterial pellets and/or culture supernatants were estimated from images of immunoblots with Image Lab (Bio-Rad). Where appropriate, we calculated the percentage of secretion as the ratio between the amounts of secreted protein (in the culture supernatant fraction) relative to the total amount of protein (in the culture supernatant and in the bacterial pellet fractions). The results from the quantifications are the average ± standard error of the mean (SEM) from at least three independent experiments. Detailed results for each protein analyzed are in Additional file 3: Table S3.
Y. enterocolitica translocation assays
Analyses of protein translocation into host cells by Y. enterocolitica were done essentially as previously described [49, 50]. In brief, Y. enterocolitica strains were grown in brain heart infusion (BHI; Scharlau) medium overnight at 26°C with continuous shaking (130 rpm). Bacteria were then diluted to an optical density at 600 nm of 0.2 in fresh BHI and cultured in the same conditions for 2 h. Subsequently, the yop regulon was induced by incubation for 30 min in a shaking water bath (130 rpm) at 37°C. Bacteria were then washed with DMEM supplemented with 10% (v/v) FBS and added to HeLa 229 cells, grown overnight in 24-well plates (1x105 cells/well), by using a multiplicity of infection of 50. The infected cells were incubated at 37°C in a humidified atmosphere of 5% (v/v) CO2. After 3 h of incubation, extracellular bacteria were killed by adding gentamicin (50 μg/ml), and the cells were incubated in the same conditions for additional 2 h. The infected cells were then harvested on ice, washed with phosphate-buffered saline (PBS), ressuspended in PBS containing 0.1% (v/v) Triton X-100 and a protease inhibitor cocktail (Sigma), and incubated for 10 min on ice. The samples were centrifuged (15,000 g for 15 min at 4°C) and Triton-soluble and Triton-insoluble HeLa cell lysates were loaded on sodium dodecyl sulfate-12% (v/v) polyacrilamide gels. After electrophoresis, the gels were processed for immunoblotting using 0.2 μm pore-size nitrocellulose membranes (BioRad).
The following antibodies were used for immunoblotting: rat monoclonal anti-HA (clone 3F10; Roche; used at 1:1000), mouse monoclonal anti-TEM-1 (QED Bioscience; 1:500), rabbit polyclonal anti-SycO (1:1000) , and mouse monoclonal anti-tubulin (clone B-5-1-2; Sigma; 1:1000). Immunoblot detection was done with horseradish peroxidase-conjugated secondary antibodies (GE Healthcare and Jackson ImmunoResearch), Western Lightning Plus-ECL (Perkin Elmer), and a ChemiDoc XRS + system (BioRad) or exposure to Amersham Hyperfilm ECL (GE Healthcare). All quantitative analyses were done with immunoblot images obtained using ChemiDoc XRS + (BioRad).
Real-time quantitative PCR
The expression of the newly identified candidate T3S substrates during the developmental cycle of C. trachomatis L2/434 was estimated by determining mRNA levels at different times post-infection by real-time quantitative PCR (RT-qPCR). These experiments were done as previously described . Primers (available upon request) were designed using Primer Express (Applied Biosystems). The RT-qPCR assays were done using the ABI 7000 SDS, SYBR green chemistry, and optical plates (Applied Biosystems), as previously described . At each time point, raw RT-qPCR data for each gene were normalized against the data obtained for the 16S rRNA transcript, as it was previously demonstrated that this is an adequate endogenous control . The final results were based on three independent experiments.
Selection of C. trachomatis proteins analyzed in this work
To search for previously unidentified T3S substrates of C. trachomatis, we first surveyed the genome of strain L2/434 for genes encoding mostly uncharacterized proteins, or with a putative biochemical activity compatible with the function of a T3S effector (e.g., proteases). Among these genes, we selected those encoding proteins not predicted to have a signal sequence characteristic of the general secretory pathway (according to Psortb v3.0) and that had not been previously analyzed experimentally for the presence of a T3S signal. This singled out 32 proteins (CT016, CT017, CT031, CT051, CT053, CT080, CT105, CT142, CT143, CT144, CT153, CT161, CT172, CT273, CT277, CT289, CT309, CT330, CT338, CT386, CT425, CT568, CT583, CT590, CT631, CT635, CT656, CT696, CT702, CT837, CT845, and CT849; we used the nomenclature of the annotated C. trachomatis D/UW3 strain ; the names of the corresponding genes as annotated for strain L2/434  can be found in Additional file 3: Table S3). Furthermore, for comparison purposes, we considered proteins that had been tested for the presence of a T3S signal using Shigella flexneri as a heterologous bacteria: eight proteins whose first ~40 amino acids of the corresponding C. pneumoniae homologs did not drive secretion of an adenylate cyclase (Cya) reporter protein by S. flexneri (CT066, CT429, GrgA/CT504, CT538, CT584, CT768, CT779, CT814), and three proteins whose N-terminal region of the C. pneumoniae homologs drove secretion of a Cya reporter protein by S. flexneri (CT203, CT577, CT863) . Please note that at the time this work was initiated GrgA/CT504 was an uncharacterized protein; however, it was recently described as a transcriptional activator . Finally, throughout this study we used as positive controls a C. trachomatis bona-fide T3S effector (CT694)  and a C. trachomatis likely T3S substrate (CT082) that we had previously identified , and which was recently independently confirmed , and as negative control a predicted ribosomal protein (RplJ/CT317).
In summary, in experiments that will be described below, we analyzed T3S signals in 46 C. trachomatis proteins (~5% of all proteins encoded by the L2/434 strain): 32 hypothetical proteins previously not analyzed experimentally for T3S signals, 11 proteins whose C. pneumoniae homologs were previously analyzed for T3S signals using S. flexneri as heterologous system, and 3 controls. In the selection of these proteins, we did not consider predictions made by any of the published in silico methods that suggest putative T3S substrates [28–30, 56].
The first 20 amino acids of C. trachomatis T3S substrates are sufficient to drive efficient secretion of TEM-1 hybrid proteins by Y. enterocolitica
Identification of T3S signals in C. trachomatis proteins
Analysis of the secretion of the newly identified candidate T3S substrates of C. trachomatis as full-length proteins
We next analyzed if the 23 C. trachomatis proteins carrying newly identified T3S signals, and also CT203 and the controls (CT082, CT694 and RplJ), were secreted as full-length proteins by Y. enterocolitica ΔHOPEMT. The rationale for these experiments was that some proteins cannot be type III secreted even with a T3S signal grafted at their N-termini [59–62], possibly because the secretion channel is too narrow (inner diameter of 2–3 nm ) to accommodate tightly folded proteins. For example, while we showed that YopE15-TEM-1 is efficiently type III secreted, hybrid proteins containing the first 15 or 16 amino acids of YopE fused to mouse dihydrofolate reductase (DHFR) are not type III secreted by Y. enterocolitica[59, 60]. This indicates that most T3S substrates must have particular folding properties that are compatible with them being type III secreted proteins. Based on this, we predicted that if the full-length version of chlamydial proteins were type III secreted by Yersinia this would be an additional indication that they can be T3S substrates. However, lack of secretion of the full-length proteins would not preclude that they could be T3S substrates, as they may require Chlamydia-specific chaperones, not present in Yersinia.
Secretion of full-length CT153-HA, CT172-HA, CT203-HA, CT386-HA or CT425-HA by Y. enterocolitica could occasionally be seen by immunoblotting (Figure 3A); however, this was not always reproducible and individual average percentage of secretion of these proteins was in all cases below 2% (Figure 3B). We did not detect significant amounts of CT273-HA, CT289-HA, CT309-HA, or CT631-HA in culture supernatants (Figure 3A and Additional file 3: Table S3), but as their levels of expression were either extremely low (CT273-HA, CT289-HA, and CT309-HA) or undetectable (CT631-HA) it was not possible to draw conclusions about secretion of these proteins. Furthermore, CT016-HA, and possibly CT696-HA (barely visible in Figure 3A), were immunodetected in the culture supernatant fraction in a form that migrated on SDS-PAGE at a molecular weight much lower than the one predicted from their amino acid sequence (27 kDa and 46 kDa, respectively), while in the bacterial pellet fraction their migration on SDS-PAGE corresponded roughly to their predicted molecular weight (Figure 3A). This suggests that the proteins could be cleaved during secretion, unstable in the culture supernatant, or their encoding genes possess internal Shine-Dalgarno sequences. Regardless of the exact reason, we could not confidently analyze whether CT016-HA and CT696-HA were secreted or not.
Overall, the full set of T3S assays revealed 10 proteins (CT053, CT105, CT142, CT143, CT144, CT161, CT338, CT429, CT656, and CT849) as newly identified likely T3S substrates of C. trachomatis, and therefore as possible effectors.
CT053, CT105, CT142, CT143, CT161, CT338, and CT429 can be translocated into host cells by Y. enterocolitica
In summary, these experiments showed that CT053-HA, CT105-HA, CT142-HA, CT143-HA, CT161-HA, CT338-HA and CT429-HA have the capacity of being translocated into infected host cells further suggesting that the endogenous C. trachomatis proteins could be effectors. The results do not preclude that CT144, CT656 or CT849 could be effectors, but the evidence is not as strong as for the other 7 proteins.
Expression of genes encoding newly identified likely T3S substrates during development of C. trachomatis
Earlier studies using heterologous systems have led to the identification of several bona-fide or putative C. trachomatis T3S effectors [13–15, 21, 22, 24–27]. While these and other analyses covered a significant portion of all C. trachomatis proteins, we hypothesized that there could be previously unidentified T3S substrates. By combining basic bioinformatics searches, exhaustive T3S assays, translocation assays, and analyses of chlamydial gene expression in infected cells, we revealed 10 C. trachomatis proteins (CT053, CT105, CT142, CT143, CT144, CT161, CT338, CT429, CT656, and CT849) as likely T3S substrates and possible effectors. In particular, CT053, CT105, CT142, CT143, CT338, and CT429 were type III secreted by Y. enterocolitica, could be translocated into host cells, and their encoding genes were clearly expressed in C. trachomatis strain L2/434. Therefore, these 6 proteins have a high likelihood of being effectors. However, additional future studies are required to show that all of these 10 proteins are indeed translocated by C. trachomatis into host cells and to show that they are bona-fide effectors, i.e., that they interfere with host cell processes.
Among the likely T3S effectors of C. trachomatis that we identified, CT105 and CT142 have been previously singled out as possible modulators of host cell functions, based on the phenotypic consequences of their ectopic expression in yeast S. cerevisiae. In addition, the genes encoding CT142, CT143, and CT144 have been shown to be markedly transcriptionally regulated by a protein (Pgp4) encoded by the Chlamydia virulence plasmid . This plasmid is present in almost all C. trachomatis clinical isolates , and studies in animal models of infection showed that it is a virulence factor in vivo[67, 68]. Additional studies are needed to understand if the putative effector function of CT142, CT143, and CT144 can partially explain the virulence role of the chlamydial plasmid. Furthermore, the predicted amino acid sequence of CT849 reveals a domain of unknown function (DUF720) that can only be found in Chlamydia proteins. In C. trachomatis, besides CT849, a DUF720 domain is found in CT847, a T3S effector that interacts with human G rap2 c yclin D-i nteracting p rotein (GCIP) , and in CT848, which has been indicated as a T3S substrate using S. flexneri as a heterologous system . Therefore, this further supports a possible role of CT849 as an effector. In contrast with CT105, CT142, CT143, CT144 or CT849, no significant information is available or could be retrieved about CT053, CT338, CT429, or CT656.
CT161 is a possible T3S substrate and effector, but we could not detect significant levels of ct161 mRNA during the developmental cycle of strain L2/434. The ct161 gene is localized within the “plasticity zone”, a chromosomal region of rare high genetic diversity among C. trachomatis strains. In fact, although C. trachomatis includes strains showing remarkably different tropisms (strains that can spread into lymph nodes and cause lymphogranuloma venereum [LGV], such as L2/434, and strains causing infections usually restricted to the mucosa of the conjunctiva and genitals), their genomes are all highly similar . Preliminary data indicate that, contrarily to what is seen in LGV strains, the ct161 seems to be more expressed in some ocular and urogenital isolates (data not shown). We are currently investigating the possibility that ct161 is a pseudogene in LGV strains, perhaps inactivated by a mutation in its promoter region. Interestingly, CT161 has been shown by yeast two-hybrid to bind CT274 (a possible chlamydial T3S chaperone) . Another feature of this protein is that part of its amino acid sequence (residues 40–224, out of 246) shows 28% of identity to a region of Lda2/CT163 (residues 167–361, out of 548), a known C. trachomatis translocated protein .
Among the proteins for which we found a secretion signal but could not demonstrate their T3S as full-length proteins, we highlight CT153 and GrgA/CT504. Regarding CT153, this protein possesses a membrane attack complex/perforin (MACPF) domain , and there is previous evidence that it may be translocated by C. trachomatis, which is consistent with our data. The ct504 gene has been recently shown to encode a transcriptional activator, GrgA . Therefore, T3S of CT50420-TEM-1 could be false a positive. However, if GrgA is a T3S substrate, as our data suggests, it could have a function within the host cell or, more likely and similarly to what has been described in the T3SSs of Yersinia or Shigella[74, 75], it could be discarded by secretion once its intra-bacterial regulatory activity needs to be shut down.
We found T3S signals in 56% proteins analyzed (26 out of 46, including controls). This high percentage of proteins showing a T3S signal suggests that some should be false positives. It is conceivable that within a single bacterium non-secreted proteins possess T3S signals but are not targeted to the T3SS machinery because they also carry signals (e.g. DNA-, membrane-, or protein-binding) that preferentially direct them to other location within the bacterial cell. To help differentiating between true or false positives among chlamydial proteins carrying a T3S signal we analyzed their secretion as full-length proteins. This is because, as explained above in the Results section, not all proteins have folding characteristics compatible with T3S [59–62]. However, we cannot exclude that some of the C. trachomatis full-length proteins that were not type III secreted by Yersinia have a T3S chaperone that maintains them in a secretion-competent state  and enables their secretion during infection by C. trachomatis. Intriguingly, CT082 or CT694 have dedicated T3S chaperones, CT584 and Slc1, respectively , and, in agreement with what we previously observed , they were both secreted as full-length proteins in the absence of the chaperones. Considering that T3S chaperones have various functions [76, 77], the chaperone role of CT584 or Slc1 should be different from maintaining their substrates in a secretion-competent state.
Eleven of the Chlamydia proteins that we analyzed have been previously studied for T3S using S. flexneri has a heterologous system . In the majority of the cases the outcome of the experiments was identical; however, differently from what was shown in Shigella, we detected a T3S signal in the N-terminal of CT429 (which was also secreted as a full-length protein, and could be translocated into HeLa cells), GrgA/CT504, and CT779 and we did not detect a T3S signal in CT577. Evidence for a T3S signal in only one of the heterologous systems may suggest a false positive. However, there is a myriad of possible explanations for these discrepancies, when considering that different heterologous systems (Shigella and Yersinia) and reporter proteins (Cya and TEM-1) were used, and that the N-terminal regions in the hybrid proteins consisted in different lengths of amino acids and were in some cases from different Chlamydia species.
We compared the data from our T3S assays (including the controls, CT082, CT694, and RplJ) with predictions of T3S substrates by in silico methods (Effective T3S , SIEVE , Modlab , and T3_MM ) using resources available in the Web (Effective T3S, Modlab and T3_MM) and Table three in reference  (SIEVE), as detailed in Additional file 3: Table S3. When considering the analysis of T3S signals in TEM-1 hybrids, the vast majority of proteins (60%; 12 out of 20) in which we did not find a T3S signal were also predicted not to be secreted by each of the in silico methods. In contrast, the vast majority of proteins (58%; 15 out of 26) in which we detected a T3S signal were also predicted to be secreted by at least one of the in silico methods. The correlation between our experimental data and the in silico predictions was more striking when considering the T3S of full-length proteins. Among the 16 full-length proteins for which we could not find definitive evidence of T3S, 10 (i.e., 62.5%) were also predicted not to be secreted by each of the in silico methods, but among the 11 proteins that we showed or confirmed to be T3S substrates, 10 (i.e., 83%) were also predicted to be secreted by at least one of the in silico methods. Overall, this indicates some correlation between our experimental data and the in silico methods that predict T3S substrates. However, for many proteins, each of these in silico methods generates different predictions (see Additional file 3: Table S3). It is possible that the quantitative data on T3S such as the one we generated in this and in a previous study , can be used to normalize and improve the predictive value of such methods.
We found 10 C. trachomatis proteins (CT053, CT105, CT142, CT143, CT144, CT161, CT338, CT429, CT656, and CT849) with a high likelihood of being T3S substrates, and therefore possible effectors delivered by the bacteria into host cells. For 6 of these proteins (CT053, CT105, CT142, CT143, CT338, and CT429), the hypothesis that they could be effectors was supported by their capacity of being translocated into host cells and by the expression of their encoding genes by C. trachomatis. The identification of all C. trachomatis effectors is a crucial step towards a comprehensive understanding of the mechanisms by which this pathogen subverts host cells. The recently developed methods for genetic manipulation of Chlamydia indicate that it should be possible to ectopically express candidate effectors in C. trachomatis[17, 78], which would facilitate the analysis of their translocation into host cells. Our work highlights C. trachomatis proteins that should be prioritized in such studies, thus aiding the future identification of chlamydial effectors. Furthermore, the quantitative analysis of T3S of TEM-1 hybrid proteins that we carried out could help to further develop the in silico methods for identification of T3S substrates [28–30, 56].
This work was supported by Fundação para a Ciência e a Tecnologia (FCT) through grants ERA-PTG/0005/2010 (in the frame of ERA-NET PathoGenoMics) to LJM, ERA-PTG/0004/2010 (in the frame of ERA-NET PathoGenoMics) to JPG, and PEst-OE/EQB/LA0004/2011; by the European Commission through a Marie Curie European Re-integration Grant (PERG03-GA-2008-230954) to LJM; and by a European Society for Clinical Microbiology and Infectious Diseases (ESCMID) research grant to LJM. MdC, FA, and VB hold PhD fellowships SFRH/BD/62728/2009, SFRH/BD/73545/2010, and SFRH/BD/68527/2010, respectively, from FCT.
- Wright HR, Turner A, Taylor HR: Trachoma. Lancet. 2008, 371 (9628): 1945-1954.View ArticlePubMedGoogle Scholar
- Bebear C, de Barbeyrac B: Genital Chlamydia trachomatis infections. Clin Microbiol Infect. 2009, 15 (1): 4-10.View ArticlePubMedGoogle Scholar
- Abdelrahman YM, Belland RJ: The chlamydial developmental cycle. FEMS Microbiol Rev. 2005, 29 (5): 949-959.View ArticlePubMedGoogle Scholar
- Betts HJ, Wolf K, Fields KA: Effector protein modulation of host cells: examples in the Chlamydia spp. arsenal. Curr Opin Microbiol. 2009, 12 (1): 81-87.View ArticlePubMedGoogle Scholar
- Valdivia RH: Chlamydia effector proteins and new insights into chlamydial cellular microbiology. Curr Opin Microbiol. 2008, 11 (1): 53-59.View ArticlePubMedGoogle Scholar
- Jewett TJ, Miller NJ, Dooley CA, Hackstadt T: The conserved Tarp actin binding domain is important for chlamydial invasion. PLoS Pathog. 2010, 6 (7): e1000997-PubMed CentralView ArticlePubMedGoogle Scholar
- Lane BJ, Mutchler C, Al Khodor S, Grieshaber SS, Carabeo RA: Chlamydial entry involves TARP binding of guanine nucleotide exchange factors. PLoS Pathog. 2008, 4 (3): e1000014-PubMed CentralView ArticlePubMedGoogle Scholar
- Lutter EI, Barger AC, Nair V, Hackstadt T: Chlamydia trachomatis inclusion membrane protein CT228 recruits elements of the myosin phosphatase pathway to regulate release mechanisms. Cell Rep. 2013, 3 (6): 1921-1931.PubMed CentralView ArticlePubMedGoogle Scholar
- Scidmore MA, Hackstadt T: Mammalian 14-3-3beta associates with the Chlamydia trachomatis inclusion membrane via its interaction with IncG. Mol Microbiol. 2001, 39 (6): 1638-1650.View ArticlePubMedGoogle Scholar
- Delevoye C, Nilges M, Dehoux P, Paumet F, Perrinet S, Dautry-Varsat A, Subtil A: SNARE protein mimicry by an intracellular bacterium. PLoS Pathog. 2008, 4 (3): e1000022-PubMed CentralView ArticlePubMedGoogle Scholar
- Rzomp KA, Moorhead AR, Scidmore MA: The GTPase Rab4 interacts with Chlamydia trachomatis inclusion membrane protein CT229. Infect Immun. 2006, 74 (9): 5362-5373.PubMed CentralView ArticlePubMedGoogle Scholar
- Mital J, Miller NJ, Fischer ER, Hackstadt T: Specific chlamydial inclusion membrane proteins associate with active Src family kinases in microdomains that interact with the host microtubule network. Cell Microbiol. 2010, 12 (9): 1235-1249.PubMed CentralView ArticlePubMedGoogle Scholar
- Chellas-Gery B, Linton CN, Fields KA: Human GCIP interacts with CT847, a novel Chlamydia trachomatis type III secretion substrate, and is degraded in a tissue-culture infection model. Cell Microbiol. 2007, 9 (10): 2417-2430.View ArticlePubMedGoogle Scholar
- Hower S, Wolf K, Fields KA: Evidence that CT694 is a novel Chlamydia trachomatis T3S substrate capable of functioning during invasion or early cycle development. Mol Microbiol. 2009, 72 (6): 1423-1437.View ArticlePubMedGoogle Scholar
- Pennini ME, Perrinet S, Dautry-Varsat A, Subtil A: Histone methylation by NUE, a novel nuclear effector of the intracellular pathogen Chlamydia trachomatis. PLoS Pathog. 2010, 6 (7): e1000995-PubMed CentralView ArticlePubMedGoogle Scholar
- Derre I, Swiss R, Agaisse H: The lipid transfer protein CERT interacts with the Chlamydia inclusion protein IncD and participates to ER-Chlamydia inclusion membrane contact sites. PLoS Pathog. 2011, 7 (6): e1002092-PubMed CentralView ArticlePubMedGoogle Scholar
- Wang Y, Kahane S, Cutcliffe LT, Skilton RJ, Lambden PR, Clarke IN: Development of a transformation system for Chlamydia trachomatis: restoration of glycogen biosynthesis by acquisition of a plasmid shuttle vector. PLoS Pathog. 2011, 7 (9): e1002258-PubMed CentralView ArticlePubMedGoogle Scholar
- Gerard HC, Mishra MK, Mao G, Wang S, Hali M, Whittum-Hudson JA, Kannan RM, Hudson AP: Dendrimer-enabled DNA delivery and transformation of Chlamydia pneumoniae. Nanomedicine. 2013, 9 (7): 996-1008.View ArticlePubMedGoogle Scholar
- Sisko JL, Spaeth K, Kumar Y, Valdivia RH: Multifunctional analysis of Chlamydia-specific genes in a yeast expression system. Mol Microbiol. 2006, 60 (1): 51-66.View ArticlePubMedGoogle Scholar
- Ho TD, Starnbach MN: The Salmonella enterica serovar Typhimurium-encoded type III secretion systems can translocate Chlamydia trachomatis proteins into the cytosol of host cells. Infect Immun. 2005, 73 (2): 905-911.PubMed CentralView ArticlePubMedGoogle Scholar
- Subtil A, Delevoye C, Balana ME, Tastevin L, Perrinet S, Dautry-Varsat A: A directed screen for chlamydial proteins secreted by a type III mechanism identifies a translocated protein and numerous other new candidates. Mol Microbiol. 2005, 56 (6): 1636-1647.View ArticlePubMedGoogle Scholar
- Muschiol S, Boncompain G, Vromman F, Dehoux P, Normark S, Henriques-Normark B, Subtil A: Identification of a family of effectors secreted by the type III secretion system that are conserved in pathogenic Chlamydiae. Infect Immun. 2011, 79 (2): 571-580.PubMed CentralView ArticlePubMedGoogle Scholar
- Furtado AR, Essid M, Perrinet S, Balana ME, Yoder N, Dehoux P, Subtil A: The chlamydial OTU domain-containing protein ChlaOTU is an early type III secretion effector targeting ubiquitin and NDP52. Cell Microbiol. 2013, 15 (12): 2064-2079.View ArticlePubMedGoogle Scholar
- Fields KA, Hackstadt T: Evidence for the secretion of Chlamydia trachomatis CopN by a type III secretion mechanism. Mol Microbiol. 2000, 38 (5): 1048-1060.View ArticlePubMedGoogle Scholar
- Clifton DR, Fields KA, Grieshaber SS, Dooley CA, Fischer ER, Mead DJ, Carabeo RA, Hackstadt T: A chlamydial type III translocated protein is tyrosine-phosphorylated at the site of entry and associated with recruitment of actin. Proc Natl Acad Sci U S A. 2004, 101 (27): 10166-10171.PubMed CentralView ArticlePubMedGoogle Scholar
- Pais SV, Milho C, Almeida F, Mota LJ: Identification of novel type III secretion chaperone-substrate complexes of Chlamydia trachomatis. PLoS ONE. 2013, 8 (2): e56292-PubMed CentralView ArticlePubMedGoogle Scholar
- Hovis KM, Mojica S, McDermott JE, Pedersen L, Simhi C, Rank RG, Myers GS, Ravel J, Hsia RC, Bavoil PM: Genus-optimized strategy for the identification of chlamydial type III secretion substrates. Pathog Dis. 2013, 69 (3): 213-222.View ArticlePubMedGoogle Scholar
- Arnold R, Brandmaier S, Kleine F, Tischler P, Heinz E, Behrens S, Niinikoski A, Mewes HW, Horn M, Rattei T: Sequence-based prediction of type III secreted proteins. PLoS Pathog. 2009, 5 (4): e1000376-PubMed CentralView ArticlePubMedGoogle Scholar
- Samudrala R, Heffron F, McDermott JE: Accurate prediction of secreted substrates and identification of a conserved putative secretion signal for type III secretion systems. PLoS Pathog. 2009, 5 (4): e1000375-PubMed CentralView ArticlePubMedGoogle Scholar
- Lower M, Schneider G: Prediction of type III secretion signals in genomes of gram-negative bacteria. PLoS One. 2009, 4 (6): e5917-PubMed CentralView ArticlePubMedGoogle Scholar
- Fling SP, Sutherland RA, Steele LN, Hess B, D’Orazio SE, Maisonneuve J, Lampe MF, Probst P, Starnbach MN: CD8+ T cells recognize an inclusion membrane-associated protein from the vacuolar pathogen Chlamydia trachomatis. Proc Natl Acad Sci U S A. 2001, 98 (3): 1160-1165.PubMed CentralView ArticlePubMedGoogle Scholar
- Hobolt-Pedersen AS, Christiansen G, Timmerman E, Gevaert K, Birkelund S: Identification of Chlamydia trachomatis CT621, a protein delivered through the type III secretion system to the host cell cytoplasm and nucleus. FEMS Immunol Med Microbiol. 2009, 57 (1): 46-58.PubMed CentralView ArticlePubMedGoogle Scholar
- Kumar Y, Cocchiaro J, Valdivia RH: The obligate intracellular pathogen Chlamydia trachomatis targets host lipid droplets. Curr Biol. 2006, 16 (16): 1646-1651.View ArticlePubMedGoogle Scholar
- Li Z, Chen C, Chen D, Wu Y, Zhong Y, Zhong G: Characterization of fifty putative inclusion membrane proteins encoded in the Chlamydia trachomatis genome. Infect Immun. 2008, 76 (6): 2746-2757.PubMed CentralView ArticlePubMedGoogle Scholar
- Lei L, Qi M, Budrys N, Schenken R, Zhong G: Localization of Chlamydia trachomatis hypothetical protein CT311 in host cell cytoplasm. Microb Pathog. 2011, 51 (3): 101-109.PubMed CentralView ArticlePubMedGoogle Scholar
- Gong S, Lei L, Chang X, Belland R, Zhong G: Chlamydia trachomatis secretion of hypothetical protein CT622 into host cell cytoplasm via a secretion pathway that can be inhibited by the type III secretion system inhibitor compound 1. Microbiology. 2011, 157 (Pt 4): 1134-1144.PubMed CentralView ArticlePubMedGoogle Scholar
- Qi M, Lei L, Gong S, Liu Q, DeLisa MP, Zhong G: Chlamydia trachomatis secretion of an immunodominant hypothetical protein (CT795) into host cell cytoplasm. J Bacteriol. 2011, 193 (10): 2498-2509.PubMed CentralView ArticlePubMedGoogle Scholar
- Lu C, Lei L, Peng B, Tang L, Ding H, Gong S, Li Z, Wu Y, Zhong G: Chlamydia trachomatis GlgA Is Secreted into Host Cell Cytoplasm. PLoS ONE. 2013, 8 (7): e68764-PubMed CentralView ArticlePubMedGoogle Scholar
- Li Z, Chen D, Zhong Y, Wang S, Zhong G: The chlamydial plasmid-encoded protein pgp3 is secreted into the cytosol of Chlamydia-infected cells. Infect Immun. 2008, 76 (8): 3415-3428.PubMed CentralView ArticlePubMedGoogle Scholar
- Lei L, Dong X, Li Z, Zhong G: Identification of a novel nuclear localization signal sequence in Chlamydia trachomatis-secreted hypothetical protein CT311. PLoS ONE. 2013, 8 (5): e64529-PubMed CentralView ArticlePubMedGoogle Scholar
- Misaghi S, Balsara ZR, Catic A, Spooner E, Ploegh HL, Starnbach MN: Chlamydia trachomatis-derived deubiquitinating enzymes in mammalian cells during infection. Mol Microbiol. 2006, 61 (1): 142-150.View ArticlePubMedGoogle Scholar
- Le Negrate G, Krieg A, Faustin B, Loeffler M, Godzik A, Krajewski S, Reed JC: ChlaDub1 of Chlamydia trachomatis suppresses NF-kappaB activation and inhibits IkappaBalpha ubiquitination and degradation. Cell Microbiol. 2008, 10 (9): 1879-1892.View ArticlePubMedGoogle Scholar
- Scidmore MA: Cultivation and laboratory maintenance of Chlamydia trachomatis. Curr Protoc Microbiol. 2005, 00:11A.1.1–11A.1.25Google Scholar
- Iriarte M, Cornelis GR: YopT, a new Yersinia Yop effector protein, affects the cytoskeleton of host cells. Mol Microbiol. 1998, 29 (3): 915-929.View ArticlePubMedGoogle Scholar
- Almeida F, Borges V, Ferreira R, Borrego MJ, Gomes JP, Mota LJ: Polymorphisms in Inc proteins and differential expression of inc genes among Chlamydia trachomatis strains correlate with invasiveness and tropism of lymphogranuloma venereum isolates. J Bacteriol. 2012, 194 (23): 6574-6585.PubMed CentralView ArticlePubMedGoogle Scholar
- Sorg I, Wagner S, Amstutz M, Muller SA, Broz P, Lussi Y, Engel A, Cornelis GR: YscU recognizes translocators as export substrates of the Yersinia injectisome. EMBO J. 2007, 26 (12): 3015-3024.PubMed CentralView ArticlePubMedGoogle Scholar
- Charpentier X, Oswald E: Identification of the secretion and translocation domain of the enteropathogenic and enterohemorrhagic Escherichia coli effector Cif, using TEM-1 beta-lactamase as a new fluorescence-based reporter. J Bacteriol. 2004, 186 (16): 5486-5495.PubMed CentralView ArticlePubMedGoogle Scholar
- Marenne MN, Journet L, Mota LJ, Cornelis GR: Genetic analysis of the formation of the Ysc-Yop translocation pore in macrophages by Yersinia enterocolitica: role of LcrV, YscF and YopN. Microb Pathog. 2003, 35 (6): 243-258.View ArticlePubMedGoogle Scholar
- Denecker G, Totemeyer S, Mota LJ, Troisfontaines P, Lambermont I, Youta C, Stainier I, Ackermann M, Cornelis GR: Effect of low- and high-virulence Yersinia enterocolitica strains on the inflammatory response of human umbilical vein endothelial cells. Infect Immun. 2002, 70 (7): 3510-3520.PubMed CentralView ArticlePubMedGoogle Scholar
- Grosdent N, Maridonneau-Parini I, Sory MP, Cornelis GR: Role of Yops and adhesins in resistance of Yersinia enterocolitica to phagocytosis. Infect Immun. 2002, 70: 4165-4176.PubMed CentralView ArticlePubMedGoogle Scholar
- Letzelter M, Sorg I, Mota LJ, Meyer S, Stalder J, Feldman M, Kuhn M, Callebaut I, Cornelis GR: The discovery of SycO highlights a new function for type III secretion effector chaperones. EMBO J. 2006, 25 (13): 3223-3233.PubMed CentralView ArticlePubMedGoogle Scholar
- Borges V, Ferreira R, Nunes A, Nogueira P, Borrego MJ, Gomes JP: Normalization strategies for real-time expression data in Chlamydia trachomatis. J Microbiol Methods. 2010, 82 (3): 256-264.View ArticlePubMedGoogle Scholar
- Stephens RS, Kalman S, Lammel C, Fan J, Marathe R, Aravind L, Mitchell W, Olinger L, Tatusov RL, Zhao Q, et al: Genome sequence of an obligate intracellular pathogen of humans: Chlamydia trachomatis. Science. 1998, 282 (5389): 754-759.View ArticlePubMedGoogle Scholar
- Thomson NR, Holden MT, Carder C, Lennard N, Lockey SJ, Marsh P, Skipp P, O’Connor CD, Goodhead I, Norbertzcak H, et al: Chlamydia trachomatis: genome sequence analysis of lymphogranuloma venereum isolates. Genome Res. 2008, 18 (1): 161-171.PubMed CentralView ArticlePubMedGoogle Scholar
- Bao X, Nickels BE, Fan H: Chlamydia trachomatis protein GrgA activates transcription by contacting the nonconserved region of sigma66. Proc Natl Acad Sci U S A. 2012, 109 (42): 16870-16875.PubMed CentralView ArticlePubMedGoogle Scholar
- Wang Y, Sun M, Bao H, White AP: T3_MM: a Markov model effectively classifies bacterial type III secretion signals. PLoS ONE. 2013, 8 (3): e58173-PubMed CentralView ArticlePubMedGoogle Scholar
- Sory MP, Boland A, Lambermont I, Cornelis GR: Identification of the YopE and YopH domains required for secretion and internalization into the cytosol of macrophages, using the cyaA gene fusion approach. Proc Natl Acad Sci U S A. 1995, 92 (26): 11998-12002.PubMed CentralView ArticlePubMedGoogle Scholar
- Lloyd SA, Norman M, Rosqvist R, Wolf-Watz H: Yersinia YopE is targeted for type III secretion by N-terminal, not mRNA, signals. Mol Microbiol. 2001, 39 (2): 520-531.View ArticlePubMedGoogle Scholar
- Feldman MF, Muller S, Wuest E, Cornelis GR: SycE allows secretion of YopE-DHFR hybrids by the Yersinia enterocolitica type III Ysc system. Mol Microbiol. 2002, 46 (4): 1183-1197.View ArticlePubMedGoogle Scholar
- Lee VT, Schneewind O: Yop fusions to tightly folded protein domains and their effects on Yersinia enterocolitica type III secretion. J Bacteriol. 2002, 184 (13): 3740-3745.PubMed CentralView ArticlePubMedGoogle Scholar
- Akeda Y, Galan JE: Chaperone release and unfolding of substrates in type III secretion. Nature. 2005, 437 (7060): 911-915.View ArticlePubMedGoogle Scholar
- Sorg JA, Miller NC, Marketon MM, Schneewind O: Rejection of impassable substrates by Yersinia type III secretion machines. J Bacteriol. 2005, 187 (20): 7090-7102.PubMed CentralView ArticlePubMedGoogle Scholar
- Cornelis GR: The type III secretion injectisome. Nat Rev Microbiol. 2006, 4 (11): 811-825.View ArticlePubMedGoogle Scholar
- Stebbins CE, Galan JE: Maintenance of an unfolded polypeptide by a cognate chaperone in bacterial type III secretion. Nature. 2001, 414 (6859): 77-81.View ArticlePubMedGoogle Scholar
- Song L, Carlson JH, Whitmire WM, Kari L, Virtaneva K, Sturdevant DE, Watkins H, Zhou B, Sturdevant GL, Porcella SF, et al: Chlamydia trachomatis plasmid-encoded Pgp4 is a transcriptional regulator of virulence-associated genes. Infect Immun. 2013, 81 (3): 636-644.PubMed CentralView ArticlePubMedGoogle Scholar
- Rockey DD: Unraveling the basic biology and clinical significance of the chlamydial plasmid. J Exp Med. 2011, 208 (11): 2159-2162.PubMed CentralView ArticlePubMedGoogle Scholar
- Kari L, Whitmire WM, Olivares-Zavaleta N, Goheen MM, Taylor LD, Carlson JH, Sturdevant GL, Lu C, Bakios LE, Randall LB, et al: A live-attenuated chlamydial vaccine protects against trachoma in nonhuman primates. J Exp Med. 2011, 208 (11): 2217-2223.PubMed CentralView ArticlePubMedGoogle Scholar
- Olivares-Zavaleta N, Whitmire W, Gardner D, Caldwell HD: Immunization with the attenuated plasmidless Chlamydia trachomatis L2(25667R) strain provides partial protection in a murine model of female genitourinary tract infection. Vaccine. 2010, 28 (6): 1454-1462.PubMed CentralView ArticlePubMedGoogle Scholar
- Harris SR, Clarke IN, Seth-Smith HM, Solomon AW, Cutcliffe LT, Marsh P, Skilton RJ, Holland MJ, Mabey D, Peeling RW, et al: Whole-genome analysis of diverse Chlamydia trachomatis strains identifies phylogenetic relationships masked by current clinical typing. Nat Genet. 2012, 44 (4): 413-419. S411PubMed CentralView ArticlePubMedGoogle Scholar
- Spaeth KE, Chen YS, Valdivia RH: The Chlamydia type III secretion system C-ring engages a chaperone-effector protein complex. PLoS Pathog. 2009, 5 (9): e1000579-PubMed CentralView ArticlePubMedGoogle Scholar
- Ponting CP: Chlamydial homologues of the MACPF (MAC/perforin) domain. Curr Biol. 1999, 9 (24): R911-R913.View ArticlePubMedGoogle Scholar
- Taylor LD, Nelson DE, Dorward DW, Whitmire WM, Caldwell HD: Biological characterization of Chlamydia trachomatis plasticity zone MACPF domain family protein CT153. Infect Immun. 2010, 78 (6): 2691-2699.PubMed CentralView ArticlePubMedGoogle Scholar
- Pettersson J, Nordfelth R, Dubinina E, Bergman T, Gustafsson M, Magnusson KE, Wolf-Watz H: Modulation of virulence factor expression by pathogen target cell contact. Science. 1996, 273 (5279): 1231-1233.View ArticlePubMedGoogle Scholar
- Parsot C, Ageron E, Penno C, Mavris M, Jamoussi K, d’Hauteville H, Sansonetti P, Demers B: A secreted anti-activator, OspD1, and its chaperone, Spa15, are involved in the control of transcription by the type III secretion apparatus activity in Shigella flexneri. Mol Microbiol. 2005, 56 (6): 1627-1635.View ArticlePubMedGoogle Scholar
- Botteaux A, Sory MP, Biskri L, Parsot C, Allaoui A: MxiC is secreted by and controls the substrate specificity of the Shigella flexneri type III secretion apparatus. Mol Microbiol. 2009, 71 (2): 449-460.View ArticlePubMedGoogle Scholar
- Feldman MF, Cornelis GR: The multitalented type III chaperones: all you can do with 15 kDa. FEMS Microbiol Lett. 2003, 219 (2): 151-158.View ArticlePubMedGoogle Scholar
- Parsot C, Hamiaux C, Page AL: The various and varying roles of specific chaperones in type III secretion systems. Curr Opin Microbiol. 2003, 6 (1): 7-14.View ArticlePubMedGoogle Scholar
- Agaisse H, Derre I: A C. trachomatis cloning vector and the generation of C. trachomatis strains expressing fluorescent proteins under the control of a C. trachomatis promoter. PLoS ONE. 2013, 8 (2): e57090-PubMed CentralView ArticlePubMedGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly credited. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.