The evolution of mitochondrial DNA.

Although the massive sequencing of mitochondrial DNA from various organisms, together with studies of a different nature, has contributed enormously to the knowledge of the organization and function of this cytoplasmic genome, many issues, mainly the relationships with the nuclear genome, remain unsolved. This review critically evaluates the most recent advances in research on the evolution of the mitochondrial DNA from a qualitative and quantitative point of view, underlining the multiplicity of structures and genetic organization of this genome, which contrasts with its reduced, but rather constant, information content in various organisms. It also highlights the role that mitochondrial DNA is now playing, particularly in metazoans, in different disciplines and application fields. Among these, particular attention is focused on the discovery of the mitochondrial origin of several diseases affecting primarily the neuromuscular system.

Further analysis of the type differences of rat-liver mitochondrial DNA.

1. The sequence of the small Hind III fragment F of rat-liver mitochondrial DNA (mtDNA) type A and type B was determined in order to investigate the nature of the differences between the two types of mtDNA and to determine its position in the Hind III fragment map. 2. The three differences found were point mutations. No deletion or insertion and no modification was observed. Two of the three differences affect the sequences which are recognition sites for Eco RI, Alu I and Taq I in type A, but not in type B, mtDNA. 3. The presence of an Eco RI restriction site only within the Hind III fragment F of type A shows that the fragment is situated in between the Hind III fragments A and E. 4. In one of the six reading frames, the Hind III fragment F contains the code for a carboxyl-terminal end of a polypeptide in which the three mutations do not lead to alterations in the possible aminoacid sequence. 5. The restriction sites for Taq I and a number of the sites for Alu I and Hae III were mapped. 6. The positions of the Hap II fragment J, and of a Hind III fragment G on the mtDNA were determined.

Transport of proteins into mitochondria.

There is still much that is obscure concerning the transport of proteins into or through the mitochondrial membrane systems. In addition, as pointed out previously, it is unlikely that the details of the process are the same for proteins destined for different compartments of the organelle. A brief summary of the process for matrix proteins might be as follows: The proteins are synthesized on free polysomes as precursors of higher molecular weight than the native forms. These precursors are liberated into the cell cytosol and subsequently translocated into the mitochondria. This timing might be different in yeast under some circumstances, synthesis being completed in association with the mitochondria. The precursors interact with a receptor in the outer mitochondrial membrane interaction being mediated by the presequences of the precursors. The presequences therefore act as addressing signals as well as possibly playing a role in one or all of (a) solubilization of precursors, (b) prevention of premature assembly into multimeric structures, or (c) maintenance of nonnative configurations required for transport. Interaction occurs with a second receptor, this time in the inner membrane of the mitochondria, interaction being with multiple sites in the polypeptide chain. Transport across the inner membrane then occurs, this transport depending on a transmembrane electrochemical gradient of which the proton component is the essential part. Transport is accompanied or followed by proteolysis of the prepiece, and formation of the native structure. While steps 1 and 2 of this sequence can be considered well established, the remaining steps are still poorly understood or purely hypothetical. Nevertheless, this sequence of events is consistent with known facts about the process and provides a framework for future investigations.

Structural conservation and variation in the D-loop-containing region of vertebrate mitochondrial DNA.

The nucleotide sequences of the D-loop-containing regions of three rat mitochondrial DNAs (mtDNAs), two from the species Rattus norvegicus and one from R. rattus, were determined. Comparisons made among these sequences and with the mouse sequence showed that, on the basis of both base composition and frequency of nucleotide alterations, three domains could be defined within the D-loop-containing region: a central conserved segment, poor in L-strand adenine, flanked by two divergent, adenine-rich regions. Deletions and insertions were found to occur at an unexpectedly high frequency in these sequences and the conserved sequence block called CSB-1 was found not to be intact in the R. rattus sequence. Although in comparisons of more distantly related mtDNAs the D-loop region is the most divergent on the molecule, it does not diverge more than typical protein genes between R. norvegicus and R. rattus, and its central conserved domain appears to be one of the molecule's most conserved regions. The most variable domain borders the tRNAPhe gene and contains the L and H-strand promoters and the 5' terminus for H-strand DNA synthesis. Within this region we have found sequences in all the mtDNAs we have examined, including those of human, two artiodactyls and Xenopus, that are capable of folding into cloverleaf structures. In the other divergent domain of the same mtDNAs, we find sequences capable of assuming similar secondary structural configurations at or near the sites for the termination of D-loop DNA synthesis. The evolutionary preservation of the potential to form such structures despite the high primary-structural divergence of the regions they occur in, suggests the structures are of principal importance for some processes occurring in the D-loop-containing region.

