NOGEA: A Network-oriented Gene Entropy Approach for Dissecting Disease Comorbidity and Drug Repositioning

Rapid development of high-throughput technologies has permitted the identification of an increasing number of disease-associated genes (DAGs), which are important for understanding disease initiation and developing precision therapeutics. However, DAGs often contain large amounts of redundant or false positive information, leading to difficulties in quantifying and prioritizing potential relationships between these DAGs and human diseases. In this study, a network-oriented gene entropy approach (NOGEA) is proposed for accurately inferring master genes that contribute to specific diseases by quantitatively calculating their perturbation abilities on directed disease-specific gene networks. In addition, we confirmed that the master genes identified by NOGEA have a high reliability for predicting disease-specific initiation events and progression risk. Master genes may also be used to extract the underlying information of different diseases, thus revealing mechanisms of disease comorbidity. More importantly, approved therapeutic targets are topologically localized in a small neighborhood of master genes in the interactome network, which provides a new way for predicting drug-disease associations. Through this method, 11 old drugs were newly identified and predicted to be effective for treating pancreatic cancer and then validated by in vitro experiments. Collectively, the NOGEA was useful for identifying master genes that control disease initiation and co-occurrence, thus providing a valuable strategy for drug efficacy screening and repositioning. NOGEA codes are publicly available at https://github.com/guozihuaa/NOGEA.     

Hypertension and Type 2 Diabetes Comorbidity in Adults in the United States: Risk of Overall and Regional Adiposity

Objective: To evaluate the impact of generalized, abdominal, and truncal fat deposits on the risk of hypertension and/or diabetes and to determine whether ethnic differences in these fat patterns are independently associated with increased risk for the hypertension–diabetes comorbidity (HDC). Research Methods and Procedures: Data (n = 7075) from the Third U.S. National Health and Nutrition Examination Survey were used for this investigation. To assess risks of hypertension and/or diabetes that were due to different fat patterns, odds ratios of men and women with various cut‐points of adiposities were compared with normal subjects in logistic regression models, adjusting for age, smoking, and alcohol intake. To evaluate the contribution of ethnic differences in obesity to the risks of HDC, we compared blacks and Hispanics with whites. Results: Generalized and abdominal obesities were independently associated with increased risk of hypertension, diabetes and HDC in white, black, and Hispanic men and women. The risk of HDC due to generalized, truncal, and abdominal obesities tended to be higher in whites than blacks and Hispanics. In men, the contribution of black and Hispanic ethnicities to the increased risk of HDC due to the various obesity phenotypes was ～73% and ～61%, respectively. The corresponding values for black and Hispanic women were ～115% and ～125%, respectively. Conclusions: In addition to advocating behavioral lifestyles to curb the epidemic of obesity among at‐risk populations in the United States, there is also the need for primary health care practitioners to craft their advice to the degree and type of obesity in these at‐risk groups.     

A systems‐level approach to mapping the telomere length maintenance gene circuitry

The ends of eukaryotic chromosomes are protected by telomeres, nucleoprotein structures that are essential for chromosomal stability and integrity. Understanding how telomere length is controlled has significant medical implications, especially in the fields of aging and cancer. Two recent systematic genome‐wide surveys measuring the telomere length of deleted mutants in the yeast Saccharomyces cerevisiae have identified hundreds of telomere length maintenance (TLM) genes, which span a large array of functional categories and different localizations within the cell. This study presents a novel general method that integrates large‐scale screening mutant data with protein–protein interaction information to rigorously chart the cellular subnetwork underlying the function investigated. Applying this method to the yeast telomere length control data, we identify pathways that connect the TLM proteins to the telomere‐processing machinery, and predict new TLM genes and their effect on telomere length. We experimentally validate some of these predictions, demonstrating that our method is remarkably accurate. Our results both uncover the complex cellular network underlying TLM and validate a new method for inferring such networks.     

Mood induction in depressive patients: a comparative multidimensional approach

