Summary
Parkinson's disease has been linked to altered mitochondrial function. Mutations in parkin (park), the Drosophila ortholog of a human gene responsible for many familial cases of Parkinson's disease, shorten lifespan, abolish fertility, and disrupt mitochondrial structure. However, the role played by Park in mitochondrial function remains unclear. Here we describe a novel Drosophila gene, clueless (clu) that encodes a highly conserved tetratricopeptide repeat protein closely related to the CluA protein of Dictyostelium, Clu1 of Sacchromyces cerevisiae, and to similar proteins in diverse metazoan eukaryotes from Arabidopsis to humans. Like its orthologs, loss of Drosophila clu causes mitochondria to cluster within cells. We find that strong clu mutations resemble park mutations in their effects on mitochondrial function, and that the two genes interact genetically. Conversely, mitochondria in park homozygotes become highly clustered. We propose that Clu functions in a novel pathway that positions mitochondria within the cell based on their physiological state. Disruption of the Clu pathway may enhance oxidative damage, alter gene expression, cause mitochondria to cluster at microtubule plus ends, and lead eventually to mitochondrial failure.
Introduction
Parkinson's disease results from the catastrophic loss of dopaminergic neurons within the substantia nigra of the brain. Significant declines in mitochondrial function are now recognized as a frequent feature of both heritable and sporadic Parkinson's disease (reviewed in Abou-Sleiman et al. 2006). Mutations or environmental factors associated with the disease are thought to increase oxidative stress within susceptible neurons by a wide variety of mechanisms, leading eventually to mitochondrial failure, accelerated cell death and the onset of symptoms. Mutations in parkin (Kitada et al. 1998) are found in about 50% of heritable cases of autosomal recessive, juvenile Parkinsonism (Betarbet et al. 2005), and the protein is found within lesions in the affected neurons of patients with sporadic disease (Schlossmacher et al. 2002). Although less frequent, pink1 mutations (Valente et al. 2004) are also responsible for a significant number of cases of familial Parkinson's disease.
Studies of the Drosophila parkin (park) and pink1 orthologs (Greene et al. 2003; Wang et al. 2006; reviewed in Pallanck and Greenamyre, 2006) strongly support the view that mitochondrial defects are central to Parkinson's etiology. Flies mutant for park contain defective mitochondria in many tissues and prematurely lose dopaminergic neurons in the brain. Park is predicted to function as an E3 ubiquitin ligase and might act to remove damaged mitochondrial proteins via ubiquitinylation and degradation in the proteasome. Pink1, a mitochondrially targeted ser/thr kinase, is also required for normal mitochondrial morphology and function, possibly by facilitating mitochondrial fission/fusion (Clark et al. 2006; Wang et al. 2006; Yang et al. 2008). Overexpression of Park can rescue pink1 mutant flies, suggesting that these genes function in a common pathway (Clark et al. 2006; Park et al. 2006). Interactions of pink1 and park with genes controlling mitochondrial dynamics suggest that the mitochondrial fission pathway may be involved (Poole et al. 2008). However, the specific molecular mechanisms connecting Park and Pink1 to mitochondrial function remain imperfectly understood.
The Dictyostelium cluA gene encodes a highly conserved protein that affects mitochondrial localization within cells, but its relationship to other mitochondrial pathways is unclear (Zhu et al. 1997). Mitochondria aggregate into a single cluster in cluA mutant cells, and similar clustering occurs in Baker's yeast mutant for clu1, an orthologous gene that can functionally substitute for cluA (Fields et al. 1998, Dimmer et al. 2002). Metazoan eukaryotes contain CluA orthologs, such as friendly mitochondria (fmt) in Arabidopsis, a gene that is also required to prevent mitochondria from clustering (Logan et al. 2003). The C. elegans gene clu-1, and a human ortholog, KIAA0664 (“human Clu”), exist but little information is currently available on their function.
Drosophila oogenesis (for review see Spradling, 1993) represents a highly favorable system for studying the functional importance of mitochondrial subcellular localization. Developing follicles grow extensively, facilitating the visualization of organelles within the oocyte and its 15 interconnected nurse cells. Moreover, mitochondria undergo a series of developmentally regulated behaviors as germ cells progress from stem cell to completed egg (Cox and Spradling, 2003). These include programmed fission after the stem cell stage, Balbiani body formation at the time of follicle formation, and extensive replication within growing follicles. Microtubule-based transport mediated by Dynein (Dhc), Kinesin (Khc) and the adaptor protein Milton (Milt) play a critical role in these events (Cox and Spradling, 2006). By studying how mitochondria redistribute in ovaries bearing Dhc, Khc and milt mutations, the predominant orientation of the microtubules used to transport mitochondria could be deduced (Cox and Spradling, 2006).
Here we describe the gene “clueless (clu)” encoding a Drosophila CluA ortholog that is 53% identical to human Clu, and show that clu mutations cause mitochondrial dysfunction and clustering. clu null mutant adults display a shortened lifespan, male and female sterility, as well as severe mitochondrial abnormalities in flight muscle, similar to parkin mutants. Furthermore, park mutations cause clustering of mitochondria in follicle cells and nurse cells during oogenesis. clu interacts genetically with park because park/+; clu/+ trans-heterozygotes show enhanced mitochondrial clustering. Based on these findings we propose that Clu and Park participate in a physiological and oxidative damage control pathway(s) that link mitochondrial localization and function.
Results
clu encodes a Drosophila CluA ortholog
The Drosophila gene CG8443 produces a single annotated transcript encoding a 1448 aa predicted protein that is closely related to CluA from Dictyostelium and its orthologs (Fig. 1A, B). Consequently, we named this gene clueless (clu). Clu shows 53% identity to human Clu (KIAA0664) across the entire protein. Both proteins contain a region of tetratricopeptide repeats (TPR), as well as a putative “Clu domain” between amino acids 424-666 of Drosophila Clu, that is even more strongly conserved in metazoan Clu proteins (85% identity between Clu and KIAA0664, Fig. S1). N- and C-termini protein segments of Drosophila Clu were expressed and used to raise antibodies in guinea pigs (Fig. 1B, Methods). Both antibodies recognize a single band of apparent molecular weight 160 kd on Western blots (Fig. 1C) whose abundance is reduced in clu mutants (see below).
Figure 1. clu structure and expression.

