Adjusting Phenotypes by Noise Control

Genetically identical cells can show phenotypic variability. This is often caused by stochastic events that originate from randomness in biochemical processes involving in gene expression and other extrinsic cellular processes. From an engineering perspective, there have been efforts focused on theory and experiments to control noise levels by perturbing and replacing gene network components. However, systematic methods for noise control are lacking mainly due to the intractable mathematical structure of noise propagation through reaction networks. Here, we provide a numerical analysis method by quantifying the parametric sensitivity of noise characteristics at the level of the linear noise approximation. Our analysis is readily applicable to various types of noise control and to different types of system; for example, we can orthogonally control the mean and noise levels and can control system dynamics such as noisy oscillations. As an illustration we applied our method to HIV and yeast gene expression systems and metabolic networks. The oscillatory signal control was applied to p53 oscillations from DNA damage. Furthermore, we showed that the efficiency of orthogonal control can be enhanced by applying extrinsic noise and feedback. Our noise control analysis can be applied to any stochastic model belonging to continuous time Markovian systems such as biological and chemical reaction systems, and even computer and social networks. We anticipate the proposed analysis to be a useful tool for designing and controlling synthetic gene networks

Fast and accurate imputation of summary statistics enhances evidence of functional enrichment

Imputation using external reference panels is a widely used approach for increasing power in GWAS and meta-analysis. Existing HMM-based imputation approaches require individual-level genotypes. Here, we develop a new method for Gaussian imputation from summary association statistics, a type of data that is becoming widely available. In simulations using 1000 Genomes (1000G) data, this method recovers 84% (54%) of the effective sample size for common (>5%) and low-frequency (1-5%) variants (increasing to 87% (60%) when summary LD information is available from target samples) versus 89% (67%) for HMM-based imputation, which cannot be applied to summary statistics. Our approach accounts for the limited sample size of the reference panel, a crucial step to eliminate false-positive associations, and is computationally very fast. As an empirical demonstration, we apply our method to 7 case-control phenotypes from the WTCCC data and a study of height in the British 1958 birth cohort (1958BC). Gaussian imputation from summary statistics recovers 95% (105%) of the effective sample size (as quantified by the ratio of $\chi^2$ association statistics) compared to HMM-based imputation from individual-level genotypes at the 227 (176) published SNPs in the WTCCC (1958BC height) data. In addition, for publicly available summary statistics from large meta-analyses of 4 lipid traits, we publicly release imputed summary statistics at 1000G SNPs, which could not have been obtained using previously published methods, and demonstrate their accuracy by masking subsets of the data. We show that 1000G imputation using our approach increases the magnitude and statistical evidence of enrichment at genic vs. non-genic loci for these traits, as compared to an analysis without 1000G imputation. Thus, imputation of summary statistics will be a valuable tool in future functional enrichment analyses.Comment: 32 pages, 4 figure

Statistical mechanics for metabolic networks during steady-state growth

Which properties of metabolic networks can be derived solely from stoichiometric information about the network's constituent reactions? Predictive results have been obtained by Flux Balance Analysis (FBA), by postulating that cells set metabolic fluxes within the allowed stoichiometry so as to maximize their growth. Here, we generalize this framework to single cell level using maximum entropy models from statistical physics. We define and compute, for the core metabolism of Escherichia coli, a joint distribution over all fluxes that yields the experimentally observed growth rate. This solution, containing FBA as a limiting case, provides a better match to the measured fluxes in the wild type and several mutants. We find that E. coli metabolism is close to, but not at, the optimality assumed by FBA. Moreover, our model makes a wide range of predictions: (i) on flux variability, its regulation, and flux correlations across individual cells; (ii) on the relative importance of stoichiometric constraints vs. growth rate optimization; (iii) on quantitative scaling relations for singe-cell growth rate distributions. We validate these scaling predictions using data from individual bacterial cells grown in a microfluidic device at different sub-inhibitory antibiotic concentrations. Under mild dynamical assumptions, fluctuation-response relations further predict the autocorrelation timescale in growth data and growth rate adaptation times following an environmental perturbation.Comment: 12 pages, 4 figure

Maternal urinary metabolic signatures of fetal growth and associated clinical and environmental factors in the INMA study

