NETWORK TOPOLOGY IN HUMAN PROTEIN INTERACTION DATA PREDICTS FUNCTIONAL ASSOCIATION

High-throughput assays, such as yeast two-hybrid system, have generated a huge amount of protein-protein interaction (PPI) data in the past decade. This tremendously increases the need for developing reliable methods to systematically and automatically suggest protein functions and relationships between them. With the available PPI data, it is now possible to study the functions and relationships in the context of a large-scale network. To data, several network-based schemes have been provided to effectively annotate protein functions on a large scale. However, due to those inherent noises in high-throughput data generation, new methods and algorithms should be developed to increase the reliability of functional annotations. Previous work in a yeast PPI network (Samanta and Liang, 2003) has shown that the local connection topology, particularly for two proteins sharing an unusually large number of neighbors, can predict functional associations between proteins, and hence suggest their functions. One advantage of the work is that their algorithm is not sensitive to noises (false positives) in high-throughput PPI data. In this study, we improved their prediction scheme by developing a new algorithm and new methods which we applied on a human PPI network to make a genome-wide functional inference. We used the new algorithm to measure and reduce the influence of hub proteins on detecting functionally associated proteins. We used the annotations of the Gene Ontology (GO) and the Kyoto Encyclopedia of Genes and Genomes (KEGG) as independent and unbiased benchmarks to evaluate our algorithms and methods within the human PPI network. We showed that, compared with the previous work from Samanta and Liang, our algorithm and methods developed in this study improved the overall quality of functional inferences for human proteins. By applying the algorithms to the human PPI network, we obtained 4,233 significant functional associations among 1,754 proteins. Further comparisons of their KEGG and GO annotations allowed us to assign 466 KEGG pathway annotations to 274 proteins and 123 GO annotations to 114 proteins with estimated false discovery rates of <21% for KEGG and <30% for GO. We clustered 1,729 proteins by their functional associations and made pathway analysis to identify several subclusters that are highly enriched in certain signaling pathways. Particularly, we performed a detailed analysis on a subcluster enriched in the transforming growth factor β signaling pathway (P<10-50) which is important in cell proliferation and tumorigenesis. Analysis of another four subclusters also suggested potential new players in six signaling pathways worthy of further experimental investigations. Our study gives clear insight into the common neighbor-based prediction scheme and provides a reliable method for large-scale functional annotations in this post-genomic era.

The interorganizational implementation of state health policy: A network approach

This is an implementation analysis of three consecutive state health policies whose goal was to improve access to maternal and child health services in Texas from 1983 to 1986. Of particular interest is the choice of the unit of analysis, the policy subsystem, and the network approach to analysis. The network approach analyzes and compares the structure and decision process of six policy subsystems in order to explain program performance. Both changes in state health policy as well as differences in implementation contexts explain evolution of the program administrative and service unit, the policy subsystem. And, in turn, the evolution of the policy subsystem explains changes in program performance. ^

Extended kalman filter for estimation of parameters in nonlinear state-space models of biochemical networks.

It is system dynamics that determines the function of cells, tissues and organisms. To develop mathematical models and estimate their parameters are an essential issue for studying dynamic behaviors of biological systems which include metabolic networks, genetic regulatory networks and signal transduction pathways, under perturbation of external stimuli. In general, biological dynamic systems are partially observed. Therefore, a natural way to model dynamic biological systems is to employ nonlinear state-space equations. Although statistical methods for parameter estimation of linear models in biological dynamic systems have been developed intensively in the recent years, the estimation of both states and parameters of nonlinear dynamic systems remains a challenging task. In this report, we apply extended Kalman Filter (EKF) to the estimation of both states and parameters of nonlinear state-space models. To evaluate the performance of the EKF for parameter estimation, we apply the EKF to a simulation dataset and two real datasets: JAK-STAT signal transduction pathway and Ras/Raf/MEK/ERK signaling transduction pathways datasets. The preliminary results show that EKF can accurately estimate the parameters and predict states in nonlinear state-space equations for modeling dynamic biochemical networks.

TIME SERIES ANALYSIS AS INPUT FOR PREDICTIVE MODELING: PREDICTING CARDIAC ARREST IN A PEDIATRIC INTENSIVE CARE UNIT