A novel method for estimating substitution rate variation among sites in a large dataset of homologous DNA sequences.

We present here a novel method to estimate the site-specific relative variability in large sets of homologous sequences. It is based on the simple idea that the more closely related are the compared sequences, the higher the probability of observing nucleotide changes at rapidly evolving sites. A simulation study has been carried out to support the reliability of the method, which has been applied also to analyzing the site variability of all available human sequences corresponding to the two hypervariable regions of the mitochondrial D-loop.

Evolution of T-cell receptor gamma and delta constant region and other T-cell-related proteins in the human-rodent-artiodactyl triplet.

In this paper we report a detailed comparative and evolutionary analysis of the sequences of constant T-cell receptor (Tcr) C gamma delta genes of artiodactyls compared to the homologous sequences of rodents and primates. Because of the frequency and physiological distribution of gamma delta T-cells in different animals, rodents and humans are defined as "gamma delta low" species and ruminants as "gamma delta high" species. Such a characteristic seems to be due to an adaptive role of gamma delta T-cell function. By analyzing the ruminant gene phylogeny of Tcr C gamma we were able to estimate the distance between cattle and sheep at 18 million years ago, a time that is in agreement with other nonmolecular estimates. For Tcr C gamma delta genes a peculiar phylogenetic relationship was found, with human and mouse clustering together and leaving artiodactyls apart. By using appropriate outgroups, the same phylogenetic pattern was obtained with other T-cell related sequences: namely, Tcr C alpha chain, CD3 gamma and delta invariant subunits. Interleukin-2. Interleukin-2 receptor alpha chain and Interleukin-1 beta with the exception of Tcr C beta chain and Interleukin-1 alpha. In contrast, the analysis of all other T-cell nonrelated genes, available in primary databases reveals a different tree, where primates and artiodactyls are sister taxa and rodents are apart in accordance with the current view of mammalian phylogeny. These data are relevant to important evolutionary issues. They show how misleading a phylogeny based on a single or on a few homologous genes may be. In addition they demonstrate that genes with correlated functions may evolve in a lineage specific manner probably in relation to environmental conditions.

Evolutionary origin of nonuniversal CUGSer codon in some Candida species as inferred from a molecular phylogeny.

CUG, a universal leucine codon, has been reported to be read as serine in various yeast species belonging to the genus Candida. To gain a deeper insight into the origin of this deviation from the universal genetic code, we carried out a phylogenetic analysis based on the small-subunit ribosomal RNA genes from some Candida and other related Hemiascomycetes. Furthermore, we determined the phylogenetic relationships between the tRNA(Ser)CAG, responsible for the translation of CUG, from some Candida species and the other serine and leucine isoacceptor tRNAs in C. cylindracea. We demonstrate that the group of Candida showing the genetic code deviation is monophyletic and that this deviation could have originated more than 150 million years ago. We also describe how phylogenetic analysis can be used for genetic code predictions.

The peculiar evolution of apolipoprotein(a) in human and rhesus macaque.

