Transcription of rat mitochondrial NADH-dehydrogenase subunits. Presence of antisense and precursor RNA species

Transcriptome analysis of Edwardsiella piscicida during intracellular infection reveals excludons are involved with the activation of a mitochondrion-like energy generation program

Phylogenetic evidence suggests a shared ancestry between mitochondria and modern Proteobacteria, a phylum including several genera of intracellular pathogens. Studying these diverse pathogens, particularly during intracellular infection of their hosts, can reveal characteristics potentially representative of the mitochondrial-Proteobacterial ancestor by identifying traits shared with mitochondria. While transcriptomic approaches can provide global insights into intracellular acclimatization by pathogens, they are often limited by excess host RNAs in extracts. Here, we developed a method employing magnetic nanoparticles to enrich RNA from an intracellular Gammaproteobacterium, Edwardsiella piscicida, within zebrafish, Danio rerio, fin fibroblasts, enabling comprehensive exploration of the bacterial transcriptome. Our findings revealed that the intracellular E. piscicida transcriptome reflects a mitochondrion-like energy generation program characterized by the suppression of glycolysis and sugar transport, coupled with upregulation of the tricarboxylic acid (TCA) cycle and alternative import of simple organic acids that directly flux into TCA cycle intermediates or electron transport chain donors. Additionally, genes predicted to be members of excludons, loci of gene pairs antagonistically co-regulated by overlapping antisense transcription, are significantly enriched in the set of all genes with perturbed sense and antisense transcription, suggesting a general but important involvement of excludons with intracellular acclimatization. Notably, genes involved with the activation of the mitochondrion-like energy generation program, specifically with metabolite import and glycolysis, are also members of predicted excludons. Other intracellular Proteobacterial pathogens appear to employ a similar mitochondrion-like energy generation program, suggesting a potentially conserved mechanism for optimized energy acquisition from hosts centered around the TCA cycle.IMPORTANCEPhylogenetic evidence suggests that mitochondria and Proteobacteria, a phylum encompassing various intracellular pathogens, share a common ancestral lineage. In this study, we developed a novel method employing magnetic nanoparticles to explore the transcriptome of an aquatic Gammaproteobacterium, Edwardsiella piscicida, during intracellular infection of host cells. We show that the strategy E. piscicida uses to generate energy strikingly mirrors the function of mitochondria-energy generators devoid of glycolytic processes. Notably, several implicated genes are members of excludons-gene pairs antagonistically co-regulated by overlapping antisense transcription. Other intracellular Proteobacterial pathogens appear to adopt a similar mitochondrion-like energy generation program, indicating a possibly conserved strategy for optimized energy acquisition from hosts centered around the tricarboxylic acid cycle.

Hypertension-associated mitochondrial DNA 4401A>G mutation caused the aberrant processing of tRNAMet, all 8 tRNAs and ND6 mRNA in the light-strand transcript

Mitochondrial tRNA processing defects were associated with human diseases but their pathophysiology remains elusively. The hypertension-associated m.4401A>G mutation resided at a spacer between mitochondrial tRNAMet and tRNAGln genes. An in vitro processing experiment revealed that the m.4401A>G mutation caused 59% and 69% decreases in the 5' end processing efficiency of tRNAGln and tRNAMet precursors, catalyzed by RNase P, respectively. Using human umbilical vein endothelial cells-derived cybrids, we demonstrated that the m.4401A>G mutation caused the decreases of all 8 tRNAs and ND6 and increases of longer and uncleaved precursors from the Light-strand transcript. Conversely, the m.4401A>G mutation yielded the reduced levels of tRNAMet level but did not change the levels of other 13 tRNAs, 12 mRNAs including ND1, 12S rRNA and 16S rRNA from the Heavy-strand transcript. These implicated the asymmetrical processing mechanisms of H-strand and L-strand polycistronic transcripts. The tRNA processing defects play the determined roles in the impairing mitochondrial translation, respiratory deficiency, diminishing membrane potential, increasing production of reactive oxygen species and altering autophagy. Furthermore, the m.4401A>G mutation altered the angiogenesis, evidenced by aberrant wound regeneration and weaken tube formation in mutant cybrids. Our findings provide new insights into the pathophysiology of hypertension arising from mitochondrial tRNA processing defects.

Variation and evolution of polyadenylation profiles in sauropsid mitochondrial mRNAs as deduced from the high-throughput RNA sequencing

