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.     

Evolution of major histocompatibility complex class I genes in Cetartiodactyls

Previous studies of cattle MHC have suggested the presence of at least four classical class I loci. Analysis of haplotypes showed that any combination of one, two or three genes may be expressed, although no gene is expressed consistently. The aim of this study was to examine the evolutionary relationships among these genes and to study their phylogenetic history in Cetartiodactyl species, including cattle and their close relatives. A secondary aim was to determine whether recombination had occurred between any of the genes. MHC class I data sets were generated from published sequences or by polymerase chain reaction from cDNA. Phylogenetic analysis revealed that MHC class I sequences from Cetartiodactyl species closely related to cattle were distributed among the main cattle gene "groups", while those from more distantly related species were either scattered (sheep, deer) or clustered in a species-specific manner (sitatunga, giraffe). A comparison between gene and species trees showed a poor match, indicating that divergence of the MHC sequences had occurred independently from that of the hosts from which they were obtained. We also found two clear instances of interlocus recombination among the cattle MHC sequences. Finally, positive natural selection was documented at positions throughout the alpha 1 and 2 domains, primarily on those amino acids directly involved in peptide binding, although two positions in the alpha 3 domain, a region generally conserved in other species, were also shown to be undergoing adaptive evolution.     

Group I introns within the nuclear-encoded small-subunit rRNA gene of three green algae.

Four group I introns from the nuclear-encoded (18S) rRNA genes of three chlorophycean green algae are described; two are in Dunaliella parva, and one each is in D. salina and Characium saccatum. The introns within the gene in the latter two organisms are located at the sites equivalent to the 5' and 3' introns of D. parva, respectively. All four introns lack open reading frames and are relatively small, 381-447 bp. Both primary- and secondary-structural features place these introns within subgroup IC1 described by Michel and Westhof. Phylogenetic relationships of the three intron-containing taxa and their relatives, as inferred from comparisons of 18S rDNA sequences, suggest that inheritance of the introns along with the gene can account for their present distribution. The discovery of these four introns, in addition to two others known to exist in other chlorophycean green algae, suggests that group I introns within the 18S rRNA gene may be relatively common in the green algae.     

Sequence variability at three MHC loci of finless porpoises (<Emphasis Type="Italic">Neophocaena phocaenoides</Emphasis>)

Major histocompatibility complex (MHC) class II DQB and DRA genes and class I gene of finless porpoises (Neophocaena phocaenoides) were investigated by single-strand conformation polymorphism and sequence analysis. The DRA, DQB, and MHC-I loci each contained 5, 14, and 34 unique sequences, respectively, and considerable sequence variation was found at the MHC-I and DQB loci. Gene duplication was manifested as three to five distinct sequences at each of the DQB and MHC-I loci from some individuals, and these sequences at each of the two loci separately clustered into four groups (cluster A, B, C, and D) based on the phylogenetic trees. Phylogenetic reconstruction revealed a trans-species pattern of evolution. Relatively high rates of non-synonymous (dN) vs synonymous (dS) substitution in the peptide-binding region (PBR) suggested balancing selection for maintaining polymorphisms at the MHC-I and DQB loci. In contrast, one single locus with little sequence variation was detected in the DRA gene, and no non-synonymous substitutions in the PBR indicated no balancing selection on this gene.     

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.     

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.     

The R1 resistance gene cluster contains three groups of independently evolving, type I R1 homologues and shows substantial structural variation among haplotypes of Solanum demissum

Cultivated and wild potatoes contain a major disease-resistance cluster on the short arm of chromosome V, including the R1 resistance (R) gene against potato late blight. To explore the functional and evolutionary significance of clustering in the generation of novel disease-resistance genes, we constructed three approximately 1 Mb physical maps in the R1 gene region, one for each of the three genomes (haplotypes) of allohexaploid Solanum demissum, the wild potato progenitor of the R1 locus. Totals of 691, 919 and 559 kb were sequenced for each haplotype, and three distinct resistance-gene families were identified, one homologous to the potato R1 gene and two others homologous to either the Prf or the Bs4 R-gene of tomato. The regions with R1 homologues are highly divergent among the three haplotypes, in contrast to the conserved flanking non-resistance gene regions. The R1 locus shows dramatic variation in overall length and R1 homologue number among the three haplotypes. Sequence comparisons of the R1 homologues show that they form three distinct clades in a distance tree. Frequent sequence exchanges were detected among R1 homologues within each clade, but not among those in different clades. These frequent sequence exchanges homogenized the intron sequences of homologues within each clade, but did not homogenize the coding sequences. Our results suggest that the R1 homologues represent three independent groups of fast-evolving type I resistance genes, characterized by chimeric structures resulting from frequent sequence exchanges among group members. Such genes were first identified among clustered RGC2 genes in lettuce, where they were distinguished from slow-evolving type II R-genes. Our findings at the R1 locus in S. demissum may indicate that a common or similar mechanism underlies the previously reported differentiation of type I and type II R-genes and the differentiation of type I R-genes into distinct groups, identified here.     

Independent evolution of functional MHC class II DRB genes in New World bat species

