A comparison of approximation techniques for variance-based sensitivity analysis of biochemical reaction systems

<p>Abstract</p> <p>Background</p> <p>Sensitivity analysis is an indispensable tool for the analysis of complex systems. In a recent paper, we have introduced a thermodynamically consistent variance-based sensitivity analysis approach for studying the robustness and fragility properties of biochemical reaction systems under uncertainty in the standard chemical potentials of the activated complexes of the reactions and the standard chemical potentials of the molecular species. In that approach, key sensitivity indices were estimated by Monte Carlo sampling, which is computationally very demanding and impractical for large biochemical reaction systems. Computationally efficient algorithms are needed to make variance-based sensitivity analysis applicable to realistic cellular networks, modeled by biochemical reaction systems that consist of a large number of reactions and molecular species.</p> <p>Results</p> <p>We present four techniques, derivative approximation (DA), polynomial approximation (PA), Gauss-Hermite integration (GHI), and orthonormal Hermite approximation (OHA), for <it>analytically </it>approximating the variance-based sensitivity indices associated with a biochemical reaction system. By using a well-known model of the mitogen-activated protein kinase signaling cascade as a case study, we numerically compare the approximation quality of these techniques against traditional Monte Carlo sampling. Our results indicate that, although DA is computationally the most attractive technique, special care should be exercised when using it for sensitivity analysis, since it may only be accurate at low levels of uncertainty. On the other hand, PA, GHI, and OHA are computationally more demanding than DA but can work well at high levels of uncertainty. GHI results in a slightly better accuracy than PA, but it is more difficult to implement. OHA produces the most accurate approximation results and can be implemented in a straightforward manner. It turns out that the computational cost of the four approximation techniques considered in this paper is orders of magnitude smaller than traditional Monte Carlo estimation. Software, coded in MATLAB^®, which implements all sensitivity analysis techniques discussed in this paper, is available free of charge.</p> <p>Conclusions</p> <p>Estimating variance-based sensitivity indices of a large biochemical reaction system is a computationally challenging task that can only be addressed via approximations. Among the methods presented in this paper, a technique based on orthonormal Hermite polynomials seems to be an acceptable candidate for the job, producing very good approximation results for a wide range of uncertainty levels in a fraction of the time required by traditional Monte Carlo sampling.</p

Complete Sequencing and Pan-Genomic Analysis of Lactobacillus delbrueckii subsp. bulgaricus Reveal Its Genetic Basis for Industrial Yogurt Production

Lactobacillus delbrueckii subsp. bulgaricus (Lb. bulgaricus) is an important species of Lactic Acid Bacteria (LAB) used for cheese and yogurt fermentation. The genome of Lb. bulgaricus 2038, an industrial strain mainly used for yogurt production, was completely sequenced and compared against the other two ATCC collection strains of the same subspecies. Specific physiological properties of strain 2038, such as lysine biosynthesis, formate production, aspartate-related carbon-skeleton intermediate metabolism, unique EPS synthesis and efficient DNA restriction/modification systems, are all different from those of the collection strains that might benefit the industrial production of yogurt. Other common features shared by Lb. bulgaricus strains, such as efficient protocooperation with Streptococcus thermophilus and lactate production as well as well-equipped stress tolerance mechanisms may account for it being selected originally for yogurt fermentation industry. Multiple lines of evidence suggested that Lb. bulgaricus 2038 was genetically closer to the common ancestor of the subspecies than the other two sequenced collection strains, probably due to a strict industrial maintenance process for strain 2038 that might have halted its genome decay and sustained a gene network suitable for large scale yogurt production

CCL2-driven inflammation increases mammary gland stromal density and cancer susceptibility in a transgenic mouse model.

