The Reproducibility of an Inferred Tree and the Diploidization of Gene Segregation after Genome Duplication.

We previously introduced a numerical quantity called the stability (Ps) of an inferred tree and showed that for the tree to be reliable this stability as well as the reliability of the tree, which is usually computed as the bootstrap probability (Pb), must be high. However, if genome duplication occurs in a species, a gene family of the genome also duplicates, and for this reason alone some Ps values can be high in a tree of the duplicated gene families. In addition, the topology of the duplicated gene family can be similar to that of the original gene family if such gene families are identifiable. After genome duplication, however, the gene families are often partially deleted or partially duplicated, and the duplicated gene family may not show the same topology as that of the original family. It is therefore necessary to compute the similarity of the topologies of the duplicated and the original gene families. In this paper, we introduce another quantity called the reproducibility (Pr) for measuring the similarity of the two gene families. To show how to compute the Pr values, we first compute the Pb and Ps values for each of the MHC class II α and ß chain gene families, which were apparently generated by genome duplication. We then compute the Pr values for the MHC class II α and ß chain gene families. The Pr values for the α and ß chain gene families are now low, and this suggests that the diploidization of gene segregation has occurred after the genome duplication. Currently higher animals, defined as animals with complex phenotypic characters, generally have a higher genome size, and this increase in genome size appears to have been caused by genome duplication and diploidization of gene segregation after genome duplication.

The Reliability and Stability of an Inferred Phylogenetic Tree from Empirical Data.

The reliability of a phylogenetic tree obtained from empirical data is usually measured by the bootstrap probability (Pb) of interior branches of the tree. If the bootstrap probability is high for most branches, the tree is considered to be reliable. If some interior branches show relatively low bootstrap probabilities, we are not sure that the inferred tree is really reliable. Here, we propose another quantity measuring the reliability of the tree called the stability of a subtree. This quantity refers to the probability of obtaining a subtree (Ps) of an inferred tree obtained. We then show that if the tree is to be reliable, both Pb and Ps must be high. We also show that Ps is given by a bootstrap probability of the subtree with the closest outgroup sequence, and computer program RESTA for computing the Pb and Ps values will be presented.

Efficiencies of the NJp, Maximum Likelihood, and Bayesian Methods of Phylogenetic Construction for Compositional and Noncompositional Genes.

At the present time it is often stated that the maximum likelihood or the Bayesian method of phylogenetic construction is more accurate than the neighbor joining (NJ) method. Our computer simulations, however, have shown that the converse is true if we use p distance in the NJ procedure and the criterion of obtaining the true tree (Pc expressed as a percentage) or the combined quantity (c) of a value of Pc and a value of Robinson-Foulds' average topological error index (dT). This c is given by Pc (1 - dT/dTmax) = Pc (m - 3 - dT/2)/(m - 3), where m is the number of taxa used and dTmax is the maximum possible value of dT, which is given by 2(m - 3). This neighbor joining method with p distance (NJp method) will be shown generally to give the best data-fit model. This c takes a value between 0 and 1, and a tree-making method giving a high value of c is considered to be good. Our computer simulations have shown that the NJp method generally gives a better performance than the other methods and therefore this method should be used in general whether the gene is compositional or it contains the mosaic DNA regions or not.

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 Jost's D, which may be appropriate for data with high variability, and accepts genotype data in GENEPOP format as an input.

Evolution of the sex-lethal gene in insects and origin of the sex-determination system in Drosophila.