BACKGROUND

Genes encoded in vertebrate mitochondrial DNAs are transcribed as a polycistronic transcript for both strands, which is later processed into individual mRNAs, rRNAs and tRNAs, followed by modifications, such as polyadenylation at the 3' end of mRNAs. Although mechanisms of the mitochondrial transcription and RNA processing have been extensively studied using some model organisms, structural variability of mitochondrial mRNAs across different groups of vertebrates is poorly understood. We conducted the high-throughput RNA sequencing to identify major polyadenylation sites for mitochondrial mRNAs in the Japanese grass lizard, Takydromus tachydromoides and compared the polyadenylation profiles with those identified similarly for 23 tetrapod species, featuring sauropsid taxa (reptiles and birds).

RESULTS

As compared to the human, a major polyadenylation site for the NADH dehydrogenase subunit 5 mRNA of the grass lizard was located much closer to its stop codon, resulting in considerable truncation of the 3' untranslated region for the mRNA. Among the other sauropsid taxa, several distinct polyadenylation profiles from the human counterpart were found for different mRNAs. They included various truncations of the 3' untranslated region for NADH dehydrogenase subunit 5 mRNA in four taxa, bird-specific polyadenylation of the light-strand-transcribed NADH dehydrogenase subunit 6 mRNA, and the combination of the ATP synthase subunit 8/6 mRNA with a neighboring mRNA into a tricistronic mRNA in the side-necked turtle Pelusios castaneus. In the last case of P. castaneus, as well as another example for NADH dehydrogenase subunit 1 mRNAs of some birds, the association between the polyadenylation site change and the gene overlap was highlighted. The variations in the polyadenylation profile were suggested to have arisen repeatedly in diverse sauropsid lineages. Some of them likely occurred in response to gene rearrangements in the mitochondrial DNA but the others not.

CONCLUSIONS

These results demonstrate structural variability of mitochondrial mRNAs in sauropsids. The efficient and comprehensive characterization of the mitochondrial mRNAs will contribute to broaden our understanding of their structural and functional evolution.

A mitochondria-specific isoform of FASTK is present in mitochondrial RNA granules and regulates gene expression and function

The mitochondrial genome relies heavily on post-transcriptional events for its proper expression, and misregulation of this process can cause mitochondrial genetic diseases in humans. Here, we report that a novel translational variant of Fas-activated serine/threonine kinase (FASTK) co-localizes with mitochondrial RNA granules and is required for the biogenesis of ND6 mRNA, a mitochondrial-encoded subunit of the NADH dehydrogenase complex (complex I). We show that ablating FASTK expression in cultured cells and mice results specifically in loss of ND6 mRNA and reduced complex I activity in vivo. FASTK binds at multiple sites along the ND6 mRNA and its precursors and cooperates with the mitochondrial degradosome to ensure regulated ND6 mRNA biogenesis. These data provide insights into the mechanism and control of mitochondrial RNA processing within mitochondrial RNA granules.

Tissue-specific mtDNA abundance from exome data and its correlation with mitochondrial transcription, mass and respiratory activity

Eukaryotic cells contain a population of mitochondria, variable in number and shape, which in turn contain multiple copies of a tiny compact genome (mtDNA) whose expression and function is strictly coordinated with the nuclear one. mtDNA copy number varies between different cell or tissues types, both in response to overall metabolic and bioenergetics demands and as a consequence or cause of specific pathological conditions. Here we present a novel and reliable methodology to assess the effective mtDNA copy number per diploid genome by investigating off-target reads obtained by whole-exome sequencing (WES) experiments. We also investigate whether and how mtDNA copy number correlates with mitochondrial mass, respiratory activity and expression levels. Analyzing six different tissues from three age- and sex-matched human individuals, we found a highly significant linear correlation between mtDNA copy number estimated by qPCR and the frequency of mtDNA off target WES reads. Furthermore, mtDNA copy number showed highly significant correlation with mitochondrial gene expression levels as measured by RNA-Seq as well as with mitochondrial mass and respiratory activity. Our methodology makes thus feasible, at a large scale, the investigation of mtDNA copy number in diverse cell-types, tissues and pathological conditions or in response to specific treatments.

The post-transcriptional life of mammalian mitochondrial RNA

