Identification of primary transcriptional start sites of mouse mitochondrial DNA: accurate in vitro initiation of both heavy- and light-strand transcripts

Comparative mitogenomics of Clupeoid fish provides insights into the adaptive evolution of mitochondrial oxidative phosphorylation (OXPHOS) genes and codon usage in the heterogeneous habitats

Clupeoid fish can be considered excellent candidates to understand the role of mitochondrial DNA in adaptive evolution, as they have colonized different habitats (marine, brackish, freshwater, tropical and temperate regions) over millions of years. Here, we investigate patterns of tRNA location, codon usage bias, and lineage-specific diversifying selection signals to provide novel insights into how evolutionary improvements of mitochondrial metabolic efficiency have allowed clupeids to adapt to different habitats. Based on whole mitogenome data of 70 Clupeoids with a global distribution we find that purifying selection was the dominant force acting and that the mutational deamination pressure in mtDNA was stronger than the codon/amino acid constraints. The codon usage pattern appears evolved to achieve high translational efficiency (codon/amino acid-related constraints), as indicated by the complementarity of most codons to the GT-saturated tRNA anticodon sites (retained by deamination-induced pressure) and usage of the codons of the tRNA genes situated near to the control region (fixed by deamination pressure) where transcription efficiency was high. The observed shift in codon preference patterns between marine and euryhaline/freshwater Clupeoids indicates possible selection for improved translational efficiency in mitochondrial genes while adapting to low-salinity habitats. This mitogenomic plasticity and enhanced efficiency of the metabolic machinery may have contributed to the evolutionary success and abundance of Clupeoid fish.

Complete sequences of the mitochondrial DNA of the wild Gracilariopsis lemaneiformis and two mutagenic cultivated breeds (Gracilariaceae, Rhodophyta)

The complete mitochondrial DNA (mtDNA) of Gracilariopsis lemaneiformis was sequenced (25883 bp) and mapped to a circular model. The A+T composition was 72.5%. Forty six genes and two potentially functional open reading frames were identified. They include 24 protein-coding genes, 2 rRNA genes, 20 tRNA genes and 2 ORFs (orf60, orf142). There is considerable sequence synteny across the five red algal mtDNAs falling into Florideophyceae including Gr. lemaneiformis in this study and previously sequenced species. A long stem-loop and a hairpin structure were identified in intergenic regions of mt genome of Gr. lemaneiformis, which are believed to be involved with transcription and replication. In addition, the mtDNAs of two mutagenic cultivated breeds (“981” and “07-2”) were also sequenced. Compared with the mtDNA of wild Gr. lemaneiformis, the genome size and gene length and order of three strains were completely identical except nine base mutations including eight in the protein-coding genes and one in the tRNA gene. None of the base mutations caused frameshift or a premature stop codon in the mtDNA genes. Phylogenetic analyses based on mitochondrial protein-coding genes and rRNA genes demonstrated Gracilariopsis andersonii had closer phylogenetic relationship with its parasite Gracilariophila oryzoides than Gracilariopsis lemaneiformis which was from the same genus of Gracilariopsis.

Transfer RNA gene arrangement and codon usage in vertebrate mitochondrial genomes: a new insight into gene order conservation

Background

Mitochondrial (mt) gene arrangement has been highly conserved among vertebrates from jawless fishes to mammals for more than 500 million years. It remains unclear, however, whether such long-term persistence is a consequence of some constraints on the gene order.

Results

Based on the analysis of codon usage and tRNA gene positions, we suggest that tRNA gene order of the typical vertebrate mt-genomes may be important for their translational efficiency. The vertebrate mt-genome encodes 2 rRNA, 22 tRNA, and 13 transmembrane proteins consisting mainly of hydrophobic domains. We found that the tRNA genes specifying the hydrophobic residues were positioned close to the control region (CR), where the transcription efficiency is estimated to be relatively high. Using 47 vertebrate mt-genome sequences representing jawless fishes to mammals, we further found a correlation between codon usage and tRNA gene positions, implying that highly-used tRNA genes are located close to the CR. In addition, an analysis considering the asymmetric nature of mtDNA replication suggested that the tRNA loci that remain in single-strand for a longer time tend to have more guanine and thymine not suffering deamination mutations in their anticodon sites.

Conclusions

