Cloning and sequencing of cDNA encoding baboon liver alcohol dehydrogenase: evidence for a common ancestral lineage with the human alcohol dehydrogenase beta subunit and for class I ADH gene duplications predating primate radiation

The natural history of class I primate alcohol dehydrogenases includes gene duplication, gene loss, and gene conversion

BACKGROUND

Gene duplication is a source of molecular innovation throughout evolution. However, even with massive amounts of genome sequence data, correlating gene duplication with speciation and other events in natural history can be difficult. This is especially true in its most interesting cases, where rapid and multiple duplications are likely to reflect adaptation to rapidly changing environments and life styles. This may be so for Class I of alcohol dehydrogenases (ADH1s), where multiple duplications occurred in primate lineages in Old and New World monkeys (OWMs and NWMs) and hominoids.

METHODOLOGY/PRINCIPAL FINDINGS

To build a preferred model for the natural history of ADH1s, we determined the sequences of nine new ADH1 genes, finding for the first time multiple paralogs in various prosimians (lemurs, strepsirhines). Database mining then identified novel ADH1 paralogs in both macaque (an OWM) and marmoset (a NWM). These were used with the previously identified human paralogs to resolve controversies relating to dates of duplication and gene conversion in the ADH1 family. Central to these controversies are differences in the topologies of trees generated from exonic (coding) sequences and intronic sequences.

CONCLUSIONS/SIGNIFICANCE

We provide evidence that gene conversions are the primary source of difference, using molecular clock dating of duplications and analyses of microinsertions and deletions (micro-indels). The tree topology inferred from intron sequences appear to more correctly represent the natural history of ADH1s, with the ADH1 paralogs in platyrrhines (NWMs) and catarrhines (OWMs and hominoids) having arisen by duplications shortly predating the divergence of OWMs and NWMs. We also conclude that paralogs in lemurs arose independently. Finally, we identify errors in database interpretation as the source of controversies concerning gene conversion. These analyses provide a model for the natural history of ADH1s that posits four ADH1 paralogs in the ancestor of Catarrhine and Platyrrhine primates, followed by the loss of an ADH1 paralog in the human lineage.

Alcohol dehydrogenases: gene multiplicity and differential functions of five classes of isozymes

Mammalian alcohol dehydrogenases (ADHs) constitute an enzyme family of multiple forms (isozymes) which are differentially distributed throughout the body. Subunit types alpha, beta and gamma in dimeric combinations constitute the isozymes of human liver class I ADH, and are >94% homologous in structure. Human pi and chi subunits form homodimeric Class II and III ADH isozymes. pi-ADH is liver specific whereas chi-ADH is widely distributed throughout the body. A sixth human ADH subunit (designated mu or sigma), forming a new dimeric human stomach ADH, has been recently reported as Class IV ADH. Evidence for a seventh human ADH subunit has also been described, designated as Class V, the transcripts having been reported in the stomach and liver. All five classes of ADH represent isozymes which are homologous but exhibit at least 30% sequence differences in primary srtructure. Kinetic analyses of four of these classes of ADH indicated differential functions, serving either in the oxidative or reductive mode. Studies from various laboratories indicate the following respective functions: oxidation of aliphatic and aromatic alcohols-liver Class I and Class II, and stomach Class IV ADHs; reduction of peroxidic aldehydes-Classes I, II and IV; 'biogenic' alcohol oxidation-Classes I and II; and glutathione-dependent formaldehyde dehydrogenase-Class III.

Opossum alcohol dehydrogenases: Sequences, structures, phylogeny and evolution: evidence for the tandem location of ADH genes on opossum chromosome 5

