Sympatric Drosophila simulans flies with distinct mtDNA show difference in mitochondrial respiration and electron transport

The role of mitochondrial DNA (mtDNA) in mitochondrial metabolism is understudied yet humans harboring specific mtDNA types age at dissimilar rates, are unequally susceptible to various diseases, and differentially adapt to various environmental conditions. This study compares mitochondrial respiration, proton leak and electron transport of Drosophila simulans males with distinct mtDNA haplogroups (siII and -III) that were collected in sympatry in Kenya. Despite the large divergence among haplogroups there is very low intrahaplogroup variation and no correlated variation in the nuclear genome has been detected. We show that repeatable bioenergetic differences exist between 11d old males harboring siII and siIII mtDNA. Males with siIII mtDNA showed higher (i) state 3 respiration rates from isolated mitochondria for both complex I and complex III based substrates, and (ii) complex IV (cytochrome c oxidase) activity. Males harboring siIII mtDNA had lower (i) hydrogen peroxide formation by both complexes I and III, (ii) proton leak from isolated mitochondria, (iii) mitochondrial ATPase activity, and (iv) mitochondrial cytochrome content. In combination, the results suggest that mitochondria isolated from siIII mtDNA harboring males have more efficient metabolism than siII mtDNA harboring males.     

Sympatric Drosophila simulans flies with distinct mtDNA show age related differences in mitochondrial metabolism

The primary causes of age-related changes in mitochondrial metabolism are not known. The goal of this study is to document the influence of naturally occurring mtDNA variation on age dependent changes in mitochondrial respiration, hydrogen peroxide (H(2)O(2)) generation and antioxidant defenses in the fly Drosophila simulans. Possible changes include an increase in rates of reactive oxygen species production with age and/or an age dependent decrease in antioxidant response. For this study we have used flies harboring distinct siII and siIII mtDNA types. Previously we have shown that males harboring siII mtDNA had higher rates of mitochondrial H(2)O(2) production from complex III at 11d compared to males with the siIII mtDNA type. Here, we corroborate those results and show that Drosophila harboring the siII and siIII mtDNA types exhibit significantly different patterns of pro-oxidant and antioxidant activities as they age. Flies harboring siII mtDNA had higher rates of mitochondrial H(2)O(2) production and manganese superoxide dismutase activity at 11 and 18d of age than siIII mtDNA harboring flies. Copper-zinc superoxide dismutase activity increased from 11 to 25d in siII flies while the accumulation of oxidized glutathione did not change between 11 and 25d. In contrast, siIII harboring flies showed an age dependent increase in H(2)O(2) production, reaching higher production rates on day 25 than that observed in siII flies. Copper-zinc superoxide dismutase activities did not change between 11 and 25d while the oxidized glutathione accumulation increased with age. The results show antioxidant levels correlate with pro-oxidant levels in siII but not siIII flies. These results demonstrate our ability to correlate mtDNA variation with differences in whole mitochondrial physiology and individual complex biochemistry.     

Thermal sensitivity of mitochondrial functions in permeabilized muscle fibers from two populations of Drosophila simulans with divergent mitotypes

In ectotherms, the external temperature is experienced by the mitochondria, and the mitochondrial respiration of different genotypes is likely to change as a result. Using high-resolution respirometry with permeabilized fibers (an in situ approach), we tried to identify differences in mitochondrial performance and thermal sensitivity of two Drosophila simulans populations with two different mitochondrial types (siII and siIII) and geographical distributions. Maximal state 3 respiration rates obtained with electrons converging at the Q junction of the electron transport system (ETS) differed between the mitotypes at 24°C. Catalytic capacities were higher in flies harboring siII than in those harboring siIII mitochondrial DNA (2,129 vs. 1,390 pmol O(2)·s(-1)·mg protein(-1)). The cytochrome c oxidase activity was also higher in siII than siIII flies (3,712 vs. 2,688 pmol O(2)·s(-1)·mg protein(-1)). The higher catalytic capacity detected in the siII mitotype could provide an advantage in terms of intensity of aerobic activity, endurance, or both, if the intensity of exercise that can be aerobically performed is partly dictated by the aerobic capacity of the tissue. Moreover, thermal sensitivity results showed that even if temperature affects the catalytic capacity of the different enzymes of the ETS, both mitotypes revealed high tolerance to temperature variation. Previous in vitro study failed to detect any consistent functional mitochondrial differences between the same mitotypes. We conclude that the in situ approach is more sensitive and that the ETS is a robust system in terms of functional and regulatory properties across a wide range of temperatures.     

