Mesial temporal lobe epilepsy with childhood febrile seizure

OBJECTIVES

To evaluate the demographic and clinical manifestations of patients with mesial temporal sclerosis and temporal lobe epilepsy (MTS-TLE) with childhood febrile seizure (FS) and establishing the potential differences as compared to those without FS. We also investigated the surgery outcome in these two groups of patients.

MATERIALS AND METHODS

In this retrospective study, all patients with a clinical diagnosis of drug-resistant TLE due to mesial temporal sclerosis, who underwent epilepsy surgery at Jefferson Comprehensive Epilepsy Center, were recruited. Patients were prospectively registered in a database from 1986 through 2014. Postsurgical outcome was classified into two groups; seizure-free or relapsed. Clinical manifestations and outcome were compared between patients with MTS-TLE with FS and those without FS.

RESULTS

Two hundred and sixty-two patients were eligible for this study. One hundred and seventy patients (64.9%) did not have FS in their childhood, while 92 patients (35.1%) reported experiencing FS in their childhood. Demographic and clinical characteristics of these two groups of patients were not different. Postoperative seizure outcome was not statistically different between these two groups of patients (P = 0.19).

CONCLUSIONS

When MTS is the pathological substrate of TLE, clinical manifestations and response to surgical treatment of patients are very similar in patients with history of febrile seizure in their childhood compared to those without such an experience. In other words, when the subgroup of patients with MTS-TLE and drug-resistant seizures is examined history of childhood febrile seizure loses its value as a distinguishing factor in characteristics or predictive factor for surgery outcome.

Seasonality of birth in patients with temporal lobe epilepsy

OBJECTIVES

We investigated the seasonal pattern in births of patients with temporal lobe epilepsy and mesial temporal sclerosis. We hypothesized that the seasonal pattern in births of these patients is different from that in the general population.

MATERIALS AND METHODS

In this retrospective study, all patients who were evaluated for epilepsy surgery at Jefferson Comprehensive Epilepsy Center, Thomas Jefferson, Philadelphia, USA, between 1986 and 2014 and had a diagnosis of mesial temporal sclerosis (made by definite imaging findings of atrophy and/or sclerosis) were included. The seasonality in births of patients was compared with the seasonal pattern in the live births of the general population from Pennsylvania State.

RESULTS

Two hundred and eighty-two patients (146 females and 136 males) were studied. The seasonality pattern in birth of patients was not statistically different from that in the general population.

CONCLUSIONS

The observed contradictory findings among various studies indicate the need for further studies to elucidate whether season of birth brings the possibility of acquiring various epilepsy syndromes in the future. To investigate any possible association between season of birth and epilepsy, we suggest avoid pooling all patients with epilepsy together.

Genetic drift and estimation of effective population size