Anhedonia, reduced positive affect and enhanced negative affect are integral characteristics of major depressive disorder (MDD). Emotion dysregulation, e.g. in terms of different emotion processing deficits, has consistently been reported. The aim of the present study was to investigate mood changes in depressive patients using a multidimensional approach for the measurement of emotional reactivity to mood induction procedures. Experimentally, mood states can be altered using various mood induction procedures. The present study aimed at validating two different positive mood induction procedures in patients with MDD and investigating which procedure is more effective and applicable in detecting dysfunctions in MDD. The first procedure relied on the presentation of happy vs. neutral faces, while the second used funny vs. neutral cartoons. Emotional reactivity was assessed in 16 depressed and 16 healthy subjects using self-report measures, measurements of electrodermal activity and standardized analyses of facial responses. Positive mood induction was successful in both procedures according to subjective ratings in patients and controls. In the cartoon condition, however, a discrepancy between reduced facial activity and concurrently enhanced autonomous reactivity was found in patients. Relying on a multidimensional assessment technique, a more comprehensive estimate of dysfunctions in emotional reactivity in MDD was available than by self-report measures alone and this was unsheathed especially by the mood induction procedure relying on cartoons. The divergent facial and autonomic responses in the presence of unaffected subjective reactivity suggest an underlying deficit in the patients' ability to express the felt arousal to funny cartoons. Our results encourage the application of both procedures in functional imaging studies for investigating the neural substrates of emotion dysregulation in MDD patients. Mood induction via cartoons appears to be superior to mood induction via faces and autobiographical material in uncovering specific emotional dysfunctions in MDD.     

Evolving modular genetic regulatory networks with a recursive,top-down approach

Being able to design genetic regulatory networks (GRNs) to achieve a desired cellular function is one of the main goals of synthetic biology. However, determining minimal GRNs that produce desired time-series behaviors is non-trivial. In this paper, we propose a ‘top-down’ approach to evolving small GRNs and then use these to recursively boot-strap the identification of larger, more complex, modular GRNs. We start with relatively dense GRNs and then use differential evolution (DE) to evolve interaction coefficients. When the target dynamical behavior is found embedded in a dense GRN, we narrow the focus of the search and begin aggressively pruning out excess interactions at the end of each generation. We first show that the method can quickly rediscover known small GRNs for a toggle switch and an oscillatory circuit. Next we include these GRNs as non-evolvable subnetworks in the subsequent evolution of more complex, modular GRNs. Successful solutions found in canonical DE where we truncated small interactions to zero, with or without an interaction penalty term, invariably contained many excess interactions. In contrast, by incorporating aggressive pruning and the penalty term, the DE was able to find minimal or nearly minimal GRNs in all test problems.     

Autoregressive models for capture-recapture data: a Bayesian approach

In this article, we incorporate an autoregressive time-series framework into models for animal survival using capture-recapture data. Researchers modeling animal survival probabilities as the realization of a random process have typically considered survival to be independent from one time period to the next. This may not be realistic for some populations. Using a Gibbs sampling approach, we can estimate covariate coefficients and autoregressive parameters for survival models. The procedure is illustrated with a waterfowl band recovery dataset for northern pintails (Anas acuta). The analysis shows that the second lag autoregressive coefficient is significantly less than 0, suggesting that there is a triennial relationship between survival probabilities and emphasizing that modeling survival rates as independent random variables may be unrealistic in some cases. Software to implement the methodology is available at no charge on the Internet.     

Comorbidity between major depressive disorder and physical diseases: a comprehensive review of epidemiology,mechanisms and management

