Modeling genomic regulatory networks with big data

A generic parallel framework for inferring large-scale gene regulatory networks from expression profiles: application to Alzheimer's disease network

The inference of large-scale gene regulatory networks is essential for understanding comprehensive interactions among genes. Most existing methods are limited to reconstructing networks with a few hundred nodes. Therefore, parallel computing paradigms must be leveraged to construct large networks. We propose a generic parallel framework that enables any existing method, without re-engineering, to infer large networks in parallel, guaranteeing quality output. The framework is tested on 15 inference methods (not limited to) employing in silico benchmarks and real-world large expression matrices, followed by qualitative and speedup assessment. The framework does not compromise the quality of the base serial inference method. We rank the candidate methods and use the top-performing method to infer an Alzheimer's Disease (AD) affected network from large expression profiles of a triple transgenic mouse model consisting of 45,101 genes. The resultant network is further explored to obtain hub genes that emerge functionally related to the disease. We partition the network into 41 modules and conduct pathway enrichment analysis, revealing that a good number of participating genes are collectively responsible for several brain disorders, including AD. Finally, we extract the interactions of a few known AD genes and observe that they are periphery genes connected to the network's hub genes. Availability: The R implementation of the framework is downloadable from https://github.com/Netralab/GenericParallelFramework.

Unraveling Plant-Pathogen Interactions in Cereals Using RNA-seq

Over the past two decades, there have been significant advancements in the realm of transcriptomics, or the study of genes and their expression. Modern RNA sequencing technologies and high-performance computing are creating a "big data" revolution that provides new opportunities to explore the interactions between cereals and pathogens that affect grain yield and food safety. These data are being used to annotate genes and gene variants, as well as identify differentially expressed genes and create global gene co-expression networks. Moreover, these data can unravel the complex interactions between pathogen and host and identify genes and pathways involved in these interactions. This information can then be used for disease mitigation and the development of crops with superior resistance.

Python for gene expression

Genome biology shows substantial progress in its analytical and computational part in the last decades. Differential gene expression is one of many computationally intense areas; it is largely developed under R programming language. Here we explain possible reasons for such dominance of R in gene expression data. Next, we discuss the prospects for Python to become competitive in this area of research in coming years. We indicate that Python can be used already in a field of a single cell differential gene expression. We pinpoint still missing parts in Python and possibilities for improvement.

Emergence of co-expression in gene regulatory networks

Transcriptomes are known to organize themselves into gene co-expression clusters or modules where groups of genes display distinct patterns of coordinated or synchronous expression across independent biological samples. The functional significance of these co-expression clusters is suggested by the fact that highly coexpressed groups of genes tend to be enriched in genes involved in common functions and biological processes. While gene co-expression is widely assumed to reflect close regulatory proximity, the validity of this assumption remains unclear. Here we use a simple synthetic gene regulatory network (GRN) model and contrast the resulting co-expression structure produced by these networks with their known regulatory architecture and with the co-expression structure measured in available human expression data. Using randomization tests, we found that the levels of co-expression observed in simulated expression data were, just as with empirical data, significantly higher than expected by chance. When examining the source of correlated expression, we found that individual regulators, both in simulated and experimental data, fail, on average, to display correlated expression with their immediate targets. However, highly correlated gene pairs tend to share at least one common regulator, while most gene pairs sharing common regulators do not necessarily display correlated expression. Our results demonstrate that widespread co-expression naturally emerges in regulatory networks, and that it is a reliable and direct indicator of active co-regulation in a given cellular context.

Emergence of informative higher scales in biological systems: a computational toolkit for optimal prediction and control

The biological sciences span many spatial and temporal scales in attempts to understand the function and evolution of complex systems-level processes, such as embryogenesis. It is generally assumed that the most effective description of these processes is in terms of molecular interactions. However, recent developments in information theory and causal analysis now allow for the quantitative resolution of this question. In some cases, macro-scale models can minimize noise and increase the amount of information an experimenter or modeler has about "what does what." This result has numerous implications for evolution, pattern regulation, and biomedical strategies. Here, we provide an introduction to these quantitative techniques, and use them to show how informative macro-scales are common across biology. Our goal is to give biologists the tools to identify the maximally-informative scale at which to model, experiment on, predict, control, and understand complex biological systems.

Brief introduction of medical database and data mining technology in big data era