Sex-lethal (Sxl) functions as the switch gene for sex-determination in Drosophila melanogaster by engaging a regulatory cascade. Thus far the origin and evolution of both the regulatory system and SXL protein's sex-determination function have remained largely unknown. In this study, we explore systematically the Sxl homologs in a wide range of insects, including the 12 sequenced Drosophila species, medfly, blowflies, housefly, Megaselia scalaris, mosquitoes, butterfly, beetle, honeybee, ant, and aphid. We find that both the male-specific and embryo-specific exons exist in all Drosophila species. The homologous male-specific exon is also present in Scaptodrosophila lebanonensis, but it does not have in-frame stop codons, suggesting the exon's functional divergence between Drosophila and Scaptodrosophila after acquiring it in their common ancestor. Two motifs closely related to the exons' functions, the SXL binding site poly(U) and the transcription-activating motif TAGteam, surprisingly exhibit broader phylogenetic distributions than the exons. Some previously unknown motifs that are restricted to or more abundant in Drosophila and S. lebanonensis than in other insects are also identified. Finally, phylogenetic analysis suggests that the SXL's novel sex-determination function in Drosophila is more likely attributed to the changes in the N- and C-termini rather than in the RNA-binding region. Thus, our results provide a clearer picture of the phylogeny of the Sxl's cis-regulatory elements and protein sequence changes, and so lead to a better understanding of the origin of sex-determination in Drosophila and also raise some new questions regarding the evolution of Sxl.

Random fluctuation of selection coefficients and the extent of nucleotide variation in human populations.

It is well known that the selection coefficient of a mutant allele varies from generation to generation, and the effect of this factor on genetic variation has been studied by many theoreticians. However, no consensus has been reached. One group of investigators believes that fluctuating selection has an effect of enhancing genetic variation, whereas the other group contends that it has a diversity-reducing effect. In recent years, it has become possible to study this problem by using single nucleotide polymorphisms (SNPs) as well as exome sequence data. In this article we present the theoretical distributions of mutant nucleotide frequencies for the two models of fluctuating selection and then compare the distributions with the empirical distributions obtained from SNP and exome sequence data in human populations. Interestingly, both SNP and exome sequence data showed that the neutral mutation model fits the empirical distribution quite well. Furthermore, the mathematical models with diversity-enhancing and diversity-reducing effects also fit the empirical distribution reasonably well. This result implies that there is no need of distinguishing between the diversity-enhancing and diversity-reducing models of fluctuating selection and the nucleotide polymorphism in human populations can be explained largely by neutral mutations when long-term evolution is considered.

Origins and evolution of microRNA genes in plant species.

MicroRNAs (miRNAs) are among the most important regulatory elements of gene expression in animals and plants. However, their origin and evolutionary dynamics have not been studied systematically. In this paper, we identified putative miRNA genes in 11 plant species using the bioinformatic technique and examined their evolutionary changes. Our homology search indicated that no miRNA gene is currently shared between green algae and land plants. The number of miRNA genes has increased substantially in the land plant lineage, but after the divergence of eudicots and monocots, the number has changed in a lineage-specific manner. We found that miRNA genes have originated mainly by duplication of preexisting miRNA genes or protein-coding genes. Transposable elements also seem to have contributed to the generation of species-specific miRNA genes. The relative importance of these mechanisms in plants is quite different from that in Drosophila species, where the formation of hairpin structures in the genomes seems to be a major source of miRNA genes. This difference in the origin of miRNA genes between plants and Drosophila may be explained by the difference in the binding to target mRNAs between plants and animals. We also found that young miRNA genes are less conserved than old genes in plants as well as in Drosophila species. Yet, nearly half of the gene families in the ancestor of flowering plants have been lost in at least one species examined. This indicates that the repertoires of miRNA genes have changed more dynamically than previously thought during plant evolution.

Roles of mutation and selection in speciation: from Hugo de Vries to the modern genomic era.

