Introduction

Ubiquitination is a conserved and ubiquitous post-translational modification (PTM) of proteins involving the conjugation of the carboxy-terminal glycine residue of ubiquitin (Ub) to lysine residues of target proteins. Poly-ubiquitination of substrates involves further conjugation of a Ub internal lysine residue or amino-terminal methionine (M1) with a second Ub molecule. Seven internal lysines in Ub (K6, K11, K27, K29, K33, K48, K63) and M1 are utilized by Ub conjugating enzymes to form homogeneous, branched, or mixed poly-ubiquitin (polyUb) chains (Komander and Rape., 2012). PolyUb chains with different linkage types exhibit distinct structures and functions. For example, K48- and K11-linked polyUb chains exhibit a compact conformation and are substrates for 26S proteasome-mediated degradation (Varadan et al., 2004; Tenno et al., 2004; Eddins et al., 2007; Bremm et al., 2010; Saeki., 2017). In contrast, K63-linked polyUb conjugates adopt more open conformations that enable the recruitment of multiprotein complexes that regulate the function of the target protein by proteolytic independent events (Komander et al., 2009b; Weeks et al., 2009; Datta et al., 2009; Komander and Rape., 2012). Mixed and branched polyUb chains are also emerging as important regulators of physiological functions (Swatek and Komander., 2016; Ohtake and Tsuchiya., 2017).

Protein ubiquitination regulates numerous eukaryotic cell processes including protein degradation, signal transduction, cell cycle regulation, selective autophagy, the DNA damage response, and programmed cell death. Ub also plays key roles in modulating host innate immune responses to bacterial infection (Li et al., 2016), bacterial proteins, pathogen-containing vacuoles, and bacteria themselves by targeting them for Ub-mediated degradation by proteasomal or autophagic machineries (Li et al., 2016). Because the Ub system is critical for pathogen containment, many pathogens have evolved mechanisms to counteract the impact of this PTM (Vozandychova et al., 2021). For instance, bacterial deubiquitinases (DUBs) can remove Ub from ubiquitinated substrates thereby dampening inflammatory and cell-autonomous defense mechanisms (Kubori et al., 2019). Many DUBs are cysteine proteases with a catalytic Cys, a nearby His and an Asn/Asp (Komander et al., 2009a). DUBs are typically dedicated to the removal of Ub moieties and are unable to hydrolyze other Ub-like (Ubl) modifications such as SUMO or NEDD8. However, the CE clan of Ubl proteases (ULPs) can catalyze the removal of both SUMO and NEDD8 (Ronau et al., 2016). Bacterial pathogens also encode CE clan enzymes that function as DUBs, ULPs or both. For instance, Salmonella Typhimurium SseL, Escherichia coli ElaD, and Shigella flexneri ShiCE function as Ub specific proteases (Rytkönen et al., 2007; Catic et al., 2007; Pruneda et al., 2016) while RickCE from Rickettsia belli functions as a protease directed towards both Ub and NEDD8 as does SidE from Legionella pneumophila which displays mixed activities towards Ub, NEDD8, and ISG15 (Sheedlo et al., 2015; Pruneda et al., 2016). Similarly, XopD from Xanthamonas campestris and LotB from L. pneumophila are isopeptidases exhibiting cross reactivity towards both Ub and SUMO (Pruneda et al., 2016; Schubert et al., 2020). Some CE clan bacterial effectors display acetyltransferase activity. L. pneumophila LegCE, S. Typhimurium AvrA, and YopJ from Yersinia pestis function exclusively as acetyltransferases (Mittal et al., 2006; Mukherjee et al., 2006; Jones et al., 2008; Pruneda et al., 2016). In contrast, the Chlamydia trachomatis (Ct) effector Cdu1/ChlaDUB1 is a CE clan protein that exhibits both acetyltransferase and deubiquitinating activities (Misaghi et al., 2006; Pruneda et al., 2016; Fischer et al., 2017; Pruneda et al., 2018).

Ct is an obligate intracellular bacterial pathogen responsible for human diseases of significant clinical and public health importance (Haggerty et al., 2010). Ct has a biphasic developmental cycle in which the Ct infectious propagule or elementary body (EB) invades the target host cell. Upon internalization the EB transitions to the reticulate body (RB). RBs replicate by binary fission within a pathogenic vacuole (“inclusion”) and asynchronously differentiate back to EBs. In cell culture, starting at around 48 hours-post infection (hpi) Ct will egress after lysis of the host cell or by a process termed extrusion, wherein the entire inclusion is exocytosed from the infected cell (Moulder., 1991; Abdelrahman and Belland., 2005; Hybiske and Stephens., 2007; Lee et al., 2018). The effector Cdu1 was originally identified as a deneddylating and deubiquitinating enzyme and subsequently shown to exhibit in-vitro isopeptidase activity towards both Lys48 and Lys63 linked di-Ub substrates (Misaghi et al., 2006; Claessen et al., 2013; Pruneda et al., 2016; Fischer et al., 2017). Cdu1 is unique among CE clan enzymes in that it also functions as a bona fide lysine acetylase with both is acetylase and DUB activities catalyzed by the same catalytic active site (Pruneda et al., 2018). Intriguingly, Cdu1 autoacetylation is directed towards lysines unlike other CE clan acetylases that predominantly target serine and threonine residues (Pruneda et al., 2018). In infected cells, Cdu1 localizes to the inclusion membrane where it functions to stabilize the anti-apoptotic protein Mcl-1 and to promote the repositioning of Golgi ministacks around the Ct inclusion (Fischer et al., 2017; Wang et al., 2018; Pruneda et al., 2018; Kunz et al., 2019; Auer et al., 2020). However, the mechanism by which Cdu1 promotes redeployment of Golgi ministacks and any additional roles that Cdu1 may play during Ct infection of epithelial cells remains unknown.

In this study, we show that Cdu1 protects itself and three secreted Ct effectors, InaC, IpaM, and CTL0480 from targeted ubiquitination and proteasomal degradation. InaC, IpaM, and CTL0480 are members of a larger family of bacterial proteins embedded within the inclusion membrane (Inc proteins) (Bannantine et al., 2000; Rockey et al., 2002; Chen et al., 2006; Li et al., 2008; Alzhanov et al., 2009; Dehoux et al., 2011; Lutter et al., 2012; Lutter et al., 2013; Kokes et al., 2015; Weber et al., 2015). We show that Cdu1-mediated protection from degradation is independent from its DUB activity but relies upon its Lys acetylase activity. We show that Cdu1 protects InaC to promote repositioning of Golgi ministacks and formation of actin scaffolds around the Ct inclusion, and CTL0480 to promote recruitment of myosin phosphatase target subunit 1 (MYPT1) to the inclusion. In addition, we determined that Cdu1 and Cdu1-protected Incs are required for optimal extrusion of inclusions from host cells at the late stages of infection.

Results

The C. trachomatis inclusion membrane proteins InaC, IpaM, and CTL0480 are differentially ubiquitinated in the absence of Cdu1

Cdu1 is required for Golgi repositioning around the Ct inclusion (Pruneda et al., 2018; Auer et al., 2020). To understand how Cdu1 promotes Golgi redistribution, we first generated a cdu1 null strain in a Ct L2 (LGV L2 434 Bu) background by TargeTron mediated insertional mutagenesis (pDFTT3-aadA) (Lowden et al., 2015) (S. Figure 1). Loss of Cdu1 expression in the resulting L2 cdu1::GII aadA (cdu1::GII) strain was verified by western blot analysis and by indirect immunofluorescence with antibodies raised against Cdu1 (S. Figures 2A and 2B). Because cdu2 resides directly downstream of the cdu1 locus and encodes a second Ct DUB (Cdu2/ChlaDUB2) (Misaghi et al., 2006), we first determined whether the disruption of cdu1 impacted the expression of cdu2. We detected cdu1 and cdu2 transcripts in HeLa cells infected with Ct L2 but not for the juncture between cdu1 and cdu2 (S. Figure 2C). In cells infected with Ct cdu1::GII we only detected cdu2 transcripts (S. Figure 2C) confirming that cdu1 and cdu2 are not co-expressed as part of an operon in accordance with previous observations (Albrecht et al., 2010).

We hypothesized that Cdu1s’ DUB activity promoted Golgi redistribution around inclusions and that we could identify potential targets by comparing the protein ubiquitination profile of cells infected with WT or cdu1::GII strains by quantitative mass spectrometry (MS). HeLa cell were mock infected or infected with either WT L2 or cdu1::GII strains. At 24 hpi, poly-ubiquitinated proteins were enriched from lysed cells using Tandem Ubiquitin Binding Entities (TUBEs) (LifeSensors). TUBEs consist of concatenated Ub binding associated domains (UBAs) that bind to polyUb-modified proteins with nanomolar affinities. Poly-ubiquitinated proteins of both human and Ct origin were enriched and identified by quantitative LC-MS/MS analysis.

Over 2,000 non-ubiquitinated proteins co-precipitated with TUBE 1 bound proteins across all three conditions (mock, L2, and cdu1::GII infected HeLa cells) and 3 biological replicates (S. Table 1). Among these, 47 human proteins were significantly enriched in mock infected HeLa cells and 50 human proteins were significantly enriched during Ct infection (L2 and cdu1::GII) (S. Table 3, S. Figure 3) Pathway enrichment analysis revealed that proteins involved in RNA metabolism were overrepresented among co-precipitating proteins from mock infected cells (S. Table 4, S. Figure 4) while no biological pathways or processes were overrepresented in proteins enriched from infected cells (S. Figure 4). We also identified 8 TUBE1 co-precipitating Ct proteins in HeLa cells infected with L2 and cdu1::GII (S. Table 5, S. Figure 3).