Data mining technology can search for potentially valuable knowledge from a large amount of data, mainly divided into data preparation and data mining, and expression and analysis of results. It is a mature information processing technology and applies database technology. Database technology is a software science that researches manages, and applies databases. The data in the database are processed and analyzed by studying the underlying theory and implementation methods of the structure, storage, design, management, and application of the database. We have introduced several databases and data mining techniques to help a wide range of clinical researchers better understand and apply database technology.

A Recursive Approach to Quantized H∞ State Estimation for Genetic Regulatory Networks Under Stochastic Communication Protocols

This paper deals with the finite-horizon quantized H_∞ state estimation problem for a class of discrete time-varying genetic regulatory networks with quantization effects under stochastic communication protocols (SCPs). To better reflect the data-driven flavor of today's biological research, the network measurements (typically gigabytes in size by high-throughput sequencing technologies) are transmitted to a remote state estimator via two independent communication networks of limited bandwidths. To lighten the communication loads and avoid undesired data collisions, the measurement outputs are quantized and then transmitted under two SCPs introduced to schedule the large-scale data transmissions. The purpose of this paper is to design a time-varying state estimator such that the error dynamics of the state estimation satisfies a prescribed H_∞ performance requirement over a finite horizon in the presence of nonlinearities, quantization effects, and SCPs. By utilizing the completing-the-square technique, sufficient conditions are derived to ensure the H_∞ estimation performance and the parameters of the state estimator are designed by solving coupled backward recursive Riccati difference equations. A numerical example is given to illustrate the effectiveness of the design scheme of the proposed state estimator.

Intrinsic-overlapping co-expression module detection with application to Alzheimer's Disease

Genes interact with each other and may cause perturbation in the molecular pathways leading to complex diseases. Often, instead of any single gene, a subset of genes interact, forming a network, to share common biological functions. Such a subnetwork is called a functional module or motif. Identifying such modules and central key genes in them, that may be responsible for a disease, may help design patient-specific drugs. In this study, we consider the neurodegenerative Alzheimer's Disease (AD) and identify potentially responsible genes from functional motif analysis. We start from the hypothesis that central genes in genetic modules are more relevant to a disease that is under investigation and identify hub genes from the modules as potential marker genes. Motifs or modules are often non-exclusive or overlapping in nature. Moreover, they sometimes show intrinsic or hierarchical distributions with overlapping functional roles. To the best of our knowledge, no prior work handles both the situations in an integrated way. We propose a non-exclusive clustering approach, CluViaN (Clustering Via Network) that can detect intrinsic as well as overlapping modules from gene co-expression networks constructed using microarray expression profiles. We compare our method with existing methods to evaluate the quality of modules extracted. CluViaN reports the presence of intrinsic and overlapping motifs in different species not reported by any other research. We further apply our method to extract significant AD specific modules using CluViaN and rank them based the number of genes from a module involved in the disease pathways. Finally, top central genes are identified by topological analysis of the modules. We use two different AD phenotype data for experimentation. We observe that central genes, namely PSEN1, APP, NDUFB2, NDUFA1, UQCR10, PPP3R1 and a few more, play significant roles in the AD. Interestingly, our experiments also find a hub gene, PML, which has recently been reported to play a role in plasticity, circadian rhythms and the response to proteins which can cause neurodegenerative disorders. MUC4, another hub gene that we find experimentally is yet to be investigated for its potential role in AD. A software implementation of CluViaN in Java is available for download at https://sites.google.com/site/swarupnehu/publications/resources/CluViaN Software.rar.

Gene Coexpression Network Comparison via Persistent Homology

Persistent homology, a topological data analysis (TDA) method, is applied to microarray data sets. Although there are a few papers referring to TDA methods in microarray analysis, the usage of persistent homology in the comparison of several weighted gene coexpression networks (WGCN) was not employed before to the very best of our knowledge. We calculate the persistent homology of weighted networks constructed from 38 Arabidopsis microarray data sets to test the relevance and the success of this approach in distinguishing the stress factors. We quantify multiscale topological features of each network using persistent homology and apply a hierarchical clustering algorithm to the distance matrix whose entries are pairwise bottleneck distance between the networks. The immunoresponses to different stress factors are distinguishable by our method. The networks of similar immunoresponses are found to be close with respect to bottleneck distance indicating the similar topological features of WGCNs. This computationally efficient technique analyzing networks provides a quick test for advanced studies.