Populations with common physical diseases – such as cardiovascular diseases, cancer and neurodegenerative disorders – experience substantially higher rates of major depressive disorder (MDD) than the general population. On the other hand, people living with MDD have a greater risk for many physical diseases. This high level of comorbidity is associated with worse outcomes, reduced adherence to treatment, increased mortality, and greater health care utilization and costs. Comorbidity can also result in a range of clinical challenges, such as a more complicated therapeutic alliance, issues pertaining to adaptive health behaviors, drug-drug interactions and adverse events induced by medications used for physical and mental disorders. Potential explanations for the high prevalence of the above comorbidity involve shared genetic and biological pathways. These latter include inflammation, the gut microbiome, mitochondrial function and energy metabolism, hypothalamic-pituitary-adrenal axis dysregulation, and brain structure and function. Furthermore, MDD and physical diseases have in common several antecedents related to social factors (e.g., socioeconomic status), lifestyle variables (e.g., physical activity, diet, sleep), and stressful live events (e.g., childhood trauma). Pharmacotherapies and psychotherapies are effective treatments for comorbid MDD, and the introduction of lifestyle interventions as well as collaborative care models and digital technologies provide promising strategies for improving management. This paper aims to provide a detailed overview of the epidemiology of the comorbidity of MDD and specific physical diseases, including prevalence and bidirectional risk; of shared biological pathways potentially implicated in the pathogenesis of MDD and common physical diseases; of socio-environmental factors that serve as both shared risk and protective factors; and of management of MDD and physical diseases, including prevention and treatment. We conclude with future directions and emerging research related to optimal care of people with comorbid MDD and physical diseases.     

In silico insights into the identification of potential novel angiogenic inhibitors against human VEGFR-2: a new SAR-based hierarchical clustering approach

Abstract

In this study, binding efficiency of new pyrrolopyrimidine structural analogs against human vascular endothelial growth factor receptor-2 (VEGFR-2) were elucidated using integrated in silico methods. Optimized high-resolution model of VEGFR-2 was generated and adopted for structure-based virtual screening approaches. Pyrrolopyrimidine inhibitor (CP15) associated compounds were screened from PubChem database and subjected to virtual screening and comparative docking methods against the receptor ligand-binding domain. Accordingly, high efficient compounds were clustered with similarity indices through PubChem structure cluster module using single-linkage algorithm. Moreover, pharmacokinetics including drug metabolism activities of high-binding leads under investigation was portrayed using ADMET and similarity ensemble analysis. Optimal energy orientations of the selected protein model have been shown to be reliable, and highly recommended for screening and docking studies. Docking and clustering strategies were shown that nineteen candidates as most effective binders for VEGFR-2 than CP15, and are grouped as three classes. Lys⁸⁶⁸, Glu⁸⁸⁵, Cys⁹¹⁹, His¹⁰²⁶, Arg¹⁰²⁷, Asp¹⁰⁴⁶, and Gly¹⁰⁴⁸ residues were predicted as novel hotspot residues, and participate in H-bonds, π-cation, π-stacking, halogen bonds, and salt-bridges formation with ligands. These additional bonds are contributing extent stability that holds the receptor structure at flexible state, this make difficult to any further conformational changes for evoking angiogenic signals. The ADMET and similarity ensemble analysis results were strongly indicated that thirteen candidates as best ligands for angiogenesis targets. Altogether, these findings indicate potential angiogenic templates and their binding levels with VEGFR-2; sorted viewpoints could be useful as a promising way to describe potential angiogenesis inhibitors with related molecular targets.     

Binary regression with misclassified response and covariate subject to measurement error: a bayesian approach

We consider a Bayesian analysis for modeling a binary response that is subject to misclassification. Additionally, an explanatory variable is assumed to be unobservable, but measurements are available on its surrogate. A binary regression model is developed to incorporate the measurement error in the covariate as well as the misclassification in the response. Unlike existing methods, no model parameters need be assumed known. Markov chain Monte Carlo methods are utilized to perform the necessary computations. The methods developed are illustrated using atomic bomb survival data. A simulation experiment explores advantages of the approach.     

Evaluating interventions in health: a reconciliatory approach

Health-related Quality of Life measures have recently been attacked from two directions, both of which criticize the preference-based method of evaluating health states they typically incorporate. One attack, based on work by Daniel Kahneman and others, argues that 'experience' is a better basis for evaluation. The other, inspired by Amartya Sen, argues that 'capability' should be the guiding concept. In addition, opinion differs as to whether health evaluation measures are best derived from consultations with the general public, with patients, or with health professionals. And there is disagreement about whether these opinions should be solicited individually and aggregated, or derived instead from a process of collective deliberation. These distinctions yield a wide variety of possible approaches, with potentially differing policy implications. We consider some areas of disagreement between some of these approaches. We show that many of the perspectives seem to capture something important, such that it may be a mistake to reject any of them. Instead we suggest that some of the existing 'instruments' designed to measure HR QoLs may in fact successfully already combine these attributes, and with further refinement such instruments may be able to provide a reasonable reconciliation between the perspectives.     

