Assessing sequence comparison methods with reliable structurally identified distant evolutionary relationships

Pairwise sequence comparison methods have been assessed using proteins whose relationships are known reliably from their structures and functions, as described in the scop database [Murzin, A. G., Brenner, S. E., Hubbard, T. & Chothia C. (1995) J. Mol. Biol. 247, 536–540]. The evaluation tested the programs blast [Altschul, S. F., Gish, W., Miller, W., Myers, E. W. & Lipman, D. J. (1990). J. Mol. Biol. 215, 403–410], wu-blast2 [Altschul, S. F. & Gish, W. (1996) Methods Enzymol. 266, 460–480], fasta [Pearson, W. R. & Lipman, D. J. (1988) Proc. Natl. Acad. Sci. USA 85, 2444–2448], and ssearch [Smith, T. F. & Waterman, M. S. (1981) J. Mol. Biol. 147, 195–197] and their scoring schemes. The error rate of all algorithms is greatly reduced by using statistical scores to evaluate matches rather than percentage identity or raw scores. The E-value statistical scores of ssearch and fasta are reliable: the number of false positives found in our tests agrees well with the scores reported. However, the P-values reported by blast and wu-blast2 exaggerate significance by orders of magnitude. ssearch, fasta ktup = 1, and wu-blast2 perform best, and they are capable of detecting almost all relationships between proteins whose sequence identities are >30%. For more distantly related proteins, they do much less well; only one-half of the relationships between proteins with 20–30% identity are found. Because many homologs have low sequence similarity, most distant relationships cannot be detected by any pairwise comparison method; however, those which are identified may be used with confidence.

Estimation of evolutionary distances under stationary and nonstationary models of nucleotide substitution

Estimation of evolutionary distances has always been a major issue in the study of molecular evolution because evolutionary distances are required for estimating the rate of evolution in a gene, the divergence dates between genes or organisms, and the relationships among genes or organisms. Other closely related issues are the estimation of the pattern of nucleotide substitution, the estimation of the degree of rate variation among sites in a DNA sequence, and statistical testing of the molecular clock hypothesis. Mathematical treatments of these problems are considerably simplified by the assumption of a stationary process in which the nucleotide compositions of the sequences under study have remained approximately constant over time, and there now exist fairly extensive studies of stationary models of nucleotide substitution, although some problems remain to be solved. Nonstationary models are much more complex, but significant progress has been recently made by the development of the paralinear and LogDet distances. This paper reviews recent studies on the above issues and reports results on correcting the estimation bias of evolutionary distances, the estimation of the pattern of nucleotide substitution, and the estimation of rate variation among the sites in a sequence.

Patterns of variance in stage-structured populations: Evolutionary predictions and ecological implications

Variability in population growth rate is thought to have negative consequences for organism fitness. Theory for matrix population models predicts that variance in population growth rate should be the sum of the variance in each matrix entry times the squared sensitivity term for that matrix entry. I analyzed the stage-specific demography of 30 field populations from 17 published studies for pattern between the variance of a demographic term and its contribution to population growth. There were no instances in which a matrix entry both was highly variable and had a large effect on population growth rate; instead, correlations between estimates of temporal variance in a term and contribution to population growth (sensitivity or elasticity) were overwhelmingly negative. In addition, survivorship or growth sensitivities or elasticities always exceeded those of fecundity, implying that the former two terms always contributed more to population growth rate. These results suggest that variable life history stages tend to contribute relatively little to population growth rates because natural selection may alter life histories to minimize stages with both high sensitivity and high variation.

Substitution rate calibration of small subunit ribosomal RNA identifies chlorarachniophyte endosymbionts as remnants of green algae.

