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.

Morphological and molecular heterogeneity within nonmicrosatellite instability-high colorectal cancer

Colorectal cancer (CRC) has traditionally been classified into two groups: microsatellite stable/low-level instability (MSS/MSI-L) and high-level MSI (MSI-H) groups on the basis of multiple molecular and clinicopathologic criteria. Using methylated in tumor (MINT) markers 1, 2,12, and 31, we stratified 77 primary CRCs into three groups: MINT++ (>2), MINT+ (1-2), and MINT- (0 markers methylated). The MSS/MSI-L/ MINT++ group was indistinguishable from the MSI-H/MINT++ group with respect to methylation of p16(INK4a), p14(ARF), and RIZ1, and multiple morphological features. The only significant difference between MSI-H and non-MSI-H MINT++ cancers was the higher frequency of K-ras mutation (P < 0.004) and lower frequency of hMLH1 methylation (P < 0.001) in the latter. These data demonstrate that the separation of CRC into two nonoverlapping groups (MSI-H versus MSS/MSI-L) is a misleading oversimplification.

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.

Read-only-memory-based quantum computation: Experimental explorations using nuclear magnetic resonance and future prospects

Read-only-memory-based (ROM-based) quantum computation (QC) is an alternative to oracle-based QC. It has the advantages of being less magical, and being more suited to implementing space-efficient computation (i.e., computation using the minimum number of writable qubits). Here we consider a number of small (one- and two-qubit) quantum algorithms illustrating different aspects of ROM-based QC. They are: (a) a one-qubit algorithm to solve the Deutsch problem; (b) a one-qubit binary multiplication algorithm; (c) a two-qubit controlled binary multiplication algorithm; and (d) a two-qubit ROM-based version of the Deutsch-Jozsa algorithm. For each algorithm we present experimental verification using nuclear magnetic resonance ensemble QC. The average fidelities for the implementation were in the ranges 0.9-0.97 for the one-qubit algorithms, and 0.84-0.94 for the two-qubit algorithms. We conclude with a discussion of future prospects for ROM-based quantum computation. We propose a four-qubit algorithm, using Grover's iterate, for solving a miniature real-world problem relating to the lengths of paths in a network.

Molecular rates parallel diversification contrasts between carnivorous plant sister lineages

In the carnivorous plant family Lentibulariaceae, the bladderwort lineage (Utricularia and Genlisea) is substantially more species-rich and morphologically divergent than its sister lineage, the butterworts (Pinguicula). Bladderworts have a relaxed body plan that has permitted the evolution of terrestrial, epiphytic, and aquatic forms that capture prey in intricately designed suction bladders or corkscrew-shaped lobster-pot traps. In contrast, the flypaper-trapping butterworts maintain vegetative structures typical of angiosperms. We found that bladderwort genomes evolve significantly faster across seven loci (the trnL intron, the second trnL exon, the trnL-F intergenic spacer, the rps16 intron, rbcL, coxI, and 5.8S rDNA) representing all three genomic compartments. Generation time differences did not show a significant association. We relate these findings to the contested speciation rate hypothesis, which postulates a relationship between increased nucleotide substitution and increased cladogenesis. (C) 2002 The Willi Hennig Society.

Rapid isolation of total RNA and genomic DNA from Hakea actities

Tissues of the Australian native plant species Hakea actities (Proteaceae) contain numerous metabolites and structural compounds that hinder the isolation of nucleic acids. Separate RNA and genomic DNA extraction procedures were developed to isolate high quality nucleic acids from H. actities. Total RNA was extracted from leaves, roots and cluster roots of H. actities grown in low nutrient levels. Cluster root formation in H. actities only occurs when the plants are grown in low nutrient concentrations. However, under these conditions, nucleic acid extraction becomes increasingly difficult. The new procedures are faster than many of the published nucleic acid extraction protocols, and avoid the use of hazardous chemicals. The RNA extraction method was used successfully on another Australian species and a crop species, suggesting that the procedure is useful for molecular studies of a broad range of plants.

Quantum correlated twin atomic beams via photodissociation of a molecular Bose-Einstein condensate

We study the process of photodissociation of a molecular Bose-Einstein condensate as a potential source of strongly correlated twin atomic beams. We show that the two beams can possess nearly perfect quantum squeezing in their relative numbers.

A molecular phylogeny of rainbow skinks (Scincidae : Carlia): taxonomic and biogeographic implications

The phylogenetic relationships amongst 29 species of Carlia and Lygisaurus were estimated using a 726-base-pair segment of the protein-coding mitochondrial ND4 gene. Results do not support the recent resurrection of the genus Lygisaurus. Although most Lygisaurus species formed a single clade, this clade is nested within Carlia and includes Carlia parrhasius. Due to this new molecular evidence, and the paucity of diagnostic morphological characters separating the genera, Lygisaurus de Vis 1884 is re-synonymised with Carlia Gray 1845. Our analysis is also inconsistent with a previous suggestion that Lygisaurus timlowi should be removed to Menetia, a genus that is distantly related relative to outgroups used here. Intraspecific variation in Carlia is, in several instances, greater than interspecific distance. The most strikingly divergent lineages are found within C. rubrigularis, which appears to be paraphyletic, with southern populations more closely related to C. rhomboidalis than to northern populations of C. rubrigularis. The two C. rubrigularis-C. rhomboidalis lineages form part of a major polytomy at an intermediate level of divergence. Lack of resolution at this level, however, does not appear to be due to saturation or loss of phylogenetic signal. Rather, the polytomy probably reflects a period of relatively rapid diversification that occurred sometime during the Miocene.

Magnetic-field-induced superconductivity in layered organic molecular crystals with localized magnetic moments

l-(BETS)2FeCl4 undergoes transitions from an antiferromagnetic insulator to a metal and then to a superconductor as a magnetic field is increased. We use a Hubbard-Kondo model to clarify the role of the Fe31 magnetic ions in these phase transitions. In the high-field regime, the magnetic field acting on the electron spins is compensated by the exchange field He due to the magnetic ions. We show how He can be extracted from the observed splitting of the Shubnikov–de Haas frequencies. We predict the field range for field-induced superconductivity in other materials.

Molecular archaeology

Advances in technologies such as mass spectrometry and capillary electrophoresis have encouraged the study of ancient lipids and other ancient biomolecules. Now, microarray technology looks set to revolutionise the study of ancient DNA, perhaps with as much impact as that of PCR.

From molecular complexes to zwitterions and final products. Reactions between C3O2 and amines

The formation of molecular complexes (prereactive intermediates) between C3O2 and amines (ammonia, dimethylamine, trimethylamine, and 4-(dimethylamino)pyridine) as well as the subsequent transformation of the complexes into C3O2-amine zwitterions in cryogenic matrixes (ca. 40 K) has been observed. In the case of dimethylamine, the formation of tetramethylmalonamide has also been documented. Calculations using density functional theory (B3LYP/6-31G(2d, p)) are used to assign all above species and are in excellent agreement with the IR spectra.

Dopant extraction from scanning capacitance microscopy measurements of p-n junctions using combined inverse modeling and forward simulation

This article proposes a more accurate approach to dopant extraction using combined inverse modeling and forward simulation of scanning capacitance microscopy (SCM) measurements on p-n junctions. The approach takes into account the essential physics of minority carrier response to the SCM probe tip in the presence of lateral electric fields due to a p-n junction. The effects of oxide fixed charge and interface state densities in the grown oxide layer on the p-n junction samples were considered in the proposed method. The extracted metallurgical and electrical junctions were compared to the apparent electrical junction obtained from SCM measurements. (C) 2002 American Institute of Physics.