Search PubMed⌕ Search

SEARCH · Search PubMed

Results for “diversification rate”

Search indexed PubMed citations on genomics, clinical trials, systematic reviews and public health. Explore titles, authors and supplied subject terms, then open the PubMed record.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 235 records · Page 13Linked to original sources
PubMed1993

Population structure in Daphnia obtusa: quantitative genetic and allozymic variation.

Quantitative genetic analyses for body size and for life history characters within and among populations of Daphnia obtusa reveal substantial genetic variance at both hierarchical levels for all traits measured. Simultaneous allozymic analysis on the same population samples indicate a moderate degree of differentiation: GST = 0.28. No associations between electrophoretic genotype and phenotypic characters were found, providing support for the null hypothesis that the allozymic variants are effectively neutral. Therefore, GST can be used as the null hypothesis that neutral phenotypic evolution within populations led to the observed differentiation for the quantitative traits, which I call QST. The results of this study provide evidence that natural selection has promoted diversification for body size among populations, and has impeded diversification for relative fitness. Analyses of population differentiation for clutch size, age at reproduction, and growth rate indicate that neutral phenotypic evolution cannot be excluded as the cause.

Analysis of Variance↗
PubMed2006

The effects of switching behavior on the evolutionary diversification of generalist consumers.

Mathematical models of consumer-resource systems explore the evolution of a morphological trait that determines two resource acquisition rates in a generalist consumer. The consumer also has the ability to adjust its relative consumption of the two resources via behavioral (or developmental) plasticity subject to a trade-off. The analysis examines both stable systems and those with sustained fluctuations in abundance. In both cases, it seeks to determine how the behavioral choice affects the evolution of the morphological characters. The presence of adaptive switching behavior transforms the shape of the relationship between the morphological character and fitness in a manner that usually leads to evolution of two or more morphological types. As in models without switching, the presence of sustained cycles in resource densities often allows the evolution of a generalist as well as two specialists. However, switching expands and shifts the parameter regions yielding this outcome and in some cases allows the evolution and coexistence of at least two generalists as well as the two specialists. This level of diversity supported by only two resources is not seen in the absence of behavioral choice and resource cycles. The results suggest major roles for both behavior and environmental variation in adaptive radiation.

Adaptation, Physiological↗
PubMed2005

Adaptive diversification of vomeronasal receptor 1 genes in rodents.

The vomeronasal receptor 1 (V1R) are believed to be pheromone receptors in rodents. Here we used computational methods to identify 95 and 62 new putative V1R genes from the draft rat and mouse genome sequence, respectively. The rat V1R repertoire consists of 11 subfamilies, 10 of which are shared with the mouse, while rat appears to lack the H and I subfamilies found in mouse and possesses one unique subfamily (M). The estimations of the relative divergence times suggest that many subfamilies originated after the split of rodents and primates. The analysis also reveals that these clusters underwent an expansion very close to the split of mouse and rat. In addition, maximum likelihood analysis showed that the nonsynonymous and synonymous rate ratio for most of these clusters was much higher than one, suggesting the role of positive selection in the diversification of these duplicated V1R genes. Because V1R are thought to mediate the process of signal transduction in response to pheromone detection, we speculate that the V1R genes have evolved under positive Darwinian selection to maintain the ability to discriminate between large and complex pheromonal mixtures.

Adaptation, Physiological↗
PubMed1994

Molecular evolution of the plant R regulatory gene family.

Anthocyanin pigmentation patterns in different plant species are controlled in part by members of the myc-like R regulatory gene family. We have examined the molecular evolution of this gene family in seven plant species. Three regions of the R protein show sequence conservation between monocot and dicot R genes. These regions encode the basic helix-loop-helix domain, as well as conserved N-terminal and C-terminal domains; mean replacement rates for these conserved regions are 1.02 x 10(-9) nonsynonymous nucleotide substitutions per site per year. More than one-half of the protein, however, is diverging rapidly, with nonsynonymous substitution rates of 4.08 x 10(-9) substitutions per site per year. Detailed analysis of R homologs within the grasses (Poaceae) confirm that these variable regions are indeed evolving faster than the flanking conserved domains. Both nucleotide substitutions and small insertion/deletions contribute to the diversification of the variable regions within these regulatory genes. These results demonstrate that large tracts of sequence in these regulatory loci are evolving at a fairly rapid rate.

Amino Acid Sequence↗
PubMed2006