Apo(a) is a low density lipoprotein homologous to plasminogen and has been shown to be involved in coronary atherosclerosis. In the present paper we will try to analyze the interesting evolutionary pattern of Apo(a). The plasminogen gene contains 5 cysteine-rich sequences, called kringles, followed by a protease domain. Apo(a), probably arisen by duplication of an ancestral plasminogen gene, contains many tandemly repeated copies of a sequence domain similar to the fourth kringle of plasminogen, 37 in human and at least 10 in the partially sequenced gene of rhesus, and the protease domain. We have found that the upstream kringles of apo(a) undergo Molecular Drive-like processes that produce high intraspecies similarity, whereas the downstream kringles evolve in a molecular clock-like manner and show an high interspecies sequence similarity. The latter regions are obviously suitable for dating the duplication event by which Apo(a) arose from plasminogen, but only if they evolve at the same rate in the two genes. Thus, we propose a "Molecular Clock Test" for assessing whether the comparison of two paralogous genes (or gene regions) can give reliable information on the dating of their origin by duplication. Applying this test to the kringle-4 domain of apo(a) and plasminogen gene, we demonstrate that the separation between the two genes by duplication dates back at about 90 Mya immediately before the radiation of mammals.

Characterization of p53 mutations in colorectal liver metastases and correlation with clinical parameters.

The presence and type of mutations of the p53 tumor suppressor gene were determined in 40 patients undergoing curative hepatic resection for metastatic colorectal carcinoma. This represents the largest series in the literature on the screening of p53 mutations for liver metastases. The analysis was performed in exons 5-9 by denaturing gradient gel electrophoresis followed by direct sequencing. Forty-five percent of tumors showed mutation in p53, and this was observed only in exons 5-8. Mutations at codon positions 167, 196, 204, 213, 245, 281, 282, 286, and 306; deletion of codon 251 and of the first nucleotide of codon 252; and Leu residue (CTC) insertion downstream codon 252 are reported for the first time in colorectal liver metastasis. Mutations at codon positions 163, 248, and 273 have been reported previously. Correlation of p53 status with clinical parameters showed that patients with mutated p53 had a statistically higher number of lesions when compared with patients with wild-type p53 (P<0.050). In particular, of patients with mutated p53, 41% had three or more metastases compared with 14% of patients with wild-type p53. Synchronous metastases were present in 70% of the patients with p53 mutations and in only 29% of patients with wild-type p53 (P<0.025). In addition, patients with p53 mutations are more likely to develop recurrence (73%) compared with patients with wild-type p53 (33%; P<0.001). Other factors considered, including preoperative carcinoembryonic antigen level, bilobar distribution, and size of the lesion(s), did not show significant correlation with p53 status. These results suggest that p53 status might be an important prognostic indicator to predict the pattern and likelihood of treatment failure after hepatic resection.

Linguistic approaches to the analysis of sequence information.

Biological macromolecules have many features that resemble modern languages. Thus, linguistic approaches to the analysis of sequence information are becoming powerful tools for deciphering genetic texts. The methodologies used, to date, to determine the global parameters of the genetic language and meaningful patterns within it are described.

Long-branch attraction phenomenon and the impact of among-site rate variation on rodent phylogeny.

The phylogenetic relationships among major lineages of rodents is one of the issues most debated by both paleontologists and molecular biologists. In the present study, we have analyzed all complete mammalian mitochondrial genomes available in the databases, including five rodent species (rat, mouse, dormouse, squirrel and guinea-pig). Phylogenetic analyses were performed on H-strand amino acid sequences by means of maximum-likelihood and on H-strand protein-coding and ribosomal genes by means of distance methods. Also, log-likelihood ratio tests were applied to different tree topologies under the assumption of rodent monophyly, paraphyly or polyphyly. The analyses significantly rejected rodent monophyly and showed the existence of two differentiated clades, one containing non-murids (dormouse, squirrel and guinea-pig) and the other containing murids (rat and mouse). Long-branch attraction between murids and the outgroups could not be responsible for the existence of two different rodent clades, as no significant differences in evolutionary rate have been observed, except in the case of the squirrel, which shows a lower rate. The impact of among-site rate variation models on the phylogeny of rodents has been evaluated using the gamma distribution model. Results have shown that relationships among rodents remained unchanged, and the general topology of the tree was not affected, even though some branches were not properly resolved, most likely due to a lack of fit between estimated and real rate heterogeneity parameters.