Abstract Background Macrophages play diverse roles in mammary gland development and breast cancer. CC-chemokine ligand 2 (CCL2) is an inflammatory cytokine that recruits macrophages to sites of injury. Although CCL2 has been detected in human and mouse mammary epithelium, its role in regulating mammary gland development and cancer risk has not been explored. Methods Transgenic mice were generated wherein CCL2 is driven by the mammary epithelial cell-specific mouse mammary tumour virus 206 (MMTV) promoter. Estrous cycles were tracked in adult transgenic and non-transgenic FVB mice, and mammary glands collected at the four different stages of the cycle. Dissected mammary glands were assessed for cyclical morphological changes, proliferation and apoptosis of epithelium, macrophage abundance and collagen deposition, and mRNA encoding matrix remodelling enzymes. Another cohort of control and transgenic mice received carcinogen 7,12-Dimethylbenz(a)anthracene (DMBA) and tumour development was monitored weekly. CCL2 protein was also quantified in paired samples of human breast tissue with high and low mammographic density. Results Overexpression of CCL2 in the mammary epithelium resulted in an increased number of macrophages, increased density of stroma and collagen and elevated mRNA encoding matrix remodelling enzymes lysyl oxidase (LOX) and tissue inhibitor of matrix metalloproteinases (TIMP)3 compared to non-transgenic controls. Transgenic mice also exhibited increased susceptibility to development of DMBA-induced mammary tumours. In a paired sample cohort of human breast tissue, abundance of epithelial-cell-associated CCL2 was higher in breast tissue of high mammographic density compared to tissue of low mammographic density. Conclusions Constitutive expression of CCL2 by the mouse mammary epithelium induces a state of low level chronic inflammation that increases stromal density and elevates cancer risk. We propose that CCL2-driven inflammation contributes to the increased risk of breast cancer observed in women with high mammographic density

Fine Mapping of the NRG1 Hirschsprung's Disease Locus

The primary pathology of Hirschsprung's disease (HSCR, colon aganglionosis) is the absence of ganglia in variable lengths of the hindgut, resulting in functional obstruction. HSCR is attributed to a failure of migration of the enteric ganglion precursors along the developing gut. RET is a key regulator of the development of the enteric nervous system (ENS) and the major HSCR-causing gene. Yet the reduced penetrance of RET DNA HSCR-associated variants together with the phenotypic variability suggest the involvement of additional genes in the disease. Through a genome-wide association study, we uncovered a ∼350 kb HSCR-associated region encompassing part of the neuregulin-1 gene (NRG1). To identify the causal NRG1 variants contributing to HSCR, we genotyped 243 SNPs variants on 343 ethnic Chinese HSCR patients and 359 controls. Genotype analysis coupled with imputation narrowed down the HSCR-associated region to 21 kb, with four of the most associated SNPs (rs10088313, rs10094655, rs4624987, and rs3884552) mapping to the NRG1 promoter. We investigated whether there was correlation between the genotype at the rs10088313 locus and the amount of NRG1 expressed in human gut tissues (40 patients and 21 controls) and found differences in expression as a function of genotype. We also found significant differences in NRG1 expression levels between diseased and control individuals bearing the same rs10088313 risk genotype. This indicates that the effects of NRG1 common variants are likely to depend on other alleles or epigenetic factors present in the patients and would account for the variability in the genetic predisposition to HSCR

Quantitative High-Throughput Screening Identifies 8-Hydroxyquinolines as Cell-Active Histone Demethylase Inhibitors

Small molecule modulators of epigenetic processes are currently sought as basic probes for biochemical mechanisms, and as starting points for development of therapeutic agents. N(epsilon)-Methylation of lysine residues on histone tails is one of a number of post-translational modifications that together enable transcriptional regulation. Histone lysine demethylases antagonize the action of histone methyltransferases in a site- and methylation state-specific manner. N(epsilon)-Methyllysine demethylases that use 2-oxoglutarate as co-factor are associated with diverse human diseases, including cancer, inflammation and X-linked mental retardation; they are proposed as targets for the therapeutic modulation of transcription. There are few reports on the identification of templates that are amenable to development as potent inhibitors in vivo and large diverse collections have yet to be exploited for the discovery of demethylase inhibitors

Prevalence and attributable health burden of chronic respiratory diseases, 1990–2017: A systematic analysis for the global burden of disease study 2017