Influence of two Wolbachia strains on population structure of East African Drosophila simulans

Drosophila simulans is hypothesized to have originated in continental East Africa or Madagascar. In this study, we investigated evolutionary forces operating on mitochondrial DNA (mtDNA) in populations of D. simulans from Zimbabwe, Malawi, Tanzania, and Kenya. Variation in mtDNA may be affected by positive selection, background selection, demographic history, and/or any maternally inherited factor such as the bacterial symbiont Wolbachia. In East Africa, the wRi and wMa Wolbachia strains associate with the siII or siIII mitochondrial haplogroups, respectively. To ask how polymorphism relates to Wolbachia infection status, we sequenced 1776 bp of mitochondrial DNA and 1029 bp of the X-linked per locus from 79 lines. The two southern populations were infected with wRi and exhibited significantly reduced mtDNA variation, while Wolbachia-uninfected siII flies from Tanzania and Kenya showed high levels of mtDNA polymorphism. These are the first known populations of D. simulans that do not exhibit reduced mtDNA variation. We observed no mitochondrial variation in the siIII haplogroup regardless of Wolbachia infection status, suggesting positive or background selection. These populations offer a unique opportunity to monitor evolutionary dynamics in ancestral populations that harbor multiple strains of Wolbachia.     

Sequential evolution of a symbiont inferred from the host: Wolbachia and Drosophila simulans

This study aims to unravel the biogeography of a model symbiont/host system by exploiting the prediction that a symbiont will leave a signature of infection on the host. Specifically, a global sample of 1,442 Drosophila simulans from 33 countries and 64 sampling localities was employed to infer the phylogeography of the maternally inherited alpha-proteobacteria Wolbachia. Phylogenetic analyses, from three symbiont genes and 24 mtDNA genomes (excluding the A + T-rich region), showed that each of four Wolbachia strains infected D. simulans once. The global distribution and abundance of the Wolbachia strains and the three mtDNA haplogroups (D. simulans siI, siII and siIII) was then determined. Finally, network analyses of variable regions within siI (584 bp from seven additional lines) and siII (1,701 bp from 383 lines) facilitated a detailed biogeographic discussion. There is little variation in siIII and the haplogroup is restricted in its distribution. These data show how the history of an infection can be mapped by combining data from the symbiont and the host. They say little about the organismal history of the host because the mtDNA genome is a biased representation of the whole genome.     

Mitochondrial DNA heteroplasmy maintained in natural populations of Drosophila simulans in Réunion