Sequence analysis and compositional properties of untranslated regions of human mRNAs.

A detailed computer analysis of the untranslated regions, 5'-UTR and 3'-UTR, of human mRNA sequences is reported. The compositional properties of these regions, compared with those of the corresponding coding regions, indicate that 5'-UTR and 3'-UTR are less affected by the isochore compartmentalization than the corresponding third codon positions of mRNAs. The presence of higher functional constraints in 5'-UTR is also reported. Dinucleotide analysis shows a depletion of CpG and TpA in both sequences. A search for significant sequence motifs using the WORDUP algorithm reveals the patterns already known to have a functional role in the mRNA UTR, and several other motifs whose functional roles remain to be demonstrated. This type of analysis may be particularly useful for guiding site-directed mutagenesis experiments. In addition, it can be used for assessing the nature of anonymous sequences now produced in large amounts in megabase sequencing projects.

Structural and compositional features of untranslated regions of eukaryotic mRNAs.

The important role of 5' and 3' untranslated regions of eukaryotic mRNAs in gene regulation and expression is now widely accepted. In order to study the general structural and compositional features of these sequences we developed UTRdb, a specialized database of 5' and 3'-UTR sequences from seven different taxonomic groups of eukaryotic mRNAs cleaned of redundancy. The analysis of the UTR sequences contained in this database showed that 5'-UTR sequences, on average 200 nucleotides long, are 1.5-3 times shorter than the corresponding 3'-UTR sequences in the various taxonomic groups considered here. As to their compositional properties on average 5'-UTR sequences resulted in all cases GC richer than 3'-UTR sequences, and significant correlations were found between the GC content of 5' and 3'-UTR sequences and the GC content of the third silent codon positions of the corresponding protein coding genes. The dinucleotide analysis showed a differential depletion of CpG in vertebrate 5' and 3'-UTR, with 5'-UTR sequences being more CpG-rich, and a generalized depletion of TpA in both 5' and 3'-UTR was observed in all eukaryotic sequence collections.

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.

A 67-kDa protein binding to the curved DNA in the main regulatory region of the rat mitochondrial genome.

We have purified, by sequence-specific affinity chromatography, a mitochondrial (mt) matrix protein which binds to the curved DNA located between the replication origin (ori) of the leading strand (ori-H) and the two transcription promoters in the rat mt genome. The protein was characterized by gel electrophoresis as a 67-kDa polypeptide and seems to be involved in the DNA contact on the mt light strand. This protein differs (in the size and location of its DNA-binding site) from other DNA-binding proteins studied so far in animal mt systems. We suggest a role for the 67-kDa protein, assisted by other proteins, in regulating the initiation of leading-strand replication.

Mammalian mitochondrial D-loop region structural analysis: identification of new conserved sequences and their functional and evolutionary implications.

This paper reports the first comprehensive analysis of Displacement loop (D-loop) region sequences from ten different mammalian orders. It represents a systematic evolutionary study at the molecular level on regulatory homologous regions in organisms belonging to a well defined class, mammalia, which radiated about 150 million years ago (Mya). We have aligned and analyzed 26 complete D-loop region sequences available in the literature and the fat dormouse sequence, recently determined in our laboratory. The novelty of our alignment consists of the extensive manual revision of the preliminary output obtained by computer program to optimize sequence similarity, particularly for the two peripheral domains displaying heterogeneity in length and the presence of repeated sequences. The multialignment is available at the WWW site: http://www.ba.cnr.it/dloop.html. Our comparative study has allowed us to identify new conserved sequence blocks present in all the species under consideration and events of insertion/deletion which have important implications in both functional and evolutionary aspects. In particular we have detected two blocks, about 60 bp long, extended termination associated sequences (ETAS1 and ETAS2) conserved in all the organisms considered. Evaluation against experimental work suggests a possible functional role of ETAS1 and ETAS2 in the regulation of replication and transcription and targeted experimental approaches. The analyses on conserved sequence blocks (CSBs) clearly indicate that CSB1 is the only very essential element, common to all mammalian mt genomes, while CSB2 and CSB3 could be involved in different though related functions, probably species specific, and thus more linked to nuclear mitochondrial coevolutionary processes. Our hypothesis on the different functional implications of the conserved elements, CSBs and TASs, reported so far as main regulatory signals, would explain the different conservation of these elements in evolution. Moreover the intra-order comparison of the D-loop regions highlights peculiar features useful to define the evolutionary dynamics of this region in closely related species.