© 2020 The Author(s). Published by Elsevier Ltd. This is an Open Access article under the CC BY 4.0 license Background: Previous attempts to characterise the burden of chronic respiratory diseases have focused only on specific disease conditions, such as chronic obstructive pulmonary disease (COPD) or asthma. In this study, we aimed to characterise the burden of chronic respiratory diseases globally, providing a comprehensive and up-to-date analysis on geographical and time trends from 1990 to 2017. Methods: Using data from the Global Burden of Diseases, Injuries, and Risk Factors Study (GBD) 2017, we estimated the prevalence, morbidity, and mortality attributable to chronic respiratory diseases through an analysis of deaths, disability-adjusted life-years (DALYs), and years of life lost (YLL) by GBD super-region, from 1990 to 2017, stratified by age and sex. Specific diseases analysed included asthma, COPD, interstitial lung disease and pulmonary sarcoidosis, pneumoconiosis, and other chronic respiratory diseases. We also assessed the contribution of risk factors (smoking, second-hand smoke, ambient particulate matter and ozone pollution, household air pollution from solid fuels, and occupational risks) to chronic respiratory disease-attributable DALYs. Findings: In 2017, 544·9 million people (95% uncertainty interval [UI] 506·9–584·8) worldwide had a chronic respiratory disease, representing an increase of 39·8% compared with 1990. Chronic respiratory disease prevalence showed wide variability across GBD super-regions, with the highest prevalence among both males and females in high-income regions, and the lowest prevalence in sub-Saharan Africa and south Asia. The age-sex-specific prevalence of each chronic respiratory disease in 2017 was also highly variable geographically. Chronic respiratory diseases were the third leading cause of death in 2017 (7·0% [95% UI 6·8–7·2] of all deaths), behind cardiovascular diseases and neoplasms. Deaths due to chronic respiratory diseases numbered 3 914 196 (95% UI 3 790 578–4 044 819) in 2017, an increase of 18·0% since 1990, while total DALYs increased by 13·3%. However, when accounting for ageing and population growth, declines were observed in age-standardised prevalence (14·3% decrease), age-standardised death rates (42·6%), and age-standardised DALY rates (38·2%). In males and females, most chronic respiratory disease-attributable deaths and DALYs were due to COPD. In regional analyses, mortality rates from chronic respiratory diseases were greatest in south Asia and lowest in sub-Saharan Africa, also across both sexes. Notably, although absolute prevalence was lower in south Asia than in most other super-regions, YLLs due to chronic respiratory diseases across the subcontinent were the highest in the world. Death rates due to interstitial lung disease and pulmonary sarcoidosis were greater than those due to pneumoconiosis in all super-regions. Smoking was the leading risk factor for chronic respiratory disease-related disability across all regions for men. Among women, household air pollution from solid fuels was the predominant risk factor for chronic respiratory diseases in south Asia and sub-Saharan Africa, while ambient particulate matter represented the leading risk factor in southeast Asia, east Asia, and Oceania, and in the Middle East and north Africa super-region. Interpretation: Our study shows that chronic respiratory diseases remain a leading cause of death and disability worldwide, with growth in absolute numbers but sharp declines in several age-standardised estimators since 1990. Premature mortality from chronic respiratory diseases seems to be highest in regions with less-resourced health systems on a per-capita basis. Funding: Bill & Melinda Gates Foundation

Northern Eurasia Future Initiative (NEFI): facing the challenges and pathways of global change in the 21st century