Mitochondrial genome sequence evolution in Chlamydomonas.

The mitochondrial genomes of the Chlorophyta exhibit significant diversity with respect to gene content and genome compactness; however, quantitative data on the rates of nucleotide substitution in mitochondrial DNA, which might help explain the origin of this diversity, are lacking. To gain insight into the evolutionary forces responsible for mitochondrial genome diversification, we sequenced to near completion the mitochondrial genome of the chlorophyte Chlamydomonas incerta, estimated the evolutionary divergence between Chlamydomonas reinhardtii and C. incerta mitochondrial protein-coding genes and rRNA-coding regions, and compared the relative evolutionary rates in mitochondrial and nuclear genes. Synonymous and nonsynonymous substitution rates do not differ significantly between the mitochondrial and nuclear protein-coding genes. The mitochondrial rRNA-coding regions, however, are evolving much faster than their nuclear counterparts, and this difference might be explained by relaxed functional constraints on the mitochondrial translational apparatus due to the small number of proteins synthesized in Chlamydomonas mitochondria. Substitution rates at synonymous sites in a nonstandard mitochondrial gene (rtl) and at intronic and synonymous sites in nuclear genes expressed at low levels suggest that the mutation rate is similar in these two genetic compartments. Potential evolutionary forces shaping mitochondrial genome evolution in Chlamydomonas are discussed.

Animals↗
PubMed2004

Historical climate change and speciation: neotropical seasonally dry forest plants show patterns of both tertiary and quaternary diversification.

Historical climate changes have had a major effect on the distribution and evolution of plant species in the neotropics. What is more controversial is whether relatively recent Pleistocene climatic changes have driven speciation, or whether neotropical species diversity is more ancient. This question is addressed using evolutionary rate analysis of sequence data of nuclear ribosomal internal transcribed spacers in diverse taxa occupying neotropical seasonally dry forests, including Ruprechtia (Polygonaceae), robinioid legumes (Fabaceae), Chaetocalyx and Nissolia (Fabaceae), and Loxopterygium (Anacardiaceae). Species diversifications in these taxa occurred both during and before the Pleistocene in Central America, but were primarily pre-Pleistocene in South America. This indicates plausibility both for models that predict tropical species diversity to be recent and that invoke a role for Pleistocene climatic change, and those that consider it ancient and implicate geological factors such as the Andean orogeny and the closure of the Panama Isthmus. Cladistic vicariance analysis was attempted to identify common factors underlying evolution in these groups. In spite of the similar Mid-Miocene to Pliocene ages of the study taxa, and their high degree of endemism in the different fragments of South American dry forests, the analysis yielded equivocal, non-robust patterns of area relationships.

Base Sequence↗
PubMed2005

Comparative analysis of protein coding sequences from human, mouse and the domesticated pig.

BACKGROUND: The availability of abundant sequence data from key model organisms has made large scale studies of molecular evolution an exciting possibility. Here we use full length cDNA alignments comprising more than 700,000 nucleotides from human, mouse, pig and the Japanese pufferfish Fugu rubrices in order to investigate 1) the relationships between three major lineages of mammals: rodents, artiodactyls and primates, and 2) the rate of evolution and the occurrence of positive Darwinian selection using codon based models of sequence evolution. RESULTS: We provide evidence that the evolutionary splits among primates, rodents and artiodactyls happened shortly after each other, with most gene trees favouring a topology with rodents as outgroup to primates and artiodactyls. Using an unrooted topology of the three mammalian species we show that since their diversification, the pig and mouse lineages have on average experienced 1.44 and 2.86 times as many synonymous substitutions as humans, respectively, whereas the rates of non-synonymous substitutions are more similar. The analysis shows the highest average dN/dS ratio in the human lineage, followed by the pig and then the mouse lineages. Using codon based models we detect signals of positive Darwinian selection in approximately 5.3%, 4.9% and 6.0% of the genes on the human, pig and mouse lineages respectively. Approximately 16.8% of all the genes studied here are not currently annotated as functional genes in humans. Our analyses indicate that a large fraction of these genes may have lost their function quite recently or may still be functional genes in some or all of the three mammalian species. CONCLUSIONS: We present a comparative analysis of protein coding genes from three major mammalian lineages. Our study demonstrates the usefulness of codon-based likelihood models in detecting selection and it illustrates the value of sequencing organisms at different phylogenetic distances for comparative studies.

Animals↗
PubMed1998