A comprehensive spatial-temporal transcriptomic analysis of differentiating nascent mouse lens epithelial and fiber cells

Elucidation of both the molecular composition and organization of the ocular lens is a prerequisite to understand its development, function, pathology, regenerative capacity, as well as to model lens development and disease using in vitro differentiation of pluripotent stem cells. Lens is comprised of the anterior lens epithelium and posterior lens fibers, which form the bulk of the lens. Lens fibers differentiate from lens epithelial cells through cell cycle exit-coupled differentiation that includes cellular elongation, accumulation of crystallins, cytoskeleton and membrane remodeling, and degradation of organelles within the central region of the lens. Here, we profiled spatiotemporal expression dynamics of both mRNAs and non-coding RNAs from microdissected mouse nascent lens epithelium and lens fibers at four developmental time points (embryonic [E] day 14.5, E16.5, E18.5, and P0.5) by RNA-seq. During this critical time window, multiple complex biosynthetic and catabolic processes generate the molecular and structural foundation for lens transparency. Throughout this developmental window, 3544 and 3518 genes show consistently and significantly greater expression in the nascent lens epithelium and fibers, respectively. Comprehensive data analysis confirmed major roles of FGF-MAPK, Wnt/β-catenin, PI3K/AKT, TGF-β, and BMP signaling pathways and revealed significant novel contributions of mTOR, EIF2, EIF4, and p70S6K signaling in lens formation. Unbiased motif analysis within promoter regions of these genes with consistent expression changes between epithelium and fiber cells revealed an enrichment for both established (e.g. E2Fs, Etv5, Hsf4, c-Maf, MafG, MafK, N-Myc, and Pax6) transcription factors and a number of novel regulators of lens formation, such as Arntl2, Dmrta2, Stat5a, Stat5b, and Tulp3. In conclusion, the present RNA-seq data serves as a comprehensive reference resource for deciphering molecular principles of normal mammalian lens differentiation, mapping a full spectrum of signaling pathways and DNA-binding transcription factors operating in both lens compartments, and predicting novel pathways required to establish lens transparency.

NetControl4BioMed: a pipeline for biomedical data acquisition and analysis of network controllability

BACKGROUND

Network controllability focuses on discovering combinations of external interventions that can drive a biological system to a desired configuration. In practice, this approach translates into finding a combined multi-drug therapy in order to induce a desired response from a cell; this can lead to developments of novel therapeutic approaches for systemic diseases like cancer.

RESULT

We develop a novel bioinformatics data analysis pipeline called NetControl4BioMed based on the concept of target structural control of linear networks. Our pipeline generates novel molecular interaction networks by combining pathway data from various public databases starting from the user's query. The pipeline then identifies a set of nodes that is enough to control a given, user-defined set of disease-specific essential proteins in the network, i.e., it is able to induce a change in their configuration from any initial state to any final state. We provide both the source code of the pipeline as well as an online web-service based on this pipeline http://combio.abo.fi/nc/net_control/remote_call.php .

CONCLUSION

The pipeline can be used by researchers for controlling and better understanding of molecular interaction networks through combinatorial multi-drug therapies, for more efficient therapeutic approaches and personalised medicine.

Gene Regulatory Network Inference from Perturbed Time-Series Expression Data via Ordered Dynamical Expansion of Non-Steady State Actors

The reconstruction of gene regulatory networks from gene expression data has been the subject of intense research activity. A variety of models and methods have been developed to address different aspects of this important problem. However, these techniques are narrowly focused on particular biological and experimental platforms, and require experimental data that are typically unavailable and difficult to ascertain. The more recent availability of higher-throughput sequencing platforms, combined with more precise modes of genetic perturbation, presents an opportunity to formulate more robust and comprehensive approaches to gene network inference. Here, we propose a step-wise framework for identifying gene-gene regulatory interactions that expand from a known point of genetic or chemical perturbation using time series gene expression data. This novel approach sequentially identifies non-steady state genes post-perturbation and incorporates them into a growing series of low-complexity optimization problems. The governing ordinary differential equations of this model are rooted in the biophysics of stochastic molecular events that underlie gene regulation, delineating roles for both protein and RNA-mediated gene regulation. We show the successful application of our core algorithms for network inference using simulated and real datasets.

