Fast evolution strategies

Evolution strategies are a class of general optimisation algorithms which are applicable to functions that are multimodal, nondifferentiable, or even discontinuous. Although recombination operators have been introduced into evolution strategies, the primary search operator is still mutation. Classical evolution strategies rely on Gaussian mutations. A new mutation operator based on the Cauchy distribution is proposed in this paper. It is shown empirically that the new evolution strategy based on Cauchy mutation outperforms the classical evolution strategy on most of the 23 benchmark problems tested in this paper. The paper also shows empirically that changing the order of mutating the objective variables and mutating the strategy parameters does not alter the previous conclusion significantly, and that Cauchy mutations with different scaling parameters still outperform the Gaussian mutation with self-adaptation. However, the advantage of Cauchy mutations disappears when recombination is used in evolution strategies. It is argued that the search step size plays an important role in determining evolution strategies' performance. The large step size of recombination plays a similar role as Cauchy mutation.

An improved seeded region growing algorithm

Recently Adams and Bischof (1994) proposed a novel region growing algorithm for segmenting intensity images. The inputs to the algorithm are the intensity image and a set of seeds - individual points or connected components - that identify the individual regions to be segmented. The algorithm grows these seed regions until all of the image pixels have been assimilated. Unfortunately the algorithm is inherently dependent on the order of pixel processing. This means, for example, that raster order processing and anti-raster order processing do not, in general, lead to the same tessellation. In this paper we propose an improved seeded region growing algorithm that retains the advantages of the Adams and Bischof algorithm fast execution, robust segmentation, and no tuning parameters - but is pixel order independent. (C) 1997 Elsevier Science B.V.

Evolution of a mate recognition system after hybridization between two Drosophila species

I investigated the genetic relationship between male and female components of the mate recognition system and how this relationship influenced the subsequent evolution of the two traits, in a series of replicate populations of interspecific hybrids. Thirty populations of hybrids between Drosophila serrata and Drosophila birchii were established and maintained for 24 generations. At the fifth generation after hybridization, the mating success of hybrid individuals with the D. serrata parent was determined. The genetic correlation between male and female components of the male recognition system, as a consequence of pleiotropy or tight physical linkage, was found to be significant but low (r = 0.388). This result suggested that pleiotropy may play only a minor role in the evolution of mate recognition in this system. At the twenty-fourth generation after hybridization, the mating success of the hybrids was again determined. The evolution of male and female components was investigated by analyzing the direction of evolution of each hybrid line with respect to its initial position in relation to the genetic regression. Male and female components appeared to converge on a single equilibrium point, rather than evolving along trajectories with slope equal to the genetic regression, toward a line of equilibria.

Prediction of MHC class II-binding peptides using an evolutionary algorithm and artificial neural network

Motivation: Prediction methods for identifying binding peptides could minimize the number of peptides required to be synthesized and assayed, and thereby facilitate the identification of potential T-cell epitopes. We developed a bioinformatic method for the prediction of peptide binding to MHC class II molecules. Results: Experimental binding data and expert knowledge of anchor positions and binding motifs were combined with an evolutionary algorithm (EA) and an artificial neural network (ANN): binding data extraction --> peptide alignment --> ANN training and classification. This method, termed PERUN, was implemented for the prediction of peptides that bind to HLA-DR4(B1*0401). The respective positive predictive values of PERUN predictions of high-, moderate-, low- and zero-affinity binder-a were assessed as 0.8, 0.7, 0.5 and 0.8 by cross-validation, and 1.0, 0.8, 0.3 and 0.7 by experimental binding. This illustrates the synergy between experimentation and computer modeling, and its application to the identification of potential immunotheraaeutic peptides.

Saturation mutagenesis of the beta subunit of the human granulocyte-macrophage colony-stimulating factor receptor shows clustering of constitutive mutations, activation of ERK MAP kinase and STAT pathways, and differential beta subunit tyrosine phosphoryl

The high-affinity receptors for human granulocyte-macrophage colony-stimulating factor (GM-CSF), interleukin-1 (IL-3), and IL-5 are heterodimeric complexes consisting of cytokine-specific alpha subunits and a common signal-transducing beta subunit (h beta c). We have previously demonstrated the oncogenic potential of this group of receptors by identifying constitutively activating point mutations in the extracellular and transmembrane domains of h beta c. We report here a comprehensive screen of the entire h beta c molecule that has led to the identification of additional constitutive point mutations by virtue of their ability to confer factor independence on murine FDC-P1 cells. These mutations were clustered exclusively in a central region of h beta c that encompasses the extracellular membrane-proximal domain, transmembrane domain, and membrane-proximal region of the cytoplasmic domain. Interestingly, most h beta c mutants exhibited cell type-specific constitutive activity, with only two transmembrane domain mutants able to confer factor independence on both murine FDC-P1 and BAF-B03 cells. Examination of the biochemical properties of these mutants in FDC-P1 cells indicated that MAP kinase (ERK1/2), STAT, and JAK2 signaling molecules were constitutively activated. In contrast, only some of the mutant beta subunits were constitutively tyrosine phosphorylated. Taken together; these results highlight key regions involved in h beta c activation, dissociate h beta c tyrosine phosphorylation from MAP kinase and STAT activation, and suggest the involvement of distinct mechanisms by which proliferative signals can be generated by h beta c. (C) 1998 by The American Society of Hematology.