Virus dynamics: the effect of target cell limitation and immune responses on virus evolution.

Earlier models of virus evolution during single infections do not include target cell limitation. Here we extend the basic mathematical framework of these theories and study the effect of target cell limitation on the evolution of antigenic variation, increasing replication rates and increasing virus load. We find that target cell limitation provides a selection pressure against antigenic diversification, and can therefore provide a limit to diversity. Antigenic diversity increases virus load; at the maximum level of antigenic diversity, virus load, too, is at a maximum.

Antigenic Variation↗
PubMed1999

An evaluation of elongation factor 1 alpha as a phylogenetic marker for eukaryotes.

Elongation factor 1 alpha (EF-1 alpha) is a highly conserved ubiquitous protein involved in translation that has been suggested to have desirable properties for phylogenetic inference. To examine the utility of EF-1 alpha as a phylogenetic marker for eukaryotes, we studied three properties of EF-1 alpha trees: congruency with other phyogenetic markers, the impact of species sampling, and the degree of substitutional saturation occurring between taxa. Our analyses indicate that the EF-1 alpha tree is congruent with some other molecular phylogenies in identifying both the deepest branches and some recent relationships in the eukaryotic line of descent. However, the topology of the intermediate portion of the EF-1 alpha tree, occupied by most of the protist lineages, differs for different phylogenetic methods, and bootstrap values for branches are low. Most problematic in this region is the failure of all phylogenetic methods to resolve the monophyly of two higher-order protistan taxa, the Ciliophora and the Alveolata. JACKMONO analyses indicated that the impact of species sampling on bootstrap support for most internal nodes of the eukaryotic EF-1 alpha tree is extreme. Furthermore, a comparison of observed versus inferred numbers of substitutions indicates that multiple overlapping substitutions have occurred, especially on the branch separating the Eukaryota from the Archaebacteria, suggesting that the rooting of the eukaryotic tree on the diplomonad lineage should be treated with caution. Overall, these results suggest that the phylogenies obtained from EF-1 alpha are congruent with other molecular phylogenies in recovering the monophyly of groups such as the Metazoa, Fungi, Magnoliophyta, and Euglenozoa. However, the interrelationships between these and other protist lineages are not well resolved. This lack of resolution may result from the combined effects of poor taxonomic sampling, relatively few informative positions, large numbers of overlapping substitutions that obscure phylogenetic signal, and lineage-specific rate increases in the EF-1 alpha data set. It is also consistent with the nearly simultaneous diversification of major eukaryotic lineages implied by the "big-bang" hypothesis of eukaryote evolution.

Algorithms↗
PubMed2005

Phylogenetic analysis of Myriapoda using three nuclear protein-coding genes.

We assessed the ability of three nuclear protein-encoding genes-elongation factor-1alpha (EF-1alpha), RNA polymerase II (Pol II), and elongation factor-2 (EF-2)-from 59 myriapod and 12 non-myriapod species to resolve phylogenetic relationships among myriapod classes and orders. In a previous study using EF-1alpha and Pol II (2134 nt combined) from 34 myriapod taxa, Regier and Shultz recovered widely accepted classes, orders, and families but failed to resolve interclass and interordinal relationships. The result was attributed to heterogenous rates of cladogenesis (specifically, the inability of the slowly evolving sequences to capture phylogenetic signal during rapid phylogenetic diversification) but the possibility of inadequate taxon sampling or limited sequence information could not be excluded. In the present study, the myriapod taxon sample was increased by 25 taxa (73%) and sequence length per taxon was effectively doubled through addition of EF-2 (4318 nt combined). Parsimony and Bayesian analyses of the expanded data set recovered a monophyletic Myriapoda, all four myriapod classes and all multiply sampled orders, often with high node support. However, except for three diplopod clades (Colobognatha, Helminothomorpha, and a subgroup of Pentazonia), few interordinal relationships and no interclass relationships were well supported. These results are similar to those of the earlier study by Regier and Shultz, which indicates that taxon sample and sequence length alone do not readily explain the weakly supported resolution in the earlier study. We review recent paleontological evidence to further develop our proposal that heterogeneity in phylogenetic signal provided by our slowly evolving sequences is due to heterogeneity in the temporal structure of myriapod diversification.

Animals↗
PubMed1998

Retracing the evolution of amino acid specificity in glutaminyl-tRNA synthetase.