Mammalian mitochondria contain their own genome that encodes mRNAs for thirteen essential subunits of the complexes performing oxidative phosphorylation as well as the RNA components (two rRNAs and 22 tRNAs) needed for their translation in mitochondria. All RNA species are produced from single polycistronic precursor RNAs, yet the relative concentrations of various RNAs differ significantly. This underscores the essential role of post-transcriptional mechanisms that control the maturation, stability and translation of mitochondrial RNAs. The present review provides a detailed summary on the role of RNA maturation in the regulation of mitochondrial gene expression, focusing mainly on messenger RNA polyadenylation and stability control. Furthermore, the role of mitochondrial ribosomal RNA stability, processing and modifications in the biogenesis of the mitochondrial ribosome is discussed.

Characterization of mitochondrial mRNAs in codfish reveals unique features compared to mammals

Expression and processing of mitochondrial gene transcripts are fundamental to mitochondrial function, but information from early vertebrates like teleost fishes is essentially lacking. We have analyzed mitogenome sequences of ten codfishes (family Gadidae), and provide complete sequences from three new species (Saithe, Pollack and Blue whiting). Characterization of the mitochondrial mRNAs in Saithe and Atlantic cod identified a set of ten poly(A) transcripts, and six UAA stop codons are generated by posttranscriptional polyadenylation. Structural assessment of poly(A) sites is consistent with an RNaseP cleavage activity 5' of tRNA acceptor-like stems. COI, ND5 and ND6 mRNAs were found to harbor 3' UTRs with antisense potential extending into neighboring gene regions. While the 3' UTR of COI mRNA is complementary to the tRNA(Ser UCN) and highly similar to that detected in human mitochondria, the ND5 and ND6 3' UTRs appear more heterogenic. Deep sequencing confirms expression of all mitochondrial mRNAs and rRNAs, and provides information about the precise 5' ends in mature transcripts. Our study supports an overall evolutionary conservation in mitochondrial RNA processing events among vertebrates, but reveals some unique 5' and 3' end characteristics in codfish mRNAs with implications to antisense regulation of gene expression.

Unequal and genotype-dependent expression of mitochondrial genes in larvae of the pacific oyster Crassostrea gigas

Mitochondria are essential for regulation of energy metabolism, but little is known about patterns of mitochondrial genome expression in invertebrates. To explore the association of mitochondrial expression with differential growth of Crassostrea gigas, the Pacific oyster, we crossed two inbred lines to produce inbred and hybrid larvae, which grew at different rates under the same environmental conditions. Using high-throughput cloning and sequencing methods, we identified 1.1 million expressed sequence tags from the mitochondrial genome, 96.7% of which were perfect matches to genes targeted by the method. Expression varied significantly among genes, ranging over nearly four orders of magnitude, from mt:lRNA, which constituted 21% of all transcripts, to mt:CoII, which constituted less than 0.02% of all transcripts. Variable expression of genes coding for subunits of macromolecular complexes (e.g., mt:CoI and mt:CoII) implies that stoichiometry in these complexes must be regulated post-transcriptionally. Surprisingly, the mitochondrial transcriptome contained non-coding transcripts, which may play a role in the regulation of mitochondrial function. Finally, mitochondrial expression depended strongly on maternal factors and nuclear-cytoplasmic interactions, which may explain previously observed growth differences between reciprocal hybrids. Differences in mitochondrial gene expression could provide a biochemical index for the metabolic basis of genetically determined differences in larval growth.

Mitochondrial DNA transcription and diseases: past, present and future

The transcription of mitochondrial DNA has been studied for 30 years. However, many of the earlier observations are still unsolved. In this review we will recall the basis of mitochondrial DNA transcription, established more than twenty years ago, will include some of the recent progress in the understanding of this process and will suggest hypotheses for some of the unexplained topics. Moreover, we will show some examples of mitochondrial pathology due to altered transcription and RNA metabolism.

Oxygen tension regulates mitochondrial DNA-encoded complex I gene expression

