On the "unreasonable" effects of ELF magnetic fields upon a system of ions.

A recent experiment on a physical, nonbiological system of ions at room temperature has proved that microscopic ion currents can be induced by applying simultaneously two parallel magnetic fields, one rather weak static field, (-->)B(0) and one much weaker alternating field, (-->) B(ac),[B(ac) approximately 10(-3) B(0)] whose frequency coincides with the cyclotron frequency v = qB(0)/2pim of the selected ion. As a result, ionic bursts lasting up to 20 s and with amplitude up to 10 nA arise. The much larger exchanges of energy induced by thermal agitation (the "kT-problem") appear to play no role whatsoever. We have analyzed this problem in the framework of coherent quantum electrodynamics, reaching the following conclusions: (a) as has been shown in previous articles, water molecules in the liquid and solute ions are involved in their ground state in coherent ordered configurations; (b) ions are able to move without collisions among themselves in the interstices between water coherence domains; (c) because of coherence, ions can follow classical orbits in the magnetic fields. A full quantitative understanding of the experiments is thus reached.

The evolution of the mitochondrial D-loop region and the origin of modern man.

The origin of modern man is a highly debated issue that has recently been tackled by using mitochondrial DNA sequences. The limited genetic variability of human mtDNA has been explained in terms of a recent common genetic ancestry, thus implying that all modern-population mtDNAs originated from a single woman who lived in Africa less than 0.2 Mya. This divergence time is based on both the estimation of the rate of mtDNA change and its calibration date. Because different estimates of the rate of mtDNA evolution can completely change the scenario of the origin of modern man, we have reanalyzed the available mitochondrial sequence data by using an improved version of the statistical model, the "Markov clock," devised in our laboratory. Our analysis supports the African origin of modern man, but we found that the ancestral female from which all extant human mtDNAs originated lived in a time span of 0.3-0.8 Mya. Pushing back the date of the deepest root of the human implies that the earliest divergence would have been in the Homo erectus population.

A statistical method for detecting regions with different evolutionary dynamics in multialigned sequences.

We describe a stochastic method for tracing the evolutionary pattern of multialigned sequences. This method allows us to detect gene regions with distinct evolutionary dynamics, e.g., regions that significantly deviate from the expected behavior. Accurate detection of hypervariable or hyperconstrained regions may provide useful information on the structure/function relationship of biosequences. This information can help localize functional constraints. In addition, the selection of distinct evolutionary dynamics may assist in the correct use of biosequences as reliable molecular clocks.

Glutamine synthetase gene evolution: a good molecular clock.

Glutamine synthetase (EC 6.3.1.2) gene evolution in various animals, plants, and bacteria was evaluated by a general stationary Markov model. The evolutionary process proved to be unexpectedly regular even for a time span as long as that between the divergence of prokaryotes from eukaryotes. This enabled us to draw phylogenetic trees for species whose phylogeny cannot be easily reconstructed from the fossil record. Our calculation of the times of divergence of the various organelle-specific enzymes led us to hypothesize that the pea and bean chloroplast genes for these enzymes originated from the duplication of nuclear genes as a result of the different metabolic needs of the various species. Our data indicate that the duplication of plastid glutamine synthetase genes occurred long after the endosymbiotic events that produced the organelles themselves.

DNA microenvironments and the molecular clock.

A few years ago we presented a stationary Markov model of gene evolution according to which only homologous genes from not too divergent species obeying the condition of being stationary may behave as reliable molecular clocks. A compartmentalized model of the nuclear genome in which the genes are distributed in compartments, the isochores, defined by their G + C content has been proposed recently. We have found that only homologous gene pairs that are stationary, and belong to the same isochore, can be used consistently for the determination of phylogeny and base substitution rate. In particular, for the rodent-human couple, only about half of the homologous gene pairs are stationary. Stationary genes evolve at the third silent codon position with the same velocity independent of the genes and base composition. By contrast, nonstationary genes display apparent rate values (pseudovelocities) that are significantly higher. Our results cast doubt upon recent claims of a large acceleration in the rate of molecular evolution in rodents.

A simple quantitative model of the molecular clock.

We present the ideas, and their motivation, at the basis of a simple model of nucleic acid evolution: the stationary Markov process, or Markov clock. After a brief review of its relevant mathematical properties, the Markov clock is applied to nucleotide sequences from mitochondrial and nuclear genes of different species. Particular emphasis is given to the necessity of carrying out a correct statistical analysis, which allows us to check quantitatively the applicability of our model. We find evidence that the Markov clock ticks in many different processes, and that its limitations can be understood in terms of a simple idea that we call the "base-drift" hypothesis. This hypothesis correlates the deviations from the stationarity of the Markov process to the evolutionary distance dAB(p) of two species A and B, relative to the process P. We conclude by discussing the implications of our findings for future work.

Transition and transversion rate in the evolution of animal mitochondrial DNA.

We present a further application of the stochastic model previously described (Lanave et al., 1984, 1985) for measuring the nucleotide substitution rate in the mammalian evolution of the mitochondrial DNA (mtDNA). The applicability of this method depends on the validity of "stationarity conditions" (equal nucleotide frequencies at first, second and third silent codon positions in homologous protein coding genes). In the comparison of homologous sequences satisfying the stationarity condition at the silent sites, only the four codon families (quartets) for which both transitions and transversions are silent at the third position are considered here. This has allowed us to estimate the transition and transversion rates for any pair of species. We have analyzed the third silent codon position of the triplet rat-mouse-cow, of a series of slightly divergent primates and of two Drosophila species. In terms of two external dating input we have then determined the phylogenetic trees for rat, mouse, and cow as well as for a number of primates including man. The phylogenetic tree that we have derived for the triplet rat, mouse and cow agrees with that we had previously determined by analyzing the first, second and third silent codon positions (in both duets and quartets) of mt genes (Lanave et al., 1985). For primates our method leads to the following branching order from the oldest to the most recent: Gibbon, Orangutan, Gorilla, Chimpanzee and Man. In absolute time, fixing the distance Chimpanzee-Man as 5 million years (Myr) we estimate the dating of the divergence nodes as: Gorilla 7 Myr; Orangutan 16 Myr; Gibbon 20 Myr. In all cases analyzed, the transition rate has been found to be substantially higher than the transversion rate. Moreover we have found that the transition/transversion ratio is different in the various lineages. We suggest that this fact is probably related to the nucleotide frequencies at the third silent codon position.

Multisequence comparisons in protein coding genes. Search for functional constraints.

A very powerful method for detecting functional constraints operative in biological macromolecules is presented. This method entails performing a base permanence analysis of protein coding genes at each codon position simultaneously in different species. It calculates the degree of permanence of subregions of the gene by dividing it into segments, c codons long, counting how many sites remain unchanged in each segment among all species compared. By comparing the base permanence among several sequences with the expectations based on a stochastic evolutionary process, gene regions showing different degrees of conservation can be selected. This means that wherever the permanence deviates significantly from the expected value generated by the simulation, the corresponding regions are considered "constrained" or "hypervariable". The constrained regions are of two types: alpha and beta. The alpha regions result from constraints at the amino acid level, whereas the beta regions are those probably involved in "control" processing. The method has been applied to mitochondrial genes coding for subunit 6 of the ATPase and subunit 1 of the cytochrome oxidase in four mammalian species: human, rat, mouse, and cow. In the two mitochondrial genes a few regions that are highly conserved in all codon positions have been identified. Among these regions a sequence, common to both genes, that is complementary to a strongly conserved region of 12S rRNA has been found. This method can also be of great help in studying molecular evolution mechanisms.