Origins and divergence times of mammalian class II MHC gene clusters

The class I and II major histocompatibility complex (MHC) genes are apparently subject to evolution by a birth-and-death process. The rate of gene turnover is much slower in the latter genes than in the former. In placental mammals, the class II region can be subdivided into different orthologous subregions or gene clusters (DR, DQ, DO, and DN), but the origins and evolutionary relationships of these gene clusters are not well established. Here we report the results of our study of the times of origin and evolutionary relationships of these gene clusters in mammals. Our analysis suggests that both class II alpha-chain and beta-chain gene clusters are shared by placental mammals and marsupials, but the gene clusters from nonmammalian species are paralogous to mammalian gene clusters. We estimated the times of divergence between gene clusters in placental mammals using the linearized tree and distance regression methods. Our results indicate that most gene clusters originated 170-200 million years (MY) ago, but that DO beta-chain genes diverged from the other beta-chain gene clusters approximately 210-260 MY ago. The phylogenetic trees for the alpha- and beta-chain genes were not congruent, suggesting that the evolutionary history of the class II gene clusters is more complex than previously thought.     

Maximum-likelihood approach for gene family evolution under functional divergence

According to the observed alignment pattern (i.e., amino acid configuration), we studied two basic types of functional divergence of a protein family. Type I functional divergence after gene duplication results in altered functional constraints (i.e., different evolutionary rate) between duplicate genes, whereas type II results in no altered functional constraints but radical change in amino acid property between them (e.g., charge, hydrophobicity, etc.). Two statistical approaches, i.e., the subtree likelihood and the whole-tree likelihood, were developed for estimating the coefficients of (type I or type II) functional divergence. Numerical algorithms for obtaining maximum-likelihood estimates are also provided. Moreover, a posterior-based site-specific profile is implemented to predict critical amino acid residues that are responsible for type I and/or type II functional divergence after gene duplication. We compared the current likelihood with a fast method developed previously by examples; both show similar results. For handling altered functional constraints (type I functional divergence) in the large gene family with many member genes (clusters), which appears to be a normal case in postgenomics, the subtree likelihood provides a solution that is computationally feasible and robust against the uncertainty of the phylogeny. The cost of this feasibility is the approximation when frequencies of amino acids are very skewed. The potential bias and correction are discussed.     

Rapid expansion of killer cell immunoglobulin-like receptor genes in primates and their coevolution with MHC Class I genes

The gene family of killer cell immunoglobulin-like receptors (KIRs) in primates provides the first line of defense against virus infection and tumor transformation. Interacting with MHC class I molecules, KIRs can regulate the cytotoxic activity of natural killer (NK) cells and distinguish the tumor and virus infected cells from normal body cells. Phylogenetic analysis and comparison of domain structures identified three major groups of KIR genes (group I, II, and III genes). These groups of KIR genes, generated by a series of gene duplications, have acquired different MHC-binding specificity. Inference of ancestral KIR sequences suggested that the functional divergence of group I genes from group II genes occurred by positive selection at the MHC-binding sites after duplication. Our evolutionary study has shown that group I genes diverged from group II genes about 17 million years ago (Mya) apparently after separation of hominoids from Old World (OW) monkeys. Around the same time, gene duplication generating the class I MHC-C locus appears to have occurred. These findings suggest that KIR and MHC class I genes have coevolved as an interacting system. The KIR gene family has experienced a rapid expansion in primate species. The rate of expansion of this gene family seems to be one of the highest among all hominoid gene families. The KIR gene family is also subject to birth-and-death evolution.     

载脂蛋白多基因家族分子进化的研究