Oxygen is a major regulator of nuclear gene expression. However, although mitochondria consume almost all of the O2 available to the cells, little is known about how O2 tension influences the expression of the mitochondrial genome. We show in O2-sensitive excitable rat PC12 cells that, among the mtDNA-encoded genes, hypoxia produced a specific down-regulation of the transcripts encoding mitochondrial complex I NADH dehydrogenase (ND) subunits, particularly ND4 and ND5 mRNAs and a stable mRNA precursor containing the ND5 and cytochrome b genes. This unprecedented effect of hypoxia was fast (developed in <30 min) and fairly reversible and occurred at moderate levels of hypoxia (O2 tensions in the range of 20-70 mm Hg). Hypoxic down-regulation of the mitochondrial complex I genes was paralleled by the reduction of complex I activity and was retarded by iron chelation, suggesting that an iron-dependent post-transcriptional mechanism could regulate mitochondrial mRNA stability. It is known that cell respiration is under tight control by the amount of proteins in mitochondrial complexes of the electron transport chain. Therefore, regulation of the expression of the mitochondrial (mtDNA)-encoded complex I subunits could be part of an adaptive mechanism to adjust respiration rate to the availability of O2 and to induce fast adaptive changes in hypoxic cells.

Polyadenylation and degradation of human mitochondrial RNA: the prokaryotic past leaves its mark

RNA polyadenylation serves a purpose in bacteria and organelles opposite from the role it plays in nuclear systems. The majority of nucleus-encoded transcripts are characterized by stable poly(A) tails at their mature 3' ends, which are essential for stabilization and translation initiation. In contrast, in bacteria, chloroplasts, and plant mitochondria, polyadenylation is a transient feature which promotes RNA degradation. Surprisingly, in spite of their prokaryotic origin, human mitochondrial transcripts possess stable 3'-end poly(A) tails, akin to nucleus-encoded mRNAs. Here we asked whether human mitochondria retain truncated and transiently polyadenylated transcripts in addition to stable 3'-end poly(A) tails, which would be consistent with the preservation of the largely ubiquitous polyadenylation-dependent RNA degradation mechanisms of bacteria and organelles. To this end, using both molecular and bioinformatic methods, we sought and revealed numerous examples of such molecules, dispersed throughout the mitochondrial genome. The broad distribution but low abundance of these polyadenylated truncated transcripts strongly suggests that polyadenylation-dependent RNA degradation occurs in human mitochondria. The coexistence of this system with stable 3'-end polyadenylation, despite their seemingly opposite effects, is so far unprecedented in bacteria and other organelles.

Energy biogenesis: one key for coordinating two genomes

In metazoan organisms, energy production is the only example of a process that is under dual genetic control: nuclear and mitochondrial. We used a genomic approach to examine how energy genes of both the nuclear and mitochondrial genomes are coordinated, and discovered a novel genetic regulatory circuit in Drosophila melanogaster that is surprisingly simple and parsimonious. This circuit is based on a single DNA regulatory element and can explain both intra- and inter-genomic coordinated expression of genes involved in energy production, including the full complement of mitochondrial and nuclear oxidative phosphorylation genes, and the genes involved in the Krebs cycle.

Transcript mapping and genome annotation of ascidian mtDNA using EST data

Mitochondrial transcripts of two ascidian species were reconstructed through sequence assembly of publicly available ESTs resembling mitochondrial DNA sequences (mt-ESTs). This strategy allowed us to analyze processing and mapping of the mitochondrial transcripts and to investigate the gene organization of a previously uncharacterized mitochondrial genome (mtDNA). This new strategy would greatly facilitate the sequencing and annotation of mtDNAs. In Ciona intestinalis, the assembled mt-ESTs covered 22 mitochondrial genes ( approximately 12,000 bp) and provided the partial sequence of the mtDNA and the prediction of its gene organization. Such sequences were confirmed by amplification and sequencing of the entire Ciona mtDNA. For Halocynthia roretzi, for which the mtDNA sequence was already available, the inferred mt transcripts allowed better definition of gene boundaries (16S rRNA, ND1, ATP6, and tRNA-Ser genes) and the identification of a new gene (an additional Phe-tRNA). In both species, polycistronic and immature transcripts, creation of stop codons by polyadenylation, tRNA signal processing, and rRNA transcript termination signals were identified, thus suggesting that the main features of mitochondrial transcripts are conserved in Chordata.

Non-coding ribonucleic acids--a class of their own?

Non-coding ribonucleic acids (RNAs) do not contain a peptide-encoding open reading frame and are therefore not translated into proteins. They are expressed in all phyla, and in eukaryotic cells they are found in the nucleus, cytoplasm, and mitochondria. Non-coding RNAs either can exert structural functions, as do transfer and ribosomal RNAs, or they can regulate gene expression. Non-coding RNAs with regulatory functions differ in size ranging from a few nucleotides to over 100 kb and have diverse cell- or development-specific functions. Some of the non-coding RNAs associate with human diseases. This chapter summarizes the current knowledge about regulatory non-coding RNAs.

