Naoko Takezaki

POPTREE2: Software for Constructing Population Trees from Allele Frequency Data and Computing Other Population Statistics with Windows Interface

Currently, there is a demand for software to analyze polymorphism data such as microsatellite DNA and single nucleotide polymorphism with easily accessible interface in many fields of research. In this article, we would like to make an announcement of POPTREE2, a computer program package, that can perform evolutionary analyses of allele frequency data. The original version (POPTREE) was a command-line program that runs on the Command Prompt of Windows and Unix. In POPTREE2 genetic distances (measures of the extent of genetic differentiation between populations) for constructing phylogenetic trees, average heterozygosities (H) (a measure of genetic variation within populations) and G(ST) (a measure of genetic differentiation of subdivided populations) are computed through a simple and intuitive Windows interface. It will facilitate statistical analyses of polymorphism data for researchers in many different fields. POPTREE2 is available at http://www.med.kagawa-u.ac.jp/ approximately genomelb/takezaki/poptree2/index.html.

The origin and age of haplochromine fishes in Lake Victoria, east Africa.

According to a widely held view, the more than 300 species of haplochromine cichlid fishes in Lake Victoria (LV), East Africa, originated from a single founder species in less than 12 000 years. This view, however, does not follow from the published geological and molecular evidence. The former does indeed suggest that the LV basin dried out less than 15 000 years ago, but it does not provide any information about the species that re–colonized the new lake or that remained in the rivers draining the area. The molecular evidence is inconclusive with respect to the origin of the LV haplochromines because cichlids from critical regions around LV were not adequately sampled; and as far as the age of the LV haplochromines is concerned, it in fact led to an estimate of 250 000–750 000 years old. In the present study, mitochondrial DNA (control region) variation was determined by heteroduplex and sequencing analyses of more than 670 specimens collected at widely distributed East African riverine and lacustrine localities. The analyses revealed the existence of seven haplogroups (I–VII) distinguishable by characteristic substitutions. All endemic LV samples tested fell into one of these haplogroups (V) which, however, was also found to be present at various other localities, both riverine and lacustrine, outside LV . Within this haplogroup, four subgroups (VA through VD) could be distinguished, two of which (VB and VC) were represented in LV and at other localities. The great majority of the LV haplochromine species could be classified as belonging to the VC subgroup, which was found only in LVand in the rivers draining into it. Hence, while the endemic haplochromine species of LVcould not have originated from a single founding population, the lake does harbour a large species flock which probably arose in situ.

Empirical Tests of the Reliability of Phylogenetic Trees Constructed With Microsatellite DNA

Microsatellite DNA loci or short tandem repeats (STRs) are abundant in eukaryotic genomes and are often used for constructing phylogenetic trees of closely related populations or species. These phylogenetic trees are usually constructed by using some genetic distance measure based on allele frequency data, and there are many distance measures that have been proposed for this purpose. In the past the efficiencies of these distance measures in constructing phylogenetic trees have been studied mathematically or by computer simulations. Recently, however, allele frequencies of 783 STR loci have been compiled from various human populations. We have therefore used these empirical data to investigate the relative efficiencies of different distance measures in constructing phylogenetic trees. The results showed that (1) the probability of obtaining the correct branching pattern of a tree (PC) is generally highest for DA distance; (2) FST*, standard genetic distance (DS), and \batchmode \documentclass[fleqn,10pt,legalpaper]{article} \usepackage{amssymb} \usepackage{amsfonts} \usepackage{amsmath} \pagestyle{empty} \begin{document} \(F_{\mathrm{ST}}/(1{-}F_{\mathrm{ST}})\) \end{document} give similar PC-values, FST* being slightly better than the other two; and (3) (δμ)2 shows PC-values much lower than the other distance measures. To have reasonably high PC-values for trees similar to ours, at least 30 loci with a minimum of 15 individuals are required when DA distance is used.

POPTREEW: Web version of POPTREE for constructing population trees from allele frequency data and computing some other quantities

