A simulated annealing algorithm for finding consensus sequences

Motivation: A consensus sequence for a family of related sequences is, as the name suggests, a sequence that captures the features common to most members of the family. Consensus sequences are important in various DNA sequencing applications and are a convenient way to characterize a family of molecules. Results: This paper describes a new algorithm for finding a consensus sequence, using the popular optimization method known as simulated annealing. Unlike the conventional approach of finding a consensus sequence by first forming a multiple sequence alignment, this algorithm searches for a sequence that minimises the sum of pairwise distances to each of the input sequences. The resulting consensus sequence can then be used to induce a multiple sequence alignment. The time required by the algorithm scales linearly with the number of input sequences and quadratically with the length of the consensus sequence. We present results demonstrating the high quality of the consensus sequences and alignments produced by the new algorithm. For comparison, we also present similar results obtained using ClustalW. The new algorithm outperforms ClustalW in many cases.

Evolution of colorectal cancer: Change of pace and change of direction

This review compiles evidence for an alternative to the classical adenoma-carcinoma sequence in the evolution of colorectal cancer. It is suggested that between 30 and 50% of colorectal cancers are not initiated by mutation of the tumor suppressor gene APC, but through the epigenetic silencing of genes implicated in the control of differentiation, cell cycle control and DNA repair proficiency. The precursor polyps are often characterized by a serrated architecture, and include hyperplastic polyps, admixed polyps and serrated adenomas. The alternative pathway is heterogeneous and may culminate in cancers showing low or high level DNA microsatellite instability (MSI-L and MSI-H, respectively), and in cancers that are microsatellite stable (MSS). Cancers showing DNA MSI may be characterized by an accelerated evolution. Cancers in hereditary non-polyposis colorectal cancer show features of both classical (adenoma and APC mutation) and alternative pathways (rapid evolution, MSI-H and lack of chromosomal instability). (C) 2001 Blackwell Science Asia Pty Ltd.

Differential expression and distribution of syndecan-1 and-2 in periodontal wound healing of the rat

Cell-surface proteoglycans participate in several biological functions including interactions with adhesion molecules, growth factors and a variety of other effector molecules. Accordingly, these molecules play a central role in various aspects of cell-cell and cell-matrix interactions. To investigate the expression and distribution of the cell surface proteoglycans, syndecan-1 and -2, during periodontal wound healing, immunohistochemical analyses were carried out using monoclonal antibodies against syndecan-1, or -2 core proteins. Both syndecan-1 and -2 were expressed and distributed differentially at various stages of early inflammatory cell infiltration, granulation tissue formation, and tissue remodeling in periodontal wound healing. Expression of syndecan-1 was noted in inflammatory cells within and around the fibrin clots during the earliest stages of inflammatory cell infiltration. During granulation tissue formation it was noted in fibroblast-like cells and newly formed blood vessels. Syndecan-1 was not seen in newly formed bone or cementum matrix at any of the time periods studied. Syndecan-1 expression was generally less during the late stages of wound healing but was markedly expressed in cells that were close to the repairing junctional epithelium. In contrast, syndecan-2 expression and distribution was not evident at the early stages of inflammatory cell infiltration. During the formation of granulation tissue and subsequent tissue remodeling, syndecan-2 was expressed extracellularly in the newly formed fibrils which were oriented toward the root surface. Syndecan-2 was found to be significantly expressed on cells that were close to the root surface and within the matrix of repaired cementum covering root dentin as well as at the alveolar bone edge. These findings indicate that syndecan-1 and -2 may have distinctive functions during wound healing of the periodontium. The appearance of syndecan-1 may involve both cell-cell and cell-matrix interactions, while syndecan-2 showed a predilection to associate with cell-matrix interactions during hard tissue formation.

Interaction between natural and sexual selection during the evolution of mate recognition

The interaction between natural and sexual selection is central to many theories of how mate choice and reproductive isolation evolve, but their joint effect on the evolution of mate recognition has not, to my knowledge, been investigated in an evolutionary experiment. Natural and sexual selection were manipulated in interspecific hybrid populations of Drosophila to determine their effects on the evolution of a mate recognition system comprised of cuticular hydrocarbons (CHCs). The effect of natural selection in isolation indicated that CHCs were costly for males and females to produce. The effect of sexual selection in isolation indicated that females preferred males with a particular CHC composition. However, the interaction between natural and sexual selection had a greater effect on the evolution of the mate recognition system than either process in isolation. When natural and sexual selection were permitted to operate in combination, male CHCs became exaggerated to a greater extent than in the presence of sexual selection alone, and female CHCs evolved against the direction of natural selection. This experiment demonstrated that the interaction between natural and sexual selection is critical in determining the direction and magnitude of the evolutionary response of the mate recognition system.

Molecular clocks in reptiles: Life history influences rate of molecular evolution