BLAT (BLAST-Like Alignment Tool) analyses and interrogations of the recently published opossum genome were undertaken using previously reported rat ADH amino acid sequences. Evidence is presented for six opossum ADH genes localized on chromosome 5 and organized in a comparable ADH gene cluster to that reported for human and rat ADH genes. The predicted amino acid sequences and secondary structures for the opossum ADH subunits and the intron-exon boundaries for opossum ADH genes showed a high degree of similarity with other mammalian ADHs, and four opossum ADH classes were identified, namely ADH1, ADH3, ADH6 and ADH4 (for which three genes were observed: ADH4A, ADH4B and ADH4C). Previous biochemical analyses of opossum ADHs have reported the tissue distribution and properties for these enzymes: ADH1, the major liver enzyme; ADH3, widely distributed in opossum tissues with similar kinetic properties to mammalian class 3 ADHs; and ADH4, for which several forms were localized in extrahepatic tissues, especially in the digestive system and in the eye. These ADHs are likely to perform similar functions to those reported for other mammalian ADHs in the metabolism of ingested and endogenous alcohols and aldehydes. Phylogenetic analyses examined opossum, human, rat, chicken and cod ADHs, and supported the proposed designation of opossum ADHs as class I (ADH1), class III (ADH3), class IV (ADH4A, ADH4B and ADH4C) and class VI (ADH6). Percentage substitution rates were examined for ADHs during vertebrate evolution which indicated that ADH3 is evolving at a much slower rate to that of the other ADH classes.

Conservative evolution in duplicated genes of the primate Class I ADH cluster

Humans have seven alcohol dehydrogenase genes (ADH) falling into five classes. Three out of the seven genes (ADH1A, ADH1B and ADH1C) belonging to Class I are expressed primarily in liver and code the main enzymes catalyzing ethanol oxidization. The three genes are tandemly arrayed within the ADH cluster on chromosome 4 and have very high nucleotide similarity to each other (exons: >90%; introns: >70%), suggesting the genes have been generated by duplication event(s). One explanation for maintaining similarity of such clustered genes is homogenization via gene conversion(s). Alternatively, recency of the duplications or some other functional constraints might explain the high similarities among the genes. To test for gene conversion, we sequenced introns 2, 3, and 8 of all three Class I genes (total>15.0 kb) for five non-human primates--four great apes and one Old World Monkey (OWM)--and compared them with those of humans. The phylogenetic analysis shows each intron sequence clusters strongly within each gene, giving no evidence for gene conversion(s). Several lines of evidence indicate that the first split was between ADH1C and the gene that gave rise to ADH1A and ADH1B. We also analyzed cDNA sequences of the three genes that have been previously reported in mouse and Catarrhines (OWMs, chimpanzee, and humans) and found that the synonymous and non-synonymous substitution (dN/dS) ratios in all pairs are less than 1 representing purifying selection. This suggests that purifying selection is more important than gene conversion(s) in maintaining the overall sequence similarity among the Class I genes. We speculate that the highly conserved sequences on the three duplicated genes in primates have been achieved essentially by maintaining stability of the hetero-dimer formation that might have been related to dietary adaptation in primate evolution.

Highly abundant genes in the transcriptosome of human and baboon CD34 antigen-positive bone marrow cells

Nonhuman primates are useful large animal model systems for the in vivo study of hematopoietic stem cell biology. To better understand the degree of similarity of the hematopoietic systems between humans and baboons, and to explore the relevance of such studies in nonhuman primates to humans, this study was designed to compare the global gene expression profile of bone marrow CD34(+) cells isolated from these 2 species. Human complementary DNA (cDNA) filter arrays containing 25 920 human cDNAs were surveyed for this purpose. The expression pattern and relative gene abundance of the 2 RNA sources were similar, with a correlation coefficient of 0.87. A total of 15 970 of these cDNAs were expressed in human CD34(+) cells, of which the majority (96%) varied less than 3-fold in their relative level of expression between human and baboon. Reverse transcriptase-polymerase chain reaction analysis of selected genes confirmed that expression was comparable between the 2 species. No species-restricted transcripts have been identified, further reinforcing the high degree of similarity between the 2 populations. A subset of 1554 cDNAs, which are expressed at levels 100-fold and greater than background, is described, which includes 959 expressed sequence tags and uncharacterized cDNAs, and 595 named genes, including many that are clearly involved in hematopoiesis. The cDNAs reported here represent a selection of some of the most highly abundant genes in hematopoietic cells and provide a starting point to develop a profile of the transcriptosome of CD34(+) cells.