Our analyses imply the existence of translational constraint acting on the vertebrate mt-gene arrangement. Such translational constraint, together with the deamination-related constraint, may have contributed to long-term maintenance of gene order.

PCR-based cloning of the complete mouse mitochondrial genome and stable engineering in Escherichia coli

We have devised a method for cloning an entire mammalian mitochondrial genome (mtDNA) in Escherichia coli using PCR-based amplification and sequential ligation. Here we test this approach by cloning the complete mouse mtDNA. The mtDNA was divided into four to five fragments based on unique restriction enzyme sites and amplified by high-fidelity long-range DNA polymerase. The synthesized fragments were cloned individually to test their toxicity in the E. coli host and then combined sequentially into a vector containing the E. coli R6K origin of DNA replication. The synthetic complete mouse mtDNA clones were replicated stably and faithfully in E. coli when maintained at moderately low copy numbers per cell. The sequence integrity of the synthetic mouse mtDNA clones was confirmed by nucleotide sequencing; no mutations or rearrangements in the genome were found. This approach can facilitate the cloning of entire mammalian mitochondrial genomes in E. coli and assist in the introduction of desired modifications into the mitochondrial genome.

3D model of RNA polymerase and bidirectional transcription

In the in vitro mitochondrial (mt) transcription initiation system with mt RNA polymerase fraction and mt lysate, the transcription initiation products were shown to be synthesized bidirectionally from the only H-strand-promoter (HSP)/L-strand-promoter region (LSP) of the mitochondrial D-loop genome segment. These transcription products ranged between >100 and >800 bp with the purified mitochondrial RNA polymerase fraction, but were larger (>2030-4000 bp) in size with the mitochondrial lysate in both human and mouse. In this brief report, an in vitro reconstituted mitochondrial transcription system purified by affinity chromatography (heparin-Sepharose) from mouse hypotetraploid letter Ehrlich ascites tumor cell mitochondria was shown to initiate transcription bidirectionally from the mitochondrial D-loop region (HSP/LSP), as evidenced by in vitro generated transcription products. The in vitro generated transcription products were separated by sequencing gel. But this in vitro reconstituted transcription system was not studied beyond the D-loop region. A 3D model of the enzyme RNA polymerase was docked with both ATP and CTP.

Transcription factor activator protein-1 expressed by kainate treatment can bind to the non-coding region of mitochondrial genome in murine hippocampus

We have demonstrated previously that the transcription factor activator protein-1 (AP-1) complex is translocated into mitochondria into the nucleus in murine hippocampus after systemic kainate injection (Ogita et al. [2002] J. Neurosci. 22:2561-2570). The present study investigates whether the mitochondrial AP-1 complex translocated in response to kainate treatment binds to AP-1-like sites located at the non-coding region of the mitochondrial genome in mouse hippocampus. There are 10 sites with sequences similar to the nuclear AP-1 site in the non-coding region. Of 10 pieces (MT-1-MT-10) of synthesized double-stranded oligonucleotides, each containing a mitochondrial AP-1-like site, MT-3, MT-4, and MT-9 were effective in inhibiting mitochondrial AP-1 DNA binding enhanced by kainate. Electrophoresis mobility shift analysis using radiolabeled MT-3 and MT-9 probes demonstrated marked enhancement with binding of these 2 probes in hippocampal mitochondrial extracts prepared 2-6 hr after kainate treatment. Unlabeled AP-1 probe was more potent than unlabeled MT-9 probe in inhibiting the mitochondrial MT-9 binding. Supershift analysis revealed participation of particular Fos/Jun family proteins, such as c-Fos, Fos-B, c-Jun, Jun-B, and Jun-D, in MT-9 binding in hippocampal mitochondrial extracts prepared 4 hr after kainate treatment. Immunoprecipitation analysis using anti-c-Fos antibody demonstrated that c-Fos associated with the mitochondrial genome in hippocampal mitochondria prepared from kainate-treated animals. These results suggest that the AP-1 complex expressed by in vivo kainate treatment would bind to AP-1-like sites in the non-coding region of the mitochondrial genome after translocation into mitochondria from murine hippocampus.

Functional domains of chicken mitochondrial transcription factor A for the maintenance of mitochondrial DNA copy number in lymphoma cell line DT40