The first manuscript, entitled "Time-Series Analysis as Input for Clinical Predictive Modeling: Modeling Cardiac Arrest in a Pediatric ICU" lays out the theoretical background for the project. There are several core concepts presented in this paper. First, traditional multivariate models (where each variable is represented by only one value) provide single point-in-time snapshots of patient status: they are incapable of characterizing deterioration. Since deterioration is consistently identified as a precursor to cardiac arrests, we maintain that the traditional multivariate paradigm is insufficient for predicting arrests. We identify time series analysis as a method capable of characterizing deterioration in an objective, mathematical fashion, and describe how to build a general foundation for predictive modeling using time series analysis results as latent variables. Building a solid foundation for any given modeling task involves addressing a number of issues during the design phase. These include selecting the proper candidate features on which to base the model, and selecting the most appropriate tool to measure them. We also identified several unique design issues that are introduced when time series data elements are added to the set of candidate features. One such issue is in defining the duration and resolution of time series elements required to sufficiently characterize the time series phenomena being considered as candidate features for the predictive model. Once the duration and resolution are established, there must also be explicit mathematical or statistical operations that produce the time series analysis result to be used as a latent candidate feature. In synthesizing the comprehensive framework for building a predictive model based on time series data elements, we identified at least four classes of data that can be used in the model design. The first two classes are shared with traditional multivariate models: multivariate data and clinical latent features. Multivariate data is represented by the standard one value per variable paradigm and is widely employed in a host of clinical models and tools. These are often represented by a number present in a given cell of a table. Clinical latent features derived, rather than directly measured, data elements that more accurately represent a particular clinical phenomenon than any of the directly measured data elements in isolation. The second two classes are unique to the time series data elements. The first of these is the raw data elements. These are represented by multiple values per variable, and constitute the measured observations that are typically available to end users when they review time series data. These are often represented as dots on a graph. The final class of data results from performing time series analysis. This class of data represents the fundamental concept on which our hypothesis is based. The specific statistical or mathematical operations are up to the modeler to determine, but we generally recommend that a variety of analyses be performed in order to maximize the likelihood that a representation of the time series data elements is produced that is able to distinguish between two or more classes of outcomes. The second manuscript, entitled "Building Clinical Prediction Models Using Time Series Data: Modeling Cardiac Arrest in a Pediatric ICU" provides a detailed description, start to finish, of the methods required to prepare the data, build, and validate a predictive model that uses the time series data elements determined in the first paper. One of the fundamental tenets of the second paper is that manual implementations of time series based models are unfeasible due to the relatively large number of data elements and the complexity of preprocessing that must occur before data can be presented to the model. Each of the seventeen steps is analyzed from the perspective of how it may be automated, when necessary. We identify the general objectives and available strategies of each of the steps, and we present our rationale for choosing a specific strategy for each step in the case of predicting cardiac arrest in a pediatric intensive care unit. Another issue brought to light by the second paper is that the individual steps required to use time series data for predictive modeling are more numerous and more complex than those used for modeling with traditional multivariate data. Even after complexities attributable to the design phase (addressed in our first paper) have been accounted for, the management and manipulation of the time series elements (the preprocessing steps in particular) are issues that are not present in a traditional multivariate modeling paradigm. In our methods, we present the issues that arise from the time series data elements: defining a reference time; imputing and reducing time series data in order to conform to a predefined structure that was specified during the design phase; and normalizing variable families rather than individual variable instances. The final manuscript, entitled: "Using Time-Series Analysis to Predict Cardiac Arrest in a Pediatric Intensive Care Unit" presents the results that were obtained by applying the theoretical construct and its associated methods (detailed in the first two papers) to the case of cardiac arrest prediction in a pediatric intensive care unit. Our results showed that utilizing the trend analysis from the time series data elements reduced the number of classification errors by 73%. The area under the Receiver Operating Characteristic curve increased from a baseline of 87% to 98% by including the trend analysis. In addition to the performance measures, we were also able to demonstrate that adding raw time series data elements without their associated trend analyses improved classification accuracy as compared to the baseline multivariate model, but diminished classification accuracy as compared to when just the trend analysis features were added (ie, without adding the raw time series data elements). We believe this phenomenon was largely attributable to overfitting, which is known to increase as the ratio of candidate features to class examples rises. Furthermore, although we employed several feature reduction strategies to counteract the overfitting problem, they failed to improve the performance beyond that which was achieved by exclusion of the raw time series elements. Finally, our data demonstrated that pulse oximetry and systolic blood pressure readings tend to start diminishing about 10-20 minutes before an arrest, whereas heart rates tend to diminish rapidly less than 5 minutes before an arrest.

