Toxoplasma Infection Induces Release of Extracellular Vesicles from Host Cells that Decrease NE
Extracellular vesicles (EVs) have emerged as an important facet of host-pathogen interactions and were investigated for their relationship to the NE suppression (30, 31). EVs were purified from infected catecholamine-producing dopaminergic/noradrenergic cell cultures (henceforth termed noradrenergic cells) by step-wise ultracentrifugation followed by purification by sucrose-gradient fractionation (32, 33) (Fig. 1A). The EVs appear to be of host cell origin based on size, morphology, sedimentation rate, and exosome markers (34). The infected culture EVs contained mammalian exosomal protein markers CD81 and EpCAM (Fig. S1 and proteomic analysis, data not shown) as previously found with PC12 cells, as well as CD63, ICAM, FLOT-1, and TSG101 being detectable. The protein profiles were similar for EVs from infected and uninfected cultures (Fig. S1) and the size consistent with exosomes (Figs. 1B, 1C). Based on these properties, in this paper we have termed these Toxoplasma-induced neuronal host-derived extracellular vesicles (TINEV). Similar yields of released EVs were isolated regardless of infectious status in our experiments (5.6±1.4 µg/ml and 3.9±0.49 µg/ml, respectively; p = 0.438) Infection naive noradrenergic cells were treated with TINEVs and expression of DBH was measured. Noradrenergic cells in medium with commercial exosome-depleted serum served as control. DBH expression was significantly down-regulated (190 ± 67-fold, p = 0.006 in Fig. 1D) in cultures treated with TINEVs.
Next the in vivo effect of the TINEVs was examined. Intracranial injection of TINEVs into the locus coeruleus (LC), the central region of noradrenergic neurons in the brain, of adult rats was performed and DBH expression was measured by RT-qPCR (Fig. 1E, Fig. S2). TINEV injection induced a decrease in DBH mRNA (62 ± 43-fold, p = 0.0079) in the pons/LC region (Fig. 1F) compared to treatment with EVs from uninfected cultures. No difference was found in DBH mRNA levels in the mid-brain and prefrontal cortex for rats injected TINEVs, although low numbers of noradrenergic neurons are found in these regions. Health parameters (body weight, appearance and food intake) were normal in treated animals and expression of a housekeeping gene was unaltered in the brain with treatment (Fig. S3). Further, DBH expression in the adrenal glands was unaffected by the treatments (data not shown). The decrease in DBH mRNA observed with intracranial TINEV injection was similar to that observed in chronic infections (11).
Levels of NE and DBH mRNA have been found decreased with chronic T. gondii infection (11). A direct correlation between DBH expression and NE level (Fig. S4; correlation coefficient 0.81, p = 0.014) was observed in rodent brains. Hence DBH expression can be used as a correlate of brain NE, as found in other studies (35). Based on the above findings, T. gondii infection induces cells to release TINEVs that are able to down-regulate DBH expression and cause a widespread decrease in brain NE.
Strategy that Identified Transcriptional Gene Silencing and Epigenetic Change through Paracrine Signalling as the Mechanism of DBH Silencing
Initially, it was considered that several different factors could explain the disproportionate decrease in DBH expression (and hence NE) relative to the small percent of parasitised cells during chronic infection (e.g. neuroimmune responses). As decreased NE and DBH down-regulation (relative to other genes) has been observed in vitro with noradrenergic cell lines, mechanisms other than the host immune system are involved. A nuclear run-on assay was performed to assess whether the DBH down-regulation is at the transcriptional or post-transcriptional level. De novo transcription in nuclei isolated from infected cell cultures was measured by immunocapture of incorporated biotin-UTP. Lower amounts of nascent DBH mRNA (relative to standards) were found in host cell nuclei from infected than uninfected noradrenergic cell cultures (21 ± 1.6-fold, p = 0.00062) (Fig. 2A). Similarly, de novo transcription of DBH was downregulated in infected human noradrenergic cells (21 ± 1.7-fold, p = 0.032) (Fig. 2B). Hence, infection induced transcriptional gene silencing (TGS) of DBH in human and rat noradrenergic cells.
We examined the uninfected cells in cultures to determine whether those cells exposed to infected cells were also suppressed in NE as this could help explain the large change in expression observed. A transwell system was used that permits uninfected cells to be exposed to T. gondii-infected cell products. This method was chosen because it differentiates between diffusible signals and parasites injecting components into cells without invasion, as was observed with a Cre/loxP assay in infected mouse brains (36). DBH expression was measured in uninfected rat noradrenergic cells in the bottom reservoir of the transwell system with the top reservoir containing an infected culture (Fig. 2C). DBH expression in the cells exposed to infected cultures was found to be down-regulated (7.9 ± 2.8-fold, p = 0.02) suggesting that a transmissible factor was released from infected cells that was subsequently identified as EVs. In contrast, noradrenergic cells that were exposed to cultures containing heat-killed T. gondii (‘mock-infected’) were unchanged in DBH expression (Fig. 2D). Exposure of a human neuronal cell line to infected cells in the transwell system induced a larger decrease in DBH expression (37 ± 16-fold, p = 0.0038) than the rat cell line (Fig. 2E). The observed DBH down-regulation is likely to be a minimal baseline as it remains possible that vesicles may stick to the transwell membrane and EV passage restricted. Parasite restriction to the upper reservoir of the 0.4 µm filter transwells was confirmed by inoculating standard T. gondii cultures with media removed from upper and lower reservoirs and monitoring propagation (data not shown). Transwells were set up with infected fibroblasts in the top reservoir and uninfected noradrenergic cells in the bottom reservoir to assess whether the down-regulation was cell-type specific (Fig. 2F). The noradrenergic cells exposed to T. gondii-infected fibroblast cultures were unchanged in DBH gene expression.
A further indication that EVs were the permeable effector responsible for the DBH down-regulation was finding that the insoluble components, separated from soluble factors by ultracentrifugation, contained the DBH down-regulating activity in preliminary tests (data not shown). This provided the rationale for EV isolation and testing.
Epigenetic Changes with DBH Down-Regulation during Infection
As our findings indicated that TGS was responsible for DBH expression changes, the epigenetic state of the DBH gene was investigated (37). Methylation Sensitive Restriction Enzyme qPCR (MSRE-qPCR) was used to monitor DNA methylation levels in the DBH gene’s upstream region where the majority of CpGs are clustered (Fig. 3A). Methylation in the DBH upstream region rose from 16±4.6% to 66±3.8% in infected cultures of noradrenergic cells during the course of the infection (Fig. 3B; p = 0.00072). As the noradrenergic cells are sensitive to pH changes (ie. neurotransmitter synthesis and synaptic transmission affected) (38, 39), alkaline-shocked tachyzoites were used for in vitro infections, as in prior studies (40). This procedure elevated expression of bradyzoite markers BAG1 and SAG4 (Fig. S5). DBH methylation was also increased in infected human noradrenergic cultures (2.8-fold; range 2.3-4.8-fold; p = 0.0011) with a time-dependent increase in methylation in the region profiled (Fig. 3C).
As an indicator of EV involvement in the epigenetic changes, cells were treated with an inhibitor of EV biogenesis in the transwell system that restricts parasites and parasite-infected cells from uninfected cell culture. GW4869 inhibits sphingomyelinase which is required for vesicle budding in endosomal formation. Addition of GW4869 to the transwell system abrogated the DBH hypermethylation observed in uninfected noradrenergic cells exposed to the infected culture (Fig. 3D).
In order to examine the epigenetic effects of T. gondii on DBH in neurons, ex vivo experiments were performed with infection of organotypic brain slices of the prefrontal cortex, nucleus accumbens and ventral tegmental areas and the TGS in the neuron population measured, as the percentage of infected neurons during chronic infection is 0.002–0.14% (29, 41). Methylation of the DBH upstream region rose from 29 ± 2.7% in the brain tissue slices to 74 ± 4.6% in the infected slices (p = 0.000051) (Fig. 4A). DNA methylation levels were then analysed in the brains of chronically-infected mice. Neurons were purified by FACS and DBH methylation measured. The DBH gene in infected animals was 53 ± 7.7% methylated compared to 6.3 ± 2.0% in uninfected mice (p = 0.000045, Fig. 4B). For comparison, levels of total genomic DNA methylation was measured (Fig. 4C). No change in global DNA methylation was observed, as has previously been found with T. gondii infection (42). Hypermethylation of CpG residues upstream of the DBH gene were also found in NGS genomic bisulfite sequencing of infected cultures (Fig. S6).
The mechanism responsible for the parasite-induced DNA methylation and chromatin remodelling during DBH TGS was investigated. Experiments examined the role of DNA methyltransferase (DNMT) and histone deacetylation in this process. Infected noradrenergic cells were treated with the DNMT inhibitors RG108 and 5-azacytidine (5-AC). Both compounds disrupted the parasite-induced DBH down-regulation and DNA hypermethylation in rat noradrenergic cells, relative to marker (Figs. 4D, 4F). The DNMT inhibitors similarly blocked the TGS and epigenetic changes in human neuronal cells (Figs. 4E, 4G). This was not due to parasite sensitivity to the inhibitors since the inhibitors were non-toxic to parasites at concentrations tested (data not shown) and T. gondii lacks 5-methylcytosine (43). This implies that the DBH silencing involves DNMT. In contrast, the histone deacetylase inhibitor trichostatin A did not abrogate the DBH down-regulation or hypermethylation in infected cultures (Fig. 4D-G). Hence the mechanism may not involve histone deacetylation or activation was not captured within the experimental timeframe of 5 days and histone deacetylase is active at a different time in the epigenetic modification pathway.
To directly compare the epigenetic changes in infected versus uninfected (but infection-exposed) cells in the same culture, populations were enriched for cells containing GFP-expressing T. gondii (GFP+) and GFP- cells. DBH hypermethylation was found in both populations of cells (Fig. 5A, 5B; ANOVA test, p = 0.0089). Indeed, the GFP- cells exhibited DNA methylation levels at least equal to the GFP + cells. Hence, we conclude that direct infection was not required for hypermethylation of the DBH gene in a cell.
In order to identify the region of the DBH promoter involved in TGS in this model, the methylation status of the DBH upstream region in noradrenergic cells was measured by MSRE qPCR. Methylation of the DBH 5’ region was 34 ± 4.8% greater (p = 0.0013) in rat cells exposed to infected cells compared to exposure to uninfected cultures (Fig. 5C). Human noradrenergic cells exposed to infected cells also had increased methylation (49 ± 1.9%, p = 0.0003) of the DBH upstream region (Fig. 5D). Hence, both transcriptional down-regulation and promoter hypermethylation are induced in human and rat noradrenergic cells exposed to T. gondii cell cultures.
Based on the above findings, the isolated TINEVs ability to induce epigenetic changes was next investigated. DBH gene promoter methylation, examined using the MSRE-qPCR assay, increased from 16 ± 3.9% to 81 ± 7.9% over the course of the experiment. In a time course of exposure to TINEVs, methylation was not significantly changed after 24 or 48 hours of exposure but hypermethylation was observed at 96 hours of incubation (Fig. 5E; ANOVA p = 0.0008). Hence DBH transcriptional down-regulation and chromatin remodelling are induced by EVs from infected cells.
Toxoplasma-induced Neuronal EVs Contain an Antisense DBH lncRNA
We investigated the molecular components involved in the TINEV-induced TGS and epigenetic modifications. Preliminarily, the sensitivity of the TINEVs to ultraviolet (UV) irradiation was tested to provide an indication of whether an RNA component could be involved (44). TINEVs treated with UV radiation no longer induced hypermethylation of the DBH promoter (Fig. S7). It is possible that UV radiation inactivated proteins.
With the specificity of the differential gene expression and the role of non-coding RNAs (ncRNA) in regulating gene expression, we explored the potential presence of a long non-coding RNA (lncRNA). LncRNA has roles in neuronal gene expression such as highly specific antisense transcripts (45–48). EVs have been found to contain miRNAs and lncRNAs, although their functional significance is still unclear. With the specificity of the DBH down-regulation in T. gondii infection and its magnitude with TINEV treatment, the identification of a long non-coding RNA (lncRNA) was investigated (45, 49). A panel of primers walking upstream stepwise from the DBH coding region were used to screen for an antisense RNA (37). These identified a lncRNA in infected cultures (Fig. S8). RNA purified from TINEV preparations were found to contain the DBH antisense lncRNA (Fig. 6A-C, Fig. S9). Although the lncRNA location as internal EV cargo is unconfirmed, it is unclear of the functional significance of the location (50). The lncRNA is complementary to the DBH upstream region containing cis-regulatory elements and crosses the transcription start site. As a positive control for samples, miR-21 miRNA served to confirm RNA quality purified from the TINEVs (51). Intriguingly, the timeframe observed for the DBH hypermethylation in this study (Fig. 5E) is similar to dynamic DNA methylation changes observed in cells treated with promoter targeted antisense RNA (52). The observations above represent the first study to show a specific and functionally-relevant lncRNA to be transmitted from one neuronal cell to another regulating neurotransmission (53).
The DBH antisense lncRNA was found to contain the conserved U1 snRNP recognition site that is a consensus sequence identified in chromatin-bound lncRNAs (54). Further, predicted secondary structure analysis places the U1 snRNP recognition site in a loop region in minimum free energy secondary structures of DBH antisense lncRNA (Fig. 6D). In the published study, the recognition site, importantly to be in a loop region, allowed binding of U1 snRNP to the lncRNA to promote chromatin interaction. Future experiments will further investigate this novel induced lncRNA.