Nuclear and mitochondrial (mt) forms of chicken mt transcription factor A (c-TFAM) generated by alternative splicing of a gene (c-tfam) were cloned. c-tfam mapped at 6q1.1-q1.2 has similar exon/intron organization as mouse tfam except that the first exons encoding the nuclear and the mt form-specific sequences were positioned oppositely. When cDNA encoding the nuclear form was transiently expressed in chicken lymphoma DT40 cells after tagging at the C terminus with c-Myc, the product was localized into nucleus, whereas the only endogenous mt form of DT40 cells was immunostained exclusively within mitochondria. c-TFAM is most similar to Xenopus (xl-) TFAM in having extended C-terminal regions in addition to two high mobility group (HMG) boxes, a linker region between them, and a C-terminal tail, also found in human and mouse TFAM. Similarities between c- and xl-TFAM are higher in linker and C-terminal regions than in HMG boxes. Disruption of both tfam alleles in DT40 cells prevented proliferation. The tfam+/tfam- cells showed a 50 and 40-60% reduction of mtDNA and its transcripts, respectively. Expression of exogenous wild type c-tfam cDNA in the tfam+/tfam- cells increased mtDNA up to 4-fold in a dose-dependent manner, whereas its transcripts increased only marginally. A deletion mutant lacking the first HMG box lost this activity, whereas only marginal reduction of the activity was observed in a deletion mutant at the second HMG box. Despite the essential role of the C-terminal tail in mtDNA transcription demonstrated in vitro, deletion of c-TFAM at this region reduced the activity of maintenance of the mtDNA level only by 50%. A series of deletion mutant at the tail region suggested stimulatory and suppressive sequences in this region for the maintenance of mtDNA level.

Characterization of the red knot (Calidris canutus) mitochondrial control region

We sequenced the complete mitochondrial control regions of 11 red knots (Calidris canutus). The control region is 1168 bp in length and is flanked by tRNA glutamate (glu) and the gene ND6 at its 5' end and tRNA phenylalanine (phe) and the gene 12S on its 3' end. The sequence possesses conserved sequence blocks F, E, D, C, CSB-1, and the bird similarity box (BSB), as expected for a mitochondrial copy. Flanking tRNA regions show correct secondary structure, and a relative rate test indicated no significant difference between substitution rates in the sequence we obtained versus the known mitochondrial sequence of turnstones (Charadriiformes: Scolopacidae). These characteristics indicate that the sequence is mitochondrial in origin. To confirm this, we sequenced the control region of a single individual using both purified mitochondrial DNA and genomic DNA. The sequences were identical using both methods. The sequence and methods presented in this paper may now serve as a reference for future studies using knot and other avian control regions. Furthermore, the discovery of five variable sites in 11 knots towards the 3' end of the control region, and the variability of this region in contrast to the more conserved central domain in the alignment between knots and other Charadriiformes, highlights the importance of this area as a source of variation for future studies in knots and other birds.

Structure and evolution of the avian mitochondrial control region

The structural and evolutionary characteristics of the mitochondrial control region were studied by using control region sequences of 68 avian species. The distribution of the variable nucleotide positions within the control region was found to be genus specific and not dependant on the level of divergence, as suggested before. Saturation was shown to occur at the level of divergence of 10% in pairwise comparisons of the control region sequences, as has also been reported for the third codon positions in ND2 and cytochrome b genes of mtDNA. The ratio of control region vs cytochrome b divergence in pairwise comparisons of the sequences was shown to vary from 0.13 to 21.65, indicating that the control region is not always the most variable region of the mtDNA, but also that there are differences in the rate of divergence among the lineages. Only two of the conserved sequence blocks localized earlier for other species, D box and CSB-1, were found to show a considerable amount of sequence conservation across the avian and mammalian sequences. Additionally, a novel avian-specific sequence block was found.

Structure and patterns of sequence variation in the mitochondrial DNA control region of the great cats

Mitochondrial DNA control region structure and variation were determined in the five species of the genus Panthera. Comparative analyses revealed two hypervariable segments, a central conserved region, and the occurrence of size and sequence heteroplasmy. As observed in the domestic cat, but not commonly seen in other animals, two repetitive sequence arrays (RS-2 with an 80-bp motif and RS-3 with a 6-10-bp motif) were identified. The 3' ends of RS-2 and RS-3 were highly conserved among species, suggesting that these motifs have different functional constraints. Control region sequences provided improved phylogenetic resolution grouping the sister taxa lion (Panthera leo) and leopard (Panthera pardus), with the jaguar (Panthera onca).