A. Map of the clu locus showing sites of insertion mutation (triangles) and the translation start (arrow). Black boxes are exons. B. Drosophila Clu compared with human Clu and Dictyostelium (Dicty) CluA. The position and % amino acid identity of the Clu and TPR domains are shown. Segments used to prepare N-terminal and C-terminal-specific antibodies are indicated by heavy bars. C. Western blot of adult extracts probed with N-terminal anti-Clu antibody showing the effect of clu mutations (listed at top) on the 160 kd Clu protein.
Flies mutant for clu display abnormal mitochondrial structure and function
To analyze clu function, we identified multiple strains bearing single insertions located within the clu transcription unit (Table 1, Fig. 1A). Flies homozygous for clud08713 or cluf04554 lack detectable Clu protein on Western blots, while clud00713 or cluEP0969 homozygotes contain reduced amounts (Fig. 1C). These clu mutant flies showed striking defects. Putative null clud08713 (Fig. 2A) or cluf04554 homozygous flies (data not shown) are small, uncoordinated, sterile, and live only 3-7 days. They move slowly, do not fly, and frequently hold their wings up or down, a sign of flight muscle defects (Fig. 2A-B). The larval fat body fails to turn over within the first two days of adult life, another indication of a metabolic defect (Fig. 2C). Mutant flies bearing the weaker clud00713 or cluEP0969 alleles are viable but male sterile. Flies bearing combinations of clu alleles or trans-heterozygous with a deficiency for the clu region, failed to complement, indicating that these defects are due to disruption of clu expression. Finally, the Clu protein trap line CA06604 (Buszczak et al. 2007), produced a Clu-fusion protein of the expected size (Fig. 1C) that is likely to be functional because CA06604 homozygotes are viable and fertile.
TABLE 1. clueless alleles.
| Allele | Adult viability | Fertility | Comments |
|---|---|---|---|
| d00713 | 44% | Male semi-sterile | |
| EP0969 | 39% | Male sterile | |
| d08713 | 4% | Sterile | Short lived; Small, uncoordinated, wings up/down |
| f04554 | n.d. | Sterile | Short lived; Small, uncoordinated, wings up/down |
| CA06604 | Viable | Fertile | |
| G00271 | Viable | Fertile |
Figure 2. Mitochondrial defects in clu mutants.

A. Wild type adult female. B. clud08713 adult female, showing reduced size and abnormal wing position. C. Light micrograph showing persistent larval fat body cells from abdominal hemolymph of 4 day old clud08713 mutant female. D. Electron micrograph of normal flight muscle showing mitochondria (arrow, inset). E. clud08713 flight muscle showing abnormal, swollen mitochondria with extensive vacuolization of the inner membrane (arrow). F. Normal mitochondria from clud08713 / CyO ovarian germ cells. G. Abnormal mitochondria from clud08713 ovaries showing partial vacuolization of the inner membrane. H-J. CVα expression (green) labeling mitochondrial derivatives in elongating spermatids from cluEP0969/CyO (H), cluEP0969 (I) and park1/Df (J). In contrast to the uniform, paired mitochondrial derivatives seen normally (H, arrows), the mitochondria in both mutant animals (I, J) are uneven with swollen regions (arrowheads). nuclei (DAPI, blue). Scale bars = 1 μM, G for D-G; = 5 μM, J for H-J.
Flies mutant for park, pink1, or rho7, genes which encode mitochondrial proteins, exhibit many similar defects (Greene et al. 2003; Wang et al. 2006; McQuibban et al. 2006). To determine if changes in mitochondrial function underlie the abnormalities in clu mutant flies, we carried out electron microscopy on the flight muscles of 3-day old clud08713 adults. Nearly all the flight muscle mitochondria were swollen and showed severe abnormalities, especially vesiculation/loss of cristae (Fig. 2D, E). Mitochondria in the ovary showed similar but less severe effects, including reduced cristae and swelling (Fig. 2F, G). Mature sperm from sterile clu males were completely immobile (data not shown). In addition, the mitochondrial derivative in elongated spermatids was irregular and swollen (I) compared to wild type (Fig. 2H). We conclude that clu is required for normal mitochondrial structure and function.
To examine whether Drosophila clu is also needed for mitochondrial positioning, we analyzed mitochondria during oogenesis. Drosophila ovarian follicles arise from stem cells in the germarium of each ovariole, and subsequently develop in sequence (Fig. 3A). Changes in mitochondrial localization were observed in clu mutant germ cells beginning as early as the germline stem cells (GSCs). Normally, mitochondria in GSCs are enriched anteriorly, but are also found throughout the cytoplasm (Fig. 3B). In clud08713 mutant GSCs nearly all mitochondria occupy the anterior region of the cell (Fig. 3C). This increased clustering raises the question of whether during GSC division the daughter cystoblast and newly formed GSC inherit equal amounts of mitochondria. At mitosis in wild type GSCs and cystoblasts, mitochondria spread uniformly across the outer spindle and are inherited equally by the daughter cells (Fig. 3D; Cox and Spradling, 2003). In contrast, mitochondria remained clustered during mitosis in cluEP0969 early germ cells (Fig. 3E, arrow), suggesting that inheritance may be unequal, and that one cell may have to make up for a deficit of mitochondria by increased mitochondrial replication (Cox and Spradling, 2006).
Figure 3. clu mutations cause mitochondrial aggregation.