Chlorarachniophytes are amoeboid algae with chlorophyll a and b containing plastids that are surrounded by four membranes instead of two as in plants and green algae. These extra membranes form important support for the hypothesis that chlorarachniophytes have acquired their plastids by the ingestion of another eukaryotic plastid-containing alga. Chlorarachniophytes also contain a small nucleus-like structure called the nucleomorph situated between the two inner and the two outer membranes surrounding the plastid. This nucleomorph is a remnant of the endosymbiont's nucleus and encodes, among other molecules, small subunit ribosomal RNA. Previous phylogenetic analyses on the basis of this molecule provided unexpected and contradictory evidence for the origin of the chlorarachniophyte endosymbiont. We developed a new method for measuring the substitution rates of the individual nucleotides of small subunit ribosomal RNA. From the resulting substitution rate distribution, we derived an equation that gives a more realistic relationship between sequence dissimilarity and evolutionary distance than equations previously available. Phylogenetic trees constructed on the basis of evolutionary distances computed by this new method clearly situate the chlorarachniophyte nucleomorphs among the green algae. Moreover, this relationship is confirmed by transversion analysis of the Chlorarachnion plastid small subunit ribosomal RNA.

Chaperonin-facilitated protein folding: optimization of rate and yield by an iterative annealing mechanism.

We develop a heuristic model for chaperonin-facilitated protein folding, the iterative annealing mechanism, based on theoretical descriptions of "rugged" conformational free energy landscapes for protein folding, and on experimental evidence that (i) folding proceeds by a nucleation mechanism whereby correct and incorrect nucleation lead to fast and slow folding kinetics, respectively, and (ii) chaperonins optimize the rate and yield of protein folding by an active ATP-dependent process. The chaperonins GroEL and GroES catalyze the folding of ribulose bisphosphate carboxylase at a rate proportional to the GroEL concentration. Kinetically trapped folding-incompetent conformers of ribulose bisphosphate carboxylase are converted to the native state in a reaction involving multiple rounds of quantized ATP hydrolysis by GroEL. We propose that chaperonins optimize protein folding by an iterative annealing mechanism; they repeatedly bind kinetically trapped conformers, randomly disrupt their structure, and release them in less folded states, allowing substrate proteins multiple opportunities to find pathways leading to the most thermodynamically stable state. By this mechanism, chaperonins greatly expand the range of environmental conditions in which folding to the native state is possible. We suggest that the development of this device for optimizing protein folding was an early and significant evolutionary event.

The challenge of paleoecological stasis: reassessing sources of evolutionary stability.

The paleontological record of the lower and middle Paleozoic Appalachian foreland basin demonstrates an unprecedented level of ecological and morphological stability on geological time scales. Some 70-80% of fossil morphospecies within assemblages persist in similar relative abundances in coordinated packages lasting as long as 7 million years despite evidence for environmental change and biotic disturbances. These intervals of stability are separated by much shorter periods of ecological and evolutionary change. This pattern appears widespread in the fossil record. Existing concepts of the evolutionary process are unable to explain this uniquely paleontological observation of faunawide coordinated stasis. A principle of evolutionary stability that arises from the ecosystem is explored here. We propose that hierarchical ecosystem theory, when extended to geological time scales, can explain long-term paleoecological stability as the result of ecosystem organization in response to high-frequency disturbance. The accompanying stability of fossil morphologies results from "ecological locking," in which selection is seen as a high-rate response of populations that is hierarchically constrained by lower-rate ecological processes. When disturbance exceeds the capacity of the system, ecological crashes remove these higher-level constraints, and evolution is free to proceed at high rates of directional selection during the organization of a new stable ecological hierarchy.

Evolutionary importance for the membrane enhancement of the production of vitamin D3 in the skin of poikilothermic animals.

The photoproduction of vitamin D in the skin was essential for the evolutionary development of terrestrial vertebrates. During exposure to sunlight, previtamin D3 formed in the skin is isomerized to vitamin D3 (calciol) by a temperature-dependent process. Since early land vertebrates were poikilothermic, the relatively slow conversion of previtamin D3 to vitamin D3 at ambient temperature put them at serious risk for developing vitamin D deficiency, thus leading to a poorly mineralized skeleton that could have ultimately halted further evolutionary development of vertebrates on land. We evaluated the rate of isomerization of previtamin D3 to vitamin D3 in the skin of iguanas and found the isomerization rate was enhanced by 1100% and 1700% at 25 degrees C and 5 degrees C, respectively. It is likely that the membrane entrapment of previtamin D3 in its s-cis,s-cis conformation is responsible for the markedly enhanced conversion of previtamin D3 to vitamin D3. The membrane-enhanced production of vitamin D3 ensures the critical supply of vitamin D3 to poikilothermic animals such as iguanas.