与脂质运输有关的载脂蛋白基因构成一个复杂的多基因家族。为探讨这种演化时间长的基因家族的进化规律,本文首先建立了一种在非均衡进化速率条件下计算系统发生树中任意分支长度的简易方法,并可在此基础上算出无根分支系统树中分歧年代的期望值。进一步对本文科１０个种属共２６种载脂蛋白的系统演作作了实际分析,结果提示：①ＡｐｏＡ－Ｉ＇ＡｐｏＡ－ＩＶ,ＡｐｏＥ及ＡｐｏＡ－ＩＩ的共同祖先可能在奥陶纪水生脊椎动物中就已存     

Concerted evolution of the immunoglobulin VH gene family

With the aim of understanding the concerted evolution of the immunoglobulin VH multigene family, a phylogenetic tree for the DNA sequences of 16 mouse and five human germ line genes was constructed. This tree indicates that all genes in this family have undergone substantial evolutionary divergence. The most closely related genes so far identified in the mouse genome seem to have diverged about 6 million years (MY) ago, whereas the most distantly related genes diverged about 300 MY ago. This suggests that gene duplication caused by unequal crossing-over or gene conversion occurs very slowly in this gene family. The rate of occurrence of gene duplication in the VH gene family has been estimated to be 5 x 10(-7) per gene per year, which seems to be at least about 100 times lower than that for the rRNA gene family. This low rate of concerted evolution in the VH gene family helps retain intergenic genetic variability that in turn contributes to antibody diversity. Because of accumulation of destructive mutations, however, about one-third of the mouse and human VH genes seem to have become nonfunctional. Many of these pseudogenes have apparently originated recently, but some of them seem to have existed in the genome for more than 10 MY. The rate of nucleotide substitution for the complementarity-determining regions (CDRs) is as high as that of pseudogenes. This suggests that there is virtually no purifying selection operating in the CDRs and that germ line mutations are effectively used for generating antibody diversity.      

Mitochondrial DNA phylogeography of the red fox (Vulpes vulpes) in northern Japan

Mitochondrial DNA variation in the cytochrome b (cyt b) gene and the control region was examined in the red fox Vulpes vulpes from Japan, with special focus on the population divergence between Hokkaido and northern Honshu. Resultant haplotypes from Hokkaido were subdivided into two distinct groups (I and II), with an average genetic distance of 0.027 for cyt b. Divergence time is roughly estimated to be 1-2 million years ago, given that the conventional divergence rate of the mammalian cyt b gene is 2% per million years. Notably, Group II was only found in Hokkaido, whereas Group I comprised haplotypes from Honshu, Kyushu (Japan), eastern Russia, and Europe, as indicated by a comparison of our own data to the literature. On the other hand, judging from constructed trees, Group I haplotypes from Hokkaido appeared to differ from those from other parts of Japan, i.e., Honshu and Kyushu. This implies that Blakiston's Line, which demarcates the boundary between Hokkaido and Honshu, has been an effective barrier and has allowed the structuring of genetic variation in maternal lineages. Thus, these results suggest that the Hokkaido population, which is sometimes referred to as the distinct subspecies V. v. schrencki, has its own genetic background with multiple migration events and differs from the parapatric subspecies V. v. japonica found in Honshu and Kyushu.     

The phylogenetics of mycotoxin and sclerotium production in Aspergillus flavus and Aspergillus oryzae

Aspergillus flavus is a common filamentous fungus that produces aflatoxins and presents a major threat to agriculture and human health. Previous phylogenetic studies of A. flavus have shown that it consists of two subgroups, called groups I and II, and morphological studies indicated that it consists of two morphological groups based on sclerotium size, called "S" and "L." The industrially important non-aflatoxin-producing fungus A. oryzae is nested within group I. Three different gene regions, including part of a gene involved in aflatoxin biosynthesis (omt12), were sequenced in 33 S and L strains of A. flavus collected from various regions around the world, along with three isolates of A. oryzae and two isolates of A. parasiticus that were used as outgroups. The production of B and G aflatoxins and cyclopiazonic acid was analyzed in the A. flavus isolates, and each isolate was identified as "S" or "L" based on sclerotium size. Phylogenetic analysis of all three genes confirmed the inference that group I and group II represent a deep divergence within A. flavus. Most group I strains produced B aflatoxins to some degree, and none produced G aflatoxins. Four of six group II strains produced both B and G aflatoxins. All group II isolates were of the "S" sclerotium phenotype, whereas group I strains consisted of both "S" and "L" isolates. Based on the omt12 gene region, phylogenetic structure in sclerotium phenotype and aflatoxin production was evident within group I. Some non-aflatoxin-producing isolates of group I had an omt12 allele that was identical to that found in isolates of A. oryzae.     