TUBE 1 affinity capture lead to the identification of 43 ubiquitinated proteins (35 human proteins and 8 Ct proteins across all 3 conditions and replicates) based on the presence of peptides containing a di-glycine remnant motif (Peng et al., 2003) (S. Tables 6-8). The lack of widespread poly-ubiquitination of either human or Ct proteins in response to Ct infection (Figure 1) was surprising given that wholesale changes in protein ubiquitination has been reported during infection of HeLa cells by intracellular pathogens like S. Typhimurium (Fiskin et al., 2016). Only two human ubiquitinated proteins (ZC3H7A and DDIT4) were found to be significantly enriched in response to WT L2 infection (Figures 1A and 1C, S. Table 8) while only one human protein (MGC3121) was preferentially ubiquitinated in HeLa cells infected with the cdu1::GII mutant strain (Figures 1B and 1C, S. Table 8). In contrast three Ct proteins, InaC (K104, K107, and K149), IpaM (K29), and CTL0480 (K115) were ubiquitinated at Lys residues in the absence of Cdu1 (Figures 1B, 1C, and 1D, S. Table 8). InaC, IpaM, and CTL0480 are Ct effector proteins that localize to the inclusion membrane (Chen et al., 2006; Alzhanov et al., 2009; Lutter et al., 2013; Kokes et al., 2015). These Type 3 secretion substrates belong to a family of over 36 inclusion membrane proteins that contain a signature bi-lobal hydrophobic transmembrane domain (Bannantine et al., 2000; Rockey et al., 2002; Li et al., 2008; Dehoux et al., 2011; Lutter et al., 2012). Because Cdu1 also localizes to the inclusion membrane (Fischer et al., 2017; Wang et al., 2018; Pruneda et al., 2018; Kunz et al., 2019) we postulated that Cdu1 directly protects InaC, IpaM, and CTL0480 from ubiquitination.

The C. trachomatis inclusion membrane proteins InaC, IpaM, and CTL0480 are ubiquitinated in the absence of Cdu1.

(A-C) Volcano plots (pairwise comparisons) of the relative abundance of human and Ct ubiquitinated peptides between mock infected HeLa cells versus HeLa cells infected with L2 (24 hpi) (A), mock infected HeLa cells versus HeLa cells infected with a cdu1 null strain (cdu1::GII) (24 hpi) (B), and HeLa cells infected with L2 (24 hpi) versus HeLa cells infected with a cdu1 null strain (24 hpi) (C). All infections were performed for 24 hours. Significance values were interpolated from 3 independent biological replicates. Ubiquitinated proteins were enriched with magnetic TUBE 1 beads (binds to polyubiquitinated proteins) and peptides identified by quantitative LC MS/MS analysis. Three Ct proteins, InaC, IpaM, and CTL0480 were differentially ubiquitinated in the absence of Cdu1. (D) InaC was ubiquitinated at K104, K107, and K149, IpaM at K290, and CTL0480 at K115 in the absence of Cdu1. TMD: Transmembrane domain. Numbering corresponds to amino acids in the protein sequence of each respective inclusion membrane protein.

Cdu1 associates with InaC, IpaM, and CTL0480

We first determined if Cdu1 co-localized with InaC, IpaM, and CTL0480 at the inclusion membrane. HeLa cells were infected for 24 hours with WT L2 or L2 expressing CTL0480-Flag from its endogenous promoter and immunostained with antibodies against Cdu1, InaC, IpaM, or the Flag epitope. Both InaC and CTL0480-Flag localized throughout the inclusion membrane while IpaM was restricted to discrete microdomains as previously reported (Alzhanov et al., 2009; Dumoux et al., 2015) (Figure 2A). Cdu1 co-localized with InaC and CTL0480-Flag and with IpaM at microdomains (Figures 2A and 2B). All four antibodies specifically recognized their corresponding antigens since immunostaining for Cdu1, InaC, IpaM or the Flag epitope was not observed in Ct strains lacking Cdu1 (cdu1::GII), InaC (M407, Kokes et al., 2015) IpaM (ipaM::GII, Meier et al., 2022), or a strain that does not express CTL0480-Flag (Figure 2A).

Cdu1 associates with InaC, IpaM, and CTL0480.

(A) Co-localization of Cdu1(magenta) with endogenous InaC (green), IpaM (green), and ectopically expressed CTL0480-Flag (green) at the L2 inclusion membrane during infection of HeLa cells at 24 hpi. HeLa cells infected with an inaC null strain (M407), an ipaM null strain (ipaM::GII), and L2 pBOMB were used as antibody specificity controls. DNA stained with Hoechst is shown in blue. Scale bar: 10μm. Images are representative of multiple images captured across three independent replicates. (B) Line scan profiles of fluorescent signal intensities displayed in (A) showing co-localization of fluorescence intensities for endogenous Cdu1 with endogenous InaC, and IpaM and with CTL0480-Flag along the L2 inclusion membrane. (C) Schematic of Cdu1-GFP(C) (L2) fusion (Cdu1-GFP) and Cdu1-GFP variants used in co-transfections of HEK 293 cells. GFP: Green fluorescent protein. TMD: Transmembrane domain. PRD: Proline rich domain. CD: Catalytic domain. FL: Full length. TMD-: Cdu1-GFP variant lacking TMD domain. CD-: Cdu1-GFP variant lacking CD domain. (D-G) Western blot analysis of GFP immunoprecipitates from HEK 293 cells co-transfected with mammalian plasmids expressing: Cdu1-GFP variants and (D) truncated 3XFlag(N)-InaC (D/UW-3/CX CT813, amino acids 96-264), (E) V5(N)-IpaM (L2, full length), (F) V5(N)-CTL0480 (L2, full length), and (G) V5(N)-CpoS (L2, full length). Vec1: Empty pOPINN-GFP vector. Vec2: Empty pDEST53 vector. Vec3: Empty pcDNATM3.1/nV5-DEST vector. Western blot images are representative from two independent experiments.

We next determined if Cdu1 can interact with InaC, IpaM, and CTL0480 by co-transfecting HEK 293 cells with vectors expressing either full length Cdu1-GFP or truncated versions of Cdu1-GFP lacking transmembrane or catalytic domains (Figure 2C), and vectors expressing Flag-InaC, V5-IpaM, and V5-CTL0480. Transfected cells were lysed and Cdu1-GFP was immunoprecipitated with antibodies against GFP. Western blot analysis of the immunoprecipitates showed that Flag-InaC, V5-IpaM, and V5-CTL0480 co-precipitated with Cdu1-GFP (Figures 2D-F). Moreover, the transmembrane domain of Cdu1 was necessary for Cdu1-GFP to interact with all three Incs (Figures 2D-F). The interaction between Cdu1-GFP and the three tagged Incs was specific since we did not detect interactions between Cdu1-GFP and a V5-tagged version of the inclusion membrane protein CpoS (Sixt et al., 2017) (Figure 2G). We conclude from these results that Cdu1 can interact with all three Incs even in the absence of infection and that interactions are facilitated by the transmembrane domain of Cdu1.

Cdu1 protects InaC, IpaM, and CTL0480 proteins from degradation during infection

We next assessed if Cdu1 was required to stabilize endogenous InaC, IpaM, and CTL0480 in infected cells. HeLa cells were infected with either WT L2 or cdu1::GII strains and at various time points in the Ct infectious cycle crude cell lysates were analyzed by western blot to assess the relative abundance of Inc proteins. At 24 hpi, we observed a decrease in InaC protein levels in cells infected with cdu1::GII relative to cells infected with WT L2 (Figure 3A). In contrast, IpaM protein levels did not decrease until 36 hpi in cdu1::GII infected cells (Figure 3B). For CTL0480, we were not able to detect CTL0480 by western blot but were successful in following CTL0480 expression by indirect immunofluorescence (Figure 3C). The relative abundance of CTL0480 at inclusion membranes, like IpaM, was not affected in cdu1::GII inclusions at 24 hpi. However, at 36 hpi, a subpopulation of cells lost CTL0480 immunoreactivity and by 48 hpi, inclusion membranes of the cdu1::GII strain were devoid of CTL0480 while CTL0480 was prominently detected at the inclusion membranes of WT L2 (Figures 3C and 3D). As controls for the specificity of antibodies used for western blots and for indirect immunofluorescence, we included cells infected with Ct lacking ipaM (ipaM::GII, Meier et al., 2022), with an inaC nonsense mutant (M407, Kokes et al., 2015), and with Ct lacking ctl0480 (ctl0480::GII, Shaw et al., 2018). Overall, our results indicate that steady state protein levels of InaC and IpaM and CTL0480 localization at the inclusion membrane, especially at late stages of infection, are dependent on Cdu1, and that Cdu1 acts at different stages in the infection cycle.

Cdu1 stabilizes InaC, IpaM, and CTL0480.

(A) Western blot analysis of InaC in HeLa cells infected for 24 hours with L2, inaC null, and cdu1 null strains. Ct Slc1 and human alpha tubulin were used to determine Ct burdens and equal loading of protein extracts respectively. Western blot images are representative of 3 independent experiments. (B) Western blot analysis of endogenous IpaM in HeLa cells infected with L2, cdu1 null, and ipaM null strains for 24, 36, and 48 hours. Western blot images are representative of 2 independent experiments. (C) Localization of CTL0480 during Ct infection of HeLa cells at 24, 36, and 48 hpi. CTL0480 signal (green) co-localizes with the inclusion membrane protein IncA (magenta) at the L2 inclusion membrane. Arrowheads highlight cdu1 null inclusions lacking CTL0480 at 36 hpi. DNA stained with Hoechst is shown in blue. Scale bar: 10μm. (D) Quantification of CTL0480 localization at Ct inclusion membranes as shown in (C). CTL0480 fluorescent signal was used to determine the percent of inclusions displaying CTL0480 at the inclusion membranes. Representative images and quantification of CTL0480 positive inclusions are derived from Ct inclusions imaged in 10 fields across 3 independent biological replicates. p values were determined by a student paired t-test.

The acetylase activity of Cdu1 is required for Cdu1 to protect itself, InaC, and IpaM from polyubiquitination and proteasomal degradation

The crystal structures of Cdu1 bound to Ub or Coenzyme A indicated that the adenosine and phosphate groups of Coenzyme A make contact with a helix in variable region 3 (VR-3) of Cdu1 while the Ile36-patch of Ub binds to the opposite face of the same helix (Pruneda et al., 2018). Although Cdu1 catalyzes both of its DUB and acetylase (Act) activities with the same active site (Pruneda et al., 2018) the two activities of Cdu1 can be uncoupled by the amino acid substitution K268E in VR-3 which disrupts Coenzyme A binding required for Act activity and by the amino acid substitution I225A in the Ub-binding region of VR-3 required for DUB activity (Pruneda et al., 2018). These substitutions allowed us to test which of Cdu1’s enzymatic functions are required for the observed effects on protein stability. We generated Ct shuttle plasmids (pBOMB4-MCI backbone), expressing WT Cdu1, a catalytically inactive variant of Cdu1 lacking both DUB and Act activities (Cdu1C345A) (Pruneda et al., 2018), a Cdu1 DUB-deficient variant (Cdu1I225A), and a Cdu1 Act-deficient variant (Cdu1K268E). All Cdu1 constructs were expressed as 3X Flag epitope-tagged versions and from the Cdu1 endogenous promoter.