A. Diagram of a Drosophila ovariole, showing the nurse cells (grey) follicle cells (green) and progressively maturing follicles. The location of the germarium and a stage 5 (st5) follicle are indicated. In the germarium, germ line stem cells (GSCs) contain a specialized organelle called the spectrosome, that grows and interconnects the 16-cell germline cysts where it is called the fusome (magenta).
B. Wild type (WT) GSC showing dense mitochondria (green) at the anterior near the spectrosome (magenta circle) as well as throughout the cytoplasm. Note: the nucleus occupies the center of the cell. C. clud08713 GSCs showing increased mitochondrial clustering around spectrosomes. D-E. Mitochondria disperse across the spindle of a wild type single germ cell (D), but remain clustered in a similar mitotic clud08713 germ cell (E, arrow). DAPI (blue) stains the metaphase plate. F-H. Mitochondria disperse throughout the cytoplasm of early wild type cysts (F), but cluster at the ends of the fusome in clud08713 (G) and cluEP0969 (H) mutant cysts. I-J. Mitochondria are dispersed throughout the cytoplasm of stage 5 nurse cells in wild type (I), but cluster (arrows) in clud08713 nurse cells (J). K-L. Mitochondria in germarial follicle cells are dispersed in wild type (K), but clump extensively in cluEP0969 (L). White dashed line outlines: germ cells (B-E), or germline cysts (F-J). Fusome, lateral membranes (1B1, magenta); mitochondria (CVα, green). Scale bars = 10 μM, H for B-H; J for I, J; L for K, L.
Enhanced clustering was also observed in older germ cells in the absence of normal Clu function. Mitochondria are dispersed in young wild type germline cysts (Fig. 3F) but were highly clustered in similarly aged cysts homozygous for clud08713 or cluEP0969 where they associate with the tips of the fusome (Fig. 3G, H), an organelle that interconnects all the germ cells in a cyst. In growing ovarian follicles, mitochondria normally spread evenly throughout the nurse cell cytoplasm (Fig. 3I), but in similarly aged clud08713 follicles nearly all mitochondria cluster in a giant aggregate in each nurse cell (Fig 3J). clu mutations also caused mitochondrial clustering in some somatic follicle cells in the germarium (Fig. 3K, L).
Clu is cytoplasmic and may mediate mitochondrial movement along microtubules
We next examined the distribution of Clu protein to gain insight into its requirement for mitochondrial positioning and membrane integrity. Both the N-terminal and C-terminal antibodies, and a Clu-GFP protein trap gave indistinguishable results. Clu protein is found exclusively in the cytoplasm and is often present in discrete bodies or aggregates as well as throughout the cytoplasm. (Fig. 4A, B). Lower labeling is observed in somatic cells. Antibody reactivity is abolished in clu null mutations in all cells (Fig. 4C). Double labeling with anti-Clu and a mitochondrial antibody, reveals that multiple mitochondria associate with Clu particles (Fig. 4D). However, Clu protein is only observed outside of mitochondria (Fig. 4E, E′).
Figure 4. Clu protein is a cytoplasmic but not a mitochondrial protein.

A-C. Stage 5 follicles: In wild type (A), Clu antibody (magenta) stains the nurse cell cytoplasm including some particles. Lateral membranes (1B1, green). The distribution of Clu-GFP in the protein trap line CA06604 visualized with anti-GFP antibody (magenta) was indistinguishable (B). In clud08713, Clu antibody staining is absent (C). Lateral membranes (Hts, green). D-H. Stage 9 follicles: D. Double label of Clu (magenta, arrow) and mitochondria (green) shows that Clu particles associate with one or more mitochondria. Juxtaposition of mitochondria (E, green) and Clu protein (E, E′ magenta) can be seen more clearly at higher magnification. F. Association of microtubules (green) with some Clu bodies (magenta). G, H. Microtubules in clud08713 stage 9 nurse cells (H) appear normal compared to a wild type control (G). Scale bars = 10uM, C for A-C and H for D, F, G, H; 2.5um E′, E.
Ovarian mitochondria frequently move along microtubules as a result of tethering to both Dynein and Kinesin motors via linker proteins such as Milton (Cox and Spradling, 2006). Examining ovaries that were stained for Clu and microtubules revealed that Clu particles sometimes associate with microtubules (Fig. 4F). The distribution of microtubules in the ovary did not change appreciably in clud08713 nurse cells where mitochondria are highly clustered (Fig. 4G, H).
These observations argue that mitochondrial clustering in clu mutants is not caused by changes in microtubule organization. However, comparing the location of mitochondrial clusters in clu mutants and in follicles mutant for microtubule motor protein genes or for the milton linker (Cox and Spradling, 2006), strongly suggested that Clu affects the movement of mitochondria along microtubules. In particular, sites of mitochondrial accumulation in clu mutants such as the proximal side of germline stem cells (Fig. 3C) and adjacent to the fusome in early cysts (Fig. 3G,H), correspond to sites of microtubule plus ends in these cells (Cox and Spradling, 2006). Consequently, clu mutations cause mitochondria to cluster at the sites of microtubule plus ends in several stages where the polarity of microtubules that transport mitochondria are known. In other situations, such as in nurse cells, mitochondria are less affected by microtubule motor mutations and the systems controlling their location as well as microtubule polarity have yet to be determined.
Mitochondria cluster in park mutant cells
Many of the effects of clu mutation we observed resembled those described following disruption of park (Greene et al. 2003). Consequently, we examined mitochondria in park mutant follicles to see if they also form clusters. park1 / Df females are sterile, but follicles develop far enough to reveal strong mitochondrial clustering in the nurse cells of previtellogenic follicles (Fig. 5A, B). As follicles age, mitochondrial clustering worsens (Fig. 5C, D). Examining these clustered mitochondria using light and electron microscopy revealed that some were doughnut shaped (Fig. 5E, F), and some grew larger than normal (Fig. 5F; compare Fig. 2F). park mutants also induced clustering in somatic follicle cells (Fig. 5G, H). During spermatogenesis, only one mitochondrial derivative forms in park1 / Df mutant germ cells (Fig. 5P) as previously reported (Riparbelli and Callaini, 2007). cluEP0969 forms two derivatives but they appear frayed (Fig. 5N-O). Despite this initial difference, during spermatid elongation the park1 / Df mutant mitochondrial derivative becomes irregular and swollen like those in cluEP0969 (Fig. 2I-J).
Figure 5. park mutations cause mitochondrial clustering and interact genetically with clu.