POPTREE software, including the command line (POPTREE) and the Windows (POPTREE2) versions, is available to perform evolutionary analyses of allele frequency data, computing distance measures for constructing population trees and average heterozygosity (H) (measure of genetic diversity within populations) and G(ST) (measure of genetic differentiation among subdivided populations). We have now developed a web version POPTREEW (http://www.med.kagawa-u.ac.jp/∼genomelb/takezaki/poptreew/) to provide cross-platform access to all POPTREE functions including interactive tree editing. Furthermore, new POPTREE software (POPTREE, POPTREE2, and POPTREEW) computes standardized G(ST) and Josts D, which may be appropriate for data with high variability, and accepts genotype data in GENEPOP format as an input.

Mhc class II B gene evolution in East African cichlid fishes.

Abstract A distinctive feature of essential major histocompatibility complex (Mhc) loci is their polymorphism characterized by large genetic distances between alleles and long persistence times of allelic lineages. Since the lineages often span several successive speciations, we investigated the behavior of the Mhc alleles during or close to the speciation phase. We sequenced exon 2 of the class II B locus 4 from 232 East African cichlid fishes representing 32 related species. The divergence times of the (sub)species ranged from 6000 to 8.4 million years. Two types of evolutionary analysis were used to elucidate the pattern of exon 2 sequence divergence. First, phylogenetic methods were applied to reconstruct the most likely evolutionary pathways leading from the last common ancestor of the set to the extant sequences, and to assess the probable mechanisms involved in allelic diversification. Second, pairwise comparisons of sequences were carried out to detect differences seemingly incompatible with origin by nonparallel point mutations. The analysis revealed point mutations to be the most important mechanism behind allelic divergences, with recombination playing only an auxiliary part. Comparison of sequences from related species revealed evidence of random allelic (lineage) losses apparently associated with speciation. Sharing of identical alleles could be demonstrated between species that diverged 2 million years ago. The phylogeny of the exon was incongruent with that of the flanking introns, indicating either a high degree of convergent evolution at the peptide-binding region-encoding sites, or intron homogenization.

Inconsistency of the maximum parsimony method when the rate of nucleotide substitution is constant

The inconsistency of the maximum parsimony method is known to occur even when the rate of nucleotide substitution is constant. To understand why this inconsistency occurs, a mathematical study was conducted for the cases of five, six, and seven sequences. The results obtained indicate that this inconsistency occurs because the probability of occurrence of nucleotide configurations generated by one substitution on a short interior branch is often lower than that of configurations generated by more substitutions on other longer branches. The chance of occurrence of this event—or, the inconsistency of the maximum parsimony method—apparently increases as the number of sequences increases. The inconsistency may occur even when the extent of sequence divergence is relatively small.

Y‐chromosomal Binary Haplogroups in the Japanese Population and their Relationship to 16 Y‐STR Polymorphisms

We investigated Y chromosomal binary and STR polymorphisms in 263 unrelated male individuals from the Japanese population and further examined the relationships between the two separate types of data. Using 47 biallelic markers we distinguished 20 haplogroups, four of which (D2b1/‐022457, O3/‐002611*, O3/‐LINE1 del, and O3/‐021354*) were newly defined in this study. Most haplogroups in the Japanese population are found in one of the three major clades, C, D, or O. Among these, two major lineages, D2b and O2b, account for 66% of Japanese Y chromosomes. Haplotype diversity of binary markers was calculated at 86.3%. The addition of 16 Y‐STR markers increased the number of haplotypes to 225, yielding a haplotype diversity of 99.40%. A comparison of binary haplogroups and Y‐STR type revealed a close association between certain binary haplogroups and Y‐STR allelic or conformational differences, such as those at the DXYS156Y, DYS390m, DYS392, DYS437, DYS438 and DYS388 loci. Based on our data on the relationships between binary and STR polymorphisms, we estimated the binary haplogroups of individuals from STR haplotypes and frequencies of binary haplogroups in other Japanese, Korean and Taiwanese Han populations. The present data will enable researchers to connect data from binary haplogrouping in anthropological studies and Y‐STR typing in forensic studies in East Asian populations, especially those in and around Japan.

Phylogenetic relationships among East African haplochromine fish as revealed by short interspersed elements (SINEs)

Genomic DNA libraries were prepared from two endemic species of Lake Victoria haplochromine (cichlid) fish and used to isolate and characterize a set of short interspersed elements (SINEs). The distribution and sequences of the SINEs were used to infer phylogenetic relationships among East African haplochromines. The SINE-based classification divides the fish into four groups, which, in order of their divergence from a stem lineage, are the endemic Lake Tanganyika flock (group 1); fish of the nonendemic, monotypic, widely distributed genus Astatoreochromis (group 2); the endemic Lake Malawi flock (group 3); and group 4, which contains fish from widely dispersed East African localities including Lakes Victoria, Edward, George, Albert, and Rukwa, as well as many rivers. The group 4 haplochromines are characterized by a subset of polymorphic SINEs, each of which is present in some individuals and absent in others of the same population at a given locality, the same morphologically defined species, and the same mtDNA-defined haplogroup. SINE-defined group 4 contains six of the seven previously described mtDNA haplogroups. One of the polymorphic SINEs appears to be fixed in the endemic Lake Victoria flock; four others display the presence-or-absence polymorphism within the species of this flock. These findings have implications for the origin of Lake Victoria cichlids and for their founding population sizes.

Phylogeny of antigen-processing enzymes: cathepsins of a cephalochordate, an agnathan and a bony fish.

Cathepsins are enzymes that have been cleaving peptide bonds of lysosomal proteins probably since lysosomes appeared in early eucaryotes. When the adaptive system emerged in gnathostomes, cathepsins were recruited to produce peptides for loading onto the major histocompatibility complex class II molecules and for degrading the class II‐associated invariant chain just before the loading. The circumstances under which this recruitment took place are unclear because the knowledge about vertebrate cathepsins is limited largely to mammals. To shed light on the recruitment, 10 amphioxus, one lamprey and one cichlid fish cathepsin cDNA clone were characterized and analysed phylogenetically. Disregarding cathepsin O, whose phylogenetic position is uncertain, the analysis confirms the existence of two old lines of descent, the B and the L lineages of cathepsins, which diverged from each other early in the evolution of eucaryotes. The B lineage encompasses cathepsins B, C and Z (X). The L lineage splits off sublineages encompassing cathepsins F and W before the plant–animal separation and cathepsin H early in the evolution of the metazoa. The remaining cathepsins belonging to the L lineage diverged from one another during the evolution of vertebrates: S, K and L before the emergence of bony fishes, and the group of rodent placentally expressed cathepsins [J (P), M, Q, R, 3, 6, 7 and 8] as well as the testis/ova‐expressed cathepsins (testins) probably after the divergence of rodents from primates. The part possibly played by the adaptive immune system in some of these divergences is discussed.

Macrophage migration inhibitory factor (MIF) of jawed and jawless fishes: implications for its evolutionary origin.

The macrophage migration inhibitory factor (MIF) is a cytokine produced by T lymphocytes and macrophages in response to inflammatory stimuli. We sequenced MIF cDNA clones of two jawless fishes, the sea lamprey (Petromyzon marinus) and the North Atlantic hagfish (Myxine glutinosa), as well as of the jawed (cichlid) fish Paralabidochromis chilotes. The fish MIF-encoding genes have the same exon-intron organization as the mammalian MIF genes and are present in one copy per haploid genome. Secondary and tertiary structure predictions suggest that the fish MIF proteins have a topology characteristic of the entire MIF-family of proteins. Phylogenetic analysis separates the known nematode members of the family into two groups, one having a sister group relationship with the mammalian D-dopachrome tautomerase (DDT) proteins and the other being related to vertebrate MIFs. It also reveals a high degree of convergent evolution among the members of the family. Finally, it suggests that the divergence of MIF and DDT occurred before the emergence of nematodes in metazoan evolution.