Molecular phylogenetic studies of glutaminyl-tRNA synthetase suggest that it has relatively recently evolved from the closely related enzyme glutamyl-tRNA synthetase. We have now attempted to retrace one of the key steps in this process by selecting glutaminyl-tRNA synthetase mutants displaying enhanced glutamic acid recognition. Mutagenesis of two residues proximal to the active site, Phe-90 and Tyr-240, was found to improve glutamic acid recognition 3-5-fold in vitro and resulted in the misacylation of tRNA(Gln) with glutamic acid. In vivo expression of the genes encoding these misacylating variants of glutaminyl-tRNA synthetase reduced cellular growth rates by 40%, probably as a result of an increase in translational error rates. These results provide the first biochemical evidence that glutaminyl-tRNA synthetase originated through duplication and consequent diversification of an ancestral glutamyl-tRNA synthetase-encoding gene.

Amino Acid Sequence↗
PubMed2004

Calibration of galliform molecular clocks using multiple fossils and genetic partitions.

For more than a century, members of the traditional avian order Galliformes (i.e., pheasants, partridges, junglefowl, and relatives) have been among the most intensively studied birds, but still a comprehensive timeframe for their evolutionary history is lacking. Thanks to a number of recent cladistic interpretations for several galliform fossils, candidates now exist that can potentially be used as accurate internal calibrations for molecular clocks. Here, we describe a molecular timescale for Galliformes based on cytochrome b and ND2 using nine mostly internal fossil-based anchorpoints. Beyond application of calibrations spanning the entire evolutionary history of Galliformes, care was taken to investigate the effects of calibration choice, substitution saturation, and rate heterogeneity among lineages on divergence time estimation. Results show broad consistency in time estimation with five out of the nine total calibrations. Our divergence time estimates, based on these anchorpoints, indicate that the early history of Galliformes took place in the Cretaceous, including the origin of the basal-most megapode and perhaps cracid lineages, but that the remaining morphological diversification likely started in the earliest Tertiary. The multi-calibration/multi-genetic partition approach used here highlights the importance of understanding the genetic saturation, variation, and rate constancy spectra for the accurate calculation of divergence times by use of molecular clocks.

Animals↗
PubMed2004

Temporal diversification of Mesoamerican cichlid fishes across a major biogeographic boundary.

The Mexican Neovolcanic Plateau sharply divides the vertebrate fauna of Mesoamerica where the climate of both the neotropics and temperate North America gradually blend. Only a few vertebrate groups such as the Heroine cichlids, distributed from South America to the Rio Grande in North America, are found both north and south of the Neovolcanic Plateau. To better understand the geography and temporal diversification of cichlids at this geologic boundary, we used mitochondrial DNA sequences of the cytochrome b (cyt b) gene to reconstruct the relationships of 52 of the approximately 80 species of Heroine cichlids in Mesoamerica. Our analysis suggests several cichlids in South America should be considered as part of the Mesoamerican Heroine clade because they and the cichlids north of the Isthmus of Panama are clearly supported as monophyletic with respect to all other Neotropical cichlids. We also recovered a group containing species in Paratheraps+Paraneetroplus+Vieja as the sister clade to Herichthys. Herichthys is the only cichlid clade north of the Mexican Plateau and it is monophyletic. Non-parametric rate smoothing of cichlid cyt b sequence resulted in an estimated divergence time of approximately 6 million years for Herichthys. This temporal diversification is concordant with divergence times estimated for anurans in the genus Bufo, a group that exhibits a similar geographic distribution. Our results indicate the 5-million-year-old extension of the Mexican Neovolcanic Plateau to the Gulf Coast of Mexico has strongly influenced the current transition between the vertebrate faunas of the Neotropics and Nearctic.

Animals↗
PubMed2007

Antibodies to the superantigenic site of HIV-1 gp120: hydrolytic and binding activities of the light chain subunit.

