@article{1774, abstract = {A number of superconducting qubits, such as the transmon or the phase qubit, have an energy level structure with small anharmonicity. This allows for convenient access of higher excited states with similar frequencies. However, special care has to be taken to avoid unwanted higher-level populations when using short control pulses. Here we demonstrate the preparation of arbitrary three level superposition states using optimal control techniques in a transmon. Performing dispersive readout, we extract the populations of all three levels of the qutrit and study the coherence of its excited states. Finally we demonstrate full quantum state tomography of the prepared qutrit states and evaluate the fidelities of a set of states, finding on average 95%.}, author = {Bianchetti, R and Filipp, Stefan and Baur, Matthias P and Johannes Fink and Lang, C and Steffen, L. Kraig and Boissonneault, Maxime and Blais, Alexandre and Wallraff, Andreas}, journal = {Physical Review Letters}, number = {22}, publisher = {American Physical Society}, title = {{Control and tomography of a three level superconducting artificial atom}}, doi = {10.1103/PhysRevLett.105.223601}, volume = {105}, year = {2010}, } @article{7703, abstract = {By combining gene expression profiling with image registration, Tomer et al. (2010) find that the mushroom body of the segmented worm Platynereis dumerilii shares many features with the mammalian cerebral cortex. The authors propose that the mushroom body and cortex evolved from the same structure in the common ancestor of vertebrates and invertebrates.}, author = {Sweeney, Lora Beatrice Jaeger and Luo, Liqun}, issn = {0092-8674}, journal = {Cell}, number = {5}, pages = {679--681}, publisher = {Elsevier}, title = {{‘Fore brain: A hint of the ancestral cortex}}, doi = {10.1016/j.cell.2010.08.024}, volume = {142}, year = {2010}, } @article{8472, abstract = {Characterization of protein dynamics by solid-state NMR spectroscopy requires robust and accurate measurement protocols, which are not yet fully developed. In this study, we investigate the backbone dynamics of microcrystalline ubiquitin using different approaches. A rotational-echo double-resonance type (REDOR-type) methodology allows one to accurately measure 1H−15N order parameters in highly deuterated samples. We show that the systematic errors in the REDOR experiment are as low as 1% or even less, giving access to accurate data for the amplitudes of backbone mobility. Combining such dipolar-coupling-derived order parameters with autocorrelated and cross-correlated 15N relaxation rates, we are able to quantitate amplitudes and correlation times of backbone dynamics on picosecond and nanosecond time scales in a residue-resolved manner. While the mobility on picosecond time scales appears to have rather uniform amplitude throughout the protein, we unambiguously identify and quantitate nanosecond mobility with order parameters S2 as low as 0.8 in some regions of the protein, where nanosecond dynamics has also been revealed in solution state. The methodology used here, a combination of accurate dipolar-coupling measurements and different relaxation parameters, yields details about dynamics on different time scales and can be applied to solid protein samples such as amyloid fibrils or membrane proteins.}, author = {Schanda, Paul and Meier, Beat H. and Ernst, Matthias}, issn = {0002-7863}, journal = {Journal of the American Chemical Society}, number = {45}, pages = {15957--15967}, publisher = {American Chemical Society}, title = {{Quantitative analysis of protein backbone dynamics in microcrystalline ubiquitin by solid-state NMR spectroscopy}}, doi = {10.1021/ja100726a}, volume = {132}, year = {2010}, } @article{8473, abstract = {β2-microglobulin (β2m), the light chain of class I major histocompatibility complex, is responsible for the dialysis-related amyloidosis and, in patients undergoing long term dialysis, the full-length and chemically unmodified β2m converts into amyloid fibrils. The protein, belonging to the immunoglobulin superfamily, in common to other members of this family, experiences during its folding a long-lived intermediate associated to the trans-to-cis isomerization of Pro-32 that has been addressed as the precursor of the amyloid fibril formation. In this respect, previous studies on the W60G β2m mutant, showing that the lack of Trp-60 prevents fibril formation in mild aggregating condition, prompted us to reinvestigate the refolding kinetics of wild type and W60G β2m at atomic resolution by real-time NMR. The analysis, conducted at ambient temperature by the band selective flip angle short transient real-time two-dimensional NMR techniques and probing the β2m states every 15 s, revealed a more complex folding energy landscape than previously reported for wild type β2m, involving more than a single intermediate species, and shedding new light into the fibrillogenic pathway. Moreover, a significant difference in the kinetic scheme previously characterized by optical spectroscopic methods was discovered for the W60G β2m mutant.}, author = {Corazza, Alessandra and Rennella, Enrico and Schanda, Paul and Mimmi, Maria Chiara and Cutuil, Thomas and Raimondi, Sara and Giorgetti, Sofia and Fogolari, Federico and Viglino, Paolo and Frydman, Lucio and Gal, Maayan and Bellotti, Vittorio and Brutscher, Bernhard and Esposito, Gennaro}, issn = {0021-9258}, journal = {Journal of Biological Chemistry}, keywords = {Cell Biology, Biochemistry, Molecular Biology}, number = {8}, pages = {5827--5835}, publisher = {American Society for Biochemistry & Molecular Biology}, title = {{Native-unlike long-lived intermediates along the folding pathway of the amyloidogenic protein β2-Microglobulin revealed by real-time two-dimensional NMR}}, doi = {10.1074/jbc.m109.061168}, volume = {285}, year = {2010}, } @inbook{8506, author = {Hunt, Brian R. and Kaloshin, Vadim}, booktitle = {Handbook of Dynamical Systems}, isbn = {9780444531414}, issn = {1874-575X}, pages = {43--87}, publisher = {Elsevier}, title = {{Prevalence}}, doi = {10.1016/s1874-575x(10)00310-3}, volume = {3}, year = {2010}, } @inproceedings{8507, abstract = {We study a Cr nearly integrable Hamiltonian system defined on 𝕋3 × ℝ3. Let and µΣ1 be the restriction of Lebesgue measure on 𝕋3 × ℝ3 to ∑. We prove there is a perturbation , and an orbit (q(t), p(t)): ℝ → 𝕋3 × ℝ3 of the Hamiltonian equation such that .}, author = {Kaloshin, Vadim and ZHANG, KE and ZHENG, YONG}, booktitle = {XVIth International Congress on Mathematical Physics}, isbn = {9789814304627}, location = {Prague, Czech Republic}, pages = {314--322}, publisher = {World Scientific}, title = {{Almost dense orbit on energy surface}}, doi = {10.1142/9789814304634_0017}, year = {2010}, } @article{857, abstract = {The need to maintain the structural and functional integrity of an evolving protein severely restricts the repertoire of acceptable amino-acid substitutions. However, it is not known whether these restrictions impose a global limit on how far homologous protein sequences can diverge from each other. Here we explore the limits of protein evolution using sequence divergence data. We formulate a computational approach to study the rate of divergence of distant protein sequences and measure this rate for ancient proteins, those that were present in the last universal common ancestor. We show that ancient proteins are still diverging from each other, indicating an ongoing expansion of the protein sequence universe. The slow rate of this divergence is imposed by the sparseness of functional protein sequences in sequence space and the ruggedness of the protein fitness landscape: 98 per cent of sites cannot accept an amino-acid substitution at any given moment but a vast majority of all sites may eventually be permitted to evolve when other, compensatory, changes occur. Thus, 3.5 × 10 9 yr has not been enough to reach the limit of divergent evolution of proteins, and for most proteins the limit of sequence similarity imposed by common function may not exceed that of random sequences.}, author = {Povolotskaya, Inna and Fyodor Kondrashov}, journal = {Nature}, number = {7300}, pages = {922 -- 926}, publisher = {Nature Publishing Group}, title = {{Sequence space and the ongoing expansion of the protein universe}}, doi = {10.1038/nature09105}, volume = {465}, year = {2010}, } @article{862, abstract = {A long-standing controversy in evolutionary biology is whether or not evolving lineages can cross valleys on the fitness landscape that correspond to low-fitness genotypes, which can eventually enable them to reach isolated fitness peaks1-9. Here we study the fitness landscapes traversed by switches between different AU and GC Watson-Crick nucleotide pairs at complementary sites of mitochondrial transfer RNA stem regions in 83 mammalian species. We find that such Watson-Crick switches occur 30-40 times more slowly than pairs of neutral substitutions, and that alleles corresponding to GU and AC non-Watson-Crick intermediate states segregate within human populations at low frequencies, similar to those of non-synonymous alleles. Substitutions leading to a Watson-Crick switch are strongly correlated, especially in mitochondrial tRNAs encoded on the GT-nucleotide-rich strand of the mitochondrial genome. Using these data we estimate that a typical Watson-Crick switch involves crossing a fitness valley of a depth of about 10-3 or even about 10-2, with AC intermediates being slightly more deleterious than GU intermediates. This compensatory evolution must proceed through rare intermediate variants that never reach fixation. The ubiquitous nature of compensatory evolution in mammalian mitochondrial tRNAs and other molecules implies that simultaneous fixation of two alleles that are individually deleterious may be a common phenomenon at the molecular level.}, author = {Meer, Margarita V and Kondrashov, Alexey S and Artzy-Randrup, Yael and Fyodor Kondrashov}, journal = {Nature}, number = {7286}, pages = {279 -- 282}, publisher = {Nature Publishing Group}, title = {{Compensatory evolution in mitochondrial tRNAs navigates valleys of low fitness}}, doi = {10.1038/nature08691}, volume = {464}, year = {2010}, } @article{872, abstract = {The rate of spontaneous mutation in natural populations is a fundamental parameter for many evolutionary phenomena. Because the rate of mutation is generally low, most of what is currently known about mutation has been obtained through indirect, complex and imprecise methodological approaches. However, in the past few years genome-wide sequencing of closely related individuals has made it possible to estimate the rates of mutation directly at the level of the DNA, avoiding most of the problems associated with using indirect methods. Here, we review the methods used in the past with an emphasis on next generation sequencing, which may soon make the accurate measurement of spontaneous mutation rates a matter of routine.}, author = {Fyodor Kondrashov and Kondrashov, Alexey S}, journal = {Philosophical Transactions of the Royal Society of London. Series B, Biological Sciences}, number = {1544}, pages = {1169 -- 1176}, publisher = {Royal Society, The}, title = {{Measurements of spontaneous rates of mutations in the recent past and the near future}}, doi = {10.1098/rstb.2009.0286}, volume = {365}, year = {2010}, } @article{884, abstract = {Background: Divergence of two independently evolving sequences that originated from a common ancestor can be described by two parameters, the asymptotic level of divergence E and the rate r at which this level of divergence is approached. Constant negative selection impedes allele replacements and, therefore, is routinely assumed to decelerate sequence divergence. However, its impact on E and on r has not been formally investigated.Results: Strong selection that favors only one allele can make E arbitrarily small and r arbitrarily large. In contrast, in the case of 4 possible alleles and equal mutation rates, the lowest value of r, attained when two alleles confer equal fitnesses and the other two are strongly deleterious, is only two times lower than its value under selective neutrality.Conclusions: Constant selection can strongly constrain the level of sequence divergence, but cannot reduce substantially the rate at which this level is approached. In particular, under any constant selection the divergence of sequences that accumulated one substitution per neutral site since their origin from the common ancestor must already constitute at least one half of the asymptotic divergence at sites under such selection.Reviewers: This article was reviewed by Drs. Nicolas Galtier, Sergei Maslov, and Nick Grishin.}, author = {Kondrashov, Alexey S and Povolotskaya, Inna and Ivankov, Dmitry N and Fyodor Kondrashov}, journal = {Biology Direct}, publisher = {BioMed Central}, title = {{Rate of sequence divergence under constant selection}}, doi = {10.1186/1745-6150-5-5}, volume = {5}, year = {2010}, } @article{89, abstract = {We demonstrate the operation of a device that can produce chitosan nanoparticles in a tunable size range from 50-300 nm with small size dispersion. A piezoelectric oscillator operated at megahertz frequencies is used to aerosolize a solution containing dissolved chitosan. The solvent is then evaporated from the aerosolized droplets in a heat pipe, leaving monodisperse nanoparticles to be collected. The nanoparticle size is controlled both by the concentration of the dissolved polymer and by the size of the aerosol droplets that are created. Our device can be used with any polymer or polymer/therapeutic combination that can be prepared in a homogeneous solution and vaporized.}, author = {Wright, Ian and Higginbotham, Andrew P and Baker, Shenda and Donnelly, Tom}, journal = {ACS Applied Materials and Interfaces}, number = {8}, pages = {2360 -- 2364}, publisher = {American Chemical Society}, title = {{Generation of nanoparticles of controlled size using ultrasonic piezoelectric oscillators in solution}}, doi = {10.1021/am100375w}, volume = {2}, year = {2010}, } @article{891, abstract = {Gene duplications and their subsequent divergence play an important part in the evolution of novel gene functions. Several models for the emergence, maintenance and evolution of gene copies have been proposed. However, a clear consensus on how gene duplications are fixed and maintained in genomes is lacking. Here, we present a comprehensive classification of the models that are relevant to all stages of the evolution of gene duplications. Each model predicts a unique combination of evolutionary dynamics and functional properties. Setting out these predictions is an important step towards identifying the main mechanisms that are involved in the evolution of gene duplications.}, author = {Innan, Hideki and Fyodor Kondrashov}, journal = {Nature Reviews Genetics}, number = {2}, pages = {97 -- 108}, publisher = {Nature Publishing Group}, title = {{The evolution of gene duplications: Classifying and distinguishing between models}}, doi = {10.1038/nrg2689}, volume = {11}, year = {2010}, } @article{901, abstract = {Background: Surveying deleterious variation in human populations is crucial for our understanding, diagnosis and potential treatment of human genetic pathologies. A number of recent genome-wide analyses focused on the prevalence of segregating deleterious alleles in the nuclear genome. However, such studies have not been conducted for the mitochondrial genome.Results: We present a systematic survey of polymorphisms in the human mitochondrial genome, including those predicted to be deleterious and those that correspond to known pathogenic mutations. Analyzing 4458 completely sequenced mitochondrial genomes we characterize the genetic diversity of different types of single nucleotide polymorphisms (SNPs) in African (L haplotypes) and non-African (M and N haplotypes) populations. We find that the overall level of polymorphism is higher in the mitochondrial compared to the nuclear genome, although the mitochondrial genome appears to be under stronger selection as indicated by proportionally fewer nonsynonymous than synonymous substitutions. The African mitochondrial genomes show higher heterozygosity, a greater number of polymorphic sites and higher frequencies of polymorphisms for synonymous, benign and damaging polymorphism than non-African genomes. However, African genomes carry significantly fewer SNPs that have been previously characterized as pathogenic compared to non-African genomes.Conclusions: Finding SNPs classified as pathogenic to be the only category of polymorphisms that are more abundant in non-African genomes is best explained by a systematic ascertainment bias that favours the discovery of pathogenic polymorphisms segregating in non-African populations. This further suggests that, contrary to the common disease-common variant hypothesis, pathogenic mutations are largely population-specific and different SNPs may be associated with the same disease in different populations. Therefore, to obtain a comprehensive picture of the deleterious variability in the human population, as well as to improve the diagnostics of individuals carrying African mitochondrial haplotypes, it is necessary to survey different populations independently.Reviewers: This article was reviewed by Dr Mikhail Gelfand, Dr Vasily Ramensky (nominated by Dr Eugene Koonin) and Dr David Rand (nominated by Dr Laurence Hurst).}, author = {Breen, Michael S and Fyodor Kondrashov}, journal = {Biology Direct}, publisher = {BioMed Central}, title = {{Mitochondrial pathogenic mutations are population-specific}}, doi = {10.1186/1745-6150-5-68}, volume = {5}, year = {2010}, } @article{9012, abstract = {In this Letter, we characterize experimentally the diffusiophoretic motion of colloids and λ-DNA toward higher concentration of solutes, using microfluidic technology to build spatially and temporally controlled concentration gradients. We then demonstrate that segregation and spatial patterning of the particles can be achieved from temporal variations of the solute concentration profile. This segregation takes the form of a strong trapping potential, stemming from an osmotically induced rectification mechanism of the solute time-dependent variations. Depending on the spatial and temporal symmetry of the solute signal, localization patterns with various shapes can be achieved. These results highlight the role of solute contrasts in out-of-equilibrium processes occurring in soft matter.}, author = {Palacci, Jérémie A and Abécassis, Benjamin and Cottin-Bizonne, Cécile and Ybert, Christophe and Bocquet, Lydéric}, issn = {10797114}, journal = {Physical Review Letters}, number = {13}, publisher = {American Physical Society}, title = {{Colloidal motility and pattern formation under rectified diffusiophoresis}}, doi = {10.1103/physrevlett.104.138302}, volume = {104}, year = {2010}, } @article{9013, abstract = {In this Letter, we investigate experimentally the nonequilibrium steady state of an active colloidal suspension under gravity field. The active particles are made of chemically powered colloids, showing self propulsion in the presence of an added fuel, here hydrogen peroxide. The active suspension is studied in a dedicated microfluidic device, made of permeable gel microstructures. Both the microdynamics of individual colloids and the global stationary state of the suspension under gravity are measured with optical microscopy. This yields a direct measurement of the effective temperature of the active system as a function of the particle activity, on the basis of the fluctuation-dissipation relationship. Our work is a first step in the experimental exploration of the out-of-equilibrium properties of active colloidal systems.}, author = {Palacci, Jérémie A and Cottin-Bizonne, Cécile and Ybert, Christophe and Bocquet, Lydéric}, issn = {10797114}, journal = {Physical Review Letters}, number = {8}, publisher = {American Physical Society }, title = {{Sedimentation and effective temperature of active colloidal suspensions}}, doi = {10.1103/physrevlett.105.088304}, volume = {105}, year = {2010}, } @article{9145, abstract = {We have found a new way to express the solutions of the RSM (Reynolds Stress Model) equations that allows us to present the turbulent diffusivities for heat, salt and momentum in a way that is considerably simpler and thus easier to implement than in previous work. The RSM provides the dimensionless mixing efficiencies Γα (α stands for heat, salt and momentum). However, to compute the diffusivities, one needs additional information, specifically, the dissipation ε. Since a dynamic equation for the latter that includes the physical processes relevant to the ocean is still not available, one must resort to different sources of information outside the RSM to obtain a complete Mixing Scheme usable in OGCMs. As for the RSM results, we show that the Γα’s are functions of both Ri and Rρ (Richardson number and density ratio representing double diffusion, DD); the Γα are different for heat, salt and momentum; in the case of heat, the traditional value Γh = 0.2 is valid only in the presence of strong shear (when DD is inoperative) while when shear subsides, NATRE data show that Γh can be three times as large, a result that we reproduce. The salt Γs is given in terms of Γh. The momentum Γm has thus far been guessed with different prescriptions while the RSM provides a well defined expression for Γm(Ri, Rρ). Having tested Γh, we then test the momentum Γm by showing that the turbulent Prandtl number Γm/Γh vs. Ri reproduces the available data quite well. As for the dissipation ε, we use different representations, one for the mixed layer (ML), one for the thermocline and one for the ocean’s bottom. For the ML, we adopt a procedure analogous to the one successfully used in PB (planetary boundary layer) studies; for the thermocline, we employ an expression for the variable εN−2 from studies of the internal gravity waves spectra which includes a latitude dependence; for the ocean bottom, we adopt the enhanced bottom diffusivity expression used by previous authors but with a state of the art internal tidal energy formulation and replace the fixed Γα = 0.2 with the RSM result that brings into the problem the Ri, Rρ dependence of the Γα; the unresolved bottom drag, which has thus far been either ignored or modeled with heuristic relations, is modeled using a formalism we previously developed and tested in PBL studies. We carried out several tests without an OGCM. Prandtl and flux Richardson numbers vs. Ri. The RSM model reproduces both types of data satisfactorily. DD and Mixing efficiency Γh(Ri, Rρ). The RSM model reproduces well the NATRE data. Bimodal ε-distribution. NATRE data show that ε(Ri < 1) ≈ 10ε(Ri > 1), which our model reproduces. Heat to salt flux ratio. In the Ri ≫ 1 regime, the RSM predictions reproduce the data satisfactorily. NATRE mass diffusivity. The z-profile of the mass diffusivity reproduces well the measurements at NATRE. The local form of the mixing scheme is algebraic with one cubic equation to solve.}, author = {Canuto, V.M. and Howard, A.M. and Cheng, Y. and Muller, Caroline J and Leboissetier, A. and Jayne, S.R.}, issn = {1463-5003}, journal = {Ocean Modelling}, keywords = {Computer Science (miscellaneous), Geotechnical Engineering and Engineering Geology, Atmospheric Science, Oceanography}, number = {3-4}, pages = {70--91}, publisher = {Elsevier}, title = {{Ocean turbulence, III: New GISS vertical mixing scheme}}, doi = {10.1016/j.ocemod.2010.04.006}, volume = {34}, year = {2010}, } @article{9146, abstract = {The factors governing the rate of change in the amount of atmospheric water vapor are analyzed in simulations of climate change. The global-mean amount of water vapor is estimated to increase at a differential rate of 7.3% K − 1 with respect to global-mean surface air temperature in the multi-model mean. Larger rates of change result if the fractional change is evaluated over a finite change in temperature (e.g., 8.2% K − 1 for a 3 K warming), and rates of change of zonal-mean column water vapor range from 6 to 12% K − 1 depending on latitude. Clausius–Clapeyron scaling is directly evaluated using an invariant distribution of monthly-mean relative humidity, giving a rate of 7.4% K − 1 for global-mean water vapor. There are deviations from Clausius–Clapeyron scaling of zonal-mean column water vapor in the tropics and mid-latitudes, but they largely cancel in the global mean. A purely thermodynamic scaling based on a saturated troposphere gives a higher global rate of 7.9% K − 1. Surface specific humidity increases at a rate of 5.7% K − 1, considerably lower than the rate for global-mean water vapor. Surface specific humidity closely follows Clausius–Clapeyron scaling over ocean. But there are widespread decreases in surface relative humidity over land (by more than 1% K − 1 in many regions), and it is argued that decreases of this magnitude could result from the land/ocean contrast in surface warming.}, author = {O’Gorman, P A and Muller, Caroline J}, issn = {1748-9326}, journal = {Environmental Research Letters}, keywords = {Renewable Energy, Sustainability and the Environment, Public Health, Environmental and Occupational Health, General Environmental Science}, number = {2}, publisher = {IOP Publishing}, title = {{How closely do changes in surface and column water vapor follow Clausius–Clapeyron scaling in climate change simulations?}}, doi = {10.1088/1748-9326/5/2/025207}, volume = {5}, year = {2010}, } @article{920, abstract = {Most eukaryotic cells sense and respond to the mechanical properties of their surroundings. This can strongly influence their collective behavior in embryonic development, tissue function, and wound healing. We use a deformable substrate to measure collective behavior in cell motion due to substrate mediated cell-cell interactions. We quantify spatial and temporal correlations in migration velocity and substrate deformation, and show that cooperative cell-driven patterns of substrate deformation mediate long-distance mechanical coupling between cells and control collective cell migration.}, author = {Angelini, Thomas and Hannezo, Edouard B and Trepat, Xavier and Fredberg, Jeffrey and Weitz, David}, journal = {Physical Review Letters}, number = {16}, publisher = {American Physical Society}, title = {{Cell migration driven by cooperative substrate deformation patterns}}, doi = {10.1103/PhysRevLett.104.168104}, volume = {104}, year = {2010}, } @article{9452, abstract = {Eukaryotic cytosine methylation represses transcription but also occurs in the bodies of active genes, and the extent of methylation biology conservation is unclear. We quantified DNA methylation in 17 eukaryotic genomes and found that gene body methylation is conserved between plants and animals, whereas selective methylation of transposons is not. We show that methylation of plant transposons in the CHG context extends to green algae and that exclusion of histone H2A.Z from methylated DNA is conserved between plants and animals, and we present evidence for RNA-directed DNA methylation of fungal genes. Our data demonstrate that extant DNA methylation systems are mosaics of conserved and derived features, and indicate that gene body methylation is an ancient property of eukaryotic genomes.}, author = {Zemach, Assaf and McDaniel, Ivy E. and Silva, Pedro and Zilberman, Daniel}, issn = {1095-9203}, journal = {Science}, keywords = {Multidisciplinary}, number = {5980}, pages = {916--919}, publisher = {American Association for the Advancement of Science}, title = {{Genome-wide evolutionary analysis of eukaryotic DNA methylation}}, doi = {10.1126/science.1186366}, volume = {328}, year = {2010}, } @article{9485, abstract = {Cytosine methylation silences transposable elements in plants, vertebrates, and fungi but also regulates gene expression. Plant methylation is catalyzed by three families of enzymes, each with a preferred sequence context: CG, CHG (H = A, C, or T), and CHH, with CHH methylation targeted by the RNAi pathway. Arabidopsis thaliana endosperm, a placenta-like tissue that nourishes the embryo, is globally hypomethylated in the CG context while retaining high non-CG methylation. Global methylation dynamics in seeds of cereal crops that provide the bulk of human nutrition remain unknown. Here, we show that rice endosperm DNA is hypomethylated in all sequence contexts. Non-CG methylation is reduced evenly across the genome, whereas CG hypomethylation is localized. CHH methylation of small transposable elements is increased in embryos, suggesting that endosperm demethylation enhances transposon silencing. Genes preferentially expressed in endosperm, including those coding for major storage proteins and starch synthesizing enzymes, are frequently hypomethylated in endosperm, indicating that DNA methylation is a crucial regulator of rice endosperm biogenesis. Our data show that genome-wide reshaping of seed DNA methylation is conserved among angiosperms and has a profound effect on gene expression in cereal crops.}, author = {Zemach, Assaf and Kim, M. Yvonne and Silva, Pedro and Rodrigues, Jessica A. and Dotson, Bradley and Brooks, Matthew D. and Zilberman, Daniel}, issn = {1091-6490}, journal = {Proceedings of the National Academy of Sciences}, number = {43}, pages = {18729--18734}, publisher = {National Academy of Sciences}, title = {{Local DNA hypomethylation activates genes in rice endosperm}}, doi = {10.1073/pnas.1009695107}, volume = {107}, year = {2010}, }