State Estimation for Discrete Time-Delayed Genetic Regulatory Networks With Stochastic Noises Under the Round-Robin Protocols

This paper investigates the problem of state estimation for discrete time-delayed genetic regulatory networks with stochastic process noises and bounded exogenous disturbances under the Round-Robin (RR) protocols. The network measurement outputs obtained by two groups of sensors are transmitted to two remote sub-estimators via two independent communication channels, respectively. To lighten the communication loads of the networks and reduce the occurrence rate of data collisions, two RR protocols are utilized to orchestrate the transmission orders of sensor nodes in two groups, respectively. The error dynamics of the state estimation is governed by a switched system with periodic switching parameters. By constructing a transmission-order-dependent Lyapunov-like functional and utilizing the up-to-date discrete Wirtinger-based inequality together with the reciprocally convex approach, sufficient conditions are established to guarantee the exponentially ultimate boundedness of the estimation error dynamics in mean square with a prescribed upper bound on the decay rate. An asymptotic upper bound of the outputs of the estimation errors in mean square is derived and the estimator parameters are then obtained by minimizing such an upper bound subject to linear matrix inequality constraints. The repressilator model is utilized to illustrate the effectiveness of the designed estimator.

An efficient method to transcription factor binding sites imputation via simultaneous completion of multiple matrices with positional consistency

Transcription factors (TFs) are DNA-binding proteins that have a central role in regulating gene expression. Identification of DNA-binding sites of TFs is a key task in understanding transcriptional regulation, cellular processes and disease. Chromatin immunoprecipitation followed by high-throughput sequencing (ChIP-seq) enables genome-wide identification of in vivo TF binding sites. However, it is still difficult to map every TF in every cell line owing to cost and biological material availability, which poses an enormous obstacle for integrated analysis of gene regulation. To address this problem, we propose a novel computational approach, TFBSImpute, for predicting additional TF binding profiles by leveraging information from available ChIP-seq TF binding data. TFBSImpute fuses the dataset to a 3-mode tensor and imputes missing TF binding signals via simultaneous completion of multiple TF binding matrices with positional consistency. We show that signals predicted by our method achieve overall similarity with experimental data and that TFBSImpute significantly outperforms baseline approaches, by assessing the performance of imputation methods against observed ChIP-seq TF binding profiles. Besides, motif analysis shows that TFBSImpute preforms better in capturing binding motifs enriched in observed data compared with baselines, indicating that the higher performance of TFBSImpute is not simply due to averaging related samples. We anticipate that our approach will constitute a useful complement to experimental mapping of TF binding, which is beneficial for further study of regulation mechanisms and disease.

The Molecular Revolution in Cutaneous Biology: Era of Next-Generation Sequencing

Like any true conceptual revolution, next-generation sequencing (NGS) has not only radically changed research and clinical practice, it has also modified scientific culture. With the possibility to investigate DNA contents of any organism and in any context, including in somatic disorders or in tissues carrying complex microbial populations, it initially seemed as if the genetic underpinning of any biological phenomenon could now be deciphered in an almost streamlined fashion. However, over the past recent years, we have once again come to understand that there is no such a thing as great opportunities without great challenges. The steadily expanding use of NGS and related applications is now facing biologists and physicians with novel technological obstacles, analytical hurdles and increasingly pressing ethical questions.

Computational inference of gene regulatory networks: Approaches, limitations and opportunities

Gene regulatory networks lie at the core of cell function control. In E. coli and S. cerevisiae, the study of gene regulatory networks has led to the discovery of regulatory mechanisms responsible for the control of cell growth, differentiation and responses to environmental stimuli. In plants, computational rendering of gene regulatory networks is gaining momentum, thanks to the recent availability of high-quality genomes and transcriptomes and development of computational network inference approaches. Here, we review current techniques, challenges and trends in gene regulatory network inference and highlight challenges and opportunities for plant science. We provide plant-specific application examples to guide researchers in selecting methodologies that suit their particular research questions. Given the interdisciplinary nature of gene regulatory network inference, we tried to cater to both biologists and computer scientists to help them engage in a dialogue about concepts and caveats in network inference. Specifically, we discuss problems and opportunities in heterogeneous data integration for eukaryotic organisms and common caveats to be considered during network model evaluation. This article is part of a Special Issue entitled: Plant Gene Regulatory Mechanisms and Networks, edited by Dr. Erich Grotewold and Dr. Nathan Springer.