Standardizing the analysis of conditioned fear in rodents: a multidimensional software approach

Data comparability between different laboratories strongly depends on the individually applied analysis method. This factor is often a critical source of variation in rodent phenotyping and has never been systematically investigated in Pavlovian fear conditioning paradigms. In rodents, fear is typically quantified in terms of freezing duration via manual observation or automated systems. While manual analysis includes biases such as tiredness or inter‐personal scoring variability, computer‐assisted systems are unable to distinguish between freezing and immobility. Consequently, the novel software called MOVE follows a semi‐automatized approach that prefilters video sequences of interest for the final human judgment. Furthermore, MOVE allows integrating additional data sources (e.g. force‐sensitive platform, EEG) to reach the most accurate and precise results. MOVE directly supports multi‐angle video recordings with webcams or standard laboratory equipment. The integrated manual key logger and internal video player complement this all‐in‐one software solution. Calculating the interlaboratory variability of manual freezing evaluation revealed significantly different freezing scores in two out of six laboratories. This difference was minimized when all experiments were analyzed with MOVE. Applied to a genetically modified mouse model, MOVE revealed higher fear responses of CB1 deficient mice compared to their wild‐type littermates after foreground context fear conditioning. Multi‐angle video analysis compared to the single‐camera approach reached up to 15% higher accuracy and two fold higher precision. Multidimensional analysis provided by integration of additional data sources further improved the overall result. We conclude that the widespread usage of MOVE could substantially improve the comparability of results from different laboratories.     

In silico model‐based inference: A contemporary approach for hypothesis testing in network biology

Inductive inference plays a central role in the study of biological systems where one aims to increase their understanding of the system by reasoning backwards from uncertain observations to identify causal relationships among components of the system. These causal relationships are postulated from prior knowledge as a hypothesis or simply a model. Experiments are designed to test the model. Inferential statistics are used to establish a level of confidence in how well our postulated model explains the acquired data. This iterative process, commonly referred to as the scientific method, either improves our confidence in a model or suggests that we revisit our prior knowledge to develop a new model. Advances in technology impact how we use prior knowledge and data to formulate models of biological networks and how we observe cellular behavior. However, the approach for model‐based inference has remained largely unchanged since Fisher, Neyman and Pearson developed the ideas in the early 1900s that gave rise to what is now known as classical statistical hypothesis (model) testing. Here, I will summarize conventional methods for model‐based inference and suggest a contemporary approach to aid in our quest to discover how cells dynamically interpret and transmit information for therapeutic aims that integrates ideas drawn from high performance computing, Bayesian statistics, and chemical kinetics. © 2014 American Institute of Chemical Engineers Biotechnol. Prog., 30:1247–1261, 2014     

Elderly complete denture wearers: a social approach to tooth loss

doi: 10.1111/j.1741‐2358.2011.00550.x Elderly complete denture wearers: a social approach to tooth loss Objectives: To correlate emotional reactions to tooth loss with denture satisfaction attributes in elderly complete denture wearers. Background: Total tooth loss is a serious life event, and poor oral health has an impact on daily life. Edentulism treated by rehabilitation with dentures can have a positive effect on patients’ self‐image and social behaviour. Methods: A group of 80 edentulous subjects undergoing routine prosthetic care in a Greek Department of Prosthetic Dentistry were interviewed using two structured questionnaires. The first questionnaire explored reactions to tooth loss, whereas the second measured their subjective experience of complete dentures. The responses to both questionnaires were compared using the statistical package spss v.17. Results: The results showed significant correlation between aspects of tooth loss experience and complete denture satisfaction. Despite the fact that a substantial proportion of patients were satisfied with their complete dentures, some patients experienced increased social and psychological problems related to their edentulousness and the wearing of complete dentures. The aesthetic and functional aspects of complete dentures affected both patients’ social behaviour and self‐confidence. Conclusions: Total tooth loss was not only reflected in patient’s social behaviour and self‐image, but it had a complex and multifaceted impact on satisfaction from complete dentures.     