One of the most important problems in evolutionary biology is to understand how new species are generated in nature. In the past, it was difficult to study this problem because our lifetime is too short to observe the entire process of speciation. In recent years, however, molecular and genomic techniques have been developed for identifying and studying the genes involved in speciation. Using these techniques, many investigators have already obtained new findings. At present, however, the results obtained are complex and quite confusing. We have therefore attempted to understand these findings coherently with a historical perspective and clarify the roles of mutation and natural selection in speciation. We have first indicated that the root of the currently burgeoning field of plant genomics goes back to Hugo de Vries, who proposed the mutation theory of evolution more than a century ago and that he unknowingly found the importance of polyploidy and chromosomal rearrangements in plant speciation. We have then shown that the currently popular Dobzhansky-Muller model of evolution of reproductive isolation is only one of many possible mechanisms. Some of them are Oka's model of duplicate gene mutations, multiallelic speciation, mutation-rescue model, segregation-distorter gene model, heterochromatin-associated speciation, single-locus model, etc. The occurrence of speciation also depends on the reproductive system, population size, bottleneck effects, and environmental factors, such as temperature and day length. Some authors emphasized the importance of natural selection to speed up speciation, but mutation is crucial in speciation because reproductive barriers cannot be generated without mutations.

MEGA5: molecular evolutionary genetics analysis using maximum likelihood, evolutionary distance, and maximum parsimony methods.

Comparative analysis of molecular sequence data is essential for reconstructing the evolutionary histories of species and inferring the nature and extent of selective forces shaping the evolution of genes and species. Here, we announce the release of Molecular Evolutionary Genetics Analysis version 5 (MEGA5), which is a user-friendly software for mining online databases, building sequence alignments and phylogenetic trees, and using methods of evolutionary bioinformatics in basic biology, biomedicine, and evolution. The newest addition in MEGA5 is a collection of maximum likelihood (ML) analyses for inferring evolutionary trees, selecting best-fit substitution models (nucleotide or amino acid), inferring ancestral states and sequences (along with probabilities), and estimating evolutionary rates site-by-site. In computer simulation analyses, ML tree inference algorithms in MEGA5 compared favorably with other software packages in terms of computational efficiency and the accuracy of the estimates of phylogenetic trees, substitution parameters, and rate variation among sites. The MEGA user interface has now been enhanced to be activity driven to make it easier for the use of both beginners and experienced scientists. This version of MEGA is intended for the Windows platform, and it has been configured for effective use on Mac OS X and Linux desktops. It is available free of charge from http://www.megasoftware.net.

Evolutionary changes of the target sites of two microRNAs encoded in the Hox gene cluster of Drosophila and other insect species.

MicroRNAs (miRs) are noncoding RNAs that regulate gene expression at the post-transcriptional level. In animals, the target sites of a miR are generally located in the 3' untranslated regions (UTRs) of messenger RNAs. However, how the target sites change during evolution is largely unknown. MiR-iab-4 and miR-iab-4as are known to regulate the expression of two Hox genes, Abd-A and Ubx, in Drosophila melanogaster. We have therefore studied the evolutionary changes of these two miR genes and their target sites of the Hox genes in Drosophila, other insect species, and Daphnia. Our homology search identified a single copy of each miR gene located in the same genomic position of the Hox gene cluster in all species examined. The seed nucleotide sequence was also the same for all species. Searching for the target sites in all Hox genes, we found several target sites of miR-iab-4 and miR-iab-4as in Antp in addition to Abd-A and Ubx in most insect species examined. Our phylogenetic analysis of target sites in Abd-A, Ubx, and Antp showed that the old target sites, which existed before the divergence of the 12 Drosophila species, have been well maintained in most species under purifying selection. By contrast, new target sites, which were generated during Drosophila evolution, were often lost in some species and mostly located in unalignable regions of the 3' UTRs. These results indicate that these regions can be a potential source of generating new target sites, which results in multiple target genes for each miR in animals.

Origins and evolution of microRNA genes in Drosophila species.

