Author Correction: Self-organization of a human organizer by combined Wnt and Nodal signalling.

Ref. 7 from Benvenisty and colleagues was inadvertently omitted; this has now been cited in the text and added to the reference list, and subsequent references have been renumbered. The Letter has been corrected online.

Self-organization of a human organizer by combined Wnt and Nodal signalling.

In amniotes, the development of the primitive streak and its accompanying 'organizer' define the first stages of gastrulation. Although these structures have been characterized in detail in model organisms, the human primitive streak and organizer remain a mystery. When stimulated with BMP4, micropatterned colonies of human embryonic stem cells self-organize to generate early embryonic germ layers 1 . Here we show that, in the same type of colonies, Wnt signalling is sufficient to induce a primitive streak, and stimulation with Wnt and Activin is sufficient to induce an organizer, as characterized by embryo-like sharp boundary formation, markers of epithelial-to-mesenchymal transition and expression of the organizer-specific transcription factor GSC. Moreover, when grafted into chick embryos, human stem cell colonies treated with Wnt and Activin induce and contribute autonomously to a secondary axis while inducing a neural fate in the host. This fulfils the most stringent functional criteria for an organizer, and its discovery represents a milestone in human embryology.

Forced periodic expression of G1 cyclins phase-locks the budding yeast cell cycle.

Phase-locking (frequency entrainment) of an oscillator, in which a periodic extrinsic signal drives oscillations at a frequency different from the unperturbed frequency, is a useful property for study of oscillator stability and structure. The cell cycle is frequently described as a biochemical oscillator; however, because this oscillator is tied to key biological events such as DNA replication and segregation, and to cell growth (cell mass increase), it is unclear whether phase locking is possible for the cell cycle oscillator. We found that forced periodic expression of the G(1) cyclin CLN2 phase locks the cell cycle of budding yeast over a range of extrinsic periods in an exponentially growing monolayer culture. We characterize the behavior of cells in a pedigree using a return map to determine the efficiency of entrainment to the externally controlled pulse. We quantify differences between mothers and daughters and how synchronization of an expanding population differs from synchronization of a single oscillator. Mothers only lock intermittently whereas daughters lock completely and in a different period range than mothers. We can explain quantitative features of phase locking in both cell types with an analytically solvable model based on cell size control and how mass is partitioned between mother and daughter cells. A key prediction of this model is that size control can occur not only in G(1), but also later in the cell cycle under the appropriate conditions; this prediction is confirmed in our experimental data. Our results provide quantitative insight into how cell size is integrated with the cell cycle oscillator.

Cross-species comparison of Drosophila male accessory gland protein genes.

Drosophila melanogaster males transfer seminal fluid proteins along with sperm during mating. Among these proteins, ACPs (Accessory gland proteins) from the male's accessory gland induce behavioral, physiological, and life span reduction in mated females and mediate sperm storage and utilization. A previous evolutionary EST screen in D. simulans identified partial cDNAs for 57 new candidate ACPs. Here we report the annotation and confirmation of the corresponding Acp genes in D. melanogaster. Of 57 new candidate Acp genes previously reported in D. melanogaster, 34 conform to our more stringent criteria for encoding putative male accessory gland extracellular proteins, thus bringing the total number of ACPs identified to 52 (34 plus 18 previously identified). This comprehensive set of Acp genes allows us to dissect the patterns of evolutionary change in a suite of proteins from a single male-specific reproductive tissue. We used sequence-based analysis to examine codon bias, gene duplications, and levels of divergence (via dN/dS values and ortholog detection) of the 52 D. melanogaster ACPs in D. simulans, D. yakuba, and D. pseudoobscura. We show that 58% of the 52 D. melanogaster Acp genes are detectable in D. pseudoobscura. Sequence comparisons of ACPs shared and not shared between D. melanogaster and D. pseudoobscura show that there are separate classes undergoing distinctly dissimilar evolutionary dynamics.

Regulatory element detection using correlation with expression.

We present here a new computational method for discovering cis-regulatory elements that circumvents the need to cluster genes based on their expression profiles. Based on a model in which upstream motifs contribute additively to the log-expression level of a gene, this method requires a single genome-wide set of expression ratios and the upstream sequence for each gene, and outputs statistically significant motifs. Analysis of publicly available expression data for Saccharomyces cerevisiae reveals several new putative regulatory elements, some of which plausibly control the early, transient induction of genes during sporulation. Known motifs generally have high statistical significance.

Diffusion in inhomogeneous media: theory and simulations applied to whole cell photobleach recovery.

A continuum description for diffusion in a simple model for an inhomogeneous but isotropic media is derived and implemented numerically. The locally averaged density of diffusible marker is input from experiment to define the sample. Then a single additional parameter, the effective diffusion constant, permits the quantitative simulation of diffusive relaxation from any initial condition. Using this simulation, it is possible to model the recovery of a fluorescently tagged protein in the endoplasmic reticulum (ER) after photobleaching a substantial region of a live cell, and fit an effective diffusion constant which is a property both of the geometry of the ER and the marker. Such quantitative measurements permit inferences about the topology and internal organization of this organelle.