Genes of the major histocompatibility complex (MHC) play a pivotal role in the vertebrate immune system and are attractive markers for functional, fitness-related, genetic variation. Although bats (Chiroptera) represent the second largest mammalian order and are prone to various emerging infectious diseases, little is known about MHC evolution in bats. In the present study, we examined expressed MHC class II DRB sequences (exons 1 to 4) of New World bat species, Saccopteryx bilineata, Carollia perspicillata, Noctilio albiventris and Noctilio leporinus (only exon 2). We found a wide range of copy number variation of DRB loci with one locus detected in the genus Noctilio and up to ten functional loci observed in S. bilineata. Sequence variation between alleles of the same taxa was high with evidence for positive selection. We found statistical support for recombination or gene conversion events among sequences within the same but not between bat species. Phylogenetic relationships among DRB alleles provided strong evidence for independent evolution of the functional MHC class II DRB genes in the three investigated species, either by recent gene duplication, or homogenization of duplicated loci by frequent gene conversion events. Phylogenetic analysis of all available chiropteran DRB exon 2 sequences confirmed their monophyletic origin within families, but revealed a possible trans-species mode of evolution pattern in congeneric bat species, e.g. within the genera Noctilio and Myotis. This is the first study investigating phylogenetic relationships of MHC genes within bats and therefore contributes to a better understanding of MHC evolution in one of the most dominant mammalian order.     

Sequence diversification of 45S rRNA ITS, trnH-psbA spacer, and matK genic regions in several Allium species

Allium is a very diverse genus with over 600 species distributed worldwide. Haplotype analyses of 45S rRNA ITS, trnH-psbA spacer, and matK gene sequences in 9 Allium species were carried out, subsequent to which phylogenetic relations of the nine species were also analyzed. Of the three genes, the nuclear 45S rRNA ITS sequences showed the highest variation with one haplotype in each species. The other two chloroplast genes revealed that more than one haplotype was present in each species, and each haplotype was present in several of the species. In the matK gene, EcoRI restriction revealed heteroplasmy in which the functional gene retains the EcoRI recognition site while the nonfunctional, pseudogene does not. Phylogenetic patterns were not consistent among the haplotypes of the 45 rRNA ITS, trnH-psbA spacer, and matK genic regions. This phylogenetic incongruency might be due to the presence of multiple haplotypes in each of the chloroplast genes. However, the inconsistency of the phylogenetic relationships, based on the 45S rRNA ITS sequences makes a strong case for further analysis.     

Analysis of MHC class I genes across horse MHC haplotypes

The genomic sequences of 15 horse major histocompatibility complex (MHC) class I genes and a collection of MHC class I homozygous horses of five different haplotypes were used to investigate the genomic structure and polymorphism of the equine MHC. A combination of conserved and locus-specific primers was used to amplify horse MHC class I genes with classical and nonclassical characteristics. Multiple clones from each haplotype identified three to five classical sequences per homozygous animal and two to three nonclassical sequences. Phylogenetic analysis was applied to these sequences, and groups were identified which appear to be allelic series, but some sequences were left ungrouped. Sequences determined from MHC class I heterozygous horses and previously described MHC class I sequences were then added, representing a total of ten horse MHC haplotypes. These results were consistent with those obtained from the MHC homozygous horses alone, and 30 classical sequences were assigned to four previously confirmed loci and three new provisional loci. The nonclassical genes had few alleles and the classical genes had higher levels of allelic polymorphism. Alleles for two classical loci with the expected pattern of polymorphism were found in the majority of haplotypes tested, but alleles at two other commonly detected loci had more variation outside of the hypervariable region than within. Our data indicate that the equine major histocompatibility complex is characterized by variation in the complement of class I genes expressed in different haplotypes in addition to the expected allelic polymorphism within loci.     

Shift of phylogenic position in megalocytiviruses based on three different genes

Major capsid protein (MCP), the adenosine triphosphatase (ATPase), and the PstI fragment genes from five Japanese and three Korean megalocytivirus isolates were sequenced and phylogenetically analyzed with known megalocytiviruses. Phylogenetic trees formed three major clusters (M1, M2, and M3 or P1, P2, and P3), and genogroup I was divided into two minor clusters (M1a/M1b and P1a/P1b) using three target genes. Sequence identity was >97% within each cluster, except cluster II of the PstI fragment (>94% of sequence identity). Interestingly, different genotyping patterns were observed for the same isolates depending on the gene analyzed. The JPN-YelTail and JPN-BfTuna isolates located in the minor M1a cluster, based on MCP and ATPase nucleotide sequences, appeared in the minor P1b cluster based on the PstI fragment, suggesting a shift of phylogenic position in megalocytiviruses. Further study will be conducted to compare the viral antigenicity and pathogenicity between the two isolates showing the shift of phylogenic position and the other isolates clustered within genogroup I.     

Proceedings of the SMBE Tri-National Young Investigators' Workshop 2005. MHC Class I genes in the Tuatara (Sphenodon spp.): evolution of the MHC in an ancient reptilian order