Antibodies (Abs) to the superantigenic determinant of HIV gp120 (gp120(SAg)) are potential protective agents against HIV infection. We report that the light chain subunits of Abs cloned from lupus patients using phage library methods bind and hydrolyze gp120(SAg) independent of the heavy chain. Unlike frequent gp120(SAg) recognition by intact Abs attributable to V(H) domain structural elements, the isolated light chains expressed this activity rarely. Four light chains capable of gp120(SAg) recognition were identified by fractionating phage displayed light chains using peptide probes containing gp120 residues 421-433, a gp120(SAg) component. Three light chains expressed non-covalent gp120(SAg) binding and one expressed gp120(SAg) hydrolyzing activity. The hydrolytic light chain was isolated by covalent phage fractionation using an electrophilic analog of residues 421-433. This light chain hydrolyzed a reporter gp120(SAg) substrate and full-length gp120. Other peptide substrates and proteins were hydrolyzed at lower rates or not at all. Consistent with the expected nucleophilic mechanism of hydrolysis, the light chain reacted selectively and covalently with the electrophilic gp120(SAg) peptide analog. The hydrolytic reaction entailed a fast initial step followed by a slower rate limiting step, suggesting rapid substrate acylation and slow deacylation. All four gp120(SAg)-recognizing light chains contained sequence diversifications relative to their germline gene counterparts. These observations indicate that in rare instances, the light chain subunit can bind and hydrolyze gp120(SAg) without the participation of the heavy chain. The pairing of such light chains with heavy chains capable of gp120(SAg) recognition represents a potential mechanism for generating protective Abs with enhanced HIV binding strength and anti-viral proteolytic activity.

Amino Acid Sequence↗
PubMed2002

HLA class II allele and haplotype frequencies in Koreans based on 107 families.

We have investigated the distribution of HLA class II alleles and haplotypes in 107 Korean families (207 parents and 291 children) for the HLA-DRB1, DRB3/B4/B5, DQA1, DQB1 and DPB1 loci. Numbers of alleles observed for each locus were DRB1: 25, DQA1: 14, DQB1: 15, and DPB1: 13. Only two to three alleles were observed for the DRB3 (*0101, *0202, *0301), DRB4 (*0103, * 0103102 N), and DRB5 (*0101, *0102) loci. These alleles showed strong associations with DRB1 alleles: DRB3*0101 with DRB1*1201, *1301 and *1403; DRB3*0301 with DRB1*1202 and *1302; DRB3*0202 with DRB1*0301, *1101, *1401 and *1405; DRB5*0101 and *0102 were exclusively associated with DRB1*1501 and *1502, respectively. The seven most common DRB1-DQB1 haplotypes of frequencies > 0.06 accounted for 52% of the total haplotypes. These haplotypes were exclusively related with the seven most common DRB1-DRB3/B4/B5-DQA1-DQB1 haplotypes: DRB1*1501-DRB5*0101-DQA1*0102-DQB1*0602 (0.085), DRB1*0405-DRB4*0103-DQA1*0303-DQB1*0401 (0.082), DRB1*09012-DRB4*0103-DQA1*0302-DQB1*03032 (0.082), DRB1*0101-DQA1*0101-DQB1*0501 (0.075), DRB1*0701-DRB4*0103-DQA1*0201-DQB1*0202 (0.065), DRB1*0803-DQA1*0103-DQB1*0601 (0.065), and DRB1*1302-DRB3*0301-DQA1*0102-DQB1*0604 (0.065). When these haplotypes were extended to the DPB1 locus, much diversification of haplotypes was observed and only one haplotype remained with a frequency of > 0.06: DRB1*0405-DRB4*0103-DQA1*0303-DQB1*0401-DPB1*0501 (0.062). Such diversification would have resulted from cumulated events of recombination within the HLA class II region, and the actual recombination rate observed between the HLA-DQB1 and DPB1 loci was 2.3% (10/438 informative meioses, including 2 recombinants informative by analysis of TAP genes). Comparison of the distribution of DRB1-DQB1 haplotypes with other populations revealed that Koreans are closest to Japanese people. However, Koreans share a few haplotypes with white people and Africans, which are rare in Japanese: DRB1*0701-DQB1*0202 and DRB1*1302-DQB1*0609. The results obtained in this study will provide useful information for anthropology, organ transplantation and disease association studies.

Alleles↗
PubMed2006

Arrival and diversification of caviomorph rodents and platyrrhine primates in South America.