Plasmids expressing each Cdu1 variant were transformed into the cdu1::GII mutant and the resulting strains used to infect HeLa cells for 36 and 48 hours. The levels of endogenous InaC and IpaM in cell extracts of infected cells were monitored by western blot analysis. At 36 hpi, InaC protein levels drastically decreased in cells infected with cdu1 null strains transformed with empty vector or expressing the catalytic inactive variant of Cdu1 (Cdu1C345A-Flag) (Figure 4A). Likewise, IpaM protein levels diminished at 48 hpi during infection with the same strains (Figure 4B). Both InaC and IpaM protein levels were restored to wild type levels in cdu1 null strains complemented with wild type Cdu1-Flag (Figures 4A and 4B). Unexpectedly, cells infected with a cdu1 null strain ectopically expressing the DUB deficient Cdu1 variant (Cdu1I225A-Flag) displayed wild type levels of InaC and IpaM while the Act deficient variant (Cdu1K268E-Flag) did not (Figures 4A and 4B). These results suggest that the acetylase activity of Cdu1 rather than its DUB activity is required for Cdu1’s ability to stabilize InaC and IpaM proteins.

The acetylase activity of Cdu1 is required to stabilize Cdu1, InaC, and IpaM.

(A) Western blot analysis of endogenous InaC and Cdu1-Flag catalytic variants expressed from a plasmid. HeLa cells were infected for 36 hours with L2, a cdu1 null strain, and cdu1 null strains expressing wild type Cdu1-Flag and the Cdu1 variants C345A-Flag (catalytic inactive), I225A-Flag (DUB deficient), and K268E-Flag (Act deficient). Protein lysates from HeLa cells infected with M407 (inaC null) were used to control for the specificity of anti-InaC antibodies. Western blot images are representative of two independent experiments. (B) Western blot analysis of endogenous IpaM and Cdu1-FLAG variants in crude extracts of HeLa cells infected for 48 hours with the same strains as describe in (A). Infection of HeLa cells with ipaM::GII was used to test for the specificity of the anti-IpaM antibody. Western blot images are representative of two independent experiments. (C) Western blot analysis of Cdu1C345A-Flag (catalytic inactive) and Cdu1K268E-Flag (Act deficient) expressed in a cdu1 null strain or wild type L2 background after infection of HeLa cells for 24 hours. Both Cdu1 variants are stabilized by endogenous Cdu1. (D) Western blot analysis of Cdu1-Flag and Cdu1C345A-Flag expressed in a cdu1 null strain and Cdu1C345A-Flag expressed in wild type L2 following immunoprecipitation (anti-Flag) from HeLa cell extracts after infection for 24 hours and treatment with MG132 (25 μM, 5 hours). Western blot image is a representative blot from at least three independent experiments. (E) Western blot analysis of endogenous InaC and IpaM following immunoprecipitation of proteins acetylated at lysines from inclusion membrane subcellular fractions derived from HeLa cells infected with WT L2 for 24 hours. Western blot image is representative of two independent experiments. (F) Western blot analysis (Flag WB) of acetylated lysine immunoprecipitates generated from inclusion membrane subcellular fractions derived from HeLa cells infected with WT L2 strains expressing Cdu1-Flag, CTL0480-Flag, and IpaM-Flag at 40 hpi. Western blot image is representative of two independent experiments. (G) Flag WB of acetylated lysine immunoprecipitates generated from inclusion membrane subcellular fractions derived from HeLa cells infected with Rif-R L2 (M407 parental strain) and with M407 (inaC null) expressing InaC-Flag at 24 hpi. Western blot image is representative of two independent experiments. (H) The initiator methionine, Lys10, Lys126, and Lys263 of Cdu1 are acetylated during Ct infection of HeLa cells (24 hpi). One tyrosine (Y) residue and multiple serine (S) and threonine (T) residues in Cdu1 are also phosphorylated during Ct infection of HeLa cells (24 hpi). Three PX(S/T)P MAPK phosphorylation consensus sequence motifs were identified in the proline rich domain of Cdu1. Modified residues were identified by quantitative LC MS/MS analysis of immunoprecipitated Cdu1-Flag across 3 independent biological replicates. TMD: Transmembrane domain. PR: Proline rich. PSTP: PX(S/T)P motifs. Ac: Acetylation. P: Phosphorylation. Numbering corresponds to amino acids in Cdu1 protein sequence.

When we monitored the stability of each Flag-tagged Cdu1 variant, we observed that the catalytically inactive variant of Cdu1 (C345A-Flag) was destabilized (Figures 4A and 4B-Flag WB). We reasoned that Cdu1 also protects itself from being targeted for degradation in infected cells. As with InaC and IpaM, the acetylase but not the DUB activity of Cdu1 was required for Cdu1’s stability (Figures 4A and 4B). Although the C345A (DUB-, Act-) and K268E (Act-) amino acid substitutions in Cdu1 do not destabilize Cdu1 expressed in E. coli (Pruneda et al., 2018), it was possible that these substitutions impacted the expression and/or folding of Cdu1 in Chlamydia. To determine if these Cdu1 mutants were inherently unstable, we expressed Cdu1C345A-Flag and Cdu1K268E-Flag in WT L2 in the presence of endogenous Cdu1. We found that each Flag-tagged variant was stabilized (Figure 4C), indicating that endogenous Cdu1 protected the catalytically-deficient Cdu1 variants in trans. In addition, western blot analysis of immunoprecipitated Cdu1-Flag and Cdu1C345A-Flag expressed in a cdu1::GII mutant showed that while wild-type Cdu1-Flag was not modified by Lys48-linked poly-ubiquitination, Cdu1C345A-Flag was robustly modified by Lys48-linked polyUb in the presence of the proteasome inhibitor MG132 (Figure 4D). Moreover, endogenous Cdu1 protected Cdu1C345A-Flag from Lys48-linked polyUb when Cdu1C345A-Flag was expressed in a wild type L2 background (Figure 4D). These results indicate that the loss of Cdu1 activity likely leads to its Lys48-linked poly-ubiquitination and subsequent proteasome-dependent degradation.

Because Cdu1 autoacetylates itself and its Act activity is directed towards lysines (Pruneda et al., 2018) we postulated that Cdu1 may stabilize proteins from degradation by acetylating lysine residues that are potential targets of ubiquitination. We tested this hypothesis by assessing whether Cdu1, InaC, CTL0480, and IpaM are acetylated at lysines during infection. Fractions enriched for inclusion membranes were isolated by sub-cellular fractionation from HeLa cells infected with wild-type L2 (24 hpi), and proteins acetylated at lysines were immunoprecipitated. Western blot analysis of acetyl-lysine immunoprecipitates indicated that InaC and IpaM, but not the Inc protein IncA, were acetylated (Figure 4E). Western blot analysis of anti-acetyl-lysine immunoprecipitates of inclusion membrane-enriched membrane fractions of HeLa cells infected with L2 expressing Cdu1-Flag (24 hpi), CTL0480-Flag (40 hpi), and IpaM-Flag (40 hpi) also showed that all three Flag tagged effectors were acetylated at lysines (Figure 4F). We also determined that Flag-tagged InaC expressed in an M407 (inaC null) background was acetylated at lysines as determined by western blot analysis of anti-acetyl-lysine immunoprecipitates (Figure 4G). In addition, we identified acetylated forms of Cdu1 from mass spectrometric analysis of Flag immunoprecipitates derived from extracts of HeLa cells infected with L2 expressing Cdu1-Flag, (24 hpi) (Figure 4H).

We reasoned that the lysine acetylase activity of Cdu1 is a prominent mechanism by which Cdu1 protects client proteins (at 24 hpi), since loss of Cdu1 or expression of the acetylase deficient variant of Cdu1 (Cdu1K268E-Flag) leads to a marked increase in Ub immunostaining at or near the periphery of cdu1::GII inclusions (> 80% of inclusions) compared to HeLa cells infected with cdu1::GII strains complemented with wild type or DUB deficient (I225A-Flag) Cdu1 strains (Figures 5A and 5B). Given that Cdu1 appears to localize exclusively at inclusion membranes we predicted that its activity would be spatially restricted to the inclusion periphery. We tested this premise by co-infecting HeLa cells with an incA null strain (M923 (IncA R197*), Kokes et al., 2015) and the cdu1::GII strain. IncA mediates homotypic fusion of inclusion membranes and loss of IncA results in the accumulation of multiple unfused inclusions in cells infected at high MOIs (Hackstadt et al., 1999; Suchland et al., 2000; Pannekoek et al., 2005) (Figure 5C). As expected, incA mutants which retain Cdu1 activity did not accumulate Ub at or near the periphery of inclusion membranes. In HeLa cells bearing both cdu1::GII and M923 inclusions, cdu1::GII Cdu1C345A and M923 inclusions, or cdu1::GII Cdu1K268E and M923 inclusions, only the M923 inclusions were protected from ubiquitination (Figures 5C and 5D). Based on these observations we conclude that the acetylase activity of Cdu1 protects proteins in cis and that this activity is constrained to the membrane of the pathogenic vacuole consistent with previous reports (Auer et al., 2020).

The DUB activity of Cdu1 is not required for blocking the ubiquitination of inclusion membranes.

(A) Loss of Cdu1’s acetylase activity results in the accumulation of Ub (magenta signal) at the periphery of the inclusion. HeLa cells were infected for 24 hours. Antisera against the membrane protein Cap1 (green) was used to mark inclusion membranes. DNA stained with Hoechst is shown in blue. Scale bar: 10μm. (B) Quantification of ubiquitinated inclusions as shown in (A). The Ub fluorescent signal was used to determine the number of infected cells with Ub decorated inclusions. The total number of ubiquitinated inclusions was divided by the total number of inclusions analyzed (defined by Cap1 staining). 87%, 86%, and 91% of inclusions were decorated with Ub in HeLa cells infected with a cdu1 null strain and cdu1 null strains expressing Cdu1C345A (DUB-, Act-), and Cdu1K268E (Act-) variants respectively. Representative images in (A) and quantification of ubiquitinated inclusions in (B) were obtained from inclusions imaged in 10 fields across 3 independent biological replicates for each strain. p values were determined by a student paired t-test. (C) HeLa cells co-infected with cdu1::GII (cdu1⊗) and incA null (incA⊗, M923) strains at 24 hpi. IncA-deficient inclusions do not fuse with other inclusions. In co-infected cells, Cdu1 present on the inclusion membranes of incA⊗ strains did not block ubiquitination events at or near the inclusion membranes of neighboring cdu1⊗ strains (mCherry signal from pBOMB4-MCI plasmid). DNA stained with Hoechst is shown in blue. Scale bar: 10μm. (D) Quantification of cdu1⊗ inclusions as shown in (C) in which ubiquitination events are observed in regions of cdu1⊗ inclusion membranes that are in direct apposition to incA⊗ inclusion membranes (junctions). The total number of cdu1⊗ inclusions ubiquitinated at inclusion junctions was divided by the total number of inclusions analyzed (Cap1 staining). 66% and 61% of cdu1⊗ inclusions were decorated with Ub at junctions in HeLa cells co-infected with incA⊗ and cdu1⊗ strains or with incA⊗ and a cdu1⊗ strain ectopically expressing Cdu1K268E (Act-) respectively. Representative images in (C) and quantification of cdu1⊗ inclusions ubiquitinated at junctions in (D) are derived from inclusions imaged in 6 fields across 3 independent biological replicates for each condition. p values were determined by a student paired t-test.