The genetic relationships of two subspecies of striped field mice, Apodemus agrarius coreae and Apodemus agrarius chejuensis

We obtained 282 base pairs of sequence for the mitochondrial control region of 70 individuals of Korean striped field mice Apodemus agrarius coreae and Apodemus agrarius chejuensis to determine the levels of genetic divergence between these morphologically distinct taxa. The DNA sequences showed more genetic diversity (pi) in A. agrarius coreae (2.98%) than in A. agrarius chejuensis (1.86%). Our data do not support the current concept that the two morphotypes are different species, but phylogenetic analysis indicates that animals of A. agrarius coreae with large body size from Wan Island cluster with the large-bodied A. agrarius chejuensis, and should be included in that taxon. As currently accepted A. agrarius coreae is not strictly monophyletic, because the large-bodied samples cluster within the range of mitochondrial variation of A. agrarius chejuensis. The fact that the two morphotypes do not share mitochondrial haplotypes (chi2=66, P < 0.001) suggests that there is little gene flow between them. A molecular clock estimate suggests that the two subspecies might have been isolated at the time of separation of the islands from the mainland.

Replicative advantage and tissue-specific segregation of RR mitochondrial DNA between C57BL/6 and RR heteroplasmic mice

To investigate the interactions between mtDNA and nuclear genomes, we produced heteroplasmic maternal lineages by transferring the cytoplasts between the embryos of two mouse strains, C57BL/6 (B6) and RR. A total of 43 different nucleotides exist in the displacement-loop (D-loop) region of mtDNA between B6 and RR. Heteroplasmic embryos were reconstructed by electrofusion using a blastomere from a two-cell stage embryo of one strain and an enucleated blastomere from a two-cell stage embryo of the other strain. Equivalent volumes of both types of mtDNAs were detected in blastocyst stage embryos. However, the mtDNA from the RR strain became biased in the progeny, regardless of the source of the nuclear genome. The RR mtDNA population was very high in most of the tissues examined but was relatively low in the brain and the heart. An age-related increase of RR mtDNA was also observed in the blood. The RR mtDNAs in the reconstructed embryos and in the embryos collected from heteroplasmic mice showed a different segregation pattern during early embryonic development. These results suggest that the RR mtDNA has a replicative advantage over B6 mtDNA during embryonic development and differentiation, regardless of the type of nuclear genome.

Human mitochondrial transcription factor A and promoter spacing integrity are required for transcription initiation

The two major promoters for transcription of the human mitochondrial genome are located near each other in the displacement-loop region of the molecule. Previous work has localized these promoters to regions of < 100 nucleotides each; the DNA sequence at the transcription start site is stringently required, as is the region from -10 to -40 base pairs upstream of each respective start site. Each upstream site is recognized and bound by human mitochondrial transcription factor A (h-mtTFA), an event previously shown to be important for transcriptional activation. We report here results using recombinant h-mtTFA that demonstrate the dependence of transcription initiation of h-mtTFA. In addition, altering the distance between the h-mtTFA binding site and the transcription start site greatly impairs transcription initiation efficiency. The decrease in transcription initiation efficiency was shown to be a consequence of altering the position of h-mtTFA binding as opposed to the strength of h-mtTFA binding, as judged by DNA footprinting ability. Analysis of a chimeric yeast-human promoter revealed that the yeast mtTFA homologue cannot substitute for the human protein, even when bound at an appropriate position upstream of the human transcription start site.

Identification of a stable RNA encoded by the H-strand of the mouse mitochondrial D-loop region and a conserved sequence motif immediately upstream of its polyadenylation site