The circumstellar environment and evolutionary status of the supergiant B[e] star Wd1-9

Context. Historically, supergiant (sg)B[e] stars have been difficult to include in theoretical schemes for the evolution of massive OB stars. Aims. The location of Wd1-9 within the coeval starburst cluster Westerlund 1 means that it may be placed into a proper evolutionary context and we therefore aim to utilise a comprehensive multiwavelength dataset to determine its physical properties and consequently its relation to other sgB[e] stars and the global population of massive evolved stars within Wd1. Methods. Multi-epoch R- and I-band VLT/UVES and VLT/FORS2 spectra are used to constrain the properties of the circumstellar gas, while an ISO-SWS spectrum covering 2.45−45μm is used to investigate the distribution, geometry and composition of the dust via a semi-analytic irradiated disk model. Radio emission enables a long term mass-loss history to be determined, while X-ray observations reveal the physical nature of high energy processes within the system. Results. Wd1-9 exhibits the rich optical emission line spectrum that is characteristic of sgB[e] stars. Likewise its mid-IR spectrum resembles those of the LMC sgB[e] stars R66 and 126, revealing the presence of equatorially concentrated silicate dust, with a mass of ~10−4M⊙. Extreme historical and ongoing mass loss (≳ 10−4M⊙yr−1) is inferred from the radio observations. The X-ray properties of Wd1-9 imply the presence of high temperature plasma within the system and are directly comparable to a number of confirmed short-period colliding wind binaries within Wd1. Conclusions. The most complete explanation for the observational properties of Wd1-9 is that it is a massive interacting binary currently undergoing, or recently exited from, rapid Roche-lobe overflow, supporting the hypothesis that binarity mediates the formation of (a subset of) sgB[e] stars. The mass loss rate of Wd1-9 is consistent with such an assertion, while viable progenitor and descendent systems are present within Wd1 and comparable sgB[e] binaries have been identified in the Galaxy. Moreover, the rarity of sgB[e] stars - only two examples are identified from a census of ~ 68 young massive Galactic clusters and associations containing ~ 600 post-Main Sequence stars - is explicable given the rapidity (~ 104yr) expected for this phase of massive binary evolution.

An evolutionary economic perspective on technical change and adjustment in cane harvesting systems in the Australian sugar industry

Australian sugar-producing regions have differed in terms of the extent and rate of incorporation of new technology into harvesting systems. The Mackay sugar industry has lagged behind most other sugar-producing regions in this regard. The reasons for this are addressed by invoking an evolutionary economics perspective. The development of harvesting systems, and the role of technology in shaping them, is mapped and interpreted using the concept of path dependency. Key events in the evolution of harvesting systems are identified, which show how the past has shaped the regional development of harvesting systems. From an evolutionary economics perspective, the outcomes observed are the end result of a specific history.

Do harvest refuges buffer kangaroos against evolutionary responses to selective harvesting?

There is a wealth of literature documenting a directional change of body size in heavily harvested populations. Most of this work concentrates on aquatic systems, but terrestrial populations are equally at risk. This paper explores the capacity of harvest refuges to counteract potential effects of size-selective harvesting on the allele frequency,of populations. We constructed a stochastic, individual-based model parameterized with data on red kangaroos. Because we do not know which part of individual growth would change in the course of natural selection, we explored the effects of two alternative models of individual growth in which alleles affect either the growth rate or the maximum size. The model results show that size-selective harvesting can result in significantly smaller kangaroos for a given age when the entire population is subject to harvesting. In contrast, in scenarios that include dispersal from harvest refuges, the initial allele frequency remains virtually unchanged.