Background Maternal metabolism during pregnancy is a major determinant of the intra-uterine environment and fetal outcomes. Herein, we characterize the maternal urinary metabolome throughout pregnancy to identify maternal metabolic signatures of fetal growth in two subcohorts and explain potential sources of variation in metabolic profiles based on lifestyle and clinical data. Methods We used 1H nuclear magnetic resonance (NMR) spectroscopy to characterize maternal urine samples collected in the INMA birth cohort at the first (n = 412 and n = 394, respectively, in Gipuzkoa and Sabadell cohorts) and third trimesters of gestation (n = 417 and 469). Metabolic phenotypes that reflected longitudinal intra- and inter-individual variation were used to predict measures of fetal growth and birth weight. Results A metabolic shift between the first and third trimesters of gestation was characterized by 1H NMR signals arising predominantly from steroid by-products. We identified 10 significant and reproducible metabolic associations in the third trimester with estimated fetal, birth, and placental weight in two independent subcohorts. These included branched-chain amino acids; isoleucine, valine, leucine, alanine and 3 hydroxyisobutyrate (metabolite of valine), which were associated with a significant fetal weight increase at week 34 of up to 2.4 % in Gipuzkoa (P < 0.005) and 1 % in Sabadell (P < 0.05). Other metabolites included pregnancy-related hormone by-products of estrogens and progesterone, and the methyl donor choline. We could explain a total of 48–53 % of the total variance in birth weight of which urine metabolites had an independent predictive power of 12 % adjusting for all other lifestyle/clinical factors. First trimester metabolic phenotypes could not predict reproducibly weight at later stages of development. Physical activity, as well as other modifiable lifestyle/clinical factors, such as coffee consumption, vitamin D intake, and smoking, were identified as potential sources of metabolic variation during pregnancy. Conclusions Significant reproducible maternal urinary metabolic signatures of fetal growth and birth weight are identified for the first time and linked to modifiable lifestyle factors. This novel approach to prenatal screening, combining multiple risk factors, present a great opportunity to personalize pregnancy management and reduce newborn disease risk in later life

Phenotypic flexibility and the evolution of organismal design

Evolutionary biologists often use phenotypic differences between species and between individuals to gain an understanding of organismal design. The focus of much recent attention has been on developmental plasticity – the environmentally induced variability during development within a single genotype. The phenotypic variation expressed by single reproductively mature organisms throughout their life, traditionally the subject of many physiological studies, has remained underexploited in evolutionary biology. Phenotypic flexibility, the reversible within-individual variation, is a function of environmental conditions varying predictably (e.g. with season), or of more stochastic fluctuations in the environment. Here, we provide a common framework to bring the different categories of phenotypic plasticity together, and emphasize perspectives on adaptation that reversible types of plasticity might provide. We argue that better recognition and use of the various levels of phenotypic variation will increase the scope for phenotypic experimentation, comparison and integration.

Accurate estimation of SNP-heritability from biobank-scale data irrespective of genetic architecture.

SNP-heritability is a fundamental quantity in the study of complex traits. Recent studies have shown that existing methods to estimate genome-wide SNP-heritability can yield biases when their assumptions are violated. While various approaches have been proposed to account for frequency- and linkage disequilibrium (LD)-dependent genetic architectures, it remains unclear which estimates reported in the literature are reliable. Here we show that genome-wide SNP-heritability can be accurately estimated from biobank-scale data irrespective of genetic architecture, without specifying a heritability model or partitioning SNPs by allele frequency and/or LD. We show analytically and through extensive simulations starting from real genotypes (UK Biobank, N = 337 K) that, unlike existing methods, our closed-form estimator is robust across a wide range of architectures. We provide estimates of SNP-heritability for 22 complex traits in the UK Biobank and show that, consistent with our results in simulations, existing biobank-scale methods yield estimates up to 30% different from our theoretically-justified approach

Population Structure and Cryptic Relatedness in Genetic Association Studies

We review the problem of confounding in genetic association studies, which arises principally because of population structure and cryptic relatedness. Many treatments of the problem consider only a simple ``island'' model of population structure. We take a broader approach, which views population structure and cryptic relatedness as different aspects of a single confounder: the unobserved pedigree defining the (often distant) relationships among the study subjects. Kinship is therefore a central concept, and we review methods of defining and estimating kinship coefficients, both pedigree-based and marker-based. In this unified framework we review solutions to the problem of population structure, including family-based study designs, genomic control, structured association, regression control, principal components adjustment and linear mixed models. The last solution makes the most explicit use of the kinships among the study subjects, and has an established role in the analysis of animal and plant breeding studies. Recent computational developments mean that analyses of human genetic association data are beginning to benefit from its powerful tests for association, which protect against population structure and cryptic kinship, as well as intermediate levels of confounding by the pedigree.Comment: Published in at http://dx.doi.org/10.1214/09-STS307 the Statistical Science (http://www.imstat.org/sts/) by the Institute of Mathematical Statistics (http://www.imstat.org