Recommended nomenclature for the vertebrate alcohol dehydrogenase gene family

The alcohol dehydrogenase (ADH) gene family encodes enzymes that metabolize a wide variety of substrates, including ethanol, retinol, other aliphatic alcohols, hydroxysteroids, and lipid peroxidation products. Studies on 19 vertebrate animals have identified ADH orthologs across several species, and this has now led to questions of how best to name ADH proteins and genes. Seven distinct classes of vertebrate ADH encoded by non-orthologous genes have been defined based upon sequence homology as well as unique catalytic properties or gene expression patterns. Each class of vertebrate ADH shares <70% sequence identity with other classes of ADH in the same species. Classes may be further divided into multiple closely related isoenzymes sharing >80% sequence identity such as the case for class I ADH where humans have three class I ADH genes, horses have two, and mice have only one. Presented here is a nomenclature that uses the widely accepted vertebrate ADH class system as its basis. It follows the guidelines of human and mouse gene nomenclature committees, which recommend coordinating names across species boundaries and eliminating Roman numerals and Greek symbols. We recommend that enzyme subunits be referred to by the symbol "ADH" (alcohol dehydrogenase) followed by an Arabic number denoting the class; i.e. ADH1 for class I ADH. For genes we recommend the italicized root symbol "ADH" for human and "Adh" for mouse, followed by the appropriate Arabic number for the class; i.e. ADH1 or Adh1 for class I ADH genes. For organisms where multiple species-specific isoenzymes exist within a class, we recommend adding a capital letter after the Arabic number; i.e. ADH1A, ADH1B, and ADH1C for human alpha, beta, and gamma class I ADHs, respectively. This nomenclature will accommodate newly discovered members of the vertebrate ADH family, and will facilitate functional and evolutionary studies.

A dynein heavy chain homologue gene in Hexamita inflata

A gene encoding an unusually small dynein heavy chain homologue, hDYHH, was cloned from the genome of a free-living diplomonad, Hexamita inflata (Hi). The open reading frame (ORF) of hDYHH is 867bp and encodes a polypeptide of 289 amino acids (aa), hDYHH. hDYHH is homologous to the region around the third P-loop ATP-binding site of several dynein heavy chain polypeptides that are around 4000aa. Northern blot analysis showed that hDYHH is expressed in vivo and that the mRNA length (approximately 1.8kb) is consistent with the gene length (1.67kb). Southern blot analysis indicated that there are hDYHH homologues within the Hi genome, possibly including a longer dynein heavy chain gene. An hDYHH homologue was also identified in Hexamita pusilla (Hp). hDYHH is the first full-length protein-encoding gene cloned from Hexamita.

Gene expression in a young multigene family: tissue-specific differences in the expression of the human alcohol dehydrogenase genes ADH1, ADH2, and ADH3