Life history has been implicated as a determinant of variation in rate of molecular evolution amongst vertebrate species because of a negative correlation between bode size and substitution rate for many Molecular data sets. Both the generality and the cause of the negative bode size trend have been debated, and the validity of key studies has been questioned (particularly concerning the failure to account for phylogenetic bias). In this study, a comparative method has been used to test for an association between a range of life-history variables-such as body size age at maturity, and clutch size-and DNA substitution rate for three genes (NADH4, cytochrome b, and c-mos). A negative relationship between body size and rate of molecular evolution was found for phylogenetically independent pairs of reptile species spanning turtles. lizards. snakes, crocodile, and tuatara. Although this Study was limited by the number of comparisons for which both sequence and lite-history data were available, the results, suggest that a negative bode size trend in rate of molecular evloution may be a general feature of reptile molecular evolution. consistent with similar studies of mammals and birds. This observation has important implications for uncovering the mechanisms of molecular evolution and warns against assuming that related lineages will share the same substitution rate (a local molecular clock) in order to date evolutionary divergences from DNA sequences.

Evolution and developmental expression of nuclear receptor genes in the ascidian Herdmania

Nuclear receptors are a superfamily of metazoan transcription factors that have been shown to be involved in a wide range of developmental and physiological processes. A PCR-based survey of genomic DNA and developmental cDNAs from the ascidian Herdmania identifies eight members of this multigene family. Sequence comparisons and phylogenetic analyses reveal that these ascidian nuclear receptors are representative of five of the six previously defined nuclear receptor subfamilies and are apparent homologues of retinoic acid [NR1B], retinoid X [NR2B], peroxisome proliferator-activated [NR1C], estrogen related [NR3B], neuron-derived orphan (NOR) [NR4A3], nuclear orphan [NR4A], TR2 orphan [NR2C1] and COUP orphan [NR2F3] receptors. Phylogenetic analyses that include the ascidian genes produce topologically distinct trees that suggest a redefinition of some nuclear receptor subfamilies. These trees also suggest that extensive gene duplication occurred after the vertebrates split from invertebrate chordates. These ascidian nuclear receptor genes are expressed differentially during embryogenesis and metamorphosis.

Adaptive multiscale principal components analysis for online monitoring of wastewater treatment

Fault detection and isolation (FDI) are important steps in the monitoring and supervision of industrial processes. Biological wastewater treatment (WWT) plants are difficult to model, and hence to monitor, because of the complexity of the biological reactions and because plant influent and disturbances are highly variable and/or unmeasured. Multivariate statistical models have been developed for a wide variety of situations over the past few decades, proving successful in many applications. In this paper we develop a new monitoring algorithm based on Principal Components Analysis (PCA). It can be seen equivalently as making Multiscale PCA (MSPCA) adaptive, or as a multiscale decomposition of adaptive PCA. Adaptive Multiscale PCA (AdMSPCA) exploits the changing multivariate relationships between variables at different time-scales. Adaptation of scale PCA models over time permits them to follow the evolution of the process, inputs or disturbances. Performance of AdMSPCA and adaptive PCA on a real WWT data set is compared and contrasted. The most significant difference observed was the ability of AdMSPCA to adapt to a much wider range of changes. This was mainly due to the flexibility afforded by allowing each scale model to adapt whenever it did not signal an abnormal event at that scale. Relative detection speeds were examined only summarily, but seemed to depend on the characteristics of the faults/disturbances. The results of the algorithms were similar for sudden changes, but AdMSPCA appeared more sensitive to slower changes.

Differential requirements for COPI coats in formation of replication complexes among three genera of Picornaviridae

Picornavirus RNA replication requires the formation of replication complexes (RCs). consisting of virus-induced vesicles associated with viral nonstructural proteins and RNA. Brefeldin A (BFA) has been shown to strongly inhibit RNA replication of poliovirus but not of encephalomyocarditis virus (EMCV). Here, we demonstrate that the replication of parechovirus 1 (ParV1) is partly resistant to BFA, whereas echovirus 11 (EV11) replication is strongly inhibited. Since BFA inhibits COPI-dependent steps in endoplasmic reticulum (ER)-Golgi transport, we tested a hypothesis that different picornaviruses may have differential requirements for COPI in the formation of their RCs. Using immunofluorescence and cryo-immunoelectron microscopy we examined the association of a COPI component, beta-COP, with the RCs of EMCV, ParV1, and EV11 EMCV RCs did not contain beta-COP. In contrast, beta-COP appeared to be specifically distributed to the RCs of EV11 In ParV1-infected cells beta-COP was largely dispersed throughout the cytoplasm, with some being present in the RCs. These results suggest that there are differences in the involvement of COPI in the formation of the RCs of various picornaviruses, corresponding to their differential sensitivity to BFA. EMCV RCs are likely to be formed immediately after vesicle budding from the ER, prior to COPI association with membranes. ParV1 RCs are formed from COPI-containing membranes but COPI is unlikely to be directly involved in their formation, whereas formation of EV11 RCs appears to be dependent on COPI association with membranes.

