Evolutionary patterns of the mitochondrial genome in Metazoa: exploring the role of mutation and selection in mitochondrial protein coding genes

The mitochondrial genome is a fundamental component of the eukaryotic domain of life, encoding for several important subunits of the respiratory chain, the main energy production system in cells. The processes by means of which mitochondrial DNA (mtDNA) replicates, expresses itself and evolves have been explored over the years, although various aspects are still debated. In this review, we present several key points in modern research on the role of evolutionary forces in affecting mitochondrial genomes in Metazoa. In particular, we assemble the main data on their evolution, describing the contributions of mutational pressure, purifying, and adaptive selection, and how they are related. We also provide data on the evolutionary fate of the mitochondrial synonymous variation, related to the nonsynonymous variation, in comparison with the pattern detected in the nucleus.

Elevated mutational pressure characterizes the evolution of the mitochondrial synonymous variation, whereas purging selection, physiologically due to phenomena such as cell atresia and intracellular mtDNA selection, guarantees coding sequence functionality. This enables mitochondrial adaptive mutations to emerge and fix in the population, promoting mitonuclear coevolution.

Mitochondrial DNA variability of West New Guinea populations

This paper reports human mitochondrial DNA variability in West New Guinea (the least known, western side of the island of New Guinea), not yet described from a molecular perspective. The study was carried out on 202 subjects from 12 ethnic groups, belonging to six different Papuan language families, representative of both mountain and coastal plain areas. Mitochondrial DNA hypervariable region 1 (HVS 1) and the presence of the 9-bp deletion (intergenic region COII-tRNA(Lys)) were investigated. HVS 1 sequencing identified 73 polymorphic sites defining 89 haplotypes; the 9-bp deletion, which is considered a marker of Austronesian migration in the Pacific, was found to be absent in the whole West New Guinea study sample. Statistical analysis applied to the resulting haplotypes reveal high heterogeneity and an intersecting distribution of genetic variability in these populations, despite their cultural and geographic diversity. The results of subsequent phylogenetic approaches subdivide mtDNA diversity in West New Guinea into three main clusters (groups I-III), defined by sets of polymorphisms which are also shared by some individuals from Papua New Guinea. Comparisons with worldwide HVS 1 sequences stored in the MitBASE database show the absence of these patterns outside Oceania and a few Indonesian subjects, who also lack the 9-bp deletion. This finding, which is consistent with the effects of genetic drift and prolonged isolation of West New Guinea populations, lead us to regard these patterns as New Guinea population markers, which may harbor the genetic memory of the earliest human migrations to the island.

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.

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.

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.

New p53 mutations in hilar cholangiocarcinoma

BACKGROUND

Inactivation of the tumour suppressor gene, p53, is the commonest genetic abnormality in human cancer. The study of the type of p53 mutation in a given tumour may provide prognostic information, clues to aetiology and become useful for therapeutics.

MATERIALS AND METHODS

The molecular characterisation of p53 was performed by restriction analysis, denaturing gradient gel electrophoresis, and gene sequencing for exons 5-9.

RESULTS

We report, p53 mutational analysis in exons 5-9 in 29 European patients with hilar cholangiocarcinoma who underwent attempted resection. Four patients (14%) showed somatic single nucleotide substitutions with amino acid changes (146, 163, 175, 158, and 175) with one showing mutations in two different positions involving a loss of two CfoI sites. All the mutations occurred in exon 5. Three patients had a germline polymorphism (10%) with a silent substitution in codon 213 (exon 6).

CONCLUSIONS

The systematic screening for p53 mutations in European patients with hilar cholangiocarcinoma has shown that the type of mutation (except 175) is different and its incidence is much lower when compared to the pattern previously reported for intrahepatic cholangiocarcinoma in East Asian patients. A probable explanation is that the presence and type of p53 mutation is dependent on geographic and environmental factors which vary in different populations.

Lineage-specific evolutionary rate in mammalian mtDNA