By using a combination of Northern blot hybridization with strand-specific DNA probes, S1 nuclease protection, and sequencing of oligo-dT-primed cDNA clones, we have identified a 0.8 kb poly(A)-containing RNA encoded by the H-strand of the mouse mitochondrial D-loop region. The 5' end of the RNA maps to nucleotide 15417, a region complementary to the start of tRNA(Pro) gene and the 3' polyadenylated end maps to nucleotide 16295 of the genome, immediately upstream of tRNA(Phe) gene. The H-strand D-loop region encoded transcripts of similar size are also detected in other vertebrate systems. In the mouse, rat, and human systems, the 3' ends of the D-loop encoded RNA are preceded by conserved sequences AAUAAA, AAUUAA, or AACUAA, that resemble the polyadenylation signal. The steady-state level of the RNA is generally low in dividing or in vitro cultured cells, and markedly higher in differentiated tissues like liver, kidney, heart, and brain. Furthermore, an over 10-fold increase in the level of this RNA is observed during the induced differentiation of C2C12 mouse myoblast cells into myotubes. These results suggest that the D-loop H-strand encoded RNA may have yet unknown biological functions. A 20 base pair DNA sequence from the 3' terminal region containing the conserved sequence motif binds to a protein from the mitochondrial extracts in a sequence-specific manner. The binding specificity of this protein is distinctly different from the previously characterized H-strand DNA termination sequence in the D-loop or the H-strand transcription terminator immediately downstream of the 16S rRNA gene. Thus, we have characterized a novel poly(A)-containing RNA encoded by the H-strand of the mitochondrial D-loop region and also identified the putative ultimate termination site for the H-strand transcription.

Decreased mitochondrial gene expression in isolated islets of rats injected neonatally with streptozotocin

The aim of the present study was to evaluate the possible role of the expression of the mitochondrial genome for the regulation of insulin production in the pancreatic Beta cell. For this purpose, islets of Langerhans were isolated from adult control rats and rats injected neonatally with streptozotocin and the islet contents of specific mitochondrial DNAs and RNAs together with nuclear-encoded RNAs were determined. The contents of mitochondrial cytochrome b mRNA, the mitochondrial 12 S rRNA and insulin mRNA were all 30-40% lower in islets isolated from the streptozotocin-treated rats as compared to islets from control rats. In contrast, the nuclear mRNA coding for the mitochondrial adenine nucleotide translocator was not decreased in the streptozotocin-treated rats. Contents of mitochondrial DNA, as assessed by the Southern blotting technique, were markedly decreased in the streptozotocin islets. Sequence analysis of mitochondrial DNA from streptozotocin islets and control islets however, did not reveal any differences in nucleotide sequences. In control islets the contents of mitochondrial cytochrome b mRNA increased in response to a high glucose concentration during a 4-h incubation period. Serum deprivation or the addition of theophylline or 4-phorbol 12-myristate 13-acetate failed to affect the cytochrome b mRNA contents in vitro. It is concluded that islets of streptozotocin-treated rats contain low contents of mitochondrial DNA and RNA. Since a lower mitochondrial RNA content may result in a diminished oxidative capacity, it is conceivable that a deficiency of this messenger may contribute to the development of insulin deficiency.

The complete and symmetric transcription of the main non coding region of rat mitochondrial genome: in vivo mapping of heavy and light transcripts

The experiments here reported demonstrate that the main non-coding region of rat mitochondrial DNA is symmetrically transcribed. We have identified stable heavy and light transcripts, whose pattern is rather complex, in the D-loop region of rat mitochondrial DNA. Their relative concentrations have been determined. We detected heavy transcripts which encompass the whole D-loop and more abundant heavy RNA species which we interpreted as transcripts terminating downstream of the 3' end of the last coded gene (Thr-tRNA). The processed heavy RNA species contain polyA, suggesting a strict association between cleavage and polyadenylation. The pattern of light transcripts shows a long RNA, which, starting from the light strand promoter, covers the whole segment, and shorter RNA species which seems to be actively processed at the level of the conserved sequence boxes, probably acting as primers. The symmetric transcription of the D-loop containing region of rat mitochondrial DNA, and in particular the presence of stable transcripts complementary to the putative RNA primers, suggest that mechanisms mediated by interaction between complementary transcripts (antisense RNAs) might play a role in the regulation of mitochondrial DNA replication and expression.

Molecular population genetics of mtDNA size variation in crickets

