Nucleotide substitution rate of mammalian mitochondrial genomes

We present here for the first time a comprehensive study based on the analysis of closely related organisms to provide an accurate determination of the nucleotide substitution rate in mammalian mitochondrial genomes. This study examines the evolutionary pattern of the different functional mtDNA regions as accurately as possible on the grounds of available data, revealing some important "genomic laws." The main conclusions can be summarized as follows. (1) High intragenomic variability in the evolutionary dynamic of mtDNA was found. The substitution rate is strongly dependent on the region considered, and slow- and fast-evolving regions can be identified. Nonsynonymous sites, the D-loop central domain, and tRNA and rRNA genes evolve much more slowly than synonymous sites and the two peripheral D-loop region domains. The synonymous rate is fairly uniform over the genome, whereas the rate of nonsynonymous sites depends on functional constraints and therefore differs considerably between genes. (2) The commonly accepted statement that mtDNA evolves more rapidly than nuclear DNA is valid only for some regions, thus it should be referred to specific mitochondrial components. In particular, nonsynonymous sites show comparable rates in mitochondrial and nuclear genes; synonymous sites and small rRNA evolve about 20 times more rapidly and tRNAs about 100 times more rapidly in mitochondria than in their nuclear counterpart. (3) A species-specific evolution is particularly evident in the D-loop region. As the divergence times of the organism pairs under consideration are known with sufficient accuracy, absolute nucleotide substitution rates are also provided.

Update of the Human MitBASE database

