A new method for calculating evolutionary substitution rates

In 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 techniques. 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 X 10(-8) nucleotide substitutions per site per year. We have also determined the divergence time ratios between the couples mouse-cow/rat-mouse and rat-cow/rat-mouse. In both cases this value is approximately 1.4.

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.

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.

ACNUC--a portable retrieval system for nucleic acid sequence databases: logical and physical designs and usage

ACNUC is a database structure and retrieval software for use with either the GenBank or EMBL nucleic acid sequence data collections. The nucleotide and textual data furnished by both collections are each restructured into a database that allows sequence retrieval on a multi-criterion basis. The main selection criteria are: species (or higher order taxon), keyword, reference, journal, author, and organelle; all logical combinations of these criteria can be used. Direct access to sequence regions that code for a specific product (protein, tRNA or rRNA) is provided. A versatile extraction procedure copies selected sequences, or fragments of them, from the database to user files suitable to be analysed by user-supplied application programs. A detailed help mechanism is provided to aid the user at any time during the retrieval session. All software has been written in FORTRAN 77 which guarantees a high degree of transportability to minicomputers or mainframes.

The nucleotide sequence of the large ribosomal RNA gene and the adjacent tRNA genes from rat mitochondria

We have sequenced the Eco R(1) fragment D from rat mitochondrial DNA. It contains one third of the tRNA (Val) gene (the remaining part has been sequenced from the 3' end of the Eco R(1) fragment A) the complete gene for the large mt 16S rRNA, the tRNA (Leu) gene and the 5' end of an unidentified reading frame. The mt gene for the large rRNA from rat has been aligned with the homologous genes from mouse and human using graphic computer programs. Hypervariable regions at the center of the molecule and highly conserved regions toward the 3' end have been detected. The mt gene for tRNA Leu is of the conventional type and its primary structure is highly conserved among mammals. The mt gene for tRNA(Val) shows characteristics similar to those of other mt tRNA genes but the degree of homology is lower. Comparative studies confirm that AGA and AGG are read as stop codons in mammalian mitochondria.

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.

Evidence for two distinct subgroups of alfalfa mosaic virus (AMV) from france and italy and their relationships with other AMV strains Brief report

The nucleotide sequence of the putative coat protein open reading frame of seven previously uncharacterized AMV strains from Italy and France was determined and aligned with comparable sequences of other AMV strains (425 L, 425 M, YSMV, S, VRU, 15/64 and Da). The data set of AMV sequences was used to determine phylogenetic relationships by both a stochastic (stationary Markov model) and a deterministic method (maximum-parsimony) of analysis. The topology of the trees obtained with the two methods was essentially the same showing that all AMV strains clustered in two monophyletic groups. Close clustering of Italian strains in subgroup I and of French strains in subgroup II seems to suggests the effect of geographic distinctiveness of evolutionary dynamics of these AMV strains. This separation did not correlate with differences in host range or symptoms (necrotic or non necrotic) induced in tomato but rather it reflected variations in the amino acid sequence of their CP, which might be related to structural properties of virus particles. A simple and rapid procedure based on the reverse transcriptase-polymerase chain reaction (RT-PCR) followed by ezymatic digestion (RFLP) was developed to identify and classify AMV isolates into the two subgroups. The method applied to a number of other AMV isolates from Italy and France supported their division in two distinct subgroups. This RT-PCR RFLP method may be useful way to investigate the dynamics of AMV populations in nature.

Evolutionary dynamics of cucumber mosaic virus satellite RNA during natural epidemics in Italy

The evolutionary dynamics of 22 variants of cucumber mosaic virus satellite RNA (CMV satRNA) isolated in Italy during virus epidemics from 1988 to 1993 were investigated on the basis of their primary structure and biological properties. Most of the variants were amplified from total nucleic acid preparations extracted from field-infected plants, thus representing wild isolates of CMV satRNA. Eleven variants were associated with subgroup II CMV strains, 10 with subgroup I and 1 with a mixed infection by both strains. When inoculated onto tomato seedlings, the variants induced the phenotype (necrogenic or ameliorative) predicted by their nucleotide sequence. Phylogenetic relationships between the satRNA variants were determined using the stationary Markov model, a stochastic model for evolution. For each satRNA, the Markov analysis gave a good correlation between position in the phylogenetic tree and biological properties. The variants with ameliorative and necrogenic phenotypes in tomato followed two different evolutionary dynamics in nature. Tfn-satRNA, a 390-nt-long molecule, followed a third type of evolutionary dynamic far apart from that of the shorter satRNA molecules (i.e., those in the 334- to 340-nt-length class). Average values of the mean constant rate of nucleotide substitutions/site (Ksubs/site) indicated that in nature the variants tend to keep their heterogeneity unchanged from one epidemic episode to the other, even if the outbreaks occur in places very far from each other. This seems to be in agreement with the proposed maintenance of a functional molecular structure as a constraint to CMV satRNA evolution.

Evolution of TRG clusters in cattle and sheep genomes as drawn from the structural analysis of the ovine TRG2@ locus

The availability of genomic clones representative of the T-cell receptor constant gamma (TRGC) ovine genes enabled us to demonstrate, by fluorescent in situ hybridization (FISH) on cattle and sheep metaphases, the presence of two T-cell receptor gamma (TRG1@ and TRG2@) paralogous loci separated by at least five chromosomal bands on chromosome 4. Only TRG1@ is included within a region of homology with human TRG locus on chromosome 7, thus TRG2@ locus appears to be peculiar to ruminants. The structure of the entire TRG2@ locus, the first complete physical map of any ruminant animal TCR gamma locus, is reported here. The TRG2@ spans about 90 kb and consists of three clusters that we named TRG6, TRG2, and TRG4, according to the constant genes name. Phylogenetic analysis has highlighted the correlation between the grouping pattern of cattle and sheep variable gamma (TRGV) genes and the relevant TRGC; variable (V), joining (J), and constant (C) rearrange to be found together in mature transcripts. The simultaneous results on the TRG2@ locus molecular organization in sheep and on the phylogenetic analysis of cattle and sheep V expressed sequences indicate that at least six TRG clusters distributed in the two loci are present in these ruminant animals. The inferred evolution of TRG clusters in cattle and sheep genomes is consistent with a scenario where a minimal ancient cluster, containing the basic structural scheme of one V, one J, and one C gene, has undergone a process of duplication and intrachromosomal transposition.

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.

Mammalian genes as molecular clocks?

In analyzing the silent nucleotide substitutions in some mammalian mitochondrial mRNA coding genes, we had found that the frequency of each of the four nucleotides in rat, mouse, and cow, but not in humans, is the same in the silent third codon position (Lanave C, Preparata G, Saccone C, Serio G (1984) J Mol Evol 20:86-93). Because our findings for these three species were compatible with a stationary Markov process for the evolution of nucleotide sequences, we applied such a model to calculate the effective evolutionary silent substitution rate (vs) and the divergence times among the species. In this paper we have analyzed the first and second codon positions in the same mammalian mitochondrial genes. We found that in the first and second codon positions the human mitochondrial genes satisfy the stationarity conditions. This has allowed us to use the stochastic model mentioned above to calculate the divergence times among mouse, rat, cow, and human. Furthermore, we have analyzed the silent substitution rate in one nuclear gene for these four mammals. We found that in this gene the effective silent substitution rate is about 3 times lower than in mitochondrial genes, and that humans are in this case stationary with respect to the other three mammals in the third codon position as well. Application of our Markov model to this latter gene yields divergence times consistent with our previous determinations.

The nucleotide sequences of several tRNA genes from rat mitochondria: common features and relatedness to homologous species

We have determined the nucleotide sequences of thirteen rat mt tRNA genes. The features of the primary and secondary structures of these tRNAs show that those for Gln, Ser, and f-Met resemble, while those for Lys, Cys, and Trp depart strikingly from the universal type. The remainder are slightly abnormal. Among many mammalian mt DNA sequences, those of mt tRNA genes are highly conserved, thus suggesting for those genes an additional, perhaps regulatory, function. A simple evolutionary relationship between the tRNAs of animal mitochondria and those of eukaryotic cytoplasm, of lower eukaryotic mitochondria or of prokaryotes, is not evident owing to the extreme divergence of the tRNA sequences in the two groups. However, a slightly higher homology does exist between a few animal mt tRNAs and those from prokaryotes or from lower eukaryotic mitochondria.

Artiodactyl emergence is accompanied by the birth of an extensive pool of diverse germline TRDV1 genes

Molecular cloning of cDNA from gamma/delta T cells has shown that in sheep, the variable domain of the delta chain is chiefly determined by the expression of the TRDV1 subgroup, apparently composed of a large number of genes. There are three other TRDV subgroups, but these include only one gene each. To evaluate the extent and the complexity of the genomic TRDV repertoire, we screened a sheep liver genomic library from a single individual of the Altamurana breed and sheep fibroblast genomic DNA from a single individual of the Gentile di Puglia breed. We identified a total of 22 TRDV1 genes and the TRDV4 gene. A sequence comparison between germline and the rearranged genes indicates that, in sheep, the TRDV repertoire is generated by the VDJ rearrangement of at least 40 distinct TRDV1 genes. All germline TRDV1 genes present a high degree of similarity in their coding as well as in 5' and 3' flanking regions. However, a systematic analysis of the translation products reveals that these genes present a broadly different and specific repertoire in the complementarity-determining regions or recognition loops, allowing us to organize the TRDV genes into sets. We assume that selection processes operating at the level of ligand recognition have shaped the sheep TRDV germline repertoire. A phylogenetic study based on a sequence analysis of the TRDV genes from different mammalian species shows that the diversification level of these genes is higher in artiodactyl species compared to humans and mice.

Organization, structure and evolution of 41kb of genomic DNA spanning the D-J-C region of the sheep TRB locus

A genomic region of 41,045 bp encompassing the 3'-end of the sheep T cell receptor beta chain was sequenced. Extensive molecular analysis has revealed that this region retains a unique structural feature for the presence of a third D-J-C cluster, never detected in any other mammalian species examined so far. A total of 3 TRBD, 18 TRBJ and 3 substantially identical TRBC genes were identified in about 28kb. At 13kb, downstream from the last TRBC gene, in an inverted transcriptional orientation, lies a TRBV gene. Sequence comparison and phylogenetic analyses have demonstrated that the extra D-J-C cluster originated from an unequal crossing over between the two ancestral TRBC genes. Interspersed repeats spanning 22.2% of the sequence, contribute to the wider size of the sheep TRB locus with respect to the other mammalian counterparts, without modifying the general genomic architecture. The nucleotide and predicted amino acid sequences from peripheral T cells cDNA clones indicated that the genes from cluster 3 are fully implicated in the beta chain recombination machinery. Closer inspections of the transcripts have also shown that inter-cluster rearrangements and splice variants, involving the additional cluster, increase the functional diversity of the sheep beta chain repertoire.

Lectin-related resistance factors against bruchids evolved through a number of duplication events

Abundant lectin-related proteins found in common beans ( Phaseolus vulgaris L.) have been shown to confer resistance against the larvae of a number of bruchid species. Genes encoding for these proteins are members of the lectin multigene family, the most representative components being arcelins, phytohemagglutinins and alpha-amylase inhibitors. Arcelins have been described in seven variants, some of which are resistance factors against the Mexican bean weevil ( Zabrotes subfasciatus), a major bean predator. In this study the isolation and sequencing of arcelin genes from wild P. vulgaris genotypes, containing Arc3 and Arc7 variants, is reported, and similarities and evolutionary relationships among the seven known arcelins are described. The evolutionary analysis shows that arcelins 3 and 4 cluster together and are the most-ancient variants. A duplication event gave rise to two additional clusters, one comprising arcelins 1, 2 and 6 and separated from the cluster of arcelins 5 and 7. A multiple number of arcelin genes were found in arcelin 3 and 4 genotypes indicating that more than one type of arcelin gene may be present in the same locus. Some of these sequences are reminiscent of ancient duplication events in arcelin evolution demonstrating that arcelins have evolved through multiple duplications. A further aim of this paper was to better understand and describe the evolution of the entire lectin multigene family. Beside arcelins, a number of other types of sequences, such as putative lectins and sequences not easily classifiable, were found in genotypes containing Arc3 and Arc4. These results, together with the evolutionary analysis, indicate that lectin loci are quite complex and confirm their origin by multiple duplication events.

Genomic organization of the sheep TRG1@ locus and comparative analyses of Bovidae and human variable genes

gammadelta T cells commonly account for 0.5%-5% of human (gammadelta low species) circulating T cells, whereas they are very common in chickens, and they may account for >70% of peripheral cells in ruminants (gammadelta high species). We have previously reported the ovine TRG2@ locus structure, the first complete physical map of any ruminant animal TCR locus. Here we determined the TRG1@ locus organization in sheep, reported all variable (V) gamma gene segments in their germline configuration and included human and cattle sequences in a three species comparison. The TRG1@ locus spans about 140 kb and consists of three clusters named TRG5, TRG3, and TRG1 according to the constant (C) genes. The predicted tertiary structure of cattle and sheep V proteins showed a remarkably high degree of conservation between the experimentally determined human Vgamma9 and the proteins belonging to TRG5 Vgamma subgroup. However systematic comparison of primary and tertiary structure highligthed that in Bovidae the overall conformation of the gammadelta TCR, is more similar to the Fab fragment of an antibody than any TCR heterodimer. Phylogenetic analysis showed that the evolution of cattle and sheep V genes is related to the rearrangement process of V segments with the relevant C, and consequentely to the appartenence of the V genes to a given cluster. The TRG cluster evolution in cattle and sheep pointed out the existence of a TRG5 ancient cluster and the occurrence of duplications of its minimal structural scheme of one V, two joining (J), and one C.

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.

Physical relationship between satellite I and II DNA in centromeric regions of sheep chromosomes

Fluorescence in situ hybridization (FISH) with probes representing sheep satellite I and satellite II DNAs shows a different distribution of the two repetitive DNA families in the centromeric region of most chromosomes. The single signal per chromosome produced by the satellite I probe suggests close proximity of this DNA family to the primary constriction. Satellite II produces two separate signals on the sister chromatids, and large blocks of satellite II DNA constitute most of the short arm of all acrocentric chromosomes. We have isolated and sequenced a phage clone containing a junction between discrete blocks of satellite I and satellite II sequences. The junction is characterized by an abrupt juxtaposition of arrays of the two satellites. The possibility that the peculiar structural features of this junction could have a functional significance is discussed.

The complete Tirant transposable element in Drosophila melanogaster shows a structural relationship with retrovirus-like retrotransposons

We have determined the structure and organization of Tirant, a retrotransposon of Drosophila melanogaster reported in literature to be responsible for four independent mutations. Tirant is a long terminal repeat (LTR) retrotransposon 8527bp long. It possesses three open reading frames (ORF) encoding Gag, Pol and Env proteins with a strong similarity with ZAM, a recently identified member of the gypsy class of retrovirus-like mobile elements. Molecular analysis of the Tirant genomic copies present in four D. melanogaster strains revealed that most of them are defective, non-autonomous elements that differ in the position and extension of the conserved internal portion. Defective elements lacking the Gag ORF but retaining the Env ORF are abundant in heterochromatin. Four discrete Tirant transcripts are observed during embryogenesis in the strain Oregon-R, the smaller of which, 1.8kb in size, originates from the splicing of a primary transcript and leads to a subgenomic RNA coding for the Env product.

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.

Time and biosequences

In this paper we discuss and demonstrate the importance of several factors relative to the relationship between time and evolution of biosequences. In both quantitative and qualitative measurements of the genetic distances, the compositional constraints of the nucleotide sequences play a very important role. We demonstrate that when homologous sequences significantly differ in base composition we get erratic branching order and/or wrong evaluation of the evolutionary rates. We must consider that every gene may have a different evolutionary dynamic along its sequence, generally linked to its functional constraints; this too can seriously affect its clock-like behavior. We report some cases showing how these factors can affect the quantitative measurements of the genetic distances of biosequences.

Lectin and lectin-related proteins in lima bean (Phaseolus lunatus L.) seeds: biochemical and evolutionary studies

Lectin-related polypeptides are a class of defence proteins found in seeds of Phaseolus species. In Lima bean (P. lunatus), these proteins and their genes have been well characterized in the Andean morphotype, which represents one of the two gene pools of this species. To study the molecular evolution of the lectin family in Lima bean we characterized the polypeptides belonging to this multigene family and cloned the genes belonging to the Mesoamerican gene pool. The latter gene pool contains components similar to those of the Andean pool, namely: an amylase inhibitor-like (AIL), an arcelin-like (ARL) lectin and the less abundant Lima bean lectin (LBL). These proteins originate from an ancestor gene of the lectin type which duplicated to yield the lectin gene and the progenitor of ARL and AIL. In this species. ARL represents an evolutionary intermediate form that precedes AIL. Phylogenetic analysis supports an Andean origin for Lima bean. The molecular evolutionary studies were extended to the genes of common bean and demonstrated that true lectin genes and the ancestor of lectin-related genes are the result of a duplication event that occurred before speciation. Lima and common bean followed different evolutionary pathways and in the latter species a second duplication event occurred that gave rise, in Mesoamerican wild genotypes, to arcelin genes.

Genetic characterization of Pompeii and Herculaneum Equidae buried by Vesuvius in 79 AD

DNA extracted from the skeletons of five equids discovered in a Pompeii stable and of a horse found in Herculaneum was investigated. Amino acid racemization level was consistent with the presence of DNA. Post-mortem base modifications were excluded by sequencing a 146 bp fragment of the 16S rRNA mitochondrial gene. Sequencing of a 370 bp fragment of mitochondrial (mt)DNA control region allowed the construction of a phylogenetic tree that, along with sequencing of nuclear genes (epsilon globin, gamma interferon, and p53) fragments, gave us the possibility to address some questions puzzling archaeologists. What animals-donkeys, horses, or crossbreeds-were they? And, given they had been evidently assigned to one specific job, were they all akin or were they animals with different mitochondrial haplotypes? The conclusions provided by molecular analysis show that the Pompeii remains are those of horses and mules. Furthermore one of the equids (CAV5) seems to belong to a haplotype, which is either not yet documented in the GenBank or has since disappeared. As its characteristics closely recall those of donkeys, which is the out group chosen to construct the tree, that appears to have evolved within the Equidae family much earlier than horses, this assumption seems to be nearer the truth.

Exon-intron organization of TRGC genes in sheep

A series of genomic clones derived from a sheep library were used to determine the germline configuration and the exon-intron organization of TRGC2, TRGC3, and TRGC4 genes. Based on the outcomes of molecular analysis, we compared and aligned the genomic sequences with the known complete cDNA sequences of sheep and deduced the exon-intron organization of TRGC genes in this ruminant animal, EX1, corresponding to the disulfide-linked constant domain, and EX3, corresponding to the transmembrane and cytoplasmatic domains, are similar in length in all genes. Conversely, the hinge-encoding EX2A, EX2B, and EX2C exons differ in number and length between genes, and EX2A contains the TTKPP motif irrespective of whether it occurs in single or triplicate form. The molecular data also indicate that at least one additional gene is present in sheep. Phylogenetic analysis grouped the ruminant TRGC genes in two clusters that could have emerged from two ancestral forms that underwent a series of duplications giving rise to the new sequences that were selected and then fixed in the ruminant lineages. A correlation between the cluster distribution in the phylogenetic tree of TRGC genes and their expression during fetal development is discussed.