Bifurcation and singularity analysis of a molecular network for the induction of long-term memory.

Withdrawal reflexes of the mollusk Aplysia exhibit sensitization, a simple form of long-term memory (LTM). Sensitization is due, in part, to long-term facilitation (LTF) of sensorimotor neuron synapses. LTF is induced by the modulatory actions of serotonin (5-HT). Pettigrew et al. developed a computational model of the nonlinear intracellular signaling and gene network that underlies the induction of 5-HT-induced LTF. The model simulated empirical observations that repeated applications of 5-HT induce persistent activation of protein kinase A (PKA) and that this persistent activation requires a suprathreshold exposure of 5-HT. This study extends the analysis of the Pettigrew model by applying bifurcation analysis, singularity theory, and numerical simulation. Using singularity theory, classification diagrams of parameter space were constructed, identifying regions with qualitatively different steady-state behaviors. The graphical representation of these regions illustrates the robustness of these regions to changes in model parameters. Because persistent protein kinase A (PKA) activity correlates with Aplysia LTM, the analysis focuses on a positive feedback loop in the model that tends to maintain PKA activity. In this loop, PKA phosphorylates a transcription factor (TF-1), thereby increasing the expression of an ubiquitin hydrolase (Ap-Uch). Ap-Uch then acts to increase PKA activity, closing the loop. This positive feedback loop manifests multiple, coexisting steady states, or multiplicity, which provides a mechanism for a bistable switch in PKA activity. After the removal of 5-HT, the PKA activity either returns to its basal level (reversible switch) or remains at a high level (irreversible switch). Such an irreversible switch might be a mechanism that contributes to the persistence of LTM. The classification diagrams also identify parameters and processes that might be manipulated, perhaps pharmacologically, to enhance the induction of memory. Rational drug design, to affect complex processes such as memory formation, can benefit from this type of analysis.

Social network size and perceived network sufficiency as predictors of perceived health status

This study focused on the relationship between social network size (number of friends and relatives), perceived sufficiency of the network and self-rated health utilizing data from the National Survey of Personal Health Practices and Consequences, 1979. For men neither perceived sufficiency nor number of relatives were associated with self-rated health status. The number of friends was positively associated with health status. For women perceived network sufficiency was positively and significantly related to health status, independent of network size. The number of friends and relatives was not associated with self-rated health status. The sociodemographic variables accounted for most of the explained variance in health status for both males and females. Social networks may hold different meanings for women and men, and may require qualitative as well as quantitative analysis. There may have been insufficient variance in the major variables to produce meaningful results. ^

Network-based analysis of genome-wide scans of natural selection

The genomic era brought by recent advances in the next-generation sequencing technology makes the genome-wide scans of natural selection a reality. Currently, almost all the statistical tests and analytical methods for identifying genes under selection was performed on the individual gene basis. Although these methods have the power of identifying gene subject to strong selection, they have limited power in discovering genes targeted by moderate or weak selection forces, which are crucial for understanding the molecular mechanisms of complex phenotypes and diseases. Recent availability and rapid completeness of many gene network and protein-protein interaction databases accompanying the genomic era open the avenues of exploring the possibility of enhancing the power of discovering genes under natural selection. The aim of the thesis is to explore and develop normal mixture model based methods for leveraging gene network information to enhance the power of natural selection target gene discovery. The results show that the developed statistical method, which combines the posterior log odds of the standard normal mixture model and the Guilt-By-Association score of the gene network in a naïve Bayes framework, has the power to discover moderate/weak selection gene which bridges the genes under strong selection and it helps our understanding the biology under complex diseases and related natural selection phenotypes.^

Saturday Morning Television Advertisements Aired on English and Spanish Language Networks along the Texas-Mexico Border