During the past several decades, the Earth system has changed significantly, especially across Northern Eurasia. Changes in the socio-economic conditions of the larger countries in the region have also resulted in a variety of regional environmental changes that can have global consequences. The Northern Eurasia Future Initiative (NEFI) has been designed as an essential continuation of the Northern Eurasia Earth Science Partnership Initiative (NEESPI), which was launched in 2004. NEESPI sought to elucidate all aspects of ongoing environmental change, to inform societies and, thus, to better prepare societies for future developments. A key principle of NEFI is that these developments must now be secured through science-based strategies co-designed with regional decision makers to lead their societies to prosperity in the face of environmental and institutional challenges. NEESPI scientific research, data, and models have created a solid knowledge base to support the NEFI program. This paper presents the NEFI research vision consensus based on that knowledge. It provides the reader with samples of recent accomplishments in regional studies and formulates new NEFI science questions. To address these questions, nine research foci are identified and their selections are briefly justified. These foci include: warming of the Arctic; changing frequency, pattern, and intensity of extreme and inclement environmental conditions; retreat of the cryosphere; changes in terrestrial water cycles; changes in the biosphere; pressures on land-use; changes in infrastructure; societal actions in response to environmental change; and quantification of Northern Eurasia's role in the global Earth system. Powerful feedbacks between the Earth and human systems in Northern Eurasia (e.g., mega-fires, droughts, depletion of the cryosphere essential for water supply, retreat of sea ice) result from past and current human activities (e.g., large scale water withdrawals, land use and governance change) and potentially restrict or provide new opportunities for future human activities. Therefore, we propose that Integrated Assessment Models are needed as the final stage of global change assessment. The overarching goal of this NEFI modeling effort will enable evaluation of economic decisions in response to changing environmental conditions and justification of mitigation and adaptation efforts

Global, regional, and national incidence, prevalence, and years lived with disability for 354 diseases and injuries for 195 countries and territories, 1990-2017: a systematic analysis for the Global Burden of Disease Study 2017.

The Global Burden of Diseases, Injuries and Risk Factors 2017 includes a comprehensive assessment of incidence, prevalence, and years lived with disability (YLDs) for 354 causes in 195 countries and territories from 1990 to 2017. Previous GBD studies have shown how the decline of mortality rates from 1990 to 2016 has led to an increase in life expectancy, an ageing global population, and an expansion of the non-fatal burden of disease and injury. These studies have also shown how a substantial portion of the world's population experiences non-fatal health loss with considerable heterogeneity among different causes, locations, ages, and sexes. Ongoing objectives of the GBD study include increasing the level of estimation detail, improving analytical strategies, and increasing the amount of high-quality data. METHODS: We estimated incidence and prevalence for 354 diseases and injuries and 3484 sequelae. We used an updated and extensive body of literature studies, survey data, surveillance data, inpatient admission records, outpatient visit records, and health insurance claims, and additionally used results from cause of death models to inform estimates using a total of 68 781 data sources. Newly available clinical data from India, Iran, Japan, Jordan, Nepal, China, Brazil, Norway, and Italy were incorporated, as well as updated claims data from the USA and new claims data from Taiwan (province of China) and Singapore. We used DisMod-MR 2.1, a Bayesian meta-regression tool, as the main method of estimation, ensuring consistency between rates of incidence, prevalence, remission, and cause of death for each condition. YLDs were estimated as the product of a prevalence estimate and a disability weight for health states of each mutually exclusive sequela, adjusted for comorbidity. We updated the Socio-demographic Index (SDI), a summary development indicator of income per capita, years of schooling, and total fertility rate. Additionally, we calculated differences between male and female YLDs to identify divergent trends across sexes. GBD 2017 complies with the Guidelines for Accurate and Transparent Health Estimates Reporting

Population and fertility by age and sex for 195 countries and territories, 1950–2017: a systematic analysis for the Global Burden of Disease Study 2017