Regulatory element detection using a probabilistic segmentation model.

The availability of genome-wide mRNA expression data for organisms whose genome is fully sequenced provides a unique data set from which to decipher how transcription is regulated by the upstream control region of a gene. A new algorithm is presented which decomposes DNA sequence into the most probable "dictionary" of motifs or words. Identification of words is based on a probabilistic segmentation model in which the significance of longer words is deduced from the frequency of shorter words of various length. This eliminates the need for a separate set of reference data to define probabilities, and genome-wide applications are therefore possible. For the 6,000 upstream regulatory regions in the yeast genome, the 500 strongest motifs from a dictionary of size 1,200 match at a significance level of 15 standard deviations to a database of cis-regulatory elements. Analysis of sets of genes such as those up-regulated during sporulation reveals many new putative regulatory sites in addition to identifying previously known sites.

Building a dictionary for genomes: identification of presumptive regulatory sites by statistical analysis.

The availability of complete genome sequences and mRNA expression data for all genes creates new opportunities and challenges for identifying DNA sequence motifs that control gene expression. An algorithm, "MobyDick," is presented that decomposes a set of DNA sequences into the most probable dictionary of motifs or words. This method is applicable to any set of DNA sequences: for example, all upstream regions in a genome or all genes expressed under certain conditions. Identification of words is based on a probabilistic segmentation model in which the significance of longer words is deduced from the frequency of shorter ones of various lengths, eliminating the need for a separate set of reference data to define probabilities. We have built a dictionary with 1,200 words for the 6, 000 upstream regulatory regions in the yeast genome; the 500 most significant words (some with as few as 10 copies in all of the upstream regions) match 114 of 443 experimentally determined sites (a significance level of 18 standard deviations). When analyzing all of the genes up-regulated during sporulation as a group, we find many motifs in addition to the few previously identified by analyzing the subclusters individually to the expression subclusters. Applying MobyDick to the genes derepressed when the general repressor Tup1 is deleted, we find known as well as putative binding sites for its regulatory partners.

Modeling RNA folding paths with pseudoknots: application to hepatitis delta virus ribozyme.

A quantitative understanding of nucleic acid hybridization is essential to many aspects of biotechnology, such as DNA microarrays, as well as to the structure and folding kinetics of RNA. However, predictions of nucleic acid secondary structures have long been impeded by the presence of helices interior to loops, so-called pseudoknots, which impose complex three-dimensional conformational constraints. In this paper we compute the pseudoknot free energies analytically in terms of known standard parameters, and we show how the results can be included in a kinetic Monte Carlo code to follow the succession of secondary structures during quenched or sequential folding. For the hepatitis delta virus ribozyme, we predict several nonnative stems on the folding path, characterize a kinetically trapped state, interpret several experimentally characterized mutations in terms of the folding path, and suggest how hybridization with other parts of the genome inactivates the newly formed ribozyme.

Kinetic analysis of secretory protein traffic and characterization of golgi to plasma membrane transport intermediates in living cells.

Quantitative time-lapse imaging data of single cells expressing the transmembrane protein, vesicular stomatitis virus ts045 G protein fused to green fluorescent protein (VSVG-GFP), were used for kinetic modeling of protein traffic through the various compartments of the secretory pathway. A series of first order rate laws was sufficient to accurately describe VSVG-GFP transport, and provided compartment residence times and rate constants for transport into and out of the Golgi complex and delivery to the plasma membrane. For ER to Golgi transport the mean rate constant (i.e., the fraction of VSVG-GFP moved per unit of time) was 2.8% per min, for Golgi to plasma membrane transport it was 3.0% per min, and for transport from the plasma membrane to a degradative site it was 0.25% per min. Because these rate constants did not change as the concentration of VSVG-GFP in different compartments went from high (early in the experiment) to low (late in the experiment), secretory transport machinery was never saturated during the experiments. The processes of budding, translocation, and fusion of post-Golgi transport intermediates carrying VSVG- GFP to the plasma membrane were also analyzed using quantitative imaging techniques. Large pleiomorphic tubular structures, rather than small vesicles, were found to be the primary vehicles for Golgi to plasma membrane transport of VSVG-GFP. These structures budded as entire domains from the Golgi complex and underwent dynamic shape changes as they moved along microtubule tracks to the cell periphery. They carried up to 10,000 VSVG-GFP molecules and had a mean life time in COS cells of 3.8 min. In addition, they fused with the plasma membrane without intersecting other membrane transport pathways in the cell. These properties suggest that the post-Golgi intermediates represent a unique transport organelle for conveying large quantities of protein cargo from the Golgi complex directly to the plasma membrane.