Large functional range of steady-state levels of nuclear and mitochondrial transcripts coding for the subunits of the human mitochondrial OXPHOS system

We have measured, by reverse transcription and real-time quantitative PCR, the steady-state levels of the mitochondrial and nuclear transcripts encoding several subunits of the human oxidative phosphorylation (OXPHOS) system, in different normal tissues (muscle, liver, trachea, and kidney) and in cultured cells (normal fibroblasts, 143B osteosarcoma cells, 143B206 rho(0) cells). Five mitochondrial transcripts and nine nuclear transcripts were assessed. The measured amounts of these OXPHOS transcripts in muscle samples corroborated data obtained by others using the serial analysis of gene expression (SAGE) method to appraise gene expression in the same type of tissue. Steady-state levels for all the transcripts were found to range over more than two orders of magnitude. Most of the time, the mitochondrial H-strand transcripts were present at higher levels than the nuclear transcripts. The mitochondrial L-strand transcript ND6 was usually present at a low level. Cultured 143B cells contained significantly reduced amounts of mitochondrial transcripts in comparison with the tissue samples. In 143B206 rho(0) cells, fully depleted of mitochondrial DNA, the levels of nuclear OXPHOS transcripts were not modified in comparison with the parental cells. This observation indicated that nuclear transcription is not coordinated with mitochondrial transcription. We also observed that in the different tissues and cells, there is a transcriptional coregulation of all the investigated nuclear genes. Nuclear OXPHOS gene expression seems to be finely regulated.

Changes in gene expression associated with developmental arrest and longevity in Caenorhabditis elegans

Gene expression in a developmentally arrested, long-lived dauer population of Caenorhabditis elegans was compared with a nondauer (mixed-stage) population by using serial analysis of gene expression (SAGE). Dauer (152,314) and nondauer (148,324) SAGE tags identified 11,130 of the predicted 19,100 C. elegans genes. Genes implicated previously in longevity were expressed abundantly in the dauer library, and new genes potentially important in dauer biology were discovered. Two thousand six hundred eighteen genes were detected only in the nondauer population, whereas 2016 genes were detected only in the dauer, showing that dauer larvae show a surprisingly complex gene expression profile. Evidence for differentially expressed gene transcript isoforms was obtained for 162 genes. H1 histones were differentially expressed, raising the possibility of alternative chromatin packaging. The most abundant tag from dauer larvae (20-fold more abundant than in the nondauer profile) corresponds to a new, unpredicted gene we have named tts-1 (transcribed telomere-like sequence), which may interact with telomeres or telomere-associated proteins. Abundant antisense mitochondrial transcripts (2% of all tags), suggest the existence of an antisense-mediated regulatory mechanism in C. elegans mitochondria. In addition to providing a robust tool for gene expression studies, the SAGE approach already has provided the advantage of new gene/transcript discovery in a metazoan.

Differential expression of sense and antisense transcripts of the mitochondrial DNA region coding for ATPase 6 in fetal and adult porcine brain: identification of novel unusually assembled mitochondrial RNAs

The mammalian mitochondrial genome is a double-stranded circular DNA molecule, which is transcribed from both strands as polycistronic RNAs, which are further processed to yield the mature polyadenylated mRNAs, rRNAs and tRNAs. We compared the gene expression patterns of foetal and adult porcine brains and identified a sequence tag from the ATPase 6 region of the mitochondrial genome which, in adult brain, was more abundant in the sense (H-strand) form, but, in foetal brain, more abundant in the antisense form (L-strand). By means of solution hybridisation/S1 nuclease protection assay, Northern blotting, and PCR based techniques, we demonstrated that the ATPase 6 region of the porcine mitochondrial genome is transcribed as co-existing, stable sense and antisense RNAs. Furthermore, we identified sense and antisense transcripts from this region consisting of inversely assembled fragments joined together at a direct repeat of 7 nucleotides. Our results suggest that transcription and post-transcriptional processing of mitochondrial RNAs are much more complex than presently thought.

The deafness-associated mitochondrial DNA mutation at position 7445, which affects tRNASer(UCN) precursor processing, has long-range effects on NADH dehydrogenase subunit ND6 gene expression