A. Wild type stage 5 follicle lacks clustered mitochondria (green). 1B1 (magenta) stains membranes. B. Severe mitochondrial clustering in stage 5 park1/Df follicle. C. Wild type stage 9 follicle lacks clustered mitochondria. D. Extreme mitochondrial clustering in stage 9 park1/Df follicle. E. Mitochondria in parkZ472 are clustered, long (arrow) and frequently form ring structures (arrowhead). F. Electron micrograph of a parkZ472 mitochondrial cluster showing doughnut-shaped mitochondria (arrowhead). G. Wild type stage 9 follicle cells mostly lack clustered mitochondria (arrow). H. Mitochondria cluster in the cell center in most park1/Df stage 9 follicle cells (arrow). I-J. Clu protein (white) levels are present at similar levels in wild type (I) and in park1/Df (J) follicles, but large Clu particles are absent in the mutant. K-M. Normal dispersed mitochondrial distribution in cluEP969/CyO (K) and park1/TM6b (L) stage 5 follicles. M. Clustered mitochondria in stage 5 follicle from cluEP0969 / + ; park1/+ female. N-P. Mitochondrial derivative in cluEP0969/CyO (N), cluEP0969 (O), and (P) park1/Df male germ cells at the leaf blade stage. In wild type (N), the Nebenkern unfurls revealing two mitochondria (asterisks) elongating below the nucleus (dotted circle). These mitochondria are less compact and appear frayed in cluEP0969 mutants (O). park1/Df (P) have only one mitochondrion (asterisk) that looks smooth as in (N). green = mitochondria, CVα; magenta = membranes, 1B1. Scale bars = 10uM, M for A, B, G, H, K-M; D for C, D; J for I, J; = 5uM in P for N-P
clu and park interact
Because of the similarities between clu and park, we investigated whether these two genes interact genetically. park mutant ovaries still contained normal amounts of Clu protein (Fig. 5I, J), so the clustering observed in these animals is not simply due to reduced Clu expression. The only difference appeared to be a reduction in the number of Clu particles in the park mutant cells. Females heterozygous for either clu or park alone were phenotypically wild type and mitochondria did not cluster in nurse cells from stage 5 follicles (Fig. 5K, L). In contrast, mitochondria were noticeably clustered in 100% of similarly aged follicles from park1 / + ; cluEP0969 / + females (Fig. 5M). Evidence of severe mitochondrial dysfunction was not observed, however, and both male and female trans-heterozygotes were fertile. We conclude that clu and park genetically interact to promote mitochondrial clustering.
Gene expression changes in clu mutant cells
To further analyze the roles played by Clu and Park we investigated if clu or park mutant animals contain differences in major mitochondrial proteins that are encoded in the cell nucleus. Western blots revealed both pyruvate dehydrogenase and Complex V-alpha chain levels were significantly reduced in null clud08713 homozygotes. This decrease was not observed in cluEP0969 homozygotes, which retain a small amount of Clu protein. The observed reduction indicates that clud08713 mitochondria contain reduced amounts of these major mitochondrial proteins, or that the total number of mitochondria is reduced in clud08713 homozygotes.
To try and identify genes that are directly impacted by clu mutation, we analyzed stage 10B ovarian follicles bearing the weaker cluEP0969 allele using microarrays. Such follicles appear normal and successfully complete oogenesis, however, Clu protein is significantly depleted (Fig. 1C), and mitochondria cluster early in germ cell development. These studies revealed a small number of high-level changes in gene expression between wild type and cluEP0969 follicles (Table 2).
TABLE 2. Gene expression changes in cluEP0969 stage 10B follicles.
| Gene | fold change | comment |
|---|---|---|
| Top downreglated genes (<0.15 of control) | ||
| CG14036 | 0.013 | related to gametocyte specific factor 1 |
| CG18446 | 0.014 | SALL4 related; regulates survival and apoptosis in mammals |
| U26 | 0.019 | putative aminoadipate-semialdehyde dehydrogenase |
| CG5181 | 0.022 | single-stranded DNA binding protein/ DNA repair |
| SCAP | 0.028 | sterol regulatory element binding protein cleavage |
| CG12896 | 0.051 | thiol specific antioxidant |
| Glaz | 0.068 | lipid metabolic process |
| Pde11 | 0.083 | cGMP-specific phosphodiesterase |
| GstE1 | 0.097 | metabolizes products of lipid peroxidation |
| CG10962 | 0.148 | oxidoreductase activity |
| Rep1 | 0.132 | inhibits caspase activated DNAse (CAD), apoptosis |
| CG10354 | 0.135 | 5′-3′ exoribonuclease activity |
| Phr6-4 | 0.149 | (6-4)-photolyase |
| Top upregulated genes (>5X of control) | ||
| CG32912 | 18.0 | N-acetylmuramoyl-L-alanine amidase activity |
| CG12708 | 17.9 | ovary enriched protein |
| CG7650 | 11.1 | phosducin-like protein; thioredoxin_like superfamily |
| path | 10.7 | PAT-class amino acid/fatty acid transporter |
| CG8964 | 10.3 | IgG domain protein resembling human PTK7 |
| CG11880 | 9.70 | putative choline transporter/solute carrier protein |
| GstE10 | 8.33 | glutathione transferase |
| CG32364 | 8.27 | gonad-enriched protein, function unknown |
| CG14636 | 6.0 | function unknown |
| Cyp6a17 | 5.7 | cytochrome P450 protein |
| Cyp6d4 | 5.6 | cytochrome P450 protein |
Several of the most strongly altered genes protect cells from reactive oxygen species that are generated during mitochondrial metabolism. For example, both CG12896 and GstE1, genes thought to counteract reactive oxygen damage to thiols and to lipids, respectively (Sawicki et al. 2003), are strongly decreased by clu mutation. Several putative DNA repair proteins, including CG10534, a 5′3′ exonuclease, CG5181 (SSB1), a single-stranded DNA binding protein, and phr6-4, a 6-4 photolyase, are all strongly lowered. Apopotosis regulators such as the Rep1 caspase inhibitor (Mukae et al. 2000) and the putative SALL4-homolog CG18446, are reduced as much as tenfold. Several genes involved in mitochondrial metabolism are also changed, including SCAP, which regulates fatty acid production (reviewed in Rawson, 2003). In Drosophila, SCAP and its partner SREBP sense 16-carbon fatty acids and control the transcription of target genes involved in fatty acid biosynthesis in mitochondria. Perhaps in response to altered lipid metabolism, transcripts of several other candidate metabolic genes are also altered, including CG3424 encoding an amino acid/fatty acid transporter, CG11880 encoding a putative choline transporter, U26 encoding a putative aminoadipate-semialdehyde dehydrogenase, and others. Thus, the changes in transcript levels we observe in relatively mild clu mutant follicles suggest that multiple cellular functions linked to mitochondria have been compromised. The results suggest that clu mutation reduces mitochondrial function, enhances oxidative damage and decreases repair or apoptosis.
Discussion
clu encodes a highly conserved protein required for mitochondrial localization and function
Our studies identify a novel Drosophila gene, clu, that is essential for the normal cytoplasmic localization and function of cellular mitochondria. Clu-related proteins are highly conserved, especially within a novel 243 amino acid region (the “Clu domain”) which is 85% identical (0 gaps) between Drosophila and human Clu. Clu function also appears to be remarkably conserved because mutations in orthologs in Dictyostelium, Saccharomyces, and Arabidopsis all cause mitochondria to cluster into a restricted cytoplasmic region (Zhu et al. 1997; Fields et al. 1998; Dimmer et al. 2002; Logan et al. 2003). In Drosophila, the Clu pathway is also essential to maintain mitochondrial integrity, while Dictyostelium and yeast cells lacking CluA or Clu1 have no detected functional defects.
park interacts with clu
Disruptions in a conserved process of mitochondrial maintenance involving the park and pink1 genes underlies many cases of Parkinson's disease (Abou-Sleiman et al. 2006). We find that Drosophila clu mutations cause phenotypes that strongly resemble those caused by park mutations. Moreover, we showed that mitochondria in park mutants cluster in ovarian cells in a very similar manner to their behavior in clu mutants. In addition, park and clu genetically interact; compound heterozygotes have a mitochondrial clustering phenotype not seen in either heterozygote alone. These observations argue that Clu and Park function in the same or closely related pathways, and that the human Clu ortholog KIAA0664 might play a role in Parkinson's disease. Furthermore, these findings suggest that changes in mitochondrial localization may be of significance in this condition.
Mitochondria cluster at microtubule plus ends in the early germ cells of clu mutants
Mitochondria are often positioned within cells by motor-dependent movement along microtubules (reviewed in Hollenbeck and Saxton, 2005). Clu protein appears to contact mitochondria, especially in the large Clu particles, many of which are located adjacent to microtubules. Our experiments suggest that clu-induced mitochondrial aggregates are caused by changes in microtubule-based mitochondrial transport. In both the Drosophila brain and ovary, mitochondria associate with transport complexes containing the adaptor protein Milton that are linked to both plus end-directed motors such as Kinesin and minus-end directed motors such as Dynein (Glater et al. 2006; Cox and Spradling, 2006). During oogenesis, loss of Khc, or of one of the Milton isoforms that interacts with Khc, cause mitochondria to cluster in cellular regions rich in microtubule “−“ ends, presumably because countervailing movement toward “+” ends has been lost (see Fig. 5 in Cox and Spradling, 2006). Strikingly, mitochondria in clu mutant ovaries accumulate at predicted sites of microtuble “+” ends. These include the proximal region of germline stem cells (Fig. 3C) and against the fusome in early cysts (Fig. 3G, H). Thus, clu mutations may act by interfering with minus-end directed mitochondrial movement along microtubules in some situations. Whether Clu particles correspond to mitochondrial transport complexes that mediate these changes, or whether Clu acts indirectly on mitochondrial positioning could not be determined with certainty from our studies.
In other cells, such as nurse cells, clu mutations also cause pronounced mitochondrial clustering (Fig. 3J), but that clustering is not mimicked by mutations in genes encoding microtubule motor proteins or Milton (Cox and Spradling, 2006). Other systems of mitochondrial localization may predominate in such cells, but still be subject to regulation by the Clu pathway. For example, in plant cells, Clu was postulated to control mitochondrial localization by regulating the choice between microtubule-dependent and microfilament-dependent transport (Logan et al. 2003). In particular, the TPR domain of Clu was proposed to repress the interaction of mitochondria with microtubules by competing for binding between a similar TPR region in kinesin light chain, thereby allowing transport along actin to predominate. The loss of this postulated interaction via clu mutation predicts that mitochondria would accumulate at microtubule minus ends in these cells. Although the default position where mitochondria accumulate in the absence of Clu might vary depending on the particular cell type and transport systems involved, the use of a Clu-dependent pathway to control mitochondrial subcellular location appears to be widespread.
Clu and Park may function in a novel pathway affecting mitochondrial physiology
Our experiments show that mitochondria require Clu to maintain their subcellular location and structural integrity. In the presence of even a relatively mild clu mutation, transcripts from genes involved in mitochondrial function and in protection from oxidative damage are reduced. Consequently, we propose that the normal role of Clu is to function in a pathway that controls the location and activity of mitochondria within the cell. Frequently, mitochondria may need to move to a different subcellular location in order to maximize access to substrates or to mitigate damage from oxidative metabolism. For example, in neurons, microtubule-dependent mitochondrial transport is modulated based on the level of respiratory activity (Miller and Sheetz, 2004). Clu protein would participate in a pathway that senses the internal physiological state of individual mitochondria and transduces this information into homeostatic changes in its position and metabolic activity. When the Clu pathway is impaired, mitochondria would not move or operate normally, and might consequently suffer damage, a scenario consistent with the changes we observed in the inner mitochondrial membranes of Clu mutants, with their reduced levels of mitochondrial enzymes, and with the observed changes in nuclear gene expression. The greater severity of clu mutations in Drosophila compared to Dictyostelium and yeast clu might be due to an intrinsically greater requirement for dynamic mitochondrial positioning in the specialized cells of complex metazoans. This model may provide a rationale for the unexpected effects of microtubule inhibitors on mitochondrial function that were reported recently (Wagner et al. 2008).
Fig. 7 illustrates how the clu and park pathway(s) might link microtubule-based mitochondrial transport, mitochondrial physiology and oxidative damage. Under normal conditions, the clu and park pathway(s) would sense the physiological state of mitochondria and activate appropriate levels of minus-end directed mitochondrial movement along microtubules. For example, mitochondria low in respiratory substrates might move to cellular regions where these substrates are abundant. Mitochondria in need of repair might move close to the nucleus where appropriate repair genes would be induced. In conjunction with active plus end-directed motors, the result would be that mitochondria move dynamically to locations throughout the cell appropriate to their physiological state. However, without clu and park pathway(s) function or if the local level of toxic metabolic products exceeded a threshold, the mitochondria in question would cease minus-end directed transport and undergo concerted plus-end directed movement. In many cells, the major locus of microtubule minus ends is found near the nucleus so that plus end-directed movement would increase the distance between reactive oxygen production and the nuclear DNA.
Figure 7. Model for clu and park function.
The figure shows a model for the role of the clu and park pathway(s) in cells (such as early ovarian germ cells) in which mitochondrial movement is controlled by the regulated engagement of plus end-directed kinesin motors (pink) and minus-end directed dynein motors (blue). The clu and park pathway(s) are located in the cytoplasm where they are proposed to sense the physiological state of a mitochondrion and modify its movement accordingly. A. Under normal conditions, modulating the level of Dynein engagement relative to plus-end directed motor activity allows the clu and park pathway(s) to cause a mitochondrion to move in either direction (arrows). B. When the functioning of the clu and park pathway(s) is defective, or when levels of reactive oxygen species are elevated, minus-end microtubule movement is blocked, causing mitochondria to move preferentially toward microtubule plus ends (large arrow). Note that in other cells, the clu and park pathway(s) might regulate mitochondrial movement by linking to different motors and transport systems.
Clu and mitochondrial inheritance
During early germ cell development in many organisms, a subset of mitochondria within the germ cell cyst is transported and accumulates in a concerted fashion to form the Balbiani body (Cox and Spradling, 2003; Kloc et al, 2004). It has been postulated that Balbiani body formation is a selective process that rids the germline of mitochondria with defective mitochondrial genomes (Pepling et al. 1999; Cox and Spradling, 2003). In Drosophila, mitochondrial movement leading up to Balbiani body formation takes place by minus-end directed transport along microtubules (Cox and Spradling, 2003; 2006). Consequently, the Clu pathway might contribute to mitochondrial screening. However, Balbiani bodies of normal size were observed to form in clu mutants, suggesting that Clu function is not utilized or is redundant at this critical stage of germ cell development.
The Clu pathway and Parkinson's disease
Our studies have several implications for understanding Parkinson's disease. According to our model, the Clu pathway would contribute strongly to the ability of mitochondria to remain functional during aging, despite the high metabolic requirements and environmental stresses experienced by many tissues. This may be particularly important in neural cells such as those that are compromised in Parkinson's disease. In a cell whose Clu pathway is compromised, mitochondria would operate less efficiently, suffer more damage and wear out faster because they would spend more time in cellular locations that are not appropriate to their metabolic state. The number of mitochondria producing elevated levels of ROS would rise, increasing reactive oxygen damage to the nucleus and other cellular components, thereby leading to greater cell death. Consequently, the level of Clu pathway function within an individual might influence their susceptibility to sporadic Parkinson's disease and to other late-onset neurological disorders.
Mitochondrial mutations may fall into two classes based on their phenotype, and on their relationship to the Clu pathway. Some mutations, such as those affecting mitochondrial ribosomal proteins (Royden et al. 1987), or certain metabolic pathways (Zhang et al. 1999), cause bang-sensitivity, but differ phenotypically from clu mutation (Fergestad et al. 2006). In contrast, mutations in park (Greene et al. 2003), pink1 (Wang et al. 2006), rho7 (McQuibban et al. 2006), mitochondrial ATP6 (Celotto et al. 2006), and mitochondrial cytochrome oxidase (Xu et al. 2008) resemble clu mutations in causing general defects in movement, flight, muscle or nerve degeneration, male fertility and longevity. The effects of this latter class may be too severe to be compensated by Clu pathway operation, leading to mitochondrial mis-positioning and increased production of reactive oxygen species. Consistent with this, expression of anti-oxidants partially suppresses the affects of pink1 on neurodegeneration (Wang et al. 2006). Mutations in some of these genes, in addition to park, may also cause mitochondrial clustering and genetically interact with clu. Clearly, further studies of the Clu pathway will deepen our understanding of how mitochondria are maintained in cells and why they sometimes become damaged with age.
Materials and Methods
Drosophila strains
All symbols are as described in FlyBase. Park[1] (Cha et al, 2005) was obtained from Dr. Kyoung Sang Cho, and park[Z472] (Greene et al, 2003) was obtained from Leo Pallanck. WT flies were CantonS, except for Fig. 3C, which is OregonR. All other stocks were obtained from Bloomington or Harvard. All stocks were cultured at 22-25°C on standard food.
Antibody production
Antibodies were raised to amino acids 1-104 and 1301-1448 of Clu using the protocol described in Buszczak and Spradling (2006).
Immunofluorescene Microscopy
Fattened females and males were dissected and fixed following the protocol in Cox and Spradling (2003). Primary antibodies were diluted at follows: mouse Complex V, alpha subunit (1:500, MitoSciences), mouse 1B1 (1:100, Developmental Studies Hybridoma Bank), ginea pig Clu-C (1:2000). The following secondaries were used: goat anti-mouse IgG2b and goat anti-mouse AlexaFluor488, goat anti-mouse IgG1 AlexaFluor 568, goat anti-ginea pig AlexaFluor 568 (1:400, Molecular Probes). For microtubule labeling, fattened ovaries were fixed for fifteen minues in 4% paraformaldehyde in Graces with 1mM EGTA. Primary antibodies were a 1:1:1 mix of E7, AA 12.1, AA 4.3 (1:50, Developmental Studies Hybridoma Bank). For DNA labeling, DAPI was added 5 ug/ml for 10 minutes prior to adding VectoShield (Vector Labs, Inc.). Confocal analysis was done using a 63× (NA 1.32) PlanApo lens and Leica TCS SP2 or SP5 confocal microscopes, or a Zeiss LSM510 Pascal.
Electron microscopy
Adult flight muscles and fattened ovaries were dissected in Grace's media (Sigma), then processed and imaged as described in Cox and Spradling (2006).
Western blotting
Western blotting was done as described in Wilhelm et al (2000) with the following modifications. To isolate protein, adults were ground with a pestle in Laemmli buffer + BME. Sample was boiled seven minutes before loading on a 4-15% acrylamide gel. The following primary antibodies were used: goat α-Clu-N was used 1:15,000, mouse α-actin (JLA-20, 1:200, Sigma), mouse α-Pyruvate dehydrogenase E1 α-subunit (PDF, 1:500, Mitosciences), mouse α-Complex V, α-subunit (CVα, 1:10,000, Mitosciences)
Microarray analyses
RNA was extracted from 100 stage 10B ovarian follicles as described in Allen and Spradling (2008). Microarray analyses were performed in duplicate with RNA samples isolated on separate days, and greater than 90% of the genes that changed significantly were reproducible.
Supplementary Material
Figure 6. clu adults contain reduced amounts of mitochondrial protein.