Three human alcohol dehydrogenase genes, ADH1, ADH2, and ADH3, were formed by tandem duplications and have diverged in their tissue-specific and developmental expression. Their proximal promoters remain 80-84% identical in sequence, approximately the same degree of identity as at synonymous sites in the coding regions of these three genes. To understand the evolution of tissue specificity, gene expression must be studied in many different cells and tissues. A systematic comparison of their promoters reveals the effects of subtle sequence differences on the binding of nuclear proteins to their cis-acting elements. There are differences in the affinity with which some proteins are bound to altered sites including C/EBP sites, USF/MLTF sites, and the G3T site (which binds Sp1). There are also differences in the sites that are occupied, e.g. CTF/NFI-related sites. These sequence differences are reflected in differences in gene expression in three cell lines. In H4IIE-C3 hepatoma cells, the ADH1 promoter was more active than the ADH2 promoter, and the ADH3 promoter was nearly nonfunctional. In HeLa cells, both ADH1 and ADH2 promoters directed expression; again the ADH3 promoter was extremely weak. None of the three promoters had much activity in CV-1 cells. Coexpression of C/EBP alpha greatly stimulated expression of the ADH1 promoter in HeLa cells and in CV-1 cells, but only weakly stimulated expression in H4IIE-C3 cells. The stimulation of the ADH1 promoter by C/EBP alpha was comparable to that of ADH2, despite the weaker binding to the C/EBP sites that flank the TATA box in ADH1. The ADH3 promoter was not greatly stimulated by C/EBP alpha, despite good binding of C/EBP alpha. These results demonstrate that small differences in the cis-acting elements affect affinity of binding by transcription factors and the pattern of gene expression.

Linking of isozyme and class variability patterns in the emergence of novel alcohol dehydrogenase functions. Characterization of isozymes in Uromastix hardwickii

The nature of the isozyme differences in the class-I alcohol dehydrogenase structure from the lizard, Uromastix hardwickii, was determined and related to those in the human and horse enzymes, for which isozyme structures have also been established. The Uromastix isozymes differ much (at a total of 72 positions, 19%) but, in spite of this, have similar properties and were not obtained resolved. Their structures were analyzed in mixture, and the two sub-sets of peptides obtained could be distinguished by evaluation of the recovery ratios within the peptide pairs. The isozymes have class-I activities, with an ethanol dehydrogenase activity of 0.6 U/mg and no formaldehyde dehydrogenase activity, have typical class-I structures, and are composed of N-terminally acetylated 375-residue subunits (a and b). Importantly, variability patterns between the isozymes are reminiscent of those both in the other two lines with isozymes (primates and horse) and in the class distinctions of the enzyme. Hence, the variability pattern since the distant stage of class-I emergence is also visible within the more recent isozyme divergence, illustrating a continuity in the evolution of isozymes to classes (and then to enzymes). The pattern also links the different levels of multiplicity and may suggest an acceptability in common to duplications and mutations, compatible with the emergence of novel functions.

A Giardia duodenalis gene encoding a protein with multiple repeats of a toxin homologue

A gene, CRP136, from the intestinal protozoan parasite Giardia duodenalis, expressed at a high level in a metronidazole-resistant line, encodes a 136 kDa protein with 23 copies of a 40 amino acid repeat. The protein is cysteine-rich and has the typical membrane-spanning region and CXXC amino acid motifs of a family of Giardia cysteine-rich surface proteins (CRSPs). The repeat unit in CRP136, shares 57% homology with the gene encoding the precursor of the sarafotoxins, a group of snake toxins from the burrowing adder known to cause symptoms similar to those of humans acutely infected with Giardia. The sarafotoxins are low molecular weight sulphydryl cross-linked peptides which are proteolytically cleaved from a precursor polyprotein. CRP136 has homology over the entire length of the sarafotoxin precursor, and the repeats are of the same length. Thus CRP136 represents the first evidence for a potential Giardia toxin. The genomic copy number of CRP136 appears to be the same in both the parent and drug-resistant lines and expression of this gene, and at least one other, is associated with a conserved partial duplication, but not amplification, of one chromosome.

Molecular cloning of the human and monkey sperm surface protein PH-20

The guinea pig sperm surface protein PH-20 has an essential function in sperm adhesion to the zona pellucida of guinea pig eggs. Fully effective contraception has been achieved by immunizing either male or female guinea pigs with purified guinea pig PH-20. Here we report the isolation of human and cynomolgus monkey PH-20 cDNAs as a key step toward testing the function of primate PH-20 and the contraceptive efficacy of PH-20 immunization in primates. The deduced amino acid sequence of human PH-20 has 509 residues and is 59% identical with guinea pig PH-20, suggesting they may have a conserved function and immunogenicity. Southern blots show that there is a single PH-20 gene in the human genome and Northern blots of human testis poly(A)+ RNA show a 2.4-kb message. Northern blots of tissues other than testis are negative for PH-20, indicating that human PH-20 is testis-specific.