The major histocompatibility complex (MHC) is an extremely dynamic region of the genome, characterized by high polymorphism and frequent gene duplications and rearrangements. This has resulted in considerable differences in MHC organization and evolution among vertebrate lineages, particularly between birds and mammals. As nonavian reptiles are ancestral to both mammals and birds, they occupy an important phylogenetic position for understanding these differences. However, little is known about reptile MHC genes. To address this, we have characterized MHC class I sequences from the tuatara (Sphenodon spp.), the last survivor of an ancient order of reptiles, Sphenodontia. We isolated two different class I cDNA sequences, which share 93% sequence similarity with each other but are highly divergent from other vertebrate MHC genes. Southern blotting and polymerase chain reaction amplification of class I sequences from seven adult tuatara plus a family group indicate that these sequences represent at least two to three loci. Preliminary analysis of variation among individuals from an island population of tuatara indicates that these loci are highly polymorphic. Maximum likelihood analysis of reptile MHC class I sequences indicates that gene duplication has occurred within reptilian orders. However, the evolutionary relationships among sequences from different reptilian orders cannot be resolved, reflecting the antiquity of the major reptile lineages.     

Phylogenetic and evolutionary analysis of NBS-encoding genes in Rosaceae fruit crops

Phylogenetic relationships of the nucleotide binding site (NBS)-encoding resistance gene homologues (RGHs) among 12 species in five genera of Rosaceae fruit crops were evaluated. A total of 228 Rosaceous RGHs were deeply separated into two distinct clades, designated as TIR (sequences within this clade containing a Toll Interleukin-1 Receptor domain) and NonTIR (sequences lacking a TIR domain). Most Rosaceous RGH genes were phylogenetically distinct from Arabidopsis, Rice or Pine genes, except for a few Rosaceous members which grouped closely with Arabidopsis genes. Within Rosaceae, sequences from multiple species were often phylogenetically clustered together, forming heterogenous groups, however, apple- and chestnut rose-specific groups really exist. Gene duplication followed by sequence divergence were proposed as the mode for the evolution of a large number of distantly or closely related RGH genes in Rosaceae, and this mode may play a role in the generation of new resistance specificity. Positively selected sites within NBS-coding region were detected and thus nucleotide variation within NBS domain may function in determining disease resistance specificity. This study also discusses the synteny of a genomic region that encompass powdery mildew resistance locus among Malus, Prunus and Rosa, which may have potential use for fruit tree disease breeding and important gene cloning.     

Molecular evolution of tephritid fruit flies in the genus Bactrocera based on the cytochrome oxidase I gene

Fruit flies of the genus Bactrocera (Diptera: Tephritidae) are one of the major economically important insects in Asia and Australia. Little attention has been given to analyses of molecular phylogenetic relationships among Bactrocera subgenera. By using mitochondrial cytochrome oxidase I gene (COI) sequences, the phylogenetic relationships among four subgenera, Asiadacus, Bactrocera, Hemigymnodacus, and Zeugodacus, were investigated. Nucleotide diversity within subgenera ranged from 11.7 to 12.4%, and the net divergence among subgenera ranged from 11.2 to 15.7%. Phylogenetic trees calculated from both maximum parsimony and neighbor-joining phylogenetic analysis methods were highly congruent in terms of tree topologies. Phylogenetic analysis of mitochondrial COI sequences suggests that tephritid fruit fly species, which attack cucurbit plants, that is, Asiadacus, Hemigymnodacus and Zeugodacus, were more closely related to each other than to fruit fly species of the subgenus Bactrocera, which attack plants of numerous families. Our data supports previous classification of Bactrocera based on morphological characters. However, the phylogenetic tree showed the polyphyletic of fruit flies in subgenus Zeugodacus. Possible causes of speciation among fruit flies species in this genus were also discussed.     

Evolution and nomenclature of the zona pellucida gene family

Three subfamilies of genes are acknowledged within the zona pellucida (ZP) gene family. At present, these subfamilies each have two names that are used interchangeably: ZPA or ZP2, ZPB or ZP1, and ZPC or ZP3. The ZPA genes encode the longest protein sequences and the ZPC genes the shortest. Recently, several sequences, which have no clear relationship to the three subfamilies, have been identified. These sequences include two paralogous ZP genes from Xenopus laevis and a single gene from the fish Oryzias latipes. We have conducted extensive phylogenetic analyses of the known ZP genes. As well as establishing the evolutionary relationships among these genes, the analyses make it clear that the dual nomenclature system is no longer feasible, because major paralogous groups are present in the ZPB (ZP1) family of genes of amniotes. We propose a unified system of nomenclature for the ZP gene family that removes the existing ambiguities.     

Caenorhabditis elegans contains genes encoding two new members of the Zn-containing alcohol dehydrogenase family

We have characterized two cDNA clones from the nematode Caenorhabditis elegans that display similarity to the alcohol dehydrogenase (ADH) gene family. The nucleotide sequences of these cDNAs predict that they encode Zn-containing long-chain ADH enzymes. Phylogenetic analysis suggests that one is most similar to dimeric class III ADHs found in diverse taxa; the other is most similar to the tetrameric forms of ADH previously described only in fungi. Correspondence to: J.J. Collins