The optimization principle in phylogenetic analysis tends to give incorrect topologies when the number of nucleotides or amino acids used is small

In the maximum parsimony (MP) and minimum evolution (ME) methods of phylogenetic inference, evolutionary trees are constructed by searching for the topology that shows the minimum number of mutational changes required (M) and the smallest sum of branch lengths (S), respectively, whereas in the maximum likelihood (ML) method the topology showing the highest maximum likelihood (A) of observing a given data set is chosen. However, the theoretical basis of the optimization principle remains unclear. We therefore examined the relationships of M, S, and A for the MP, ME, and ML trees with those for the true tree by using computer simulation. The results show that M and S are generally greater for the true tree than for the MP and ME trees when the number of nucleotides examined (n) is relatively small, whereas A is generally lower for the true tree than for the ML tree. This finding indicates that the optimization principle tends to give incorrect topologies when n is small. To deal with this disturbing property of the optimization principle, we suggest that more attention should be given to testing the statistical reliability of an estimated tree rather than to finding the optimal tree with excessive efforts. When a reliability test is conducted, simplified MP, ME, and ML algorithms such as the neighbor-joining method generally give conclusions about phylogenetic inference very similar to those obtained by the more extensive tree search algorithms.

On the controllability of distributed systems

To “control” a system is to make it behave (hopefully) according to our “wishes,” in a way compatible with safety and ethics, at the least possible cost. The systems considered here are distributed—i.e., governed (modeled) by partial differential equations (PDEs) of evolution. Our “wish” is to drive the system in a given time, by an adequate choice of the controls, from a given initial state to a final given state, which is the target. If this can be achieved (respectively, if we can reach any “neighborhood” of the target) the system, with the controls at our disposal, is exactly (respectively, approximately) controllable. A very general (and fuzzy) idea is that the more a system is “unstable” (chaotic, turbulent) the “simplest,” or the “cheapest,” it is to achieve exact or approximate controllability. When the PDEs are the Navier–Stokes equations, it leads to conjectures, which are presented and explained. Recent results, reported in this expository paper, essentially prove the conjectures in two space dimensions. In three space dimensions, a large number of new questions arise, some new results support (without proving) the conjectures, such as generic controllability and cases of decrease of cost of control when the instability increases. Short comments are made on models arising in climatology, thermoelasticity, non-Newtonian fluids, and molecular chemistry. The Introduction of the paper and the first part of all sections are not technical. Many open questions are mentioned in the text.

The number of unidentified amacrine cells in the mammalian retina

The three largest known populations of amacrine cells in the rabbit retina were stained with fluorescent probes in whole mounts and counted at a series of retinal eccentricities. The retinas were counterstained using a fluorescent DNA-binding molecule and the total number of nuclei in the inner nuclear layer were counted in confocal sections. From the total number of inner nuclear layer cells and the known fraction of them occupied by amacrine cells, the fraction of amacrine cells made up by the stained populations could be calculated. Starburst cells made up 3%, indoleamine-accumulating cells made up 4%, and AII cells made up 11% of all amacrine cells. By referring four smaller populations of amacrine cells to the number of indoleamine-accumulating cells, they were estimated to make up 4% of all amacrine cells. Thus, 78% of all amacrine cells in the rabbit’s retina are known only from isolated examples, if at all. This proportion is similar in the retinas of the mouse, cat, and monkey. It is likely that a substantial fraction of the local circuit neurons present in other regions of the central nervous system are also invisible as populations to current techniques.

Phase identity of the maize leaf is determined after leaf initiation

The vegetative development of the maize shoot can be divided into juvenile and adult phases based on the types of leaves produced at different times in shoot development. Models for the regulation of phase change make explicit predictions about when the identity of these types of leaves is determined. To test these models, we examined the timing of leaf type determination in maize. Clones induced in transition leaf primordia demonstrated that the juvenile and adult regions of these leaves do not become clonally distinct until after the primordium is 700 μm in length, implying that these cell fates were undetermined at this stage of leaf development. Adult shoot apices were cultured in vitro to induce rejuvenation. We found that leaf primordia as large as 3 mm in length can be at least partially rejuvenated by this treatment, and the location of rejuvenated tissue is correlated with the maturation pattern of the leaf. The amount and distribution of juvenile tissue in rejuvenated leaves suggests that rejuvenation occurs nearly simultaneously in all leaf primordia. In vitro culture rejuvenated existing leaf primordia and the P0 primordium, but did not change the identity of subsequent primordia or the total number of leaves produced by the shoot. This result suggests that leaf identity can be regulated independently of the identity of the shoot apical meristem, and it implies that vegetative phase change is not initiated by a change in the identity of the shoot apical meristem.