Amino acid substitution at the Adh locus of Drosophila is facilitated by small population size

The number of amino acid replacement substitutions and that of synonymous substitutions are examined by using DNA sequences of the Adh locus of Drosophila. The ratio of replacement to synonymous substitutions is higher in sequence comparisons between species than in polymorphisms within species. The ratio for the between-species comparisons is highest in the Hawaiian group and lowest in the obscura group. These observations suggest that amino acid substitutions are facilitated by small population size. The result is in accord with the nearly neutral theory of molecular evolution.

Isozyme developments in mammalian class-I alcohol dehydrogenase. cDNA cloning, functional correlations, and lack of evidence for genetic isozymes in rabbit

Isozyme patterns differ widely among the classical type (class I) of mammalian alcohol dehydrogenases. For the rabbit enzyme, the possibility of isozymes has been reported but structural evidence is lacking. This system was now studied at both the mRNA/cDNA and protein levels. Ten cDNA clones, coding for class-I alcohol dehydrogenase, were isolated from a rabbit liver cDNA library using a human DNA fragment as probe. The cDNA spanned 1296 bp, including the entire coding region. All clones coded for the same polypeptide and Northern blots identified a single mRNA corresponding to about 1.5 kb. Comparison of two protein forms (CC and BC) by HPLC peptide fingerprinting and structural analysis revealed peptide segments identical in amino acid sequence. Consequently, direct protein analyses and Northern blots show the presence of only one primary translation product. The data suggest that lagomorphic alcohol dehydrogenase, like the rodent enzyme, is not as isozyme rich as it may appear superficially, and that secondary modifications contribute substantially to mammalian alcohol dehydrogenase multiplicity. The active center of the rabbit enzyme suggests similarities to the horse S, human gamma, and rat enzyme structures, compatible with a steroid dehydrogenase activity shown experimentally. Typical class-I properties were established by direct analysis and confirmed by structural properties (Km for cyclohexanol 0.8-1.1 mM, for ethanol 1.6-1.9 mM). The isozyme versus species differences mark the variability of class-I alcohol dehydrogenase versus class III and suggest a parallelism between rapid mutational differences and frequent duplicatory events.

Progressive sequence alignment and molecular evolution of the Zn-containing alcohol dehydrogenase family

Sequences of 47 members of the Zn-containing alcohol dehydrogenase (ADH) family were aligned progressively, and an evolutionary tree with detailed branch order and branch lengths was produced. The alignment shows that only 9 amino acid residues (of 374 in the horse liver ADH sequence) are conserved in this family; these include eight Gly and one Val with structural roles. Three residues that bind the catalytic Zn and modulate its electrostatic environment are conserved in 45 members. Asp 223, which determines specificity for NAD, is found in all but the two NADP-dependent enzymes, which have Gly or Ala. Ser or Thr 48, which makes a hydrogen bond to the substrate, is present in 46 members. The four Cys ligands for the structural zinc are conserved except in zeta-crystallin, the sorbitol dehydrogenases, and two bacterial enzymes. Analysis of the evolutionary tree gives estimates of the times of divergence for different animal ADHs. The human class II (pi) and class III (chi) ADHs probably diverged about 630 million years ago, and the newly identified human ADH6 appeared about 520 million years ago, implying that these classes of enzymes may exist or have existed in all vertebrates. The human class I ADH isoenzymes (alpha, beta, and gamma) diverged about 80 million years ago, suggesting that these isoenzymes may exist or have existed in all primates. Analysis of branch lengths shows that these plant ADHs are more conserved than the animal ones and that class III ADHs are more conserved than class I ADHs. The rate of acceptance of point mutations (PAM units) shows that selection pressure has existed for ADHs, implying that these enzymes play definite metabolic roles.