Cdu1 is required for F-actin assembly and Golgi ministack repositioning around the Ct inclusion, and for MYPT1 recruitment to Ct inclusions

InaC is required for Ct to assemble F-actin scaffolds and to reposition Golgi mini stacks around the peripheries of inclusion membranes (Kokes et al., 2015; Wesolowski et al., 2017; Haines et al., 2021). Because Cdu1 regulates InaC levels, we predicted that cdu1 mutants would phenocopy inaC mutants. We quantified the number of inclusions surrounded by F-actin cages at 40 hpi. In cells infected with Rif-R L2 (wild type parental strain of M407, Nguyen and Valdivia., 2012; Kokes et al., 2015), 25% of inclusions were surrounded by F-actin, consistent with previous observations (Chin et al., 2012; Kokes et al., 2015) (Figures 6A and 6E). The number of inclusions surrounded by F-actin decreased to 7% in cells infected with M407 (inaC null) and increased to 49% in HeLa cells infected with M407 complemented with wild type InaC (Figures 6A and 6E). A significant increase in inclusions surrounded by F-actin was also observed in host cells infected with cdu1::GII mutants expressing Cdu1-Flag and Cdu1I225A-Flag (DUB-) (52% and 46% respectively) (Figures 6A and 6E) while cdu1::GII mutants transformed with an empty plasmid or expressing Cdu1C345A-Flag (DUB-Act-) and Cdu1K268E-Flag (Act-) resulted in only 8%,13%, and 10% of inclusions enveloped by F-actin respectively (Figures 6A and 6E). From these observations we conclude that the acetylase activity of Cdu1 is required for Ct to promote assembly of F-actin around the Ct inclusion likely through the stabilization of InaC.

The DUB activity of Cdu1 is not required for assembly of F-actin or Golgi ministack repositioning around the Ct inclusion while Cdu1 is required for MYPT1 recruitment to the inclusion.

Representative images of F-actin (green, Alexa FluorTM Phalloidin) assembled around the Ct inclusion (magenta, anti Cdu1 and Cap1 staining) in HeLa cells infected for 40 hours. (B) Representative images of Golgi (anti GM130 staining, magenta) around Ct inclusions (arrowheads) in HeLa cells infected for 24 hours. HeLa cells were infected with an inaC null strain transformed with empty p2TK2 vector (M407 p2TK2), the wild type strain in which M407 was derived from (Rif-R L2) and M407 ectopically expressing full length InaC (M407 p2TK2-InaC). HeLa cells were also infected with L2 pBOMB, cdu1::GII pBOMB, and cdu1 null strains ectopically expressing Cdu1-Flag (Cdu1), and the Cdu1 variants C345A-Flag (DUB-, Act-), I225A-Flag (DUB-), and K268E-Flag (Act-). (C) Representative images of MYPT1 (magenta) at Ct inclusions (green, anti-IncA staining). Arrowheads represent cdu1 null inclusions with low MYPT1 signal. DNA stained with Hoechst is shown in blue in panels A-C. Scale bar: 10μm in panels A-C. (D) Quantification of Golgi dispersal around the Ct inclusion as shown in (B). The length of Golgi dispersed around each Ct inclusion imaged was measured and normalized to the perimeter length of each inclusion (% Golgi dispersion around the inclusion). Representative images in (B) and Golgi dispersal quantified in (D) were performed from inclusions imaged in 6 fields across 3 independent biological replicates. p values were determined by a student paired t-test. (E) Quantification of Ct inclusion surrounded by F-actin as shown in (A) normalized to total inclusions analyzed. Representative images in (A) and quantification of F-actin surrounding Ct inclusions in (E) were obtained from inclusions imaged in 6 fields across 3 independent biological replicates. p values were determined by one-way ANOVAs with a Student-Newman-Keuls post hoc test. (F) Quantification of MYPT1 at Ct inclusions as shown in (C). Representative images in (C) and quantification of MYPT1 recruitment in (F) were obtained from inclusions imaged in 6 fields across 3 independent replicates. p values were determined by a student paired t-test.

We next quantified Golgi dispersal in infected HeLa cells at 24 hpi, a process that is also dependent on InaC (Kokes et al., 2015; Wesolowski et al., 2017). In Hela cells infected with an inaC null strain (M407) Golgi dispersal was limited to only 26% of the Ct inclusion perimeter. In contrast, cells infected with either L2 Rif-R or with M407 complemented with wild type InaC, the Golgi is dispersed around 45% of the inclusion perimeter (Figures 6B and 6D). Similarly, Golgi dispersal around inclusions during infection with WT Ct L2 and in cdu1::GII mutants expressing wild type Cdu1-Flag or Cdu1I225A-Flag (DUB-) was 43%, 41%, and 43% respectively (Figures 6B and 6D). In HeLa cells infected with cdu1::GII and cdu1::GII strains expressing Cdu1C345A-Flag (DUB-Act-), and Cdu1K268E-Flag (Act-), Golgi repositioning was restricted to only 24%, 23%, and 23% of inclusion perimeters respectively (Figures 6B and 6D). These results confirm that both InaC and Cdu1 are required for efficient repositioning of the Golgi around the Ct inclusion as previously reported (Kokes et al., 2015; Wesolowski et al., 2017; Pruneda et al., 2018; Auer et al., 2020) and that this process is independent of Cdu1’s DUB activity but requires its acetylase activity. Moreover, our results suggest that Cdu1 promotes Golgi repositioning by protecting InaC-mediated redistribution of the Golgi around the Ct inclusion.

CTL0480 promotes recruitment of the myosin phosphatase subunit MYPT1 to the inclusion membrane where it regulates the extrusion of intact inclusions from host cells (Lutter et al., 2013; Shaw et al., 2018). Consistent with the gradual loss of CTL0480 from inclusions in cells infected with the cdu1::GII strain starting at 36 hpi (Figures 3C and 3D) we also observed a complete loss of MYPT1 recruitment to inclusions by 48 hpi (Figures 6C and 6F).

Cdu1, InaC, IpaM, and CTL0480 are required for optimal extrusion of Ct from host cells

In the absence of Cdu1, the levels of InaC, CTL0480, and IpaM decreased late in infection (36 hpi and 48 hpi, Figures 3 and 4) suggesting that a prominent role of Cdu1 is to protect these Incs from degradation late in infection. At the end of its developmental cycle, Chlamydia exits host cells by promoting cellular lysis or by extrusion of intact inclusions (Hybiske and Stephens., 2007). Ct host cell exit by extrusion is an active process requiring a remodeling of the actin cytoskeleton and the function of Inc proteins (Hybiske and Stephens., 2007; Chin et al., 2012; Lutter et al., 2013; Shaw et al., 2018; Nguyen et al., 2018). CTL0480 recruits MYPT1 (an inhibitor of Myosin II motor complexes) to the inclusion membrane which prevents premature extrusion of Ct inclusions and loss of CTL0480 leads to increased rates of extrusion by Ct from infected HeLa cells (Lutter et al., 2013; Shaw et al., 2018) (Figure 7). Actin polymerization is also required for Ct extrusion (Hybiske and Stephens., 2007; Chin et al., 2012) suggesting that InaC dependent recruitment of F-actin to the inclusion may also contribute to optimal Ct extrusion. IpaM localizes to microdomains in the inclusion membrane that are proposed to function as foci for extrusion (Nguyen et al., 2018). Based on these observations, we postulated that Cdu1-mediated protection of CTL0480, InaC, and IpaM regulates the extrusion of Ct inclusions. We quantified the number of extrusions released from infected HeLa cells at 52 hpi and observed a 60% reduction in the number of extrusions in HeLa cells infected with the cdu1::GII strain relative to cells infected with WT L2 (Figures 7A and 7B). Complementation of cdu1::GII with either wild type Cdu1-Flag or Cdu1I225A-Flag (DUB-) restored extrusion production to near wild type levels. Infection of HeLa cells with inaC (inaC::GII, Wesolowski et al., 2017), or ipaM (ipaM::GII, Meier et al., 2020) null strains led to a 75% and 58% reduction in extrusion production respectively (Figures 7A and 7B). Impaired extrusion of inclusions by cdu1::GII, inaC::GII, and ipaM::GII strains was not due to defects in inclusion biogenesis, since all three strains produced similar number of inclusions at 48 hpi relative to cells infected with WT L2 (Supp. Figure 5). We conclude that Cdu1, InaC, and IpaM function to promote optimal extrusion of Ct inclusions from host cells.

Cdu1, InaC, and IpaM are required for optimal extrusion of Ct inclusions from HeLa cells.

(A) Representative images of extrusions isolated from HeLa cell monolayers infected with Ct strains for 52 hours. Scale bar: 200 μm (B) Quantification of the number of extruded inclusions produced by infected HeLa cell monolayers. p values were determined by a student paired t-test. (C) Quantification of the size of extruded inclusions quantified in (B). Extruded inclusions varied in size among cells infected with wild type L2 (average :43 μm), ipaM mutants (average: 56 μm) and cdu1 mutants complemented wild type Cdu1 (average: 52 μm) and Cdu1I225A (DUB-) (average: 60 μm). p values were determined by one-way ANOVAs with a Student-Newman-Keuls post hoc test. Representative images in (A) and quantification of extrusion number in (B) and extrusion size in (C) are based on images obtained from 6 fields across 3 independent biological replicates.

