Predicting Phenotypic Diversity and the Underlying Quantitative Molecular Transitions

During development, signaling networks control the formation of multicellular patterns. To what extent quantitative fluctuations in these complex networks may affect multicellular phenotype remains unclear. Here, we describe a computational approach to predict and analyze the phenotypic diversity that is accessible to a developmental signaling network. Applying this framework to vulval development in C. elegans, we demonstrate that quantitative changes in the regulatory network can render ~500 multicellular phenotypes. This phenotypic capacity is an order-of-magnitude below the theoretical upper limit for this system but yet is large enough to demonstrate that the system is not restricted to a select few outcomes. Using metrics to gauge the robustness of these phenotypes to parameter perturbations, we identify a select subset of novel phenotypes that are the most promising for experimental validation. In addition, our model calculations provide a layout of these phenotypes in network parameter space. Analyzing this landscape of multicellular phenotypes yielded two significant insights. First, we show that experimentally well-established mutant phenotypes may be rendered using non-canonical network perturbations. Second, we show that the predicted multicellular patterns include not only those observed in C. elegans, but also those occurring exclusively in other species of the Caenorhabditis genus. This result demonstrates that quantitative diversification of a common regulatory network is indeed demonstrably sufficient to generate the phenotypic differences observed across three major species within the Caenorhabditis genus. Using our computational framework, we systematically identify the quantitative changes that may have occurred in the regulatory network during the evolution of these species. Our model predictions show that significant phenotypic diversity may be sampled through quantitative variations in the regulatory network without overhauling the core network architecture. Furthermore, by comparing the predicted landscape of phenotypes to multicellular patterns that have been experimentally observed across multiple species, we systematically trace the quantitative regulatory changes that may have occurred during the evolution of the Caenorhabditis genus

Unlimited multistability in multisite phosphorylation systems

Reversible phosphorylation on serine, threonine and tyrosine is the most widely studied posttranslational modification of proteins (1, 2). The number of phosphorylated sites on a protein (n) shows a significant increase from prokaryotes, with n less than or equal to 7 sites, to eukaryotes, with examples having n greater than or equal to 150 sites (3). Multisite phosphorylation has many roles (4, 5) and site conservation indicates that increasing numbers of sites cannot be due merely to promiscuous phosphorylation. A substrate with n sites has an exponential number (2^n) of phospho-forms and individual phospho-forms may have distinct biological effects (6, 7). The distribution of these phospho-forms and how this distribution is regulated have remained unknown. Here we show that, when kinase and phosphatase act in opposition on a multisite substrate, the system can exhibit distinct stable phospho-form distributions at steady state and that the maximum number of such distributions increases with n. Whereas some stable distributions are focused on a single phospho-form, others are more diffuse, giving the phospho-proteome the potential to behave as a fluid regulatory network able to encode information and flexibly respond to varying demands. Such plasticity may underlie complex information processing in eukaryotic cells (8) and suggests a functional advantage in having many sites. Our results follow from the unusual geometry of the steady-state phospho-form concentrations, which we show to constitute a rational algebraic curve, irrespective of n. We thereby reduce the complexity of calculating steady states from simulating 3 times 2^n differential equations to solving two algebraic equations, while treating parameters symbolically. We anticipate that these methods can be extended to systems with multiple substrates and multiple enzymes catalysing different modifications, as found in posttranslational modification 'codes' (9) such as the histone code (10, 11). Whereas simulations struggle with exponentially increasing molecular complexity, mathematical methods of the kind developed here can provide a new language in which to articulate the principles of cellular information processing (12)

Decomposition techniques with mixed integer programming and heuristics for home healthcare planning

We tackle home healthcare planning scenarios in the UK using decomposition methods that incorporate mixed integer programming solvers and heuristics. Home healthcare planning is a difficult problem that integrates aspects from scheduling and routing. Solving real-world size instances of these problems still presents a significant challenge to modern exact optimization solvers. Nevertheless, we propose decomposition techniques to harness the power of such solvers while still offering a practical approach to produce high-quality solutions to real-world problem instances. We first decompose the problem into several smaller sub-problems. Next, mixed integer programming and/or heuristics are used to tackle the sub-problems. Finally, the sub-problem solutions are combined into a single valid solution for the whole problem. The different decomposition methods differ in the way in which subproblems are generated and the way in which conflicting assignments are tackled (i.e. avoided or repaired). We present the results obtained by the proposed decomposition methods and compare them to solutions obtained with other methods. In addition, we conduct a study that reveals how the different steps in the proposed method contribute to those results. The main contribution of this paper is a better understanding of effective ways to combine mixed integer programming within effective decomposition methods to solve real-world instances of home healthcare planning problems in practical computation time

Relation between sleep quality and quantity, quality of life, and risk of developing diabetes in healthy workers in Japan: the High-risk and Population Strategy for Occupational Health Promotion (HIPOP-OHP) Study

