Transcription mapping of the Ori L region reveals novel precursors of mature RNA species and antisense RNAs in rat mitochondrial genome

We have identified new transcripts in the region surrounding the L-strand replication origin (Ori L) of rat liver mitochondrial DNA. In particular, we have detected previously unidentified intermediates of RNA processing on both the heavy and the light strands, such as precursors of the ND2 mRNA plus the Trp-tRNA and precursors of the tRNAs clustered in the Ori L region. This indicates that the mechanism of RNA processing in mitochondria proceeds step-wise producing a variety of precursors of the mature forms. The other striking finding is the detection of antisense RNA species in the region of L-strand replication. Since a variety of antisense transcripts were also found in the D-loop region of rat mitochondrial DNA, we suggest that they might play a regulatory role in the replication and expression of the mitochondrial genome.

Mitochondrial DNA detection and copy number determination in the spermatozoa of the sea urchin Arbacia lixula

The Polymerase Chain reaction technique has been used in order to detect and amplify a specific region of mtDNA, in a total DNA preparation extracted from the sperm of the sea urchin Arbacia lixula. The amplified fragment is the D-loop region which hybridizes with the homologous region extracted from the egg mtDNA. The results demonstrate that mtDNA is present in sperm cell, and, since the replication origin is present it is potentially able to replicate in the zygote. Furthermore, the technique used allowed us to estimate mtDNA copy number in sea urchin sperm, which has never been done before. Our results are that sea urchin sperm cell contains between 4 and 28 mtDNA molecules.

A simple method for global sequence comparison

A simple method of sequence comparison, based on a correlation analysis of oligonucleotide frequency distributions, is here shown to be a reliable test of overall sequence similarity. The method does not involve sequence alignment procedures and permits the rapid screening of large amounts of sequence data. It identifies those sequences which deserve more careful analysis of sequence similarity at the level of resolution of the single nucleotide. It uses observed quantities only and does not involve the adoption of any theoretical model.

The branching order of mammals: phylogenetic trees inferred from nuclear and mitochondrial molecular data

In order to clarify some controversial phylogenies such as those regarding the triplet of human, rodent, and cow and the evolutionary position of Lagomorpha with respect to other mammals, we have analyzed both nuclear and mitochondrial genes using the stationary Markov model developed in our laboratory. We found that the two sets of genes give different results. In particular the mitochondrial tree showed rabbit linked first to rodents and the rabbit-rodents branch linked to artiodactyls with human as the outgroup. The most favorite nuclear tree showed human linked first to artiodactyls and the human-artiodactyls branch linked to rabbit with rodents as the outgroup. The obvious questions, (1) which tree is the correct one, or (2) both trees can be incorrect, and (3) how can we explain such an evolutionary pattern, are discussed on the basis of our limited knowledge of factors that influence the clocklike behavior of biological macromolecules.

Mitochondrial DNA in the sea urchin Arbacia lixula: nucleotide sequence differences between two polymorphic molecules indicate asymmetry of mutations

Two polymorphic forms of mitochondrial DNA (mtDNA) extracted from Arbacia lixula eggs were cloned and the nucleotide sequences of specific regions determined. A comparison of the sequences of the sense strand of the two molecules demonstrates that all the differences are transitions and only of the A----G type. A change such as G----A (or A----G) on the sense mtDNA strand results from either a direct G----A (or A----G) mutation on that strand or a C----T (or T----C) on the complementary strand. None of the C----T (or T----C) changes were detected on the sense strand, which implies that the A----G mutation bias on the sense strand is not reversed for the other strand. Our observation indicates the existence of mechanisms acting asymmetrically on the two mtDNA strands, possibly during mtDNA replication.

Mitochondrial DNA in the sea urchin Arbacia lixula: evolutionary inferences from nucleotide sequence analysis

From the stirodont Arbacia lixula we determined the sequence of 5,127 nucleotides of mitochondrial DNA (mtDNA) encompassing 18 tRNAs, two complete coding genes, parts of three other coding genes, and part of the 12S ribosomal RNA (rRNA). The sequence confirms that the organization of mtDNA is conserved within echinoids. Furthermore, it underlines the following peculiar features of sea urchin mtDNA: the clustering of tRNAs, the short noncoding regulatory sequence, and the separation by the ND1 and ND2 genes of the two rRNA genes. Comparison with the orthologous sequences from the camarodont species Paracentrotus lividus and Strongylocentrotus purpuratus revealed that (1) echinoids have an extra piece on the amino terminus of the ND5 gene that is probably the remnant of an old leucine tRNA gene; (2) third-position codon nucleotide usage has diverged between A. lixula and the camarodont species to a significant extent, implying different directional mutational pressures; and (3) the stirodont-camarodont divergence occurred twice as long ago as did the P. lividus-S. purpuratus divergence.

