Testis-specific ATP synthase peripheral stalk subunits required for tissue-specific mitochondrial morphogenesis in Drosophila
- Eric M. Sawyer†1,
- Elizabeth C. Brunner†1,
- Yihharn Hwang1,
- Lauren E. Ivey1,
- Olivia Brown1,
- Megan Bannon1,
- Dennis Akrobetu1,
- Kelsey E. Sheaffer1,
- Oshauna Morgan1,
- Conroy O. Field1,
- Nishita Suresh1,
- M. Grace Gordon1,
- E. Taylor Gunnell1,
- Lindsay A. Regruto1,
- Cricket G. Wood2,
- Margaret T. Fuller2 and
- Karen G. Hales1, 2Email author
© The Author(s). 2017
Received: 18 November 2016
Accepted: 15 March 2017
Published: 23 March 2017
In Drosophila early post-meiotic spermatids, mitochondria undergo dramatic shaping into the Nebenkern, a spherical body with complex internal structure that contains two interwrapped giant mitochondrial derivatives. The purpose of this study was to elucidate genetic and molecular mechanisms underlying the shaping of this structure.
The knotted onions (knon) gene encodes an unconventionally large testis-specific paralog of ATP synthase subunit d and is required for internal structure of the Nebenkern as well as its subsequent disassembly and elongation. Knon localizes to spermatid mitochondria and, when exogenously expressed in flight muscle, alters the ratio of ATP synthase complex dimers to monomers. By RNAi knockdown we uncovered mitochondrial shaping roles for other testis-expressed ATP synthase subunits.
We demonstrate the first known instance of a tissue-specific ATP synthase subunit affecting tissue-specific mitochondrial morphogenesis. Since ATP synthase dimerization is known to affect the degree of inner mitochondrial membrane curvature in other systems, the effect of Knon and other testis-specific paralogs of ATP synthase subunits may be to mediate differential membrane curvature within the Nebenkern.
Mitochondrial dynamics, including fusion, division, regulated movement, quality control, and interaction with the endoplasmic reticulum, have been characterized largely in yeast and cultured mammalian cells (reviewed in [1, 2]). Studies in specialized cell types with high energy needs in other organisms, such as neurons and sperm in Drosophila melanogaster, also provided early and ongoing insight into molecular mechanisms of mitochondrial shaping [3, 4]. Drosophila spermatogenesis involves a concerted developmentally-controlled program of specialized mitochondrial morphogenesis and thus is an ideal system for exploring molecular underpinnings. Mitochondria aggregate during and after male meiosis, subsequently fusing dramatically into two giant mitochondrial derivatives that interdigitate in layers at the early round spermatid stage to generate a large, onion-like, layered spherical structure called the Nebenkern [5–7]. During subsequent spermatid elongation stages, the mitochondrial derivatives within the Nebenkern unfurl and elongate beside the growing flagellar axoneme. Testis-enriched paralogs of more broadly-expressed genes  are common and are often identified via genetic analysis as associated with these functions, since null mutations in the testis-enriched paralogs can result in male-sterile but adult-viable mutant phenotypes.
ATP synthase is a large multi-protein complex, found in prokaryotes as well as in mitochondria and chloroplasts of eukaryotes, that catalyzes ATP formation from ADP and inorganic phosphate during cellular respiration and photosynthesis (reviewed in ). In mitochondria, the ATP synthase complex typically contains at least 16 subunits divided among the F1, F0, and peripheral stalk regions. The F1 portion extends into the mitochondrial matrix, physically attaching in two ways to the F0 portion, which is embedded in the inner membrane: via the F1 central stalk that allows rotation, and via the peripheral or stator stalk, which stabilizes the complex. The F0 portion rotates in response to proton movement down the gradient established by the electron transport chain, and the resulting torque and conformational changes in the F1 portion lead to ATP formation. Defects in the ATP synthase complex in humans can lead to neuromuscular disorders (reviewed in ) due to deficiencies in ATP production and overall mitochondrial dysfunction.
As shown mainly in yeast and cultured mammalian cells, the positive membrane curvature at the tips of cristae in the inner mitochondrial membrane is mediated in part by higher-order dimerization and oligomerization of ATP synthase complexes (reviewed in [11–14]). These interactions are separable from the metabolic role of ATP synthase, as some mutants with intact single ATP synthase complexes that cannot dimerize have aberrant cristae structure but can still respire [15, 16]. Cryo-electron microscopy indicates that bovine and yeast ATP synthase complex monomers induce 43° membrane curvature, and that the predicted dimer-induced membrane curvature of 86° matches observed values well [17, 18]. Tomography studies indeed show dimer ribbons corresponding with tight membrane curvature .
A logical question is whether regulation of ATP synthase higher-order structures helps to mediate variability in inner mitochondrial membrane conformation across cell types and different conditions. In Drosophila spermatogenesis, the large mitochondrial derivatives of the Nebenkern in early round spermatids do not have prominent cristae or sharp positive curvature of the inner mitochondrial membrane; instead the inner and outer mitochondrial membranes appear closely apposed across most of their area and parallel to each other in broad concentric circles of membranes reminiscent of an onion slice (hence early round spermatids are considered to be at the onion stage) [5, 20]. In this way, spermatid wild-type inner mitochondrial membranes resemble yeast inner mitochondrial membranes in cells deficient for ATP synthase dimerization.
Here we report a role for a tissue-specific isoform of ATP synthase subunit d in mediating internal Nebenkern structure within Drosophila melanogaster spermatids. We also provide evidence that other tissue-specific ATP synthase subunit paralogs help shape spermatid mitochondria. Subunit d is a major part of the peripheral stalk that stabilizes the F0 and F1 portions of ATP synthase as they rotate with respect to each other. Subunit d was first described in mammalian mitochondrial ATP synthase  and later in yeast . It is not found in bacterial or chloroplast ATP synthase. In S. cerevisiae, subunit d is oriented with its amino terminus farthest from the membrane and the carboxy terminus oriented closer to the inner mitochondrial membrane, spanning the distance between the F0 and F1 portions; while it does not itself span the inner mitochondrial membrane, it interacts via multiple coiled-coil domains with the transmembrane subunit b of the peripheral stalk [23, 24]. We identified an uncharacteristically large ATP synthase subunit d paralog expressed in Drosophila testes that is associated with internal Nebenkern structure and that may act via alteration of ATP synthase dimerization. Our work establishes a connection between tissue-specific variants of ATP synthase subunits and tissue-specific internal mitochondrial structure.
knon mutants show aberrant mitochondrial elongation during spermatogenesis and have faulty Nebenkern internal structure
knon is CG7813
We originally mapped knon by recombination to 2.1 ± 0.7 cM distal of curved (polytene region 52D2-9 on chromosome 2R), based on 38 recombinants between curved and welt. A w + -marked P element on the knon ms(2)1400 chromosome mapped near black, approximately 30 cM away, confirming its lack of association with the spermatogenesis phenotype. Further recombination mapping placed knon 0.67 ± 0.07 cM distal to Khc k13219 and 0.36 ± 0.07 cM proximal to Cdk4 k06503 , suggesting polytene interval 53B-C as the region of interest. Df(2R)7E (generated through imprecise excision of veg k03402 ), which lacks ~134 kb distal of the veg k03402 insertion site, failed to complement knon ms(2)1400 , showing a phenotype indistinguishable from knon ms(2)1400 homozygotes. In contrast, Df(2R)BSC550  complemented knon ms(2)1400 . Together, the breakpoints of Df(2R)7E and Df(2R)BSC550 defined the knon region as including two genes, Sema2a and CG7813 . Sema2a is expressed in many tissues  and known to function in the nervous system . FlyAtlas data indicate that CG7813 is expressed highly in the testis and at near-zero levels in adult carcass (minus the germline) and all other tissues tested except larval fat body, in which the larval gonad is embedded . modENCODE RNAseq data indicate that CG7813 is expressed highly in wandering third instar larvae and adult males but not at all in adult females , consistent with expression only in the male reproductive tract. Together these data suggested CG7813 as the best candidate for knon.
The MiMIC (Minos Mediated Integration Cassette) insertion line Mi[MIC]CG7813 Mi10492 [33, 34] is a hypomorphic allele of knon (CG7813) and is henceforth referred to as knon Mi10492 . The knon Mi10492 insertion is in the coding region of CG7813 between codons 554 and 555  and is therefore predicted to truncate the last 180 amino acids of Knon. Since the element is in reverse orientation with respect to the knon coding sequence, it does not provide a GFP tag to the truncated protein. Homozygous knon Mi10492 males were fertile but at slightly reduced levels, giving rise to 85% as many progeny as wild-type males (n = 56; P = 1.29x10-6), and show incompletely penetrant mild mitochondrial elongation defects (not shown). Males with knon Mi10492 in trans to either Df(2R)7E or knon ms(2)1400 were completely sterile with a mitochondrial elongation defect in spermatids similar to that in knon ms(2)1400 homozygotes, further confirming that the original knon ms(2)1400 phenotype is due to a lesion in CG7813.
Knon is a testis-enriched ATPsynD paralog with large C-terminal extension and associates with spermatid mitochondria
The knon gene encodes a 734 amino acid protein with an amino terminal region 32% identical and 50% similar to Drosophila ATPsynD (178 amino acids), the broadly-expressed canonical d subunit of the ATP synthase complex. Subunit d is part of the peripheral (or stator) stalk connecting the F0 and F1 portions of the complex. Knon, a testis-enriched paralog of subunit d, has 532 additional amino acids at its carboxy terminus beyond the region of homology. BLASTp and tBLASTn searches revealed that this C-terminal region was not homologous to other characterized gene products . Curiously, in this C-terminal region 28% of the residues are either glutamic acid or lysine, making it highly charged, and there are no predicted transmembrane domains or other motifs. The IUPred tool  predicts that most of this region is intrinsically unstructured.
Knon-GFP expressed under the control of endogenous knon regulatory sequences rescued the mutant phenotype (confirming functionality and relevance of localization) and was associated with mitochondria during later stages of spermatogenesis (Fig. 3). Knon-GFP was not detectable in the apical testis, which contains the stem cell populations and mitotically proliferating spermatogonia. The protein was detectable starting in spermatocytes during chromosome condensation. The timing of transition from absence of Knon-GFP to abundance was consistent with established patterns of transcription in primary spermatocytes . Knon-GFP co-localized with phase-dark mitochondria throughout spermatogenesis, associating with meiotic spindles, post-meiotic aggregating mitochondria, and Nebenkerne in the onion stage. Knon-GFP persisted on mitochondrial derivatives during and after mitochondrial elongation in elongated spermatid cysts. Mature, individualized wild-type sperm also contained the protein.
Exogenous expression of Knon in flight muscle alters the ratio of ATP synthase complex dimers to monomers but does not affect gross mitochondrial shaping or function
To ask whether the Knon-mediated shift toward more ATP synthase monomers was associated with altered flight muscle mitochondrial morphology, we examined by confocal microscopy tissue from flies expressing Knon in flight muscle along with a separate mitochondrially-targed GFP (mitoGFP)  (Fig. 4c-d). The mitoGFP localization was roughly similar between flies with and without exogenously expressed Knon, suggesting no gross alteration in mitochondrial number or size. Consistent with the lack of obvious effect on mitochondrial structure, flies expressing Knon in flight muscle were viable and able to fly, though there may be subtle effects at the levels of tissue architecture and flight ability that our experiments did not detect.
Other testis-enriched paralogs of ATP synthase subunits diverged at different times within insect lineages and have roles in spermatid mitochondrial dynamics
In addition to subunit d, the D. melanogaster genome includes pairs of paralogs encoding ATP synthase subunits F6, b, and g (which are in the peripheral stalk), as well as ε and β (in the F1 portion of the complex) . Expression data from modENCODE  and FlyAtlas  consistently indicate strong testis enrichment for one gene within each paralogous pair. Similar to subunit d, for subunit F6 the testis-enriched predicted gene product (147 amino acids) is significantly larger than that of the broadly-expressed paralog (106 amino acids). Conversely, for subunit b the testis-enriched predicted gene product (152 amino acids) is significantly smaller than that of the broadly expressed paralog (243 amino acids).
The Knon ATP synthase subunit d paralog is required for normal Nebenkern morphology
We found a role for Knon, an uncharacteristically large paralog of ATP synthase subunit d, in internal shaping of mitochondrial membranes in the Nebenkern during Drosophila spermatogenesis. Male flies lacking functional Knon were sterile and showed unusual membrane curvature and faulty internal structure within the Nebenkern (Fig. 1), with subsequent entanglement of the two unfurling mitochondrial derivatives and aberrant mitochondrial elongation as the flagellar axoneme grows. Simultaneously preventing mitochondrial fusion in a knon; fzo double mutant suppresses the knon single mutant elongation defect. This result indicates that the knon phenotype is due to the formation of an aberrant Nebernkern structure that cannot elongate, rather than defects at later stages. Knon is expressed only in the testis , and the protein is detected on mitochondria from meiosis through spermatid elongation (Fig. 3), consistent with the timing of dramatic mitochondrial membrane reorganization.
The effect of knon on mitochondrial shaping is likely at the level of mitochondrial dynamics and not resulting from respiratory deficiency, since 1) the mitochondrial membrane potential remains high in knon mutant spermatids; 2) other energy-requiring events of spermatid development, such as axoneme elongation and nuclear shaping, proceed; 3) disrupting ATP synthesis by knocking down the sole ATP synthase α subunit ortholog in the testis leads to tissue degeneration; and 4) in other systems, the effects of ATP synthase on organelle morphology and respiration are separable [15, 16].
Effects of Knon on ATP synthase dimerization
In yeast and mammalian cells, ATP synthase subunits e, g, and b are known to mediate dimerization of ATP synthase complexes and contribute to the sharp curvature of inner mitochondrial membranes at cristae tips [43, 44]. S. cerevisiae cells mutant for subunit e or subunit g, normally found in the membrane-embedded F0 portion of the complex, lack higher order ATP synthase dimers and have onion-like mitochondrial structures, reminiscent of the fly Nebenkern, that lack the tight inner membrane curvature normally found in cristae [45–47]. Subunit b in the peripheral stalk is also important for dimerization of the ATP synthase complex, and it can homodimerize on its own [48, 49]. In yeast with a truncated subunit b missing the transmembrane portion, no ATP synthase dimers form, and inner mitochondrial membranes lack cristae .
Peripheral stalk components farther away from the membrane also contribute to ATP synthase dimerization. Subunit h (F6 in metazoans), in the peripheral stalk, can homodimerize and stabilize dimers in yeast even if subunits e and g are missing [50, 51]. Cryo-electron microscopy of ATP synthase dimers from the alga Polytomella indicates that peripheral stalk components from adjacent ATP synthase monomers intertwine at a position 80 angstroms away from the inner mitochondrial membrane , consistent with the positions of subunits d and h/F6. The broadly-expressed D. melanogaster subunit d (CG6030) has 178 amino acids and is expressed in all tissues, although at lower levels in testes . Orthologs in humans, mice, C. elegans, and budding yeast are similar in size, with 161, 161, 191, and 174 amino acids, respectively. In contrast, Knon has 734 amino acids as a result of a large C-terminal extension. The crystal structure of bovine subunit d suggests that its carboxy terminus extends from the middle of the peripheral stalk toward the membrane and embedded F0 portion , consistent with the idea that Knon’s extra large domain in that position may alter the interactions between ATP synthase complex monomers.
Exogenous expression of Knon in flight muscle indeed reduces dimerization of ATP synthase (Fig. 4). Expression of Knon in the absence of testis-specific versions of subunits F6, g, and b may not, however, allow for assembly dynamics of a full testis-specific ATP synthase complex. We therefore cannot rule out the possibility that in wild-type spermatids, Knon along with the other testis-specific paralogs may help mediate ATP synthase dimerization at a shallower angle. Dimerization of ATP synthase complexes at a shallower angle has been seen in Tetrahymena . ATP synthase dimerization in the testis is difficult to assess directly because of tissue heterogeneity and technical constraints.
Truncated KnonMi10492 lacking the carboxy-terminal 180 amino acids is partially functional, suggesting that 554 amino acids (plus any encoded on the inserted element before the first stop codon is encountered), are sufficient for Knon to perform its function when threshold levels of protein are present. These data are consistent with the carboxy-terminal extension either providing steric hindrance that blocks dimerization or promoting unique interactions that underlie dimerization at a different angle.
Roles for other testis-enriched ATP synthase paralogs in spermatid mitochondrial shaping, and implications of different divergence times
Three other subunits of ATP synthase in the peripheral stalk and/or at the putative dimerization interface (F6, g, and b), and two subunits in the F1 portion (β and ε), appear to be encoded by broadly-expressed and testis-specific paralogs in the Drosophila melanogaster genome. We showed via RNAi knockdown experiments that the testis-specific paralogs of subunits F6 and g are required for proper Nebenkern dynamics during spermatogenesis, and that these phenotypes differed from the outcome of disrupting ATP synthesis via knockdown of the sole α subunit ortholog. The testis-enriched F6 subunit, like the Knon subunit d, is larger than the broadly-expressed version, while the two g subunit paralogs are similar in size. The testis-enriched g subunit therefore might not provide a steric hindrance to dimerization of the complex but might instead promote assembly of the testis-enriched large d and F6 subunits into the complex, or promote interactions between adjacent ATP synthase complexes at a shallower angle at the base of the peripheral stalk.
For subunit b, the RNAi transgenic line targeting the testis-enriched paralog may not be sufficiently inducible to show an effect, or alternatively the testis-enriched subunit b might truly be not required for spermatid morphogenesis. Consistent with our data, another group recently showed that knockdown of the broadly-expressed ATPsynB results in male infertility and a lack of mature sperm, though they did not examine the mitochondrial phenotype . Dispensability of the testis-specific subunit b would be consistent with its relatively recent gene duplication (Fig. 5), such that there has been less time for specialization of the two paralogs compared to the other paralogous pairs.
The mechanism of differential gene regulation between the paralogs of peripheral stalk ATP synthase subunits is unknown. Notably, the paralogous ATP synthase β subunits were duplicated earlier than the others, and a recent study showed that the broadly-expressed Drosophila ATPsynβ is normally repressed in germ cells by the hairpin RNA pathway . Derepression of ATPsynβ via knockout of the appropriate hpRNA causes defects in germ cell differentiation and male sterility . These findings suggest that the ATPsynβ paralogs have functionally diverged and have evolved RNA-based regulatory mechanisms to favor the expression of one paralog over the other in the testis. The mechanisms of regulation of the other paralogous pairs are yet to be determined.
ATP synthase, mitochondrial shaping, and tissue differentiation
In mammalian systems the ultrastructural organization of the inner mitochondrial membrane can vary among tissue types and in response to metabolic state or developmental cues (reviewed in ). Expression of subunit e, important for dimerization, can be triggered by altered nutrient and oxygen levels [56, 57]. In differentiating cardiomyocytes, mitochondria undergo structural changes concomitant with ATP synthase dimer stabilization by an accessory factor [58, 59]. In the syncytiotrophoblast outer layer of the placenta, reduced ATP synthase dimerization accompanies a lack of cristae, despite normal ATP production . Specific mutations in the mitochondrially-encoded ATP6 subunit found in Leigh syndrome patients are associated with higher levels of ATP synthase dimerization and decreased turnover of the complex .
The Drosophila ovary requires broadly-expressed ATP synthase components for germ cell differentiation , possibly associated with mitochondrial morphology and not respiration, because knockdowns of other respiratory complexes did not show the same effect. Germline-specific ATP synthase subunits in C. elegans appear required for tissue differentiation , though no mitochondrial structural aspects were examined. In trypanosomes, different ATP synthase subunit c paralogs are encoded in the genome; it is not yet known whether these gene products influence mitochondrial function and morphology in different life cycle stages , though subunit c is not positioned to influence dimerization of the complex directly. Our findings in the Drosophila testis therefore elucidate the first known instance in which tissue-specific ATP synthase subunits underlie tissue-specific regulation of mitochondrial shaping.
We describe the first evidence connecting tissue-specific mitochondrial shaping with the expression and function of tissue-specific ATP synthase subunits. A testis-specific variant of ATP synthase subunit d in Drosophila is required for Nebenkern morphology and appears to alter dimerization of ATP synthase complexes. Testis-specific paralogs of other ATP synthase subunits also appear to have a role. Our data, showing the role of a cell biological process in the context of tissue morphogenesis, contribute to the growing knowledge of ATP synthase’s role in mediating membrane conformation, beyond its canonical role in ATP production. Future work with tissue-specific ATP synthase subunit paralogs in Drosophila may explore their sub-mitochondrial localization to test localization in areas of shallower inner mitochondrial membrane curvature; further work should also explore whether the testis-specific paralogs together form a unique complex, and whether individual subunits can be substituted between the different forms of the complex.
Fly stocks, husbandry, and fertility tests
Flies were raised on standard cornmeal molasses medium or instant potato flake medium (Ward’s Scientific) at 25 °C for all crosses except for those combining GAL4 and UAS transgenes, which were maintained at 29 °C. Multiply marked stocks, transposable element insertions, the Δ2-3 transposase-encoding stock, the Act88F-GAL4 driver stock, and all deficiencies except Df(2R)7E were obtained from the Bloomington Drosophila Stock Center. Oregon R was used as the wild-type strain. Flies with transgenic UAS-RNAi constructs were obtained from the Vienna Drosophila Resource Center. The bam-GAL4 driver stock was a gift from Yukiko Yamashita. Most fertility tests were performed by mating individual males with three Oregon R virgin females and assessing the presence of larvae after five days. For quantifying fertility of knon Mi10492 allelic combinations, 56 age-matched males of each genotype were placed individually in vials with three age-matched Oregon R virgin females; pupae and eclosed offspring from vials in which the males survived were counted after fifteen days. Data were analyzed by T test.
P[GD12291]v22566 (knon; CG7813; testis-specific subunit d)
P[KK1096241]v104818 (CG12027; testis-enriched subunit F6)
P[KK108581]v107826 (CG4412; broadly-expressed subunit F6)
P[GD12493]v35387 (CG7211; testis-enriched subunit g)
P[KK108323]v107311 (CG6105; broadly-expressed subunit g)
P[GD1688]v6521 (CG17300; testis-enriched subunit b)
P[GD6129]v14211 (CG8189; broadly-expressed subunit b)
P[GD11030]v34663 (CG3612 (bellwether); sole ortholog of subunit α)
For knockdown in the testis, each RNAi line was crossed to a bam-GAL4 driver line . Testes of F1 males were dissected and imaged as described below. For RNAi knockdown in flight muscle, each relevant UAS-RNAi line was crossed to the Act88F-GAL4 third chromosome driver line (Bloomington Drosophila Stock Center).
For light microscopy, testes were dissected in TB1 buffer (7 mM K2HPO4, 7 mM KH2PO4 (pH 6.7), 80 mM KCl, 16 mM NaCl, 5 mM MgCl2, 1% PEG-6000) and opened with forceps to allow cells to form a monolayer under pressure of a cover slip. Co-labeling of DNA and polarized mitochondria was performed by incubating intact testes in TB1 supplemented with 5% DMSO, 2.5 μg/mL Hoechst, and 10 μg/mL rhodamine 123  for ~10-15 min in the dark. The testes were washed three times with TB1 to remove residual stain and then torn open on a polylysine-coated slide prior to imaging. Phase-contrast and epifluorescence microscopy were performed on a Nikon Eclipse E600W, Olympus B201, or Olympus BX60 microscope, and images were captured with a Nikon Coolpix 4500 or Olympus DP80 camera.
For confocal microscopy, thoraces of adult males and females were dissected to obtain flight muscle tissue. Thoraces were fixed in 4% paraformaldehyde with PBS, washed 3x with PBT (0.2% TritonX), and stained overnight with tetramethylrhodamine-conjugated phalloidin in PBT at a concentration of 8 μM. Whole thoraces were mounted on a slide with mounting media and sealed with a coverslip and clear nail polish. The flight muscle extruded from the exoskeleton after pressing and sealing the coverslip. Samples were observed with a Nikon 90i confocal microscope and a Plan Apo TIRF 100X/1.45 oil objective. Images were captured and overlaid using the confocal software EZ-C1.
For transmission electron microscopy, testes were dissected in TB1 and immediately placed in fixative (2% glutaraldehyde, 1% paraformaldehyde, 0.1 M sodium phosphate or sodium cacodylate buffer (pH 7)). After overnight fixation, samples were washed in 0.1% phosphate or cacodylate buffer for 15 min and stained with 1% osmium tetroxide in the same buffer for two hours. Testes were then washed three times in water, stained for 1 h in 1% uranyl acetate, washed three times in water and dehydrated through an ethanol series (30%, 50%, 70%, 95%, 100%). After five minutes in 1:1 ethanol:propylene oxide and five minutes in propylene oxide, samples were embedded in Spurr’s resin and polymerized overnight at 60 °C. Thin sections (80-90 nm) were cut with a Reichert-Jung microtome, placed on Formvar-coated slot grids and examined on a Phillips 410 transmission electron microscope.
Recombination and deficiency mapping
Virgin females of genotype knon/ al dp b r cn vg c a px bw mi sp were crossed to al dp b pr Bl c px sp/CyO males. F1 males that carried Bl along with other markers indicating single recombination events on the second chromosome were selected for each interval and crossed individually to w; knon/CyO females to score both knon and the w +-marked P element, and to make stocks of the recombinant chromosome with CyO. For finer mapping, knon/c wt px females were crossed to al dp b pr Bl c px sp/CyO males, and F1 male recombinant progeny in the c-px interval were individually crossed to knon/CyO to score knon and make stocks of recombinant chromosomes. The welt (wt) marker was scored in homozygotes in each of the recombinant stocks. For deficiency mapping, knon/CyO flies were crossed to selected 2nd chromosome deficiency lines balanced over CyO, and the Cy+ male offspring were tested for fertility and dissected as described above to assess testis subcellular phenotype.
Generation of transgenic lines
A knon transgene with endogenous regulatory sequences was generated by amplifying a ~5 kb genomic region containing CG7813 along with ~1.8 kb upstream and ~0.5 kb downstream, with primers 5’ ATATATGGTACCGCCATCAGCCAGCGAGCTT 3’ (forward) and 5’ GCATACCGCGGGGCAACACTTTTCGAACGG 3’ (reverse), and cloning the fragment into the pCaSpeR4 vector  using restriction sites KpnI and SacII. Primers were purchased from Integrated DNA Technologies (Coralville, IA). A separate knon-GFP transgene was created to include the ~1.8 kb upstream region and knon coding region through the penultimate codon, with GFP (GFPmut3* variant with S2R, S65G, S72A ) ligated in-frame 3’ to the coding region at the SacII restriction site. To create UAS-knon and UAS-knon-GFP constructs, the knon and knon-GFP constructs described above were used as PCR template in two-step splicing by overlap extension. The single knon intron was removed by amplifying the two exons of the gene in separate PCR reactions, using Phusion High-Fidelity PCR Master Mix with HF Buffer (New England Biolabs Mo531S). The two products were then pooled in a third reaction, where they annealed and amplified to form a product with sequence identical to the cDNA. The inserts were cloned into pUAST (; obtained from the Drosophila Genomics Resource Center, Bloomington, IN) using KpnI and XbaI (New England Biolabs). DNA sequencing to check constructs was performed by Retrogen (San Diego, CA), and DNA sequence analysis was performed in ApE . Constructs were injected into w 1118 embryos by Rainbow Transgenics (Camarillo, CA). w + F1 offspring were outcrossed to Kr/CyO; D/TM6C to make stocks and map insertions. knon and knon-GFP transgenes inserted on the third chromosome were crossed into the knon mutant background.
BN-PAGE and immunoblotting of flight muscle mitochondria
Mitochondria were extracted from fly thoraces using a protocol modified from . Fifteen thoraces from flies of each genotype were crushed in 50 μL mitochondrial isolation medium (MIM: 250 mM sucrose, 10 mM Tris-HCl, pH 7.4, 0.15 mM MgCl2) using a plastic pestle, then spun twice at 500 x g for 5 min at 4 °C to remove debris. The supernatant was spun at 5,000 x g for 5 min at 4 °C to pellet mitochondria. The pellet was subjected to the Life Technologies NativePAGE Novex® Bis-Tris Gel System protein isolation protocol for BN-PAGE. Protein pellets were resuspended in 25 μL of 1x Native PAGE sample buffer (Invitrogen) containing 1% digitonin. Samples were incubated on ice for 15 min and centrifuged at 20,000 x g for 30 min at 4 °C. The supernatant was frozen at -80 °C. Using reagents and protocol from Life Technologies, NativePAGE Novex® Bis-Tris Gel System, G-250 dye was added to samples to a final concentration of 0.25%, and samples were run on a 3-12% polyacrylamide gel for 90 min at 150 V constant in the X-Cell Sure Lock Mini Cell gel apparatus by Life Technologies. Immunoblotting was performed with the X Cell II Blot Module (Invitrogen). The gel was soaked for 15 min in a 0.1% SDS solution with water and samples transferred to PVDF membrane in the X Cell transfer apparatus on ice for 1 h. Following transfer to the membrane, the membrane was fixed in 8% acetic acid and rinsed with methanol and dH2O. The primary antibody was mouse monoclonal anti-ATP5A antibody [15H4C4] by Abcam at a dilution of 1:3000 or 1:5000. The blocking buffer was StartingBlock (PBS) Blocking Buffer from Thermo Scientific with 0.1% Tween. The secondary antibody was goat anti-mouse IgG-HRP by Santa Cruz Biotechnology [sc-2005] at a dilution of 1:2500. Primary incubation occurred overnight, shaking at 4 °C and secondary incubation occurred at RT, shaking for 1 h. Antibodies were diluted with blocking buffer without Tween. Membranes were imaged and band intensity quantified using a Bio-Rad gel imager for chemiluminescence and ImageJ software. To determine the internal ratio of the dimer to monomer bands within each given lane, background was subtracted and band intensity measured on non-saturated blot images.
ATP synthase subunit paralogs were identified through FlyBase annotations and orthology tools [27, 70] supplemented with NCBI BLAST  and VectorBase  to search for ATP synthase subunits encoded by multiple paralogs in genomes outside the Drosophila genus. We used maximum likelihood phylogeny  to test if the proteins encoded by all paralogs cluster into two distinct groups. Clustering, tree drawing, and evolution rate analyses were performed in the MEGA5 suite [74, 75].
We thank Patricia Wilson for the knon ms(2)1400 allele, Bee Choo Liang and Anthony Mahowald for help with initial mapping and Fran Thomas for assistance with electron microscopy. The Davidson College Fall 2009 Genetics class assisted with deficiency mapping of knon ms(2)1400 . Sarah Pyfrom and Ben Jepson provided technical assistance. We thank Barbara Lom for valuable guidance with the confocal microscope, and Jon Lim and Riley Mangan for assistance with immunoblots. Many thanks to three anonymous reviewers and Devon Harris for constructive suggestions on the manuscript. FlyBase, the Bloomington Drosophila Stock Center, and the Vienna Drosophila Resource Center were invaluable resources for data and fly stocks.
Work in the Hales lab was supported by National Science Foundation RUI 1158024 as well as by the Davidson College Department of Biology. YH, COF, and OM were supported by Davidson Research Initiative grants. In the mid-1990s, when this project was initiated, MTF was funded by National Institutes of Health (NIH) grant 1RO1-HD29194, and KGH was funded by a Howard Hughes Medical Institute predoctoral fellowship.
Availability of data and materials
The datasets used and/or analysed during the current study are available from the corresponding author on reasonable request.
MTF and KGH conceptualized, respectively, initial and later parts of this project. EMS, ECB, YH, LEI, OB, MB, DA, KES, OM, COF, NS, MGG, ETG, LAR, CGW, and KGH designed experiments and acquired, analyzed, and interpreted data. KGH drafted the manuscript with substantial contributions from EMS and ECB, and MTF also provided editorial contributions. This work was begun in the 1990s when KGH was a predoctoral student in the laboratory of MTF, and the bulk of the work was later completed after a long hiatus, in the laboratory of KGH. All authors read and approved the final manuscript.
The authors declare that they have no competing interests.
Consent for publication
Ethics approval and consent to participate
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Open AccessThis article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
- Friedman JR, Nunnari J. Mitochondrial form and function. Nature. 2014;505:335–43.View ArticlePubMedPubMed CentralGoogle Scholar
- Mishra P, Chan DC. Mitochondrial dynamics and inheritance during cell division, development and disease. Nat Rev Mol Cell Biol. 2014;15:634–46.View ArticlePubMedPubMed CentralGoogle Scholar
- Hales KG, Fuller MT. Developmentally regulated mitochondrial fusion mediated by a conserved, novel, predicted GTPase. Cell. 1997;90:121–9.View ArticlePubMedGoogle Scholar
- Stowers RS, Megeath LJ, Gorska-Andrzejak J, Meinertzhagen IA, Schwarz TL. Axonal transport of mitochondria to synapses depends on milton, a novel Drosophila protein. Neuron. 2002;36:1063–77.View ArticlePubMedGoogle Scholar
- Tates AD. Cytodifferentiation during spermatogenesis in Drosophila melanogaster: an electron microscope study. Leiden: Rijksuniversiteit; 1971.Google Scholar
- Fuller MT. Spermatogenesis. In: Bate M, Martinez-Arias A, editors. The Development of Drosophila melanogaster. Cold Spring Harbor:Cold Spring Harbor Press; 1993. p. 71–147.Google Scholar
- Fabian L, Brill JA. Drosophila spermiogenesis: big things come from little packages. Spermatogenesis. 2012;2:197–212.View ArticlePubMedPubMed CentralGoogle Scholar
- Assis R, Bachtrog D. Neofunctionalization of young duplicate genes in Drosophila. Proc Natl Acad Sci U S A. 2013;110:17409–14.View ArticlePubMedPubMed CentralGoogle Scholar
- Devenish RJ, Prescott M, Rodgers AJ. The structure and function of mitochondrial F1F0-ATP synthases. Int Rev Cell Mol Biol. 2008;267:1–58.View ArticlePubMedGoogle Scholar
- Kucharczyk R, Zick M, Bietenhader M, Rak M, Couplan E, Blondel M, Caubet SD, di Rago JP. Mitochondrial ATP synthase disorders: molecular mechanisms and the quest for curative therapeutic approaches. Biochim Biophys Acta. 2009;1793:186–99.View ArticlePubMedGoogle Scholar
- Habersetzer J, Ziani W, Larrieu I, Stines-Chaumeil C, Giraud MF, Brethes D, Dautant A, Paumard P. ATP synthase oligomerization: from the enzyme models to the mitochondrial morphology. Int J Biochem Cell Biol. 2013;45:99–105.View ArticlePubMedGoogle Scholar
- Seelert H, Dencher NA. ATP synthase superassemblies in animals and plants: two or more are better. Biochim Biophys Acta. 1807;2011:1185–97.Google Scholar
- Wittig I, Schagger H. Supramolecular organization of ATP synthase and respiratory chain in mitochondrial membranes. Biochim Biophys Acta. 2009;1787:672–80.View ArticlePubMedGoogle Scholar
- Zick M, Rabl R, Reichert AS. Cristae formation—linking ultrastructure and function of mitochondria. Biochim Biophys Acta. 2009;1793:5–19.Google Scholar
- Arnold I, Pfeiffer K, Neupert W, Stuart RA, Schagger H. Yeast mitochondrial F1F0-ATP synthase exists as a dimer: identification of three dimer-specific subunits. EMBO J. 1998;17:7170–8.View ArticlePubMedPubMed CentralGoogle Scholar
- Forster K, Turina P, Drepper F, Haehnel W, Fischer S, Graber P, Petersen J. Proton transport coupled ATP synthesis by the purified yeast H+ -ATP synthase in proteoliposomes. Biochim Biophys Acta. 2010;1797:1828–37.View ArticlePubMedGoogle Scholar
- Baker LA, Watt IN, Runswick MJ, Walker JE, Rubinstein JL. Arrangement of subunits in intact mammalian mitochondrial ATP synthase determined by cryo-EM. Proc Natl Acad Sci U S A. 2012;109:11675–80.View ArticlePubMedPubMed CentralGoogle Scholar
- Davies KM, Anselmi C, Wittig I, Faraldo-Gomez JD, Kuhlbrandt W. Structure of the yeast F1F0-ATP synthase dimer and its role in shaping the mitochondrial cristae. Proc Natl Acad Sci U S A. 2012;109:13602–7.View ArticlePubMedPubMed CentralGoogle Scholar
- Strauss M, Hofhaus G, Schroder RR, Kuhlbrandt W. Dimer ribbons of ATP synthase shape the inner mitochondrial membrane. EMBO J. 2008;27:1154–60.View ArticlePubMedPubMed CentralGoogle Scholar
- Tokuyasu KT. Dynamics of spermiogenesis in Drosophila melanogaster. VI. Significance of "onion" nebenkern formation. J Ultrastruct Res. 1975;53:93–112.View ArticlePubMedGoogle Scholar
- Walker JE, Runswick MJ, Poulter L. ATP synthase from bovine mitochondria. The characterization and sequence analysis of two membrane-associated sub-units and of the corresponding cDNAs. J Mol Biol. 1987;197:89–100.View ArticlePubMedGoogle Scholar
- Norais N, Prome D, Velours J. ATP synthase of yeast mitochondria. Characterization of subunit d and sequence analysis of the structural gene ATP7. J Biol Chem. 1991;266:16541–9.PubMedGoogle Scholar
- Bueler SA, Rubinstein JL. Location of subunit d in the peripheral stalk of the ATP synthase from Saccharomyces cerevisiae. Biochemistry. 2008;47:11804–10.View ArticlePubMedGoogle Scholar
- Dickson VK, Silvester JA, Fearnley IM, Leslie AG, Walker JE. On the structure of the stator of the mitochondrial ATP synthase. EMBO J. 2006;25:2911–8.View ArticlePubMedPubMed CentralGoogle Scholar
- Wilson PG, Fuller MT, Borisy GG. Monastral bipolar spindles: implications for dynamic centrosome organization. J Cell Sci. 1997;110(Pt 4):451–64.PubMedGoogle Scholar
- Cook RK, Christensen SJ, Deal JA, Coburn RA, Deal ME, Gresens JM, Kaufman TC, Cook KR. The generation of chromosomal deletions to provide extensive coverage and subdivision of the Drosophila melanogaster genome. Genome Biol. 2012;13:R21.View ArticlePubMedPubMed CentralGoogle Scholar
- dos Santos G, Schroeder AJ, Goodman JL, Strelets VB, Crosby MA, Thurmond J, Emmert DB, Gelbart WM, FlyBase C. FlyBase: introduction of the Drosophila melanogaster release 6 reference genome assembly and large-scale migration of genome annotations. Nucleic Acids Res. 2015;43:D690–697.View ArticlePubMedGoogle Scholar
- Chintapalli VR, Wang J, Dow JA. Using FlyAtlas to identify better Drosophila melanogaster models of human disease. Nat Genet. 2007;39:715–20.Google Scholar
- Kolodkin AL, Matthes DJ, Goodman CS. The semaphorin genes encode a family of transmembrane and secreted growth cone guidance molecules. Cell. 1993;75:1389–99.View ArticlePubMedGoogle Scholar
- modEncode Consortium, Roy S, Ernst J, Kharchenko PV, Kheradpour P, Negre N, Eaton ML, Landolin JM, Bristow CA, Ma L, et al. Identification of functional elements and regulatory circuits by Drosophila modENCODE. Science. 2010;330:1787–97.View ArticleGoogle Scholar
- Chen D, McKearin DM. A discrete transcriptional silencer in the bam gene determines asymmetric division of the Drosophila germline stem cell. Development. 2003;130:1159–70.View ArticlePubMedGoogle Scholar
- Dietzl G, Chen D, Schnorrer F, Su KC, Barinova Y, Fellner M, Gasser B, Kinsey K, Oppel S, Scheiblauer S, et al. A genome-wide transgenic RNAi library for conditional gene inactivation in Drosophila. Nature. 2007;448:151–6.View ArticlePubMedGoogle Scholar
- Nagarkar-Jaiswal S, Lee PT, Campbell ME, Chen K, Anguiano-Zarate S, Gutierrez MC, Busby T, Lin WW, He Y, Schulze KL et al. A library of MiMICs allows tagging of genes and reversible, spatial and temporal knockdown of proteins in Drosophila. eLife. 2015;4:e05338.Google Scholar
- Venken KJ, Schulze KL, Haelterman NA, Pan H, He Y, Evans-Holm M, Carlson JW, Levis RW, Spradling AC, Hoskins RA, et al. MiMIC: a highly versatile transposon insertion resource for engineering Drosophila melanogaster genes. Nature Methods. 2011;8:737–43.Google Scholar
- Altschul SF, Gish W, Miller W, Myers EW, Lipman DJ. Basic local alignment search tool. J Mol Biol. 1990;215:403–10.View ArticlePubMedGoogle Scholar
- Dosztanyi Z, Csizmok V, Tompa P, Simon I. IUPred: web server for the prediction of intrinsically unstructured regions of proteins based on estimated energy content. Bioinformatics. 2005;21:3433–4.View ArticlePubMedGoogle Scholar
- White-Cooper H, Schafer MA, Alphey LS, Fuller MT. Transcriptional and post-transcriptional control mechanisms coordinate the onset of spermatid differentiation with meiosis I in Drosophila. Development. 1998;125:125–34.PubMedGoogle Scholar
- Brand AH, Perrimon N. Targeted gene expression as a means of altering cell fates and generating dominant phenotypes. Development. 1993;118:401–15.PubMedGoogle Scholar
- Rizzuto R, Brini M, Pizzo P, Murgia M, Pozzan T. Chimeric green fluorescent protein as a tool for visualizing subcellular organelles in living cells. Curr Biol. 1995;5:635–42.View ArticlePubMedGoogle Scholar
- Diaz-Castillo C, Ranz JM. Nuclear chromosome dynamics in the Drosophila male germ line contribute to the nonrandom genomic distribution of retrogenes. Mol Biol Evol. 2012;29:2105–8.View ArticlePubMedGoogle Scholar
- Zhang C, Gschwend AR, Ouyang Y, Long M. Evolution of gene structural complexity: an alternative-splicing-based model accounts for intron-containing retrogenes. Plant Physiol. 2014;165:412–23.View ArticlePubMedPubMed CentralGoogle Scholar
- Jacobs H, Stratmann R, Lehner C. A screen for lethal mutations in the chromosomal region 59AB suggests that bellwether encodes the alpha subunit of the mitochondrial ATP synthase in Drosophila melanogaster. Mol Gen Genet. 1998;259:383–7.View ArticlePubMedGoogle Scholar
- Arselin G, Vaillier J, Salin B, Schaeffer J, Giraud MF, Dautant A, Brethes D, Velours J. The modulation in subunits e and g amounts of yeast ATP synthase modifies mitochondrial cristae morphology. J Biol Chem. 2004;279:40392–9.View ArticlePubMedGoogle Scholar
- Soubannier V, Vaillier J, Paumard P, Coulary B, Schaeffer J, Velours J. In the absence of the first membrane-spanning segment of subunit 4(b), the yeast ATP synthase is functional but does not dimerize or oligomerize. J Biol Chem. 2002;277:10739–45.View ArticlePubMedGoogle Scholar
- Giraud MF, Paumard P, Soubannier V, Vaillier J, Arselin G, Salin B, Schaeffer J, Brethes D, di Rago JP, Velours J. Is there a relationship between the supramolecular organization of the mitochondrial ATP synthase and the formation of cristae? Biochim Biophys Acta. 2002;1555:174–80.View ArticlePubMedGoogle Scholar
- Paumard P, Vaillier J, Coulary B, Schaeffer J, Soubannier V, Mueller DM, Brethes D, di Rago JP, Velours J. The ATP synthase is involved in generating mitochondrial cristae morphology. EMBO J. 2002;21:221–30.View ArticlePubMedPubMed CentralGoogle Scholar
- Velours J, Dautant A, Salin B, Sagot I, Brethes D. Mitochondrial F1F0-ATP synthase and organellar internal architecture. Int J Biochem Cell Biol. 2009;41:1783–9.View ArticlePubMedGoogle Scholar
- Spannagel C, Vaillier J, Arselin G, Graves PV, Grandier-Vazeille X, Velours J. Evidence of a subunit 4 (subunit b) dimer in favor of the proximity of ATP synthase complexes in yeast inner mitochondrial membrane. Biochim Biophys Acta. 1998;1414:260–4.View ArticlePubMedGoogle Scholar
- Gavin PD, Prescott M, Devenish RJ. F1F0-ATP synthase complex interactions in vivo can occur in the absence of the dimer specific subunit e. J Bioenerg Biomembr. 2005;37:55–66.View ArticlePubMedGoogle Scholar
- Fronzes R, Weimann T, Vaillier J, Velours J, Brethes D. The peripheral stalk participates in the yeast ATP synthase dimerization independently of e and g subunits. Biochemistry. 2006;45:6715–23.View ArticlePubMedGoogle Scholar
- Fronzes R, Chaignepain S, Bathany K, Giraud MF, Arselin G, Schmitter JM, Dautant A, Velours J, Brethes D. Topological and functional study of subunit h of the F1F0-ATP synthase complex in yeast Saccharomyces cerevisiae. Biochemistry. 2003;42:12038–49.View ArticlePubMedGoogle Scholar
- Allegretti M, Klusch N, Mills DJ, Vonck J, Kuhlbrandt W, Davies KM. Horizontal membrane-intrinsic alpha-helices in the stator a-subunit of an F-type ATP synthase. Nature. 2015;521:237–40.View ArticlePubMedGoogle Scholar
- Balabaskaran Nina P, Dudkina NV, Kane LA, van Eyk JE, Boekema EJ, Mather MW, Vaidya AB. Highly divergent mitochondrial ATP synthase complexes in Tetrahymena thermophila. PLoS Biol. 2010;8:e1000418.View ArticlePubMedPubMed CentralGoogle Scholar
- Chen YN, Wu CH, Zheng Y, Li JJ, Wang JL, Wang YF. Knockdown of ATPsyn-b caused larval growth defect and male infertility in Drosophila. Archive Insect Biochem Physiol. 2015;88:144–54.Google Scholar
- Wen J, Duan H, Bejarano F, Okamura K, Fabian L, Brill JA, Bortolamiol-Becet D, Martin R, Ruby JG, Lai EC. Adaptive regulation of testis gene expression and control of male fertility by the Drosophila harpin RNA pathway. Molecular Cell. 2015;57:165–78.Google Scholar
- Levy FH, Kelly DP. Regulation of ATP synthase subunit e gene expression by hypoxia: cell differentiation stage-specific control. Am J Physiol. 1997;272:C457–465.Google Scholar
- Swartz DA, Park EI, Visek WJ, Kaput J. The e subunit gene of murine F1F0-ATP synthase. Genomic sequence, chromosomal mapping, and diet regulation. J Biol Chem. 1996;271:20942–8.View ArticlePubMedGoogle Scholar
- Bisetto E, Comelli M, Salzano AM, Picotti P, Scaloni A, Lippe G, Mavelli I. Proteomic analysis of F1F0-ATP synthase super-assembly in mitochondria of cardiomyoblasts undergoing differentiation to the cardiac lineage. Biochim Biophys Acta. 1827;2013:807–16.Google Scholar
- Campanella M, Casswell E, Chong S, Farah Z, Wieckowski MR, Abramov AY, Tinker A, Duchen MR. Regulation of mitochondrial structure and function by the F1F0-ATPase inhibitor protein, IF1. Cell Metab. 2008;8:13–25.View ArticlePubMedGoogle Scholar
- De los Rios Castillo D, Zarco-Zavala M, Olvera-Sanchez S, Pardo JP, Juarez O, Martinez F, Mendoza-Hernandez G, Garcia-Trejo JJ, Flores-Herrera O. Atypical cristae morphology of human syncytiotrophoblast mitochondria: role for complex V. J Biol Chem. 2011;286:23911–9.View ArticlePubMedPubMed CentralGoogle Scholar
- Cortes-Hernandez P, Vazquez-Memije ME, Garcia JJ. ATP6 homoplasmic mutations inhibit and destabilize the human F1F0-ATP synthase without preventing enzyme assembly and oligomerization. J Biol Chem. 2007;282:1051–8.View ArticlePubMedGoogle Scholar
- Teixeira FK, Sanchez CG, Hurd TR, Seifert JR, Czech B, Preall JB, Hannon GJ, Lehmann R. ATP synthase promotes germ cell differentiation independent of oxidative phosphorylation. Nat Cell Biol. 2015;17:689–96.View ArticlePubMedPubMed CentralGoogle Scholar
- Kawasaki I, Hanazawa M, Gengyo-Ando K, Mitani S, Maruyama I, Iino Y. ASB-1, a germline-specific isoform of mitochondrial ATP synthase b subunit, is required to maintain the rate of germline development in Caenorhabditis elegans. Mech Dev. 2007;124:237–51.Google Scholar
- Gulde PE, Christen L, Brown SV, Williams N. Three distinct isoforms of ATP synthase subunit c are expressed in T. brucei and assembled into the mitochondrial ATP synthase complex. PLoS One. 2013;8:e54039.View ArticlePubMedPubMed CentralGoogle Scholar
- Johnson LV, Walsh ML, Chen LB. Localization of mitochondria in living cells with rhodamine 123. Proc Natl Acad Sci U S A. 1980;77:990–4.View ArticlePubMedPubMed CentralGoogle Scholar
- Thummel C, Pirrotta V. Technical notes: new pCaSpeR P element vectors. Dros Inf Serv. 1992;71:150.Google Scholar
- Cormack BP, Valdivia RH, Falkow S. FACS-optimized mutants of the green fluorescent protein (GFP). Gene. 1996;173:33–8.View ArticlePubMedGoogle Scholar
- Davis MW. ApE: a plasmid editor. http://biologylabs.utah.edu/jorgensen/wayned/ape/. Accessed 27 Aug 2012.
- Rera M, Bahadorani S, Cho J, Koehler CL, Ulgherait M, Hur JH, Ansari WS, Lo Jr T, Jones DL, Walker DW. Modulation of longevity and tissue homeostasis by the Drosophila PGC-1 homolog. Cell Metab. 2011;14:623–34.View ArticlePubMedPubMed CentralGoogle Scholar
- Drosophila 12 Genomes Consortium, Clark AG, Eisen MB, Smith DR, Bergman CM, Oliver B, Markow TA, Kaufman TC, Kellis M, Gelbart W, et al. Evolution of genes and genomes on the Drosophila phylogeny. Nature. 2007;450:203–18.Google Scholar
- Johnson M, Zaretskaya I, Raytselis Y, Merezhuk Y, McGinnis S, Madden TL. NCBI BLAST: a better web interface. Nucleic Acids Res. 2008;36:W5–9.View ArticlePubMedPubMed CentralGoogle Scholar
- Megy K, Emrich SJ, Lawson D, Campbell D, Dialynas E, Hughes DS, Koscielny G, Louis C, Maccallum RM, Redmond SN, et al. VectorBase: improvements to a bioinformatics resource for invertebrate vector genomics. Nucleic Acids Res. 2012;40:D729–734.View ArticlePubMedGoogle Scholar
- Yang Z, Kumar S, Nei M. A new method of inference of ancestral nucleotide and amino acid sequences. Genetics. 1995;141:1641–50.PubMedPubMed CentralGoogle Scholar
- Tamura K, Peterson D, Peterson N, Stecher G, Nei M, Kumar S. MEGA5: molecular evolutionary genetics analysis using maximum likelihood, evolutionary distance, and maximum parsimony methods. Mol Biol Evol. 2011;28:2731–9.View ArticlePubMedPubMed CentralGoogle Scholar
- Hall BG. Building phylogenetic trees from molecular data with MEGA. Mol Biol Evol. 2013;30:1229–35.View ArticlePubMedGoogle Scholar
- Sievers F, Wilm A, Dineen D, Gibson TJ, Karplus K, Li W, Lopez R, McWilliam H, Remmert M, Soding J, et al. Fast, scalable generation of high-quality protein multiple sequence alignments using Clustal Omega. Mol Syst Biol. 2011;7:539.View ArticlePubMedPubMed CentralGoogle Scholar