Significance analysis of microarrays applied to the ionizing radiation response

Microarrays can measure the expression of thousands of genes to identify changes in expression between different biological states. Methods are needed to determine the significance of these changes while accounting for the enormous number of genes. We describe a method, Significance Analysis of Microarrays (SAM), that assigns a score to each gene on the basis of change in gene expression relative to the standard deviation of repeated measurements. For genes with scores greater than an adjustable threshold, SAM uses permutations of the repeated measurements to estimate the percentage of genes identified by chance, the false discovery rate (FDR). When the transcriptional response of human cells to ionizing radiation was measured by microarrays, SAM identified 34 genes that changed at least 1.5-fold with an estimated FDR of 12%, compared with FDRs of 60 and 84% by using conventional methods of analysis. Of the 34 genes, 19 were involved in cell cycle regulation and 3 in apoptosis. Surprisingly, four nucleotide excision repair genes were induced, suggesting that this repair pathway for UV-damaged DNA might play a previously unrecognized role in repairing DNA damaged by ionizing radiation.

Comprehensive copy number and gene expression profiling of the 17q23 amplicon in human breast cancer

The biological significance of DNA amplification in cancer is thought to be due to the selection of increased expression of a single or few important genes. However, systematic surveys of the copy number and expression of all genes within an amplified region of the genome have not been performed. Here we have used a combination of molecular, genomic, and microarray technologies to identify target genes for 17q23, a common region of amplification in breast cancers with poor prognosis. Construction of a 4-Mb genomic contig made it possible to define two common regions of amplification in breast cancer cell lines. Analysis of 184 primary breast tumors by fluorescence in situ hybridization on tissue microarrays validated these results with the highest amplification frequency (12.5%) observed for the distal region. Based on GeneMap'99 information, 17 known genes and 26 expressed sequence tags were localized to the contig. Analysis of genomic sequence identified 77 additional transcripts. A comprehensive analysis of expression levels of these transcripts in six breast cancer cell lines was carried out by using complementary DNA microarrays. The expression patterns varied from one cell line to another, and several overexpressed genes were identified. Of these, RPS6KB1, MUL, APPBP2, and TRAP240 as well as one uncharacterized expressed sequence tag were located in the two common amplified regions. In summary, comprehensive analysis of the 17q23 amplicon revealed a limited number of highly expressed genes that may contribute to the more aggressive clinical course observed in breast cancer patients with 17q23-amplified tumors.

Long term trends in the evolution of H(3) HA1 human influenza type A

We have studied the HA1 domain of 254 human influenza A(H3N2) virus genes for clues that might help identify characteristics of hemagglutinins (HAs) of circulating strains that are predictive of that strain’s epidemic potential. Our preliminary findings include the following. (i) The most parsimonious tree found requires 1,260 substitutions of which 712 are silent and 548 are replacement substitutions. (ii) The HA1 portion of the HA gene is evolving at a rate of 5.7 nucleotide substitutions/year or 5.7 × 10−3 substitutions/site per year. (iii) The replacement substitutions are distributed randomly across the three positions of the codon when allowance is made for the number of ways each codon can change the encoded amino acid. (iv) The replacement substitutions are not distributed randomly over the branches of the tree, there being 2.2 times more changes per tip branch than for non-tip branches. This result is independent of how the virus was amplified (egg grown or kidney cell grown) prior to sequencing or if sequencing was carried out directly on the original clinical specimen by PCR. (v) These excess changes on the tip branches are probably the result of a bias in the choice of strains to sequence and the detection of deleterious mutations that had not yet been removed by negative selection. (vi) There are six hypervariable codons accumulating replacement substitutions at an average rate that is 7.2 times that of the other varied codons. (vii) The number of variable codons in the trunk branches (the winners of the competitive race against the immune system) is 47 ± 5, significantly fewer than in the twigs (90 ± 7), which in turn is significantly fewer variable codons than in tip branches (175 ± 8). (viii) A minimum of one of every 12 branches has nodes at opposite ends representing viruses that reside on different continents. This is, however, no more than would be expected if one were to randomly reassign the continent of origin of the isolates. (ix) Of 99 codons with at least four mutations, 31 have ratios of non-silent to silent changes with probabilities less than 0.05 of occurring by chance, and 14 of those have probabilities <0.005. These observations strongly support positive Darwinian selection. We suggest that the small number of variable positions along the successful trunk lineage, together with knowledge of the codons that have shown positive selection, may provide clues that permit an improved prediction of which strains will cause epidemics and therefore should be used for vaccine production.