Evolutionary genomics in Metazoa: the mitochondrial DNA as a model system.

One of the most important aspects of mitochondrial (mt) genome evolution in Metazoa is constancy of size and gene content of mtDNA, whose plasticity is maintained through a great variety of gene rearrangements probably mediated by tRNA genes. The trend of mtDNA to maintain the same genetic structure within a phylum (e.g., Chordata) is generally accepted, although more recent reports show that a considerable number of transpositions are observed also between closely related organisms. Base composition of mtDNA is extremely variable. Genome GC content is often low and, when it increases, the two complementary bases distribute asymmetrically, creating, particularly in vertebrates, a negative GC-skew. In mammals, we have found coding strand base composition and average degree of gene conservation to be related to the asymmetric replication mechanism of mtDNA. A quantitative measurement of mtDNA evolutionary rate has revealed that each of the various components has a different evolutionary rate. Non-synonymous rates are gene specific and fall in a range comparable to that of nuclear genes, whereas synonymous rates are about 22-fold higher in mt than in nuclear genes. tRNA genes are among the most conserved but, when compared to their nuclear counterparts, they evolve 100 times faster. Finally, we describe some molecular phylogenetic reconstructions which have produced unexpected outcomes, and might change our vision of the classification of living organisms.

Analysis of oligonucleotide AUG start codon context in eukariotic mRNAs.

The AUG start codon context features have been investigated by analyzing eukaryotic mRNAs belonging to various taxonomic groups. The functional relevance of each specific position surrounding the AUG start codon has been established as a function of the measured shift between base composition observed at that particular position, and base composition averaged over all the 5'untranslated regions. A more detailed analysis carried out on human genes belonging to different isochores showed significant isochore-specific fea-tures that cannot be explained only by a mutational bias effect. The most represented heptamers spanning from position -3 to +4 with respect to the initiator AUG have been determined for mRNAs belonging to different taxonomic groups and a web page utility has been set up (http://bigarea.area.ba.cnr.it:8000/BioWWW/ATG.html) to determine the relative abundance of a user submitted oligonucleotide context in a given species or taxon.

Evolution of the mitochondrial genetic system: an overview.

Mitochondria, semi-autonomous organelles possessing their own genetic system, are commonly accepted to descend from free-living eubacteria, namely hydrogen-producing alpha-proteobacteria. The progressive loss of genes from the primitive eubacterium to the nucleus of the eukaryotic cell is strongly justified by the Muller rachet principle, which postulates that asexual genomes, like mitochondrial ones, accumulate deleterious and sublethal mutations faster than sexual genomes, like the nucleus. According to this principle, the mitochondrial genome would be doomed to death; instead, we observe that the mitochondrial genome has a variable size and structure in the different organisms, though it contains more or less the same set of genes. This is an example of genetic conservation versus structural diversity. From an evolutionary point of view the genetic system of organelles is clearly under strong selective pressure and for its survival it needs to utilize strategies to slow down or halt the ratchet. Anyway, the mitochondrial genome changes with time, and the rate of evolution is different for both diverse regions of the mtDNA and between lineages, as demonstrated in the case of mammalian mt genomes. We report here our data on the evolution of the mitochondrial DNA in mammals which demonstrate the suitability of mtDNA as a molecular tool for evolutionary analyses.

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.