MetaNetVar: Pipeline for applying network analysis tools for genomic variants analysis

Network analysis can make variant analysis better. There are existing tools like HotNet2 and dmGWAS that can provide various analytical methods. We developed a prototype of a pipeline called MetaNetVar that allows execution of multiple tools. The code is published at https://github.com/NCBI-Hackathons/Network_SNPs. A working prototype is published as an Amazon Machine Image - ami-4510312f .

Towards an Enhanced Understanding of Plant-Microbiome Interactions to Improve Phytoremediation: Engineering the Metaorganism

Phytoremediation is a promising technology to clean-up contaminated soils based on the synergistic actions of plants and microorganisms. However, to become a widely accepted, and predictable remediation alternative, a deeper understanding of the plant-microbe interactions is needed. A number of studies link the success of phytoremediation to the plant-associated microbiome functioning, though whether the microbiome can exist in alternative, functional states for soil remediation, is incompletely understood. Moreover, current approaches that target the plant host, and environment separately to improve phytoremediation, potentially overlook microbial functions and properties that are part of the multiscale complexity of the plant-environment wherein biodegradation takes place. In contrast, in situ studies of phytoremediation research at the metaorganism level (host and microbiome together) are lacking. Here, we discuss a competition-driven model, based on recent evidence from the metagenomics level, and hypotheses generated by microbial community ecology, to explain the establishment of a catabolic rhizosphere microbiome in a contaminated soil. There is evidence to ground that if the host provides the right level and mix of resources (exudates) over which the microbes can compete, then a competitive catabolic and plant-growth promoting (PGP) microbiome can be selected for as long as it provides a competitive superiority in the niche. The competition-driven model indicates four strategies to interfere with the microbiome. Specifically, the rhizosphere microbiome community can be shifted using treatments that alter the host, resources, environment, and that take advantage of prioritization in inoculation. Our model and suggestions, considering the metaorganism in its natural context, would allow to gain further knowledge on the plant-microbial functions, and facilitate translation to more effective, and predictable phytotechnologies.

Reference sequence (RefSeq) database at NCBI: current status, taxonomic expansion, and functional annotation