In contrast, infection of HeLa cells with a ctl0480::GII mutant strain led to an increase in the number of extruded inclusions as previously observed (Shaw et al., 2018) (Figures 7A and 7B). Therefore, even though the Cdu1-mediated protection of InaC and IpaM is important for the extrusion of inclusions and cdu1 mutants phenocopy the loss of InaC and IpaM, the phenotypic similarities do not extend to the increased number of extruded inclusions observed in cells infected with the ctl0480::GII mutant strain (Figures 7A and 7B). We infer from these observations that functions for both InaC and IpaM in the extrusion of inclusions are epistatic to CTL0480. Extruded inclusions produced during infection of HeLa cells also varied in size with an average diameter of 40 μm (Figures 7A and 7C). Interestingly, the loss of IpaM and over expression of Cdu1-Flag and Cdu1I225A-Flag (DUB-) shifted the size distribution of extrusions toward larger extrusions (Figures 7A and 7C) suggesting that Ct regulates the size of extruded inclusions through Cdu1.

Discussion

Several Chlamydia effector proteins localize to the inclusion membrane to regulate interactions between the Chlamydia pathogenic vacuole and the host cytoskeleton, organelles, and vesicular trafficking pathways. They also modulate host cell death programs and promote Chlamydia exit from host cells (reviewed in Bugalhão and Mota., 2019). Given the central roles that Incs play in promoting Chlamydia intracellular infection, it is not surprising that they are targeted for inactivation by host cellular defenses. In response, Chlamydia has evolved mechanisms to protect Incs. In this study we show that the acetylase activity of the effector Cdu1 protects itself as well as three independent Incs; InaC, IpaM, and CTL0480, from ubiquitination and degradation (Figure 8). Interestingly, all three Incs play prominent roles in regulating the extrusion of Chlamydia inclusions from host cells (Figure 7). Observations that the encapsulation of Chlamydia within an extruded inclusion enhances survival of Chlamydia within macrophages (Zuck et al., 2017) together with the broad conservation of extrusion as an exit strategy among Chlamydia (Zuck et al., 2016) suggests that this mechanism is important for Chlamydia pathogenesis. Thus, targeting Chlamydia Incs that regulate extrusion for Ub-mediated destruction may be advantageous for the host. In this context, we hypothesize that InaC is targeted for degradation due to its role in promoting F-actin assembly around the inclusion, since F-actin is essential for extrusion (Hybiske and Stephens., 2013; Chin et al., 2012). However, other roles of InaC like its regulation of microtubule scaffolds around the inclusion cannot be ruled out (Wesolowski et al., 2017; Haines et al., 2021). CTL0480 was previously shown to function as an inhibitor of extrusions through its role in modulating the activity of myosin light chain 2 (MLC2) (Lutter et al., 2013; Shaw et al., 2018). IpaM localizes to specialized microdomains in the inclusion membrane which are also sites of enrichment for over 9 inclusion membrane proteins including CTL0480 and MrcA, both of which are required for Chlamydia extrusion (Mital et al., 2010; Lutter et al., 2013; Nguyen et al., 2018). IpaM is a third inclusion membrane protein localizing to microdomains that is required for the optimal extrusion of inclusions and lends support to a potential role for microdomains as foci for extrusion (Nguyen et al., 2018). We also find that the loss of IpaM shifted the size distribution of extrusions towards larger inclusions indicating that Ct constraints the size of extruded inclusions (Figure 7). We speculate that heterogeneity in the size of extrusions might facilitate uptake of some extrusions by innate immune cells at infected mucosal sites. Uptake of Chlamydia by host phagocytes has been proposed to promote Chlamydia LGV dissemination to distal sites in the genital tract to avoid clearance of Chlamydia by other immune cells (Zuck et al., 2017). Protection of IpaM by Cdu1 could therefore be advantageous for Chlamydia EBs within extrusions to escape to distal sites and avoid exposure to localized inflammatory responses.

A model for acetylation mediated protection of the translocated Ct effectors InaC, IpaM, and CTL0480 from degradation.

The human cellular Ub machinery targets C. trachomatis effectors, including the inclusion membrane proteins InaC, IpaM, and CTL0480 for ubiquitination and subsequent protein degradation. C. trachomatis has evolved strategies to counter such defense mechanisms by translocating the effector Cdu1. Cdu1 utilizes its acetylase activity rather than its deubiquitinase activity to protect itself and all three inclusion membrane proteins from being targeted for ubiquitination and degradation. Cdu1 protects InaC in order to promote the recruitment of F-actin scaffolds and Golgi ministacks to the inclusion perimeter and CTL0480 to recruit the Myosin II regulator MYPT1. All three inclusion proteins and Cdu1 also participate in promoting extrusion of C. trachomatis inclusions late in infection. Image created with BioRender.com.

Effectors that modulate the activity of other translocated effectors are referred to as “metaeffectors”, a term coined by Kubori and colleagues after observing that the L. pneumophila effector LubX which functions as an E3 ligase, ubiquitinates the translocated effector SidH leading to its degradation (Kubori et al., 2010). Several other effector-metaeffector interactions have been described in L. pneumophila, Salmonella enterica, and Brucella abortus which regulate the activity of other effectors either directly or indirectly by modifying the same host target or cellular process (Kubori et al., 2010; Neunuebel et al., 2011; Jeong et al., 2015; Urbanus et al., 2016; Smith et al., 2020; Iyer and Das., 2021). In this context Cdu1 functions as an atypical metaeffector since it targets multiple effectors and unlike other effector-metaeffector pairs it is not encoded in the same locus as its targets. We also observed that Cdu1 interactions with InaC, IpaM, and CTL0480 likely occur independently from each other and that the stoichiometry and kinetics of degradation varies for each Inc during infection of host cells. For example, InaC protein levels decrease as early as 24 hpi while IpaM and CTL0480 protein levels decrease at later stages of infection (36 and 48 hpi) (Figure 3). Furthermore, in the absence of Cdu1, we observe almost complete loss of InaC at 36 hpi (Figure 4A) and of CTL0480 at 48 hpi (Figures 3C and 3D), while only a subpopulation of IpaM is destabilized throughout infection (Figures 3B and 4B).

The observation that the DUB activity of Cdu1 was not required to protect InaC, IpaM, and CTL0480 from ubiquitination was unexpected and pointed to a potential role for its acetylation activity instead. Because Cdu1 can interact with InaC, IpaM, and CTL0480 we hypothesized that Cdu1 may protect these Incs through acetylation of lysine residues that are otherwise targeted for ubiquitination. Indeed, we found that all three Incs and Cdu1 are acetylated at lysines in infected cells, however, we were unable to determine if lysine acetylation in all four proteins was dependent on Cdu1s’ Act activity or if these PTMs are protective. Why the DUB activity of Cdu1 is unable to compensate for loss of its Act activity remains unknown. It is possible that Cdu1, like other DUBs, is regulated by PTMs (Komander et al., 2009a). For instance, phosphorylation of human CYLD inhibits its DUB activity towards TRAF2 while phosphorylation of human USP8 inhibits its DUB activity toward EGFR (Reiley et al., 2005; Mizuno et al., 2007). Mass spectrometry analysis of immunoprecipitated Flag tagged Cdu1 expressed in Ct revealed that Cdu1 is phosphorylated at multiple serine and threonine residues (Figure 4H) as previously suggested (Zadora et al., 2019). We identified three PX(S/T)P MAPK phosphorylation consensus sequence motifs in the proline rich domain (PRD) of Cdu1, suggesting that MAPKs may regulate the DUB activity of Cdu1.

Cdu1 homologs are found in multiple Chlamydia species including C. trachomatis, C. muridarum, C. suis, C. psitacci, C. abortus, C. caviae, and C. felis but is notably absent in the genomes of C. pneumoniae and C. pecorum. The acquisition of a second deubiquitinase paralog (Cdu2) has also occurred in C. trachomatis, C. muridarum, and C. suis. In the genomes of all three species, cdu2 resides directly adjacent to cdu1; an arrangement that presumptively arose from a gene duplication event. Cdu2 is a dedicated ULP with deubiquitinating and deneddylating activities (Misaghi et al., 2006; Pruneda et al., 2016). New evidence suggest that both paralogues might not be functionally redundant. The crystal structure of Cdu2 has revealed differences in residues involved in substrate recognition between Cdu1 and Cdu2 and that each paralog might recognize polyUb chains differently (Hausman et al., 2020). The processivity rates for removal of terminal Ub from polyUb chains also differs between both isopeptidases with Cdu2 exhibiting limited trimming of polyUb as compared to Cdu1 (Hausman et al., 2020). Moreover, Cdu2 lacks the proline rich domain found in Cdu1 which might be important for regulation of Cdu1 enzymatic activity. The presence of Cdu2 might also explain the low incidence of human and Ct proteins that were differentially ubiquitinated in the absence of Cdu1 (Figure 1). Whereas several Chlamydia species have acquired either one or two deubiquitinase paralogs, both C. pneumoniae and C. pecorum have not. Instead, both species have acquired an unrelated deubiquitinase (ChlaOTU) belonging to the OTU family of proteases (Makarova et al., 2000; Furtado et al., 2013). Curiously, ChlaOTU is also found in C. psitacci, C. abortus, C. caviae, and C. felis all of which encode only Cdu1 and is absent in C. trachomatis, C. muridarum, and C. suis, all of which encode Cdu1 and Cdu2. It is noteworthy that Chlamydia species have independently acquired deubiquitinases multiple times (Cdu1, Cdu2, ChlaOTU) and that some of these deubiquitinases have evolved into moonlighting enzymes reflecting the diverse strategies adopted by pathogenic Chlamydia as they adapt to their particular niche.

Acknowledgements

We thank LifeSensors and the Duke Proteomics and Metabolomics Shared Resource Center for their proteomics services. We thank the Duke Light Microscopy Core Facility for microscopy services. We also thank Marcela Kokes for generating the IncA-Flag constructs used in this study. This work was supported by NIH grants AI140019 to R.J.B and AI134891 to R.H.V.

Author contributions

R.J.B and R.H.V designed the study. R.J.B wrote the manuscript with input from all listed authors. L.D generated the cdu1::GII strain. M.K verified generation of the cdu1::GII mutant strain and showed interaction of Cdu1-GFP variants with InaC, IpaM, and CTL0480 in transfected HEK cells. R.J.B performed all other experiments. B.S.S and J.N.P shared reagents. R.J.B and R.H.V proofed the manuscript.

Declaration of interests

R.H.V is a founder of Bloom Sciences (San Diego, CA), which is a microbiome therapeutics company. Findings reported in this study are unrelated to the work being performed with Bloom Sciences.

Material and Methods

Key Resource Table

Resource Availability

Materials availability

  • All newly generated materials associated with this study will be freely available upon request.

Data and code availability

  • Proteomics data will be deposited at PRIDE (Proteomics Identification Database-EMBL-EBI) and publicly available at the time of publication. Accession numbers will be listed in the key resources table.

  • Original data used for microscopy and western blots in this study will be deposited online (Mendeley data) and made available at the time of publication.

  • Any additional information required to reanalyze the data reported in this paper is available upon request.

  • This paper does not report original code.