The pathogenetic mechanism of the deafness-associated mitochondrial DNA (mtDNA) T7445C mutation has been investigated in several lymphoblastoid cell lines from members of a New Zealand pedigree exhibiting the mutation in homoplasmic form and from control individuals. We show here that the mutation flanks the 3' end of the tRNASer(UCN) gene sequence and affects the rate but not the sites of processing of the tRNA precursor. This causes an average reduction of approximately 70% in the tRNASer(UCN) level and a decrease of approximately 45% in protein synthesis rate in the cell lines analyzed. The data show a sharp threshold in the capacity of tRNASer(UCN) to support the wild-type protein synthesis rate, which corresponds to approximately 40% of the control level of this tRNA. Strikingly, a 7445 mutation-associated marked reduction has been observed in the level of the mRNA for the NADH dehydrogenase (complex I) ND6 subunit gene, which is located approximately 7 kbp upstream and is cotranscribed with the tRNASer(UCN) gene, with strong evidence pointing to a mechanistic link with the tRNA precursor processing defect. Such reduction significantly affects the rate of synthesis of the ND6 subunit and plays a determinant role in the deafness-associated respiratory phenotype of the mutant cell lines. In particular, it accounts for their specific, very significant decrease in glutamate- or malate-dependent O2 consumption. Furthermore, several homoplasmic mtDNA mutations affecting subunits of NADH dehydrogenase may play a synergistic role in the establishment of the respiratory phenotype of the mutant cells.

Mitochondrial Antisense RNA for Cytochrome C Oxidase (MARCO) Can Induce Morphologic Changes and Cell Death in Human Hematopoietic Cell Lines

To identify essential molecules capable of inducing terminal morphologic maturation and cell death of myeloid progenitor cells, we isolated cDNA clones by functional expression cloning using a library constructed from all-trans retinoic acid (ATRA)-treated human promyelocytic HL-60 cells. Clones which induced morphologic changes in HL-60 cells from blastic cells to mature neutrophilic granulocytes were selected. The isolated positive cDNA clone was demonstrated to encode an antisense RNA for cytochrome c oxidase/serine tRNA derived from a mitochondrial gene (MARCO). When MARCO was expressed in HL-60 cells with the lac switch system, blastic cell morphology became neutrophilic after 48-hour incubation with IPTG, and cell death was observed after 3 days. Also, high molecular weight DNA fragmentation was observed after 36 hours in culture. Similar results were observed using transformants from human K562 cells and CMK cells. RT-PCR analysis revealed that MARCO was transcribed in both ATRA and TNF-α systems, and also in human blood neutrophilic granulocytes. Following transfection with cytochrome c oxidase expression plasmids, TNF-α–induced high molecular weight DNA fragmentation in U937 cells and HL-60 cells was inhibited in these transformants. These results indicate that maturational changes in hematopoietic cells and the process of cell death may be induced by mitochondrial respiratory insufficiency, and also that the mitochondrial gene MARCO may be used as one of the candidates for gene supplementation therapy for the acute leukemias.

Mitochondrial Antisense RNA for Cytochrome C Oxidase (MARCO) Can Induce Morphologic Changes and Cell Death in Human Hematopoietic Cell Lines

To identify essential molecules capable of inducing terminal morphologic maturation and cell death of myeloid progenitor cells, we isolated cDNA clones by functional expression cloning using a library constructed from all-trans retinoic acid (ATRA)-treated human promyelocytic HL-60 cells. Clones which induced morphologic changes in HL-60 cells from blastic cells to mature neutrophilic granulocytes were selected. The isolated positive cDNA clone was demonstrated to encode an antisense RNA for cytochrome c oxidase/serine tRNA derived from a mitochondrial gene (MARCO). When MARCO was expressed in HL-60 cells with the lac switch system, blastic cell morphology became neutrophilic after 48-hour incubation with IPTG, and cell death was observed after 3 days. Also, high molecular weight DNA fragmentation was observed after 36 hours in culture. Similar results were observed using transformants from human K562 cells and CMK cells. RT-PCR analysis revealed that MARCO was transcribed in both ATRA and TNF-α systems, and also in human blood neutrophilic granulocytes. Following transfection with cytochrome c oxidase expression plasmids, TNF-α–induced high molecular weight DNA fragmentation in U937 cells and HL-60 cells was inhibited in these transformants. These results indicate that maturational changes in hematopoietic cells and the process of cell death may be induced by mitochondrial respiratory insufficiency, and also that the mitochondrial gene MARCO may be used as one of the candidates for gene supplementation therapy for the acute leukemias.