Cecilia Saccone

A new method for calculating evolutionary substitution rates

SummaryIn this paper we present a new method for analysing molecular evolution in homologous genes based on a general stationary Markov process. The elaborate statistical analysis necessary to apply the method effectively has been performed using Monte Carlo technqiues. We have applied our method to the silent third position of the codon of the five mitochondrial genes coding for identified proteins of four mammalian species (rat, mouse, cow and man). We found that the method applies satisfactorily to the three former species, while the last appears to be outside the scope of the present approach. The method allows one to calculate the evolutionarily effective silent substitution rate (vs) for mitochondrial genes, which in the species mentioned above is 1.4×10−8 nucleotide substitutions per site per year. We have also determined the divergence time ratios between the couples mousecow/rat-mouse and rat-cow/rat-mouse. In both cases this value is approximately 1.4.

The complete nucleotide sequence of theRattus norvegicus mitochondrial genome: Cryptic signals revealed by comparative analysis between vertebrates

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

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.

The Main Regulatory Region of Mammalian Mitochondrial DNA: Structure-Function Model and Evolutionary Pattern

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

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.

Evolutionary analysis of cytochrome b sequences in some perciformes: Evidence for a slower rate of evolution than in mammals

To obtain information relative to the phylogenesis and microevolutionary rate of fish mitochondrial DNA, the nucleotide sequence of cytochrome b gene in seven fish species belonging to the order of Perciformes was determined. Sequence analysis showed that fish mitochondrial DNA has a nucleotide compositional bias similar to that of sharks but lower compared to mammals and birds. Quantitative evolutionary analysis, carried out by using a markovian stochastic model, clarifies some phylogenetic relationships within the Perciformes order, particularly in the Scombridae family, and between Perciformes, Gadiformes, Cypriniformes, and Acipenseriformes. The molecular clock of mitochondrial DNA was calibrated with the nucleotide substitution rate of cytochrome b gene in five shark species having divergence times inferred from paleontological estimates. The results of such analysis showed that Acipenseriformes diverged from Perciformes by about 200 MY, that the Perciformes common ancestor dates back to 150 MY, and that fish mitochondrial DNA has a nucleotide substitution rate three to five times lower than that of mammals.

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

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

Influence of base composition on quantitative estimates of gene evolution.

Publisher Summary This chapter presents the theoretical basis of a stationary Markov model, its mathematical formulation, and a few experimental applications to examine the influence of base composition on quantitative estimates of gene evolution. The basic hypothesis is that in a nucleotide sequence at a given time T, any particular nucleotide is the product of a stochastic process whose properties depend on the particular position in the sequence occupied by the nucleotide under consideration. The probability of mutation and the subsequent mutation fixation of that nucleotide do not depend on the history of the nucleotide but only on the environmental conditions and on the biological constraints prevailing at the time. In protein-coding genes, for example, all the nucleotide sites at the first codon position are grouped together and their collective pattern and tempo of substitution is analyzed. Similarly, the collective behavior of nucleotide sites at the second and third, silent codon positions defines other evolutionary dynamics, all distinct from each other.

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.

Mitochondrial DNA in metazoa: degree of freedom in a frozen event.

The mitochondrial genome (mtDNA), due to its peculiar features such as exclusive presence of orthologous genes, uniparental inheritance, lack of recombination, small size and constant gene content, certainly represents a major model system in studies on evolutionary genomics in metazoan. In 800 million years of evolution the gene content of metazoan mitochondrial genomes has remained practically frozen but several evolutionary processes have taken place. These processes, reviewed here, include rearrangements of gene order, changes in base composition and arising of compositional asymmetry between the two strands, variations in the genetic code and evolution of codon usage, lineage-specific nucleotide substitution rates and evolutionary patterns of mtDNA control regions.