Genetic assignment methods for the direct, real-time estimation of migration rate: a simulation-based exploration of accuracy and power

Genetic assignment methods use genotype likelihoods to draw inference about where individuals were or were not born, potentially allowing direct, real-time estimates of dispersal. We used simulated data sets to test the power and accuracy of Monte Carlo resampling methods in generating statistical thresholds for identifying F-0 immigrants in populations with ongoing gene flow, and hence for providing direct, real-time estimates of migration rates. The identification of accurate critical values required that resampling methods preserved the linkage disequilibrium deriving from recent generations of immigrants and reflected the sampling variance present in the data set being analysed. A novel Monte Carlo resampling method taking into account these aspects was proposed and its efficiency was evaluated. Power and error were relatively insensitive to the frequency assumed for missing alleles. Power to identify F-0 immigrants was improved by using large sample size (up to about 50 individuals) and by sampling all populations from which migrants may have originated. A combination of plotting genotype likelihoods and calculating mean genotype likelihood ratios (D-LR) appeared to be an effective way to predict whether F-0 immigrants could be identified for a particular pair of populations using a given set of markers.

Sociality and the rate of molecular evolution

The molecular clock does not tick at a uniform rate in all taxa but maybe influenced by species characteristics. Eusocial species (those with reproductive division of labor) have been predicted to have faster rates of molecular evolution than their nonsocial relatives because of greatly reduced effective population size; if most individuals in a population are nonreproductive and only one or few queens produce all the offspring, then eusocial animals could have much lower effective population sizes than their solitary relatives, which should increase the rate of substitution of nearly neutral mutations. An earlier study reported faster rates in eusocial honeybees and vespid wasps but failed to correct for phylogenetic nonindependence or to distinguish between potential causes of rate variation. Because sociality has evolved independently in many different lineages, it is possible to conduct a more wide-ranging study to test the generality of the relationship. We have conducted a comparative analysis of 25 phylogenetically independent pairs of social lineages and their nonsocial relatives, including bees, wasps, ants, termites, shrimps, and mole rats, using a range of available DNA sequences (mitochondrial and nuclear DNA coding for proteins and RNAs, and nontranslated sequences). By including a wide range of social taxa, we were able to test whether there is a general influence of sociality on rates of molecular evolution and to test specific predictions of the hypothesis: (1) that social species have faster rates because they have reduced effective population sizes; (2) that mitochondrial genes would show a greater effect of sociality than nuclear genes; and (3) that rates of molecular evolution should be correlated with the degree of sociality. We find no consistent pattern in rates of molecular evolution between social and nonsocial lineages and no evidence that mitochondrial genes show faster rates in social taxa. However, we show that the most highly eusocial Hymenoptera do have faster rates than their nonsocial relatives. We also find that social parasites (that utilize the workers from related species to produce their own offspring) have faster rates than their social relatives, which is consistent with an effect of lower effective population size on rate of molecular evolution. Our results illustrate the importance of allowing for phylogenetic nonindependence when conducting investigations of determinants of variation in rate of molecular evolution.

Coadaptation: A unifying principle in evolutionary thermal biology

Over the last 50 yr, thermal biology has shifted from a largely physiological science to a more integrated science of behavior, physiology, ecology, and evolution. Today, the mechanisms that underlie responses to environmental temperature are being scrutinized at levels ranging from genes to organisms. From these investigations, a theory of thermal adaptation has emerged that describes the evolution of thermoregulation, thermal sensitivity, and thermal acclimation. We review and integrate current models to form a conceptual model of coadaptation. We argue that major advances will require a quantitative theory of coadaptation that predicts which strategies should evolve in specific thermal environments. Simply combining current models, however, is insufficient to understand the responses of organisms to thermal heterogeneity; a theory of coadaptation must also consider the biotic interactions that influence the net benefits of behavioral and physiological strategies. Such a theory will be challenging to develop because each organism's perception of and response to thermal heterogeneity depends on its size, mobility, and life span. Despite the challenges facing thermal biologists, we have never been more pressed to explain the diversity of strategies that organisms use to cope with thermal heterogeneity and to predict the consequences of thermal change for the diversity of communities.

