A purine-pyrimidine motif verifying an identical presence in almost all gene taxonomic groups

A statistical parameter identifies, with a high degree of significance, a motif which is present in protein-coding sequences of eukaryotes, prokaryotes, chloroplasts, mitochondria, viral introns, ribosomal RNA genes, and transfer RNA genes. The random probability of occurrence of such a situation is 10(-12). This motif has the following properties: (i) its significant presence in almost all present-day genes explains why it can be considered as primitive oligonucleotide, (ii) its nucleotide order is: YRY (N)6YRY, R being a purine base, Y a pyrimidine one and N any base, (iii) its length and its terminal trinucleotides YRY suggest a primordial function related to the spatial structure of the DNA sequences. This motif is found in some viral protein-coding genes, but not in eukaryotic introns.     

A model of DNA sequence evolution

Statistical studies of gene populations on the purine/pyrimidine alphabet have shown that the mean occurrence probability of thei-motif YRY(N) _i YRY (R=purine, Y=pyrimidine, N=R or Y) is not uniform by varyingi in the range [1,99], but presents a maximum ati=6 in the following populations: protein coding genes of eukaryotes, prokaryotes, chloroplasts and mitrochondria, and also viral introns, ribosomal RNA genes and transfer RNA genes (Arquès and Michel, 1987b,J. theor. Biol. 128, 457–461). From the “universality” of this observation, we suggested that the oligonucleotide YRY(N)₆ is a primitive one and that it has a central function in DNA sequence evolution (Arquès and Michel, 1987b,J. theor. Biol. 128, 457–461). Following this idea, we introduce a concept of a model of DNA sequence evolution which will be validated according to a shema presented in three parts. In the first part, using the last version of the gene database, the YRY(N)₆YRY preferential occurrence (maximum ati=6) is confirmed for the populations mentioned above and is extended to some newly analysed populations: chloroplast introns, chloroplast 5′ regions, mitochondrial 5′ regions and small nuclear RNA genes. On the other hand, the YRY(N)₆YRY preferential occurrence and periodicities are used in order to classify 18 gene populations. In the second part, we will demonstrate that several statistical features characterizing different gene populations (in particular the YRY(N)₆YRY preferential occurrence and the periodicities) can be retrieved from a simple Markov model based on the mixing of the two oligonucleotides YRY(N)₆ and YRY(N)₃ and based on the percentages of RYR and YRY in the unspecified trinucleotides (N)₃ of YRY(N)₆ and YRY(N)₃. Several properties are identified and prove in particular that the oligonucleotide mixing is an independent process and that several different features are functions of a unique parameter. In the third part, the return of the model to the reality shows a strong correlation between reality and simulation concerning the presence of large alternating purine/pyrimidine stretches and of periodicities. It also contributes to a greater understanding of biological reality, e.g. the presence or the absence of large alternating purine/pyrimidine stretches can be explained as being a simple consequence of the mixing of two particular oligonucleotides. Finally, we believe that such an approach is the first step toward a unified model of DNA sequence evolution allowing the molecular understanding of both the origin of life and the actual biological reality.     

Analytical expression of the purine/pyrimidine codon probability after and before random mutations

Recently, we proposed a new model of DNA sequence evolution (Arquès and Michel. 1990b.Bull. math. Biol. 52, 741–772) according to which actual genes on the purine/pyrimidine (R/Y) alphabet (R=purine=adenine or guanine, Y=pyrimidine=cytosine or thymine) are the result of two successive evolutionary genetic processes: (i) a mixing (independent) process of non-random oligonucleotides (words of base length less than 10: YRY(N)₆, YRYRYR and YRYYRY are so far identified; N=R or Y) leading to primitive genes (words of several hundreds of base length) and followed by (ii) a random mutation process, i.e. transformations of a base R (respectively Y) into the base Y (respectively R) at random sites in these primitive genes. Following this model the problem investigated here is the study of the variation of the 8 R/Y codon probabilities RRR,..., YYY under random mutations. Two analytical expressions solved here allow analysis of this variation in the classical evolutionary sense (from the past to the present, i.e. after random mutations), but also in the inverted evolutionary sense (from the present to the past, i.e. before random mutations). Different properties are also derived from these formulae. Finally, a few applications of these formulae are presented. They prove the proposition in Arquès and Michel (1990b.Bull. math. Biol. 52, 741–772), Section 3.3.2, with the existence of a miximal mean number of random mutations per base of the order 0.3 in the protein coding genes. They also confirm the mixing process of oligonucleotides by excluding the purine/pyrimidine contiguous and alternating tracts from the formation process of primitive genes.     