The main regulatory region of mammalian mitochondrial DNA: structure-function model and evolutionary pattern

The evolution of the main regulatory region (D-loop) of the mammalian mitochondrial genome was analyzed by comparing the sequences of eight mammalian species: human, common chimpanzee, pygmy chimpanzee, dolphin, cow, rat, mouse, and rabbit. The best alignment of the sequences was obtained by optimization of the sequence similarities common to all these species. The two peripheral left and right D-loop domains, which contain the main regulatory elements so far discovered, evolved rapidly in a species-specific manner generating heterogeneity in both length and base composition. They are prone to the insertion and deletion of elements and to the generation of short repeats by replication slippage. However, the preservation of some sequence blocks and similar cloverleaf-like structures in these regions, indicates a basic similarity in the regulatory mechanisms of the mitochondrial genome in all mammalian species. We found, particularly in the right domain, significant similarities to the telomeric sequences of the mitochondrial (mt) and nuclear DNA of Tetrahymena thermophila. These sequences may be interpreted as relics of telomeres present in ancestral linear forms of mtDNA or may simply represent efficient templates of RNA primase-like enzymes. Due to their peculiar evolution, the two peripheral domains cannot be used to estimate in a quantitative way the genetic distances between mammalian species. On the other hand the central domain, highly conserved during evolution, behaves as a good molecular clock. Reliable estimates of the times of divergence between closely and distantly related species were obtained from the central domain using a Markov model and assuming nonhomogeneous evolution of nucleotide sites.

Glutamine synthetase gene evolution: a good molecular clock

Glutamine synthetase (EC 6.3.1.2) gene evolution in various animals, plants, and bacteria was evaluated by a general stationary Markov model. The evolutionary process proved to be unexpectedly regular even for a time span as long as that between the divergence of prokaryotes from eukaryotes. This enabled us to draw phylogenetic trees for species whose phylogeny cannot be easily reconstructed from the fossil record. Our calculation of the times of divergence of the various organelle-specific enzymes led us to hypothesize that the pea and bean chloroplast genes for these enzymes originated from the duplication of nuclear genes as a result of the different metabolic needs of the various species. Our data indicate that the duplication of plastid glutamine synthetase genes occurred long after the endosymbiotic events that produced the organelles themselves.

Evolutionary analysis of the nucleus-encoded subunits of mammalian cytochrome c oxidase

The cytochrome c oxidase enzyme complex of eukaryotes is made up of three mitochondrial-coded subunits and a variable number of nuclear-coded subunits. Some nuclear-coded subunits are present in multiple forms and probably perform a tissue- or development-specific function. A detailed evolutionary analysis of the cytochrome c oxidase subunits that have been sequenced to date is reported here. We have found that gene duplication events from which the liver and heart isoforms of rat subunits VIa and subunit VIII originated can both be dated at about 240 +/- 90 million years ago, long before the radiation of mammalian lineages. Sequence divergence between the processed-type pseudogenes for the subunits IV, VIc and VIII have been estimated. Our results indicate that they arose fairly recently, thus suggesting that retroposition is a continuing process. We show that the rate of silent substitution in mitochondrial-coded subunits is 5-10 times higher than in nuclear-coded subunits; on the other hand replacement rates, although differing from gene to gene, are roughly of the same order of magnitude in both nuclear and mitochondrial genes. In the case of most of the nuclear-coded proteins we observed a slightly greater similarity between rats and cow, which agrees with the data obtained for mitochondrial-coded subunits.

Direct evidence that restriction endonucleases may under estimate the degree of divergence between molecules

We studied two polymorphic forms of mtDNA extracted from A. lixula eggs. In order to compare and to quantitate the variability, we sequenced specific regions of the two molecules. In this way, we obtained a precise measurement of the variability within two haplotypes. We also obtained a direct demonstration that some differences in nucleotide sequence can escape detection when restriction endonuclease analysis is used. Our results underline the unreliability of the use of restriction mapping to estimate divergence between relatively short and closely related DNA sequences.

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.

Sequence-dependent DNA curvature: conformational signal present in the main regulatory region of the rat mitochondrial genome

Theoretical analysis and experimental approaches by gel electrophoresis in retarding conditions allowed us to identify the presence of an intrinsic bending in the D-loop containing region of the rat mitochondrial genome. The curvature was located in the right domain of the sequence analyzed, between the origin of replication of the heavy strand and its promoter. The preliminary evidence of a specific recognition of the bent DNA with mitochondrial matrix proteins suggests a probable role of this DNA conformation in the duplication and/or expression of the mammalian mitochondrial genome.

DNA microenvironments and the molecular clock

