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 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 non-synonymous variation, in comparison with the pattern detected in the nucleus. Elevated mutational pressure characterises 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 mito-nuclear 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.

[Histiocytic necrotizing lymphadenitis or Kikuchi-Fujimoto disease. Report of a case]. / Linfoadenite istiocitica necrotizzante o malattia di Kikuchi-Fujimoto. Descrizione di un caso.

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.

[Juvenile primary refractory anemia with excess of blasts. Report of a case]. / Anemia refrattaria con eccesso di blasti primaria giovanile. Descrizione di un caso.

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.

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.

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.

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.

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.

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.

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/

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.

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.

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.

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.