The Wilhelmine E. Key 2001 Invitational Lecture. Estimation of divergence times for a few mammalian and several primate species

Statistical methods for estimating divergence times by using multiprotein gamma distances are discussed. When a large number of proteins are used, even a small degree of deviation from the molecular clock hypothesis can be detected. In this case, one may use the stem-lineage method for estimating divergence times. However, the estimates obtained by this method are often similar to those obtained by the linearized tree method. Application of these methods to a dataset of 104 proteins from several vertebrate species indicated that the divergence times between humans and mice and between mice and rats are about 96 and 33 million years (MY) ago, respectively. These estimates were obtained by assuming that birds and mammals diverged 310 MY ago. Similarly application of the methods to the protein sequence data from primate species indicated that the human lineage separated from the chimpanzee, gorilla, Old World monkeys, and New World monkeys about 6.0, 7.0, 23.0, and 33.0 MY ago, respectively. In this case the use of two calibration points, that is, the divergence time (13 MY ago) between humans and orangutans and between primates and artiodactyls (90 MY ago) gave essentially the same estimates.     

Analysis of the evolution of angiotensin II type 1 receptor gene in mammals (mouse, rat, bovine and human).

The nucleotide and amino acid sequences for mouse angiotensin II (AII) type 1A and 1B receptors were deduced from their complementary and genomic DNAs. Evolutionary analyses based on the nucleotide sequences of the coding region of AII type 1 receptor genes indicated that the duplication event of the type 1 gene occurred 24 +/- 2 million years ago before the divergence between the rat and mouse but after the divergence between rodents and the human/artiodactyls couple. This conclusion was consistent with the results of genomic Southern blot analyses, which revealed that the mouse and rat possess 2 similar but separate genes, whereas the bovine and human have only a single class gene.     

A Mesophilic,Autotrophic, Ammonia-Oxidizing Archaeon of Thaumarchaeal Group I.1a Cultivated from a Deep Oligotrophic Soil Horizon

Soil nitrification plays an important role in the reduction of soil fertility and in nitrate enrichment of groundwater. Various ammonia-oxidizing archaea (AOA) are considered to be members of the pool of ammonia-oxidizing microorganisms in soil. This study reports the discovery of a chemolithoautotrophic ammonia oxidizer that belongs to a distinct clade of nonmarine thaumarchaeal group I.1a, which is widespread in terrestrial environments. The archaeal strain MY2 was cultivated from a deep oligotrophic soil horizon. The similarity of the 16S rRNA gene sequence of strain MY2 to those of other cultivated group I.1a thaumarchaeota members, i.e., Nitrosopumilus maritimus and “Candidatus Nitrosoarchaeum koreensis,” is 92.9% for both species. Extensive growth assays showed that strain MY2 is chemolithoautotrophic, mesophilic (optimum temperature, 30°C), and neutrophilic (optimum pH, 7 to 7.5). The accumulation of nitrite above 1 mM inhibited ammonia oxidation, while ammonia oxidation itself was not inhibited in the presence of up to 5 mM ammonia. The genome size of strain MY2 was 1.76 Mb, similar to those of N. maritimus and “Ca. Nitrosoarchaeum koreensis,” and the repertoire of genes required for ammonia oxidation and carbon fixation in thaumarchaeal group I.1a was conserved. A high level of representation of conserved orthologous genes for signal transduction and motility in the noncore genome might be implicated in niche adaptation by strain MY2. On the basis of phenotypic, phylogenetic, and genomic characteristics, we propose the name “Candidatus Nitrosotenuis chungbukensis” for the ammonia-oxidizing archaeal strain MY2.     