The RefSeq project at the National Center for Biotechnology Information (NCBI) maintains and curates a publicly available database of annotated genomic, transcript, and protein sequence records (http://www.ncbi.nlm.nih.gov/refseq/). The RefSeq project leverages the data submitted to the International Nucleotide Sequence Database Collaboration (INSDC) against a combination of computation, manual curation, and collaboration to produce a standard set of stable, non-redundant reference sequences. The RefSeq project augments these reference sequences with current knowledge including publications, functional features and informative nomenclature. The database currently represents sequences from more than 55,000 organisms (>4800 viruses, >40,000 prokaryotes and >10,000 eukaryotes; RefSeq release 71), ranging from a single record to complete genomes. This paper summarizes the current status of the viral, prokaryotic, and eukaryotic branches of the RefSeq project, reports on improvements to data access and details efforts to further expand the taxonomic representation of the collection. We also highlight diverse functional curation initiatives that support multiple uses of RefSeq data including taxonomic validation, genome annotation, comparative genomics, and clinical testing. We summarize our approach to utilizing available RNA-Seq and other data types in our manual curation process for vertebrate, plant, and other species, and describe a new direction for prokaryotic genomes and protein name management.

Budding off: bringing functional genomics to Candida albicans

Candida species are the most prevalent human fungal pathogens, with Candida albicans being the most clinically relevant species. Candida albicans resides as a commensal of the human gastrointestinal tract but is a frequent cause of opportunistic mucosal and systemic infections. Investigation of C. albicans virulence has traditionally relied on candidate gene approaches, but recent advances in functional genomics have now facilitated global, unbiased studies of gene function. Such studies include comparative genomics (both between and within Candida species), analysis of total RNA expression, and regulation and delineation of protein-DNA interactions. Additionally, large collections of mutant strains have begun to aid systematic screening of clinically relevant phenotypes. Here, we will highlight the development of functional genomics in C. albicans and discuss the use of these approaches to addressing both commensalism and pathogenesis in this species.

Network measures for chemical library design

In this overview, we examine recent developments in network approaches to drug design. A brief overview of networks is followed by a discussion of how chemical similarity networks and their properties address challenges in drug design. Multiple methods used to assess or enhance chemical diversity for early-stage drug discovery are discussed, as well as methods that can be used for drug repositioning and ligand polypharmacology.

Transcriptional networks in plant immunity

Next to numerous abiotic stresses, plants are constantly exposed to a variety of pathogens within their environment. Thus, their ability to survive and prosper during the course of evolution was strongly dependent on adapting efficient strategies to perceive and to respond to such potential threats. It is therefore not surprising that modern plants have a highly sophisticated immune repertoire consisting of diverse signal perception and intracellular signaling pathways. This signaling network is intricate and deeply interconnected, probably reflecting the diverse lifestyles and infection strategies used by the multitude of invading phytopathogens. Moreover it allows signal communication between developmental and defense programs thereby ensuring that plant growth and fitness are not significantly retarded. How plants integrate and prioritize the incoming signals and how this information is transduced to enable appropriate immune responses is currently a major research area. An important finding has been that pathogen-triggered cellular responses involve massive transcriptional reprogramming within the host. Additional key observations emerging from such studies are that transcription factors (TFs) are often sites of signal convergence and that signal-regulated TFs act in concert with other context-specific TFs and transcriptional co-regulators to establish sensory transcription regulatory networks required for plant immunity.

Environmentally induced changes in correlated responses to selection reveal variable pleiotropy across a complex genetic network

Selection in novel environments can lead to a coordinated evolutionary response across a suite of characters. Environmental conditions can also potentially induce changes in the genetic architecture of complex traits, which in turn could alter the pattern of the multivariate response to selection. We describe a factorial selection experiment using the nematode Caenorhabditis remanei in which two different stress-related phenotypes (heat and oxidative stress resistance) were selected under three different environmental conditions. The pattern of covariation in the evolutionary response between phenotypes or across environments differed depending on the environment in which selection occurred, including asymmetrical responses to selection in some cases. These results indicate that variation in pleiotropy across the stress response network is highly sensitive to the external environment. Our findings highlight the complexity of the interaction between genes and environment that influences the ability of organisms to acclimate to novel environments. They also make clear the need to identify the underlying genetic basis of genetic correlations in order understand how patterns of pleiotropy are distributed across complex genetic networks.

Profiling the transcription factor regulatory networks of human cell types

Neph et al. (2012) (Circuitry and dynamics of human transcription factor regulatory networks. Cell, 150: 1274-1286) reported the transcription factor (TF) regulatory networks of 41 human cell types using the DNaseI footprinting technique. This provides a valuable resource for uncovering regulation principles in different human cells. In this paper, the architectures of the 41 regulatory networks and the distributions of housekeeping and specific regulatory interactions are investigated. The TF regulatory networks of different human cell types demonstrate similar global three-layer (top, core and bottom) hierarchical architectures, which are greatly different from the yeast TF regulatory network. However, they have distinguishable local organizations, as suggested by the fact that wiring patterns of only a few TFs are enough to distinguish cell identities. The TF regulatory network of human embryonic stem cells (hESCs) is dense and enriched with interactions that are unseen in the networks of other cell types. The examination of specific regulatory interactions suggests that specific interactions play important roles in hESCs.

Escape from the laboratory: new horizons for plant genetics

Next generation sequencing (NGS) is changing the way biologists work, as large amounts of genetic data can be easily outsourced commercially. Consequently, crucial research efforts in plant genetics can now be found outside the traditional laboratory setting, allowing for novel and more challenging scientific questions to be answered by virtual collaborative networks.

Network dynamics in the tumor microenvironment

The evolutionary path from tumor initiation to metastasis can only be fully understood by considering cancer cells as part of a multi-species ecosystem within the tumor microenvironment. This paper reviews and suggests two important recent trends. Firstly, I review arguments that interactions among diverse cells in the tumor microenvironment create a distinct cellular environment that can confer growth advantages, resist interventions, and allow tumors to remain dormant for long periods. Second, I review and highlight a trend toward data-rich, molecularly detailed, computational models of the tumor microenvironment. I argue that data-driven molecularly detailed tumor microenvironment models can now be built using data from multiple emerging high-throughput technologies, and that such models can pinpoint mechanisms of dysregulation and suggest specific drug targets and follow up experiments.