Objectives: The aim of this content analysis study is to characterize the TV advertisements aired to an at-risk child population along the Texas-Mexico border. Methods: We characterized the early Saturday morning TV advertisements aired by three broadcast network categories (U.S. English language, U.S. Spanish language, and Mexican Spanish language) in Spring 2010. The number, type (food related vs. non-food related), target audience, and persuasion tactics used were recorded. Advertised foods, based on nutrition content, were categorized as meeting or not meeting current dietary guidelines. Results: Most commercials were non-food related (82.7%, 397 of 480). The majority of the prepared foods (e.g., cereals, snacks, and drinks) advertised did not meet the current U.S. Dietary Guidelines. Additionally, nutrition content information was not available for many of the foods advertised on the Mexican Spanish language broadcast network category. Conclusions: For U.S. children at risk for obesity along the Texas-Mexico border exposure to TV food advertisements may result in the continuation of sedentary behavior as well as an increased consumption of foods of poor nutritional quality. An international regulatory effort to monitor and enforce the reduction of child-oriented food advertising is needed. Editors' Note: This article was submitted in response to the first issue of the Journal of Applied Research on Children: Latino Children.

B cell memory and idiotype network regulation

This investigation examined the clonal dynamics of B-cell expression and evaluated the role of idiotype network interactions in shaping the expressed secondary B-cell repertoire. Three interrelated experimental approaches were applied. The first approach was designed to distinguish between regulatory influences controlled by the major histocompatibility complex (MHC) and regulatory influences controlled by non-MHC factors including the idiotype network. This approach consisted of studies on the clonal dynamics and heterogeneity of the expressed IgG antibody repertoire of BALB/c mice. The second approach involved the analysis of the clonal dynamics of antibody responses of outbred rabbits. This analysis was coupled with studies to detect the occurrence and activity of constituents of the idiotype network. In the third approach the transfer of rabbit lymphocytes from immunized donors to MHC matched naive recipients was used to examine the effects of recipient non-MHC immunoregulatory influences on the expression of donor memory B-cells. Although many memory B cells were unaffected by non-MHC influences, these data show that non-MHC immunoregulatory influences can affect the expression of B-cells in the secondary response of inbred mice and outbred rabbits. The results also indicate that most IgG antibody responses are heterogeneous and are characterized by a stable group of dominant clonotypes. Clonal dominance and B-cell memory were found to be established early in an immune response. The expression of B memory clones appeared to be favored over the expression of virgin B cells. The injection of anti-tetanus antibody induced the antigen independent production of anti-tetanus antibody, probably through idiotypic mechanisms. These results demonstrate that both antibody and antigen can affect the expressed B-ceIl repertoire. Thus, idiotypic interactions are capable of influencing the expression of B-cells and these findings support the existence and function of an idiotype network with strong immunoregulatory potential. ^

Convergent and divergent modulatory pathways in {\it Aplysia\/}: Cellular and network studies

Neuromodulation is essential to many functions of the nervous system. In the simple gastropod mollusk Aplysia californica, neuromodulation of the circuits for the defensive withdrawal reflexes has been associated with several forms of learning. In the present work, the neurotransmitters and neural circuitry which contribute to the modulation of the tail-siphon withdrawal reflex were examined.^ A recently-identified neuropeptide transmitter, buccalin A was found to modulate the biophysical properties of the sensory neurons that mediate the reflex. The actions of buccalin A on the sensory neurons were compared with those of the well-characterized modulatory transmitter serotonin, and convergence and divergence in the actions of these two transmitters were evaluated. Buccalin A dramatically increased the excitability of sensory neurons and occluded further enhancement of excitability by serotonin. Buccalin A produced no significant change in spike duration, and it did not block serotonin-induced spike broadening. Voltage-clamp analysis revealed the currents that may be involved in the effects on spike duration and excitability. Buccalin A decreased an outward current similar to the S-K$\sp+$ current (I$\sb{\rm K,S}$). Buccalin A appeared to occlude further modulation of I$\sb{\rm K,S}$ by serotonin, but did not block serotonin-induced modulation of the voltage-dependent delayed rectifier K$\sp+$ current (I$\sb{\rm K,V}$). These results suggest that buccalin A converges on some, but not all, of the same subcellular modulatory pathways as serotonin.^ In order to begin to understand neuromodulation in a more physiological context for the tail-siphon withdrawal reflex, the modulatory circuitry for the tail-withdrawal circuit was examined. Mechanoafferent neurons in the J cluster of the cerebral ganglion were identified as elements of a modulatory circuit for the reflex. Excitatory and inhibitory connections were observed between the J cells and the pleural sensory neurons, the tail motor neurons, and several classes of interneurons for the tail-siphon withdrawal circuit. The J cells produced both fast and slow PSPs in these neurons. Of particular interest was the ability of the J cells to produce slow EPSPs in the pleural sensory neurons. These slow EPSPs were associated with an increase in the excitability of the sensory neurons. The J cells appear to mediate both sensory and modulatory inputs to the circuit for the tail-siphon withdrawal reflex from the anterior part of the animal. ^