The area code hypothesis revisited: Olfactory receptors and other related transmembrane receptors may function as the last digits in a cell surface code for assembling embryos

Recent evidence emerging from several laboratories, integrated with new data obtained by searching the genome databases, suggests that the area code hypothesis provides a good heuristic model for explaining the remarkable specificity of cell migration and tissue assembly that occurs throughout embryogenesis. The area code hypothesis proposes that cells assemble organisms, including their brains and nervous systems, with the aid of a molecular-addressing code that functions much like the country, area, regional, and local portions of the telephone dialing system. The complexity of the information required to code cells for the construction of entire organisms is so enormous that we assume that the code must make combinatorial use of members of large multigene families. Such a system would reuse the same receptors as molecular digits in various regions of the embryo, thus greatly reducing the total number of genes required. We present the hypothesis that members of the very large families of olfactory receptors and vomeronasal receptors fulfill the criteria proposed for area code molecules and could serve as the last digits in such a code. We discuss our evidence indicating that receptors of these families are expressed in many parts of developing embryos and suggest that they play a key functional role in cell recognition and targeting not only in the olfactory system but also throughout the brain and numerous other organs as they are assembled.

Single-molecule spectroscopy of the β2 adrenergic receptor: Observation of conformational substates in a membrane protein

Single-molecule studies of the conformations of the intact β2 adrenergic receptor were performed in solution. Photon bursts from the fluorescently tagged adrenergic receptor in a micelle were recorded. A photon-burst algorithm and a Poisson time filter were implemented to characterize single molecules diffusing across the probe volume of a confocal microscope. The effects of molecular diffusion and photon number fluctuations were deconvoluted by assuming that Poisson distributions characterize the molecular occupation and photon numbers. Photon-burst size histograms were constructed, from which the source intensity distributions were extracted. Different conformations of the β2 adrenergic receptor cause quenching of the bound fluorophore to different extents and hence produce different photon-burst sizes. An analysis of the photon-burst histograms shows that there are at least two distinct substates for the native adrenergic membrane receptor. This behavior is in contrast to one peak observed for the dye molecule, rhodamine 6G. We test the reliability and robustness of the substate number determination by investigating the application of different binning criteria. Conformational changes associated with agonist binding result in a marked change in the distribution of photon-burst sizes. These studies provide insight into the conformational heterogeneity of G protein-coupled receptors in the presence and absence of a bound agonist.

Preserved neuron number in the hippocampus of aged rats with spatial learning deficits.

Hippocampal neuron loss is widely viewed as a hallmark of normal aging. Moreover, neuronal degeneration is thought to contribute directly to age-related deficits in learning and memory supported by the hippocampus. By taking advantage of improved methods for quantifying neuron number, the present study reports evidence challenging these long-standing concepts. The status of hippocampal-dependent spatial learning was evaluated in young and aged Long-Evans rats using the Morris water maze, and the total number of neurons in the principal cell layers of the dentate gyrus and hippocampus was quantified according to the optical fractionator technique. For each of the hippocampal fields, neuron number was preserved in the aged subjects as a group and in aged individuals with documented learning and memory deficits indicative of hippocampal dysfunction. The findings demonstrate that hippocampal neuronal degeneration is not an inevitable consequence of normal aging and that a loss of principal neurons in the hippocampus fails to account for age-related learning and memory impairment. The observed preservation of neuron number represents an essential foundation for identifying the neurobiological effects of hippocampal aging that account for cognitive decline.

Enhanced green fluorescence by the expression of an Aequorea victoria green fluorescent protein mutant in mono- and dicotyledonous plant cells.

The expression of the jellyfish green fluorescent protein (GFP) in plants was analyzed by transient expression in protoplasts from Nicotiana tabacum, Arabidopsis thaliana, Hordeum vulgare, and Zea mays. Expression of GFP was only observed with a mutated cDNA, from which a recently described cryptic splice site had been removed. However, detectable levels of green fluorescence were only emitted from a small number of protoplasts. Therefore, other mutations in the GFP cDNA leading to single-amino acid exchanges in the chromophore region, which had been previously studied in Escherichia coli, were tested in order to improve the sensitivity of this marker protein. Of the mutations tested so far, the exchange of GFP amino acid tyrosine 66 to histidine (Y66H) led to detection of blue fluorescence in plant protoplasts, while the exchange of amino acid serine 65 to cysteine (S65C) and threonine (S65T) increased the intensity of green fluorescence drastically, thereby significantly raising the detection level for GFP. For GFP S65C, the detectable number of green fluorescing tobacco (BY-2) protoplasts was raised up to 19-fold, while the fluorimetricly determined fluorescence was raised by at least 2 orders of magnitude.