Biochemical genetics of alcohol dehydrogenase isozymes in the gray short-tailed opossum (Monodelphis domestica)

Polyacrylamide gel-isoelectric focusing (PAGE-IEF) methods were used to examine the multiplicity, tissue distribution, and biochemical genetics of alcohol dehydrogenase (ADH) isozymes among gray short-tailed opossums (Monodelphis domestica). Seven ADH isozymes were resolved and distinguished on the basis of their isoelectric points, tissue distributions, and substrate and inhibitor specificities. ADH1 and ADH2 exhibited Class I properties and were observed in liver (and intestine) extracts. ADH3, ADH4, and ADH5 showed "high-Km" (possibly Class IV) properties, with ADH3 and ADH4 exhibiting high activity in cornea, ear, stomach, and esophagus extracts. ADH6 and ADH7 exhibited Class III properties, including activities as formaldehyde dehydrogenases, with each showing different tissue distribution characteristics; ADH6 was widely distributed, and ADH7 was restricted to prostate extracts. An additional form of formaldehyde dehydrogenase (FDH) was observed, which was inactive with hexenol and ethanol as substrates. Isoelectric point variants were observed for ADH3 (three forms) and for ADH4 (two forms), and the inheritance of ADH3 was studied in 15 families of M. domestica. The data were consistent with codominant inheritance of two alleles (ADH3*A and ADH3*B) at a single autosomal locus (designated ADH3) and with a model involving a dimeric ADH isozyme: ADH3 (gamma 2 isozyme, forming three dimers designated gamma 1(2), gamma 1 gamma 2, and gamma 2(2) in heterozygous individuals).

Induction of the major integral membrane protein of mouse liver peroxisomes by peroxisome proliferators

Peroxisome proliferators are known to increase the volume of the peroxisomal compartment in rodent liver. We have examined the induction of the major integral membrane protein of mouse liver peroxisomes (PMP68) by a number of these agents, and compared this with their effect on the peroxisomal bifunctional protein (PBP), an enzyme of the beta-oxidation pathway which is located in the peroxisome matrix. Dietary clofibrate, di-2-(ethylhexyl)phthalate and Wy-14,643, three structurally unrelated proliferators, all increased the mRNA and protein content of PMP68 approx. 2-fold, whereas PBP was induced 8-13-fold. The kinetics and sequence of induction of PMP68 and PBP following a single dose of Wy-14,643 were compared and shown to be similar, and the effects were reversible. Another proliferator, BM 15766, caused maximal induction of PMP68 but only a low induction of PBP; further PBP induction was achieved by the administration of BM 15766 in combination with Wy-14,643. Similarly, BM 15766 and Wy-14,643 increased transcription of the PMP68 gene in vitro, whereas PBP gene transcription was increased by Wy-14,643 but not by BM 15766. Thus peroxisome proliferators enhance the expression of the genes for both the membrane protein PMP68 and the matrix protein PBP, but the regulation of this expression appears to be mediated by different mechanisms.

Multiplication of the class I alcohol dehydrogenase locus in mammalian evolution

Chromosomal DNA samples derived from various primates and other mammals (horse, sheep, rabbit, and mouse) were digested with restriction endonuclease and hybridized with a probe of the sixth exon of the human ADH gene, which is highly conserved in the class I alcohol dehydrogenase of these mammalian species. The copy number of the class I ADH gene in each species was estimated from the number of hybridized bands. Primate DNA samples showed three distinct bands in the blots of PstI digest and DraI digest. Moreover, most of the bands from primate DNA showed a similarity in size so as to allow us to assign the ADH1, ADH2, and ADH3 homologues in each species. In contrast, mouse has only one gene, and rabbit, sheep, and horse seem to have only two genes, for the class I ADH, which showed divergent hybridization bands. These results are consistent with the view that the human class I ADH gene cluster has been generated through gene multiplication events which occurred before the Catarrhini branch point in the course of primate evolution.