Evolutionary neural networks : models and applications

The scaling problems which afflict attempts to optimise neural networks (NNs) with genetic algorithms (GAs) are disclosed. A novel GA-NN hybrid is introduced, based on the bumptree, a little-used connectionist model. As well as being computationally efficient, the bumptree is shown to be more amenable to genetic coding lthan other NN models. A hierarchical genetic coding scheme is developed for the bumptree and shown to have low redundancy, as well as being complete and closed with respect to the search space. When applied to optimising bumptree architectures for classification problems the GA discovers bumptrees which significantly out-perform those constructed using a standard algorithm. The fields of artificial life, control and robotics are identified as likely application areas for the evolutionary optimisation of NNs. An artificial life case-study is presented and discussed. Experiments are reported which show that the GA-bumptree is able to learn simulated pole balancing and car parking tasks using only limited environmental feedback. A simple modification of the fitness function allows the GA-bumptree to learn mappings which are multi-modal, such as robot arm inverse kinematics. The dynamics of the 'geographic speciation' selection model used by the GA-bumptree are investigated empirically and the convergence profile is introduced as an analytical tool. The relationships between the rate of genetic convergence and the phenomena of speciation, genetic drift and punctuated equilibrium arc discussed. The importance of genetic linkage to GA design is discussed and two new recombination operators arc introduced. The first, linkage mapped crossover (LMX) is shown to be a generalisation of existing crossover operators. LMX provides a new framework for incorporating prior knowledge into GAs.Its adaptive form, ALMX, is shown to be able to infer linkage relationships automatically during genetic search.

Megtakarítási típusok - egy adaptív-evolúciós megközelítés = Types of saving - an adaptive-evolutionary approach

A közgazdaságtanban az ágensalapú modellezés egyik alkalmazási területe a makro ökonómia. Ebben a tanulmányban néhány népszerű megtakarítási szabály létét feltételezve adaptív-evolúciós megközelítésben endogén módon próbálunk következtetni e szabályok relatív életképességére. Három különböző típusú ágenst vezetünk be: egy prudens, egy rövidlátó és egy, a permanensjövedelem-elméletnek megfelelően működőt. Rendkívül erős szelekciós nyomás mellett a prudens típus egyértelműen kiszorítja a másik kettőt. A második legéletképesebbnek a rövidlátó típus tűnik, de már közepes szelekciós nyomásnál sem tűnik el egyik típus sem. Szokásos tőkehatékonyság mellett a prudens típus túlzott beruházási tendenciát visz a gazdaságba, és a gazdaság az aranykori megtakarítási rátánál magasabbat ér el. A hitelkorlátok oldása még nagyobb mértékű túlzott beruházáshoz vezethet, a hitelek mennyiségének növekedése mellett a tőketulajdonosok mintegy "kizsákmányoltatják" magukat azokkal, akiknek nincs tőkejövedelmük. A hosszú távú átlagos fogyasztás szempontjából a három típus kiegyensúlyozott aránya adja a legjobb eredményt, ugyanakkor ez jóval nagyobb ingadozással jár, mint amikor csak prudens típusú háztartások léteznek. ____ Agent-based modelling techniques have been employed for some time in macroeconomics. This paper tests some popular saving rules in an adaptive-evolutionary context of looking at their relative survival values. The three types are prudent, short-sighted, and responsive to the permanent-income hypothesis. It is found that where selection pressure is very high, only the prudent type persists. The second most resilient seems to be the short-sighted type, but all three coexist even at medium levels of selection pressure. When the efficiency of capital approaches the level usually assumed in macroeconomics, the prudent type drives the economy towards excessive accumulation of capital, i. e. a long-term savings rate that exceeds the golden rule. If credit constraints are relaxed, this tendency strengthens as credit grows and capital-owners seem to allow themselves to be exploited" by workers. From the angle of average consumption, the best outcome is obtained from a random distribution of types, although this is accompanied by higher volatility.