Platyrrhine primates and caviomorph rodents are clades of mammals that colonized South America during its period of isolation from the other continents, between 100 and 3 million years ago (Mya). Until now, no molecular study investigated the timing of the South American colonization by these two lineages with the same molecular data set. Using sequences from three nuclear genes (ADRA2B, vWF, and IRBP, both separate and combined) from 60 species, and eight fossil calibration constraints, we estimated the times of origin and diversification of platyrrhines and caviomorphs via a Bayesian relaxed molecular clock approach. To account for the possible effect of an accelerated rate of evolution of the IRBP gene along the branch leading to the anthropoids, we performed the datings with and without IRBP (3768 sites and 2469 sites, respectively). The time window for the colonization of South America by primates and by rodents is demarcated by the dates of origin (upper bound) and radiation (lower bound) of platyrrhines and caviomorphs. According to this approach, platyrrhine primates colonized South America between 37.0 +/- 3.0 Mya (or 38.9 +/- 4.0 Mya without IRBP) and 16.8 +/- 2.3 (or 20.1 +/- 3.3) Mya, and caviomorph rodents between 45.4 +/- 4.1 (or 43.7 +/- 4.8) Mya and 36.7 +/- 3.7 (or 35.8 +/- 4.3) Mya. Considering both the fossil record and these molecular datings, the favored scenarios are a trans-Atlantic migration of primates from Africa at the end of the Eocene or beginning of the Oligocene, and a colonization of South America by rodents during the Middle or Late Eocene. Based on our nuclear DNA data, we cannot rule out the possibility of a concomitant arrival of primates and rodents in South America. The caviomorphs radiated soon after their arrival, before the Oligocene glaciations, and these early caviomorph lineages persisted until the present. By contrast, few platyrrhine fossils are known in the Oligocene, and the present-day taxa are the result of a quite recent, Early Miocene diversification.

Animal Migration↗
PubMed2004

Molecular phylogenetics and biogeography of transisthmian and amphi-Atlantic needlefishes (Belonidae: Strongylura and Tylosurus): perspectives on New World marine speciation.

Phylogenetic relationships among New World and eastern Atlantic species in the belonid genera Strongylura and Tylosurus were hypothesized using 3689bp of nucleotide sequence; including the entire mitochondrial (mtDNA) ATP synthase 6 and 8 genes; partial cytochrome b; 12S and 16S ribosomal genes; and introns and exons, 2 and 3 of the nuclear-encoded creatine kinase B gene. Concordant mtDNA and nuclear genealogies permitted well-supported inference of species relationships within Strongylura and Tylosurus, and of the chronology of diversification in the two genera. Our phylogenetic hypothesis permitted an assessment of Rosen's [Syst. Zool. 24 (1975) 431] model of species diversification across the eastern Atlantic to eastern Pacific marine biogeographic track. The spatial predictions of the Rosen model were generally supported, but not the temporal predictions. Furthermore, long branches leading to terminal Belonidae indicated that many species have persisted for millions of years or that nucleotide substitution rates were elevated for some clades. Though heterogeneity of nucleotide substitution rate was indicated across some belonid lineages, molecular clock estimates were used to hypothesize biogeographic scenarios for Strongylura across the eastern Pacific and Atlantic region. Furthermore, use of a molecular clock indicated; that early diversification among contemporary Strongylura may have been initiated by changes in Atlantic Ocean circulation precipitated by closure of the Tethys Sea; and provided approximate dates for the isolation of the freshwater species on the American continents.

ATP Synthetase Complexes↗
PubMed1997

Coordinated amino acid changes in the evolution of mammalian defensins.

The mammalian defensin molecule is a short, highly cationic peptide cytotoxic to both microbial and mammalian cells which is cleaved from a precursor including a signal peptide and a highly anionic propiece. A phylogenetic analysis of 28 complete sequences from five mammalian species (mouse, rat, guinea pig, rabbit, and human) showed species-specific clusters of sequences, indicating that the genes duplicated after divergence of these species. Comparison of rates of synonymous and nonsynonymous nucleotide substitution suggested that gene duplication has often been followed by a period in which diversification of the mature defensins at the amino acid level has been selectively favored. In some comparisons, it appeared that amino acid differences in this region have appeared in a nonrandom fashion so as to change the pattern of residue charges. Because it has been hypothesized that the negative charge in the propiece serves to balance the positive charge in the mature defensin and thus to prevent cytotoxicity prior to cleavage, we used a maximum likelihood method of reconstructing ancestral states in order to test whether this balance has been maintained over evolutionary time in spite of rapid diversification of the mature defensin at the amino acid level. Reconstructed ancestral sequences always maintained a charge balance between mature defensin and propiece, and changes in the net positive charge of the mature defensin were balanced by corresponding changes in the propiece. The results support the hypothesis that, in the evolution of these proteins, amino acid changes have occurred in a coordinated fashion so as to preserve an adaptive phenotype.

Amino Acid Sequence↗