Experimental Model And Subject Details

Cell lines

Vero (CCL-81; RRID:CVCL_0030), HeLa (CCL-2; RRID:CVCL_0059), and HEK293T (CRL-3216; CVCL_0063) cells were purchased from ATCC and cultured in High Glucose Dulbecco’s Modified Eagle’s Medium supplemented with L-glutamine, sodium pyruvate (DMEM; Gibco) and 10% fetal bovine serum (FBS; Sigma-Aldrich). Cells were grown at 37°C in a 5% CO2 humidified incubator. Vero, Hela, and HEK293T cells were derived from females. All three cell lines have been authenticated by the Duke Cell Culture and DNA analysis facility.

Chlamydia strains and propagation

Chlamydia strains used in this study are listed in the Key Resources Table. Ct strains were propagated in Vero cells and harvested by osmotic lysis at 48 hours post infection. Following lysis extracts were sonicated and bacteria pelleted by centrifugation at 21,000 x g. Bacteria were resuspended in SPG storage buffer (75g/L sucrose, 0.5 g/L KH4HPO4, 1.2 g/L Na2HPO4, 0.72 g/L glutamic acid, pH 7.5) and stored as single use aliquots at -80°C.

Method Details

Chlamydia infections

Chlamydia infections were synchronized by centrifugation (2,500 x g for 30 minutes at 10°C) onto HeLa cell monolayers and incubated for the indicated times. Co-infections were performed by infecting HeLa cell monolayers at a 1:1 ratio using MOIs of 2 for each co-infecting strain.

Insertional mutagenesis of CTL0247 (cdu1)

Primer sequences for TargeTronTM mediated mutagenesis of the LGV L2 434 Bu cdu1 (CTL0247) ORF were designed at the TARGETRONICS, LLC web portal (www.targetrons.com). IBS1/2, EBS1/delta, and EBS2 primers (primer sequences are listed in S. Table 10) were used in a PCR reaction to generate homing sequences for TargeTronTM integration between nucleotides 635 and 636 of the cdu1 ORF using a TargeTronTM gene knockout system (Sigma-Aldrich; TA0100) according to the manufacturer instructions. Homing sequences were gel purified, digested with HindIII and BsrGI, and ligated into HindIII and BsrGI digested pDFTT3-aadA (Lowden et al., 2015). Ligations were transformed into E. coli DH5α, clones isolated, and cdu1 redirected pDFTT3-aadA plasmids identified by restriction digest and verified by Sanger sequencing (Eton Bioscience) using a T7-promoter specific primer. The resulting plasmid was transformed into a C. trachomatis LGV L2 434 Bu strain and transformants selected with 150 μg/mL spectinomycin and plaque purified as previously described (Kędzior and Bastidas., 2019). Insertion of the GII aadA intron at the cdu1 locus was verified by PCR analysis (S. Figure 1) using primers that amplify amplicons spanning the cdu1::GII 5’ (RBP409 and RBP436) and 3’ (RBP468 and RBP118) junctions, the cdu1 CDS (RBP409 and RBP118), and the aadA CDS (RBP512 and RBP513). Primer sequences are listed in S. Table 10. Loss of Cdu1 protein was verified by western blot and indirect immunofluorescence analysis (S. Figures 2A and 2B).

Analysis of cdu1 and cdu2 transcription by RT-PCR

Confluent HeLa cell monolayers (2.9 x 106 cells/infection) were infected with wild type L2 434 Bu or L2 cdu1::GII aadA strains. At 24 hpi, total RNA was isolated with a Qiagen RNeasy kit (Qiagen; 74004) according to the manufacturer instructions. Total RNA was treated twice with DNAse I (NEB; M0303S) and used for cDNA synthesis using a SuperScript IV Reverse Transcriptase kit (Thermo Fisher Scientific; 18090010). cDNAs synthesized with and without reverse transcriptase were used as templates for PCR analysis (S. Figure 2C) using primers that amplify amplicons spanning the cdu1 (CTL0247_F and CTL0247_R) and cdu2 ORFs (CTL0246_F and CTL0247_R), and the intergenic junction between the cdu1 and cdu2 ORFs (CTL0246-0247_F and CTL0246-0247_R). Primer sequences can be found in the S. Table 10.

TUBE1 based global ubiquitin profiling