Mitochondrial DNA (mtDNA) variation in Drosophila simulans was studied to determine whether the cytoplasmic state of mtDNA heteroplasmy persists in natural populations in Réunion. For this purpose, 172 isofemale lines, newly collected from two local populations, were examined, among which three types of mtDNA (siII, siIII and siIII') were found, based on the Hpa II restriction pattern. Ten of the lines were heteroplasmic for a combination of siII and siIII, as determined by autoradiography. The same type of heteroplasmy had been noted in one of the two local populations 8 years before (Satta et al. 1988). The present results suggest that the heteroplasmic state occurs recurrently in natural populations of D. simulans in Réunion.     

Mitochondrial DNA in the Drosophila Melanogaster Complex

Mitochondrial DNA in the complex Drosophila melanogaster was among the first studied in metazoans. The variability of the molecule was extensively studied using restriction enzymes, gene sequence and recently sequence of the whole coding region. Within the complex, seven major haplotypes have been described, one (me) in D. melanogaster, three in D. simulans (siI, siII, siIII), two in D. mauritiana (maI, maII), and one in D. sechellia (se). The molecular distance between the haplotypes is comprised between 1 and 5%, except for siII and maI, which are virtually identical. The nucleotide diversity within each of these haplotypes is very low, varying from 0 to 0.0005. Most of the cytoplasms are infected by the bacterium Wolbachia and different bacterial strains infect cytoplasms harboring different mtDNA types. mtDNA polymorphism is discussed in relation with Wolbachia, nuclear polymorphism and speciation events.     

A Wolbachia-associated fitness benefit depends on genetic background in Drosophila simulans

The alpha-proteobacteria Wolbachia infect a number of insect species and influence host reproduction to favour the spread of infected females through a population. The fitness effect of this infection is important in understanding the spread and maintenance of Wolbachia within and among host populations. However, a full elucidation of fitness effect requires careful control of host genetic background. Here, I transferred a single clone of Wolbachia (the wHa strain) into three genetically distinct isofemale lines of the fly Drosophila simulans using microinjection methodology. These lines carried one of the three described mitochondrial haplogroups (siI, siII or siIII) and differ in nuclear genome as well. Population cage assays showed that wHa-infected siIII flies enjoyed a dramatic fitness benefit compared to uninfected siIII. In contrast, wHa did not affect the fitness of siI or siII flies. This study points to the importance of host-by-symbiont interaction terms that may play an important role in organismal-fitness.     

Viability of cytochrome c genotypes depends on cytoplasmic backgrounds in Tigriopus californicus

Because of their extensive functional interaction, mitochondrial DNA (mtDNA) and nuclear genes may evolve to form coadapted complexes within reproductively isolated populations. As a consequence of coadaptation, the fitness of particular nuclear alleles may depend on mtDNA genotype. Among populations of the copepod Tigriopus californicus, there are high levels of amino acid substitutions in both the mtDNA genes encoding subunits of cytochrome c oxidase (COX) and the nuclear gene encoding cytochrome c (CYC), the substrate for COX. Because of the functional interaction between enzyme and substrate proteins, we hypothesized that the fitness of CYC genotypes would depend on mtDNA genotype. To test this hypothesis, segregation ratios for CYC and a second nuclear marker (histone H1) unrelated to mitochondrial function were scored in F2 progeny of several reciprocal interpopulation crosses. Genotypic ratios at the CYC locus (but not the H1 locus) differed between reciprocal crosses and differed from expected Mendelian ratios, suggesting that CYC genotypic fitnesses were strongly influenced by cytoplasmic (including mtDNA) background. However, in most cases the nature of the deviations from Mendelian ratios and differences between reciprocal crosses are not consistent with simple coevolution between CYC and mtDNA background. In a cross in which both newly hatched larvae and adults were sampled, only the adult sample showed deviations from Mendelian ratios, indicating that genotypic viabilities differed. In two of six crosses, large genotypic ratio differences for CYC were observed between the sexes. These results suggest that significant variation in nuclear-mtDNA coadaptation may exist between T. californicus populations and that the relative viability of specific cytonuclear allelic combinations is somehow affected by sex.     

A Candidate Complex Approach to Study Functional Mitochondrial DNA Changes: Sequence Variation and Quaternary Structure Modeling of Drosophila simulans Cytochrome c Oxidase

A problem with studying evolutionary dynamics of mitochondrial (mt) DNA is that classical population genetic techniques cannot identify selected substitutions because of genetic hitchhiking. We circumvented this problem by employing a candidate complex approach to study sequence variation in cytochrome c oxidase (COX) genes within and among three distinct Drosophila simulans mtDNA haplogroups. First, we determined sequence variation in complete coding regions for all COX mtDNA and nuclear loci and their isoforms. Second, we constructed a quaternary structure model of D. simulans COX. Third, we predicted that six of nine amino acid changes in D. simulans mtDNA are likely to be functionally important. Of these seven, genetic crosses can experimentally determine the functional significance of three. Fourth, we identified two single amino acid changes and a deletion of two consecutive amino acids in nuclear encoded COX loci that are likely to influence cytochrome c oxidase activity. These data show that linking population genetics and quaternary structure modeling can lead to functional predictions of specific mtDNA amino acid mutations and validate the candidate complex approach. Electronic supplementary material The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

Excess amino acid polymorphism in mitochondrial DNA: contrasts among genes from Drosophila, mice, and humans

Recent studies of mitochondrial DNA (mtDNA) variation in mammals and Drosophila have shown an excess of amino acid variation within species (replacement polymorphism) relative to the number of silent and replacement differences fixed between species. To examine further this pattern of nonneutral mtDNA evolution, we present sequence data for the ND3 and ND5 genes from 59 lines of Drosophila melanogaster and 29 lines of D. simulans. Of interest are the frequency spectra of silent and replacement polymorphisms, and potential variation among genes and taxa in the departures from neutral expectations. The Drosophila ND3 and ND5 data show no significant excess of replacement polymorphism using the McDonald-Kreitman test. These data are in contrast to significant departures from neutrality for the ND3 gene in mammals and other genes in Drosophila mtDNA (cytochrome b and ATPase 6). Pooled across genes, however, both Drosophila and human mtDNA show very significant excesses of amino acid polymorphism. Silent polymorphisms at ND5 show a significantly higher variance in frequency than replacement polymorphisms, and the latter show a significant skew toward low frequencies (Tajima's D = -1.954). These patterns are interpreted in light of the nearly neutral theory where mildly deleterious amino acid haplotypes are observed as ephemeral variants within species but do not contribute to divergence. The patterns of polymorphism and divergence at charge-altering amino acid sites are presented for the Drosophila ND5 gene to examine the evolution of functionally distinct mutations. Excess charge-altering polymorphism is observed at the carboxyl terminal and excess charge-altering divergence is detected at the amino terminal. While the mildly deleterious model fits as a net effect in the evolution of nonrecombining mitochondrial genomes, these data suggest that opposing evolutionary pressures may act on different regions of mitochondrial genes and genomes.      

Nuclear-mitochondrial epistasis and drosophila aging: introgression of Drosophila simulans mtDNA modifies longevity in D. melanogaster nuclear backgrounds

Under the mitochondrial theory of aging, physiological decline with age results from the accumulated cellular damage produced by reactive oxygen species generated during electron transport in the mitochondrion. A large body of literature has documented age-specific declines in mitochondrial function that are consistent with this theory, but relatively few studies have been able to distinguish cause from consequence in the association between mitochondrial function and aging. Since mitochondrial function is jointly encoded by mitochondrial (mtDNA) and nuclear genes, the mitochondrial genetics of aging should be controlled by variation in (1) mtDNA, (2) nuclear genes, or (3) nuclear-mtDNA interactions. The goal of this study was to assess the relative contributions of these factors in causing variation in Drosophila longevity. We compared strains of flies carrying mtDNAs with varying levels of divergence: two strains from Zimbabwe (<20 bp substitutions between mtDNAs), strains from Crete and the United States (approximately 20-40 bp substitutions between mtDNAs), and introgression strains of Drosophila melanogaster carrying mtDNA from Drosophila simulans in a D. melanogaster Oregon-R chromosomal background (>500 silent and 80 amino acid substitutions between these mtDNAs). Longevity was studied in reciprocal cross genotypes between pairs of these strains to test for cytoplasmic (mtDNA) factors affecting aging. The intrapopulation crosses between Zimbabwe strains show no difference in longevity between mtDNAs; the interpopulation crosses between Crete and the United States show subtle but significant differences in longevity; and the interspecific introgression lines showed very significant differences between mtDNAs. However, the genotypes carrying the D. simulans mtDNA were not consistently short-lived, as might be predicted from the disruption of nuclear-mitochondrial coadaptation. Rather, the interspecific mtDNA strains showed a wide range of variation that flanked the longevities seen between intraspecific mtDNAs, resulting in very significant nuclear x mtDNA epistatic interaction effects. These results suggest that even "defective" mtDNA haplotypes could extend longevity in different nuclear allelic backgrounds, which could account for the variable effects attributable to mtDNA haplogroups in human aging.     

Mitochondrial genotype affects fitness in Drosophila simulans

Drosophila simulans is known to harbor three distinct mitochondrial DNA (mtDNA) haplotype groups (siI, -II, and -III) with nearly 3.0% interhaplotypic divergence but <0.06% intrahaplotypic diversity. With the large amount of genetic variation in this system, the potential power to detect intraspecific fitness differences in fly lines that carry distinct haplotypes is great. We test three life-history traits on fly lines with known sequence differences in the mtDNA genome after controlling the nuclear genome by backcrossing. We find that flies with the siI haplotype are fastest developing and have the lowest probability of surviving to three experimental periods (2-6, 12-17, and 34-39 days of age). Wild-type males with siIII mtDNA were more active while disruption of specific coadapted nucleo-mitochondrial complexes caused a significant decrease in activity. These results are discussed in the context of the geographic distribution of each haplotype.     

Cytonuclear coadaptation in Drosophila: disruption of cytochrome c oxidase activity in backcross genotypes

Abstract The cytochrome c oxidase enzyme (COX) is comprised of 10 nuclear-encoded subunits and three mito-chondrial-encoded subunits in close physical association in the inner mitochondrial membrane. COX passes electrons from cytochrome c to molecular oxygen and pumps protons into the inner mitochondrial space for ATP production. Selection on nuclear-mitochondrial interactions within species should lead to coadaptation of the proteins comprising this important enzyme. Under this model, there should be relatively little disruption of COX activity when mitochondrial genomes are crossed among strains within species. A more pronounced disruption of activity is expected when the mitochondrial genome is expressed in the nuclear background of a different species. We test these hypotheses in Drosophila using hybridization and backcrossing among lines of D. simulans and D. mauritiana. Disrupted cytonuclear genotypes were constructed using backcrosses between two lines of D. simulans (siI and si II ) that introduced each divergent mitochondrial DNA (mtDNA) into each nuclear background due to maternal inheritance of mtDNA. Similar crosses were used to introduce eachD. simulans mtDNA into the D. mauritiana maI nuclear background. Reconstituted cytonuclear control genotypes were constructed by backcrossing the initial F1 females to males of the maternal genotype. COX enzyme activities were compared among these disrupted and reconstituted backcross genotypes within and between species. The disruption effect on COX activity was restricted to males of interspecific genotypes. These data support the coadaptation hypothesis and are consistent with predictions that the evolution of modifiers of male mitochondrial dysfunction is hindered by the maternal inheritance of mtDNA. New sequence data for nuclear encoded subunits of COX identified amino acids that may play a role in the disruption effect.     

Unraveling Selection in the Mitochondrial Genome of Drosophila

We examine mitochondrial DNA variation at the cytochrome b locus within and between three species of Drosophila to determine whether patterns of variation conform to the predictions of neutral molecular evolution. The entire 1137-bp cytochrome b locus was sequenced in 16 lines of Drosophila melanogaster, 18 lines of Drosophila simulans and 13 lines of Drosophila yakuba. Patterns of variation depart from neutrality by several test criteria. Analysis of the evolutionary clock hypothesis shows unequal rates of change along D. simulans lineages. A comparison within and between species of the ratio of amino acid replacement change to synonymous change reveals a relative excess of amino acid replacement polymorphism compared to the neutral prediction, suggestive of slightly deleterious or diversifying selection. There is evidence for excess homozygosity in our world wide sample of D. melanogaster and D. simulans alleles, as well as a reduction in the number of segregating sites in D. simulans, indicative of selective sweeps. Furthermore, a test of neutrality for codon usage shows the direction of mutations at third positions differs among different topological regions of the gene tree. The analyses indicate that molecular variation and evolution of mtDNA are governed by many of the same selective forces that have been shown to govern nuclear genome evolution and suggest caution be taken in the use of mtDNA as a ``neutral' molecular marker.     

Deletion of mtDNA disrupts mitochondrial function and structure, but not biogenesis

The role of nuclear DNA (nDNA)-encoded proteins in the regulation of mitochondrial fission and fusion has been documented, yet the role of mitochondrial DNA (mtDNA) and encoded proteins in mitochondrial biogenesis remains unknown. Long-term treatment of a lymphoblastoid cell line Molt-4 with ethidium bromide generated mtDNA-deficient rho0 mutants. Depletion of mtDNA in rho0 cells produced functional and morphological changes in mitochondria without affecting the nuclear genome and encoded proteins. Indeed, the gene encoding subunit II of mitochondrial cytochrome c oxidase (COX II), a prototypical mitochondrial gene, was reduced in rho0 mutants blunting the activity of mitochondrial cytochrome coxidase. Yet, the amount of the nuclear beta-actin gene and the activity of citrate synthase, a mitochondrial matrix enzyme encoded by nDNA, remained unaffected in rho0 cells. Loss of mtDNA in rho0 cells was associated with significant distortion of mitochondrial structure, decreased electron density of the matrix and disorganized inner and outer membranes, resulting in the appearance of 'ghost-like' mitochondria. However, the number of mitochondria-like structures was not significantly different between mtDNA-deficient and parental cells. Thus, we conclude that cells lacking mtDNA still generate mitochondrial scaffolds, albeit with aberrant function.     

Within‐population genetic effects of mtDNA on metabolic rate in Drosophila subobscura

A growing body of research supports the view that within‐species sequence variation in the mitochondrial genome (mtDNA) is functional, in the sense that it has important phenotypic effects. However, most of this empirical foundation is based on comparisons across populations, and few studies have addressed the functional significance of mtDNA polymorphism within populations. Here, using mitonuclear introgression lines, we assess differences in whole‐organism metabolic rate of adult Drosophila subobscura fruit flies carrying either of three different sympatric mtDNA haplotypes. We document sizeable, up to 20%, differences in metabolic rate across these mtDNA haplotypes. Further, these mtDNA effects are to some extent sex specific. We found no significant nuclear or mitonuclear genetic effects on metabolic rate, consistent with a low degree of linkage disequilibrium between mitochondrial and nuclear genes within populations. The fact that mtDNA haplotype variation within a natural population affects metabolic rate, which is a key physiological trait with important effects on life‐history traits, adds weight to the emergent view that mtDNA haplotype variation is under natural selection and it revitalizes the question as to what processes act to maintain functional mtDNA polymorphism within populations.     

Divergent evolution of life span associated with mitochondrial DNA evolution

Mitochondria play a key role in ageing. The pursuit of genes that regulate variation in life span and ageing have shown that several nuclear‐encoded mitochondrial genes are important. However, the role of mitochondrial encoded genes (mtDNA) is more controversial and our appreciation of the role of mtDNA for the evolution of life span is limited. We use replicated lines of seed beetles that have been artificially selected for long or short life for >190 generations, now showing dramatic phenotypic differences, to test for a possible role of mtDNA in the divergent evolution of ageing and life span. We show that these divergent selection regimes led to the evolution of significantly different mtDNA haplotype frequencies. Selection for a long life and late reproduction generated positive selection for one specific haplotype, which was fixed in most such lines. In contrast, selection for reproduction early in life led to both positive selection as well as negative frequency‐dependent selection on two different haplotypes, which were both present in all such lines. Our findings suggest that the evolution of life span was in part mediated by mtDNA, providing support for the emerging general tenet that adaptive evolution of life‐history syndromes may involve mtDNA.     

Molecular and ecological signs of mitochondrial adaptation: consequences for introgression?

The evolution of the mitochondrial genome and its potential adaptive impact still generates vital debates. Even if mitochondria have a crucial functional role, as they are the main cellular energy suppliers, mitochondrial DNA (mtDNA) introgression is common in nature, introducing variation in populations upon which selection may act. Here we evaluated whether the evolution of mtDNA in a rodent species affected by mtDNA introgression is explained by neutral expectations alone. Variation in one mitochondrial and six nuclear markers in Myodes glareolus voles was examined, including populations that show mtDNA introgression from its close relative, Myodes rutilus. In addition, we modelled protein structures of the mtDNA marker (cytochrome b) and estimated the environmental envelopes of mitotypes. We found that massive mtDNA introgression occurred without any trace of introgression in the analysed nuclear genes. The results show that the native glareolus mtDNA evolved under past positive selection, suggesting that mtDNA in this system has selective relevance. The environmental models indicate that the rutilus mitotype inhabits colder and drier habitats than the glareolus one that can result from local adaptation or from the geographic context of introgression. Finally, homology models of the cytochrome b protein revealed a substitution in rutilus mtDNA in the vicinity of the catalytic fraction, suggesting that differences between mitotypes may result in functional changes. These results suggest that the evolution of mtDNA in Myodes may have functional, ecological and adaptive significance. This work opens perspective onto future experimental tests of the role of natural selection in mtDNA introgression in this system.     

Mitochondrial proteome: altered cytochrome c oxidase subunit levels in prostate cancer

Laser capture microdissection was combined with reverse phase protein lysate arrays to quantitatively analyze the ratios of mitochondrial encoded cytochrome c oxidase subunits to nuclear encoded cytochrome c oxidase subunits, and to correlate the ratios with malignant progression in human prostate tissue specimens. Cytochrome c oxidase subunits I-III comprise the catalytic core of the enzyme and are all synthesized from mitochondrial DNA. The remaining subunits (IV-VIII) are synthesized from cellular nuclear DNA. A significant (P < 0.001, 30/30 prostate cases) shift in the relative concentrations of nuclear encoded cytochrome c oxidase subunits IV, Vb, and VIc compared to mitochondrial encoded cytochrome c oxidase subunits I and II was noted during the progression of prostate cancer from normal epithelium through premalignant lesions to invasive carcinoma. Significantly, this shift was discovered to begin even in the premalignant stage. Reverse phase protein lysate array-based observations were corroborated with immunohistochemistry, and extended to a few human carcinomas in addition to prostate. This analysis points to a role for nuclear DNA encoded mitochondrial proteins in carcinogenesis; underscoring their potential as targets for therapy while highlighting the need for full characterization of the mitochondrial proteome.