The existence of a lineage-specific nucleotide substitution rate in mammalian mtDNA has been investigated by analyzing the mtDNA of all available species, that is, 35 complete mitochondrial genomes from 14 mammalian orders. A detailed study of their evolutionary dynamics has been carried out on both ribosomal RNA and first and second codon positions (P12) of H-strand protein-coding genes by using two different types of relative-rate tests. Results are quite congruent between ribosomal and P12 sites. Significant rate variations have been observed among orders and among species of the same order. However, rate variation does not exceed 1.8-fold between the fastest (Proboscidea and Primates) and the slowest (Perissodactyla) evolving orders. Thus, the observed mitochondrial rate variations among taxa do not invalidate the suitability of mtDNA for drawing mammalian phylogeny. Dependence of evolutionary rate differences on variations in mutation and/or fixation rates was examined. Body size, generation time, and metabolic rate were tested, and no significant correlation was observed between them and the taxon-specific evolutionary rates, most likely because the latter might be influenced by multiple overlapping variable constraints.

[Juvenile primary refractory anemia with excess of blasts. Report of a case]

Myelodysplastic syndromes (MDS) are diseases typical of the adult age, characterized by a clonal alteration of the totipotent staminal cell which causes an inefficient hemopoiesis, reduction of bone marrow cell number, increased bone marrow cell destruction, dysplasia of at least two of the three hemopoietic cell lines and by the tendency to evolve towards acute myeloid leukemia (AML). In patients with MSD, particularly in the advanced primary form, single or multiple chromosomal abnormalities can be found frequently, which may show up and/or modify themselves in any moment of the disease (multistep pathogenetic hypothesis) and whose severity influences significantly the prognosis of MSD patients. In November 1998, a 22 year old female patient (C. C.) was admitted to the Department of Internal medicine of the Garibaldi Hospital, Catania for anemia of unknown origin. The normalities of the hematochemical tests and of other instrumental examinations, induces to perform an osteo-medullary biopsy which revealed the presence of a typical MSD, refractory anemia with excess of blasts type. This disease in its primary form is rare in youth and has a high tendency to evolve in AML. The good clinical conditions of the patient, the absence of chromosomal abnormalities, the normal levels of HbF, the short time interval of the initial diagnosis induced to proceed to bone marrow transplantation, that, as shown by the data reported in the international literature, may give her a good quoad vitam prognosis.

[Histiocytic necrotizing lymphadenitis or Kikuchi-Fujimoto disease. Report of a case]

Histiocytic Necrotising Lymphadenitis by Kikuchi-Fujimoto (NHL) is a rare disease of unknown etiopathogenesis, characterized by cervical lymphadenomegaly, fever and asthenia. It has a good prognosis with a complete functional recovery of the affected lymph nodes. In 1998 a 28 year-old patient (A. G.) was admitted in the Department of Internal Medicine, Garibaldi Hospital, Catania for fever, asthenia and cervical lymphadenopathy. Hemato-chemical tests performed during hospitalization showed only relative leukopenia and a significant increase of ESR. An initial cervical lymph node biopsy made the diagnosis of Hodgkin's lymphoma. However, the favourable clinical course and the normalization of the previous altered hematochemical tests, induced to perform a second lymph node biopsy which led to diagnosis of NHL. The patient was given antibiotics and anti-inflammatory drugs for ten days with a complete healing which persisted at a twelve-month follow-up. The conclusion is drawn that NHL of unknown etiopathogenesis and with a good prognosis is more frequent than what revealed by the data in the literature because of its insidious and aspecific clinical presentation.

UTRdb and UTRsite: specialized databases of sequences and functional elements 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, 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. All internet resources implemented for retrieval and functional analysis of 5' and 3' untranslated regions of eukaryotic mRNAs are accessible at http://bigarea.area.ba.cnr.it:8000/EmbIT/UTRH ome/

MitoNuc and MitoAln: two related databases of nuclear genes coding for mitochondrial proteins

Mitochondria, besides their central role in energy metabolism, have recently been found to be involved in a number of basic processes of cell life and to contribute to the pathogenesis of many degenerative diseases. All functions of mitochondria depend on the interaction of nuclear and organellar genomes. Mitochondrial genomes have been extensively sequenced and analysed and the data collected in several specialised databases. In order to collect information on nuclear coded mitochondrial proteins we developed MitoNuc and MitoAln, two related databases containing, respectively, detailed information on sequenced nuclear genes coding for mitochondrial proteins in Metazoa and yeast, and the multiple alignments of the relevant homologous protein coding regions. MitoNuc and MitoAln retrieval through SRS at http://bio-www.ba.cnr.it:8000/srs6/ can easily allow the extraction of sequence data, subsequences defined by specific features and nucleotide or amino acid multiple alignments.