The complete chloroplast genome of Gracilariopsis lemaneiformis (Rhodophyta) gives new insight into the evolution of family Gracilariaceae

The complete chloroplast genome of Gracilariopsis lemaneiformis was recovered from a Next Generation Sequencing data set. Without quadripartite structure, this chloroplast genome (183,013 bp, 27.40% GC content) contains 202 protein‐coding genes, 34 tRNA genes, 3 rRNA genes, and 1 tmRNA gene. Synteny analysis showed plasmid incorporation regions in chloroplast genomes of three species of family Gracilariaceae and in Grateloupia taiwanensis of family Halymeniaceae. Combined with reported red algal plasmid sequences in nuclear and mitochondrial genomes, we postulated that red algal plasmids may have played an important role in ancient horizontal gene transfer among nuclear, chloroplast, and mitochondrial genomes. Substitution rate analysis showed that purifying selective forces maintaining stability of protein‐coding genes of nine red algal chloroplast genomes over long periods must be strong and that the forces acting on gene groups and single genes of nine red algal chloroplast genomes were similar and consistent. The divergence of Gp. lemaneiformis occurred ~447.98 million years ago (Mya), close to the divergence time of genus Pyropia and Porphyra (443.62 Mya).     

Different patterns of evolution in the centromeric and telomeric regions of group A and B haplotypes of the human killer cell Ig-like receptor locus

The fast evolving human KIR gene family encodes variable lymphocyte receptors specific for polymorphic HLA class I determinants. Nucleotide sequences for 24 representative human KIR haplotypes were determined. With three previously defined haplotypes, this gave a set of 12 group A and 15 group B haplotypes for assessment of KIR variation. The seven gene-content haplotypes are all combinations of four centromeric and two telomeric motifs. 2DL5, 2DS5 and 2DS3 can be present in centromeric and telomeric locations. With one exception, haplotypes having identical gene content differed in their combinations of KIR alleles. Sequence diversity varied between haplotype groups and between centromeric and telomeric halves of the KIR locus. The most variable A haplotype genes are in the telomeric half, whereas the most variable genes characterizing B haplotypes are in the centromeric half. Of the highly polymorphic genes, only the 3DL3 framework gene exhibits a similar diversity when carried by A and B haplotypes. Phylogenetic analysis and divergence time estimates, point to the centromeric gene-content motifs that distinguish A and B haplotypes having emerged ~6 million years ago, contemporaneously with the separation of human and chimpanzee ancestors. In contrast, the telomeric motifs that distinguish A and B haplotypes emerged more recently, ~1.7 million years ago, before the emergence of Homo sapiens. Thus the centromeric and telomeric motifs that typify A and B haplotypes have likely been present throughout human evolution. The results suggest the common ancestor of A and B haplotypes combined a B-like centromeric region with an A-like telomeric region.     

Genotypic and phenotypic diversity of rhizobia isolated from Lathyrus japonicus indigenous to Japan