<p>Abstract</p> <p>Background</p> <p>The effect of sleep on the risk of developing diabetes has not been explored in an Asian population. The objective of this study is to investigate the effect of self-reported sleep duration and sleep quality on the risk of developing diabetes in a prospective cohort in Japan.</p> <p>Methods</p> <p>Data were analyzed from the cohort of participants in a High-risk and Population Strategy for Occupational Health Promotion Study (HIPOP-OHP), conducted in Japan from the year 1999 until 2004. A Cox proportional hazard model was used to evaluate the association between sleep duration or sleep quality and the risk of diabetes.</p> <p>Results</p> <p>Of 6509 participants (26.1% of women, 19–69 years of age), a total of 230 type 2 diabetes cases were reported over a median 4.2 years of follow-up. For participants who often experienced difficulty in initiating sleep, the multivariate-adjusted hazard ratios for diabetes were 1.42 (95%CI, 1.05–1.91) in participants with a medium frequency of difficulty initiating sleep, and 1.61 (95%CI, 1.00–2.58) for those with a high frequency, with a statistically significant linear trend. Significant association was not observed in the association between difficulty of maintaining sleep or duration of sleep, and risk of diabetes.</p> <p>Conclusion</p> <p>Medium and high frequencies of difficulty initiating sleep, but not difficulty in maintaining sleep or in sleep duration, are associated with higher risks of diabetes in relatively healthy Asian workers, even after adjusting for a large number of possible further factors.</p

Utilización de la yoga en la mejora del peso corporal y flujo pico espiratorio de personas con sobrepeso y obesidad

The present investigation is of the quasi-experimental type of temporal series design interrupted with a group with repeated measures of pretest and posttest. The objective of the same was the evaluation of a Yoga training program on the effects on the Peak Expiratory Flow (PEF) and on overweight-obesity of female patients who attended the Mutual Aid Groups of the Secretary of Health in Manzanillo, Mexico. The intervention group consisted of 15 women with an average age of 46 years (16) without chronic non-communicable diseases, the program was designed for 6 with two weekly sessions of 40 minutes each (80 min per week). The instruments used to verify the progress in the PEF and the BMI are: clinical history, portable flowmeters and a minimum protocol of weight and height Quetelet´s formula to determine the BMI, as well as an interview at the end of the intervention to know aspects about participation in the intervention. Among the main findings, there was a significant progress in Peak Expiratory Flow (PEF) after the intervention, the participants were able to make a deeper inhalation and a greater exhalation. &nbsp;The results show positive effects on the PEF after a 6-week program with specific yoga work to improve breathing, but not for BMI and changes in body weight. &nbsp;La presente investigación es de tipo cuasiexperimental de diseño serie temporal interrumpida con un grupo con medidas repetidas de pretest y postest. El objetivodel mismo, fue la evaluación de un programa de entrenamiento de Yoga sobre los efectosen elFlujo Pico Espiratorio (PEF) yen sobrepeso-obesidad de pacientes femeninos que asistieron a los Grupos de Ayuda Mutua de la Secretaria de Salud en Manzanillo, México. El grupo de intervención estuvo integrado por 15 mujeres con un promedio de edad de 46 años (±16) sin enfermedades crónicas no transmisibles, el programa se diseñó para 6 semanas con dos sesiones semanales de 40 minutos cada una (80 min por semana). Los instrumentos utilizados para verificar el progreso en el PEF e IMC fueron,historia clínica, flujometros portátiles y un protocolo mínimo de toma de peso y talla y la fórmula de Quetelet para determinar el IMC, así como, una entrevista al final de la intervención para conocer aspectos sobre la participación en la intervención. Entre los principales hallazgos se evidencia un progreso significativo en Flujo Pico Espiratorio (PEF) posterior a la intervención, las participantes fueron capaces de hacer una inspiración más profunda y una espiración mayor. Los resultados muestran efectos positivos en el PEF posterior a un programa de 6 semanas con trabajo de yoga específico para la mejora de respiración, pero no para disminución de IMC y cambios en el peso corporal

Genome3D: integrating a collaborative data pipeline to expand the depth and breadth of consensus protein structure annotation.