Cash, customers, and care: The experience and meaning of differential payment for high care places in aged care facilities

The introduction of new asset/income tested charges for high care residents was the 1997-98 Commonwealth government policy response to concerns about financing residential aged care. This in-depth study of residents, families, staff and managers in three aged care facilities explores issues of equity, access and empowerment arising when some residents pay more for the same level of care and amenity. The study reports little evidence of financial contributions affecting access to high care places and the delivery of care, the potential for differential access to amenities such as single rooms linked to the extra payments, and no evidence of a sense of empowerment linked to payment of the new charges. The complexity of current financial arrangements, access to appropriate financial advice at the time of entry, and the potential for an informal two tier system in relation to the allocation Of amenities are identified as developing policy issues.

Global distribution and genetic discontinuities of mangroves - emerging patterns in the evolution of Rhizophora

Mangroves are often described as a group of plants with common features and common origins based mostly on their broad distributional patterns, together with an erroneous view of comparable abilities in long-distance dispersal. However, whilst mangroves have common needs to adapt to rigorous environmental constraints associated with regular seawater inundation, individual taxa have developed different strategies and characteristics. Since mangroves are a genetically diverse group of mostly flowering plants, they may also have evolved at quite different geological periods, dispersed at different rates from different locations and developed different adaptive strategies. Current distributions of individual taxa show numerous instances of unusual extant distribution which demonstrate finite dispersal limitations, especially across open water. Our preliminary assessment of broad distribution and discontinuities reveals important patterns. Discontinuities, in the absence of current dispersal barriers, may be explained by persistent past barriers. As we learn more about discontinuities, we are beginning to appreciate their immense implications and what they might tell us about past geological conditions and how these might have influenced the distribution and evolution of mangroves. In this article, we describe emerging patterns in genetic relationships and distributions based on both current knowledge and preliminary results of our studies of molecular and morphometric characteristics of Rhizophora species in the Indo West Pacific region.

Boundary value problems for systems of difference equations associated with systems of second-order ordinary differential equations

We study difference equations which arise as discrete approximations to two-point boundary value problems for systems of second-order ordinary differential equations. We formulate conditions which guarantee a priori bounds on first differences of solutions to the discretized problem. We establish existence results for solutions to the discretized boundary value problems subject to nonlinear boundary conditions. We apply our results to show that solutions to the discrete problem converge to solutions of the continuous problem in an aggregate sense. (C) 2002 Elsevier Science Ltd. All rights reserved.

Difference equations associated with fully nonlinear boundary value problems for second order ordinary differential equations

We study the continuous problem y"=f(x,y,y'), xc[0,1], 0=G((y(0),y(1)),(y'(0), y'(1))), and its discrete approximation (y(k+1)-2y(k)+y(k-1))/h(2) =f(t(k), y(k), v(k)), k = 1,..., n-1, 0 = G((y(0), y(n)), (v(1), v(n))), where f and G = (g(0), g(1)) are continuous and fully nonlinear, h = 1/n, v(k) = (y(k) - y(k-1))/h, for k =1,..., n, and t(k) = kh, for k = 0,...,n. We assume there exist strict lower and strict upper solutions and impose additional conditions on f and G which are known to yield a priori bounds on, and to guarantee the existence of solutions of the continuous problem. We show that the discrete approximation also has solutions which approximate solutions of the continuous problem and converge to the solution of the continuous problem when it is unique, as the grid size goes to 0. Homotopy methods can be used to compute the solution of the discrete approximation. Our results were motivated by those of Gaines.

The trematodes of groupers (Serranidae : Epinephelinae) knowledge, nature and evolution

Groupers (Epinephelinae) are prominent marine fishes distributed in the warmer waters of the world. Review of the literature suggests that trematodes are known from only 62 of the 159 species and only 9 of 15 genera; nearly 90% of host-parasite combinations have been reported only once or twice. All 20 families and all but 7 of 76 genera of trematodes found in epinephelines also occur in non-epihephelines. Only 12 genera of trematodes are reported from both the Atlantic-Eastern Pacific and the Indo-West Pacific. Few (perhaps no) species are credibly cosmopolitan but some have wide distributions across the Indo-West Pacific. The hierarchical 'relatedness' of epinephelines as suggested by how they share trematode taxa (families, genera, species) shows little congruence with what is known of their phylogeny. The major determinant of relatedness appears to be geographical proximity. Together these attributes suggest that host-parasite coevolution has contributed little to the evolution of trematode communities of epinephelines. Instead, they appear to have arisen through localized episodes of host-switching, presumably both into and out of the epinephelines. The Epinephelinae may well be typical of most groups of marine fishes both in the extent to which their trematode parasites are known and in that, apparently, co-evolution has contributed little to the evolution of their communities of trematodes.