Human MitBASE is a database collecting human mtDNA variants. This database is part of a greater mitochondrial genome database (MitBASE) funded within the EU Biotech Program. The present paper reports the recent improvements in data structure, data quality and data quantity. As far as the database structure is concerned it is now fully designed and implemented. Based on the previously described structure some changes have been made to optimise both data input and data quality. Cross-references with other bio-databases (EMBL, OMIM, MEDLINE) have been implemented. Human MitBASE data can be queried with the MitBASE Simple Query System (http://www.ebi.ac.uk/htbin/Mitbase/mit base.pl) and with SRS at the EBI under the 'Mutation' section (http://srs.ebi.ac.uk/srs5/). At present the HumanMitBASE node contains approximately 5000 variants related to studies investigating population polymorphisms and pathologies.

Vertebrate MitBASE: a specialised database on vertebrate mitochondrial DNA sequences

Vertebrate MitBASE is a specialized database where all the vertebrate mitochondrial DNA entries from primary databases are collected, revised and integrated with new information emerging from the literature. Variant sequences are also analyzed, aligned and linked to reference sequences. Data related to the same species and fragment can be viewed over the WWW. The database has a flexible interface and a retrieval system to help non-expert users and contains information not currently available in the primary databases. Vertebrate MitBASE is now available through the MitBASE home page at URL: http://www.ebi.ac.uk/htbin/Mitbase/mitb ase.pl. This work is part of a larger project, MitBASE which is a network of databases covering the full panorama of knowledge on mitochondrial DNA from protists to human sequences.

UTRdb: a specialized database of 5' and 3' untranslated regions of eukaryotic mRNAs

The 5' and 3' untranslated regions of eukaryotic mRNAs may play a crucial role in the regulation of gene expression controlling mRNA localization, stability and translational efficiency. For this reason we developed UTRdb (http://bigarea.area.ba.cnr.it:8000/BioWWW/#U TRdb), a specialized database of 5' and 3' untranslated sequences of eukaryotic mRNAs cleaned from redundancy. UTRdb entries are enriched with specialized information not present in the primary databases including the presence of nucleotide sequence patterns already demonstrated by experimental analysis to have some functional role. All these patterns have been collected in the UTRsite database so that it is possible to search any input sequence for the presence of annotated functional motifs. Furthermore, UTRdb entries have been annotated for the presence of repetitive elements.

MitBASE: a comprehensive and integrated mitochondrial DNA database

MitBASE is an integrated and comprehensive database of mitochondrial DNA data which collects all available information from different organisms and from intraspecie variants and mutants. Research institutions from different countries are involved, each in charge of developing, collecting and annotating data for the organisms they are specialised in. The design of the actual structure of the database and its implementation in a user-friendly format are the care of the European Bioinformatics Institute. The database can be accessed on the Web at the following address: http://www.ebi.ac. uk/htbin/Mitbase/mitbase.pl. The impact of this project is intended for both basic and applied research. The study of mitochondrial genetic diseases and mitochondrial DNA intraspecie diversity are key topics in several biotechnological fields. The database has been funded within the EU Biotechnology programme.

Update of AMmtDB: a database of multi-aligned metazoa mitochondrial DNA sequences

The present paper describes AMmtDB, a database collecting the multi-aligned sequences of vertebrate mitochondrial genes coding for proteins and tRNAs, as well as the multiple alignment of the mammalian mtDNA main regulatory region (D-loop) sequences. The genes coding for proteins are multi-aligned based on the translated sequences and both the nucleotide and amino acid multi-alignments are provided. As far as the genes coding for tRNAs are concerned, the multi-alignments based on the primary and the secondary structures are both provided; for the mammalian D-loop multi-alignments we report the conserved regions of the entire D-loop (CSB1, CSB2, CSB3, the central region, ETAS1 and ETAS2) as defined by Sbisà et al. [ Gene (1997), 205, 125-140). A flatfile format for AMmtDB has been designed allowing its implementation in SRS (http://bio-www.ba.cnr.it:8000/BioWWW/#AMMTDB ). Data selected through SRS can be managed using GeneDoc or other programs for the management of multi-aligned data depending on the user's operative system. The multiple alignments have been produced with CLUSTALV and PILEUP programs and then carefully optimized manually.

Asymmetrical directional mutation pressure in the mitochondrial genome of mammals

The base composition of 25 complete mammalian mitochondrial (mt) genomes has been analyzed taking into account all three codon positions (P1230 and fourfold degenerate sites (P4FD) of H-strand genes. In the nontranscribed L strand, G is the less represented base and A is the most represented one in all cases, while C and T differ among species. H-strand protein-coding genes show an asymmetric distribution of the four bases between the two strands. The asymmetry indexes AT and GC skews on P4FD are much higher than those on P123, suggesting the existence of asymmetrical directional mutation pressure. Relationships between the compositional features and transcription of replication processes have been investigated in order to find a possible mechanism that could explain the origin of this asymmetry. AT and GC skews, the base composition in fourfold degenerate sites, and the number of variable sites for each gene are significantly correlated with the duration of single-stranded state of the H-stranded genes during replication. We tested different replication-related hypotheses, such as the existence of biased dNTP pools, gamma DNA polymerase mispairing, and the asymmetric replication itself. Most of them failed to explain the observed results, hydrolytic deaminations being the only one in agreement with our data. Thus, we hypothesize that one of the crucial processes for the origin of asymmetric and biased base composition of mammalian mitochondrial genomes is the spontaneous deamination of C and A in the H strand during replication.

Complete mitochondrial DNA sequence of the fat dormouse, Glis glis: further evidence of rodent paraphyly

The complete mitochondrial genome of the fat dormouse, Glis glis, has been sequenced (16,602 bp). A total of 23 complete mitochondrial mammalian genomes have been taken into account for phylogenetic reconstruction. Phylogenetic analyses were performed with parsimony, distance (stationary Markov model), and maximum-likelihood methods. In all cases, data strongly support the paraphyly of rodents, with dormouse and guinea pig in a different clade from rat and mouse, reaching bootstrap values of 95%. Rodent monophyly and the existence of Glires (Rodentia and Lagomorpha) are weakly supported, with maximum bootstrap values of 11% and 8.6%, respectively. This result agrees with the analyses of isochore patterns in the nuclear genome and the B2 and B2-like retroposons, which show a close relationship between dormice and guinea pigs rather than between dormice and rats and mice.

p53 mutation is a poor prognostic indicator for survival in patients with hepatocellular carcinoma undergoing surgical tumour ablation

Forty-two patients with hepatocellular carcinoma (HCC) were resected and their tumours were analysed for p53 mutations by GC-clamped denaturing gradient gel electrophoresis (DGGE), single-strand conformation polymorphism (SSCP) and gene sequencing. All the exons have been analysed in this study. Eight of 12 HCCs with cirrhosis due to viral hepatitis and the two patients with sarcomatoid changes displayed p53 mutations. In contrast, no mutation was observed in the fibrolamellar variant (n = 9), non-cirrhotics (n = 13) and alcoholic cirrhosis (n = 6). The mutations observed were in exons 5-8. Two mutations were observed in codons 136 and 213 as well as a T insertion between residues 156 and 157 (exon 5) and these are reported for the first time in HCC. Likewise, the silent mutation polymorphism in codon 213 was noticed in 3 of the 42 patients. Survival analysis of these patients after surgery showed the mean and median survival in patients with wild-type p53 to be 60 and 43 months respectively. In the group with p53 mutations, the mean and median survival was 15 and 12 months. The difference was statistically significant (P= 0.003).

UTRdb: a specialized database of 5'- and 3'-untranslated regions of eukaryotic mRNAs

The important role the untranslated regions of eukaryotic mRNAs may play in gene regulation and expression is now widely acknowledged. For this reason we developed UTRdb, a specialized database of 5'- and 3'-untranslated sequences of eukaryotic mRNAs cleaned from redundancy. UTRdb entries are enriched with specialized information not present in the primary databases, including the presence of functional patterns already demonstrated by experimental analysis to have some functional role. A collection of such patterns is being collected in UTRsite database (http://bio-www.ba.cnr.it:8000/srs5/) which can also be used with appropriate computational tools to detect known functional patterns contained in mRNA untranslated regions.

The mitBASE human dataset structure

MitBASE is a comprehensive and integrated mitochondrial genome database funded within the EU BIOTECH PROGRAM. It is a project for the development and implementation of an integrated and comprehensive database of mitochondrial data which will collect all available information from different organisms and from intraspecies variants and mutants. The present paper describes the structure of the Human dataset in mitBASE where human molecular data are distinguished from clinical and pathological data. MitBASE home page address is: http://www.ebi.ac.uk/htbin/Mitbase/mitb ase.pl

Update of MmtDB: a Metazoa mitochondrial DNA variants database

The present paper describes the improvements in MmtDB, a specialised database designed to collect Metazoa mitochondrial DNA variants. Priority in the data collection has been given to Metazoa for which a large amount of variants is available, e.g., for humans. Starting from the sequences available in the Nucleotide Sequence Databases, the redundant sequences have been removed and new sequences from other sources have been added. Value-added information is associated to each variant sequence, e.g., analysed region, experimental method, tissue and cell lines, population data, sex, age, family code and information about the variation events (nucleotide position, involved gene, restriction site gain or loss). Cross-references are introduced to the EMBL Data Library, as well as an internal cross-referencing among MmtDB entries according to tissual, heteroplasmic, familiar and aplotypical correlation. Furthermore MmtDB has a new section, AMmtDB: Aligned Metazoan mitochondrial biosequences. MmtDB can be accessed through the World Wide Web at URL http://WWW.ba.cnr.it/[symbol: see text]areamt08/MmtDBWWW.htm

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.

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.

Lineage-specific evolution of echinoderm mitochondrial ATP synthase subunit 8

Peculiar evolutionary properties of the subunit 8 of mitochondrial ATP synthase (ATPase8) are revealed by comparative analyses carried out between both closely and distantly related species of echinoderms. The analysis of nucleotide substitution in the three echinoids demonstrated a relaxation of amino acid functional constraints. The deduced protein sequences display a well conserved domain at the N-terminus, while the central part is very variable. At the C-terminus, the broad distribution of positively charged amino acids, which is typical of other organisms, is not conserved in the two different echinoderm classes of the sea urchins and of the sea stars. Instead, a motif of three amino acids, so far not described elsewhere, is conserved in sea urchins and is found to be very similar to the motif present in the sea stars. Our results indicate that the N-terminal region seems to follow the same evolutionary pattern in different organisms, while the maintenance of the C-terminal part in a phylum-specific manner may reflect the co-evolution of mitochondrial and nuclear genes.

Molecular evolution of mammalian aquaporin-2: further evidence that elephant shrew and aardvark join the paenungulate clade

A 328-bp sequence from exon 1 of the gene for aquaporin-2 (AQP2) was compared in 12 mammalian species, representing as many eutherian orders. This sequence encodes the N-terminal half of this kidney-specific water channel protein. Most amino acid replacements, as well as an insertion, have occurred in extracellular loops connecting the transmembrane helices, in agreement with a lower functional importance of these loops. Phylogenetic analyses were performed with parsimony, distance, and maximum-likelihood methods. The AQP2 data set, alone as well as in combination with previously published alpha A-crystallin protein sequences, strongly supports a clade consisting of elephant, hyrax, aardvark, and elephant shrew, reaching bootstrap values of 99%. This finding fully agrees with the only other presently available sequence data sets that include these taxa, those of von Willebrand factor and interphotoreceptor retinoid-binding protein, and suggests that this extended paenungulate clade is one of the most conspicuous superordinal groupings in eutherian phylogeny. Some support was obtained for an artiodactyl/perissodactyl clade, while the grouping of pholidotes with edentates was contradicted.

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.

Apolipoproteins(a): a puzzling evolutionary story

Human apolipoprotein(a), a risk factor for heart disease, has over 80% sequence identity to plasminogen. Plasminogen contains five distinct kringle domains plus a catalytic protease subunit. Human apo(a) consists of multiple copies (the number varies in individuals) of a domain resembling kringle 4, a single copy of a domain resembling kringle 5, and a protease-like domain. The recently cloned hedgehog version of apolipoprotein(a), which contains 31 nearly identical copies of plasminogen kringle 3 and lacks a protease domain, has prompted us to investigate the evolutionary history of the apolipoprotein(a) gene in mammals. Our analysis supports the nonfunctionality of the human apolipoprotein(a) protease domain, and a single (or multiple) duplication of plasminogen gene before mammal radiation, which originated apolipoprotein(a) in mammals.

MmtDB: a Metazoa mitochondrial DNA variants database

The present paper describes the structure of MmtDB-a specialized database designed to collect Metazoa mitochondrial DNA variants. Priority in the data collection is given to the Metazoa species for which a large amount of variants is available, as it is the case for human variants. Starting from the sequences available in the Nucleotide Sequence Databases, the redundant sequences are removed and new sequences from other sources are added. Value-added information are associated to each variant sequence, e.g. analysed region, experimental method, tissue and cell lines, population data, sex, age, family code and information about the variation events (nucleotide position, involved gene, restriction site's gain or loss). Cross-references are introduced to the EMBL Data Library, as well as an internal cross-referencing among MmtDB entries according to their tissual, heteroplasmic, familiar and aplotypical correlation. MmtDB can be accessed through the World Wide Web at URL [see text].

Evolution of the nad3-rps12 gene cluster in angiosperm mitochondria: comparison of edited and unedited sequences

We have analyzed the nad3-rps12 locus for eight angiosperms in order to compare the utility of mitochondrial DNA and edited mRNA sequences in phylogenetic reconstruction. The two coding regions, containing from 25 to 35 editing sites in the various plants, have been concatenated in order to increase the significance of the analysis. Differing from the corresponding chloroplast sequences, unedited mitochondrial DNA sequences seem to evolve under a quasi-neutral substitution process which undifferentiates the nucleotide substitution rates for the three codon positions. By using complete gene sequences (all codon positions) we found that genomic sequences provide a classical angiosperm phylogenetic tree with a clear-cut grouping of monocotyledons and dicotyledons with Magnoliidae at the basal branch of the tree. Conversely, owing to their low nucleotide substitution rates, edited mRNA sequences were found not to be suitable for studying phylogenetic relationships among angiosperms.

Evolutionary analysis of sea urchin mitochondrial tRNAs: folding of the molecules as suggested by the non-random occurrence of nucleotides

Comparative analyses of the mitochondrial tRNA sequences of the sea urchins Arbacia lixula, Paracentrotus lividus and Strongylocentrotus purpuratus revealed that conserved nucleotides may be involved in determining the typical L-shaped spatial conformation of tRNAs. These results shed light on the specific tertiary interactions that allow the folding of the atypical mitochondrial tRNAs into a functional form. A consensus mitochondrial tRNA secondary structure was derived. It shows the presence of nucleotides virtually conserved only in these organisms that represent a sort of molecular signature in sea urchins and suggests a possible physiological role. Finally, we speculate that the non-canonical structure of animal tRNAs, as well as the deviations from the universality of the genetic code, may be due to the reduction in size of the metazoan mitochondrial genome, with the concomitant acquisition of new functions by the mitochondrial tRNAs.

The evolution of the RNase P- and RNase MRP-associated RNAs: phylogenetic analysis and nucleotide substitution rate

We report a detailed evolutionary study of the RNase P- and RNase MRP- associated RNAs. The analyses were performed on all the available complete sequences of RNase MRP (vertebrates, yeast, plant), nuclear RNase P (vertebrates, yeast), and mitochondrial RNase P (yeast) RNAs. For the first time the phylogenetic distance between these sequences and the nucleotide substitution rates have been quantitatively measured.The analyses were performed by considering the optimal multiple alignments obtained mostly by maximizing similarity between primary sequences. RNase P RNA and MRP RNA display evolutionary dynamics following the molecular clock. Both have similar rates and evolve about one order of magnitude faster than the corresponding small rRNA sequences which have been, so far, the most common gene markers used for phylogeny. However, small rRNAs evolve too slowly to solve close phylogenetic relationships such as those between mammals. The quicker rate of RNase P and MRP RNA allowed us to assess phylogenetic relationships between mammals and other vertebrate species and yeast strains. The phylogenetic data obtained with yeasts perfectly agree with those obtained by functional assays, thus demonstrating the potential offered by this approach for laboratory experiments.

The guinea-pig is not a rodent

In 1991 Graur et al. raised the question of whether the guinea-pig, Cavia porcellus, is a rodent. They suggested that the guinea-pig and myomorph rodents diverged before the separation between myomorph rodents and a lineage leading to primates and artiodactyls. Several findings have since been reported, both for and against this phylogeny, thereby highlighting the issue of the validity of molecular analysis in mammalian phylogeny. Here we present findings based on the sequence of the complete mitochondrial genome of the guinea-pig, which strongly contradict rodent monophyly. The conclusions are based on cumulative evidence provided by orthologically inherited genes and the use of three different analytical methods, none of which joins the guinea-pig with myomorph rodents. In addition to the phylogenetic conclusions, we also draw attention to several factors that are important for the validity of phylogenetic analysis based on molecular data.

Isolation of a 25-kDa protein binding to a curved DNA upstream the origin of the L strand replication in the rat mitochondrial genome

The presence of a curved DNA sequence in the gene for the NADH-dehydrogenase subunit 2 of rat mitochondrial genome, upstream from the origin of the light strand replication have been demonstrated through theoretical analysis and experimental approaches. Gel retardation assays showed that this structure makes a complex with a protein component extracted from the mitochondrial matrix. The isolation and purification of this protein is reported. With a Sepharose CL-6B and magnetic DNA affinity chromatography a polypeptide was purified to homogeneity having 25-kDa mass as shown by gel electrophoresis. To functionally characterize this protein, its capability to bind to other sequences of the homologous or heterologous DNA and to specific riboprobes was also investigated. A role for this protein as a trans-acting agent required for the expression of the mammalian mitochondrial genome is suggested.