Genome3D (https://www.genome3d.eu) is a freely available resource that provides consensus structural annotations for representative protein sequences taken from a selection of model organisms. Since the last NAR update in 2015, the method of data submission has been overhauled, with annotations now being 'pushed' to the database via an API. As a result, contributing groups are now able to manage their own structural annotations, making the resource more flexible and maintainable. The new submission protocol brings a number of additional benefits including: providing instant validation of data and avoiding the requirement to synchronise releases between resources. It also makes it possible to implement the submission of these structural annotations as an automated part of existing internal workflows. In turn, these improvements facilitate Genome3D being opened up to new prediction algorithms and groups. For the latest release of Genome3D (v2.1), the underlying dataset of sequences used as prediction targets has been updated using the latest reference proteomes available in UniProtKB. A number of new reference proteomes have also been added of particular interest to the wider scientific community: cow, pig, wheat and mycobacterium tuberculosis. These additions, along with improvements to the underlying predictions from contributing resources, has ensured that the number of annotations in Genome3D has nearly doubled since the last NAR update article. The new API has also been used to facilitate the dissemination of Genome3D data into InterPro, thereby widening the visibility of both the annotation data and annotation algorithms

A new era in palaeomicrobiology: prospects for ancient dental calculus as a long-term record of the human oral microbiome.

The field of palaeomicrobiology is dramatically expanding thanks to recent advances in high-throughput biomolecular sequencing, which allows unprecedented access to the evolutionary history and ecology of human-associated and environmental microbes. Recently, human dental calculus has been shown to be an abundant, nearly ubiquitous, and long-term reservoir of the ancient oral microbiome, preserving not only microbial and host biomolecules but also dietary and environmental debris. Modern investigations of native human microbiota have demonstrated that the human microbiome plays a central role in health and chronic disease, raising questions about changes in microbial ecology, diversity and function through time. This paper explores the current state of ancient oral microbiome research and discusses successful applications, methodological challenges and future possibilities in elucidating the intimate evolutionary relationship between humans and their microbes

Multisite Phosphorylation Provides an Effective and Flexible Mechanism for Switch-Like Protein Degradation

Phosphorylation-triggered degradation is a common strategy for elimination of regulatory proteins in many important cell signaling processes. Interesting examples include cyclin-dependent kinase inhibitors such as p27 in human and Sic1 in yeast, which play crucial roles during the G1/S transition in the cell cycle. In this work, we have modeled and analyzed the dynamics of multisite-phosphorylation-triggered protein degradation systematically. Inspired by experimental observations on the Sic1 protein and a previous intriguing theoretical conjecture, we develop a model to examine in detail the degradation dynamics of a protein featuring multiple phosphorylation sites and a threshold site number for elimination in response to a kinase signal. Our model explains the role of multiple phosphorylation sites, compared to a single site, in the regulation of protein degradation. A single-site protein cannot convert a graded input of kinase increase to much sharper output, whereas multisite phosphorylation is capable of generating a highly switch-like temporal profile of the substrate protein with two characteristics: a temporal threshold and rapid decrease beyond the threshold. We introduce a measure termed temporal response coefficient to quantify the extent to which a response in the time domain is switch-like and further investigate how this property is determined by various factors including the kinase input, the total number of sites, the threshold site number for elimination, the order of phosphorylation, the kinetic parameters, and site preference. Some interesting and experimentally verifiable predictions include that the non-degradable fraction of the substrate protein exhibits a more switch-like temporal profile; a sequential system is more switch-like, while a random system has the advantage of increased robustness; all the parameters, including the total number of sites, the threshold site number for elimination and the kinetic parameters synergistically determine the exact extent to which the degradation profile is switch-like. Our results suggest design principles for protein degradation switches which might be a widespread mechanism for precise regulation of cellular processes such as cell cycle progression

The Magnitude and Kinetics of the Mucosal HIV-Specific CD8+ T Lymphocyte Response and Virus RNA Load in Breast Milk

BACKGROUND: The risk of postnatal HIV transmission is associated with the magnitude of the milk virus load. While HIV-specific cellular immune responses control systemic virus load and are detectable in milk, the contribution of these responses to the control of virus load in milk is unknown. METHODS: We assessed the magnitude of the immunodominant GagRY11 and subdominant EnvKY9-specific CD8+ T lymphocyte response in blood and milk of 10 A*3002+, HIV-infected Malawian women throughout the period of lactation and correlated this response to milk virus RNA load and markers of breast inflammation. RESULTS: The magnitude and kinetics of the HIV-specific CD8+ T lymphocyte responses were discordant in blood and milk of the right and left breast, indicating independent regulation of these responses in each breast. However, there was no correlation between the magnitude of the HIV-specific CD8+ T lymphocyte response and the milk virus RNA load. Further, there was no correlation between the magnitude of this response and markers of breast inflammation. CONCLUSIONS: The magnitude of the HIV-specific CD8+ T lymphocyte response in milk does not appear to be solely determined by the milk virus RNA load and is likely only one of the factors contributing to maintenance of low virus load in milk

Effects of APOE, APOB and LDLR variants on serum lipids and lack of association with xanthelasma in individuals from Southeastern Brazil

Xanthelasma might be a clinical manifestation of dyslipidemia, a recognized risk factor for coronary artery disease. We investigated the association of apolipoprotein E (APOE HhaI), apolipoprotein B (APOB XbaI and Ins/Del) and LDL receptor (LDLR AvaII and HincII) gene polymorphisms with lipid profiles in 100 Brazilians with xanthelasma and 100 controls. Allele frequencies were similar in both groups. APOE, APOB and LDLR genotypes were not correlated with differences in the serum lipid profile. In individuals with xanthelasma, the APOB D allele was associated with less chance of having increased LDL-cholesterol (O.R. = 0.16, CI95% = 0.03-0.94, p = 0.042). In the control group, the APOB X+ allele was associated with less chance of having both increased total cholesterol (O.R. = 0.16, CI95% = 0.03-0.78, p = 0.023) and increased LDL-cholesterol (O.R. = 0.10, CI95% = 0.02-0.60, p = 0.012). Moreover, there was a significantly higher frequency of control individuals (68%) with elevated serum triglyceride levels, compared to patients (48%, p = 0.008). On the other hand, triglyceride levels in controls also seemed to be influenced by all other gene polymorphisms studied, an effect that might be enhanced by environmental factors