A simulation of the genetic periodicities modulo 2 and 3 with processes of nucleotide insertions and deletions.

Recently, a new genetic process termed RNA editing has been identified showing insertions and deletions of nucleotides in particular RNA molecules. On the other hand, there are a few non-random statistical properties in genes: in particular, the periodicity modulo 3 (P3) associated with an open reading frame, the periodicity modulo 2 (P2) associated with alternating purine/pyrimidine stretches, the YRY(N)6YRY preferential occurrence (R = purine = adenine or guanine, Y = pyrimidine = cytosine or thymine, N = R or Y) representing a "code" of the DNA helix pitch, etc. The problem investigated here is whether a process of the type RNA editing can lead to the non-random statistical properties commonly observed in genes. This paper will show in particular that: The process of insertions and deletions of mononucleotides in the initial sequence [YRY(N)3]* [series of YRY(N)3] can lead to the periodicity modulo 2 (P2). The process of insertions and deletions of trinucleotides in the initial sequence [YRY(N)6]* [series of YRY(N)6] can lead to the periodicity modulo 3 (P3) and the YRY(N)6YRY preferential occurrence. Furthermore, these two processes lead to a strong correlation with the reality, namely the mononucleotide insertion/deletion process, with the 5' eukaryotic regions and the trinucleotide insertion/deletion process, with the eukaryotic protein coding genes.     

Cladistical analysis of primitive G-band sequences for the karyotype of the ancestor of the Cricetidae complex of rodents

Homologous segments identified by G-banding sequences of chromosomes of Peromyscus boylii, Neotoma micropus, Oryzomys capito, (Family Cricetidae) Rattus norvegicus, Melomys burtoni, and Apodemus sylvaticus (Family Muridae) were used to hypothesize a chromosomal condition for the cricetid ancestor. A critical assumption in proposing the primitive G-banding sequences for a given chromosome is that if the outgroup and ingroup taxa have a specific sequence, then the ancestor of the ingroup taxa also had that same sequence. Using this methodology, (chromosome numbers refer to proposed homology to the standardized karyotype for Peromyscus), we propose that: (1) the primitive banding pattern of chromosome 1 was identical to that of Neotoma; (2) the primitive patterns of chromosomes 2, 3, 4, 6, 7, 8, 9, 10, 11, and 12 were primitive banding patterns of 5 and 13 were undetermined; (4) a major portion of the banding patterns of 14 and X were present in the ancestral karyotype. Only the largest 14 autosomes and X were examined because the smaller elements had insufficient G-band definition to ensure reasonable accuracy. The karyotype ancestral to that of Peromyscus, Neotoma, and Oryzomys may be as above and the banding patterns of 5, 13, and 14 were acrocentric and identical to those shown for Peromyscus, Neotoma, and Oryzomys (Fig. 1). In the primitive karyotype, heterochromatin (C-band material) was probably limited to the centromeric regions. If the primitive karyotype is as described above, then it is possible to determine the direction, type, and magnitude of chromosomal evolution evident in the various cricetid lineages. Based on the available data, radiation from the ancestral cytotype is characterized by a nonrandom distribution of types of chromosomal changes. Within many genera, more rearrangements occur in the 14 largest autosomal chromosomes of some congeneric species than distinguish the proposed primitive conditions for the genera Peromyscus, Neotoma, and Oryzomys. It would appear that the extensive morphological radiation from the primitive cricetid ancestor as indicated by the presence of over 100 surviving genera within the family, was not accompanied by extensive karyotypic changes. The magnitude of chromosomal variation that accompanies speciation in these genera appears to range from no detectable chromosomal evolution to a radical reorganization of the genome.     