MitBASE : a comprehensive and integrated mitochondrial DNA database. The present status

MitBASE is an integrated and comprehensive database of mitochondrial DNA data which collects, under a single interface, databases for Plant, Vertebrate, Invertebrate, Human, Protist and Fungal mtDNA and a Pilot database on nuclear genes involved in mitochondrial biogenesis in Saccharomyces cerevisiae. MitBASE reports all available information from different organisms and from intraspecies variants and mutants. Data have been drawn from the primary databases and from the literature; value adding information has been structured, e.g., editing information on protist mtDNA genomes, pathological information for human mtDNA variants, etc. The different databases, some of which are structured using commercial packages (Microsoft Access, File Maker Pro) while others use a flat-file format, have been integrated under ORACLE. Ad hoc retrieval systems have been devised for some of the above listed databases keeping into account their peculiarities. The database is resident at the EBI and is available at the following site: http://www3.ebi.ac.uk/Research/Mitbase/mitbas e.pl. The impact of this project is intended for both basic and applied research. The study of mitochondrial genetic diseases and mitochondrial DNA intraspecies diversity are key topics in several biotechnological fields. The database has been funded within the EU Biotechnology programme.

Isochore specificity of AUG initiator context of human genes

The efficiency of AUG start codon recognition in translation initiation is modulated by its sequence context. Here we investigated a non-redundant set of 5914 human genes and show that this context is different in genes located in different isochores. In particular, of the two main consensus start sequences, RCCaugR is five-fold more represented than AARaugR in genes from the GC-rich H3 isochores compared to genes from the GC-poor L isochores. Furthermore, genes located in GC-rich isochores have shorter 5' UTRs and stronger avoidance of upstream AUG than genes located in GC-poor isochores. This suggests that genes requiring highly efficient translation are located in GC-rich isochores and genes requiring fine modulation of expression are located in GC-poor isochores. This is in agreement with independent data from the literature concerning the location of housekeeping and tissue-specific genes, respectively.

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.

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.

Guinea pig p53 mRNA: identification of new elements in coding and untranslated regions and their functional and evolutionary implications

We report the sequence of the guinea pig p53 cDNA. The comparative analysis of the coding and noncoding regions of p53 cDNAs of all available complete vertebrate sequences has allowed us to single out new conserved signals possibly involved in p53 functional activity. We have focused our attention on the most variable region of the protein, the proline (P)-rich domain, suggested to play a fundamental role in antiproliferative pathways. In this domain we have identified the PXXXXP repeated motif and singled out a common consensus sequence that can be considered a signature for mammalian p53: PXXXXPX{0,4}PX{0,9}PA(T,P,I,)(S,P)WPL. We have demonstrated the significance of the PXXXXP motif in SH3-binding protein and suggested its structure to be a loop. Also, the 5' and 3' untranslated regions (UTRs) of the guinea pig were sequenced, and this study represents the first detailed structural analysis of the UTRs of the p53 mRNAs available in literature. The 5' UTR of guinea pig (233 nt) can be folded into a stable secondary structure resembling that predicted in mouse. The 3' UTR of guinea pig is 771 nt long and shows higher similarity with human than with rodent sequences, having a region of about 350 nt that is deleted in rat and mouse. In the 3' UTR we have identified the presence of a mammalian-wide interspersed repeat sequence and of a cytoplasmic polyadenylation element, which could be involved in translational activation by promoting polyadenylation of mRNA, providing information about a possible mechanism of regulation of p53 expression mediated by the 3' UTR of the mRNA. The observations presented here could open new avenues to targeted mutations and experimental approaches useful in investigating new regulation mechanisms of p53 translation, activity, and stability.

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.

[Ultra-low laparoscopic rectal resection and colo-anal anastomosis. Experimental study on swine]

