Abstract
Rare genetic diseases preponderantly affect the nervous system causing neurodegeneration to neurodevelopmental disorders. This is the case for both Menkes and Wilson disease, arising from mutations in ATP7A and ATP7B, respectively. The ATP7A and ATP7B proteins localize to the Golgi and regulate copper homeostasis. We demonstrate genetic and biochemical interactions between ATP7 paralogs with the conserved oligomeric Golgi (COG) complex, a Golgi apparatus vesicular tether. Disruption of Drosophila copper homeostasis by ATP7 tissue-specific transgenic expression caused alterations in epidermis, aminergic, sensory, and motor neurons. Prominent among neuronal phenotypes was a decreased mitochondrial content at synapses, a phenotype that paralleled with alterations of synaptic morphology, transmission, and plasticity. These neuronal and synaptic phenotypes caused by transgenic expression of ATP7 were rescued by downregulation of COG complex subunits. We conclude that the integrity of Golgi-dependent copper homeostasis mechanisms, requiring ATP7 and COG, are necessary to maintain mitochondria functional integrity and localization to synapses.
SIGNIFICANCE STATEMENT Menkes and Wilson disease affect copper homeostasis and characteristically afflict the nervous system. However, their molecular neuropathology mechanisms remain mostly unexplored. We demonstrate that copper homeostasis in neurons is maintained by two factors that localize to the Golgi apparatus, ATP7 and the conserved oligomeric Golgi (COG) complex. Disruption of these mechanisms affect mitochondrial function and localization to synapses as well as neurotransmission and synaptic plasticity. These findings suggest communication between the Golgi apparatus and mitochondria through homeostatically controlled cellular copper levels and copper-dependent enzymatic activities in both organelles.
Introduction
Rare diseases include genetic maladies that disproportionately affect the nervous system with phenotypes ranging from neurodegeneration to behavioral impairments (Sanders et al., 2019; Lee et al., 2020). The vast collection of mutated human genes and neurologic symptomatology offers a wide field for the discovery of novel cellular mechanisms necessary for neuronal function. Here, we focus on rare neurologic diseases affecting ATP7A, ATP7B, and the conserved oligomeric Golgi (COG) complex subunit genes, whose products localize to the Golgi complex at steady state (Kaler, 2011; Polishchuk and Lutsenko, 2013; Climer et al., 2015). We report that these ATP7 paralogs and the COG complex converge to maintain copper homeostasis and, unexpectedly, mitochondrial distribution in neurons.
Mutations in the ATP7A cause Menkes disease, a systemic copper-depletion affliction caused by intestinal copper malabsorption. Menkes is characterized by childhood systemic and neurologic phenotypes. Menkes' brain phenotypes result from this systemic copper depletion and span from intellectual disability to widespread gray matter neurodegeneration. (Menkes, 1988, 1999; Kaler, 2011; Zlatic et al., 2015; Hartwig et al., 2019; Guthrie et al., 2020). These Menkes disease phenotypes are caused in part by defects in the activity of copper-dependent enzymes that either traverse the Golgi complex or localize to mitochondria (Lutsenko et al., 2007; Zlatic et al., 2015). In contrast, with the systemic depletion of copper observed in Menkes patients, cell-autonomous ATP7A gene defects cause a cellular overload of copper because of impaired copper efflux from cells (Camakaris et al., 1980; Morgan et al., 2019). This ATP7A cell-autonomous copper phenotype can be recapitulated systemically in Wilson disease. Wilson disease is caused by mutations in the ATP7A paralogue, ATP7B, a copper transporter gene expressed preferentially in the liver. Wilson disease leads to organismal copper overload because of defective copper excretion by the liver. This systemic copper overload causes liver damage, psychiatric symptoms, and lenticular neurodegeneration at late stages of disease (Lutsenko et al., 2007; Kaler, 2011).
ATP7A and ATP7B are P-type ATPases residing in the Golgi complex at steady state where they contribute to cellular copper homeostasis. These ATPases sequester copper into the Golgi lumen away from the cytoplasm by a mechanism that requires ATP hydrolysis to drive copper transport. Furthermore, after a copper challenge, these ATPases are translocated to the plasma membrane where they pump copper outside of the cell (Petris et al., 1996; Kaler, 2011; Polishchuk and Lutsenko, 2013; Polishchuk et al., 2014). In turn, extracellular copper is taken up by the activity of a plasma membrane copper transporter CTR1 (Kuo et al., 2001; Lutsenko et al., 2007; Kaler, 2011). The expression of ATP7A and CTR1 is controlled by the activity of the COG complex, which tethers vesicles to the Golgi complex (Comstra et al., 2017). The COG complex is an octamer localized to the Golgi apparatus and required for incoming vesicle fusion with the Golgi. Elimination of any one of the eight COG subunits leads to Golgi fragmentation because of destabilization and degradation of the whole octamer (Ungar et al., 2002; Zolov and Lupashin, 2005; Climer et al., 2015; Bailey Blackburn et al., 2016). Human mutations in seven of the eight COG subunits cause a group of diseases collectively known as congenital disorder of glycosylation Type II (Wu et al., 2004; Kranz et al., 2007; Foulquier, 2009; Climer et al., 2015). These rare disorders are characterized, in part, by neurodevelopmental pathology and behavioral phenotypes which are loosely similar to neurologic phenotypes in Menkes disease. However, whether congenital disorder of glycosylation Type II neurologic phenotypes may be linked to copper metabolism remains unanswered (Climer et al., 2018).
The convergence of neurologic, neurodevelopmental, and psychiatric phenotypes in these three diseases with proteins localizing to the Golgi complex, point to fundamental, neuronal mechanisms necessary to maintain cellular copper homeostasis. These phenotypes raise the question: how are the activities of copper transporters and the COG complex coordinated to maintain copper homeostasis in cells, tissues, and organisms? Here, we demonstrate that ATP7 paralogues and the COG complex, both Golgi localized machineries, control synapse development, neurotransmission, and the subcellular localization of mitochondria at synapses.
Material and Methods
Drosophila strains
The following strains where obtained from the Bloomington Drosophila Stock Center, Bloomington, IL: w1118 (#5905), pnr-GAL4 (#3039), Ddc-GAL4 (#7009), elavC155 -GAL4 (#458), UAS-COG1 RNAi (#38943), UAS-COG5 RNAi (#56882), UAS-COG8 RNAi (#55192), Cog1e02840 (#18089), Mito-GFP (#25748). UAS-ATP7-wt (OE) and UAS-ATP7-FLAG (OE1) were gifts of Richard Burke, Monash University, Australia. The UAS-ATP7 RNAi (#108159) stock was obtained from the Vienna Drosophila Resource Center. GAL4477, UAS-mCD8::GFP/CyO, tubP-GAL80;GAL4ppk.1.9, UAS-mCD8::GFP (CIV-GAL4 driver), GAL4477, UAS-mito-HA-GFP.AP/CyO;ppk::tdTomato, and OregonR (wild-type) were used for dendritic analysis. UAS-mito-HA-GFP.AP line was obtained from Bloomington (#8442).
Antibodies and beads
Antibodies used in Drosophila stainings were anti-HRP-FITC (MP Biomedicals, 0855977, RRID:AB_2334736), phalloidin-rhodamine (Abcam, ab235138), phalloidin-Alexa Fluor 633 (Invitrogen, A22284), anti-Cy3 (Santa Cruz Biotechnology, 166894), the primary antibody rabbit anti-GFP (Sysy 132002, RRID:AB_887725), and the secondary antibody Alexa Fluor 488 goat anti-rabbit (Invitrogen).
The following antibodies were used: rabbit anti COG7 from V. Lupashin's laboratory, rabbit anti-COG5 (HPA020300, Sigma, RRID:AB_2081406), mouse monoclonal anti-GOLGIN 97 CDF4 (A-21270, Thermo Fisher Scientific, RRID:AB_221447), mouse anti-Lamp1 (DSHB H4A3-c, Developmental Studies Hybridoma Bank, RRID:AB_528126), rabbit anti-ATP7B (AB124973, Abcam, RRID:AB_10975785), mouse anti-GAPDH (0411; sc-47 724, Santa Cruz, RRID:AB_627678), mouse anti-myc (c-myc (9E10), sc-40, Santa Cruz, RRID:AB_627268), mouse anti-FLAG M2 (F3165, Sigma, RRID:AB_259529), rabbit anti-COX11 (11 498–1-AP, Proteintech, RRID:AB_10598173), rabbit anti-COX17 (HPA048158, Sigma, RRID:AB_2680292), rabbit anti-MPC1 (HPA045119, Sigma, RRID:AB_10960421), mouse anti-MFN1 (ab126575, Abcam, RRID:AB_11141234), rabbit anti-Seladin/DHCR24 (10 471–1-AP, Proteintech, RRID:AB_2091446), rabbit anti-HSP60 (12 165, Cell Signaling, RRID:AB_2636980), rabbit anti-SDHA (11 998, Cell Signaling, RRID:AB_2750900), and mouse monoclonal F-10 anti-Tomm20 (sc-17 764, Santa Cruz, RRID:AB_628381). The following beads were used in immunoprecipitations: protein A-Sepharose 4 B Fast Flow from Staphylococcus aureus (P9424, Sigma), Protein G Sepharose four Fast Flow (17-0618-01, GE Healthcare), Dynal magnetic beads (110.31, Invitrogen).
Drosophila immunofluorescence microscopy
For experiments where larval neuromuscular junctions were analyzed, Drosophila stocks were reared in polystyrene vials on standard fly food (900 ml milli-Q water, 48 g active dry yeast, 120 g cornmeal, 9 g agar, 120 g molasses, 2.4 g tegosept, and 9 ml propionic acid). Animals were housed in a humidified, 25°C incubator (Jeio Tech Co, Inc, IB-15G) with a 12/12 h light/dark cycle.
Larval NMJ dissections, immunohistochemistry, and confocal microscopy were performed as described previously (Gokhale et al., 2015b, 2016; Mullin et al., 2015). Wandering late third-instar female larvae were dissected with a dorsal incision in a 10 mm cell culture dish in normal HL3 composed of 70 mm NaCl, 5 mm KCl, 21.5 mm MgCl2, 10 mm NaHCO3, 5 mm trehalose, 115 mm sucrose, and 5 mm BES, pH 7.2–7.3. Larvae were fixed in 4% paraformaldehyde for 1 h at room temperature (RT) and stained with anti-HRP-FITC (1:500) and anti-phalloidin-rhodamine (1:1000; data not shown, used for muscle orientation) in PBS supplemented with 0.1% Triton X-100 (PBS-T) overnight at 4°C. All steps were followed by three 10 min rinses in PBS-T and a final rinse in PBS. The larval preparation was affixed to a poly-L-lysine coated glass side with a drop of Vectashield (Vector Labs H-1000-10) and coverslipped. The coverslip was adhered to the slide with a layer of nail polish in preparation for inverted microscopy. Samples were stored at 4°C until imaged. An inverted Olympus FV1000 microscope was used for confocal imaging of synapses on muscle 6/7 in the A3 segment. Maximum intensity projections were converted to jpegs. Total NMJ length and bouton perimeter measurements were calculated in FIJI (Schindelin et al., 2012). The FIJI Cell Counter plugin was used to quantify bouton number (Rai et al., 2018). All files were blind for quantification.
For body wall mosaics, larvae were dissected, fixed, and stained as described above; however, a lateral incision was made to preserve one set of the somatic musculature. The dissection was incubated in the primary antibodies anti-phalloidin-633 (1:500), anti-HRP-Cy3 (1:1000), and rabbit anti-GFP (1:1000) overnight at 4°C and incubated with Alexa Fluor 488 goat anti-rabbit (1:1000) with additional anti-phalloidin-633 (1:500), and anti-HRP-Cy3 (1:1000) for 2 h at RT the following day before placement on a slide. Mosaics were imaged using the Multi Area Time Lapse Controller of the Olympus FV1000 confocal microscope. Image stitching was completed, and mosaics converted to jpegs with the Multi Area Time Lapse Viewer of Olympus Fluoview software.
For live confocal imaging, GAL4477, UAS-mCD8::GFP;GAL4ppk.1.9, UAS-mCD8::GFP virgin flies were crossed to males from individual transgenic fly lines, or outcrossed to OregonR males which served as control. Fluorescent neurons from wandering third-instar larvae were imaged on a Zeiss LSM780 confocal microscope as previously described (Das et al., 2017). Briefly, individual animals were placed in a microscopic slide and immersed in 1:5 (v/v) diethyl ether to halocarbon oil and covered with a 22 × 50-mm coverslip. Images were collected as z-stacks using the 20× dry objective at a step-size of 2 µm. Maximum intensity projections of the images acquired were exported using the Zen blue software. The images were then manually curated using the Flyboys software to remove non-specific background noise such as denticle belts. The ImageJ software was used to skeletonize and process these images as previously described (Iyer et al., 2013a,b), and custom Python algorithms used to compile quantitative data such as total dendritic length and number of branches. For Strahler order analysis, images were reconstructed using NeuronStudio and the branch order extracted using the centripetal branch order function (Rodriguez et al., 2008).
For mitochondrial analysis, virgin female flies from GAL4477, UAS-mito-HA.GFP.AP/Cyo;ppk-tdTomato were outcrossed to either males from OregonR (control) or individual transgenic fly lines. Images were acquired as described above and maximum intensity projections of the z stacks were exported using Zen Blue software. Using Adobe Photoshop, images were cropped to a fixed size (6.6 × 4.5 inches) in the same quadrant of the neurons and the mitochondria manually counted along the dendrites in these quadrants. The cropped images were processed as described above to obtain the total dendritic length which was then used to analyze the number of organelles per 1000 pixels dendritic length. For branch point analysis, mitochondria localized at branch points were manually counted and normalized to the total number of branch points.
Drosophila viability
For experiments measuring animal viability, flies were reared on standard Molasses Food (Genesee Scientific) at 25°C in 12/12 h light/dark cycle. The total number of balanced (control) and unbalanced (experimental) male and female adult offspring from each cross was counted every 1–2 d for 9 d with at least 200 offspring scored. Viability was observed with at least three replicate vials per experiment. Four to five flies per genotype per sex, where available, were imaged at 4× magnification (Gokhale et al., 2015a; Zlatic et al., 2018).
Drosophila copper toxicity
Male and female flies were collected at eclosion and aged 4–7 d. Flies were then starved for 3 h before being placed in empty vials that contained a filter disk soaked with 200 μl of 100 mm CuSO4. Approximately 10 flies per sex per genotype were analyzed with a minimum of six replicates. Flies were maintained at 25°C in a 12/12 h light/dark cycle and the number of dead flies was assessed daily (Gokhale et al., 2015a; Zlatic et al., 2018).
Drosophila synaptic physiology
Drosophila were raised in supplemented fly food (German Formulation, Genesee Scientific), at 2°C temperature and 80% of humidity. Recordings were performed in Modified HL3 solution (10 mm NaHCO3, 5 mm KCl, 4 mm MgCl, 5 mm HEPES, 70 mm NaCl, 5 mm trehalose, and 115 mm sucrose; pH 7.2). The final Ca2+ concentration was set to the desired level indicated in the text. Postsynaptic currents from specific genotypes were recorded at segment A3 of ventral longitudinal muscle 6 as indicated in third-instar larvae using two-electrode voltage clamp with a –80-mV holding potential as was previously described (Acharya et al., 2006; Jorquera et al., 2012; Stevens et al., 2012; Astorga et al., 2016). The motor nerves innervating the specified muscles were severed and placed into a suction electrode, so an action potential could be evoked by electrical stimulation at the indicated frequencies using a programmable stimulator (Master8, AMPI). Each stimulus was adjusted to 6-V intensity for 0.1-ms duration, allowing to recruit both motor neurons and preventing electrotonic induced hyperexcitability.
Data Acquisition and analysis was performed using the National Instrument DAQ system controlled by WinEDR software (University of Strathclyde) on a desktop computer. Electrical recordings were analyzed using Clampfit 10 software (Molecular Devices) or IgorPro 4.0 software (WaveMetrics). Graphs and statistics were generated using OriginPro 9.0 software (OriginLab Corp), and Microsoft Excel routines.
To quantify the post-tetanic plasticity in the different genotypes, we adjusted the normalized post-tetanic responses with a previously described model that adds depression in the function as P·e⋀(–x/τ) + 1 – D, were P is the post-tetanic potentiation (PTP) and D the post-tetanic depression (PTD), which decays between two states with a time constant tau. To solve this function, we used numerical solutions by repetitive iteration with solver routines minimizing the summation of the squares of the differences between the data and the model. To help the simulation, we incorporated values of potentiation and depression obtained during the tetanic stimulation. Numerical analysis of normalized post-tetanic enhancement in control and mutant animals generates values for potentiation, depression, and the time constant.
Quantitation of total copper
The measurement of total copper was conducted as previously described (McAllum et al., 2020). Briefly, quantitation of copper metal in larva were determined with a Thermo iCAP-TQ series inductively coupled plasma mass spectrometer (ICP-MS). The instrument was operated in kinetic energy discrimination (KED) mode using He as a collision gas to remove polyatomic interferences. Quantitation of copper (as m/z 63) was determined with nine-point external calibration curve (Accustandard, AG-MECAL2A-ASL-5) generated using 1000 and 100 µg l−1 stock solutions of standard with online dilutions conducted by the prepFast M5 system (Elemental Scientific) and ranged from 2 to 1000 µg l−1. The yttrium internal standard (10 µg l−1, Accustandard, AG-INT2-ASL-5) was mixed with the sample online with the prepFast system. The wash, sample dilution, and sample carrier solutions were 2% nitric acid in 18.2-Ω water. The wet weight of each larva was determined, and each larva were digested at 95°C for 10 min in 20 µl of concentrated (70%) nitric acid (metal basis grade-225711, Aldrich) followed by overnight at RT. The samples were then diluted 30× with 2% nitric acid to a final volume of 600 µl. This solution was transferred to a polypropylene 96-deep-well plate (Greiner Bio-One GmbH, 780201), and samples were measured directly from the 96-well plate with M5 prepFast autosampler using 250-µl sample loops.
Cell culture
SH-SY5Y (ATCC, RRID:CVCL_0019) and HEK293T (RRID:CVCL_0063) cells were cultured in DMEM supplemented with 10% fetal bovine serum (FBS) and 100 μg/ml penicillin and streptomycin at 37°C in 10% CO2. HEK293T cells deficient for COG subunits were generated as described (Blackburn and Lupashin, 2016; Willett et al., 2016). HepG2 cells (RRID:CVCL_0027) were grown in DMEM supplemented with 10% FBS (decomplemented at 56°C for 30 min), 2 mm L-glutamine, 100 μg/ml penicillin, and streptomycin at 37°C in 10% CO2. HepG2 copper incubations were performed by adding 200 μm CuSO4 or 400 μm BCS diluted in DMEM at 37°C in a CO2 incubator for 2 h.
SH-SY5Y (ATCC, RRID:CVCL_0019) were genome edited using gRNA and Cas9 preassembled complexes by Synthego with a knock-out efficiency of 80%. The gRNAs used was ACAGACUCCAAAGACCCUAC which targeted transcript ENST 00000341514 exon 3. Wild-type and mutant cells were cloned by clonal dilution and clones mutagenesis was confirmed by Sanger sequencing.
Immunoprecipitation of ATP7 paralogues
SH-SY5Y cells were grown to 80–90% confluency in 10-cm tissue culture dishes. On the day of the experiment, cells were washed twice in ice cold PBS followed by lysis in a buffer containing 150 mm NaCl, 10 mm HEPES, 1 mm EGTA, and 0.1 mm MgCl2, with 0.5% Triton X-100 and Complete Protease Inhibitor (Roche 11245200), pH 7.4. The lysate was sonicated and incubated on ice for 30 min followed by a clarifying spin at 16,100 × g for 10 min. The soluble supernatant was obtained and 500 μg of the protein lysate was added to 30 μl Dynal magnetic beads coated with 5 μl ATP7A antibody (75–142, NeuroMab, RRID: AB_10672736) and incubated for 2 h at 4°C in an end-to-end rotor. As negative controls, immunoprecipitations were done in the presence of a custom made antigenic ATP7A peptide (VSLEEKNATIIYDPKLQTPK, Biosynthesis) prepared in 10 mm MOPS and used at 22 μm (Comstra et al., 2017). Additionally as a non-specific antibody control, we used equivalent amounts of the monoclonal antibody against FLAG (F3165, Sigma, RRID::AB_259529). After the incubation period, magnetic beads were washed six times in a buffer containing 150 mm NaCl, 10 mm HEPES, 1 mm EGTA, and 0.1 mm MgCl2, pH 7.4 with 0.1% Triton X-100. Proteins were eluted from the beads with Laemli sample buffer. Samples were resolved by SDS-PAGE and contents analyzed by immunoblot.
For immunoprecipitations using HepG2 cells, cells were solubilized at 4°C for 10 min at RT in a lysis buffer containing 0.5% Triton X-100, 20 mm Tris/HCl (pH 7.4), 150 mm NaCl, 1 mm EDTA (pH 8), 0.5% NP-40, and 10% glycerol, supplemented with 1× protease inhibitor cocktail (Sigma). The mixture was placed into a Microfuge tube, kept on ice for 10 min, and then centrifuged at 16,100 × g for 15 min at 4°C. Cell lysates were incubated with anti-ATP7B antibody overnight or with anti-FLAG and anti-Myc antibodies as controls. Then protein A or G Sepharose beads were added to each specimen for 4 h, and immune complexes were collected by centrifugation. The beads were then washed and immunoprecipitated proteins were eluted with sample buffer, separated by SDS-PAGE and analyzed by Western blotting.
RNA interference
Small interfering RNAs (siRNAs) targeting COG5 were purchased from Sigma-Aldrich. The following siRNA were used: CCAUUCAUAGUUGGCAGUU and GGACUUUGAAGGAUACUA.
Scrambled siRNAs were used as a control. HepG2 cells were transfected with siRNAs using Dharmafect4 (Dharmacon T-2004), according to manufacturer instructions.
RNA preparation and qRT-PCR
Control cells and COG5-silenced cells were lyzed and total RNA purified by QIAshredder (QIAGEN) and extracted with RNeasy Protect Mini kit (QIAGEN) using standard conditions. Total RNA (1 μg) was reverse-transcribed by QuantiTect Reverse Transcription kit (QIAGEN) according to the manufacturer's instructions. qRT-PCR experiments were performed using Light Cycler 480 Syber Green MasterMix (Roche) for cDNA amplification and LightCycler 480 II (Roche) for signal detection. qRT-PCR results were analyzed using the comparative Ct method normalized against housekeeping gene β-Actin. The following primers were used: β-ACTIN forward (5′-AAGAGCTACGAGCTGCCTGA-3′), β-ACTIN reverse (5′-GACTCCATGCCCAGGAAGG-3′), COG5 forward (5′-AGAGCCCGACTTGAAGTGGAAA-3′), and COG5 reverse (5′-ACCTGAAGAGCTGTTCCGACTT-3′).
Immunofluorescence
For HEPG2 cells, control and COG5-silenced cells were fixed (10 min, 4% paraformaldehyde in 0.2 m HEPES) directly or 2 h after incubation with 200 μm CuSO4. Fixed cells were then incubated with blocking/permeabilizing solution (0.5% bovine serum albumin, 0.1% saponin, 50 mm NH4Cl in PBS) for 20–30 min. Primary and secondary antibodies were diluted in blocking/permeabilizing solution and added to the cells for 1 h and 45 min, respectively. Samples were examined with a Zeiss LSM 700 confocal microscope equipped with a 63× 1.4 NA oil objective.
For HEK293 cells, control and COG1 mutant cells were plated onto Matrigel coated coverslips and incubated at 37°C, 5% CO2 overnight in DMEM (Corning, 10-013-CV) supplemented with 10% FBS (97068-085, VWR). Cells were washed twice in 37°C warmed PBS (21-040-CV, Corning) supplemented with 1 mm MgCl2 and 100 μm CaCl2 and fixed at 37°C for 20 min in warmed, 4% paraformaldehyde (P6148, Sigma) in PBS/1 mm MgCl2/100 μm CaCl2. Fixed cells were blocked and permeablized for 30 min at RT in block solution containing, 2% BSA (03116956001, Roche), 1% fish gelatin (G7765, Sigma), 0.02% saponin (S4521, Sigma), and 15% horse serum (SH30074.03, HyClone) in PBS. Cells were incubated with primary antibody the mouse monoclonal anti-Tomm20 antibody diluted 1:200 in block solution at 37°C for 30 min. Coverslips were washed three times in Block Solution and incubated with secondary antibody [Alexa Fluor 555 goat anti-mouse IgG (H + L), A21422, Invitrogen, RRID:AB_2535844] diluted 1:1000 in block solution at 37°C for 30 min. Coverslips were then washed twice in Blocking Solution, once in PBS/1 mm MgCl2/100 μm CaCl2 and mounted on slides with Dapi Fluoromount-G (0100-20, Southern Biotech).
Coverslips were imaged on a Nikon Ti2-E inverted microscope 60× oil immersion Apo objective, NA 1.40 with a Galvano scanner. DAPI acquisition settings include excitation/emission wavelength of 405.0/450.0, laser scan confocal with PMT detector. Alexa Fluor 555 acquisition settings include excitation/emission wavelength or 561.0/595.0 Laser Scan Confocal with a GaAsP detector. Images were captured with a Nikon A1 LFOV camera. NIS-Elements C Imaging v5.20.01 and Imaris software were used for data collection, analysis, and 3D reconstruction.
Extracellular flux analysis of tissue culture cells
All oxygen consumption rates were measured on the Seahorse XFe96 Analyzer (Seahorse Bioscience) following manufacturer recommendations. XFe96 extracellular flux assay kit probes (Seahorse Bioscience 102601-100) were incubated with the included calibration solution overnight at 37°C under non-CO2-injected conditions. Cells were trypsinized, counted (Bio-Rad TC20 automated Cell Counter), and seeded into XFe96 cell culture microplates (Seahorse Bioscience 101085-004) at 20,000 cells/well and incubated at 37°C with 10% CO2 in complete culture media for 24 h before initialization of the stress test. The next day, wells were washed twice in Seahorse stress test media consisting of Seahorse XF base media (Seahorse Bioscience 102353-100), 2 mm L-glutamine (HyClone SH30034.01), 1 mm sodium pyruvate (Sigma S8636), and 10 mm D-glucose (Sigma G8769), pH7.4 before final volume added and then incubated at 37°C in non-CO2-injected conditions for 1 h before stress test. Seahorse injection ports were filled with 10-fold concentrated solution of disulfiram (Sigma 86 720), oligomycin A (Sigma 75 351), carbonyl cyanide-4-(trifluoromethoxy)phenylhydrazone (FCCP; Sigma C2920), and rotenone (Sigma R8875)/antimycin A (Sigma A8674) for final testing conditions of disulfiram (0.1 μm), oligomycin (1.0 μm), FCCP (0.125 μm), rotenone (0.5 μm), and antimycin A (0.5 μm). The flux analyzer protocol conditions consisted of three basal read cycles, 40 reads following disulfiram injection, and three reads following each injection of oligomycin A, FCCP, and rotenone plus antimycin A. Each read cycle consisted of a 3-min mix cycle, followed by a 3-min read cycle where oxygen and pH levels were determined over time. The Seahorse Wave Software version 2.2.0.276 was used for data analysis of oxygen consumption rates. Experiments were repeated in quadruplicate and aliquots of each batch of stress test media were sent for inherent copper concentration analysis by ICP mass spectrometry performed by the Center for Applied Isotope Studies (CAIS) at the University of Georgia. Non-mitochondrial respiration was determined as the lowest oxygen consumption rate following injection of rotenone plus antimycin A. Basal respiration was calculated from the oxygen consumption rate just before oligomycin injection minus the non-mitochondrial respiration. Acute response was calculated from the difference in oxygen consumption rates just before disulfiram injection to just before oligomycin injection. Oligomycin A sensitivity was calculated as the difference in oxygen consumption rates just before oligomycin injection to the minimum oxygen consumption rate following oligomycin injection but before FCCP injection.
Evolutionary rate covariation analysis
We used the web engine https://csb.pitt.edu/erc_analysis/ (Clark et al., 2012, 2013; Findlay et al., 2014) using as inputs the CORUM annotated complexes database (Ruepp et al., 2010).
Statistical analysis
Experimental conditions were compared using Synergy KaleidaGraph (version 4.1.3) or Aabel NG2 v5 x64 (Gigawiz) as specified in each figure. Metal determinations were analyzed by estimation statistics using the web engine https://www.estimationstats.com/#/ (Ho et al., 2019).
Results
P-ATPase copper transporters interact with the COG complex
ATP7A and ATP7B are paralogs present in chordates, which originated from a single ATP7 gene found in the last common ancestor to eumetazoans. The Drosophila genus retained this single ATP7 gene (http://www.treefam.org/family/TF300460). We previously demonstrated that human ATP7A and the COG complex biochemically interact in cells (Comstra et al., 2017). Moreover, mutations of COG complex subunits reduced the expression of ATP7A and CTR1, resulting in cellular copper depletion (Comstra et al., 2017). As ATP7B shares 52.5% identity with ATP7A (alignment of Uniprot entries ATP7A_HUMAN and ATP7B_HUMAN), then ATP7B may also interact with the COG complex. To test this hypothesis, we immunoprecipitated endogenous ATP7A or ATP7B from neuroblastoma or hepatoma human cell lines, respectively, and assessed the presence of COG complex subunits by immunoblot with antibodies against COG5 and COG7 (Fig. 1). Both ATP7A and ATP7B immunoprecipitated COG complex subunits (Fig. 1A, lane 4, B, lanes 3 and 4). The selectivity of these ATP7A/B coprecipitations with COG subunits was determined with control immunoprecipitations with unrelated antibodies (Fig. 1A, lane 3, B, lanes 3′ and 4′), out-competing immune complexes with the ATP7A antigenic peptide (Fig. 1A, lane 5), or by blotting ATP7B immunoprecipitates with GAPDH antibodies (Fig. 1B). Similar to ATP7A (Comstra et al., 2017), the association of ATP7B with the COG complex remained after challenging cells with excess extracellular copper (Fig. 1B, compare lanes 3 and 4). Downregulation of COG5 in hepatoma cells (Fig. 1C) disrupted Golgi morphology and decreased ATP7B expression, both at Golgi and lysosomal compartments after copper-induced translocation, as revealed by co-labeling cells with ATP7B and GOLGIN95 or LAMP1 antibodies, respectively (Fig. 1D). These findings demonstrate a conserved interaction between the COG complex with ATP7 paralogs.
Biochemical and genetic interactions between COG complex subunits and ATP7 paralogs. SH-SY5Y (A, E, F) and HepG2 (B) cells were incubated in normal media or media supplemented with a 200 μm copper excess (B, even lanes). Detergent soluble extracts from (A) SH-SY5Y and (B) HepG2 cells were immunoprecipitated with antibodies against ATP7A or ATP7B, respectively. The presence of COG complex in the immunoprecipitates was detected with antibodies against COG5 and COG7. Controls were performed with unrelated antibodies (A, lane 3 or B, lanes 3′−4′) or with the ATP7A antigenic peptide (A, lane 5). Asterisk marks non-specific band. A, Lane 1 corresponds to input. B, Lanes 1,2,1′ and 2′ correspond to inputs. C, siRNA knock-down of COG5 in HepG2 cells. Relative expression corresponds to the ratio between COG5 and ACTB (n = 3). D, Immunomicroscopy analysis of the expression and localization of ATP7B in Golgi (untreated cells) and lysosomes (CuSO4-treated cells). These compartments were detected with antibodies against GOLGIN97 and LAMP1, respectively. Scale bar: 10 µm. E, ATP7A control (lanes 1–2) and null clone (lanes 3–4) cell extracts were immunoblotted with the indicated antibodies. F, Quantification of blots presented in E; n = 4 Wilcoxon–Mann–Whitney rank-sum test.
Since mutations or RNAi downregulation of COG complex subunits reduced the expression of ATP7A or ATP7B (Comstra et al., 2017), we next asked whether the expression of COG subunits and GPP130/GOLIM4, an established COG-dependent Golgi resident protein, were sensitive to elimination of ATP7A expression. To address this question, we engineered SH-SY5Y neuroblastoma clones lacking ATP7A by CRISPR-Cas9 mutagenesis. We chose these cells because their ATP7B expression was undetectable by mass spectrometry analysis (V. Faundez' unpublished data). These ATP7A null clones robustly and significantly increased the expression of the COG complex subunits COG5 and seven as well as a Golgi resident and COG-interacting protein GPP130/GOLIM4 (Fig. 1E,F; Oka et al., 2004; Zolov and Lupashin, 2005; Sohda et al., 2010; Comstra et al., 2017). These findings demonstrate bidirectional genetic control of ATP7 and COG complex machinery expression.
Synaptic development requires ATP7-dependent and COG-dependent copper homeostasis
We tested the phylogenetic conservation of the interactions between the COG complex and copper transporters in adult and developing Drosophila neurons. We took advantage of the fact that the expression of mammalian ATP7A, ATP7B, and CTR1 is reduced in cells that lack COG complex subunits (Fig. 1C,D; Comstra et al., 2017). We reasoned that Drosophila neurons would help us to genetically test whether the COG complex balances the expression of Drosophila ATP7 and CTR1 to control copper homeostasis. We predicted that neuronal cell autonomous and systemic changes in transporter expression should differentially and predictably affect copper content in developing tissues and animals (Fig. 2A). For example, neuronal cell autonomous ATP7 overexpression would decrease cellular copper by localizing this transporter to the plasma membrane, thus constituting a model of copper depletion by metal efflux (Fig. 2A4; Norgate et al., 2006; Burke et al., 2008; Binks et al., 2010). In contrast, copper content after systemic ATP7 overexpression at the plasma membrane would result from a balance between two effects (Fig. 2A2). First, an increased enterocyte copper absorption resulting in increased copper in the hemolymph bathing cells in the organism (Fig. 2A2). This enterocyte-hemolymph effect would be counterbalanced by increased cellular copper efflux from internal organ cells, such as neurons (Fig. 2A2; Bahadorani et al., 2010). Conversely, neuronal cell autonomous ATP7 knock-down would increase cellular copper because of defective copper efflux from neurons (Fig. 2A6; Norgate et al., 2006; Burke et al., 2008; Binks et al., 2010). However, systemic ATP7 inhibition would be dominated by enterocyte copper absorption defects, hemolymph copper deficiency, and organ copper depletion, thus mimicking Menkes disease (Fig. 2A3). We tested these copper content predictions by systemic overexpression of ATP7 or downregulation of ATP7 in all tissues using the Actin-GAL4 driver (Fig. 2B). Systemic overexpression of ATP7 did not significantly modified the larval copper content per mg of larva (Fig. 2B). In contrast, ubiquitously expressed ATP7 RNAi transgene decreased copper content per mg of larva although these larvae were significantly smaller (Fig. 2B). This ATP7 RNAi copper content phenotype partially mimics the organismal copper depletion phenotype observed in Menkes disease (Figs. 2A3,B).
Model of copper homeostasis by the COG complex and copper transporters and copper content in Drosophila. A, Model of predicted larval and brain copper status in different genotypes. Brain and gut are depicted in systemic genetic models, while only brain is depicted in the neuronal cell autonomous models. Color intensities of organs show predicted copper content (color code in the left of the figure). Magnitude of brain copper influx and efflux are depicted in the size of blue arrows in the bottom row. Note that changes in copper content in brain and whole larvae can be independently modulated by neuronal cell autonomous (C155>) or systemic (actin driver, Act>) modifications of ATP7 gene expression. B, ICP-MS determinations of total larval copper content in UAS-ATP7 or UAS-ATP7 RNAi transgenes expressed from a ubiquitous Actin-GAL4 driver. The mean difference for two comparisons against the shared control Act are shown in the above Cumming estimation plot. The raw data are plotted on the upper axes. On the lower axes, mean differences are plotted as bootstrap sampling distributions. Each mean difference is depicted as a dot. Each 95% confidence interval is indicated by the ends of the vertical error bars. The p value of the two-sided permutation t test is indicated. The number of larvae analyzed is depicted by italic figures in parentheses.
To determine whether systemic ATP7 overexpression affects synaptic structure, we analyzed neuromuscular junction morphology and found that overexpression of ATP7 had no effect on synapse size (Figs. 2A2, 3A,B). In contrast, flies heterozygous for COG1 had increased synapse size measured as the total length of all branches (Fig. 3A,B, COG1e02840/+), a phenotype rescued by overexpression of ATP7 (Fig. 3A,B). This result suggests that the effects of reduced COG1 expression are due in part to the COG effects on endogenous ATP7 expression (Fig. 3A,B). If loss of COG1 results in reduced ATP7 expression, then phenotypes caused by neuronal cell autonomous ATP7 overexpression should be reverted by either systemic or neuronal-selective downregulation of the COG1 subunit (compare Fig. 2A4 and A5). We find that neuronal-specific overexpression of ATP7 using C155-GAL4 results in decreased synapse size (Fig. 3C,D), and this phenotype was rescued by either COG1 heterozygosity or neuronal-specific knock-down of COG1 (Fig. 3C,D). Second, we manipulated the expression of CTR1 and the COG complex cell autonomously in larval motor neurons in which ATP7 expression is downregulated (Figs. 2A7, 3E–G). Neuron-specific downregulation of ATP7 did not change the neuromuscular junction size (Fig. 3E,F) but increased the size of synaptic boutons (Figs. 2A6, 3E,G, insert). This phenotype was rescued by either neuronal-specific COG1 or CTR1 knock-down, indicating that COG1 downregulation modulates the entry of copper from the plasma membrane phenocopying CTR1 downregulation (Figs. 2A7, 3E,G). These results demonstrate that the COG complex regulates cellular copper homeostasis by controlling the function of ATP7 and CTR1 in developing motor neurons.
Copper homeostasis by the COG complex and copper transporters modulate glutamatergic synapse morphology. Muscle 6–7 neuromuscular junctions from wild-type animals and animals expressing ATP7, COG, and CTR1 UAS-transgenes stained with anti-HRP antibodies (A, C, E). C, Synapses from animals carrying the UAS-ATP7 transgene expressed under the systemic driver Actin-GAL4, the COG1e02840/+ genotype, each in isolation or their combination. Scale bar: 50 µm. B, Total branch length quantifications of animals in A. C, Synapses from animals carrying the UAS-ATP7 transgene, UAS-COG1 RNAi expressed under the C155-GAL4 driver, or the COG1e02840/+ genotype, each in isolation or their combinations. D, Quantifications of animals in C. E, Synapses from animals carrying the UAS-ATP7 RNAi, UAS-CTR1 RNAi expressed with the C155-GAL4 driver, or their combinations. F, G, Total branch length and bouton size area quantifications, respectively, of animals in E. Number of animals per genotype (n) are presented by the number at the base each column. One-way ANOVA followed by two-tailed Fisher's least significant difference comparison.
To assess whether ATP7 and COG genetic interactions modulating synapses were just restricted to developing motor neurons, we measured dendritic arborization (da) complexity in Drosophila third-instar larvae class IV (C-IV) da epidermal sensory neurons (Fig. 4). The C-IV da neuron possesses an elaborated and stereotypic dendritic arbor whose complexity is responsive to gene perturbations (Iyer et al., 2013a; Singhania and Grueber, 2014; Das et al., 2017). We overexpressed ATP7 in a neuronal cell-autonomous manner, a model of copper efflux (Fig. 2A4), using the ppk-GAL4 driver (Yang et al., 2009). ATP7 expression significantly decreased the complexity of the dendritic arbor, as determined by the total dendritic length and number of branches (Fig. 4A–C). COG1 haploinsufficiency with the COG1e02840/+ allele did not significantly change dendritic arbors (Fig. 4A–C). However, this COG1 gene dosage reduction was enough for suppressing the ATP7 overexpression dendritic phenotype (Figs. 2A5, 4A–C). The dendritic phenotype caused by ATP7 overexpression and its rescue were more prominent in terminal branches, as indicated by the reversed Strahler analysis, in which distal branches of the 1 and 2 categories were the most affected as compared with branches closer to the cell body (Fig. 4D). These results demonstrate that the COG complex balances the expression of copper influx and efflux transport activities from cells (compare Fig. 2A4, A5). In turn, these changes in copper homeostasis modulate dendritic and presynaptic morphology in diverse developing neuronal types.
The COG complex and ATP7 modulate sensory C-IV da neuron dendrite architecture. A, Representative live confocal images of C-IV da neurons of the specified genotypes labeled by GFP driven by CIV-GAL4 (477;ppk). B–D, Quantitative analysis of total dendritic length and total number of branches in the specified genotypes. OE and OE1 represent different transgenic animal strains overexpressing ATP7 (Norgate et al., 2006). Scale bar: 200 µm. D, Reversed Strahler branch order analysis. Terminal branches correspond to Strahler numbers 1 and 2. Average ± SEM n of animals per genotype are presented by the number at the base each column. For B, C, comparisons were made with one-way ANOVA followed by two-tailed Fisher's least significant difference comparison. D, Analysis was made with two-way ANOVA followed by Tukey's multiple comparisons test.
The COG complex genetically interacts with ATP7 transporters in adult Drosophila tissues
We applied the same experimental logic used for larval synapses to adult Drosophila catecholaminergic and serotoninergic neurons. These cells are only sensitive to ATP7 expression manipulation when the animals are exposed to excess dietary copper (Comstra et al., 2017). We focused on these cells because copper-dependent molecules and transporters are highly expressed in vertebrate catecholaminergic cells (Xiao et al., 2018). As we previously demonstrated, copper feeding control animals increased mortality close to a 100% (Fig. 5A, column 1; Comstra et al., 2017). We used the Ddc-GAL4 driver to overexpress ATP7 in catecholaminergic and serotoninergic neurons. However, this driver is also expressed in larval epidermis (Li et al., 2000; Landgraf et al., 2003). ATP7 overexpression using the Ddc-GAL4 driver decreased copper-induced mortality close to 30% and 60% in adult females and males, respectively (Fig. 5A, columns 1 and 2). We interpret this protective effect as an increased copper efflux because of overexpressed ATP7 (Fig. 2A4). A similar protective outcome was achieved by downregulation of the COG1, COG5, or COG8 subunits (Fig. 5A, column 1 compare to 3, 5, and 7). Here, the protective effect of COG downregulation likely uncovers a decreased CTR1 activity necessary for copper influx into aminergic neurons (Comstra et al., 2017). To define how the COG complex balances the expression of ATP7 and a CTR1 activity, we downregulated COG complex subunits in the ATP7 overexpression background and find that this leads to increased mortality as compared with the COG subunit RNAi alone (Fig. 5A, compare columns 3–4 and 7–8), suggesting that the COG complex preferentially affects the expression of ATP7 over a copper uptake transporter activity (compare Fig. 2A4 and A5). If the COG complex has a stronger influence on ATP7 than CTR1 in aminergic neurons, then knock-down of ATP7, a model of neuronal copper accumulation (Fig. 2A6), should not protect animals after a copper feeding challenge. Additionally, the ATP7 knock-down phenotype should be suppressed by inhibition of the COG complex because of disruption of copper influx mechanisms (Fig. 2A7; Comstra et al., 2017). As predicted, ATP7 RNAi-driven mortality was indistinguishable from wild-type animals (Fig. 5B, compare columns 1 and 2). However, the effect of copper feeding on mortality were suppressed by inhibition of the COG subunits alone or in combination with inhibition of ATP7 (Fig. 5B, compare columns 3–4, 5–6, and 7–8). These results demonstrate that in neurons expressing ATP7, COG complex deficiency preferentially suppresses the effects of ATP7 overexpression (Fig. 2A5). These findings demonstrate the COG complex balances cellular copper homeostasis by modulating the expression of ATP7 and copper influx activity, presumably through CTR1.
Interactions between ATP7 and the COG complex are conserved in Drosophila adult dopaminergic neurons. Male and female adult animals expressing UAS-ATP7A (A), UAS-ATP7 RNAi (B), or RNAi against COG1, COG5, or COG8 alone or in combination in dopaminergic neurons using the ddc-GAL4 driver (A, B). Mortality was induced by copper feeding for 2–3 d. For experiments assessing UAS-ATP7 in both sexes, we used 10 experiments with 86–96 animals in total per genotype and sex, except for UAS-COG8 RNAi experiments were n = 7–8 with 46–77 animals used in total. For experiments assessing UAS-ATP7 RNAi in both sexes, we used 12–38 experiments with 80–336 animals in total per genotype and sex; p values calculated with Kruskal–Wallis test followed by pairwise comparisons with two-tailed Mann–Whitney rank-sum test. See Extended Data Figure 5-1.
Extended Data Figure 5-1
Interactions between ATP7 and the COG complex are conserved in Drosophila adult epidermis. A, Phenotype of transgenic animals expressing UAS-ATP7 or COG subunit RNAi, alone or their combinations, under the control of pnr-GAL4. B, Mortality of adult transgenic animals expressing UAS-ATP7, UAS-COG1, UAS-COG5, or UAS-COG8 RNAi alone or their combinations using the epidermal pnr-GAL4 driver. N = 498 total number of animals used in these experiments; p values and z-scores were estimated as two-tailed differences between two proportions. C, Model of the preferential genetic interactions between the COG complex and copper transporters in epidermis or neurons. Arrow size indicates the relative effect of the COG complex in the activity of the copper transporter. Download Figure 5-1, TIF file.
We validated these genetic interactions and their interpretations in the Drosophila adult epidermis where the interplay of ATP7 and CTR1 in controlling copper homeostasis has been precisely established. In this tissue, overexpression of ATP7 with the epidermal pnr-GAL4 driver depletes cellular copper, a phenomenon rescued by increased expression of the plasma membrane copper uptake transporter CTR1 (Norgate et al., 2006; Burke et al., 2008; Binks et al., 2010; Hwang et al., 2014; Mercer et al., 2016). We reasoned that, in contrast with neurons, phenotypes caused by epidermal overexpressed ATP7 would be additive with epidermal COG complex subunit RNAi because of downregulation of a postulated CTR1 copper influx activity. Thus, supporting the concept that copper homeostasis is chiefly controlled through a putative CTR1 activity in the epidermis (Extended Data Fig. 5-1). Overexpression of ATP7 caused incomplete dorsal closure and depigmentation of the thorax and abdomen of adult female and males flies (Extended Data Fig. 5-1A). These phenotypes were maladaptive, as demonstrated by decreased animal survival to ∼38% and ∼65% in adult females and males, respectively (Extended Data Fig. 5-1B, column 2). These ATP7-dependent epidermal phenotypes were absent in animals where the expression of the COG complex was reduced by RNAi against any one of three COG subunits (Extended Data Fig. 5-1B, column 3). However, simultaneous expression of ATP7 in animals where the COG complex was downregulated by RNAi significantly decreased survival to ∼22% and ∼35% in adult females and males, respectively (Extended Data Fig. 5-1B, column 4). This doubling in mortality represents an additive phenotype. These results suggest that the COG complex modulates copper homeostasis by preferentially controlling the expression of a copper influx activity, possibly CTR1, in epidermal cells (Extended Data Fig. 5-1C). These epidermal findings further validate our genetic strategy to examine whether an intact COG complex is required for balancing copper influx and efflux in neurons.
Mitochondrial respiration requires a COG-dependent and copper-dependent mechanisms
A possible target of copper dyshomeostasis in motor and sensory Drosophila neurons is the mitochondria as both genetic defects in ATP7A and ATP7B perturb mitochondrial morphology and function (Yoshimura and Kudo, 1983; Yamano and Suzuki, 1985; Onaga et al., 1987; Yamano et al., 1988; Gu et al., 2000; Zischka et al., 2011; Bhattacharjee et al., 2016; Lichtmannegger et al., 2016; Guthrie et al., 2020). Similarly, COG complex mutants decrease tetrazolium salt metabolization mediated by the NAD(P)H-dependent oxidoreductase and dehydrogenase enzymes, which localize both to the cytoplasm and mitochondria (Berridge et al., 2005; Comstra et al., 2017). This suggests that the COG complex could modulate mitochondrial function in a copper-dependent manner. We tested this hypothesis using a combination of phylogenetic-bioinformatic analyses, measuring mitochondrial respiration with Seahorse oximetry, and biochemical approaches.
We used evolutionary rate covariation to determine functional associations between either the ATP7A and ATP7B proteins or the COG complex with subunits of each one of the five mitochondrial respiratory chain complexes. This approach allows us to examine coevolutionary forces that act on genes even when their products do not physically interact or may reside in different compartments (Fig. 6A,B, ERC; Clark et al., 2012, 2013; Findlay et al., 2014). The phylogenic correlation between ATP7A and ATP7B was significant with mitochondrial Complexes I, III, and IV, the latter, a copper-containing respiratory chain complex (Fig. 6A; Horn and Barrientos, 2008). Moreover, the eight subunits of the COG complex also covaried significantly with Complex I (Fig. 6A). We analyzed the evolutionary rate covariation among all pairs of a gene group integrated by ATP7A, ATP7B, the COG complex, and the mitochondrial Complex I subunits. Gene group analysis of the evolutionary rate covariation among these 48 genes indicated that this gene set is significantly correlated phylogenetically (Fig. 6B, p = 0.00130). Also, pairwise analyses showed that ATP7A coevolved with ATP7B (evolutionary rate covariation scores >0.3, p < 0.05), and multiple Complex I subunits (Fig. 6B,C). Conversely, COG complex subunits significantly coevolved with ATP7B, and Complex I subunits (Fig. 6B,C). These results support a model where ATP7A, ATP7B, and the COG complex participate in a pathway shared with mitochondria.
The COG complex activity is required for mitochondrial respiration. A–C, Evolutionary rate covariation (ERC) between gene groups constituted by human ATP7A and ATP7B or the COG complex subunits with subunits of the five mitochondrial respiratory complexes defined by the CORUM database. A, Group probability of evolutionary covariation. B, ERC for pairs of genes in a gene group that includes ATP7A and ATP7B, the eight COG complex subunits, and the 38 nuclear encoded subunits of the mitochondrial respiratory Complex I (CORUM complex ID:388). ERC threshold >0.3 identifies covariated gene pairs. C, p values for genes pairs presented in B; α <0.05 Wilcoxon rank-sum test against 100,000 permutations. ERC and p values are represented by circle color and size. D–H, COG null cells have defective respiration that can be ameliorated with a copper ionophore. Extracellular Seahorse flow oxymetry was used to measure oxygen consumption rate (OCR) in wild-type, COG1Δ/Δ, and COG8Δ/Δ HEK293 cells with or without the addition of 100 nm disulfiram in the presence of 180 nm copper. (D) OCR each time points over 300 min. E, Average basal OCR. F, Acute response. Mitochondrial respiration was isolated by the addition of 0.1 μm oligomycin (G), and non-mitochondrial respiration was measured by the addition of 0.5 μm rotenone and antimycin A (H). N = 4, non-parametric Kruskal–Wallis test followed by pairwise Mann–Whitney U test comparisons. See Extended Data Figure 6-1.
Extended Data Figure 6-1
Biochemical and morphological mitochondrial phenotypes in COG-deficient cells. A, Immunoblots of wild-type, COG1, and COG8 mutant cell extracts probed with antibodies against mitochondrial markers listed. HSP60 and SDHA immunoblots were used as loading controls. B, Confocal immunomicroscopy of wild type and COG1 stained with antibodies against the mitochondrial antigen TOMM20. 3D rendering was performed with Imaris. Download Figure 6-1, TIF file.
We tested whether COG genetic defects could impair mitochondrial function by measuring oxygen consumption rates with SeahorseTM technology (Divakaruni et al., 2014). We used HEK293 cells where both genomic copies of the COG subunits 1 (COG1Δ/Δ) or 8 (COG8Δ/Δ) genes were edited by CRISPR-Cas9 (Bailey Blackburn et al., 2016; Blackburn and Lupashin, 2016). These COG null cells phenocopy established COG complex deficiency phenotypes, which include defective Golgi-dependent glycosylation of membrane proteins and their degradation in lysosomal compartments (Shestakova et al., 2006; Bailey Blackburn et al., 2016; Blackburn and Lupashin, 2016). COG null HEK293 cells had a lower basal oxygen consumption rate compared with wild-type controls (Fig. 6D,E). To determine whether this reduced basal respiration is because of impaired copper availability, we treated cells with disulfiram, a drug with several targets and also a copper ionophore (Allensworth et al., 2015; https://go.drugbank.com/drugs/DB00822). We find that the basal respiration rate in wild-type cells in the presence of 180 nm copper was not changed even after 4 h (Fig. 6D,F, acute response). In contrast, the addition of copper-disulfiram to COG1Δ/Δ and COG8Δ/Δ HEK293 cells significantly improved the basal respiration rate, a response that was evident at 2 h and progressively increased up to 4 h of incubation with copper-ionophore (Fig. 6D,F, acute response). This selective improvement in oxygen consumption by COG null cells was sensitive to the mitochondrial ATPase inhibitor oligomycin (Divakaruni et al., 2014), thus demonstrating the mitochondrial origin of the copper-disulfiram rescue of the respiration phenotype in COG-null HEK293 cells (Fig. 6D,G). Non-mitochondrial sources of oxygen consumption were not affected by genotype or copper ionophore addition, as determined by measuring oxygen consumption in the presence of antimycin and rotenone, inhibitors of cytochrome C reductase and Complex I, respectively (Fig. 6D,H; Divakaruni et al., 2014). These results demonstrate that mitochondrial respiration defects observed in COG-deficient cells can be acutely corrected by circumventing impaired copper transport with a copper ionophore.
We next asked whether these mitochondrial functional defects in COG mutant cells correlated with changes in mitochondria composition and/or morphology. A proteome of COG mutant cell membranes identified ∼13% of mitochondria annotated proteins among all significant COG-sensitive hits (V. Lupashin's unpublished data). We used this dataset to inform a selection of mitochondria-localized proteins to assess their expression in COG1 and COG8 mutant cells. These selected mitochondrial proteins include copper-delivery factors to Complex IV in the electron transport chain (COX11 and COX17; Baker et al., 2017), mitofusin (MFN1; Detmer and Chan, 2007), the pyruvate carrier MPC1 (Bricker et al., 2012), and a mitochondrial protein involved in cholesterol metabolism (DHCR24). These five proteins were annotated to mitochondria by Mitocarta 2.0 (Calvo et al., 2016). Moreover, these five proteins increased their content in COG1 and COG8 null cells as compared with wild-type cells (Extended Data Fig. 6-1A). Despite the Seahorse and biochemical mitochondrial phenotypes in COG null cells, we could not detect alterations in the morphology of mitochondria as revealed by immunofluorescence microscopy of wild-type and COG1 null cells (Extended Data Fig. 6-1B). Thus, phylogenetic, functional, and biochemical results show that COG complex-dependent copper homeostasis is necessary for mitochondrial integrity.
ATP7-dependent and COG-dependent copper homeostasis sustain mitochondrial synaptic content
We sought to answer whether ATP7-dependent and COG-dependent copper homeostasis was required for normal Drosophila mitochondria at motor neuron synapses. To address this question, we generated animals carrying both UAS-ATP7 and a UAS-mitochondria-targeted GFP reporter (Fig. 7, MitoGFP), under the control of the C155-Gal-4 neuronal driver. We focused on the neuronal ATP7 overexpression/copper efflux paradigm because ATP7 overexpression synaptic phenotypes are robustly rescued just by a single gene copy loss of COG1 (Figs. 2A4,A5, 3C,D). We measured synapse morphology on neuromuscular junctions double-stained with antibodies against HRP and GFP to label the neuronal plasma membrane and mitochondria, respectively. Mitochondria were present in neuromuscular junctions, including muscles 6 and 7, a standard synaptic model (Fig. 7A). A single copy loss of the COG1 gene is not sufficient to alter mitochondrial content and distribution in neurons (Fig. 7A–D) even when the COG1 haploinsufficiency increased synapse branching (COG1e02840/+; Fig. 7A,B). In contrast, ATP7 overexpression in neurons resulted in a nearly complete depletion of mitochondria from synapses (Figs. 2A4, 7A,C,D) with a concomitant increase in mitochondrial signal at ventral nerve cords, a site where motor neuron cell bodies reside (Fig. 7E). This pronounced mitochondrial synapse depletion phenotype was rescued in animals where the UAS-ATP7 transgene is expressed in the COG1e02840/+ background (compare Figs. 2A4 and 2A5, 7A,C,D). The COG1e02840/+-dependent rescue of the ATP7 overexpression mitochondrial phenotype also resulted in a normalization of the mitochondrial signal in the ventral nerve cord (Fig. 7E).
The COG complex and ATP7 are required to maintain mitochondrial content in glutamatergic synapses. A, B, Muscle 6–7 neuromuscular junctions from third-instar larvae from wild-type animals and animals expressing ATP7 transgene, carrying the COG1e02840/+ allele or the combination were crossed to the mitochondrial UAS-mitochondrial-GFP reporter. Reporter genes were expressed with the neuronal C155-GAL4 driver. Dissected larvae were stained with anti-HRP antibodies and anti-GFP. Scale bar: 50 µm. B–D, Quantification of synaptic branches and mitochondria content expressed as HRP or GFP area and their ratios. Average ± SEM n of animals per genotype are presented by the number at the base each column. One-way ANOVA followed by two-tailed Fisher's least significant difference comparison. E, Third-instar larvae from wild-type animals and animals expressing ATP7 transgene, carrying the COG1e02840/+ allele, or their combination were crossed to the mitochondrial UAS-GFP reporter. Transgenes were expressed with the neuronal C155-GAL4 driver. Dissected ventral cords were stained with anti-GFP and imaged by confocal microscopy. Scale bar: 25 µm.
In parallel, we conducted analyses in Drosophila third-instar larvae C-IV da neurons by overexpressing ATP7A, the plasma membrane marker CD4-tdTomato, and MitoGFP in a cell-autonomous manner using the 477-GAL4 and ppk-GAL4 drivers (Fig. 8). ATP7 overexpression decreased the mitochondrial content in dendrites (Figs. 2A4, 8A–C). This phenotype was partially rescued by the COG1e02840/+ allele (Figs. 2A5, 8A), as indicated by a restoration of the mitochondria numbers at branch points of the dendritic tree (Fig. 8C). These results demonstrate that ATP7, COG, and their interactions modulate dendrite and presynaptic mitochondrial content in motor and sensory neurons. These ATP7 and COG-dependent synaptic branching and mitochondrial phenotypes predict intricate alterations in synaptic transmission and plasticity in neurons with copper dyshomeostasis.
The COG complex and ATP7 are required to maintain mitochondria in sensory C-IV da neuron dendrites. A, Representative live confocal images of C-IV da neurons of the specified genotypes expressing a mitochondria-targeted GFP UAS-transgene and a UAS-plasma membrane marker (CD4-td-Tomato) under the control of the 477-GAL4 driver. B, C, Quantification of mitochondria content as organelles per branch length unit or as organelles at branching points. Scale bar: 100 µm. Average ± SEM n of animals per genotype are presented by the number at the base each column. One-way ANOVA followed by two-tailed Fisher's least significant difference comparison.
The COG complex and ATP7 modulate basal release probability
We analyzed the effects of single or combined changes in the expression of ATP7 and the COG complex with the COG1e02840/+ allele on neurotransmission (Figs. 2A4, A5, 9–10). We measured synaptic transmission at the neuromuscular junction of Drosophila larvae synapses on muscles 6–7 under two-electrode voltage-clamp. We recorded nerve evoked junctional currents (nEJCs) in the different genotypes at low extracellular Ca2+, a condition where basal changes in release probability are easily distinguishable (Acharya et al., 2006; Jorquera et al., 2012; Stevens et al., 2012). Control strains carrying the neuronal GAL4 driver (elav-C155) or the UAS-ATP7 transgene were similar in their nEJCs amplitude (Fig. 9A). However, COG1e02840/+ haploinsufficiency and ATP7 neuronal overexpression/copper efflux displayed increased nEJC amplitudes (Fig. 9A). The nEJCs were significantly larger in COG1 haploinsufficiency and ATP7 neuronal overexpression by 3-fold and 4-fold, respectively (Fig. 9A,B). This nEJC amplitude phenotype was rescued in COG1e02840/+ animals overexpressing ATP7 in neurons (Fig. 9A,B). The phenotype of large amplitude in nEJCs in COG1e02840/+ animals suggests an increased release probability, as synaptic size and bouton number directly correlate with the amplitude of the overall synaptic responses in healthy synapses (Zhang et al., 2001; Menon et al., 2013). However, in pathologic conditions, diverse homeostatic mechanisms of compensation can modify the release probability (Wondolowski and Dickman, 2013; Davis and Müller, 2015; Frank et al., 2020), evident in ATP7A overexpressing neurons, where the amplitude inversely correlates with the synaptic growth phenotype (Figs. 3C,D, 9A,B).
The COG complex and ATP7 modifies nerve-evoked synaptic transmission and fast plasticity. Representative nEJCs in each genotype at 0.2 mm extracellular Ca2+ (A) and the average amplitude (B). Overlapping of the average normalized nEJCs in each genotype magnified in a new time window (C) and the skew factor of kurtosis in each genotype (D). Facilitation in each genotype induced by two pulses spaced at 50 ms (E) and the average (F). Representative recording of nEJCs during short-term tetanic facilitation induced with nerve stimulation at 20 Hz (G) and the average tetanic facilitation (H); n of animals per genotype are presented by the number at the base each column. One-way ANOVA followed by Tukey's HSD post hoc test.
A synaptic response can be prolonged with increases in release probability, therefore we looked at the nEJCs kinetics. Normalizing average responses during nerve stimulation at low-frequency (Fig. 9C) and estimating a skew values index of kurtosis dividing the charge transferred by their peak amplitude (Fig. 9D), show the slower timing in COG1-haploinsufficiency and ATP7 neuronal overexpression by 2.5 and 2 ms, respectively. This phenotype is rescued partially in animals carrying COG1e02840/+ and expressing the UAS-ATP7 transgene in neurons, indicating that the ATP7 levels modulate amplitude and kinetics of synaptic responses. The neurotransmission phenotype in COG1e02840/+ and neuronal expressed ATP7 animals agree with increased release probability, we further tested this hypothesis measuring facilitation using paired-pulse and tetanic stimulation. To determine the facilitation index in these genotypes, we applied paired-pulse stimulation spaced by 50 ms (Fig. 9E,F). Control strains facilitate, COG1e02840/+, and ATP7 neuronal overexpression facilitate similarly (Fig. 9E,F). However, tetanic facilitation after a brief train of stimuli at 20 Hz is decreased significantly in animals overexpressing ATP7 in neurons (Fig. 9G,H), indicating a quick saturation in the release probability during nerve activity. These phenotypes were rescued by COG1 haploinsufficiency with the COG1e02840/+ allele (Fig. 9G,H). Altogether, our results demonstrate that ATP7 levels in motor neurons increase the overall initial release probability.
The COG complex and ATP7 modify tetanic and post-tetanic short-term synaptic plasticity
Neuronal transmission undergoes short-term changes in synaptic efficacy during high-frequency nerve activation. After the conditioning episode ceases, these changes in synaptic efficacy persist and are progressively restored to their preconditioned values over time, a phenomenon known as short-term synaptic memory. Short-term synaptic memory is key for learning and working memory, as determined in Drosophila mutants with altered Pavlovian memory retention where neuromuscular synaptic modifications correlate with memory retention phenotypes (Zhong and Wu, 1991). To scrutinize the effects of tetanic and post-tetanic plasticity, we analyzed the nEJCs during a protocol that induces PTP under low release probability conditions (Fig. 10). In all genotypes, test stimulation for 40 s at 0.5 Hz evoked stable nEJCs (Fig. 10A,B). However, the tetanic episode, induced by 25 s of nerve stimulation at 20 Hz, dramatically enhance the synaptic transmission decaying in control to the pretetanic average values within 50 s of test stimulation after the conditioning episode ceased (Fig. 10A,B). However, COG1e02840/+ and neuronal overexpressed ATP7 synapses cannot sustain this increase (Fig. 10A,B). Besides, synaptic transmission decayed after a tetanic episode below the test control in both genotypes, suggesting the presence of a tetanic depression and PTD. The phenotype observed in synapses overexpressing ATP7 was rescued by the COG1e02840/+ allele (Fig. 10A,B) indicating that the genetic interaction between ATP7 and the COG complex is necessary to sustain short-term plasticity.
The COG complex and ATP7 alter short-term synaptic-plasticity in conditioning and postconditioning episodes. A, Average amplitude of nEJCs during the application of a short-term synaptic memory protocol in each genotype, and representative electrophysiological recordings of synaptic transmission (inset). The average instantaneous skewness in each condition is displayed at the top. Test and conditioning stimulation were conducted at 0.5 and 20 Hz, respectively. B, Average responses during the low-frequency test stimulation (0.5 Hz) before and after the tetanic episode normalized by the pretetanic value, and the average increase in synaptic efficacy during tetanic activity in each genotype. C, The average estimation for tetanic potentiation, PTP, and depression, with the time constant of transition. All graphs in C have the same n of animals per genotype. These are presented by the number at the base the tau columns. One-way ANOVA followed by Tukey's HSD post hoc test.
To evaluate the presynaptic sources that could modify neurotransmission in these genotypes, we tested the capacity to change their synaptic efficacy by looking in detail at the tetanic and post-tetanic plasticity parameters (Fig. 10C). Normalized tetanic responses indicate that all genotypes display different profiles of enhancement after the initial tetanic facilitation phase (Fig. 10B). Nerve activity induces the characteristic phase of augmentation within 10 s and it reaches the maximal potentiation in synaptic enhancement in control strains. Synapses from animals carrying the COG1e02840/+ allele took less time during the augmentation phase, also reaching lower potentiation levels (Fig. 10B,C). In turn, ATP7 neuronal overexpression rapidly progressed to depression (Fig. 10A–C). These data are consistent with an increase in release probability and abnormal synaptic vesicle availability during periods of high demand. Interestingly, this synaptic behavior resembles the phenotype observed in Drosophila drp1, mutants devoid of mitochondria in presynapses (Verstreken et al., 2005). However, the synaptic transmission behavior occurring during postconditioning episodes, such as short-term synaptic memory, is undetermined in this kind of synapses (Verstreken et al., 2005). Pretetanic and post-tetanic response analysis indicate that control strains induced the characteristic PTP decaying to pretetanic values with a time constant of 45 s (Fig. 10A–C). However, the COG1e02840/+ allele displayed less PTP which progressed to PTD (Fig. 10A–C). This PTP phenotype was even more pronounced in ATP7 overexpressing neurons, which supports the concept of abnormal synaptic vesicle availability during high demand activity. Importantly, these PTP phenotypes observed in either COG1e02840/+ or ATP7 overexpressing animals were rescued in animals simultaneously carrying the COG1e02840/+ allele and expressing ATP7 in neurons (Fig. 10A–C). Our data demonstrate that genetic interactions between ATP7 and the COG complex, that maintain neuronal copper homeostasis, are necessary for synaptic function, including basal probability of vesicle release, synaptic vesicle availability, and short-term plasticity during conditioning and postconditioning episodes.
Discussion
We demonstrate phylogenetically conserved genetic interactions between ATP7 paralogs and a Golgi complex tether, the COG complex. Disruption of these interactions alters presynaptic function, mitochondrial content, and plasticity. Drosophila phenotypes caused by changes in the expression of ATP7 are modulated by downregulation of COG complex subunits in diverse neurons. These cell types are affected in Menkes disease by mechanisms yet to be fully explored but postulated to include defective activity of cuproenzymes in mitochondria (Zlatic et al., 2015; Hartwig et al., 2019; Guthrie et al., 2020). We conclude that genetic and biochemical interactions between ATP7 and COG are necessary to maintain cellular copper homeostasis, which in turn is necessary for mitochondria functional integrity and localization.
The steady state expression of ATP7A, CTR1, and COG complex subunits requires the mutual integrity of these genes to maintain their balanced expression. The ablation of the COG complex downregulates these copper transporters in mammalian cells in a cell-autonomous manner (Comstra et al., 2017). This observation led us to design a set of genetic tests of copper homeostasis predictions in Drosophila tissues (Fig. 2A). We scrutinized these copper homeostasis predictions using both systemic and cell-autonomous changes in the expression of ATP7, CTR1, and COG complex subunits (Fig. 2A). We used ATP7 neuronal cell-autonomous overexpression to induce cytoplasmic copper efflux from cells, thus mimicking one of the Menkes-like states which are characterized by both reduced cytoplasmic and Golgi copper (Fig. 2A4). Conversely, ATP7 RNAi was used to increase cellular copper (Fig. 2A6). These are established paradigms to study copper homeostasis in Drosophila (Norgate et al., 2006; Burke et al., 2008; Hwang et al., 2014; Mercer et al., 2016). We conclude that ATP7-dependent and CTR1-dependent copper homeostasis is under control of the COG complex in neurons. The following genetic evidence supports this assertion: (1) the rescue of diverse developmental neuronal cell-autonomous ATP7 overexpression/copper efflux phenotypes in neuronal tissues by COG gene defects (Figs. 3C, 4A, 6–10); (2) the COG-dependent rescue of cell-autonomous ATP7 RNAi/copper accumulation phenotypes in neurons which we propose represents a reduced CTR transport activity at the plasma membrane (Figs. 3E, 5B); (3) the rescue of copper induced mortality phenotype by cell-autonomous downregulation of three COG complex subunits in aminergic neurons (Fig. 5B); (4) the rescue of a synaptic haploinsufficiency COG phenotype by systemic overexpression of ATP7 (Fig. 2A); (5) the similar rescue of cell-autonomous ATP7 RNAi/copper accumulation synaptic phenotypes granted by either neuronal-specific COG subunit or CTR1 downregulation (Fig. 2E); and (6) finally, the rescue of mitochondrial respiration defects by copper ionophores in COG-null human cells (Fig. 6).
While our genetic and biochemical analyses point to a contribution of copper in COG deficient phenotypes in flies and mammalian cells, it is reasonable to expect the contribution of other receptors and membrane transporter to neuronal phenotypes described here. COG genetic defects characteristically impair the activity of glycosyltransferases localized to the Golgi complex (Ungar et al., 2002; Zolov and Lupashin, 2005; Climer et al., 2015; Bailey Blackburn et al., 2016; Blackburn et al., 2018; D'Souza et al., 2019). These enzymes are no longer retained in the Golgi apparatus of COG-deficient cells, leading to a global decrease in the glycosylation of proteins traversing the secretory pathway. This defective glycosylation destabilizes membrane proteins residing in the endomembrane system in mammals and Drosophila (Ungar et al., 2002; Zolov and Lupashin, 2005; Shestakova et al., 2006; Pokrovskaya et al., 2011; Climer et al., 2015; Bailey Blackburn et al., 2016; Frappaolo et al., 2017). These COG-dependent biochemical defects would determine phenotypes in human COG mutations, also known as congenital disorders of glycosylation Type II. However, which neuronal proteins undergo defective glycosylation and destabilizations in this disorder is still unclear (Barone et al., 2014; Climer et al., 2015, 2018). We speculate that neurotrophic factors and receptors, synaptic adhesion molecules, calcium and ion channels, and/or synaptic membrane proteins may contribute to the synaptic phenotypes observed in COG deficient Drosophila neurons (Kadmon et al., 1990; Han et al., 2004; Kleene and Schachner, 2004; Dong et al., 2008; Kwon and Chapman, 2012; Weiss et al., 2013; Lazniewska and Weiss, 2017).
A prominent phenotype observed in neurons where copper efflux is increased by ATP7 overexpression (Fig. 2A4), and a focus of our studies, is a depletion of mitochondria at neuromuscular synapses and dendrites. This mitochondrial phenotype is concomitant with an accumulation of mitochondria in the ventral nerve cord where motor neuron cell bodies reside (Fig. 7A,E). These Drosophila mitochondrial phenotypes are reminiscent of mitochondrial defects in Menkes disease (Yoshimura and Kudo, 1983; Onaga et al., 1987; Yamano et al., 1988; Zlatic et al., 2015). We attribute the mitochondrial phenotype in ATP7 overexpressing neurons to defective copper loading in mitochondria, yet this idea needs a direct experimental test as the one performed in a mouse Menkes model (Guthrie et al., 2020). Our observation that COG-null human cells are defective in cellular copper and mitochondrial respiration, the later rescued by copper ionophores, support a model of copper depletion-dependent damage of mitochondria (Fig. 6). Irrespective of whether mitochondria are impaired solely by copper-dependent mechanisms or other additional mechanisms that require copper transporter activity, such as oxidative stress; we postulate that functional mitochondria defects may either prevent their delivery to nerve terminals by impairing mitochondria attachment to kinesin motors, prevent their anchoring to nerve terminals, and/or induce their engulfment by mitophagy, thus trapping them in cell bodies. Copper-dependent damage of mitochondria alters their morphology, uncouples the respiratory chain depolarizing mitochondria, and induces mitophagy (Yoshimura and Kudo, 1983; Yamano and Suzuki, 1985; Onaga et al., 1987; Yamano et al., 1988; Gu et al., 2000; Zischka et al., 2011; Bhattacharjee et al., 2016; Lichtmannegger et al., 2016; Polishchuk et al., 2019). Mitochondrial uncoupling and depolarization change mitochondria kinesin-dependent movement and induce mitophagy in cell bodies (Overly et al., 1996; Miller and Sheetz, 2004; Roberts et al., 2008; Verburg and Hollenbeck, 2008; Wang et al., 2011; Cai et al., 2012; Sheng and Cai, 2012; Ashrafi et al., 2014; Mishra and Chan, 2016; Misgeld and Schwarz, 2017). This evidence suggests that these mitochondrial distribution mechanisms may be perturbed or engaged in neurons whose copper homeostasis is disrupted by defects in ATP7 and/or the COG complex.
Synapse morphology, mitochondrial content, evoked neurotransmission, and short-term plasticity are modified by changes in the expression of ATP7 and COG complex subunits (Figs. 3–4, 7–10). How these phenotypes are intertwined into cause and effect relationships remains to be explored. We speculate, however, that the absence of mitochondria may account in part for the modified neurotransmission observed in ATP7 overexpression/copper efflux synapses. Mitochondria could contribute to synapses: their calcium buffering capacity (Vos et al., 2010), their ATP synthase activity (Rangaraju et al., 2014), their Krebs cycle metabolites capable of modulating neurotransmission (Ugur et al., 2017), or enzymatic activities necessary for glutamate metabolism (Kovacevic and McGivan, 1983). In addition, mitochondria could offer metabolic flexibility diversifying sources of carbon entering the Krebs cycle for energy production at the synapse (pyruvate, glutamate, or fatty acids; Olson et al., 2016). As an example, synaptic calcium regulation by mitochondria may explain some of the neurotransmission phenotypes in ATP7 overexpression/copper efflux brains. Short episodes of high-frequency nerve-stimulation induce plasticity in the form of facilitation, augmentation, and potentiation. These are all short-term increases in the efficacy of synaptic transmission, also known as synaptic enhancement (Zucker and Regehr, 2002). Facilitation depends on the residual Ca2+ after influx; in turn, augmentation relies on the management of the resting Ca2+ during nerve activity, which promotes the exocytosis of the readily releasable pool of vesicles, whereas potentiation depends mostly on the exo/endocytic balance (Zucker et al., 1991; Kamiya and Zucker, 1994; Regehr et al., 1994). Post-tetanic facilitation and augmentation phenomena are related to the decay of residual calcium (Edelman et al., 1987; Zucker et al., 1991; Zucker and Regehr, 2002). However, the resting calcium explains the PTP, indicating that during activity the rising calcium drives phosphorylation-dependent signaling activation, which lasts longer than the calcium transient itself (Macleod et al., 2002; Lnenicka et al., 2006). In that sense, mitochondria could be recruited during the development of PTP (Tang and Zucker, 1997; Zhong et al., 2001). Thus, mitochondria could maintain or modulate calcium levels at the nerve terminal to sustain diverse forms of synaptic plasticity by buffering the tetanic activity critical for postconditioning episodes like short-term synaptic memory. Our data show that tetanic and post-tetanic synaptic enhancements are altered in ATP7 overexpressing/copper efflux synapses, a phenotype that can be reverted by COG complex downregulation. We do not have direct evidence that mitochondria are solely responsible for the synaptic phenotypes in ATP7 overexpression. However, the parallel rescue of both mitochondrial content and neurotransmission phenotypes in ATP7 overexpression synapses after COG complex downregulation strongly support the hypothesis of a mitochondrial-dependent mechanisms account for the neurotransmission phenotypes in these terminals. We propose that copper dyshomeostasis synaptic phenotypes in the fly provide a foundation to understand behavioral and neurologic phenotypes characteristic of Menkes and Wilson disease.
Footnotes
This work was supported by National Institutes of Health (NIH) Grants 1RF1AG060285 (to V.F.), R01NS108778 and R01NS108778-01S1 (to R.A.J.), R15AR070505 (to A.D.V.-M.), Telethon TIGEM-CBDM9 (to R.S.P.), R01NS086082 (to D.N.C.), R01GM083144 (to V.L.), and 5K12GM000680-19 (to C.H.). Stocks obtained from the Bloomington Drosophila Stock Center (NIH Grant P40OD018537) were used in this study.
The authors declare no competing financial interests.
Acknowledgements: We thank the Faundez lab members for their comments.
- Correspondence should be addressed to Ramon A. Jorquera at ramon.jorquera{at}unab.cl or Victor Faundez at vfaunde{at}emory.edu