A speculation on the origin of protein synthesis.

It is suggested that protein sythesis may have begun without even a primitive ribosome if the primitive tRNA could take up two configuration and could bind to the messenger RNA with five base-pairs instead of the present three. This idea would impose base sequence restriction on the early messages and on the early genetic code such that the first four amino acids coded were glycine, serine, aspartic acid and aspargine. A possible mechanism is suggested for the polymerization of the early message.     

Genomic analysis of Hox clusters in the sea lamprey Petromyzon marinus

The sea lamprey Petromyzon marinus is among the most primitive of extant vertebrates. We are interested in the organization of its Hox gene clusters, because, as a close relative of the gnathostomes, this information would help to infer Hox cluster organization at the base of the gnathostome radiation. We have partially mapped the P. marinus Hox clusters using phage, cosmid, and P1 artificial chromosome libraries. Complete homeobox sequences were obtained for the 22 Hox genes recovered in the genomic library screens and analyzed for cognate group identity. We estimate that the clusters are somewhat larger than those of mammals (roughly 140 kbp vs. 105 kbp) but much smaller than the single Hox cluster of the cephalochordate amphioxus (at more than 260 kb). We never obtained more than three genes from any single cognate group from the genomic library screens, although it is unlikely that our screen was exhaustive, and therefore conclude that P. marinus has a total of either three or four Hox clusters. We also identify four highly conserved non-coding sequence motifs shared with higher vertebrates in a genomic comparison of Hox 10 genes.     

A speculation on the origin of protein synthesis

It is suggested that protein synthesis may have begun without even a primitive ribosome if the primitive tRNA could take up two configurations and could bind to the messenger RNA with five base-pairs instead of the present three. This idea would impose base sequence restriction on the early messages and on the early genetic code such that the first four amino acids coded were glycine, serine, aspartic acid and aspargine. A possible mechanism is suggested for the polymerization of the early message.This paper is dedicated to the memory of Dr. Aharon Katzir.     

Studies of Pollen Morphology and Its Systematic Position in the Order Piperales

The pollen morphology of 25 species and 10 genera of Piperales (Chloranthaceae, Piperaceae and Saururaceae) has been examined under light microscope, of which 7 species were observed under scanning electron microscope and 1 species, Hedyosmum orentale Merr. & Chun transmision electron microseope. Three principal types of pollen were found: anasulcate (mostly) (sometime trichotomosulcate), inaperturate (partly) and multicolpoidate (partly). The present article has discussed the palynological data mainly in relation to the classification and the systematic position of Cbloranthaceae and also deals with the systematic position of the order Piperales. The present author agrees to put the family Chloranthaceae into the order Piperales. Because this family differs from Piperaceae and Saururaceae in pollen morphology, therefore, Chloranthaceae should raise to the level of suborder. Among three families of the order Piperales, the present author considers Chloranthaceae to be the most primitive family, on account of the following reasons: 1. The family Chloranthaceae shows the characteristics of primitive entomophilous plants in the sculpture of exine, while in the other two families, Piperaceae and Saururaceae, their exine is almost smooth and represents wind-pollenated plants; 2. Pollen of the family Chloranthaceae are larger than those of Piperaceae and Saururaceae; 3. The fossil pollen Clavatipollenites has been proved to be one of the most primitive angiosperms on the earth, that it is known, it occurred in the early Cretaceous, and at that time ferns and gymnosperms were predominant, while the Chloran- thaceae has already existed at that time; 4. Sarcandra of Chloranthaceae possesses the characters of a vesselless secondary xylem and a delayed development of embryo. Thus, Chloranthaceae would be considered as the most primitive family in the order Piperales. The systematic position of the order Piperales is also discussed. Itutehinson makes a point that order Ranales is more primitive than Piperales, and his system is arranged in the following order: Ranales → Piperales → to climax family Chloranthaceae. This view-point, however, is net supported by the palynological data. Pollen morphology shows that Piperales is more primitive than Ranales, because the pollen in Piperales possess the ancient aperture type of Pteridospermes, i.e., the type of anasulcate aperture is prevailing in Piperales, moreover, pollen grains of Ranales are mainly tricolpate type, and tricolpate pollen is a characteristic of typical angiosperms. In addition, the Piperales possesses a series of characters that are common among monocots, but rare among dicots. As the divergence between dicer and monocot took place in the early Cretaceous, their ancestor possesses common chararcters both of dicots and monocots while the extant Piperales still possess many characters of monocots that indicate it is much nearer to the point of divergence, and it explains that the Piperales is closely related to the ancestor of monocots and dicers Piperales, therefore, is more primitive than Ranales.     