Background: Population estimates underpin demographic and epidemiological research and are used to track progress on numerous international indicators of health and development. To date, internationally available estimates of population and fertility, although useful, have not been produced with transparent and replicable methods and do not use standardised estimates of mortality. We present single-calendar year and single-year of age estimates of fertility and population by sex with standardised and replicable methods. Methods: We estimated population in 195 locations by single year of age and single calendar year from 1950 to 2017 with standardised and replicable methods. We based the estimates on the demographic balancing equation, with inputs of fertility, mortality, population, and migration data. Fertility data came from 7817 location-years of vital registration data, 429 surveys reporting complete birth histories, and 977 surveys and censuses reporting summary birth histories. We estimated age-specific fertility rates (ASFRs; the annual number of livebirths to women of a specified age group per 1000 women in that age group) by use of spatiotemporal Gaussian process regression and used the ASFRs to estimate total fertility rates (TFRs; the average number of children a woman would bear if she survived through the end of the reproductive age span [age 10–54 years] and experienced at each age a particular set of ASFRs observed in the year of interest). Because of sparse data, fertility at ages 10–14 years and 50–54 years was estimated from data on fertility in women aged 15–19 years and 45–49 years, through use of linear regression. Age-specific mortality data came from the Global Burden of Diseases, Injuries, and Risk Factors Study (GBD) 2017 estimates. Data on population came from 1257 censuses and 761 population registry location-years and were adjusted for underenumeration and age misreporting with standard demographic methods. Migration was estimated with the GBD Bayesian demographic balancing model, after incorporating information about refugee migration into the model prior. Final population estimates used the cohort-component method of population projection, with inputs of fertility, mortality, and migration data. Population uncertainty was estimated by use of out-of-sample predictive validity testing. With these data, we estimated the trends in population by age and sex and in fertility by age between 1950 and 2017 in 195 countries and territories. Findings: From 1950 to 2017, TFRs decreased by 49\ub74% (95% uncertainty interval [UI] 46\ub74–52\ub70). The TFR decreased from 4\ub77 livebirths (4\ub75–4\ub79) to 2\ub74 livebirths (2\ub72–2\ub75), and the ASFR of mothers aged 10–19 years decreased from 37 livebirths (34–40) to 22 livebirths (19–24) per 1000 women. Despite reductions in the TFR, the global population has been increasing by an average of 83\ub78 million people per year since 1985. The global population increased by 197\ub72% (193\ub73–200\ub78) since 1950, from 2\ub76 billion (2\ub75–2\ub76) to 7\ub76 billion (7\ub74–7\ub79) people in 2017; much of this increase was in the proportion of the global population in south Asia and sub-Saharan Africa. The global annual rate of population growth increased between 1950 and 1964, when it peaked at 2\ub70%; this rate then remained nearly constant until 1970 and then decreased to 1\ub71% in 2017. Population growth rates in the southeast Asia, east Asia, and Oceania GBD super-region decreased from 2\ub75% in 1963 to 0\ub77% in 2017, whereas in sub-Saharan Africa, population growth rates were almost at the highest reported levels ever in 2017, when they were at 2\ub77%. The global average age increased from 26\ub76 years in 1950 to 32\ub71 years in 2017, and the proportion of the population that is of working age (age 15–64 years) increased from 59\ub79% to 65\ub73%. At the national level, the TFR decreased in all countries and territories between 1950 and 2017; in 2017, TFRs ranged from a low of 1\ub70 livebirths (95% UI 0\ub79–1\ub72) in Cyprus to a high of 7\ub71 livebirths (6\ub78–7\ub74) in Niger. The TFR under age 25 years (TFU25; number of livebirths expected by age 25 years for a hypothetical woman who survived the age group and was exposed to current ASFRs) in 2017 ranged from 0\ub708 livebirths (0\ub707–0\ub709) in South Korea to 2\ub74 livebirths (2\ub72–2\ub76) in Niger, and the TFR over age 30 years (TFO30; number of livebirths expected for a hypothetical woman ageing from 30 to 54 years who survived the age group and was exposed to current ASFRs) ranged from a low of 0\ub73 livebirths (0\ub73–0\ub74) in Puerto Rico to a high of 3\ub71 livebirths (3\ub70–3\ub72) in Niger. TFO30 was higher than TFU25 in 145 countries and territories in 2017. 33 countries had a negative population growth rate from 2010 to 2017, most of which were located in central, eastern, and western Europe, whereas population growth rates of more than 2\ub70% were seen in 33 of 46 countries in sub-Saharan Africa. In 2017, less than 65% of the national population was of working age in 12 of 34 high-income countries, and less than 50% of the national population was of working age in Mali, Chad, and Niger. Interpretation: Population trends create demographic dividends and headwinds (ie, economic benefits and detriments) that affect national economies and determine national planning needs. Although TFRs are decreasing, the global population continues to grow as mortality declines, with diverse patterns at the national level and across age groups. To our knowledge, this is the first study to provide transparent and replicable estimates of population and fertility, which can be used to inform decision making and to monitor progress. Funding: Bill &amp; Melinda Gates Foundation