A few years ago we presented a stationary Markov model of gene evolution according to which only homologous genes from not too divergent species obeying the condition of being stationary may behave as reliable molecular clocks. A compartmentalized model of the nuclear genome in which the genes are distributed in compartments, the isochores, defined by their G + C content has been proposed recently. We have found that only homologous gene pairs that are stationary, and belong to the same isochore, can be used consistently for the determination of phylogeny and base substitution rate. In particular, for the rodent-human couple, only about half of the homologous gene pairs are stationary. Stationary genes evolve at the third silent codon position with the same velocity independent of the genes and base composition. By contrast, nonstationary genes display apparent rate values (pseudovelocities) that are significantly higher. Our results cast doubt upon recent claims of a large acceleration in the rate of molecular evolution in rodents.

The complete nucleotide sequence, gene organization, and genetic code of the mitochondrial genome of Paracentrotus lividus

The 15,697-nucleotide sequence of Paracentrotus lividus mitochondrial DNA is reported. This genome codes for 2 rRNAs, 22 tRNAs, and 12 mRNAs which specify 13 subunits of the mitochondrial inner membrane respiratory complexes. The gene arrangement differs from that of other animal species. The two ribosomal genes 16 S and 12 S are separated by a stretch of about 3.3 kilobase pairs which contains the ND1 and ND2 genes and a cluster of 15 tRNA genes. The ND4L coding sequence is not contained in the ND4 mRNA but has its own mRNA which maps between the tRNA(Arg) and the Co II genes. The main noncoding region, located in the tRNA gene cluster, is only 132 nucleotides long, but contains sequences homologous to the mammalian displacement loop. Other short noncoding sequences are interspersed in the genome: they contain a conserved AT consensus which probably has a role in transcription or RNA processing. As regards the mitochondrial genetic code, the codons AGA and AGG specify serine and are recognized by a tRNA with a GCU anticodon, whereas AUA and AAA code for isoleucine and asparagine rather than for methionine and lysine. Except for ND4L which starts with AUC and ATPase 8 which starts with GUG, AUG is used as the initiation codon. In 11 out of 13 cases the genes terminate with the canonical stop codons UAA or UAG. These observations suggest that during invertebrate evolution each lineage developed its own mechanism of mitochondrial DNA replication and transcription and of RNA processing and translation.

The complete nucleotide sequence of the Rattus norvegicus mitochondrial genome: cryptic signals revealed by comparative analysis between vertebrates

This paper reports the nucleotide sequence of rat mitochondrial DNA, only the fourth mammalian mitochondrial genome to be completely sequenced. Extensive comparative studies performed with similar genomes from other organisms revealed a number of interesting features. 1) Messenger RNA genes: the codon strategy is mainly dictated by the base compositional constraints of the corresponding codogenic DNA strand. The usage of the initiation and termination codons follows well-established rules. In general the canonical initiator, ATG, and terminators, TAA and TAG (in rat, only TAA), are always present when there is gene overlapping or when the mRNAs possess untranslated nucleotides at the 5' or 3' ends. 2) Transfer RNA genes: a number of features suggest the peculiar evolutionary behavior of this class of genes and confirm their role in the duplication and rearrangement processes that took place in the evolution of the animal mitochondrial genome. 3) Ribosomal RNA genes: accurate sequence analysis revealed a number of significant examples of complementarity between ribosomal and messenger RNAs. This suggests that they might play an important role in the regulation of mitochondrial translation and transcription mechanisms. The properties revealed by our work shed new light on the organization and evolution of the vertebrate mitochondrial genome and more importantly open up the way to clearly aimed experimental studies of the regulatory mechanisms in mitochondria.

Mitochondrial genome in animal cells. Structure, organization, and evolution

In the past decade, the development of new DNA, RNA, and protein technologies has greatly incremented the knowledge about the organization and expression of mitochondrial DNA. The complete base sequence of mitochondrial DNA of several animals is known and many data are rapidly accumulating on the mitochondrial genomes of other systems. Here we discuss the results so far obtained that disclosed unexpected features of mitochondrial genetics. Furthermore, mitochondrial DNA has become established as a powerful tool for evolutionary studies in animals. Evidences are presented demonstrating that the evolution of mitochondrial DNA has proceeded in different ways in the various taxonomic groups. Data on heteroplasmic animals, which demonstrate the rapid evolution of mitochondrial DNA, are also presented.

Clustering of tRNA genes in Paracentrotus lividus mitochondrial DNA