MicroRNAs (miRs) regulate gene expression at the posttranscriptional level. To obtain some insights into the origins and evolutionary patterns of miR genes, we have identified miR genes in the genomes of 12 Drosophila species by bioinformatics approaches and examined their evolutionary changes. The results showed that the extant and ancestral Drosophila species had more than 100 miR genes and frequent gains and losses of miR genes have occurred during evolution. Although many miR genes appear to have originated from random hairpin structures in intronic or intergenic regions, duplication of miR genes has also contributed to the generation of new miR genes. Estimating the rate of nucleotide substitution of miR genes, we have found that newly arisen miR genes have a substitution rate similar to that of synonymous nucleotide sites in protein-coding genes and evolve almost neutrally. This suggests that most new miR genes have not acquired any important function and would become inactive. By contrast, old miR genes show a substitution rate much lower than the synonymous rate. Moreover, paired and unpaired nucleotide sites of miR genes tend to remain unchanged during evolution. Therefore, once miR genes acquired their functions, they appear to have evolved very slowly, maintaining essentially the same structures for a long time.

The neutral theory of molecular evolution in the genomic era.

The neutral theory of molecular evolution has been widely accepted and is the guiding principle for studying evolutionary genomics and the molecular basis of phenotypic evolution. Recent data on genomic evolution are generally consistent with the neutral theory. However, many recently published papers claim the detection of positive Darwinian selection via the use of new statistical methods. Examination of these methods has shown that their theoretical bases are not well established and often result in high rates of false-positive and false-negative results. When the deficiencies of these statistical methods are rectified, the results become largely consistent with the neutral theory. At present, genome-wide analyses of natural selection consist of collections of single-locus analyses. However, because phenotypic evolution is controlled by the interaction of many genes, the study of natural selection ought to take such interactions into account. Experimental studies of evolution will also be crucial.

Analysis of the immunoglobulin light chain genes in zebra finch: evolutionary implications.

All jawed vertebrates produce immunoglobulins (IGs) as a defense mechanism against pathogens. Typically, IGs are composed of two identical heavy chains (IGH) and two identical light chains (IGL). Most tetrapod species encode more than one isotype of light chains. Chicken is the only representative of birds for which genomic information is currently available and is an exception to the above rule because it encodes only a single IGL isotype (i.e., lambda). Here, we show that the genome of zebra finch, another bird species, encodes a single IGL isotype, that is, lambda, like the chicken. These results strongly suggest that the second isotype (i.e., kappa) present in both reptiles and mammals was lost in a very early stage of bird evolution. Furthermore, we show that both chicken and zebra finch contain a single set of functional variable, joining, and constant region genes and multiple variable region pseudogenes. The latter finding suggests that this type of genomic organization was already present in the common ancestor of these bird species and remained unchanged over a long evolutionary time. This conservation is in contrast with the high levels of variation observed in the mammalian IGL loci. The presence of a single functional variable region gene followed by multiple variable pseudogenes in zebra finch suggest that this species may be generating antibody diversity by a gene conversion-like mechanism like the chicken.

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.

Genomic organization and evolution of immunoglobulin kappa gene enhancers and kappa deleting element in mammals.