Prospects for a complete molecular map of the human genome

Linkage maps are limited by the number of recombinations that can be scored in human pedigrees to a resolution of ca. 1 centimorgan (relative distance between genes on a chromosome having a crossover value of 1%) which is estimated to be about 10 megabases. Molecular maps can be formed at any resolution down to the base sequence. To complement the linkage approach, the most useful molecular map would be one that helped to locate disease loci, by using restriction fragment length polymorphisms (RFLPS) and accurate localization of recombinations, and which then helped to find candidate genes in this region, by providing the positions of coding sequences. This paper discusses the appropriate form and scale of such a map, how it can be produced with methods now available, and the most efficient strategy for building the map, based on present knowledge of the organization of the human genome.     

Molten earth and the origin of prebiological molecules.

Evidence for the molten Earth at its accretion time has been accumulated through the geochemical investigations and the observations of the surfaces of planets by space probes such as Venera 8, Mariner 9, Surveyor, Luna, and Apollo. The primitive terrestrial atmosphere might have been derived from the volcanic gases, as suggested by Rubey, but of a higher temperature than so far assumed. A thermochemical calculation of the composition of the volcanic gas suggests the following possibilities: (1) Large amounts of H2 and CO were present in the primitive atmosphere. This gives a theoretical basis for the HCN-production experiment by Abelson. (2) HCHO and NH3 existed in the primitive oceans, of the amount comparable with the weight of the present biosphere. (3) Plenty of NO3-, SO4, and PO4 were expected in the primitive oceans. The NO3- ions might have been useful for the nitrate respiration advocated by Egami. In an appendix, it is argued, on ;he basis of the observational evidence of the exospheric temperatures of planets by space probes, that a highly reducing atmosphere would (if it existed on the primitive Earth) have disappeared very quickly due to the thermal escape of hydrogen from its exosphere.     

Ciona intestinalis cDNA projects: expressed sequence tag analyses and gene expression profiles during embryogenesis

Ascidians are primitive chordates. Their fertilized egg develops quickly into a tadpole-type larva, which consists of a small number but distinct types of cells, including those of epidermis, central nervous system with two sensory organs, endoderm and mesenchyme in the trunk, and notochord and muscle in the tail. This configuration of the ascidian tadpole is thought to represent the most simplified and primitive chordate body plan. In addition, the free-swimming and non-feeding larvae metamorphose into sessile and filter-feeding adults. The genome size of Ciona intestinalis is estimated to be about 160 Mb, and the number of genes approximately 15,500. The present Ciona cDNA projects focused on gene expression profiles of fertilized eggs, 32-110-cell stage embryos, tailbud embryos, larvae, and young adults. Expressed sequence tags (ESTs) of the 5'-most end and 3'-most end of more than 3000 clones were determined at each developmental stage, and the clones were categorized into independent clusters using the 3'-end sequences. Nearly 1000 clusters of them were then analyzed in detail of their sequences against a BLASTX search. This analysis demonstrates that, on average, half of the clusters showed proteins with sequence similarities to known proteins and the other half did not show sequence similarities to known proteins. Genes with sequence similarities were further categorized into three major subclasses, depending on their functions. Furthermore, the expression profiles of all of the clusters were analyzed by whole-mount in situ hybridization. This analysis highlights gene expression patterns characteristic to each developmental stage. As a result, the present study provides many new molecular markers for each of the tissues and/or organs that constitutes the Ciona tailbud embryo. This sequence information will be used for further comparative genome studies to explore molecular mechanisms involved in the formation of one of the most primitive chordate body plans. All of the data fully characterized may be viewed at the web site http://ghost.zool.kyoto-u.ac.jp.     