Functional and motor neuronal analysis of defensive siphon responses in {\it Aplysia californica}

One of the central goals of neuroscience research is to determine how networks of neurons control and modify behavior. One of the most influential model systems for this kind of analysis is the siphon and gill withdrawal reflex of the marine mollusc A. californica. In response to tactile stimulation, the siphon displays 3 different responses: (1) a posterior pointing and leveling (flaring) of the siphon in response to tail stimulation (the siphon T response), (2) constriction and anterior pointing to head stimulation (the siphon H response) and (3) constriction and withdrawal between the animal's parapodia (the siphon S response). The siphon S response is pseudoconditioned by a noxious tail stimulus to resemble the siphon T response. Behavioral and combined behavioral/intracellular studies were conducted to determine the motor neuronal control of these behaviors and to search for mechanisms of siphon response transformation following pseudoconditioning. The present studies have found that the flaring component of pseudoconditioned siphon S responses occurs during mantle pumping (MP) triggered by noxious tail stimulation. Siphon stimulation also triggers MP, as recorded in neurons of the Interneuron II pattern generator which commands MP. The 4 LF$\rm\sb{SB}$ siphon motor neurons (SMNs) were found necessary and sufficient for the siphon T response, while SMNs RD$\rm\sb S$ and LD$\rm\sb{S1}$ were found necessary and sufficient for the siphon H response. Following pseudoconditioning, there is an increase in the number of evoked spikes to the test stimulus for the LF$\rm\sb{SB}$ cells and a decreased number for RD$\rm\sb S.$ Siphon flaring occurring during the pseudoconditioned response correlates with increased LF$\rm\sb{SB}$ activity during triggered MP cycles. This suggests that psuedoconditioning is in part due to reconfiguration of the motor outputs of the Interneuron II network. These results suggest that these defensive responses are controlled and patterned by a well-defined, finite set of motor neurons and interneurons (Interneuron II) that are dedicated to specific behavioral functions, but also have parallel distributed properties. ^

Achieving effective biosurveillance with the Nationwide Health Information Network

The federal government is currently developing the Nationwide Health Information Network (NHIN). Described as a “network of networks,” the NHIN seeks to provide a nationwide, interoperable health information infrastructure that will securely connect consumers with those involved in health care. As part of the national health information technology (HIT) agenda, the NHIN aims to improve individual and population health by enabling health information to follow the consumer, be available for clinical decision-making, and support important public health measures such as biosurveillance. While the NHIN promises to improve clinical care to individuals and to reduce U.S. health care system costs overall, this electronic environment presents novel challenges for protecting individually identifiable health information. A major barrier to achieving public trust in the NHIN is the development of, and adherence to, a consistent and coordinated approach to privacy and security of health information. This paper will analyze the policy framework for electronic health information exchange with the NHIN. This exercise will demonstrate that the current policy is an effective framework for achieving effective biosurveillance with the NHIN. ^

Relating one's social support network to one's mid term grade