Dynamical quantum noise in trapped Bose-Einstein condensates

We introduce the study of dynamical quantum noise in Bose-Einstein condensates through numerical simulation of stochastic partial differential equations obtained using phase-space representations. We derive evolution equations for a single trapped condensate in both the positive-P and Wigner representations and perform simulations to compare the predictions of the two methods. The positive-P approach is found to be highly susceptible to the stability problems that have been observed in other strongly nonlinear, weakly damped systems. Using the Wigner representation, we examine the evolution of several quantities of interest using from a variety of choices of initial stare for the condensate and compare results to those for single-mode models. [S1050-2947(98)06612-8].

mRNA differential display of 2-amino-1-methyl-6-phenylinlidazo[4,5-b]pyridine-induced rat mammary gland tumors

The mRNA differential display technique was used to compare mRNAs between normal mammary gland and turner-derived epithelial cells from female Sprague-Dawley rat mammary gland tumors induced by the heterocyclic amine 2-amino-1-methyl-6-phenylimidazo[4,5-b]pyridine (PhIP) and promoted by a high-fat diet (23.5% corn oil). Two genes, beta-casein and transferrin, were identified as differentially expressed. The expression of these genes was examined across a bank of rat mammary gland tumors derived from animals on a low-fat diet (5% corn oil) or the high-fat diet. Carcinomas had over a 10- and 50-fold lower expression of beta-casein and transferrin, respectively than normal mammary gland. In addition, carcinomas from animals on the high-fat diet showed on average a 5-fold higher expression of beta-casein, and transferrin than carcinomas from animals on the low-fat diet. The results indicate the process of mammary gland tumorigenesis alters the expression of certain genes in the mammary gland, and that the level of dietary fat further modulates the expression of these genes.

Genetic structure and evolution of species in the mangrove genus Avicennia (Avicenniaceae) in the Indo-West Pacific