Polymer models of meiotic and mitotic chromosomes.

Polymers tied together by constraints exhibit an internal pressure; this idea is used to analyze physical properties of the bottle-brush-like chromosomes of meiotic prophase that consist of polymer-like flexible chromatin loops, attached to a central axis. Using a minimal number of experimental parameters, semiquantitative predictions are made for the bending rigidity, radius, and axial tension of such brushes, and the repulsion acting between brushes whose bristles are forced to overlap. The retraction of lampbrush loops when the nascent transcripts are stripped away, the oval shape of diplotene bivalents between chiasmata, and the rigidity of pachytene chromosomes are all manifestations of chromatin pressure. This two-phase (chromatin plus buffer) picture that suffices for meiotic chromosomes has to be supplemented by a third constituent, a chromatin glue to understand mitotic chromosomes, and explain how condensation can drive the resolution of entanglements. This process resembles a thermal annealing in that a parameter (the affinity of the glue for chromatin and/or the affinity of the chromatin for buffer) has to be tuned to achieve optimal results. Mechanical measurements to characterize this protein-chromatin matrix are proposed. Finally, the propensity for even slightly chemically dissimilar polymers to phase separate (cluster like with like) can explain the apparent segregation of the chromatin into A + T- and G + C-rich regions revealed by chromosome banding.

Driving proteins off DNA using applied tension.

Proteins that bind DNA so as to reduce its end-to-end length can be dissociated by application of force. The thermodynamics of this process are discussed, with special attention to the case of histones bound to DNA (i.e., a string of nucleosomes, or chromatin fiber). The histone octamer is predicted to be driven off chromatin fiber for tensions >2 piconewtons.

Nuclear membrane dynamics and reassembly in living cells: targeting of an inner nuclear membrane protein in interphase and mitosis.

The mechanisms of localization and retention of membrane proteins in the inner nuclear membrane and the fate of this membrane system during mitosis were studied in living cells using the inner nuclear membrane protein, lamin B receptor, fused to green fluorescent protein (LBR-GFP). Photobleaching techniques revealed the majority of LBR-GFP to be completely immobilized in the nuclear envelope (NE) of interphase cells, suggesting a tight binding to heterochromatin and/or lamins. A subpopulation of LBR-GFP within ER membranes, by contrast, was entirely mobile and diffused rapidly and freely (D = 0. 41 +/- 0.1 microm2/s). High resolution confocal time-lapse imaging in mitotic cells revealed LBR-GFP redistributing into the interconnected ER membrane system in prometaphase, exhibiting the same high mobility and diffusion constant as observed in interphase ER membranes. LBR-GFP rapidly diffused across the cell within the membrane network defined by the ER, suggesting the integrity of the ER was maintained in mitosis, with little or no fragmentation and vesiculation. At the end of mitosis, nuclear membrane reformation coincided with immobilization of LBR-GFP in ER elements at contact sites with chromatin. LBR-GFP-containing ER membranes then wrapped around chromatin over the course of 2-3 min, quickly and efficiently compartmentalizing nuclear material. Expansion of the NE followed over the course of 30-80 min. Thus, selective changes in lateral mobility of LBR-GFP within the ER/NE membrane system form the basis for its localization to the inner nuclear membrane during interphase. Such changes, rather than vesiculation mechanisms, also underlie the redistribution of this molecule during NE disassembly and reformation in mitosis.

Electrostatics of lipid bilayer bending.

The electrostatic contribution to spontaneous membrane curvature is calculated within Poisson-Boltzmann theory under a variety of assumptions and emphasizing parameters in the physiological range. Asymmetrical surface charges can be fixed with respect to bilayer midplane area or with respect to the lipid-water area, but induce curvatures of opposite signs. Unequal screening layers on the two sides of a vesicle (e.g., multivalent cationic proteins on one side and monovalent salt on the other) also induce bending. For reasonable parameters, tubules formed by electrostatically induced bending can have radii in the 50-100-nm range, often seen in many intracellular organelles. Thus membrane associated proteins may induce curvature and subsequent budding, without themselves being intrinsically curved. Furthermore, we derive the previously unexplored effects of respecting the strict conservation of charge within the interior of a vesicle. The electrostatic component of the bending modulus is small under most of our conditions and is left as an experimental parameter. The large parameter space of conditions is surveyed in an array of graphs.

Fluctuations and supercoiling of DNA.

Frequently, DNA in vivo is organized into loops that are partially underwound and consequently form interwound helical supercoils. Methods from polymer statistical mechanics are used to show how the competition between entropy (thermal fluctuations) and elastic energy determines supercoil radius and pitch, in good agreement with recent experiments and simulations. Supercoil reorganization by means of slithering (reptation) of the DNA along the supercoil is argued to be a slow process. Extension of supercoiled DNA by an applied force shows a number of unexpected features, including coexistence of interwound and helical states.