Microglial trogocytosis and the complement system regulate axonal pruning in vivo
Abstract
Partial phagocytosis—called trogocytosis—of axons by microglia has been documented in ex vivo preparations but has not been directly observed in vivo. The mechanisms that modulate microglial trogocytosis of axons and its function in neural circuit development remain poorly understood. Here, we directly observe axon trogocytosis by microglia in vivo in the developing Xenopus laevis retinotectal circuit. We show that microglia regulate pruning of retinal ganglion cell axons and are important for proper behavioral response to dark and bright looming stimuli. Using bioinformatics, we identify amphibian regulator of complement activation 3, a homolog of human CD46, as a neuronally expressed synapse-associated complement inhibitory molecule that inhibits trogocytosis and axonal pruning. Using a membrane-bound complement C3 fusion protein, we demonstrate that enhancing complement activity enhances axonal pruning. Our results support the model that microglia remodel axons via trogocytosis and that neurons can control this process through expression of complement inhibitory proteins.
Introduction
Microglia, the immune cells of the CNS, are vital for the maintenance and development of a healthy brain. Constantly surveilling the brain (Nimmerjahn et al., 2005; Wake et al., 2009), these highly phagocytic cells are thought to contribute to developmental synaptic remodeling by phagocytosing inappropriate or supernumerary synapses, a hypothesis that has derived considerable support from histological and immunohistochemical evidence identifying synaptic components within endocytic compartments in microglia (Paolicelli et al., 2011; Schafer et al., 2012; Stevens et al., 2007; Tremblay et al., 2010). This hypothesis is further supported by numerous studies demonstrating that microglia depletion leads to exuberant axonal outgrowth (Pont-Lezica et al., 2014; Squarzoni et al., 2014), impaired pruning of excess synapses (Ji et al., 2013; Milinkeviciute et al., 2019), and increased spine density (Wallace et al., 2020) during development.
Studies depleting microglia or disrupting microglial function provide indirect evidence to support the hypothesis that microglia remodel synapses through phagocytic mechanisms. An inherent weakness of indirect approaches is that the source of the synaptic material within microglia is unknown. For example, it is possible that synaptic components may be found within microglia due to clearance of apoptotic neurons rather than synaptic pruning. More direct approaches are required to verify whether microglia collect synaptic material from living neurons. Currently, direct evidence of complete elimination of synapses by microglial engulfment remains elusive. However, instead of removing entire synapses, microglia have been documented engaging in trogocytosis, or partial elimination, of axons and presynaptic boutons in ex vivo organotypic cultures and in fixed brain tissue using electron microscopy (Weinhard et al., 2018), although it remains to be seen whether microglial trogocytosis of axons is a phenomenon that occurs in vivo.
Even if we accept the hypothesis that microglia trogocytose the axonal compartment, many questions remain. What impact does partial elimination of presynaptic structures have on circuit remodeling? It is unclear if this phenomenon is required for proper wiring of neuronal circuits. Does microglial trogocytosis of axons affect axon morphology? While disrupting microglial function enhances axon tract outgrowth (Pont-Lezica et al., 2014; Squarzoni et al., 2014), it is unknown if this result is because of a disruption in microglial trogocytosis, or whether non-phagocytic mechanisms are in play. Is axonal trogocytosis by microglia mechanistically similar to complement-mediated synaptic pruning? There is extensive evidence demonstrating that the complement system regulates synaptic pruning by microglia via the complement protein C3 (Paolicelli et al., 2011; Schafer et al., 2012; Stevens et al., 2007). However, knockout (KO) mice lacking complement receptor type 3 (CR3), the receptor for activated C3, do not exhibit a deficit in microglial trogocytosis (Weinhard et al., 2018), raising the possibility that microglial trogocytosis of axons is mechanistically distinct from complement-mediated synaptic pruning.
In this study, we addressed these questions and directly observed in vivo trogocytosis of retinal ganglion cell (RGC) axons by individual microglial cells in real-time using the developing Xenopus laevis retinotectal circuit. We then developed an assay to monitor microglial trogocytosis of axons among the population of microglia in the optic tectum. To investigate the functional role of microglial trogocytosis, we depleted microglial cells and found that this enhanced axon arborization and reversed the behavioral responses to dark and bright looming stimuli. We identified amphibian regulator of complement activation 3 (aRCA3) (Oshiumi et al., 2009), a homolog of mammalian CD46, as an endogenously expressed, synapse-associated, complement inhibitory molecule in Xenopus laevis RGC neurons. Overexpression of aRCA3 inhibited trogocytosis and axonal pruning. Conversely, expression of a membrane-bound complement C3 fusion protein in RGCs enhanced axonal pruning. Our findings provide direct in vivo evidence supporting the hypothesis that microglia trogocytose presynaptic axonal elements and supports a model in which microglial trogocytosis regulates axonal pruning to promote proper neural wiring during development (Schafer et al., 2012; Weinhard et al., 2018). In this model, neurons exert local control of microglial trogocytosis and axonal pruning by expressing complement regulatory proteins.
Results
Microglia in Xenopus laevis tadpoles resemble neonatal mammalian microglia
Microglia in albino Xenopus laevis tadpoles were labeled with IB4-isolectin conjugated fluorophores for in vivo imaging. IB4-isolectin binds the RET receptor tyrosine kinase on microglial cells (Boscia et al., 2013), but does not lead to production of tumor necrosis factor-α or alterations in microglial morphology (Grinberg et al., 2011). Injection of fluorophore conjugated IB4-isolectin into the third ventricle of Xenopus laevis tadpoles labels highly mobile cells (Figure 1A and Video 1A) that have both ameboid-like and primitive ramified-like morphologies, resembling the morphology of embryonic microglial cells. Ameboid microglia are round or irregular cells with filopodia and/or pseudopodia, while primitive ramified microglia have scantly developed, poorly branching processes (Dalmau et al., 1997; Dalmau et al., 1998). Interestingly, IB4-isolectin-labeled cells can be observed switching back and forth between ameboid-like and primitive ramified-like morphologies (Figure 1B and Video 1B), suggesting that microglial morphology during development may be more dynamic than previously thought. Microglial cell bodies in the adult brain demonstrate low mobility (migration), while microglia in developing brain have highly mobile cell bodies (Smolders et al., 2019). Microglial mobility was determined to be 1.7 ± 0.3 µm/min (mean ± SD) (Figure 1C), comparable to what has been observed in neonatal rodent slice culture (1.4 μm/min at E15.5) (Smolders et al., 2017). To provide functional confirmation that IB4-isolectin-labeled cells are microglial cells, a laser irradiation injury was performed on the neuropil region. As expected of microglial cells, IB4-isolectin-labeled cells responded to injury by surrounding the damaged zone and removing injured tissue (Figure 1D and Video 2).
In developing zebrafish larvae, microglia primarily localize to the cell body layer of the optic tectum and are largely excluded from the tectal neuropil (Svahn et al., 2013). Conversely, microglia in developing mammalian models are found in neuropil regions (Dalmau et al., 1997; Hoshiko et al., 2012; Tremblay et al., 2010). In the Xenopus laevis retinotectal circuit, RGC axons project to the contralateral optic tectum, where they arborize and synapse on tectal neurons in the neuropil region (Figure 2A). To examine whether Xenopus microglia interact with the tectal neuropil, RGC axons innervating the neuropil were labeled by bulk electroporation with a plasmid encoding pH-stable green fluorescent protein (pHtdGFP). In vivo live imaging revealed that, microglia in developing Xenopus laevis can be found in both the cell body layer and the neuropil region (Video 3A). Microglia move in and out of the neuropil region from the cell body layer (Figure 2B and Video 3B) and move freely through the neuropil (Figure 2C and Video 3B). Additionally, microglia extend processes into the neuropil to contact axons, with interactions ranging from minutes to hours in duration (Figure 2D and Video 3C).
In vivo imaging of RGC neurons reveals microglial trogocytosis of axons and presynaptic structures
We then sought to examine whether microglia cells engage in trogocytosis of RGC axons in the developing Xenopus laevis retinotectal circuit. As endosomal organelles are typically acidified (Casey et al., 2010), pH-stability of dyes and fluorescent proteins is an important consideration when performing live imaging of trogocytosis (Shinoda et al., 2018). To reduce quenching of fluorescence, we used pHtdGFP (pKa = 4.8) (Roberts et al., 2016) which is more pH-stable than EGFP (pKa = 6.15). We expressed pHtdGFP in RGC axons by electroporation in the eye and labeled microglia with Alexa 594-conjugated IB4-isolectin. Two-photon live imaging revealed that the amount of green fluorescence associated with individual microglial cells increased following interactions with pHtdGFP-labeled axons. In the example shown in Video 4A, a microglial cell increases in green fluorescence by twofold following interaction with a pHtdGFP-labeled axon (Figure 2—figure supplement 1A). Similarly, in the example shown in Figure 2E and Video 4B, the green fluorescence in the microglial cell increased threefold following interaction with a pHtdGFP-labeled axon (Figure 2—figure supplement 1B). The real-time increase in microglial green fluorescence suggests direct transfer of fluorescent protein from the pHtdGFP-labeled axon and provides direct in vivo evidence of microglial trogocytosis of the presynaptic RGC axon.
Because microglia in developing Xenopus laevis tadpoles are highly mobile, the possibility that they may leave the imaging field complicates measuring trogocytosis of axons by individual microglia. To quantify microglial trogocytosis, we instead took the approach of measuring microglial green fluorescence across the population of microglia within the optic tectum. A greater number of pHtdGFP-labeled axons in the optic tectum is expected to lead to more frequent opportunities for trogocytotic interactions between microglia and pHtdGFP-labeled axons, resulting in greater amounts of green fluorescent material found within the microglial population. Based on this principle, we developed an assay to measure microglial trogocytosis of RGC axons in Xenopus laevis tadpoles. At developmental stage 39/40, RGC axons were labeled by retinal electroporation with plasmid encoding pHtdGFP, and microglia were labeled by intraventricular injection of Alexa 594-conjugated IB4-isolectin (Figure 3A). By 2 days post-labeling, axons begin expressing pHtdGFP and innervate the optic tectum. This also corresponds to the period when microglia begin to extensively colonize the optic tectum. At day 4 and day 5 post-labeling, the optic tectum was imaged by two-photon microscopy. The number of pHtdGFP-labeled axons present in the optic tectum was counted and the green fluorescence within the population of microglia in the imaging field was quantified using 3D masking with the IB4-isolectin channel (Figure 3B). To control for the possibility of RGC apoptosis, data were excluded if axonal blebbing was observed, or if the number of axons decreased from day 4 to day 5. Even when electroporation yields no labeled axons, some basal green fluorescence in microglial cells is still observed (Figure 3C and Figure 3—figure supplement 1). This is because microglia have high levels of autofluorescent molecules such as lipofuscin, bilirubin, and porphyrins (Mitchell et al., 2010). At day 4, there is a weak positive correlation between microglial green fluorescence and the number of pHtdGFP-labeled axons in the optic tectum, a relationship which is significantly strengthened on day 5 (Figure 3D). The change in green fluorescence associated with microglia from day 4 to day 5 significantly correlates with the number of pHtdGFP-labeled axons in the optic tectum (Figure 3—figure supplement 2A) suggesting that microglial cells accumulate pHtdGFP from intact axons between day 4 and day 5 post-labeling.
To determine if synaptic material is being trogocytosed by microglia, we generated a synaptophysin-pHtdGFP fusion protein (SYP-pHtdGFP). SYP is a presynaptic vesicle protein (Valtorta et al., 2004), and SYP fusion proteins are commonly used as synaptic vesicle markers (Nakata et al., 1998; Ruthazer et al., 2006). Expressing SYP-pHtdGFP in RGC neurons yielded axons with pHtdGFP puncta along the length of their terminal arbors (Figure 3E). When SYP-pHtdGFP is expressed in RGC axons, on day 4, no correlation is observed between microglial green fluorescence and the number of SYP-pHtdGFP-labeled axons in the optic tectum. However, by day 5, a correlation between microglial green fluorescence and the number of SYP-pHtdGFP-labeled axons is observed (Figure 3F). The change in microglial green fluorescence from day 4 to day 5 is proportionate to the number of SYP-pHtdGFP-labeled axons in the optic tectum (Figure 3—figure supplement 2B), suggesting that microglial cells trogocytose and accumulate presynaptic elements from RGC axons over the period from day 4 to day 5.
Depletion of microglial cells enhances RGC axon branching and reverses the profile of behavioral responses to dark and bright looming stimuli
To assess the functional roles of microglial trogocytosis, we first depleted microglia using PLX5622. PLX5622 is an inhibitor of colony stimulating factor one receptor (CSF1R), which is a tyrosine kinase receptor essential for microglia survival (Elmore et al., 2014; Erblich et al., 2011). Animals reared in 10 μM PLX5622 had significantly reduced microglia numbers in the optic tectum compared to vehicle-treated animals (Figure 4A and Figure 4B). Morphological analysis of surviving microglia also revealed a reduction in the number of processes per microglial cell, indicating a less ramified morphology (Figure 4C).
Next, we interrogated the effect of microglial depletion on the morphology of single RGC axons. Axons were followed for several days in control and microglia-depleted animals (Figure 4D). Microglial depletion with PLX5622 did not affect axon arbor length (Figure 4E) but significantly increased the number of axon branches (Figure 4F), suggesting that microglia negatively regulate axonal arborization during development.
We then sought to delineate the functional effects of microglial depletion on the development of the retinotectal circuit. Previous reports in Xenopus laevis and in zebrafish have shown that the retinotectal circuit is a vital processing and decision-making center for the visual detection of looming objects (Dong et al., 2009; Khakhalin et al., 2014). Therefore, we developed a free-swimming looming stimulus assay in Xenopus tadpoles (Figure 5A) to determine the functional outcomes of disrupting microglial function in this circuit. Tadpole behavioral responses to dark looming stimuli or bright looming stimuli (Figure 5B) were recorded and custom computer vision software was used to track the locomotor response of the tadpole (Figure 5C).
Figure 5D shows contrails from representative animals in response to dark looming stimuli. In vehicle-treated animals, dark looming stimuli evoked stereotypical defensive escape behavior (Video 5A), whereas microglia-depleted animals were less likely to make defensive responses to dark looming stimuli (Video 5B). Presentation of dark looming stimuli elicited an increase in velocity in control animals (Figure 5E), with a peak in instantaneous velocity at 1.1 ± 0.1 s (mean ± SD, n = 8) post-stimulus. The distance traveled by control tadpoles over the 3 s immediately following presentation of a dark looming stimulus increased compared to the period before the stimulus was presented (Figure 5F). This increase was absent in microglia-depleted animals. We categorized tadpole responses to looming stimuli as exhibiting defensive behavior, absence of defensive behavior or undeterminable (excluded from response rate calculation). Microglia-depletion significantly reduced the response rate to dark looming stimuli (Figure 5G).
Figure 5H shows contrails from representative animals in response to bright looming stimuli. Surprisingly, bright looming stimuli rarely evoked defensive escape behavior in control animals (Video 6A) but often evoked robust responses in microglia-depleted animals (Video 6B). Presentation of bright looming stimuli elicited an increase in velocity in microglia-depleted animals (Figure 5I), with a peak in instantaneous velocity at 1.9 ± 0.6 s (mean ± SD, n = 8) after the onset of stimulus presentation. The distance traveled by microglia-depleted tadpoles over the 3 s immediately following presentation of a bright looming stimulus increased compared to the period before the stimulus was presented (Figure 5J). This increase was absent in vehicle-treated animals. Furthermore, microglia-depletion significantly increased the response rate to bright looming stimuli (Figure 5K).
In trials where animals exhibited escape behavior to dark or bright looming stimuli, microglia depletion did not significantly alter escape distance, maximum escape velocity, or escape angle, suggesting that microglial depletion did not disrupt motor function (Figure 5—figure supplement 1).
Bioinformatic identification of amphibian regulator of complement activation 3 (aRCA3), a neuronally expressed membrane-bound complement inhibitory molecule and homolog of human CD46
Microglial depletion enhanced axonal arborization and perturbed functional development of the retinotectal circuit. As microglia are known to secrete trophic factors and cytokines (Parkhurst et al., 2013; Solek et al., 2018), signaling mechanisms unrelated to trogocytosis may be at work. To better discern a role of trogocytosis in modulating axon morphology, we searched for a molecule that would inhibit axon trogocytosis. Because complement C3 is a known mediator of phagocytosis (Brown and Neher, 2014) and has been proposed to accumulate at synapses and tag them for removal by microglial phagocytosis (Schafer et al., 2012; Stevens et al., 2007), we carried out a bioinformatics screen to identify neuronally expressed proteins that have potential inhibitory activity against complement C3.
To search for proteins that inhibit complement C3, we first used the STRING Protein-Protein Association Network to identify proteins which interact with human complement C3 (Figure 6A). To select for complement inhibitory proteins, our search was then refined by focusing on proteins containing CCP (domains abundant in complement control proteins) motifs, which are highly conserved modules abundant in proteins responsible for negative regulation of the complement system (Norman et al., 1991). The top CCP-containing proteins identified by STRING were queried on the Allen Brain Institute human multiple cortical areas RNA-seq dataset (Allen Institute for Brain Science, 2015; Hodge et al., 2019). Of the CCP-containing proteins identified by STRING, only CD46 is highly expressed by human CNS neurons (Figure 6B).
Using NCBI protein-protein BLAST, the Xenopus laevis protein with the highest similarity to human CD46 was identified as amphibian Regulator of Complement Activation 3 (aRCA3) (Table 1). aRCA3 is a protein uncharacterized in Xenopus laevis but characterized as a membrane-associated complement regulatory protein in Xenopus tropicalis (Oshiumi et al., 2009). Like human CD46, Xenopus laevis aRCA3 is a type I transmembrane protein with many extracellular CCP domains and a stop-transfer anchor sequence near the C-terminus of the protein (Figure 6C). As a comparison, the second and third most similar Xenopus laevis proteins, complement receptor two and complement component four binding protein, do not share a protein architecture with CD46. The tertiary structure of Xenopus laevis aRCA3 resembles that of human CD46 (Video 7A and Video 7B). Despite having twice the number of CCP domains as human CD46, aRCA3 is a likely homolog of CD46. The number of CCP domains of CD46 homologs varies across species—for example, human (Lublin et al., 1988) and rodent CD46 have four CCP domains (Miwa et al., 1998), the Gallus gallus homolog (CREM) has six CCP or CCP-like domains (Oshiumi et al., 2005), and homologs in Danio rerio (rca2.1) (Tsujikura et al., 2015), Equus caballus (MCP), and Coturnix japonica (MCP-like) each have five CCP domains.
To investigate whether RGC neurons endogenously express aRCA3, retinal sections were probed by fluorescence RNAscope in situ hybridization (Wang et al., 2012). Probes against aRCA3, the positive control transcript RNA polymerase II subunit A (polr2a.L), and the negative control transcript bacterial dihydrodipicolinate reductase (DapB) were hybridized on PFA fixed retinal sections (Figure 6D). As expected, DapB transcripts were not detected in the retina, whereas the housekeeping polr2a transcript was ubiquitously expressed. In contrast, aRCA3 transcripts localize to the ganglion cell layer (Figure 6E), revealing that the aRCA3 gene is endogenously expressed in RGC neurons.
aRCA3 associates with synapses, enhances axonal arborization and inhibits microglial trogocytosis
It is unknown whether aRCA3 localizes to the correct cellular compartments to protect axons from complement attack. To investigate this, we examined the localization of aRCA3 by tagging it with mCherry (Figure 7A). This aRCA3-mCherry fusion protein was expressed together with SYP-pHtdGFP, which localizes to synaptic puncta (Figure 3E). SYP-pHtdGFP and aRCA3-mCherry colocalize when co-expressed in the same axon (Figure 7B). The high degree of colocalization between SYP-pHtdGFP and aRCA3-mCherry suggests that aRCA3 distributes to similar subcellular compartments as the presynaptic marker SYP (Figure 7C).
Next, we examined the effect of aRCA3 overexpression on microglial trogocytosis. pHtdGFP was expressed in RGC neurons with or without overexpression of aRCA3, and microglial trogocytosis was quantified by measuring microglial green fluorescence on day 4 and day 5 post-labeling. Expression of pHtdGFP in RGC neurons resulted in a significant correlation between microglial green fluorescence and the number of pHtdGFP-labeled axons on day 5 post-labeling (Figure 7D), a relationship that was not observed when aRCA3 was overexpressed (Figure 7E). Overexpression of aRCA3 significantly reduced the correlation coefficient between the change in microglial green fluorescence from day 4 to day 5 and the number of pHtdGFP-labeled axons (Figure 7—figure supplement 1), suggesting that overexpression of aRCA3 reduces microglial trogocytosis.
We then examined the effect of aRCA3 overexpression on axon morphology in single RGC axons over several days (Figure 7F). Overexpression of aRCA3 resulted in significantly greater axon arbor length compared to control axons (Figure 7G). Similarly, axon branch number was greater in axons overexpressing aRCA3 compared to control axons (Figure 7H). These results show that aRCA3 overexpression causes axons to have larger and more highly branched arbors.
While CD46 is best known for its ability to inactivate complement C3 and complement C4, CD46 can also signal through intracellular tyrosine kinase activity under certain conditions (Riley-Vargas et al., 2004). This raises the potential caveat that the increased axonal arborization induced by overexpression of aRCA3 may result from aberrant intracellular signaling. Using a bioinformatic tool, NetPhos3.1 with cutoff scores > 0.6 (Blom et al., 2004), we did not find a predicted tyrosine kinase phosphorylation site on the cytoplasmic region of aRCA3. Nonetheless, it is possible that aRCA3 may exert some of its effects through an intracellular signaling pathway. To provide further functional validation that the complement system affects axonal arborization, we examined the effects of enhancing complement activity on axon morphology.
Expression of a membrane-bound complement C3 fusion protein reduces RGC axon size and branching
If aRCA3 affects axon morphology by inhibiting complement activity, enhancing complement activity on single axons should produce effects opposite to that of aRCA3 overexpression. To explore this possibility, we designed an axon surface-localized complement C3 fusion protein to enhance complement activity on individual axons (Figure 8A). Synaptobrevin, also known as vesicle-associated membrane protein 2 (VAMP2), is concentrated in synaptic vesicles, though a significant fraction of VAMP2 is also present on the axon surface (Ahmari et al., 2000; Sankaranarayanan and Ryan, 2000). We cloned Xenopus laevis complement C3 and fused the N-terminus to the extracellular C-terminus of Xenopus laevis VAMP2. This design was chosen in favor of a GPI anchor design that modifies the C-terminus C345C domain of complement C3 as this domain undergoes major rearrangement during activation and proteolysis (32° rotation, 10 Å translation) (Janssen et al., 2005). In contrast, the N-terminus of complement C3 is exposed on the surface of the protein and located on the MG1 domain, a domain that does not undergo marked confirmational changes upon complement C3 activation and proteolysis (3° rotation, 1 Å translation). Thus, expression of VAMP2-C3 in RGC neurons results in axons tagged with extracellular membrane-bound complement C3. While we used complement C3 precursor to generate the VAMP2-C3 fusion protein, complement C3 undergoes spontaneous, low-level activation through the alternative complement pathway (Pangburn et al., 1981).
VAMP2-C3 was co-expressed with pHtdGFP in RGC neurons, and axons were monitored over several days (Figure 8B). To control for the possibility that VAMP2 overexpression may affect axon morphology, we also overexpressed VAMP2 alone in RGC axons. Expression of VAMP2-C3 in RGC axons significantly reduced axon arbor length and axon branch number when compared to control or VAMP2 overexpression (Figure 8C and Figure 8D), demonstrating that enhancing complement activity negatively regulates axonal arborization.
Discussion
In vivo evidence of microglial trogocytosis of axons during healthy development
Previous real-time imaging experiments in ex vivo slice culture have shown that microglia trogocytose presynaptic elements (Weinhard et al., 2018). Our real-time imaging results now add in vivo support to the hypothesis that microglia engulf synaptic material via trogocytosis. Additionally, labeling of axons with pHtdGFP or SYP-pHtdGFP resulted in an increase in fluorescent label within microglia between days 4 and 5 post-labeling. As SYP-pHtdGFP is primarily localized to synaptic vesicles, this suggests that the axonal material that microglia engulf contains presynaptic vesicles and is consistent with correlative electron microscopy studies that have demonstrated putative presynaptic vesicles within microglia (Weinhard et al., 2018), although we did not use correlative electron microscopy to confirm this in our study.
Numerous past studies have shown that during development microglia have synaptic components within phagocytic compartments (Paolicelli et al., 2011; Schafer et al., 2012; Stevens et al., 2007; Tremblay et al., 2010). However, apoptosis is a prominent feature of neural development (Nijhawan et al., 2000), and it is unclear whether the source of synaptic material within microglia is because of engulfment of apoptotic neuronal components or if it is due to engulfment of synaptic material from live neurons. When investigating trogocytosis, we minimized the effect of apoptosis by excluding measurements when the number of axons between imaging sessions decreased or when axonal blebbing was observed, providing in vivo direct evidence that microglia accumulate fluorescent label from living axons.
Expression of complement inhibitory or complement enhancing molecules in neurons regulates axon morphology and trogocytosis by microglia
Activation of complement C3 exposes a reactive thioester bond that covalently attaches to amine or carbohydrate groups on cell surfaces (Sahu et al., 1994). Microglia express CR3, a receptor for complement C3 (Ling et al., 1990), and binding of CR3 to its ligand induces phagocytosis (Newman et al., 1985). Complement C3 localizes to synapses during development and tags synapses for removal by microglia (Stevens et al., 2007). The importance of complement C3 in synaptic pruning has been studied by disrupting the C3 pathway using C3 KO mice (Stevens et al., 2007), CR3 KO mice (Schafer et al., 2012), or recently, with the exogenous neuronal expression of the complement inhibitory proteins Crry (Werneburg et al., 2020) and CD55 (Wang et al., 2020). Here we show that enhancing the C3 pathway with a membrane-bound VAMP2-C3 fusion protein increases axonal pruning at the single axon level. Conversely, overexpression of an endogenous membrane-bound complement inhibitory molecule inhibits axonal pruning and trogocytosis (Figure 9A).
While VAMP2-C3 clearly produced effects on axon morphology, we did not rule out the possibility that the expression of VAMP2-C3 in RGC neurons may be producing effects on axons that are not mediated by the complement system. Additionally, we were unable to determine whether VAMP2-C3 expression enhanced microglial trogocytosis as electroporation of RGC neurons with VAMP2-C3 rarely yielded labeled axons. We speculate that expression of VAMP2-C3 in RGC neurons reduces viability, as previous studies have shown that the complement cascade mediates phagocytosis of RGC neurons during development (Anderson et al., 2019). Furthermore, our study indirectly examined the effects of VAMP2-C3 and aRCA3 on microglial trogocytosis. Future experiments could directly examine the effects of VAMP2-C3 expression and aRCA3 overexpression on trogocytosis rates by real-time imaging.
It is hypothesized that neurons endogenously express membrane-bound complement inhibitory molecules to protect synapses from phagocytosis (Stephan et al., 2012; Stevens et al., 2007), and our data support this hypothesis. Recently, sushi repeat protein X-linked 2 (SRPX2) has been identified as a endogenous neuronal complement inhibitor in the mammalian system which protects synapses from complement C1 tagging (Cong et al., 2020). SRPX2 is a secreted protein, and while secreted complement inhibitory molecules can act via cell-autonomous mechanisms to protect the neuron that produced it, such molecules also have non-autonomous actions due to the ability to diffuse through the extracellular environment. Conversely, a membrane-bound complement inhibitory molecule acts solely through cell-autonomous mechanisms and allows for greater spatiotemporal control of local complement protection. While the identity of a membrane-bound complement inhibitory molecule endogenous to mammalian neurons remains elusive, in this study we characterized Xenopus laevis aRCA3, an endogenously expressed, synaptic vesicle-associated, complement inhibitory molecule. aRCA3 is the most similar amphibian homolog of human CD46. CD46 is a membrane-bound complement inhibitory molecule that cleaves activated complement C3 and C4 (Barilla-LaBarca et al., 2002). Using the Allen Brain Institute human multiple cortical areas RNA-seq dataset (Allen Institute for Brain Science, 2015), we report that CD46 transcripts are enriched in human neurons. Interestingly, CD46 associates directly with β1-integrins (Lozahic et al., 2000), an adhesion molecule present on neuronal surfaces (Neugebauer and Reichardt, 1991) and enriched in synaptosomes (Chan et al., 2003), suggesting that it too may be localized to axon surfaces and synapses. In our study, we show that aRCA3 expression inhibits both axon trogocytosis by microglia and axonal pruning—we speculate that CD46 may perform similar functions in mammals. Excessive synaptic pruning is thought to be one of the underlying causes of schizophrenia (Sellgren et al., 2019), and three large-scale genetic susceptibility studies have identified the CD46 gene as a significant Schizophrenia-risk locus (Håvik et al., 2011; Kim et al., 2020; Ripke et al., 2014). Clearly, the role of CD46 in neurodevelopment warrants further study.
Microglia actively suppress exuberant arborization at the single axon level
Disrupting microglial function by depletion increases axon tract innervation in prenatal models (Pont-Lezica et al., 2014; Squarzoni et al., 2014). As depleting microglia increases the number of neural progenitor cells (Cunningham et al., 2013) and the number of RGC neurons in the embryonic retina (Anderson et al., 2019), it was unclear whether the increased axon tract innervation that occurs following microglial depletion resulted from a deficit of microglial-mediated axonal pruning, or whether it was because of an increase in the overall number of axons. To address this, we examined the effect of microglial depletion on the morphology of single axons, showing that microglia depletion significantly increased axon branch number. Thus, our study supports the model that the increase in axon tract innervation resulting from microglial depletion is because of increased arborization of individual axons (Figure 9B).
The observations that microglial depletion and overexpression of the complement inhibitor aRCA3 in RGC neurons enhanced axonal arborization suggests that both phenomena may act through a common mechanism. These results support the hypothesis that microglial trogocytosis suppresses both axon branching and axon growth through the removal of extra branches and axonal material. While this explanation is at odds with the observation that microglia-depletion did not increase axon arbor length, it is important to note that with aRCA3 overexpression, a single axon gains a relative growth advantage over other axons—this is distinct from the case of microglial depletion, where all axons profit from the same increased growth. If all axons are growing larger at the same pace, competitive mechanisms (Gosse et al., 2008; Ruthazer et al., 2003) might be expected to constrain axon arbor sizes. Another interesting possibility is that microglia preferentially engage in trogocytosis at short, dynamic branches over longer, established, mature branches. In this case, microglial depletion results in a greater number of short branches, which over the short time course in this study, do not significantly impact axon branch length.
It is worth considering that there may be unintended side-effects of microglial depletion with PLX5622 and that many microglial functions beyond trogocytosis of axons may be altered. In our study, we found that CSF1R inhibition reduced process branching of remaining microglial cells, an observation that has also been seen in rodents (Elmore et al., 2014). Decreased microglial process complexity often accompanies release of cytokines and other neuromodulatory molecules (Karperien et al., 2013; Hanisch, 2002). While we did not study the effect of microglial depletion on cytokine and chemokine production in this study, past studies have shown that cytokines and chemokines in non-disease mouse models are not altered by PLX5622 (Reshef et al., 2017). A further potential confound is that low expression of CSF1R has been reported in some neurons in mice hippocampus (Luo et al., 2013). It is possible that CSF1R is present in Xenopus laevis RGC neurons and we did not exclude the possibility that PLX5622 may be acting directly on RGC axons in our study.
Microglia contribute to proper wiring of the retinotectal circuit in developing Xenopus laevis
A feed-forward network drives visually evoked escape behavior in Xenopus tadpoles (Khakhalin et al., 2014). In the retina, photoreceptors act on bipolar cells, which act on RGCs. Some RGCs are tuned to detect looming (Dunn et al., 2016; Münch et al., 2009), and the population of RGCs that respond to dark looming stimuli are distinct from the population of RGCs that respond to bright looming stimuli (Temizer et al., 2015). In zebrafish, these distinct RGC populations project to different arborization fields in the optic tectum (Robles et al., 2014; Temizer et al., 2015) where the execution of looming computation and behavioral decision-making occurs (Barker and Baier, 2015; Fotowat and Gabbiani, 2011). In highly predated animals such as crabs, zebrafish, and mice (Oliva et al., 2007; Temizer et al., 2015; Yilmaz and Meister, 2013), dark looming stimuli—which signal imminent threat such as an oncoming object or predator—reliably induce escape responses, whereas bright looming stimuli—which may occur as the animal exits a tunnel or traverses its environment—are less effective.
Microglial depletion does not appear to alter motor functionality, suggesting that the hindbrain and motor circuitry remain intact. Instead, the effects of microglia depletion on looming-evoked escape behavior must occur further upstream in the retina, the optic tectum, or in the projections from the tectal neurons to the hindbrain. As microglia depletion induces exuberant RGC axonal arborization, one explanation for why microglia depletion disrupts defensive behavior to dark looming stimuli is that the dark looming-sensitive RGC axons cannot effectively wire with their tectal neuron counterparts, resulting in a misclassification of threatening visual stimuli when axonal pruning is disrupted. When microglia are depleted, axons form more errant connections, and dark looming-sensitive RGC axons may instead wire together with tectal neurons which are not associated with threat classification, reducing the probability that dark looming stimuli elicit escape responses. We predict that the reverse is true for bright looming-sensitive RGC axons, whereby microglia depletion results in increased likelihood of errant wiring between bright looming RGC axons and tectal neurons that classify threatening visual stimuli, leading to enhancement of escape behavior to bright looming stimuli. Conversely, it is possible that dark and bright looming-sensitive RGC axons are differentially affected by microglial depletion or the circuitry of the retina upstream of RGCs is remodeled following microglial depletion, both of which could produce the abnormal behavioral outcomes we observed. Indeed, the mechanism linking microglial depletion to improper response to looming stimuli is unclear and warrants future investigation.
Synaptic pruning and trogocytosis are two sides of the same coin: axonal pruning
Classically, synaptic pruning has been described as the engulfment and elimination of synapses, a phenomenon dependent on the complement pathway (Paolicelli et al., 2011; Perry and O'Connor, 2008; Schafer et al., 2012; Stevens et al., 2007). While no deficit in microglial trogocytosis has been observed in ex vivo cultures from CR3 KO mice (Weinhard et al., 2018), one interpretation for this result is that synaptic pruning and trogocytosis are unique phenomena mediated by different mechanisms. However, our study suggests that complement-mediated synaptic pruning and trogocytosis are mechanistically related. We find that the complement system influences both axonal trogocytosis and axonal arborization, and it is possible that in CR3 KO mice, compensation through one of the other receptors for activated C3 occurs. For example, microglia express complement receptor type 4 (CR4) (Allen Institute for Brain Science, 2015; Hodge et al., 2019) which can bind activated complement C3 to induce phagocytosis (Ross et al., 1992). Our data supports the hypothesis that neurons control microglial-mediated circuit remodeling through the expression of endogenous membrane-bound complement inhibitory molecules to regulate microglial trogocytosis.
Materials and methods
Lead contact and materials availability
Request a detailed protocolPlasmids generated in this study will be made available upon request. Further information and requests for resources and reagents should be directed to the Lead Contact, Edward Ruthazer (edward.ruthazer@mcgill.ca).
Xenopus laevis tadpoles
Request a detailed protocolAdult albino Xenopus laevis frogs (RRID:XEP_Xla200) were maintained and bred at 18°C. Female frogs were primed by injection of 50 IU pregnant mare serum gonadotropin (ProSpec-Tany TechnoGene Ltd., Ness-Ziona, Isreal). After 3 days, male and primed female frogs were injected with 150 IU and 400 IU of human chorionic gonadotropin (Sigma-Aldrich, Oakville, CA) into the dorsal lymph sac, respectively. The injected male and female frogs were placed in isolated tanks for mating. Embryos were collected the following day and maintained in Modified Barth’s Saline with HEPES (MBSH) in an incubator at 20°C with LED illumination set to a 12 hr/12 hr day-night cycle and staged according to Nieuwkoop and Faber (NF) developmental stages (Nieuwkoop and Faber, 1994). All experiments were conducted according to protocol application number 2015–7728 approved by The Animal Care Committee of the Montreal Neurological Institute and in accordance with Canadian Council on Animal Care guidelines. Xenopus laevis sex cannot be determined visually pre-metamorphosis, and thus the sex of experimental animals was unknown.
Labeling of microglia, axons, and brain structures
Request a detailed protocolStage 39/40 tadpoles were anesthetized with 0.02% MS-222 in 0.1X MBSH. For labeling of microglia, tadpoles received an intracerebroventricular (icv) injection to the third ventricle of 1 mg/ml Alexa 594 conjugated IB4-isolectin. A minimum of 48 hr was allowed to pass before live imaging studies were commenced to allow for binding and update of IB4-isolectin by microglial cells. For concurrent labeling of brain structures, 1 mM CellTracker Green BODIPY in 10% DMSO was injected icv after isolectin injection. For labeling of RGC axons, electroporation of plasmid DNA into RGC cells was performed as described previously (Ruthazer et al., 2006; Ruthazer et al., 2013a). In brief, stage 39/40 tadpoles were anesthetized with 0.02% MS-222 in 0.1X MBSH. A glass micropipette was back filled with endotoxin free maxi-prep plasmid solution (2–3 μg/µl) and fast green to visualize the injection. The micropipette was advanced into the eye, and DNA solution pressure injected. The micropipette was then withdrawn, and parallel platinum electrodes were placed to bracket the eye. For bulk labeling of RGCs, a Grass Instruments SD9 electrical stimulator was used to apply 4–6 pulses of 2.4 ms duration at 36 V with a 3 µF capacitor was placed in parallel to obtain exponential decay current pulses. For single-cell labeling of RGCs, 4–6 pulses of 1.6 ms duration at 36 V was applied.
Microglial depletion
Request a detailed protocolPLX5622 was dissolved in DMSO at 20 mM, aliquoted, and stored at −20°C. Thawed aliquots of PLX5622 were sonicated briefly for 1 min before dilution in polyethylene glycol (PEG) 400. PEG 400 solution was then further diluted in a solution of 0.1X MBSH containing non-ionic surfactants poloxamer 407 and D-α-tocopherol polyethylene glycol 1000 succinate (TPGS) to form a mixed micelle drug delivery vehicle (Guo et al., 2013). Final vehicle rearing solution consisted of 2.5% PEG 400, 0.04% poloxamer 407, 0.01% TPGS, and 0.05% DMSO in 0.1X MBSH. Initially, animals were reared in 0.1X MBSH. At stage 35–40, animals were transferred to vehicle rearing solution with or without 10 μM PLX5622 for rearing. Rearing solutions were refreshed daily with newly prepared drug and vehicle solutions.
Two-photon live imaging
Request a detailed protocolTwo-photon live imaging of axons and microglia cells was performed as described previously (Ruthazer et al., 2013b). Tadpoles were placed in a Sylgard 184 silicone imaging chamber with their dorsal side facing up and were covered with a #1 thickness glass coverslip. One μm interval z-stacks of the optic tectum were collected on a custom-built two-photon microscope (Olympus BX61WI with Olympus FV300 confocal scan head) outfitted with an Olympus 1.0 NA 60x water-immersion objective, or a Thorlabs multiphoton resonant scanner imaging system outfitted with an Olympus 1.0 NA 20x water-immersion objective. Excitation was produced using a Spectra-Physics InSight3X femtosecond pulsed laser. Images were collected at 512 × 512 pixels with Fluoview 5.0 or ThorImage LS.
For real-time imaging studies, stage 46–48 tadpoles were immobilized by brief (2–8 min) immersion in freshly thawed 2 mM pancuronium bromide in 0.1X MBSH. Animals were then maintained in 0.1X MBSH for imaging. Z-stacks were collected at 6 min intervals. For daily imaging studies, stage 43–47 tadpoles were anesthetized by immersion in 0.02% MS-222. Z-stacks were collected, and animals were returned to 0.1X MBSH rearing solution.
When imaging microglia (Alexa 594 conjugated IB4-isolectin) and RGC axons (eGFP or pHtdGFP) concurrently, excitation wavelength was set at 830 nm and a 565 nm emission dichroic was used in conjunction with green (500–550 nm) and red (584–676 nm) filters for fluorescence emission detection on separate photomultiplier tubes. When imaging microglia (Alexa 594 conjugated IB4-isolectin) and brain structures (CellTracker Green BODIPY), imaging was done sequentially. Alexa 594 was first imaged with excitation wavelength set to 810 nm and fluorescence emission detection through the red filter. Imaging of CellTracker Green BODIPY was then performed at excitation wavelength 710 nm and the fluorescence emission detected through the green filter. When imaging RGC axons alone (pHtdGFP), excitation wavelength was set to 910 nm and fluorescence emission detected through the green filter. When imaging aRCA3-mCherry-Myc and SYP-pHtdGFP, images were captured concurrently at wavelength 990 nm through the red and green filters, respectively.
Two-photon laser-induced injury
Request a detailed protocolLaser irradiation injury was carried out by repeatedly scanning an approximately 10 × 10 μm region under high laser power at wavelength 710 nm for 10 s.
Microscopy image processing and analysis
For visualization purposes, the 3D object counter (Bolte and Cordelières, 2006) and 3D ROI manager (Ollion et al., 2013) plugins were used on z-stacks for segmentation and effacement of melanophores on the dorsal dermal surface of the animal before maximum intensity projection. Analysis was performed by an experimenter blind to treatment group.
Axon morphology analysis
Request a detailed protocolTwo-photon z-stacks were captured daily from 2 to 5 days post-electroporation. Z-stacks were denoised with CANDLE (Coupé et al., 2012), and were manually traced using Imaris six software.
Microglia quantification
Request a detailed protocolThe number of microglia and microglia process numbers were counted manually within 150-μm-thick z-stacks of the optic tectum.
Microglia mobility
Request a detailed protocolFor microglial mobility experiments, microglia in the optic tectum were imaged over 2 hr at 6 min intervals and manual tracking of microglia was carried out using the TrackMate plugin (Tinevez et al., 2017).
Real-time imaging of trogocytosis of axons by individual microglial cells
Request a detailed protocolCollected z-stacks were denoised with the 3D hybrid median filter plugin. Time series were registered by descriptor based series registration (Preibisch et al., 2010). Microglia 3D region of interests (ROIs) were generated using the 3D object counter plugin (Bolte and Cordelières, 2006) using the same threshold across all time points. When microglia 3D ROIs were overlapping, they were split manually. Green fluorescence intensity within microglia 3D ROIs was measured using the 3D ROI manager plugin (Ollion et al., 2013) and background subtraction was carried out with a rotated or mirrored z-stack. Baseline values were defined as the average of the first 15 frames, except in the case of microglial interaction, wherein any frames after an interaction were excluded from baseline calculation. Movies were generated using Imaris software or ImageJ.
Trogocytosis assay
Request a detailed protocolMicroglia were labeled with Alexa 594-conjugated isolectin and axons with pHtdGFP or SYP-pHtdGFP as described above. On day 4 and day 5 post-labeling, 150-µM-thick z-stacks were collected from the optic tectum. Laser power and photomultiplier tube voltage was kept constant on both days. The number of labeled axons was counted manually. In studies when axons were labeled with pHtdGFP, data were excluded if axonal blebbing was detected, or if the number of labeled axons fell from day 4 to day 5. This data exclusion step was not performed when axons were labeled with SYP-pHtdGFP. Using the red microglia channel, microglia 3D ROIs were automatically generated using the 3D object counter plugin (Bolte and Cordelières, 2006). Microglia segmentation threshold was set to the mean pixel intensity of the red channel plus two standard deviations. A minimum threshold size of 1500 pixels was used to exclude background from analysis. The 3D ROI manager plugin (Ollion et al., 2013) was then used to measure the mean gray value of the green channel within microglial cell ROIs. To reduce background signal in the green channel, the mode of the green pixel intensity was calculated from the z-stack and subtracted. Non-microglial ROIs such as melanophores as well as microglia which had overlapping pixels with the axon were excluded manually and the average microglial green fluorescence intensity for the population of microglia in the z-stack was determined.
Colocalization analysis
Request a detailed protocolTwo-photon z-stacks were captured on day three post-electroporation. aRCA3-mCherry-Myc and SYP-pHtdGFP channels were denoised using the 3D hybrid median filter ImageJ plugin. A 3D mask of the axon was generated by summing the red and green channels and using an intensity threshold for segmentation of the axon. Pearson’s correlation coefficient (Manders et al., 1993) of the masked axon was calculated using Imaris software, with thresholds determined automatically using the method of Costes (Costes et al., 2004). The ScatterJ plugin (Zeitvogel et al., 2016) was used to generate scatter plots.
Looming stimulus behavioral assay
Experimental setup
Request a detailed protocolStage 47 tadpoles were placed within a closed 60 mm petri dish and allowed to swim freely. The petri dish was placed on the bottom of a large shallow (20 cm diameter, 8 cm deep) glass culture dish filled completely with 0.1X MBSH. The large shallow glass culture dish was placed on a purple 3D printed stage to allow for automated segmentation of albino tadpoles. A webcam was placed above the tadpole to record tadpole behavior, while a 2000 lumens projector customized by 3D printing to shorten the focal distance, was used to project visual stimuli onto a piece of 96 white 20 lb bond copy paper taped to the side of the large glass culture dish.
Stimulus presentation and recording
Request a detailed protocolDark or bright looming stimuli were presented using custom code written in Python 3. An expanding circle was projected onto the culture dish (800 × 600 pixels or 10.6 × 8 cm) using the PsychoPy library (Peirce et al., 2019). For dark looming stimuli, a black (29 cd/m2) expanding circle was shown over a white (208 cd/m2) background. For bright looming stimuli, a white expanding circle was shown over a black background. During looming, the diameter of the circle expanded exponentially at 10% per frame at 60 frames per second, from a diameter of 54 pixels (7 mm) until it encompassed the entire screen 0.5 s later. After another 0.8 s, the screen was reset and a 10 s refractory period commenced. In parallel to stimulus presentation, 480 p video was recorded by webcam using the OpenCV library (Bradski, 2000). Ten looming stimulus trials per animal were recorded.
Automated tracking of tadpole behavior
Request a detailed protocolCustom computer vision tadpole tracking code was written using Python three and OpenCV (Bradski, 2000). Feature detection on the petri dish was used to automatically determine the scale of video data. Background was subtracted and segmentation of the tadpole was carried out by thresholding. The resulting mask was fit to an ellipse to extract location and speed, as well as directional data and escape angle. Location data from 0 to 2 s following the onset of the looming stimulus was summarized and displayed as a contrail. Instantaneous velocity for the 3 s before and after the looming stimulus was also extracted and compared by area under curve analysis to obtain distance traveled.
Discrimination of positive and negative responses to looming stimuli
Request a detailed protocolAn automated python script was used to randomize tadpole videos and a user blinded to treatment categorized tadpole responses to looming stimuli as defensive escape behavior (positive response), absence of defensive escape behavior (negative response), or undeterminable (excluded from response rate calculation). Tadpole responses were categorized as undeterminable if the tadpole was moving quickly when the looming stimulus was presented. The response rate was calculated as the number of positive responses divided by the combined number of positive and negative responses.
Assessment of motor responses
Request a detailed protocolPositive escape responses were further analyzed to compare motor data. Escape distance and maximum velocity over 2 s from the onset of the loom was calculated. Escape angle was also measured from heading data. In dark looming stimuli trials, the escape angle was defined as the absolute value of the change in heading from 0 to 0.6 s. For bright looming stimuli trials, the escape angle was defined as the absolute value of the change in heading from 0 to 1.2 s. This difference is because dark looming stimuli evoked escape movement beginning approximately 0.5 s after looming onset, whereas bright looming stimuli evoked escape movement beginning approximately 1.0 s after looming onset.
Bioinformatics analysis
Identification of complement inhibitory proteins
Request a detailed protocolTo identify candidate complement inhibitory proteins, the STRING database (Szklarczyk et al., 2019) was used to look for functional interactions between complement C3 and other proteins. Human complement C3 was queried using a medium confidence (0.4) minimum interaction score, and a limit of 10 and 20 interactors for the first and second shell, respectively. Sources was limited to textmining, experiments and databases.
Expression of complement inhibitory proteins in human cortical cells
Request a detailed protocolThe Allen Brain Map Human Transcriptomics Cell Types Database (Allen Institute for Brain Science, 2015) was used to examine the expression of candidate proteins identified from the STRING query. Cell taxonomy and hierarchical clustering was retained from previous analysis (Hodge et al., 2019). Gene expression of complement inhibitory proteins across distinct brain cell clusters was examined and plotted using R software (R Development Core Team, 2018) and the Bioconductor (Huber et al., 2015), ComplexHeatmap (Gu et al., 2016), and ColorBrewer (Neuwirth, 2014) packages.
Homology search
Request a detailed protocolThe protein sequence of human CD46 (NCBI accession #: NP_002380.3) was queried in the Xenopus laevis genome with NCBI protein-protein BLAST (Altschul et al., 1990). The top three hits were predicted proteins XP_018102062.1, XP_018100169.1 and XP_018102058.1. The identity of the predicted proteins was obtained by cross-referencing with the Xenopus tropicalis genome using BLAST.
Conserved modular architecture analysis and transmembrane topology
Request a detailed protocolThe SMART protein domain research tool (Letunic and Bork, 2018; Schultz et al., 1998) was used to predict conserved modular architecture and the Phobius transmembrane topology tool (Käll et al., 2004; Käll et al., 2007) was used to predict transmembrane topology.
Tertiary protein structure prediction
Request a detailed protocolDe novo tertiary protein structure was predicted using trRosetta (Yang et al., 2020). The extracellular region, transmembrane region, and intracellular regions of human CD46 and Xenopus laevis aRCA3 were modeled separately by trRosetta and joined using pyMOL (Schrödinger, LLC, 2020) The orientation of the protein in the phospholipid bilayer was predicted using the protein bilayer lipid membrane orientation package (Kimmett et al., 2014).
Kinase prediction
Request a detailed protocolThe NetPhos 3.1 tool (Blom et al., 2004) was used to make predictions of tyrosine phosphorylation sites on the intracellular C-terminus of aRCA3. A cutoff score of 0.6 was used.
RNAscope in situ hybridization
Tissue preparation
Request a detailed protocolThe RNAscope (Wang et al., 2012) fixed-frozen tissue sample preparation and pretreatment protocol provided by the manufacturer was modified for Xenopus laevis tadpoles. Stage 46 tadpoles were euthanized in 0.2% MS-222 in 0.1X MBSH. Tadpoles were then fixed in 4% PFA at 4°C for 24 hr on a laboratory rocker. For cryoprotection, tadpoles were then moved sequentially through 10%, 20% and 30% sucrose in 1X PBS, until samples sunk to the bottom of the container. Cryoprotected tadpoles were then embedded in OCT blocks on dry ice. OCT blocks were sectioned at 8 μm thickness on a cryostat and mounted on superfrost plus slides. To enhance tissue adhesion, slides were air dried for 2 hr at −20°C and baked at 60°C for 30 min. Slides were then post-fixed by immersion in 4% PFA for 15 min, and then dehydrated with 50%, 70% and 100% ethanol. Slides were treated with hydrogen peroxide to block endogenous peroxidases. For target retrieval, using a hot plate and beaker, slides were boiled in RNAscope target retrieval reagent and then treated with RNAscope Protease III.
RNAscope assay
Request a detailed protocolThe RNAscope Multiplex Fluorescent 2.0 Assay was performed according to manufacturer’s protocols using the HybEZ oven. In brief, probes were applied to slides, and three amplification steps were carried out. Opal 570 dye was applied to slides along with Hoechst counterstain. Coverslips were mounted with Prolong gold antifade mountant. Slides were imaged with a Leica TCS SP8 confocal.
Molecular biology
The following plasmid constructs were constructed and used in this study: pEF1α-pHtdGFP; pEF1α-SYP-pHtdGFP; pEF1α-aRCA3-Myc-P2A-pHtdGFP; pEF1α-aRCA3-mCherry-Myc-P2A-SYP-pHtdGFP; pEF1α-pHtdGFP-P2A-Myc-VAMP2-C3; pEF1α-pHtdGFP-P2A-VAMP2.
cDNA library
Request a detailed protocolStage 40 Tadpoles were euthanized in 0.2% MS-222 in 0.1X MBSH. Animals were transferred to TRIzol reagent (Invitrogen, Carlsbad, CA) and homogenized by sonication and processed according to manufacturer’s protocols to isolate mRNA. Superscript IV reverse transcriptase (Invitrogen) was then used according to manufacturer’s protocol to generate whole tadpole cDNA.
Cloning and plasmid isolation
Request a detailed protocolCloning of Xenopus laevis cDNA was carried out using primers flanking genes of interest. Primers were generated according to mRNA sequences for aRCA3 (XM_018246573.1), VAMP2.S (NM_001087474.1), and c3.L (XM_018253729.1) predicted by NCBI (NCBI Resource Coordinators, 2018) and Xenbase (Karimi et al., 2018). Synaptophysin was subcloned from mouse synaptophysin (Ruthazer et al., 2006). pHtdGFP was subcloned from a pFA6a plasmid (Roberts et al., 2016). PCR was performed with Phusion High-Fidelity DNA polymerase (Thermo Scientific, Waltham, MA) and DNA fragments were cut with NEB enzymes and ligated into a plasmid with an EF1α promotor, ampicillin resistance, and a P2A self-cleaving construct or Myc tag when appropriate. DH5α bacteria were transformed and single clone colonies containing inserts were isolated. Endotoxin-free plasmid preparations of high yield and purity for electroporation was prepared by maxi-prep (Qiagen, Hilden, Germany).
Design of fusion proteins
Request a detailed protocolThe SYP-pHtdGFP fusion protein was created by in-frame fusion of pHtdGFP to the C-terminus of SYN between a three amino acid linker (QGT). The aRCA3-mCherry-Myc fusion protein was created by an in-frame fusion of mCherry-Myc to the C-terminus of aRCA3 between a two amino acid linker (AC). The Myc-VAMP2-C3 fusion protein was created by an in-frame fusion of Myc to the N-terminus of VAMP2 between a four amino acid linker (PGKI) and an in-frame fusion of C3 to the C-terminus of VAMP2 between an 18 amino acid linker (ASIKSPVQPLSAHSPVCI). The long linker for VAMP2-C3 was chosen to ensure that the C-terminus transmembrane domain of VAMP2 would not sterically hinder proper folding of C3.
Statistical analysis
Statistical analysis was performed using Graphpad Prism nine software. All data is presented as mean ± SEM, unless otherwise noted. Data was assessed for normality by D’Agostino-Pearson omnibus test and examination of quantile-quantile plots and assessed for homoscedasticity by examination of residual plots. All statistical tests were two-sided, and differences were considered significant for p<0.05.
Pairwise comparisons
Request a detailed protocolWhen normality and homoscedasticity assumptions were met, differences between two groups were tested by unpaired t-test. Welch’s t-test was used if the assumption of homogeneity of variance was not met. Non-normally distributed pairwise comparisons were tested by Mann-Whitney test.
Group data
Request a detailed protocolIf no datapoints were missing, two-way repeated measures ANOVA was used to test group data. If an interaction was found, the interaction was reported and Sidak’s multiple comparison test was used to test between groups or time points as appropriate. If no interaction was found, then the main effect was reported. Conversely, if the data contained missing values (due to loss of animals during a time course), the mixed effect restricted maximum likelihood model was used.
Correlation
Request a detailed protocolCorrelations were calculated using Pearson’s correlation coefficient. Trendlines were plotted using simple linear regression. Significance testing of the difference between two correlation coefficients was carried out using Fisher’s z transformation.
Non-linear curve fitting
Request a detailed protocolLog-transformed histograms were fit by Poisson regression with no weighting. The means of the best-fit curves were compared by likelihood-ratio test.
Data and code availability
Request a detailed protocolThe quantitative data used for analyses and figure generation in this paper are publicly available at https://doi.org/10.6084/m9.figshare.12895172.v1. The documented source code and user guide for the looming stimulus presentation and tadpole tracking software module, XenLoom (Beta), developed in this study are available at https://github.com/tonykylim/XenLoom_beta.
Data availability
Data used for analysis and figure generation are included in the manuscript and supporting files. Source data files have been provided for Figures 1–8.
References
-
Assembly of presynaptic active zones from cytoplasmic transport packetsNature Neuroscience 3:445–451.https://doi.org/10.1038/74814
-
Basic local alignment search toolJournal of Molecular Biology 215:403–410.https://doi.org/10.1016/S0022-2836(05)80360-2
-
Complement targets newborn retinal ganglion cells for phagocytic elimination by microgliaThe Journal of Neuroscience 39:2025–2040.https://doi.org/10.1523/JNEUROSCI.1854-18.2018
-
Role of membrane cofactor protein (CD46) in regulation of C4b and C3b deposited on cellsThe Journal of Immunology 168:6298–6304.https://doi.org/10.4049/jimmunol.168.12.6298
-
Sensorimotor decision making in the zebrafish tectumCurrent Biology 25:2804–2814.https://doi.org/10.1016/j.cub.2015.09.055
-
A guided tour into subcellular colocalization analysis in light microscopyJournal of Microscopy 224:213–232.https://doi.org/10.1111/j.1365-2818.2006.01706.x
-
The isolectin IB4 binds RET receptor tyrosine kinase in microgliaJournal of Neurochemistry 126:428–436.https://doi.org/10.1111/jnc.12209
-
Microglial phagocytosis of live neuronsNature Reviews Neuroscience 15:209–216.https://doi.org/10.1038/nrn3710
-
Sensors and regulators of intracellular pHNature Reviews Molecular Cell Biology 11:50–61.https://doi.org/10.1038/nrm2820
-
Integrin requirement for hippocampal synaptic plasticity and spatial memoryThe Journal of Neuroscience 23:7107–7116.https://doi.org/10.1523/JNEUROSCI.23-18-07107.2003
-
A CANDLE for a deeper in vivo insightMedical Image Analysis 16:849–864.https://doi.org/10.1016/j.media.2012.01.002
-
Microglia regulate the number of neural precursor cells in the developing cerebral cortexJournal of Neuroscience 33:4216–4233.https://doi.org/10.1523/JNEUROSCI.3441-12.2013
-
Development of microglia in the prenatal rat HippocampusThe Journal of Comparative Neurology 377:70–84.https://doi.org/10.1002/(SICI)1096-9861(19970106)377:1<70::AID-CNE7>3.0.CO;2-G
-
Visual avoidance in Xenopus Tadpoles Is correlated with the maturation of visual responses in the optic tectumJournal of Neurophysiology 101:803–815.https://doi.org/10.1152/jn.90848.2008
-
Collision detection as a model for sensory-motor integrationAnnual Review of Neuroscience 34:1–19.https://doi.org/10.1146/annurev-neuro-061010-113632
-
The applications of vitamin E TPGS in drug deliveryEuropean Journal of Pharmaceutical Sciences 49:175–186.https://doi.org/10.1016/j.ejps.2013.02.006
-
Orchestrating high-throughput genomic analysis with bioconductorNature Methods 12:115–121.https://doi.org/10.1038/nmeth.3252
-
A combined transmembrane topology and signal peptide prediction methodJournal of Molecular Biology 338:1027–1036.https://doi.org/10.1016/j.jmb.2004.03.016
-
Advantages of combined transmembrane topology and signal peptide prediction--the phobius web serverNucleic Acids Research 35:W429–W432.https://doi.org/10.1093/nar/gkm256
-
Xenbase: a genomic, epigenomic and transcriptomic model organism databaseNucleic Acids Research 46:D861–D868.https://doi.org/10.1093/nar/gkx936
-
Quantitating the subtleties of microglial morphology with fractal analysisFrontiers in Cellular Neuroscience 7:3.https://doi.org/10.3389/fncel.2013.00003
-
Excitation and inhibition in recurrent networks mediate collision avoidance in Xenopus tadpolesEuropean Journal of Neuroscience 40:2948–2962.https://doi.org/10.1111/ejn.12664
-
ProBLM web server: protein and membrane placement and orientation packageComputational and Mathematical Methods in Medicine 2014:1–7.https://doi.org/10.1155/2014/838259
-
20 years of the SMART protein domain annotation resourceNucleic Acids Research 46:D493–D496.https://doi.org/10.1093/nar/gkx922
-
SoftwareXenLoom_beta, version swh:1:rev:b487791d1d91a5950eebfd1e7640e0c3db761cf5Software Heritage.
-
CD46 (membrane cofactor protein) associates with multiple beta1 integrins and tetraspansEuropean Journal of Immunology 30:900–907.https://doi.org/10.1002/1521-4141(200003)30:3<900::AID-IMMU900>3.0.CO;2-X
-
Colony-stimulating factor 1 receptor (CSF1R) signaling in injured neurons facilitates protection and survivalJournal of Experimental Medicine 210:157–172.https://doi.org/10.1084/jem.20120412
-
Measurement of co-localization of objects in dual-colour confocal imagesJournal of Microscopy 169:375–382.https://doi.org/10.1111/j.1365-2818.1993.tb03313.x
-
Microglia regulate pruning of specialized synapses in the auditory brainstemFrontiers in Neural Circuits 13:55.https://doi.org/10.3389/fncir.2019.00055
-
Technical advance: autofluorescence as a tool for myeloid cell analysisJournal of Leukocyte Biology 88:597–603.https://doi.org/10.1189/jlb.0310184
-
Approach sensitivity in the retina processed by a multifunctional neural circuitNature Neuroscience 12:1308–1316.https://doi.org/10.1038/nn.2389
-
Visualization of the dynamics of synaptic vesicle and plasma membrane proteins in living axonsJournal of Cell Biology 140:659–674.https://doi.org/10.1083/jcb.140.3.659
-
Database resources of the national center for biotechnology informationNucleic Acids Research 46:D8–D13.https://doi.org/10.1093/nar/gkx1095
-
Phagocytosis by receptors for C3b (CR1), iC3b (CR3), and IgG (Fc) on human peritoneal macrophagesJournal of Leukocyte Biology 38:267–278.https://doi.org/10.1002/jlb.38.2.267
-
BookNormal Table of Xenopus Laevis (Daudin): A Systematical and Chronological Survey of the Development From the Fertilized Egg Till the End of MetamorphosisGarland Pub.
-
Apoptosis in neural development and diseaseAnnual Review of Neuroscience 23:73–87.https://doi.org/10.1146/annurev.neuro.23.1.73
-
Three-dimensional structure of a complement control protein module in solutionJournal of Molecular Biology 219:717–725.https://doi.org/10.1016/0022-2836(91)90666-T
-
Escape behavior and neuronal responses to looming stimuli in the crab Chasmagnathus granulatus (Decapoda: grapsidae)Journal of Experimental Biology 210:865–880.https://doi.org/10.1242/jeb.02707
-
PsychoPy2: experiments in behavior made easyBehavior Research Methods 51:195–203.https://doi.org/10.3758/s13428-018-01193-y
-
C1q: the perfect complement for a synaptic feast?Nature Reviews Neuroscience 9:807–811.https://doi.org/10.1038/nrn2394
-
Microglia shape corpus callosum axon tract fasciculation: functional impact of prenatal inflammationEuropean Journal of Neuroscience 39:1551–1557.https://doi.org/10.1111/ejn.12508
-
SoftwareR: A Language and Environment for Statistical ComputingR Foundation for Statistical Computing, Vienna, Austria.
-
CD46: expanding beyond complement regulationTrends in Immunology 25:496–503.https://doi.org/10.1016/j.it.2004.07.004
-
Identification and characterisation of a pH-stable GFPScientific Reports 6:1–9.https://doi.org/10.1038/srep28166
-
Stabilization of axon branch dynamics by synaptic maturationJournal of Neuroscience 26:3594–3603.https://doi.org/10.1523/JNEUROSCI.0069-06.2006
-
Bulk electroporation of retinal ganglion cells in live Xenopus tadpolesCold Spring Harbor Protocols 2013:prot076471.https://doi.org/10.1101/pdb.prot076471
-
In vivo Time-Lapse imaging of neuronal development in XenopusCold Spring Harbor Protocols 2013:top077156.https://doi.org/10.1101/pdb.top077156
-
Specificity of the thioester-containing reactive site of human C3 and its significance to complement activationBiochemical Journal 302 ( Pt 2:429–436.https://doi.org/10.1042/bj3020429
-
Fiji: an open-source platform for biological-image analysisNature Methods 9:676–682.https://doi.org/10.1038/nmeth.2019
-
Fluorescent proteins for investigating biological events in acidic environmentsInternational Journal of Molecular Sciences 19:1548.https://doi.org/10.3390/ijms19061548
-
Microglia: brain cells on the moveProgress in Neurobiology 178:101612.https://doi.org/10.1016/j.pneurobio.2019.04.001
-
Maternal immune activation in neurodevelopmental disordersDevelopmental Dynamics 247:588–619.https://doi.org/10.1002/dvdy.24612
-
Microglia modulate wiring of the embryonic forebrainCell Reports 8:1271–1279.https://doi.org/10.1016/j.celrep.2014.07.042
-
The complement system: an unexpected role in synaptic pruning during development and diseaseAnnual Review of Neuroscience 35:369–389.https://doi.org/10.1146/annurev-neuro-061010-113810
-
Development of ramified microglia from early macrophages in the zebrafish optic tectumDevelopmental Neurobiology 73:60–71.https://doi.org/10.1002/dneu.22039
-
A visual pathway for Looming-Evoked escape in larval zebrafishCurrent Biology 25:1823–1834.https://doi.org/10.1016/j.cub.2015.06.002
-
RNAscope: a novel in situ RNA analysis platform for formalin-fixed, paraffin-embedded tissuesThe Journal of Molecular Diagnostics : JMD 14:22–29.https://doi.org/10.1016/j.jmoldx.2011.08.002
-
Rapid innate defensive responses of mice to looming visual stimuliCurrent Biology 23:2011–2015.https://doi.org/10.1016/j.cub.2013.08.015
-
ScatterJ: an ImageJ plugin for the evaluation of analytical microscopy datasetsJournal of Microscopy 261:148–156.https://doi.org/10.1111/jmi.12187
Article and author information
Author details
Funding
Canadian Institutes of Health Research (FDN-143238)
- Edward S Ruthazer
Fonds de Recherche du Québec - Santé (31036)
- Edward S Ruthazer
Canadian Institutes of Health Research (Postdoctoral fellowship)
- Tony KY Lim
Faculty of Medicine, McGill (Postdoctoral Fellowship)
- Tony KY Lim
The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.
Acknowledgements
We thank Dr Wayne Sossin, Dr Jean-Francois Cloutier and the MNI Microscopy Core Facility for sharing access to equipment. We also thank Dr Philip Kesner and Anne Schohl for technical advice, Dr Larissa Ferguson for technical assistance, and all members of the Ruthazer lab for helpful discussions. The pFA6a-pHtdGFP plasmid was generously provided by Dr Joerg Stelling. PLX5622 was kindly provided by Plexxikon, Inc TKL was supported by a CIHR Postdoctoral Fellowship and the McGill Faculty of Medicine McLaughlin Postdoctoral Fellowship. This work was supported by grants to ESR from FRQS (31036) and CIHR (FDN-143238).
Ethics
Animal experimentation: All experiments were conducted according to protocol application number 2015-7728 approved by The Animal Care Committee of the Montreal Neurological Institute and in accordance with Canadian Council on Animal Care guidelines.
Copyright
© 2021, Lim and Ruthazer
This article is distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use and redistribution provided that the original author and source are credited.
Metrics
-
- 3,992
- views
-
- 416
- downloads
-
- 77
- citations
Views, downloads and citations are aggregated across all versions of this paper published by eLife.
Download links
Downloads (link to download the article as PDF)
Open citations (links to open the citations from this article in various online reference manager services)
Cite this article (links to download the citations from this article in formats compatible with various reference manager tools)
Further reading
-
- Immunology and Inflammation
Chronic antigenic stimulation can trigger the formation of interleukin 10 (IL-10)-producing T-regulatory type 1 (TR1) cells in vivo. We have recently shown that murine T-follicular helper (TFH) cells are precursors of TR1 cells and that the TFH-to-TR1 cell transdifferentiation process is characterized by the progressive loss and acquisition of opposing transcription factor gene expression programs that evolve through at least one transitional cell stage. Here, we use a broad range of bulk and single-cell transcriptional and epigenetic tools to investigate the epigenetic underpinnings of this process. At the single-cell level, the TFH-to-TR1 cell transition is accompanied by both, downregulation of TFH cell-specific gene expression due to loss of chromatin accessibility, and upregulation of TR1 cell-specific genes linked to chromatin regions that remain accessible throughout the transdifferentiation process, with minimal generation of new open chromatin regions. By interrogating the epigenetic status of accessible TR1 genes on purified TFH and conventional T-cells, we find that most of these genes, including Il10, are already poised for expression at the TFH cell stage. Whereas these genes are closed and hypermethylated in Tconv cells, they are accessible, hypomethylated, and enriched for H3K27ac-marked and hypomethylated active enhancers in TFH cells. These enhancers are enriched for binding sites for the TFH and TR1-associated transcription factors TOX-2, IRF4, and c-MAF. Together, these data suggest that the TR1 gene expression program is genetically imprinted at the TFH cell stage.
-
- Genetics and Genomics
- Immunology and Inflammation
Systemic lupus erythematosus (SLE) is an autoimmune disease, the pathophysiology and genetic basis of which are incompletely understood. Using a forward genetic screen in multiplex families with SLE, we identified an association between SLE and compound heterozygous deleterious variants in the non-receptor tyrosine kinases (NRTKs) ACK1 and BRK. Experimental blockade of ACK1 or BRK increased circulating autoantibodies in vivo in mice and exacerbated glomerular IgG deposits in an SLE mouse model. Mechanistically, NRTKs regulate activation, migration, and proliferation of immune cells. We found that the patients’ ACK1 and BRK variants impair efferocytosis, the MERTK-mediated anti-inflammatory response to apoptotic cells, in human induced pluripotent stem cell (hiPSC)-derived macrophages, which may contribute to SLE pathogenesis. Overall, our data suggest that ACK1 and BRK deficiencies are associated with human SLE and impair efferocytosis in macrophages.