Phylogenetic relationships and the primitive X chromosome inferred from chromosomal and satellite DNA analysis in Bovidae

The early phylogeny of the 137 species in the Bovidae family is difficult to resolve; knowledge of the evolution and relationships of the tribes would facilitate comparative mapping, understanding chromosomal evolution patterns and perhaps assist breeding and domestication strategies. We found that the study of the presence and organization of two repetitive DNA satellite sequences (the clone pOaKB9 from sheep, a member of the 1.714 satellite I family and the pBtKB5, a 1.715 satellite I clone from cattle) on the X and autosomal chromosomes by in situ hybridization to chromosomes from 15 species of seven tribes, was informative. The results support a consistent phylogeny, suggesting that the primitive form of the X chromosome is acrocentric, and has satellite I sequences at its centromere. Because of the distribution of the ancient satellite I sequence, the X chromosome from the extant Tragelaphini (e.g. oryx), rather than Caprini (sheep), line is most primitive. The Bovini (cow) and Tragelaphini tribes lack the 1.714 satellite present in the other tribes, and this satellite is evolutionarily younger than the 1.715 sequence, with absence of the 1.714 sequence being a marker for the Bovini and Tragelaphini tribes (the Bovinae subfamily). In the other tribes, three (Reduncini, Hippotragini and Aepycerotini) have both 1.714 and 1.715 satellite sequences present on both autosomes and the X chromosome. We suggest a parallel event in two lineages, leading to X chromosomes with the loss of 1.715 satellite from the Bovini, and the loss of both 1.714 and 1.715 satellites in a monophyletic Caprini and Alcelaphini lineage. The presence and X chromosome distribution of these satellite sequences allow the seven tribes to be distributed to four groups, which are consistent with current diversity estimates, and support one model to resolve points of separation of the tribes.     

Haemonchus contortus: evaluation of two signal sequence trapping systems for detection of secreted molecules

Given that signal sequences between secreted proteins of different species can be interchanged, it is reasonable to expect that both mammalian and yeast signal sequence trapping (SST) systems would secrete Haemonchus contortus proteins with similar efficiency and quality. To determine if H. contortus cDNAs that contain a signal sequence could re-establish secretion of a reporter protein, mammalian and yeast SST vectors were designed, 10 H. contortus genes selected, and their respective cDNAs cloned into these two SST vectors. The selected molecules included genes known to code for excretory/secretory or membrane-bound proteins as potential test 'positives', and genes known to code for non-secreted proteins as test 'negatives'. While differentiation between secretion and non-secretion was evident in both systems, the results indicated greater efficiency was achieved when the mammalian system was used. Therefore, mammalian SST using COS cells would be a more useful tool to screen H. contortus cDNA libraries for potential secreted and type-1 integral membrane proteins than yeast SST.     

Progress and challenges in studies of the evolution of development

Plant evolutionary developmental genetics (EDG) has made considerable progress over the last decade. This is in part due to the accumulation of large amounts of sequence data that have provided robust organismal phylogenies and, increasingly, broad assessments of molecular evolution. Attempts to use primary sequence data to identify genes that have changed function in evolutionary time have not been as successful as initially hoped. The coding sequences of most genes, which are more amenable to statistical analysis than are regulatory sequences, are generally under purifying selection, as would be expected if much evolutionary change is the result of changes in cis-regulatory sequences. Sequence-based analysis of the regulatory sequences themselves remains difficult. Comparative studies of gene expression have been useful to identify genes whose developmental role may have changed in evolutionary time and will be critical to the future development of EDG. Such studies can be used to test hypotheses of gene function. Transformation experiments are often illuminating, but can be hard to interpret, particularly if genes from multiple species are all placed into a single heterologous system such as Arabidopsis. The ideal experiment would be a gene swap or promoter swap between two species, but this awaits development of good transformation systems. The immediate need for EDG is studies of gene expression on a massive scale, far broader than any studies undertaken to date.