We have determined the base sequence of the restriction fragment Bam1-2 (3,593) of Paracentrotus lividus (sea urchin) mtDNA. This fragment contains, in addition to genes previously identified (part of the 12S rRNA, ND1 and part of the ND2 mRNA), a cluster of 15 tRNA genes located between the 12S and ND1 genes. Also to be found in the tRNA gene cluster, between the tRNA(Thr) and tRNA(Pro) genes, is a sequence of 134 bp which constitutes the only non-coding region of this DNA so far identified. The distinctive organization of the tRNA genes and the extreme size reduction of the non-coding region suggest the existence of unique mechanisms for the regulation of gene expression in this organism.

Structural elements highly preserved during the evolution of the D-loop-containing region in vertebrate mitochondrial DNA

A detailed comparative study of the regions surrounding the origin of replication in vertebrate mitochondrial DNA (mtDNA) has revealed a number of interesting properties. This region, called the D-loop-containing region, can be divided into three domains. The left (L) and right (R) domains, which have a low G content and contain the 5' and the 3' D-loop ends, respectively, are highly variable for both base sequence and length. They, however, contain thermodynamically stable secondary structures which include the conserved sequence blocks called CSB-1 and TAS which are associated with the start and stop sites, respectively, for D-loop strand synthesis. We have found that a "mirror symmetry" exists between the CSB-1 and TAS elements, which suggests that they can act as specific recognition sites for regulatory, probably dimeric, proteins. Long, statistically significant repeats are found in the L and R domains. Between the L and R domains we observed in all mtDNA sequences a region with a higher G content which was apparently free of complex secondary structure. This central domain, well preserved in mammals, contains an open reading frame of variable length in the organisms considered. The identification of common features well preserved in evolution despite the high primary structural divergence of the D-loop-containing region of vertebrate mtDNA suggests that these properties are of prime importance for the mitochondrial processes that occur in this region and may be useful for singling out the sites on which one should operate experimentally in order to discover functionally important elements.

A simple quantitative model of the molecular clock

We present the ideas, and their motivation, at the basis of a simple model of nucleic acid evolution: the stationary Markov process, or Markov clock. After a brief review of its relevant mathematical properties, the Markov clock is applied to nucleotide sequences from mitochondrial and nuclear genes of different species. Particular emphasis is given to the necessity of carrying out a correct statistical analysis, which allows us to check quantitatively the applicability of our model. We find evidence that the Markov clock ticks in many different processes, and that its limitations can be understood in terms of a simple idea that we call the "base-drift" hypothesis. This hypothesis correlates the deviations from the stationarity of the Markov process to the evolutionary distance dAB(p) of two species A and B, relative to the process P. We conclude by discussing the implications of our findings for future work.

Duplication and remoulding of tRNA genes during the evolutionary rearrangement of mitochondrial genomes

During the evolution of sea urchins, a transfer RNA gene lost its tRNA function and became part of a protein-coding gene. This functional loss of a tRNA with specificity for one group of leucine codons (CUN, where N is any base) was accompanied by the gain of a new tRNA with that specificity. The new tRNA gene for CUN codons appears to have evolved by duplication and divergence from a tRNA gene specific for another group of leucine codons (UUR, where R is a purine). These proposals account for (1) the strong sequence resemblance between the modern tRNA genes for CUN and UUR codons in Paracentrotus, (2) the altered location of the CUN gene in mitochondrial DNA of this urchin, and (3) the persistence of a 72-base pair sequence containing a trace of the old CUN gene at its original location. The old CUN gene now codes for an extra 24 amino acids at the amino end of subunit 5 in NADH dehydrogenase. Besides giving clues about the mechanisms by which tRNA genes move during mitochondrial DNA evolution, this finding leads us to propose a pathway relating the arrangements of other genes in mitochondrial DNAs from four animal phyla.

A novel gene order in the Paracentrotus lividus mitochondrial genome

The mitochondrial DNA (mtDNA) from Paracentrotus lividus (sea urchin) eggs, a circular molecule of about 15,500 bp, has been cloned in plasmid vectors after cleavage with various restriction enzymes. By a combination of Northern blot hybridization and nucleotide sequence analysis we have characterized most of the P. lividus mitochondrial transcripts and determined the basic gene organization of the mtDNA. The nucleotide sequence of a gene for one NADH dehydrogenase (ND) subunit, ND4L, has also been determined. Our results show the existence of a novel gene order. The 12S and 16S rRNA genes are not contiguous but are separated from each other by ND1 and ND2 genes. The ND4L gene is not adjacent to ND4 but is located between the tRNAArg gene and the gene for subunit II of cytochrome oxidase (CoII). The tRNA genes are reshuffled and contrary to all vertebrate mitochondrial genomes studied so far, there are no intergenic regions between the tRNAPhe and the cytochrome b genes. These characteristics suggest a peculiar mechanism for the regulation of gene expression in this organism and provide information on the evolution of the mitochondrial genetic system in animal cells.