Confluent HeLa cell monolayers (5.04 x 106 cells/infection) were mock infected or separately infected with WT LGV L2 434 Bu or a L2 cdu1::GII aadA strain at MOIs of 3. At 24 hpi cells were collected and spun down (700 x g for 10 minutes), frozen at -80°C and shipped on dry ice to LifeSensors (Malvern, PA) for quantitative TUBE1-based Mass Spectrometry Analysis. Cell pellets from three independent biological replicates were sent to LifeSensors for analysis. Cell were subsequently lysed in lysis buffer (50 mM Tris-HCL, pH 7.5, 150 mM NaCl, 2 mM EDTA, 1% NP-40, 10% glycerol, 1% Sodium Deoxycholate) supplemented with a protease inhibitor cocktail, the DUB inhibitor PR-619 (Sigma-Aldrich; SML0430), and the proteasomal inhibitor MG-132 (Sigma-Aldrich; 474791). Lysates were clarified by high-speed centrifugation (14,000 x g, 10 minutes, 4°C) and supernatants containing 5 mg of protein were equilibrated with magnetic TUBE-1 (LifeSensors; UM401M) and incubated overnight at 4°C under rotation. TUBEs were isolated with a magnetic stand and washed sequentially with PBST and TUBE wash buffer. Poly-ubiquitinated and associated proteins were eluted with TUBE elution buffer. Eluted supernatants were neutralized with neutralization buffer and loaded onto SDS-gels and run until SDS buffer reached 0.5 cm into the gel. Gels were stained with Coomassie Blue and lanes excised, reduced with TCEP, alkylated with iodoacetamide, and digested with Trypsin (Trypsin Gold, Mass Spectrometry Grade) (Promega; V5280). Tryptic digests were analyzed using a 150 min LC run on a Thermo ScientificTM Q Exactive HF OrbitrapTM LC-MS/MS system. MS data was searched against the UniProt human database (UniProt; Proteome ID: UP000005640) and the Chlamydia trachomatis L2 434 Bu reference database (NCBI:txid47472) using MaxQuant 1.6.2.3 (https://www.maxquant.org). Proteins, peptides, and site identification was set to a false discovery rate of 1%. N-terminal acetylation, Met oxidation, and diGly remnant on lysine residues was also identified. All peptides and proteins identified can be found in S. Table 1. The intensities (sum of all peptide MS peak areas for a protein or ubiquitinated peptide) for each protein and ubiquitinated peptide across all three biological replicates were used to determine mean intensities and to calculate p-values based on one-way student t-tests. Volcano plots of mean intensities vs. p-values were generated with VolcaNoseR (Goedhart and Lujsterburg., 2020) (https://huygens.science.uva.nl/VolcaNoseR/) and used to identify significantly enriched proteins and ubiquitinated peptides. Data used to generate each Volcano plot can be found in S. Table 2. Pathway enrichment analysis was performed with Metascape (Zhou et al., 2019) (https://metascape.org/gp/index.html#/main/step1) and DAVID bioinformatic resources (Huang et al., 2009a; Huang et al., 2009b) (https://david.ncifcrf.gov/tools.jsp).

Inhibitors, antibodies, and western blots

MG-132 (25 μM) (Sigma-Aldrich; 474791) was added to infected cell monolayers 5 hours prior to extract preparations. Recombinant Cdu1 protein (LGV L2 434 Bu, amino acids 71-401) was generated as previously described (Pruneda et al., 2018) and kindly provided by Jonathan Pruneda (Oregon Health and Science University, OR). Recombinant Cdu1 protein was used to generate antibodies in immunized New Zealand White rabbits. Cdu1 antisera was pre-adsorbed with crude cell extracts from HeLa cells infected with a cdu1::GII aadA strain. Pre-adsorbed antisera was used for western blot analysis at a 1:500 dilution in a solution containing 5% BSA supplemented with crude extracts from HeLa cells infected with a cdu1::GII aadA strain (0.1mg/mL total protein). Antibodies, antibody dilutions, and antibody diluents used in this study are listed in S. Table 9. For western blot analysis, lysates from infected HeLa cell monolayers (2.4 x 106 cells) were prepared by incubating cell monolayers with boiling hot 1% SDS lysis buffer (1% SDS, 100 mM NaCl, 50 mM Tris, pH 7.5). Lysates were collected, briefly sonicated, and total protein concentration measured with a DCTM Protein Assay Kit (BIO-RAD; 5000111). 8 μg (Slc1 and alpha Tubulin blots) and 25 μg of total protein lysates (all other blots) were loaded onto 4-15% Mini-PROTEAN and TGX Stain FreeTM Protein Gels (BIO-RAD; 4568084), transferred to PVDF membranes (BIO-RAD; 1620177), blocked with 5% Milk/TBSt, and incubated with primary antibodies overnight at 4°C. Protein signals were detected with Goat anti mouse (H+L) IgG (ThermoFisher scientific; 31430) or Goat anti rabbit (H+L) IgG HRP (ThermoFisher scientific; 31460) conjugated secondary antibodies (1:1000 in 5% Milk/TBSt) and SuperSignal West Femto HRP substrate (ThermoFisher scientific; 34096). Antibody bound membranes were imaged with a LI-COR imaging system (LI-COR, Inc.). Varying amounts of protein extracts were used to determine the linear range of detection for InaC, IpaM, and Slc1 antibodies prior to quantification of western blot images with LI-COR imaging software (data not shown).

Immunofluorescence microscopy

HeLa cells were grown on coverslips to 50% confluency (0.1 x 105 cells) and infected at MOIs of 0.6. At indicated times, infected cells were separately fixed with ice cold Methanol or with warm PBS containing 4% formaldehyde for 20 minutes. After fixative removal, cells were washed with PBS and formaldehyde fixed cells were incubated either in 5% BSA/PBS supplemented with 0.1% Triton X-100 or in 5% BSA/PBS supplemented with 0.05% Saponin for 30 minutes with gentle rocking. Following washing with PBS, Methanol fixed cells were incubated with primary antibodies diluted in 5% BSA/PBS and formaldehyde fixed cells were incubated with primary antibodies diluted in 5% BSA/PBS supplemented with 0.1% Triton X-100 or 0.05% Saponin for 1 hour with gentle rocking. Dilutions for each antibody used can be found in S. Table 9. Methanol fixed cells were washed with PBS and incubated with secondary antibodies diluted in 5% BSA/PBS and supplemented with Hoechst 33342 (2 μg/mL) (ThermoFisher Scientific; H3570). Formaldehyde fixed cells were washed and incubated with 5% BSA/PBS supplemented with 0.1% Triton X-100 and Hoechst or 0.05% Saponin and Hoechst for 1 hour protected from light and with gentle rocking. For detection of F-actin, Phalloidin conjugated to Alexa FluorTM 488 (1:5000) (Act-Stain 488 Phalloidin; Cytoskeleton Inc.; PHDG1) was added for the last 20 minutes of incubation with the secondary antibodies. Coverslips were transferred to glass slides, mounted with 10 μL of Vectashield (Vector Labs; H-1000) and incubated over night at room temperature prior to imaging. Secondary antibodies used were goat anti-mouse (H+L) IgG (ThermoFisher scientific; A-11001 and A-21235) and goat anti-rabbit (H+L) IgG (ThermoFisher scientific; A-11008 and A-21244) conjugated to Alexa FluorTM 488 and Alexa FluorTM 647. All of the antibodies used for indirect immunofluorescence analysis were analyzed under all three staining conditions (Methanol, Formaldehyde/Triton X-100, and Formaldehyde/Saponin).

Representative images were acquired with an inverted confocal laser scanning microscope (Zeiss 880) equipped with an Airyscan detector (Hamamatsu) and with diode (405 nm), argon ion (488 nm), double solid-state (561 nm) and helium-neon (633) lasers. Images were acquired with a 63x C-Apochromatic NA 1.2 oil-objective (Zeiss). Images acquired in Airyscan mode were deconvoluted using automatic Airyscan processing in Zen software (Zeiss). Image acquisition was performed at the Light Microscopy Core Facility at Duke University. Images used for quantification were captured in an inverted microscope (Ti2-Nikon instruments) equipped with an ORCA Flash 4.0 V3 sCMOS camera (Hamamatsu) and a SOLA solid-state white light illuminator (Lumencro). Images were acquired using a 60x Plan Apochromatic NA 1.40 oil objective. All images were opened with ImageJ (Schneider et al., 2012) and only linear adjustments were made to fluorescence intensity for the entire image. Images were exported as TIFFs and compiled with Adobe suite software (Illustrator).

Vector construction and C. trachomatis transformation

Constructs used in co-transfection experiments

Mammalian vectors expressing Cdu1-GFP constructs were kindly provided by Jonathan Pruneda (Oregon Health and Science University, OR). Briefly, geneblocks encoding full length Cdu1 (LGV L2 434 Bu, CTL0247) (amino acids 1-401), Cdu1 lacking its transmembrane domain (amino acids 71-401), and Cdu1 lacking its catalytic domain (amino acids 1-130) were generated and inserted into the pOPIN-GFP vector (Berrow et al., 2007) by In-FusionTM cloning (Takara Bio; 638947), resulting in Cdu1 constructs with a C-terminal eGFP-His tag preceded by a 3C protease cleavage site. The Flag-InaC mammalian expression vector was derived from a GatewayTM entry clone containing the C. trachomatis Serovar D/UW-3/CX CT813 (inaC) ORF (amino acids 41-264) obtained from a C. trachomatis ORFeome library (Roan et al., 2006). The entry vector was used as a donor plasmid for GatewayTM based transfer into a modified pcDNATM DEST53 (ThermoFisher Scientific; 12288015) vector in which the cycle 3 GFP ORF was removed. A NEB Q5-Site Directed Mutagenesis Kit (New England Biolabs; E0554S) was used to introduce a 3X Flag epitope tag at the N-terminus of the CT813 ORF and a stop codon at the end of the CT813 ORF. L2 ipaM (CTL0476), L2 CTL0480, and L2 cpoS (CTL0481) ORFs were PCR amplified from cell lysates derived from Vero cells infected with wild type L2 LGV 434 Bu with primers containing attB sequences (primers ipaM forward, ipaM reverse, CTL0480 forward, CTL0480 reverse, cpoS forward, and cpoS reverse). Primer sequences can be found in the S. Table 10. PCR amplicons were used as donors for GatewayTM BP ClonaseTM based transfers into the donor vector pDONRTM221 (ThermoFisher Scientific; 12536017) to generate entry plasmids. Entry plasmids were used to transfer ipaM, CTL0480, and cpoS into the GatewayTM destination vector pcDNATM3.1/nV5-DEST (ThermoFisher Scientific; 12290010) by GatewayTM LR ClonaseTM based reactions. The resulting mammalian expression vectors express IpaM, CTL0480, and CpoS with V5-epitopes fused to their N-terminus.

pBOMB4-MCI based plasmids

CTL0480 and ipaM ORFs were amplified by PCR from cell extracts derived from Vero cells infected with wild type LGV L2 434 Bu. The CTL0480 ORF, 149 b.p of upstream sequence, and a 3X FLAG epitope was amplified by PCR using primers RBP628 and RBP629. The CTL0476 (ipaM) ORF, 400 b.p of upstream sequence, and a 3X FLAG epitope was amplified with primers RBP623 and RBP624. CTL0480 and ipaM amplicons were digested with Not1 and Pst1 and cloned into Not1 and Pst1 digested pBOMB4-MCI (Bauler and Hackstadt., 2014) to generate pBOMB4-MCI_CTL0480-3X Flag and pBOMB4-MCI_IpaM-3X Flag plasmids respectively. pBOMB4-MCI_Cdu1-3XFlag plasmids were generated by PCR amplification of 175 b.p of genomic sequence directly upstream of the L2 434 Bu CTL0247 (cdu1) ORF and the entire cdu1 ORF tagged with a C-terminal 3X Flag epitope tag (primers RBP460 and RBP461). PCR amplicons were generated from gradient purified LGV L2 434 Bu EBs and cloned into a pCRTM-Blunt II TOPO vector using a Zero BluntTM TOPOTM PCR cloning Kit (ThermoFisher Scientific; K2800J10) according to the manufacturer instructions. Cdu1 catalytic variants were generated with a NEB Q5-Site Directed Mutagenesis Kit (New England Biolabs; E0554S) using the cdu1p-cdu1-3X Flag construct cloned into pCRTM-Blunt II TOPO as a template and following the manufacturer instructions. Primers for introducing base pair changes were designed on the NEBaseChanger website (https://nebasechanger.neb.com). The Cdu1C345A variant was generated by changing the TGC codon located at positions 1033-1035 in the cdu1 ORF to GCT (primers RBP525 and RBP526). The Cdu1I225A variant was generated by substituting the ATC codon located at positions 673-675 for GCT (primers RBP527and RBP528). The Cdu1K268E variant was generated by introducing an A802G base pair substitution (primers RBP529 and RBP530). Wild type cdu1p-cdu1-3XFLAG and all three cdu1 variants were digested with Not1 and Pst1 and ligated into Not1 and Pst1 digested pBOMB4-MCI (Bauler and Hackstadt., 2014). Primer sequences can be found in the S. Table 10.

p2TK2_SW2-inaC-3XFlag

The CTL0184 (inaC) ORF and 250 b.p of upstream sequence was amplified by PCR from cell extracts derived from Vero cells infected with L2 434 Bu and cloned into the p2TK2_SW2 vector (Agaisse and Derré., 2013). A NEB Q5-Site Directed Mutagenesis Kit (New England Biolabs; E0554S) was used to insert a 3X FLAG epitope sequence at the C-terminus (stop codon removed) of the CTL0184 ORF to generate the p2TK2_SW2-InaC-3XFlag plasmid.

pBOMB4-MCI based plasmids and p2TK2_SW2-InaC-3X Flag plasmids were transformed into corresponding Chlamydia strains, and transformants were selected with 10 U/mL Penicillin G and plaque purified as previously described (Kędzior and Bastidas., 2019). All primer sequences and plasmids generated in this study are listed in S. Table 10 and Key Resources Table.

Subcellular fractionation

HeLa cells (2.16 x 107 cells/strain) seeded in six well plates were mock infected or infected with Chlamydia L2 434 Bu strains. At indicated time points cells were washed with ice-cold PBS, collected in ice-cold PBS with a cell scraper, and transferred to 15 mL conical tubes. Cell suspensions were centrifuged at 500 x g for 5 minutes at 4°C and cell pellets were resuspended in 400 μL of ice-cold subcellular fractionation buffer (20 mM HEPES (pH 7.4), 10 mM KCl, 2mM MgCl2, 1mM EDTA, 1mM EGTA) supplemented with 1mM DTT and a 1x cOmplete Mini-EDTA free protease inhibitor cocktail (Sigma-Aldrich; 11836170001). Cells were incubated on ice for 20 minutes and lysed with 30 strokes of a Dounce homogenizer. Cell lysates were sequentially centrifuged twice at 720 x g for 5 minutes at 4°C to remove intact nuclei. Supernatants were centrifuged at 10,000 x g for 5 minutes at 4°C and the heavy membrane (inclusion) fraction was recovered and resuspended in IP lysis buffer (25 mM Tris-Hcl, 150 mM NaCl, 1% NP-40, 5% Glycerol) supplemented with 1mM PMSF and a 1x cOmplete Mini-EDTA free protease inhibitor cocktail (Sigma-Aldrich; 11836170001).

Immunoprecipitations (IPs)

Transfections and GFP-immunoprecipitations: HEK 293T cell monolayers (1.32 x 107 cells/transfection) seeded in 10 cm cell culture dishes pre coated with poly-L-Lysine (Sigma-Aldrich; P4707) were grown to 50% confluency and transfected with 10 μg of each plasmid used per co-transfection in a 1.5 to 1 ratio of jetOPTIMUS (Polyplus; 101000051) transfection reagent to total plasmid DNA according to the manufacturer instructions. At 24 hpi, transfected cells were lysed in IP lysis buffer (described above) supplemented with 1mM PMSF and 1x cOmplete Mini-EDTA free protease inhibitor cocktail (Sigma-Aldrich; 11836170001). Lysates were transferred to Eppendorf tubes, sonicated, and cleared by centrifugation (21,000 x g, 15 minutes, 4°C). Supernatants containing 2 mg of total protein were incubated with magnetic GFP-Trap® agarose (Proteintech; gta) for 1 hour at 4°C with rotation. Beads were washed according to manufacturer instructions and immunoprecipitated proteins eluted with 2X Laemmli sample buffer.

Flag immunoprecipitations

Mock and infected HeLa cell monolayers (1.44 x 107 cells) grown in six well plates were lysed in IP lysis buffer (described above) and transferred to Eppendorf tubes. Lysates were sonicated and cleared by centrifugation (21,000 x g, 15 minutes, 4°C). Supernatants containing 2 mg of total protein were pre-cleared by incubating with Protein A/G PLUS-Agarose (Santa Cruz Biotechnology; sc-2203) for 30 minutes at 4°C followed by sedimentation of agarose resins by centrifugation. Supernatants were incubated with M2-anti Flag mouse mAb (1:400) (Sigma-Aldrich; F1804) overnight at 4°C with rotation followed by incubation with Protein A/G PLUS-Agarose for 3 hours at 4°C with rotation. Agarose resins were sedimented and washed according to manufacturer instructions and immunoprecipitated proteins eluted with 50 μL of 100 μg/mL 3xFLAG peptides (APExBIO; A6001).

Acetylated lysine immunoprecipitations

Heavy membrane (inclusion) subcellular fractions isolated from infected HeLa cells and containing 1 mg of total protein were pre-cleared by incubating with Protein A/G PLUS-Agarose (Santa Cruz Biotechnology; sc-2203) for 30 minutes at 4°C followed by sedimentation of agarose resins by centrifugation. Supernatants were incubated with an anti-acetylated lysine rabbit antibody (Cell signaling; #9441) (1:100) overnight at 4°C with rotation followed by incubation with Protein A/G PLUS-Agarose for 3 hours at 4°C with rotation. Agarose resins were sedimented and washed according to manufacturer instructions and immunoprecipitated proteins eluted with 50 μL of 2X Laemmli sample buffer.

Input (8 μg of total protein for Slc1 and alpha Tubulin blots, and 25 μg for all other blots) and immunoprecipitates (GFP, Flag, proteins acetylated at lysines) were loaded onto 4-15% Mini-PROTEAN and TGX Stain FreeTM Protein Gels (BIO-RAD; 4568084), transferred to PVDF membranes (BIO-RAD; 1620177), blocked with 5% Milk/TBSt, and incubated with primary antibodies overnight at 4°C. Protein signals were detected with Goat anti mouse (H+L) IgG (ThermoFisher scientific; 31430) or Goat anti rabbit (H+L) IgG HRP (ThermoFisher scientific; 31460) conjugated secondary antibodies (1:1000 in 5% Milk/TBSt) and SuperSignal West Femto HRP substrate (ThermoFisher scientific; 34096). Antibody bound membranes were imaged with a LI-COR imaging system (LI-COR, Inc.).

Identification Cdu1 lysine acetylation and phosphorylation sites by LC-MS/MS

HeLa cell monolayers (8.64 x107 cells/strain) grown in six well plates were infected with a wild type L2 strain transformed with empty pBOMB4-MCI (Bauler and Hackstadt., 2014) plasmid or with a wild type L2 strain transformed with a pBOMB4-MCI_cdu1-3X Flag plasmid. At 24 hpi infected cells were lysed and Flag tagged Cdu1 was immunoprecipitated as described above. Flag eluates from 3 independent biological replicates were sent to the Proteomics and Metabolomics Shared Resource Facility at Duke University for quantitative LC-MS/MS analysis. Samples were spiked with undigested casein, reduced with 10 mM dithiothreitol, and alkylated with 20 mM iodoacetamide. Eluates were then supplemented with 1.2% phosphoric acid and S-Trap (Protifi) binding buffer (90% Methanol, 100 mM TEAB). Proteins were trapped on the S-Trap, digested with 20 ng/μL Trypsin (Trypsin Gold, Mass Spectrometry Grade) (Promega; V5280), and eluted with 50 mM TEAB, 0.2% FA, and 50% ACN/0.2% FA. Samples were lyophilized and resuspended in 1% TFA/2% acetonitrile containing 12.5 fmol/μL yeast alcohol dehydrogenase. Quantitative LC/MS/MS was performed using a nanoAcquitiy UPLC system (Waters Corp) coupled to Thermo ScientificTM OrbitrapTM Fusion Lumos high resolution accurate mass tandem mass spectrometer via a nanoelectrospray ionization source. Data was analyzed with Proteome Discoverer 2.3 (Thermo Fisher ScientificTM.) and MS/MS data searched against the Chlamydia trachomatis LGV L2 434 Bu reference database (NCBI:txid47472). Cdu1-Flag MS/MS data was analyzed with Mascot software (Matrix Science) using Trypsin/P specificity for N-terminal acetylation, lysine acetylation, lysine Ub, and S/T/Y phosphorylation identification. Analysis identified multiple acetylated and phosphorylated Cdu1 peptides and no Cdu1 ubiquitinated peptides. Data was viewed in Scaffold with Scaffold PTM (Scaffold Software).

Isolation and imaging of Chlamydia inclusion extrusions

HeLa cell monolayers (1.2 x 106 cells) were infected with Ct strains at MOIs of 0.8. At 48 hpi, infected monolayers were imaged using an EVOS FL Cell Imaging System (ThermoFisher Scientific) equipped with a 20x/0.4 NA objective and a CCD camera. Following imaging, growth media was removed, cell monolayers washed with fresh growth media, and monolayers incubated for an additional 4 hours at 37°C. At 52 hpi growth supernatants were collected and transferred to Eppendorf tubes. Extrusions were enriched by centrifugation (1,500 rpm, 5 min) and pellets (not always visible) containing extrusions were resuspended in 30 μL of 4% Formaldehyde/PBS supplemented with Hoechst (2 μg/mL) and 0.2% Trypan Blue Solution (Gibco, 0.4%). Extrusions were analyzed by plating 10 μL drops on a glass side (without coverslips) and immediately imaged using an EVOS FL Cell Imaging System (ThermoFisher Scientific) equipped with a 20x/0.4 NA objective and a CCD camera. Intact extrusions were identified based on morphology, lacking nuclei, and being impermeable to trypan blue. Images were opened in ImageJ (Schneider et al., 2012) and enumeration of inclusions and extrusions was performed manually. The sizes of individual extrusions and inclusions were determined by manually tracing a line around the perimeter of each extrusion and inclusion in ImageJ and measuring perimeter length. All measurements were exported to Microsoft Excel. Data plots and statistical analyses were done with Prism 9 (GraphPad) software. Datasets were analyzed for significance using a paired student t-test.

Image analysis

Line scan profiles of Cdu1 co-localization with inclusion membrane proteins was performed with ImageJ (Schneider et al., 2012) by tracing a line through regions of interest and plotting fluorescent signal intensities with the Plot Profile function. Localization of CTL0480, recruitment of MYPT1, and association of Ub with Ct inclusions was performed manually from maximum projections in ImageJ. Assessment of F-actin recruitment to Ct inclusions was performed manually in ImageJ by projecting four to five sections in order to capture the entire inclusion. Redistribution of Golgi around the Ct inclusion was measured in ImageJ from maximum projections. The perimeters of individual inclusions were manually traced, and lengths measured. The length of dispersed Golgi was measured by tracing and measuring the length of the GM130 signal directly adjacent to each inclusion. Dispersed Golgi length was divided by inclusion perimeter length. All measurements were exported to Microsoft Excel for quantification. Data plots and statistical analyses were done with Prism 9 (GraphPad).

Quantification And Statistical Analysis

Quantifications were generated from three independent experiments and measurements derived from blinded images. Data plots and statistical analyses were done with Prism 9 (GraphPad) software. Datasets were analyzed for significance using a paired student t-test or one-way ANOVAs with a Student-Newman-Keuls post hoc test. Data graphs show means and error bars represent standard error. p-values less than 0.05 are defined as statistically significant. The indicated statistical test for each experiment can be found in the figure legends.

Figures

TargeTron mediated disruption of the L2 cdu1 ORF.

(A) Depiction of the cdu1 ORFs in WT L2 and cdu1::GII strains with corresponding diagnostic PCR amplicons. (B) PCR analysis of the disrupted cdu1 ORF in the cdu1::GII strain. *Amplification of cdu1 ORF in cdu1::GII strain was unsuccessful due to size of amplicon.

Generation of a cdu1 null strain in C. trachomatis (L2).

(A) Western blot analysis of protein lysates derived from HeLa cells infected with a Ct L2 434 Bu parental strain (L2) and an L2 cdu1::GII aadA (cdu1::GII) strain for 24 hours. Blots were probed with antibodies raised against recombinant Cdu1 (amino acids 71-401), Ct RpoB (bacterial RNA polymerase), and human α tubulin. (B) HeLa cells infected with L2 transformed with pBOMB4-MCI empty vector (L2 pBOMB) and cdu1::GII transformed with empty pBOMB4-MCI (cdu1::GII pBOMB) were fixed at 24 hpi and stained with Cdu1 antisera. Cdu1 signal is depicted in green, Ct cells expressing mCherry (encoded in pBOMB4-MCI vector) are shown in red, and DNA stained with Hoechst in blue. Scale bar: 10μm. (C) Total RNA isolated from HeLa cells infected with L2 and cdu1::GII strains for 24 hours was used to synthesize cDNAs which served as templates for PCR analysis of cdu1, cdu2, and cdu1-cdu2 bicistronic expression. RT+: Total RNA + reverse transcriptase. RT-: Total RNA without reverse transcriptase. gDNA: Genomic DNA.

Proteins co-precipitating (non ubiquitinated) with human and Ct proteins enriched by TUBE 1 pulldowns.

(A) Volcano plots (pairwise comparisons) of human proteins (non-ubiquitinated) that co-precipitate with TUBE1 bound proteins in mock infected HeLa cells and HeLa cells infected with L2 or cdu1 null strains (24 hpi). (B) Volcano plots (pairwise comparisons) of Ct TUBE1 co-precipitating proteins (non-ubiquitinated) identified in mock infected HeLa cells and HeLa cells infected with L2 or cdu1 null strains (24 hpi).

Pathway enrichment analysis of human and Ct proteins that co-precipitate (non ubiquitinated) with TUBE1 bound proteins.

(A) Metascape functional enrichment analysis of human co-precipitating proteins in mock infected HeLa cells. (B) Metascape functional enrichment analysis of human co-precipitating proteins in cdu1 null infected HeLa cells. (C) DAVID pathway enrichment analysis of Ct co-precipitating proteins enriched in infected (L2 and cdu1 null) HeLa cells.

The number of inclusions and the size of inclusions in cdu1 null, inaC null, ipaM null, and ctl0480 null strains are similar across each strain

(A) Representative images of inclusions in infected HeLa cell monolayers at 48 hpi. Scale bar: 200 μm (B) Quantification of the number of inclusions in infected HeLa cell monolayers. (C) Quantification of inclusion size.