We have studied the genomic structure and evolutionary pattern of immunoglobulin kappa deleting element (KDE) and three kappa enhancers (KE5', KE3'P, and KE3'D) in eleven mammalian genomic sequences. Our results show that the relative positions and the genomic organization of the KDE and the kappa enhancers are conserved in all mammals studied and have not been affected by the local rearrangements in the immunoglobulin kappa (IGK) light chain locus over a long evolutionary time ( approximately 120 million years of mammalian evolution). Our observations suggest that the sequence motifs in these regulatory elements have been conserved by purifying selection to achieve proper regulation of the expression of the IGK light chain genes. The conservation of the three enhancers in all mammals indicates that these species may use similar mechanisms to regulate IGK gene expression. However, some activities of the IGK enhancers might have evolved in the eutherian lineage. The presence of the three IGK enhancers, KDE, and other recombining elements (REs) in all mammals (including platypus) suggest that these genomic elements were in place before the mammalian radiation.

Genomic drift and evolution of microsatellite DNAs in human populations.

In recent years, copy number variation (CNV) of DNA segments has become a hot topic in the study of genetic variation, and a large amount of CNVs has been uncovered in human populations. The CNVs involving the smallest units of DNA segments are microsatellite DNAs, and the evolutionary change of microsatellite DNAs is believed to occur mostly by the increase or decrease of one repeat unit at a time in a more or less neutral fashion. If we note that eukaryotic genomes contain millions of microsatellite loci, this pattern of nucleotide change is expected to generate random changes of genome size, that is, genomic drift, and will provide a neutral model of CNV evolution. We therefore investigated the amount of variation of the total number of repeats (TNR) per individual concerned with 145 microsatellite loci in three human populations, Africans, Europeans, and Asians. It was shown that the TNR follows the normal distribution in all three populations and that the extent of variation of TNR is more than 50% greater in Africans than in Europeans and Asians as expected from the hypothesis of African origin of modern humans. If we consider all microsatellite loci in the human genome and compute the variation of the total number of nucleotides involved (TNN), it is possible to study the contribution of microsatellite loci to the genome size variation. This study has shown that the genome sizes of human individuals are affected considerably by genomic drift of microsatellite DNA alone. This pattern of evolution is similar to that of olfactory receptor (OR) genes previously studied in human populations and support the idea that the number of OR genes has evolved in a more or less neutral fashion. However, this conclusion does not necessarily apply to the genomewide CNVs of various DNA segments, and it appears that long variant DNA fragments are deleterious and under purifying selection.

Reliabilities of identifying positive selection by the branch-site and the site-prediction methods.

Natural selection operating in protein-coding genes is often studied by examining the ratio (omega) of the rates of nonsynonymous to synonymous nucleotide substitution. The branch-site method (BSM) based on a likelihood ratio test is one of such tests to detect positive selection for a predetermined branch of a phylogenetic tree. However, because the number of nucleotide substitutions involved is often very small, we conducted a computer simulation to examine the reliability of BSM in comparison with the small-sample method (SSM) based on Fisher's exact test. The results indicate that BSM often generates false positives compared with SSM when the number of nucleotide substitutions is approximately 80 or smaller. Because the omega value is also used for predicting positively selected sites, we examined the reliabilities of the site-prediction methods, using nucleotide sequence data for the dim-light and color vision genes in vertebrates. The results showed that the site-prediction methods have a low probability of identifying functional changes of amino acids experimentally determined and often falsely identify other sites where amino acid substitutions are unlikely to be important. This low rate of predictability occurs because most of the current statistical methods are designed to identify codon sites with high omega values, which may not have anything to do with functional changes. The codon sites showing functional changes generally do not show a high omega value. To understand adaptive evolution, some form of experimental confirmation is necessary.

Evolution of F-box genes in plants: different modes of sequence divergence and their relationships with functional diversification.

F-box proteins are substrate-recognition components of the Skp1-Rbx1-Cul1-F-box protein (SCF) ubiquitin ligases. In plants, F-box genes form one of the largest multigene superfamilies and control many important biological functions. However, it is unclear how and why plants have acquired a large number of F-box genes. Here we identified 692, 337, and 779 F-box genes in Arabidopsis, poplar and rice, respectively, and studied their phylogenetic relationships and evolutionary patterns. We found that the plant F-box superfamily can be divided into 42 families, each of which has a distinct domain organization. We also estimated the number of ancestral genes for each family and identified highly conservative versus divergent families. In conservative families, there has been little or no change in the number of genes since the divergence between eudicots and monocots approximately 145 million years ago. In divergent families, however, the numbers have increased dramatically during the same period. In two cases, the numbers of genes in extant species are >100 times greater than that in the most recent common ancestor (MRCA) of the three species. Proteins encoded by highly conservative genes always have the same domain organization, suggesting that they interact with the same or similar substrates. In contrast, proteins of rapidly duplicating genes sometimes have quite different domain structures, mainly caused by unusually frequent shifts of exon-intron boundaries and/or frameshift mutations. Our results indicate that different F-box families, or different clusters of the same family, have experienced dramatically different modes of sequence divergence, apparently having resulted in adaptive changes in function.