Nucleotide sequence analysis of a region of cricket (Gryllus firmus) mtDNA showing discrete length variation revealed tandemly repeated sequences 220 base pairs (bp) in length. The repeats consist of 206 bp sequences bounded by the dyad symmetric sequence 5'GGGGGCATGCCCCC3'. The sequence data showed that mtDNA size variation in this species is due to variation in the number of copies of tandem repeats. Southern blot analysis was used to document the frequency of crickets heteroplasmic for two or more different-sized mtDNAs. In New England populations of G. firmus and a close relative Gryllus pennsylvanicus approximately 60% of the former and 45% of the latter were heteroplasmic. From densitometry of autoradiographs the frequencies of mtDNA size classes were determined for the population samples and are shown to very different in the two species. However, in populations where hybridization between the two species has occurred, the frequencies of size classes and cytoplasmic genotypes in each species' distinct mtDNA lineage were shifted in a manner suggesting nuclear-cytoplasmic interactions. The data were applied to reported diversity indices and hierarchical statistics. The hierarchical statistics indicated that the greatest proportion of variation for mtDNA size was due to variation among individuals in their cytoplasmic genotypes (heteroplasmic or homoplasmic state). The diversity indices were used to estimate a per-generation mutation rate for size variants of 10(-4). The data are discussed in light of the relationship between genetic drift and mutation in maintaining variation for mtDNA size.

Nuclear gene products that function in mitochondrial DNA replication

Mammalian mitochondrial DNAs replicate unidirectionally from two distinct strand-specific origins. A round of replication begins at the heavy-strand origin (the D-loop) where transcripts from an upstream promoter serve as the primers for DNA synthesis. The transition from RNA to DNA synthesis occurs within short, conserved nucleotide sequence blocks and is mediated by specific endonucleolytic cleavage of the primary transcript. An enzymic component involved in the generation of primer RNA in mouse mitochondria has been identified. It is a sequence-specific endoribonuclease that cleaves single-stranded RNA substrate precisely at one of the transition sites. The other origin, that for light-strand synthesis, is located well apart on the genome and functions only when in a single-stranded template form. This origin has a defined secondary structure that is the most highly conserved sequence element in mammalian mitochondrial DNAs. Initiation of replication at this origin is by the action of a mitochondrial DNA primase, which is capable of synthesizing a short stretch of ribonucleotides before switching to DNA synthesis. Mitochondrial DNA primase appears to have an associated RNA species and the evidence to date suggests that components of both the D-loop endoribonuclease and the DNA primase are nuclear gene products.

Yeast mitochondrial RNA polymerase is homologous to those encoded by bacteriophages T3 and T7

Analysis of the nucleotide sequence of the genetic locus for yeast mitochondrial RNA polymerase (RPO41) reveals a continuous open reading frame with the coding potential for a polypeptide of 1351 amino acids, a size consistent with the electrophoretic mobility of this enzymatic activity. The transcription product from this gene spans the singular reading frame. In vivo transcript abundance reflects codon usage and growth under stringent conditions for mitochondrial biogenesis and function results in a several fold higher level of gene expression than growth under glucose repression. A comparison of the yeast mitochondrial RNA polymerase amino acid sequence to those of E. coli RNA polymerase subunits failed to demonstrate any regions of homology. Interestingly, the mitochondrial enzyme is highly homologous to the DNA-directed RNA polymerases of bacteriophages T3 and T7, especially in regions most highly conserved between the T3 and T7 enzymes themselves.

Promoter selection in human mitochondria involves binding of a transcription factor to orientation-independent upstream regulatory elements

Selective transcription of human mitochondrial DNA requires a transcription factor (mtTF) in addition to an essentially nonselective RNA polymerase. Partially purified mtTF is able to sequester promoter-containing DNA in preinitiation complexes in the absence of mitochondrial RNA polymerase, suggesting a DNA-binding mechanism for factor activity. Functional domains, required for positive transcriptional regulation by mtTF, are identified within both major promoters of human mtDNA through transcription of mutant promoter templates in a reconstituted in vitro system. These domains are essentially coextensive with DNA sequences protected from nuclease digestion by mtTF-binding. Comparison of the sequences of the two mtTF-responsive elements reveals significant homology only when one sequence is inverted; the binding sites are in opposite orientations with respect to the predominant direction of transcription. Thus mtTF may function bidirectionally, requiring additional protein-DNA interactions to dictate transcriptional polarity. The mtTF-responsive elements are arrayed as direct repeats, separated by approximately 80 bp within the displacement-loop region of human mitochondrial DNA; this arrangement may reflect duplication of an ancestral bidirectional promoter, giving rise to separate, unidirectional promoters for each strand.