Allozyme variation in species of the mangrove genus Avicennia was screened in 25 populations collected from 22 locations in the Indo-West Pacific and eastern North America using 11 loci. Several fixed gene differences supported the specific status of Avicennia alba, A. integra, A. marina, and A. rumphiana from the Indo-West Pacific, and A. germinans from the Atlantic-East Pacific. The three varieties of A. marina, var. marina, var. eucalyptifolia, and var. australasica, had higher genetic similarities (Nei's I) and no fixed gene differences, confirming their conspecific status. Strong genetic structuring was observed in A. marina, with sharp changes in gene frequencies at the geographical margins of varietal distributions. The occurrence of alleles found otherwise in only one variety, in only immediately adjacent populations of another variety, provided evidence of introgession between varieties. The varieties appear to have diverged recently in the Pleistocene and are apparently not of ancient Cretaceous origin, as suggested earlier. Despite evidence of high degrees of outcrossing, gene flow among populations was relatively low (N(e)m less than or equal to 1-2), except where populations were geographically continuous, questioning assumptions that these widespread mangrove species achieve high levels of long-distance dispersal.

Increased expression of mitochondrial genes in human alcoholic brain revealed by differential display

Polymerase chain reaction (PCR)-based differential display was used to screen for alterations in gene expression in the mesolimbic system of the human alcoholic brain. Total RNA was extracted from the nucleus accumbens of five alcoholic and five control brains. A selected subpopulation of mRNA was reverse-transcribed to cDNA and amplified by PCR. A differentially expressed cDNA fragment was recovered, cloned, and sequenced. Full sequence analysis of this 467 bp fragment revealed 98.2% homology with the human mitochondrial 12S rRNA gene. Dot-blot analysis showed increased expression of this gem in nucleus accumbens and hippocampus, but not in the superior frontal cortex, primary motor cortex, caudate, and pallidus/putamen In a total of eight human alcoholic brains, compared with seven control brains. A similar increased expression was observed by dot-blot analysis, using RNA from the cerebral cortex of rats chronically treated with alcohol vapor. Hybridization of a 16S rRNA oligonucleotide probe indicated that the expression of both rRNAs genes was significantly increased in nucleus accumbens. These results indicate that chronic alcohol consumption induces alteration in expression of mitochondrial genes in selected brain regions. The altered gene expression may reflect mitochondrial dysfunction In the alcohol-affected brain.

Crystal structure of human T cell leukemia virus type 1 gp21 ectodomain crystallized as a maltose-binding protein chimera reveals structural evolution of retroviral transmembrane proteins

Retroviral entry into cells depends on envelope glycoproteins, whereby receptor binding to the surface-exposed subunit triggers membrane fusion by the transmembrane protein (TM) subunit. We determined the crystal structure at 2.5-Angstrom resolution of the ectodomain of gp21, the TM from human T cell leukemia virus type 1. The gp21 fragment was crystallized as a maltose-binding protein chimera, and the maltose-binding protein domain was used to solve the initial phases by the method of molecular replacement. The structure of gp21 comprises an N-terminal trimeric coiled coil, an adjacent disulfide-bonded loop that stabilizes a chain reversal, and a C-terminal sequence structurally distinct from HIV type 1/simian immunodeficiency virus gp41 that packs against the coil in an extended antiparallel fashion. Comparison of the gp21 structure with the structures of other retroviral TMs contrasts the conserved nature of the coiled coil-forming region and adjacent disulfide-bonded loop with the variable nature of the C-terminal ectodomain segment. The structure points to these features having evolved to enable the dual roles of retroviral TMs: conserved fusion function and an ability to anchor diverse surface-exposed subunit structures to the virion envelope and infected cell surface. The structure of gp21 implies that the N-terminal fusion peptide is in close proximity to the C-terminal transmembrane domain and likely represents a postfusion conformation.

Some new examples for nonuniqueness of the evolution problem of harmonic maps

We find some new examples to show nonuniquence for the heat flow of harmonic maps where weak solutions satisfy the same monotonicity property.

A consistent point-searching algorithm for solution interpolation in unstructured meshes consisting of 4-node bilinear quadrilateral elements

To translate and transfer solution data between two totally different meshes (i.e. mesh 1 and mesh 2), a consistent point-searching algorithm for solution interpolation in unstructured meshes consisting of 4-node bilinear quadrilateral elements is presented in this paper. The proposed algorithm has the following significant advantages: (1) The use of a point-searching strategy allows a point in one mesh to be accurately related to an element (containing this point) in another mesh. Thus, to translate/transfer the solution of any particular point from mesh 2 td mesh 1, only one element in mesh 2 needs to be inversely mapped. This certainly minimizes the number of elements, to which the inverse mapping is applied. In this regard, the present algorithm is very effective and efficient. (2) Analytical solutions to the local co ordinates of any point in a four-node quadrilateral element, which are derived in a rigorous mathematical manner in the context of this paper, make it possible to carry out an inverse mapping process very effectively and efficiently. (3) The use of consistent interpolation enables the interpolated solution to be compatible with an original solution and, therefore guarantees the interpolated solution of extremely high accuracy. After the mathematical formulations of the algorithm are presented, the algorithm is tested and validated through a challenging problem. The related results from the test problem have demonstrated the generality, accuracy, effectiveness, efficiency and robustness of the proposed consistent point-searching algorithm. Copyright (C) 1999 John Wiley & Sons, Ltd.

Rapid evolution towards equal sex ratios in a system with heterogamety

The equal sex ratios found in many species with heterogametic sex determination may be a consequence of selection for equality or the result of the Mendelian segregation of the two sex chromosomes. A lack of genetic variation in sex ratio in species with heterogamety has been the major obstacle in distinguishing between these two hypotheses. We overcome this obstacle by generating hybrids between two species of Drosophila. The resulting hybrid lines had biased sex ratios, allowing us to observe the evolution of sex ratio in replicate populations. Sex ratio converged towards 1:1 after 16 generations of natural selection. These changes in sex ratio were not due to differences in viability between the sexes and the loci underlying the variation in sex ratio were not sex-linked. Equal sex ratios may therefore be the result of natural selection as Fisher predicted.

Electrochemistry of macrocyclic cobalt(III/II) hexaamines: Electrocatalytic hydrogen evolution in aqueous solution

The macrocyclic cobalt hexaamines [Co(trans-diammac)](3+) and [Co(cis-diammac)](3+) (diammac = 6,13-dimethyl-1,4,8,11-tetraazacyclotetradecane-6,13-diamine) are capable of reducing the overpotential for hydrogen evolution on a mercury cathode in aqueous solution. Protons are reduced in a catalytic process involving reoxidation of the Co-II species to its parent Co-III complex. The cycle is robust at neutral pH with no decomposition of catalyst. The stability of the [Co(trans-diammac)](2+) and [Co(cis-diammac)](2+) complexes depends on the pH of the solution and the coordinating properties of the supporting electrolyte. Electrochemical studies indicate that the adsorbed Co-II complex on the surface of mercury is the active catalyst for the reduction of protons to dihydrogen.