Much has been written about the relation of social support to health outcomes. Support networks were found to be predictive of health status. Not so clear was the manner in which social support helped the individual to avoid health complications. Whereas some aspects of the support network were protective, others were burdensome. Duties to one's network could serve as a stressor and duties outside one's network might stress the support system itself. Exposure to one's network was associated with certain health risks while disruption in one's social support network was associated with other health risks.^ Many factors contributed to the impact of a social support network upon the individual member: the characteristics of the individual, the individual's role or position within the network, qualities of the network and duties or indebtedness of the individual to the network. This investigation considered the possibility that performance could serve as a stressor in a fashion similar to an exposure to a health hazard.^ Because the literature includes many examples of studies in which the subjects were college students, academic progress is a performance common to most subjects. A profile of the support networks of successful students was contrasted with those of less successful students in this correlational study.^ What was uncovered in this investigation was a very complex web of interrelated constructs. Most aspects of the social support network did not significantly predict academic performance. Only a limited number of characteristics were associated with academic success: the frequency of support, student age, the existence of a 'mentor' within one' s network, and the extent to which one received a predominant source of support. Other factors had a tendency to be negatively correlated with midterm grade, suggesting those factors may impede academic performance.^ Medical status did not predict grades, but was correlated with many aspects of the network. Disruptions in particular parts of one's network were correlated with particular health categories. In fact, disruption in social support was more predictive of academic outcomes than medical complications. Whereas the individual's values were related to the contributing factors, only the individual's satisfaction with certain aspects of the support network were predictive of higher midterm grades in a psychology class. Dissatisfaction was associated with lower grades, suggesting a disruptive effect within the network. Associations among the features of support networks which predicted academic progress were considered. ^

Network-based genome-wide association study of Crohn's disease

Genome-wide association studies (GWAS) have rapidly become a standard method for disease gene discovery. Many recent GWAS indicate that for most disorders, only a few common variants are implicated and the associated SNPs explain only a small fraction of the genetic risk. The current study incorporated gene network information into gene-based analysis of GWAS data for Crohn's disease (CD). The purpose was to develop statistical models to boost the power of identifying disease-associated genes and gene subnetworks by maximizing the use of existing biological knowledge from multiple sources. The results revealed that Markov random field (MRF) based mixture model incorporating direct neighborhood information from a single gene network is not efficient in identifying CD-related genes based on the GWAS data. The incorporation of solely direct neighborhood information might lead to the low efficiency of these models. Alternative MRF models looking beyond direct neighboring information are necessary to be developed in the future for the purpose of this study.^

A comparative state policy analysis of the adoption of birth defects surveillance legislation

Birth defects are the leading cause of infant mortality in the United States and are a major cause of lifetime disability. However, efforts to understand their causes have been hampered by a lack of population-specific data. During 1990–2004, 22 state legislatures responded to this need by proposing birth defects surveillance legislation (BDSL). The contrast between these states and those that did not pass BDSL provides an opportunity to better understand conditions associated with US public health policy diffusion. ^ This study identifies key state-specific determinants that predict: (1) the introduction of birth defects surveillance legislation (BDSL) onto states' formal legislative agenda, and (2) the successful adoption of these laws. Secondary aims were to interpret these findings in a theoretically sound framework and to incorporate evidence from three analytical approaches. ^ The study begins with a comparative case study of Texas and Oregon (states with divergent BDSL outcomes), including a review of historical documentation and content analysis of key informant interviews. After selecting and operationalizing explanatory variables suggested by the case study, Qualitative Comparative Analysis (QCA) was applied to publically available data to describe important patterns of variation among 37 states. Results from logistic regression were compared to determine whether the two methods produced consistent findings. ^ Themes emerging from the comparative case study included differing budgetary conditions and the significance of relationships within policy issue networks. However, the QCA and statistical analysis pointed to the importance of political parties and contrasting societal contexts. Notably, state policies that allow greater access to citizen-driven ballot initiatives were consistently associated with lower likelihood of introducing BDSL. ^ Methodologically, these results indicate that a case study approach, while important for eliciting valuable context-specific detail, may fail to detect the influence of overarching, systemic variables, such as party competition. However, QCA and statistical analyses were limited by a lack of existing data to operationalize policy issue networks, and thus may have downplayed the impact of personal interactions. ^ This study contributes to the field of health policy studies in three ways. First, it emphasizes the importance of collegial and consistent relationships among policy issue network members. Second, it calls attention to political party systems in predicting policy outcomes. Finally, a novel approach to interpreting state data in a theoretically significant manner (QCA) has been demonstrated.^