A. Western blot of whole fly extracts probed for Clu, Pyruvate dehydrogenase (PDH), Complex V-alpha subunit (CV-α), and Actin. The level of the two mitochondrial proteins PDH and CVα were unchanged from wild type levels in heterozygotes mutant for clu or park. However, both PDH and CVα levels were reduced significantly relative to the actin control in strong clud08713 mutants (but not in cluEP0969).
Acknowledgments
We thank Drs. L. Pallanck and Kyoung Sang Cho for materials. Don Fox, Rebecca Frederick, Mary Goll and Vicki Losick made useful comments on the manuscript. Shelley Lloyd and Megan Kutzner provided technical assistance and Mike Sepanski provided technical assistance with electron microscopy. N Rusan and M. Peifer provided the protocol for labeling microtubules.
References
- Abou-Sleiman PM, Muqit MM, Wood NW. Expanding insights of mitochondrial dysfunction in Parkinson's disease. Nat Rev Neurosci. 2006;7:207–19. doi: 10.1038/nrn1868. [DOI] [PubMed] [Google Scholar]
- Betarbet R, Sherer TB, Greenamyre JT. Ubiquitin-proteasome system and Parkinson's diseases. Exp Neurol. 2005;191:S17–S27. doi: 10.1016/j.expneurol.2004.08.021. [DOI] [PubMed] [Google Scholar]
- Buszczak M, Paterno S, Lighthouse D, Bachman J, Plank J, et al. The Carnegie protein trap library: a versatile tool for Drosophila developmental studies. Genetics. 2007;175:1505–31. doi: 10.1534/genetics.106.065961. [DOI] [PMC free article] [PubMed] [Google Scholar]
- Celotto AM, Frank AC, McGrath SW, Fergestad T, Van Voorhies WA, Buttle KF, Mannella CA, Palladino MJ. Mitochondrial Encephalomyopathy in Drosophila. J Neurosci. 2006;26:810–820. doi: 10.1523/JNEUROSCI.4162-05.2006. [DOI] [PMC free article] [PubMed] [Google Scholar]
- Clark IE, Dodson MW, Jiang C, Cao JH, Huh JR, Seol JH, Yoo SJ, Hay BA, Guo M. Drosophila pink1 is required for mitochondrial function and interacts genetically with parkin. Nature. 2006;441:1162–1166. doi: 10.1038/nature04779. [DOI] [PubMed] [Google Scholar]
- Cox RT, Spradling AC. A Balbiani body and the fusome mediate mitochondrial inheritance during Drosophila oogenesis. Development. 2003;130:1579–90. doi: 10.1242/dev.00365. [DOI] [PubMed] [Google Scholar]
- Cox RT, Spradling AC. Milton controls the initial acquisition of mitochondria by Drosophila oocytes. Development. 2006;133:3371–7. doi: 10.1242/dev.02514. [DOI] [PubMed] [Google Scholar]
- Dimmer KS, Fritz S, Fuchs F, Messerschmitt M, Weinbach N, Neupert W, Westermann B. Genetic basis of mitochondrial function and morphology in Saccharomyces cerevisiae. Mol Biol Cell. 2002;13:847–53. doi: 10.1091/mbc.01-12-0588. [DOI] [PMC free article] [PubMed] [Google Scholar]
- Fergestad T, Bostwick B, Ganetzky B. Metabolic Disruption in Drosophila Bang-Sensitive Seizure Mutants. Genetics. 2006;173:1357–1364. doi: 10.1534/genetics.106.057463. [DOI] [PMC free article] [PubMed] [Google Scholar]
- Fields SD, Conrad MN, Clarke M. The S. cerevisiae CLU1 and D. discoideum cluA genes are functional homologues that influence mitochondrial morphology and distribution. J Cell Sci. 1998;111:1717–27. doi: 10.1242/jcs.111.12.1717. [DOI] [PubMed] [Google Scholar]
- Fields SD, Arana Q, Heuser J, Clarke M. Mitochondrial membrane dynamics are altered in cluA- mutants of Dictyostelium. J Muscle Res Cell Motil. 2002;23:829–38. doi: 10.1023/a:1024492031696. [DOI] [PubMed] [Google Scholar]
- Glater EE, Megeath LJ, Stowers RS, Schwarz TL. Axonal transport of mitochondria requires Milton to recruit kinesin heavy chain and is light chain independent. J Cell Bio. 2006;173:545–557. doi: 10.1083/jcb.200601067. [DOI] [PMC free article] [PubMed] [Google Scholar]
- Greene JC, Whitworth AJ, Kuo I, Andrews LA, Feany MB, Pallanck LJ. Mitochondrial pathology and apoptotic muscle degeneration in Drosophila parkin mutants. Proc Natl Acad Sci USA. 2003;100:4078–4083. doi: 10.1073/pnas.0737556100. [DOI] [PMC free article] [PubMed] [Google Scholar]
- Hollenbeck PJ, Saxton WM. The axonal transport of mitochondria. J Cell Sci. 2005;118:5411–9. doi: 10.1242/jcs.02745. [DOI] [PMC free article] [PubMed] [Google Scholar]
- Kitada T, Asakawa S, Hattori N, Matsumine H, Yamamura Y, Minoshima S, Yokochi M, Mizuno Y, Shimizu N. Mutations in the parkin gene cause autosomal recessive juvenile parkinsonism. Nature. 1998;392:605–608. doi: 10.1038/33416. [DOI] [PubMed] [Google Scholar]
- Kloc M, Bilinski S, Etkin LD. The Balbiani body and germ cell determinants: 150 years later. Curr Top Dev Biol. 2004;59:1–36. doi: 10.1016/S0070-2153(04)59001-4. [DOI] [PubMed] [Google Scholar]
- Logan DC, Scott I, Tobin AK. The genetic control of plant mitochondrial morphology and dynamics. Plant J. 2003;36:500–9. doi: 10.1046/j.1365-313x.2003.01894.x. [DOI] [PubMed] [Google Scholar]
- McQuibban GA, Lee JR, Zheng L, Juusola M, Freeman M. Normal mitochondrial dynamics requires rhomboid-7 and affects Drosophila lifespan and neuronal function. Curr Biol. 2006;16:982–989. doi: 10.1016/j.cub.2006.03.062. [DOI] [PubMed] [Google Scholar]
- Miller KE, Sheetz MP. Axonal mitochondrial transport and potential are correlated. J Cell Sci. 2004;117:2791–2804. doi: 10.1242/jcs.01130. [DOI] [PubMed] [Google Scholar]
- Mukae N, Yokoyama H, Yokokura T, Sakoyama Y, Sakahira H, Nagata S. Identification and developmental expression of inhibitor of caspase-activated DNase (ICAD) in Drosophila melanogaster. J Biol Chem. 2000;275:21402–21408. doi: 10.1074/jbc.M909611199. [DOI] [PubMed] [Google Scholar]
- Pallanck L, Greenamyre JT. Pink, parkin and the brain. Nature. 2006;441:1058. doi: 10.1038/4411058a. [DOI] [PubMed] [Google Scholar]
- Park J, Lee SB, Lee S, Kim Y, Song S, Kim S, Bae E, Kim J, Shong M, Kim JM, Chung J. Mitochondrial dysfunction in Drosophila PINK1 mutants is complemented by parkin. Nature. 2006;441:1157–1161. doi: 10.1038/nature04788. [DOI] [PubMed] [Google Scholar]
- Pepling ME, deCuevas M, Spradling AC. Germline cysts: a conserved phase of germ cell development? Trends in Cell Biology. 1999;9:257–261. doi: 10.1016/s0962-8924(99)01594-9. [DOI] [PubMed] [Google Scholar]
- Poole AC, Thomas RE, Andrews LA, McBride HM, Whitworth AJ, Pallanck LJ. The PINK1/Parkin pathway regulates mitochondrial morphology. Proc Natl Acad Sci USA. 2008;105:1638–1643. doi: 10.1073/pnas.0709336105. [DOI] [PMC free article] [PubMed] [Google Scholar]
- Rawson RB. The SREBP pathway- insights from Insigs and insects. Nat Rev Mol Cell Biol. 2003;8:631–40. doi: 10.1038/nrm1174. [DOI] [PubMed] [Google Scholar]
- Riparbelli MG, Callaini G. The Drosophila parkin homologue is required for normal mitochondrial dynamics during spermiogenesis. Dev Biol. 2007;303:108–120. doi: 10.1016/j.ydbio.2006.10.038. [DOI] [PubMed] [Google Scholar]
- Royden CS, Pirrotta V, Jan LY. The tko locus, site of a behavioral mutation in D. melanogaster, codes for a protein homologous to prokaryotic ribosomal protein S12. Cell. 1987;51:165–173. doi: 10.1016/0092-8674(87)90144-9. [DOI] [PubMed] [Google Scholar]
- Sawicki R, Singh SP, Mondal AK, Benes H, Zimniak P. Cloning, expression and biochemical characterization of one Epsilon-class (GST-3) and ten Delta-class (GST-1) glutathione S-transferases from Drosophila melanogaster, and identification of additional nine members of the Epsilon class. Biochem J. 2003;370:661–669. doi: 10.1042/BJ20021287. [DOI] [PMC free article] [PubMed] [Google Scholar]
- Schlossmacher MG, Frosch MP, Gai WP, Medina M, Sharma N, Forno L, Ochiishi T, Shimura H, Sharon R, Hattori N, Langston JW, Mizuno Y, Hyman BT, Selkoe DJ, Kosik KS. Parkin localizes to the Lewy bodies of Parkinson disease and dementia with Lewy bodies. Am J Pathol. 2002;160:1655–1667. doi: 10.1016/S0002-9440(10)61113-3. [DOI] [PMC free article] [PubMed] [Google Scholar]
- Spradling AC. In: Developmental genetics of oogenesis The Development of Drosophila melanogaster. Bate M, Martizen-Arias A, editors. Vol. 1. Cold Sprading Harbor, NY: Cold Spring Harbor Press; 1993. pp. 1–70. [Google Scholar]
- Valente EM, Abou-Sleiman PM, Caputo V, Muqit MM, Harvey K, Gispert S, Ali Z, Del Turco D, Bentivoglio AR, Healy DG, Albanese A, Nussbaum R, González-Maldonado R, Deller T, Salvi S, Cortelli P, Gilks WP, Latchman DS, Harvey RJ, Dallapiccola B, Auburger G, Wood NW. Hereditary early-onset Parkinson's disease caused by mutations in PINK1. Science. 2004;304:1158–1160. doi: 10.1126/science.1096284. [DOI] [PubMed] [Google Scholar]
- Wagner BK, et al. Large-scale chemical dissection of mitochondrial function. Nature Biotechnol. 2008;26:343–351. doi: 10.1038/nbt1387. [DOI] [PMC free article] [PubMed] [Google Scholar]
- Wang D, Qian L, Xiong H, Liu J, Neckameyer Ws, Oldham S, Xia K, Wang J, Bodmer R, Zhang Z. Antioxidants protect PINK1-dependent dopaminergic neurons in Drosophila. Proc Natl Acad Sci USA. 2006;103:13520–13525. doi: 10.1073/pnas.0604661103. [DOI] [PMC free article] [PubMed] [Google Scholar]
- Xu H, DeLuca SZ, O'Farrell PH. Manipulating the metazoan mitochondrial genome with targeted restriction enzymes. Science. 2008;321:575–7. doi: 10.1126/science.1160226. [DOI] [PMC free article] [PubMed] [Google Scholar]
- Yang Y, Gehrke S, Imai Y, Huang Z, Ouyang Y, Wang JW, Yang L, Beal MF, Vogel H, Lu B. Mitochondrial pathology and muscle and dopaminergic neuron degeneration caused by inactivation of Drosophila Pink1 is rescued by Parkin. Proc Natl Acad Sci USA. 2006;103:10793–10798. doi: 10.1073/pnas.0602493103. [DOI] [PMC free article] [PubMed] [Google Scholar]
- Yang Y, Ouyang Y, Yang L, Beal MF, McQuibban A, Vogel H, Lu B. Pink1 regulates mitochondrial dynamics through interaction with the fission/fusion machinery. Proc Natl Acad Sci USA. 2008;105:7070–5. doi: 10.1073/pnas.0711845105. [DOI] [PMC free article] [PubMed] [Google Scholar]
- Zhang YQ, Roote J, Brogna S, Davis AW, Barbash DA, Nash D, Ashburner M. stress sensitive B encodes an adenine nucleotide translocase in Drosophila melanogaster. Genetics. 1999;153:891–903. doi: 10.1093/genetics/153.2.891. [DOI] [PMC free article] [PubMed] [Google Scholar]
- Zhu Q, Hulen D, Liu T, Clarke M. The cluA- mutant of Dictyostelium identifies a novel class of proteins required for dispersion of mitochondria. Proc Natl Acad Sci USA. 1997;94:7308–13. doi: 10.1073/pnas.94.14.7308. [DOI] [PMC free article] [PubMed] [Google Scholar]
Associated Data
This section collects any data citations, data availability statements, or supplementary materials included in this article.