Global, regional, and national disability-adjusted life-years (DALYs) for 359 diseases and injuries and healthy life expectancy (HALE) for 195 countries and territories, 1990-2017: a systematic analysis for the Global Burden of Disease Study 2017.

How long one lives, how many years of life are spent in good and poor health, and how the population's state of health and leading causes of disability change over time all have implications for policy, planning, and provision of services. We comparatively assessed the patterns and trends of healthy life expectancy (HALE), which quantifies the number of years of life expected to be lived in good health, and the complementary measure of disability-adjusted life-years (DALYs), a composite measure of disease burden capturing both premature mortality and prevalence and severity of ill health, for 359 diseases and injuries for 195 countries and territories over the past 28 years. Methods We used data for age-specific mortality rates, years of life lost (YLLs) due to premature mortality, and years lived with disability (YLDs) from the Global Burden of Diseases, Injuries, and Risk Factors Study (GBD) 2017 to calculate HALE and DALYs from 1990 to 2017. We calculated HALE using age-specific mortality rates and YLDs per capita for each location, age, sex, and year. We calculated DALYs for 359 causes as the sum of YLLs and YLDs. We assessed how observed HALE and DALYs differed by country and sex from expected trends based on Socio-demographic Index (SDI). We also analysed HALE by decomposing years of life gained into years spent in good health and in poor health, between 1990 and 2017, and extra years lived by females compared with males. Findings Globally, from 1990 to 2017, life expectancy at birth increased by 7·4 years (95% uncertainty interval 7·1-7·8), from 65·6 years (65·3-65·8) in 1990 to 73·0 years (72·7-73·3) in 2017. The increase in years of life varied from 5·1 years (5·0-5·3) in high SDI countries to 12·0 years (11·3-12·8) in low SDI countries. Of the additional years of life expected at birth, 26·3% (20·1-33·1) were expected to be spent in poor health in high SDI countries compared with 11·7% (8·8-15·1) in low-middle SDI countries. HALE at birth increased by 6·3 years (5·9-6·7), from 57·0 years (54·6-59·1) in 1990 to 63·3 years (60·5-65·7) in 2017. The increase varied from 3·8 years (3·4-4·1) in high SDI countries to 10·5 years (9·8-11·2) in low SDI countries. Even larger variations in HALE than these were observed between countries, ranging from 1·0 year (0·4-1·7) in Saint Vincent and the Grenadines (62·4 years [59·9-64·7] in 1990 to 63·5 years [60·9-65·8] in 2017) to 23·7 years (21·9-25·6) in Eritrea (30·7 years [28·9-32·2] in 1990 to 54·4 years [51·5-57·1] in 2017). In most countries, the increase in HALE was smaller than the increase in overall life expectancy, indicating more years lived in poor health. In 180 of 195 countries and territories, females were expected to live longer than males in 2017, with extra years lived varying from 1·4 years (0·6-2·3) in Algeria to 11·9 years (10·9-12·9) in Ukraine. Of the extra years gained, the proportion spent in poor health varied largely across countries, with less than 20% of additional years spent in poor health in Bosnia and Herzegovina, Burundi, and Slovakia, whereas in Bahrain all the extra years were spent in poor health. In 2017, the highest estimate of HALE at birth was in Singapore for both females (75·8 years [72·4-78·7]) and males (72·6 years [69·8-75·0]) and the lowest estimates were in Central African Republic (47·0 years [43·7-50·2] for females and 42·8 years [40·1-45·6] for males). Globally, in 2017, the five leading causes of DALYs were neonatal disorders, ischaemic heart disease, stroke, lower respiratory infections, and chronic obstructive pulmonary disease. Between 1990 and 2017, age-standardised DALY rates decreased by 41·3% (38·8-43·5) for communicable diseases and by 49·8% (47·9-51·6) for neonatal disorders. For non-communicable diseases, global DALYs increased by 40·1% (36·8-43·0), although age-standardised DALY rates decreased by 18·1% (16·0-20·2)