The joint allele frequency spectrum of multiple populations: a coalescent theory approach

The allele frequency spectrum is a series of statistics that describe genetic polymorphism, and is commonly used for inferring population genetic parameters and detecting natural selection. Population genetic theory on the allele frequency spectrum for a single population has been well studied using both coalescent theory and diffusion equations. Recently, the theory was extended to the joint allele frequency spectrum (JAFS) for three populations using diffusion equations and was shown to be very useful in inferring human demographic history. In this paper, I show that the JAFS can be analytically derived with coalescent theory for a basic model of two isolated populations and then extended to multiple populations and various complex scenarios, such as those involving population growth and bottleneck, migration, and positive selection. Simulation study is used to demonstrate the accuracy and applicability of the theoretical model. The coalescent theory-based approach for the JAFS can characterize the demographic history with comprehensive statistical models as the diffusion approach does, and in addition gains several novel advantages: the computational complexity of calculating the JAFS with coalescent theory is reduced, and thus it is feasible to analytically obtain the JAFS for multiple populations; the hitchhiking effect can be efficiently modeled in coalescent theory, enabling the development of methodologies for detecting selection via multi-population polymorphism data. As an alternative to the diffusion approximation approach, the coalescent theory for the JAFS also provides a foundation for population genetic inference with the advent of large-scale genomic polymorphism data.     

Parameter estimation in a model for multidimensional recording of neuronal data: a Gibbsian approximation approach

This article proposes improved numerical procedures for estimating parameters in a spatiotemporal lattice model introduced for the analysis of cortical activities monitored from arrays of diodes. The numerical algorithms are based on approximations inspired by statistical physics. Both Gibbsian and mean-field approximations are used; they allow for computing local conditional probabilities inside the lattice. The statistical procedures rely on the computation of pseudomaximum-likelihood estimators. The estimators are evaluated on the basis of Monte Carlo simulations. These simulations show that mean-field approximations are useful for reducing the variance of estimators when the data are recorded from arrays of 144 diodes (which are in accordance with standard practice). In light of these improved methods, we give new interpretations for a data set obtained from optical recording of a Guinea pig's auditory cortex in response to pure tone stimulations.     

The epidemiology and iatrogenic transmission of hepatitis C virus in Egypt: a Bayesian coalescent approach

Hepatitis C virus (HCV) is a leading cause of liver cancer and cirrhosis, and Egypt has possibly the highest HCV prevalence worldwide. In this article we use a newly developed Bayesian inference framework to estimate the transmission dynamics of HCV in Egypt from sampled viral gene sequences, and to predict the public health impact of the virus. Our results indicate that the effective number of HCV infections in Egypt underwent rapid exponential growth between 1930 and 1955. The timing and speed of this spread provides quantitative genetic evidence that the Egyptian HCV epidemic was initiated and propagated by extensive antischistosomiasis injection campaigns. Although our results show that HCV transmission has since decreased, we conclude that HCV is likely to remain prevalent in Egypt for several decades. Our combined population genetic and epidemiological analysis provides detailed estimates of historical changes in Egyptian HCV prevalence. Because our results are consistent with a demographic scenario specified a priori, they also provide an objective test of inference methods based on the coalescent process.     

Modeling a Poisson forest in variable elevations: a nonparametric Bayesian approach

A nonparametric Bayesian formulation is given to the problem of modeling nonhomogeneous spatial point patterns influenced by concomitant variables. Only incomplete information on the concomitant variables is assumed, consisting of a relatively small number of point measurements. Residual variation, caused by other unmeasured influential factors, is modeled in terms of a spatially varying baseline intensity function. A Markov chain Monte Carlo scheme is proposed for the simultaneous nonparametric estimation of each unknown function in the model. The suggested method is illustrated by reanalysing a data set in Rathbun (1996, Biometrics 52, 226-242), and the estimated models are compared with those obtained by Rathbun.