Nowadays oncologic surgery is defining new criteria in the treatment of rectal cancer: preservation of sfincterial function, ultra-low resections with distal margin at only 2 cm distally to the tumor, role of mesorectum as preferential site of lymphatic diffusion, preservation of lombo-aortic and pelvic nerves. Laparoscopy is showing good results in bowel surgery so as previously got on biliary tract: less visceral manipulation, less stimulation of immunologic system, less pain, early resumption of peristalsis and food intake, better recovery, less hospital stay. We experimented on pig ultra-low laparoscopic resection of the rectum, with preservation of sfincterial function, and mechanical anastomosis at the upper edge of the sfincterial ring. The evaluation of surgical technique, post-operative supervision, and follow-up instrumental control (anal manometry, endo-rectal ultrasonography, sfincterial electtromanometry, Rx barium enema) show technical feasibility and confirm a better recovery with regular defecatory function.

Complete sequence of the mitochondrial DNA in the sea urchin Arbacia lixula: conserved features of the echinoid mitochondrial genome

The complete nucleotide sequence (15,719 nucleotides) of the mitochondrial DNA (mtDNA) from the sea urchin Arbacia lixula is presented. The comparison of gene arrangement between different echinoderm orders of the same class provides evidence that the gene organization is conserved within the same echinoderm class. The peculiarities of sea urchin mtDNA features, already described, are confirmed by the A. lixula mtDNA sequence. The comparison of the entire sequences of mtDNA among A. lixula, Paracentrotus lividus, and Strongylocentrotus purpuratus allowed us to detect peculiar features, common to the three sea urchin species, that can represent the molecular signature of the mt genome in the sea urchin group. Analysis of the nucleotide composition indicates that A. lixula mtDNA, in contrast with the mtDNA of other sea urchins, shows a bias in the use of T and tends to avoid the use of C, most evident in the neutral part of the molecule, such as the third codon positions. This observation indicates that the three sea urchin mtDNAs evolve under different mutation pressure. Analysis of the sequence evolution allowed us to confirm the phylogenetic tree. However, the absolute divergence time, calculated on the basis of paleontological estimates, largely diverged from the expected one.

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.

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.

Molecular classification of living organisms

Recent studies in molecular evolution have generated strong conflicts in opinion as to how world living organisms should be classified. The traditional classification of life into five kingdom has been challenged by the molecular analysis carried out mostly on rRNA sequences, which supported the division of the extant living organisms into three major groups: Archaebacteria, Eubacteria, and Eukaryota. As to the problem of placing the root of the tree of life, the analysis carried out on a few genes has provided discrepant results. In order to measure the genetic distances between species, we have carried out an evolutionary analysis of the glutamine synthetase genes, which previously have been revealed to be good molecular clocks, and of the small and large rRNA genes. All data demonstrate that archaebacteria are more closely related to eubacteria than to eukaryota, thus supporting the classical division of living organisms into two main superkingdoms, Prokaryota and Eukaryota.

Glutamine synthetase gene evolution in bacteria

The evolution of the prokaryotic glutamine synthase (GS) genes, namely the GSI and GSII isoforms, has been investigated using the second codon positions, which have previously proven to behave as a good molecular clock. Our data confirm the early divergence between prokaryotic and eukaryotic GSII before the splitting between plants and animals. The phylogenetic tree of the GSI isoforms shows Archaebacteria to be more closely related to Eubacteria than to Eukaryotes. This finding is confirmed by the phylogenetic analysis carried out on both large and small subunits of rRNA. However, differently from the rRNA analyses, Crenarchaeota and Euryarchaeota Archaebacteria, as well as high- and low-GC gram-positive bacteria, appear to be polyphyletic. We provide evidence that the observed polyphyly of Archaebacteria might be only apparent, resulting from a gene duplication event preceding the split between Archaebacteria and Eubacteria and followed by the retention of only one isoform in the extant lineages. Both gram-negative bacteria and high-GC gram-positive bacteria, which appear closely related, have GS activity regulated by an adenylylation/deadenylylation mechanism. A lateral gene transfer from Archaebacteria to low-GC eubacteria is invoked to explain the observed polyphyly of gram-positive bacteria.