Sixty-one rhizobial strains from Lathyrus japonicus nodules growing on the seashore in Japan were characterized and compared to two strains from Canada. The PCR-based method was used to identify test strains with novel taxonomic markers that were designed to discriminate between all known Lathyrus rhizobia. Three genomic groups (I, II, and III) were finally identified using RAPD, RFLP, and phylogenetic analyses. Strains in genomic group I (related to Rhizobium leguminosarum) were divided into two subgroups (Ia and Ib) and subgroup Ia was related to biovar viciae. Strains in subgroup Ib, which were all isolated from Japanese sea pea, belonged to a distinct group from other rhizobial groups in the recA phylogeny and PCR-based grouping, and were more tolerant to salt than the isolate from an inland legume. Test strains in genomic groups II and III belonged to a single clade with the reference strains of R. pisi, R. etli, and R. phaseoli in the 16S rRNA phylogeny. The PCR-based method and phylogenetic analysis of recA revealed that genomic group II was related to R. pisi. The analyses also showed that genomic group III harbored a mixed chromosomal sequence of different genomic groups, suggesting a recent horizontal gene transfer between diverse rhizobia. Although two Canadian strains belonged to subgroup Ia, molecular and physiological analyses showed the divergence between Canadian and Japanese strains. Phylogenetic analysis of nod genes divided the rhizobial strains into several groups that reflected the host range of rhizobia. Symbiosis between dispersing legumes and rhizobia at seashore is discussed.     

Tempos of gene locus deletions and duplications and their relationship to recombination rate during diploid and polyploid evolution in the Aegilops-Triticum alliance

The origin of tetraploid wheat and the divergence of diploid ancestors of wheat A and D genomes were estimated to have occurred 0.36 and 2.7 million years ago, respectively. These estimates and the evolutionary history of 3159 gene loci were used to estimate the rates with which gene loci have been deleted and duplicated during the evolution of wheat diploid ancestors and during the evolution of polyploid wheat. During diploid evolution, the deletion rate was 2.1 x 10(-3) locus(-1) MY(-1) for single-copy loci and 1.0 x 10(-2) locus(-1) MY(-1) for loci in paralogous sets. Loci were duplicated with a rate of 2.9 x 10(-3) locus(-1) MY(-1) during diploid evolution. During polyploid evolution, locus deletion and locus duplication rates were 1.8 x 10(-2) and 1.8 x 10(-3) locus(-1) MY(-1), respectively. Locus deletion and duplication rates correlated positively with the distance of the locus from the centromere and the recombination rate during diploid evolution. The functions of deleted and duplicated loci were inferred to gain insight into the surprisingly high rate of deletions of loci present apparently only once in a genome. The significance of these findings for genome evolution at the diploid and polyploid level is discussed.     

Related polymorphic F-box protein genes between haplotypes clustering in the BAC contig sequences around the S-RNase of Japanese pear

Most fruit trees in the Rosaceae exhibit self-incompatibility, which is controlled by the pistil S gene, encoding a ribonuclease (S-RNase), and the pollen S gene at the S-locus. The pollen S in Prunus is an F-box protein gene (SLF/SFB) located near the S-RNase, but it has not been identified in Pyrus and Malus. In the Japanese pear, various F-box protein genes (PpSFBB(-α-γ)) linked to the S-RNase are proposed as the pollen S candidate. Two bacterial artificial chromosome (BAC) contigs around the S-RNase genes of Japanese pear were constructed, and 649?kb around S(4)-RNase and 378?kb around S(2)-RNase were sequenced. Six and 10 pollen-specific F-box protein genes (designated as PpSFBB(4-u1-u4, 4-d1-d2) and PpSFBB(2-u1-u5,) (2-d1-d5), respectively) were found, but PpSFBB(4-α-γ) and PpSFBB(2-γ) were absent. The PpSFBB(4) genes showed 66.2-93.1% amino acid identity with the PpSFBB(2) genes, which indicated clustering of related polymorphic F-box protein genes between haplotypes near the S-RNase of the Japanese pear. Phylogenetic analysis classified 36 F-box protein genes of Pyrus and Malus into two major groups (I and II), and also generated gene pairs of PpSFBB genes and PpSFBB/Malus F-box protein genes. Group I consisted of gene pairs with 76.3-94.9% identity, while group II consisted of gene pairs with higher identities (>92%) than group I. This grouping suggests that less polymorphic PpSFBB genes in group II are non-S pollen genes and that the pollen S candidates are included in the group I PpSFBB genes.