The statistical properties of the standardized variance of gene frequency changes (a quantity equivalent to Wright's inbreeding coefficient) in a random mating population are studied, and new formulae for estimating the effective population size are developed. The accuracy of the formulae depends on the ratio of sample size to effective size, the number of generations involved (t), and the number of loci or alleles used. It is shown that the standardized variance approximately follows the chi(2) distribution unless t is very large, and the confidence interval of the estimate of effective size can be obtained by using this property. Application of the formulae to data from an isolated population of Dacus oleae has shown that the effective size of this population is about one tenth of the minimum census size, though there was a possibility that the procedure of sampling genes was improper.

Statistical Studies on Protein Polymorphism in Natural Populations. III. Distribution of Allele Frequencies and the Number of Alleles per Locus

With the aim of understanding the mechanism of maintenance of protein polymorphism, we have studied the properties of allele frequency distribution and the number of alleles per locus, using gene-frequency data from a wide range of organisms (mammals, birds, reptiles, amphibians, Drosophila and non-Drosophila invertebrates) in which 20 or more loci with at least 100 genes were sampled. The observed distribution of allele frequencies was U-shaped in all of the 138 populations (mostly species or subspecies) examined and generally agreed with the theoretical distribution expected under the mutation-drift hypothesis, though there was a significant excess of rare alleles (gene frequency, 0 approximately 0.05) in about a quarter of the populations. The agreement between the mutation-drift theory and observed data was quite satisfactory for the numbers of polymorphic (gene frequency, 0.05 approximately 0.95) and monomorphic (0.95 approximately 1.0) alleles.-The observed pattern of allele-frequency distribution was incompatible with the prediction from the overdominance hypothesis. The observed correlations of the numbers of rare alleles, polymorphic alleles and monomorphic alleles with heterozygosity were of the order of magnitude that was expected under the mutation-drift hypothesis. Our results did not support the view that intracistronic recombination is an important source of genetic variation. The total number of alleles per locus was positively correlated with molecular weight in most of the species examined, and the magnitude of the correlation was consistent with the theoretical prediction from mutation-drift hypothesis. The correlation between molecular weight and the number of alleles was generally higher than the correlation between molecular weight and heterozygosity, as expected.

Genomic drift and copy number variation of chemosensory receptor genes in humans and mice

Recent studies about the structural variation of genomic sequences have shown that there is a large amount of copy number variations (CNVs) of genes within species. Analyzing Redon et al.'s (2006) crude data on copy number variable regions (CNVRs), we previously showed that CNVs are particularly high for chemosensory receptor genes in human populations. In this paper, we reanalyzed the CNVs of these genes using more refined data by Perry et al. (2008). The results showed that the extent of CNVs is somewhat lower in this dataset than in the previous one, but that the extent is still substantial for olfactory receptor (OR), vomeronasal receptor (VR), and taste receptor (TR) genes. We also studied the CNVs for chemosensory receptor genes in mice, using CNVR data obtained from inbred strains. It was found that the extent of CNVs is quite substantial but is lower than that for human populations. However, because the mouse data came from inbred strains and might be biased, this conclusion should be regarded as tentative. Despite this reservation, the distribution of gene copy number for the OR gene family was approximately normal in both humans and mice, suggesting that genomic drift caused by random duplication and deletion of genes plays important roles in determining the evolutionary change of chemosensation.

Linear IgA bullous dermatosis of neonatal onset: case report and review of the literature

Several small blisters were noticed on the forehead and the trunk of a newborn boy on day 1. The blisters gradually enlarged and spread over the whole body including the oral mucosa. A skin biopsy was performed twice and subepidermal bullae with polymorphonuclear and mononuclear cell infiltration were demonstrated. Direct immunofluorescence showed linear IgA, IgG and C3 depositions along the basement membrane zone and this finding led to a diagnosis of linear IgA bullous dermatosis. So far, internationally, only one case has ever been reported on the disease at neonatal onset. The skin lesions spontaneously regressed and the mucosal lesions were controlled with diaminodiphenylsulfone.

CONCLUSION

In neonates with prolonged blistering, autoimmune disease such as linear IgA bullous dermatosis should be considered within the differential diagnosis and an immunofluorescence study must be performed.

Reliabilities of parsimony-based and likelihood-based methods for detecting positive selection at single amino acid sites

The reliabilities of parsimony-based and likelihood-based methods for inferring positive selection at single amino acid sites were studied using the nucleotide sequences of human leukocyte antigen (HLA) genes, in which positive selection is known to be operating at the antigen recognition site. The results indicate that the inference by parsimony-based methods is robust to the use of different evolutionary models and generally more reliable than that by likelihood-based methods. In contrast, the results obtained by likelihood-based methods depend on the models and on the initial parameter values used. It is sometimes difficult to obtain the maximum likelihood estimates of parameters for a given model, and the results obtained may be false negatives or false positives depending on the initial parameter values. It is therefore preferable to use parsimony-based methods as long as the number of sequences is relatively large and the branch lengths of the phylogenetic tree are relatively small.

MEGA2: molecular evolutionary genetics analysis software

We have developed a new software package, Molecular Evolutionary Genetics Analysis version 2 (MEGA2), for exploring and analyzing aligned DNA or protein sequences from an evolutionary perspective. MEGA2 vastly extends the capabilities of MEGA version 1 by: (1) facilitating analyses of large datasets; (2) enabling creation and analyses of groups of sequences; (3) enabling specification of domains and genes; (4) expanding the repertoire of statistical methods for molecular evolutionary studies; and (5) adding new modules for visual representation of input data and output results on the Microsoft Windows platform.

AVAILABILITY

http://www.megasoftware.net.

CONTACT

s.kumar@asu.edu

Vagal nerve stimulation: adjustments to reduce painful side effects

Vagal nerve stimulation is an approved adjunctive treatment for medically intractable epilepsy. Although it is generally well tolerated, some patients experience pain, coughing, or hoarseness during stimulation. Lowering the pulse width in these patients alleviates pain and reduces voice alteration without loss of efficacy. This allows more optimal programming of stimulation intensities.

ADAPTSITE: detecting natural selection at single amino acid sites

ADAPTSITE is a program package for detecting natural selection at single amino acid sites, using a multiple alignment of protein-coding sequences for a given phylogenetic tree. The program infers ancestral codons at all interior nodes, and computes the total numbers of synonymous (c(S)) and nonsynonymous (c(N)) substitutions as well as the average numbers of synonymous (s(S)) and nonsynonymous (s(N)) sites for each codon site. The probabilities of occurrence of synonymous and nonsynonymous substitutions are approximated by s(S) / (s(S) + s(N)) and s(N) / (s(S) + s(N)), respectively. The null hypothesis of selective neutrality is tested for each codon site, assuming a binomial distribution for the probability of obtaining c(S) and c(N).

AVAILABILITY

ADAPTSITE is available free of charge at the World-Wide Web sites http://mep.bio.psu.edu/adaptivevol.html and http://www.cib.nig.ac.jp/dda/yossuzuk/welcome.html. The package includes the source code written in C, binary files for UNIX operating systems, manual, and example files.

Evolutionary dynamics of the T-cell receptor VB gene family as inferred from the human and mouse genomic sequences

The diversity of T-cell receptors is generated primarily by the variable-region gene families, each of which is composed of a large number of member genes. The entire genomic sequence of the variable region (VB) of the T- cell receptor beta chain from humans and mice has become available. To understand the evolutionary dynamics of the VB gene family, we conducted a phylogenetic analysis of all VB genes from humans and mice, as well as a detailed analysis of internal DNA duplications in the human genomic VB region. The phylogenetic tree obtained shows that human and mouse VB genes intermingle extensively rather than forming two separate clusters and that many gene duplications occurred both before and after the divergence between primates and rodents. Analyzing the genomic maps of transposable elements (e.g., LINEs and SINEs) and relic VB genes in the VB gene region, we present evidence that a 20-kb VB region duplicated tandemly four times in the human lineage during the last 32 Myr, and 6 out of the 15 VB genes in this region have become nonfunctional during this period. Our results show that the VB gene family is subject to evolution by a birth-and-death process rather than to concerted evolution.

Estimation of divergence times from multiprotein sequences for a few mammalian species and several distantly related organisms

When many protein sequences are available for estimating the time of divergence between two species, it is customary to estimate the time for each protein separately and then use the average for all proteins as the final estimate. However, it can be shown that this estimate generally has an upward bias, and that an unbiased estimate is obtained by using distances based on concatenated sequences. We have shown that two concatenation-based distances, i.e., average gamma distance weighted with sequence length (d(2)) and multiprotein gamma distance (d(3)), generally give more satisfactory results than other concatenation-based distances. Using these two distance measures for 104 protein sequences, we estimated the time of divergence between mice and rats to be approximately 33 million years ago. Similarly, the time of divergence between humans and rodents was estimated to be approximately 96 million years ago. We also investigated the dependency of time estimates on statistical methods and various assumptions made by using sequence data from eubacteria, protists, plants, fungi, and animals. Our best estimates of the times of divergence between eubacteria and eukaryotes, between protists and other eukaryotes, and between plants, fungi, and animals were 3, 1.7, and 1.3 billion years ago, respectively. However, estimates of ancient divergence times are subject to a substantial amount of error caused by uncertainty of the molecular clock, horizontal gene transfer, errors in sequence alignments, etc.

Purifying selection and birth-and-death evolution in the ubiquitin gene family

Ubiquitin is a highly conserved protein that is encoded by a multigene family. It is generally believed that this gene family is subject to concerted evolution, which homogenizes the member genes of the family. However, protein homogeneity can be attained also by strong purifying selection. We therefore studied the proportion (p(S)) of synonymous nucleotide differences between members of the ubiquitin gene family from 28 species of fungi, plants, and animals. The results have shown that p(S) is generally very high and is often close to the saturation level, although the protein sequence is virtually identical for all ubiquitins from fungi, plants, and animals. A small proportion of species showed a low level of p(S) values, but these values appeared to be caused by recent gene duplication. It was also found that the number of repeat copies of the gene family varies considerably with species, and some species harbor pseudogenes. These observations suggest that the members of this gene family evolve almost independently by silent nucleotide substitution and are subjected to birth-and-death evolution at the DNA level.

Efficiencies of fast algorithms of phylogenetic inference under the criteria of maximum parsimony, minimum evolution, and maximum likelihood when a large number of sequences are used

In phylogenetic inference by maximum-parsimony (MP), minimum-evolution (ME), and maximum-likelihood (ML) methods, it is customary to conduct extensive heuristic searches of MP, ME, and ML trees, examining a large number of different topologies. However, these extensive searches tend to give incorrect tree topologies. Here we show by extensive computer simulation that when the number of nucleotide sequences (m) is large and the number of nucleotides used (n) is relatively small, the simple MP or ML tree search algorithms such as the stepwise addition (SA) plus nearest neighbor interchange (NNI) search and the SA plus subtree pruning regrafting (SPR) search are as efficient as the extensive search algorithms such as the SA plus tree bisection-reconnection (TBR) search in inferring the true tree. In the case of ME methods, the simple neighbor-joining (NJ) algorithm is as efficient as or more efficient than the extensive NJ+TBR search. We show that when ME methods are used, the simple p distance generally gives better results in phylogenetic inference than more complicated distance measures such as the Hasegawa-Kishino-Yano (HKY) distance, even when nucleotide substitution follows the HKY model. When ML methods are used, the simple Jukes-Cantor (JC) model of phylogenetic inference generally shows a better performance than the HKY model even if the likelihood value for the HKY model is much higher than that for the JC model. This indicates that at least in the present case, selecting of a substitution model by using the likelihood ratio test or the AIC index is not appropriate. When n is small relative to m and the extent of sequence divergence is high, the NJ method with p distance often shows a better performance than ML methods with the JC model. However, when the level of sequence divergence is low, this is not the case.

Evolution of vertebrate immunoglobulin variable gene segments

Evolution of Ig V gene segments are generally characterized by (a) evolution by "the birth and death process" and (b) diversifying selection. However, the detailed evolutionary pattern of V gene segments varies among species due to the fact that the humoral immune system itself has changed during vertebrate evolution. The change in somatic diversification system coupled with the change in lymphocyte development has imposed a significant impact on the evolution of Ig genes. In order to understand the evolution of immunological genes it is important to view it in the context of the evolution of the entire immune system itself.

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.

EKG abnormalities during partial seizures in refractory epilepsy

PURPOSE

This study assessed the frequency and character of ictal cardiac rhythm and conduction abnormalities in intractable epilepsy. Sudden unexpected death in epilepsy (SUDEP) is a major cause of excess mortality in people with refractory epilepsy, and cardiac arrhythmias during seizures may be responsible. The frequency of cardiac abnormalities during seizures in patients with refractory epilepsy must be determined.

METHODS

Fifty-one seizures in 43 patients with intractable partial epilepsy were analyzed prospectively from CCTV-EEG monitoring with one ECG channel. Arrhythmias, repolarization abnormalities, and PR and QTc intervals were determined for preictal (3 min), ictal, and postictal (3 min) periods for one or more seizures per patient. Parametric statistics were used for continuous variables, and nonparametric statistics were used for categoric variables.

RESULTS

Of the patients, 39% had one or more abnormalities of rhythm and/or repolarization during or immediately after seizures. Abnormalities included asystole (one), atrial fibrillation (one), marked or moderate sinus arrhythmia (six), supraventricular tachycardia (one), atrial premature depolarizations (APDs; eight), ventricular premature depolarizations (VPDs; two), and bundle-branch block (three). Mean seizure duration was longer in patients with abnormalities than in those without (204 vs. 71 s; p < 0.001). Generalized tonic-clonic seizures were also associated with increased occurrence of ictal ECG abnormalities (p = 0.006) as compared with complex partial seizures. There were no clinically significant differences in mean preictal and ictal/postictal PR and QTc intervals.

CONCLUSIONS

Cardiac rhythm and conduction abnormalities are common during seizures, particularly if they are prolonged or generalized, in intractable epilepsy. These abnormalities may contribute to SUDEP.

An unusual form of purifying selection in a sperm protein

Protamines are small, highly basic DNA-binding proteins found in the sperm of animals. Interestingly, the proportion of arginine residues in one type of protamine, protamine P1, is about 50% in mammals. Upon closer examination, it was found that both the total number of amino acids and the positions of arginine residues have changed considerably during the course of mammalian evolution. This evolutionary pattern suggests that protamine P1 is under an unusual form of purifying selection, in which the high proportion of arginine residues is maintained but the positions may vary. In this case, we would expect that the rate of nonsynonymous substitution is not particularly low compared with that of synonymous substitution, despite purifying selection. We would also expect that the selection for a high arginine content results in a high frequency of the nucleotide G in the coding region of this gene, because all six arginine codons contain at least one G. These expectations were confirmed in our study of mammalian protamine genes. Analysis of nonmammalian vertebrate genes also showed essentially the same patterns of evolutionary changes, suggesting that this unusual form of purifying selection has been active since the origin of bony vertebrates. The protamine gene of an insect species shows similar patterns, although its purifying selection is less intense. These observations suggest that arginine-rich selection is a general feature of protamine evolution. The driving force for arginine-rich selection appears to be the DNA-binding function of protamine P1 and an interaction with a protein kinase in the fertilized egg.

Diversity and evolution of T-cell receptor variable region genes in mammals and birds

The receptor of a T lymphocyte (TCR) recognizes nonself antigens in the company of major histocompatibility complex (MHC) molecules presented to it by the antigen-presenting cell. The variable region of TCR is encoded by either a concatenation of variable region (TCR-V), diversity region (TCR-D), and joining region (TCR-J) genes, or a concatenation of TCR-V and TCR-J genes. The TCR-V genes exist as a multigene family in vertebrate species. Here we study the evolutionary relationships of TCR-V genes from humans, sheep, cattle, rabbits, mice, and chicken. These six species can be classified into two groups according to the frequency of gamma(delta) T-cells in their peripheral T-cell populations. The "gamma(delta) low" group of species includes humans and mice, in which gamma(delta)T-cells constitute very limited portion of the T-cell population. The "gamma(delta) high" group includes sheep, cattle, rabbits, and chicken, in which gamma(delta) T-cells comprise up to 60% of the T-cell population. Here, we compiled TCR-V sequences from the six species and conducted a phylogenetic analysis. We identified various TCR-V gene subgroups based on the analysis. We found that humans and mice have representatives from nearly all of the subgroups identified, while other species have lost subgroups to different extent. Therefore, the gamma(delta) low species have a high degree of diversity of TCR-V genes, while gamma(delta) high species all have limited diversity of TCR-V genes. This pattern is similar to that found for immunoglobulin variable region (IGV) genes.

Fifty-million-year-old polymorphism at an immunoglobulin variable region gene locus in the rabbit evolutionary lineage

The Ig heavy chain variable region (V(H)) genes encode the antigen-binding regions of antibodies. The rabbit genome contains more than 100 V(H) genes, but only one (V(H)1) is preferentially used in the VDJ gene rearrangement. Three highly divergent alleles occur at this V(H)1 locus in most rabbit populations. These three V(H) alleles are also present in snowshoe hare populations, indicating that the polymorphism of the V(H)1 alleles is trans-specific. Here we report the results of a phylogenetic analysis of rabbit Ig germ-line V(H) genes (alleles) together with V(H) genes from humans and mice. We have found that all rabbit V(H) genes belong to one mammalian V(H) group (group C), which also includes various human and mouse V(H) genes. Using the rate of nucleotide substitution obtained from human and mouse V(H) sequences, we have estimated that the V(H)1 polymorphism in the rabbit lineage has been maintained for about 50 million years. This extremely long persistence of V(H)1 polymorphism is apparently caused by overdominant selection, though the real mechanism is unclear.

Locus specificity of polymorphic alleles and evolution by a birth-and-death process in mammalian MHC genes

We have conducted an extensive phylogenetic analysis of polymorphic alleles from human and mouse major histocompatibility complex (MHC) class I and class II genes. The phylogenetic tree obtained for 212 complete human class I allele sequences (HLA-A, -B, and -C) has shown that all alleles from the same locus form a single cluster, which is highly supported by bootstrap values, except for one HLA-B allele (HLA-B*7301). Mouse MHC class I loci did not show locus-specific clusters of polymorphic alleles. This was considered to be because of either interlocus genetic exchange or the confusing designation of loci in different haplotypes at the present time. The locus specificity of polymorphic alleles was also observed in human and mouse MHC class II loci. It was therefore concluded that interlocus recombination or gene conversion is not very important for generating MHC diversity, with a possible exception of mouse class I loci. According to the phylogenetic trees of complete coding sequences, we classified human MHC class I (HLA-A, -B, and -C) and class II (DRB1) alleles into three to five major allelic lineages (groups), which were monophyletic with high bootstrap values. Most of these allelic groups remained unchanged even in phylogenetic trees based on individual exons, though this does not exclude the possibility of intralocus recombination involving short DNA segments. These results, together with the previous observation that MHC loci are subject to frequent duplication and deletion, as well as to balancing selection, indicate that MHC evolution in mammals is in agreement with the birth-and-death model of evolution, rather than with the model of concerted evolution.

The EEG and prognosis in status epilepticus

PURPOSE

To examine the relation between specific EEG features and clinical outcome, determine whether a predictable sequence of EEG patterns exists during status epilepticus (SE), and examine the relation between periodic epileptiform discharges (PEDs) and SE.

METHODS

EEG records of 50 patients with SE admitted to Graduate Hospital between January 1990 and July 1995 were reviewed. Ictal EEGs were available in 72%; 28% had only postictal EEGs. Poor outcome was defined as death or persistent vegetative state, and good outcome as all others. Fisher's Exact test, chi2, and t tests were performed for data analysis.

RESULTS

Of 50 patients, 72% had a good outcome and 28%, a poor outcome. If PEDs were present at any time during or after SE, outcome tended to be worse (p = 0.053). With PEDs, eight (44%) of 18 had a poor outcome; without PEDs, six (19%) of 32 had a poor outcome. Etiologies for SE did not substantially differ in patients with or without PEDs, and structural abnormalities were not more associated with the presence of PEDs. PEDs were seen both early and late, during and after SE. Other EEG characteristics (lateralized vs. bilateral symmetric ictal EEG, discrete vs. continuous ictal activity, and postictal focal slowing) did not relate to outcome. No predictable sequence of EEG changes was found during SE.

CONCLUSIONS

PEDs are the only EEG feature related to outcome in SE and are associated with poor outcome independent of etiology.

The optimization principle in phylogenetic analysis tends to give incorrect topologies when the number of nucleotides or amino acids used is small

In the maximum parsimony (MP) and minimum evolution (ME) methods of phylogenetic inference, evolutionary trees are constructed by searching for the topology that shows the minimum number of mutational changes required (M) and the smallest sum of branch lengths (S), respectively, whereas in the maximum likelihood (ML) method the topology showing the highest maximum likelihood (A) of observing a given data set is chosen. However, the theoretical basis of the optimization principle remains unclear. We therefore examined the relationships of M, S, and A for the MP, ME, and ML trees with those for the true tree by using computer simulation. The results show that M and S are generally greater for the true tree than for the MP and ME trees when the number of nucleotides examined (n) is relatively small, whereas A is generally lower for the true tree than for the ML tree. This finding indicates that the optimization principle tends to give incorrect topologies when n is small. To deal with this disturbing property of the optimization principle, we suggest that more attention should be given to testing the statistical reliability of an estimated tree rather than to finding the optimal tree with excessive efforts. When a reliability test is conducted, simplified MP, ME, and ML algorithms such as the neighbor-joining method generally give conclusions about phylogenetic inference very similar to those obtained by the more extensive tree search algorithms.

A maternal complex partial seizure in labor can affect fetal heart rate

We describe a 33-year-old woman who had a complex partial seizure during labor. Intrauterine pressure catheter and fetal heart monitoring during the seizure revealed a strong, prolonged uterine contraction and simultaneous significant fetal heart rate deceleration for 3.5 minutes. This patient demonstrates that complex partial seizures may result in uterine hyperactivity during labor, which may result in fetal hypoxia.

Positive Darwinian selection after gene duplication in primate ribonuclease genes

Evolutionary mechanisms of origins of new gene function have been a subject of long-standing debate. Here we report a convincing case in which positive Darwinian selection operated at the molecular level during the evolution of novel function by gene duplication. The genes for eosinophil cationic protein (ECP) and eosinophil-derived neurotoxin (EDN) in primates belong to the ribonuclease gene family, and the ECP gene, whose product has an anti-pathogen function not displayed by EDN, was generated by duplication of the EDN gene about 31 million years ago. Using inferred nucleotide sequences of ancestral organisms, we showed that the rate of nonsynonymous nucleotide substitution was significantly higher than that of synonymous substitution for the ECP gene. This strongly suggests that positive Darwinian selection operated in the early stage of evolution of the ECP gene. It was also found that the number of arginine residues increased substantially in a short period of evolutionary time after gene duplication, and these amino acid changes probably produced the novel anti-pathogen function of ECP.

Evolution of immunoglobulin kappa chain variable region genes in vertebrates

The major source of immunoglobulin diversity is variation in DNA sequence among multiple copies of variable region (V) genes of the heavy- and light-chain multigene families. In order to clarify the evolutionary pattern of the multigene family of immunoglobulin light kappa chain V region (V kappa) genes, phylogenetic analyses of V kappa genes from humans and other vertebrate species were conducted. The results obtained indicate that the V kappa genes so far sequenced can be grouped into three major monophyletic clusters, the cartilaginous fish, bony fish and amphibian, and mammalian clusters, and that the cartilaginous fish cluster first separated from the rest of the V kappa genes and then the remaining two clusters diverged. The mammalian V kappa genes can further be divided into 10 V kappa groups, 7 of which are present in the human genome. Human and mouse V kappa genes from different V kappa groups are intermingled rather than clustered on the chromosome, and there are a large number of pseudogenes scattered on the chromosome. This indicates that the chromosomal locations of V kappa genes have been shuffled many times by gene duplication, deletion, and transposition in the evolutionary process and that many genes have become nonfunctional during this process. This mode of evolution is consistent with the model of birth-and-death evolution rather than with the model of concerted evolution. An analysis of duplicate V kappa functional genes and pseudogenes in the human genome has indicated that pseudogenes evolve faster than functional genes but that the rate of nonsynonymous nucleotide substitution in the complementarity-determining regions of V kappa genes has been enhanced by positive Darwinian selection.

Evolution by the birth-and-death process in multigene families of the vertebrate immune system

Concerted evolution is often invoked to explain the diversity and evolution of the multigene families of major histocompatibility complex (MHC) genes and immunoglobulin (Ig) genes. However, this hypothesis has been controversial because the member genes of these families from the same species are not necessarily more closely related to one another than to the genes from different species. To resolve this controversy, we conducted phylogenetic analyses of several multigene families of the MHC and Ig systems. The results show that the evolutionary pattern of these families is quite different from that of concerted evolution but is in agreement with the birth-and-death model of evolution in which new genes are created by repeated gene duplication and some duplicate genes are maintained in the genome for a long time but others are deleted or become nonfunctional by deleterious mutations. We found little evidence that interlocus gene conversion plays an important role in the evolution of MHC and Ig multigene families.

Color vision of ancestral organisms of higher primates

The color vision of mammals is controlled by photosensitive proteins called opsins. Most mammals have dichromatic color vision, but hominoids and Old World (OW) monkeys enjoy trichromatic vision, having the blue-, green-, and red-sensitive opsin genes. Most New World (NW) monkeys are either dichromatic or trichromatic, depending on the sex and genotype. Trichromacy in higher primates is believed to have evolved to facilitate the detection of yellow and red fruits against dappled foliage, but the process of evolutionary change from dichromacy to trichromacy is not well understood. Using the parsimony and the newly developed Bayesian methods, we inferred the amino acid sequences of opsins of ancestral organisms of higher primates. The results suggest that the ancestors of OW and NW monkeys lacked the green gene and that the green gene later evolved from the red gene. The fact that the red/green opsin gene has survived the long nocturnal stage of mammalian evolution and that it is under strong purifying selection in organisms that live in dark environments suggests that this gene has another important function in addition to color vision, probably the control of circadian rhythms.

Accuracies of ancestral amino acid sequences inferred by the parsimony, likelihood, and distance methods

Information about protein sequences of ancestral organisms is important for identifying critical amino acid substitutions that have caused the functional change of proteins in evolution. Using computer simulation, we studied the accuracy of ancestral amino acids inferred by two currently available methods (maximum-parsimony [MP] and maximum-likelihood [ML] methods) in addition to a distance method, which was newly developed in this paper. All three methods give reliable inference when the divergence of amino acid sequences is low. When the extent of sequence divergence is high, however, the ML and distance methods give more accurate results than the MP method, particularly when the phylogenetic tree includes long branches. The accuracy of inferred ancestral amino acids does not change very much when a few present-day sequences are added or eliminated. When an incorrect model of amino acid substitution is used for the ML and distance methods, the accuracy decreases, but it is still higher than that for the MP method. When the tree topology used is partially incorrect, the accuracy in the correct part of the tree is virtually unaffected. The posterior probability of inferred ancestral amino acids computed by the ML and distance methods is an unbiased estimate of the true probability when a correct substitution model is used but may become an overestimate when a simpler model is used.

Phylogenetic analysis in molecular evolutionary genetics

Recent developments of statistical methods in molecular phylogenetics are reviewed. It is shown that the mathematical foundations of these methods are not well established, but computer simulations and empirical data indicate that currently used methods such as neighbor joining, minimum evolution, likelihood, and parsimony methods produce reasonably good phylogenetic trees when a sufficiently large number of nucleotides or amino acids are used. However, when the rate of evolution varies extensively from branch to branch, many methods may fail to recover the true topology. Solid statistical tests for examining the accuracy of trees obtained by neighbor joining, minimum evolution, and least-squares method are available, but the methods for likelihood and parsimony trees are yet to be refined. Parsimony, likelihood, and distance methods can all be used for inferring amino acid sequences of the proteins of ancestral organisms that have become extinct.

Genetic distances and reconstruction of phylogenetic trees from microsatellite DNA

Recently many investigators have used microsatellite DNA loci for studying the evolutionary relationships of closely related populations or species, and some authors proposed new genetic distance measures for this purpose. However, the efficiencies of these distance measures in obtaining the correct tree topology remains unclear. We therefore investigated the probability of obtaining the correct topology (PC) for these new distances as well as traditional distance measures by using computer simulation. We used both the infinite-allele model (IAM) and the stepwise mutation model (SMM), which seem to be appropriate for classical markers and microsatellite loci, respectively. The results show that in both the IAM and SMM CAVALLI-SFORZA and EDWARDS' chord distance (DC) and NEI et al.'s DA distance generally show higher PC values than other distance measures, whether the bottleneck effect exists or not. For estimating evolutionary times, however, NEI's standard distance and GOLDSTEIN et al.'s (delta mu)2 are more appropriate than other distances. Microsatellite DNA seems to be very useful for clarifying the evolutionary relationships of closely related populations.

Central neurogenic hyperventilation in an awake patient with a pontine glioma

A 57-year-old awake man developed central neurogenic hyperventilation associated with a pontine mass. Serum pH reached as high as 7.72 with serum carbon dioxide of 6 torr. Examination of CSF during overbreathing showed that CSF pH was markedly alkaline. Pathologic study showed a well-differentiated pontine astrocytoma. The combination of alkaline CSF and an infiltrating pontine lesion supports a structural, rather than chemical, mechanism for central hyperventilation.

Efficiencies of different genes and different tree-building methods in recovering a known vertebrate phylogeny

The relative efficiencies of different protein-coding genes of the mitochondrial genome and different tree-building methods in recovering a known vertebrate phylogeny (two whale species, cow, rat, mouse, opossum, chicken, frog, and three bony fish species) was evaluated. The tree-building methods examined were the neighbor joining (NJ), minimum evolution (ME), maximum parsimony (MP), and maximum likelihood (ML), and both nucleotide sequences and deduced amino acid sequences were analyzed. Generally speaking, amino acid sequences were better than nucleotide sequences in obtaining the true tree (topology) or trees close to the true tree. However, when only first and second codon positions data were used, nucleotide sequences produced reasonably good trees. Among the 13 genes examined, Nd5 produced the true tree in all tree-building methods or algorithms for both amino acid and nucleotide sequence data. Genes Cytb and Nd4 also produced the correct tree in most tree-building algorithms when amino acid sequence data were used. By contrast, Co2, Nd1, and Nd41 showed a poor performance. In general, large genes produced better results, and when the entire set of genes was used, all tree-building methods generated the true tree. In each tree-building method, several distance measures or algorithms were used, but all these distance measures or algorithms produced essentially the same results. The ME method, in which many different topologies are examined, was no better than the NJ method, which generates a single final tree. Similarly, an ML method, in which many topologies are examined, was no better than the ML star decomposition algorithm that generates a single final tree. In ML the best substitution model chosen by using the Akaike information criterion produced no better results than simpler substitution models. These results question the utility of the currently used optimization principles in phylogenetic construction. Relatively simple methods such as the NJ and ML star decomposition algorithms seem to produce as good results as those obtained by more sophisticated methods. The efficiencies of the NJ, ME, MP, and ML methods in obtaining the correct tree were nearly the same when amino acid sequence data were used. The most important factor in constructing reliable phylogenetic trees seems to be the number of amino acids or nucleotides used.

The root of the phylogenetic tree of human populations

Although African populations have been shown to be most divergent from any other human populations, it has been difficult to establish the root of the phylogenetic tree of human populations since the rate of evolutionary change may vary from population to population owing to the fluctuation of population size and other factors. However, the root can be determined by using the chimpanzee as an outgroup and by employing proper statistical methods. Using this strategy, we constructed phylogenetic trees of human populations for five different sets of gene frequency data. The data sets used were two sets of microsatellite loci data (25 and 8 loci, respectively), restriction fragment length polymorphism (RFLP) data (79 loci), protein polymorphism data (15 loci), and Alu insertion frequency data (4 loci). All these data sets showed that the root is located in the branch connecting African and non-African populations, and in the four data sets the root was established at a significant level. These results indicate that Africans are the first group of people that split from the rest of the human populations.

Evolution of Antennapedia-class homeobox genes

Antennapedia (Antp)-class homeobox genes are involved in the determination of pattern formation along the anterior-posterior axis of the animal embryo. A phylogenetic analysis of Antp-class homeodomains of the nematode, Drosophila, amphioxus, mouse, and human indicates that the 13 cognate group genes of this gene family can be divided into two major groups, i.e., groups I and II. Group I genes can further be divided into subgroups A (cognate groups 1-2), B (cognate group 3), and C (cognate groups 4-8), and group II genes can be divided into subgroups D (cognate groups 9-10) and E (cognate groups 11-13), though this classification is somewhat ambiguous. Evolutionary distances among different amino acid sequences suggest that the divergence between group I and group II genes occurred approximately 1000 million years (MY) ago, and the five different subgroups were formed by approximately 600 MY ago, probably before the divergence of Pseudocoelomates (e.g., nematodes) and Coelomates (e.g., insects and chordates). Our results show that the genes that are phylogenetically close are also closely located in the chromosome, suggesting that the colinearity between the gene expression and gene arrangement was generated by successive tandem gene duplications and that the gene arrangement has been maintained by some sort of selection.

A new method of inference of ancestral nucleotide and amino acid sequences

A statistical method was developed for reconstructing the nucleotide or amino acid sequences of extinct ancestors, given the phylogeny and sequences of the extant species. A model of nucleotide or amino acid substitution was employed to analyze data of the present-day sequences, and maximum likelihood estimates of parameters such as branch lengths were used to compare the posterior probabilities of assignments of character states (nucleotides or amino acids) to interior nodes of the tree; the assignment having the highest probability was the best reconstruction at the site. The lysozyme c sequences of six mammals were analyzed by using the likelihood and parsimony methods. The new likelihood-based method was found to be superior to the parsimony method. The probability that the amino acids for all interior nodes at a site reconstructed by the new method are correct was calculated to be 0.91, 0.86, and 0.73 for all, variable, and parsimony-informative sites, respectively, whereas the corresponding probabilities for the parsimony method were 0.84, 0.76, and 0.51, respectively. The probability that an amino acid in an ancestral sequence is correctly reconstructed by the likelihood analysis ranged from 91.3 to 98.7% for the four ancestral sequences.

Phylogenetic test of the molecular clock and linearized trees

To estimate approximate divergence times of species or species groups with molecular data, we have developed a method of constructing a linearized tree under the assumption of a molecular clock. We present two tests of the molecular clock for a given topology: two-cluster test and branch-length test. The two-cluster test examines the hypothesis of the molecular clock for the two lineages created by an interior node of the tree, whereas the branch-length test examines the deviation of the branch length between the tree root and a tip from the average length. Sequences evolving excessively fast or slow at a high significance level may be eliminated. A linearized tree will then be constructed for a given topology for the remaining sequences under the assumption of rate constancy. We have used these methods to analyze hominoid mitochondrial DNA and drosophilid Adh gene sequences.

Molecular phylogeny and divergence times of drosophilid species

The phylogenetic relationships and divergence times of 39 drosophilid species were studied by using the coding region of the Adh gene. Four genera--Scaptodrosophila, Zaprionus, Drosophila, and Scaptomyza (from Hawaii)--and three Drosophila subgenera--Drosophila, Engiscaptomyza, and Sophophora--were included. After conducting statistical analyses of the nucleotide sequences of the Adh, Adhr (Adh-related gene), and nuclear rRNA genes and a 905-bp segment of mitochondrial DNA, we used Scaptodrosophila as the outgroup. The phylogenetic tree obtained showed that the first major division of drosophilid species occurs between subgenus Sophophora (genus Drosophila) and the group including subgenera Drosophila and Engiscaptomyza plus the genera Zaprionus and Scaptomyza. Subgenus Sophophora is then divided into D. willistoni and the clade of D. obscura and D. melanogaster species groups. In the other major drosophilid group, Zaprionus first separates from the other species, and then D. immigrans leaves the remaining group of species. This remaining group then splits into the D. repleta group and the Hawaiian drosophilid cluster (Hawaiian Drosophila, Engiscaptomyza, and Scaptomyza). Engiscaptomyza and Scaptomyza are tightly clustered. Each of the D. repleta, D. obscura, and D. melanogaster groups is monophyletic. The splitting of subgenera Drosophila and Sophophora apparently occurred about 40 Mya, whereas the D. repleta group and the Hawaiian drosophilid cluster separated about 32 Mya. By contrast, the splitting of Engiscaptomyza and Scaptomyza occurred only about 11 Mya, suggesting that Scaptomyza experienced a rapid morphological evolution. The D. obscura and D. melanogaster groups apparently diverged about 25 Mya. Many of the D. repleta group species studied here have two functional Adh genes (Adh-1 and Adh-2), and these duplicated genes can be explained by two duplication events.

Interior-branch and bootstrap tests of phylogenetic trees

We have compared statistical properties of the interior-branch and bootstrap tests of phylogenetic trees when the neighbor-joining tree-building method is used. For each interior branch of a predetermined topology, the interior-branch and bootstrap tests provide the confidence values, PC and PB, respectively, that indicate the extent of statistical support of the sequence cluster generated by the branch. In phylogenetic analysis these two values are often interpreted in the same way, and if PC and PB are high (say, > or = 0.95), the sequence cluster is regarded as reliable. We have shown that PC is in fact the complement of the P-value used in the standard statistical test, but PB is not. Actually, the bootstrap test usually underestimates the extent of statistical support of species clusters. The relationship between the confidence values obtained by the two tests varies with both the topology and expected branch lengths of the true (model) tree. The most conspicuous difference between PC and PB is observed when the true tree is starlike, and there is a tendency for the difference to increase as the number of sequences in the tree increases. The reason for this is that the bootstrap test tends to become progressively more conservative as the number of sequences in the tree increases. Unlike the bootstrap, the interior-branch test has the same statistical properties irrespective of the number of sequences used when a predetermined tree is considered. Therefore, the interior-branch test appears to be preferable to the bootstrap test as long as unbiased estimators of evolutionary distances are used. However, when the interior-branch is applied to a tree estimated from a given data set, PC may give an overestimate of statistical confidence. For this case, we developed a method for computing a modified version (P'C) of the PC value and showed that this P'C tends to give a conservative estimate of statistical confidence, though it is not as conservative as PB. In this paper we have introduced a model in which evolutionary distances between sequences follow a multivariate normal distribution. This model allowed us to study the relationships between the two tests analytically.