Ningtao Wang

Reconstructing regulatory networks from the dynamic plasticity of gene expression by mutual information

The capacity of an organism to respond to its environment is facilitated by the environmentally induced alteration of gene and protein expression, i.e. expression plasticity. The reconstruction of gene regulatory networks based on expression plasticity can gain not only new insights into the causality of transcriptional and cellular processes but also the complex regulatory mechanisms that underlie biological function and adaptation. We describe an approach for network inference by integrating expression plasticity into Shannon’s mutual information. Beyond Pearson correlation, mutual information can capture non-linear dependencies and topology sparseness. The approach measures the network of dependencies of genes expressed in different environments, allowing the environment-induced plasticity of gene dependencies to be tested in unprecedented details. The approach is also able to characterize the extent to which the same genes trigger different amounts of expression in response to environmental changes. We demonstrated the usefulness of this approach through analysing gene expression data from a rabbit vein graft study that includes two distinct blood flow environments. The proposed approach provides a powerful tool for the modelling and analysis of dynamic regulatory networks using gene expression data from distinct environments.

A Model Framework for Identifying Genes that Guide the Evolution of Heterochrony

Heterochrony, the phylogenic change in the time of developmental events or rate of development, has been thought to play an important role in producing phenotypic novelty during evolution. Increasing evidence suggests that specific genes are implicated in heterochrony, guiding the process of developmental divergence, but no quantitative models have been instrumented to map such heterochrony genes. Here, we present a computational framework for genetic mapping by which to characterize and locate quantitative trait loci (QTLs) that govern heterochrony described by four parameters, the timing of the inflection point, the timing of maximum acceleration of growth, the timing of maximum deceleration of growth, and the length of linear growth. The framework was developed from functional mapping, a dynamic model derived to map QTLs for the overall process and pattern of development. By integrating an optimality algorithm, the framework allows the so-called heterochrony QTLs (hQTLs) to be tested and quantified. Specific pipelines are given for testing how hQTLs control the onset and offset of developmental events, the rate of development, and duration of a particular developmental stage. Computer simulation was performed to examine the statistical properties of the model and demonstrate its utility to characterize the effect of hQTLs on population diversification due to heterochrony. By analyzing a genetic mapping data in rice, the framework identified an hQTL that controls the timing of maximum growth rate and duration of linear growth stage in plant height growth. The framework provides a tool to study how genetic variation translates into phenotypic innovation, leading a lineage to evolve, through heterochrony.

Towards a comprehensive picture of the genetic landscape of complex traits

The formation of phenotypic traits, such as biomass production, tumor volume and viral abundance, undergoes a complex process in which interactions between genes and developmental stimuli take place at each level of biological organization from cells to organisms. Traditional studies emphasize the impact of genes by directly linking DNA-based markers with static phenotypic values. Functional mapping, derived to detect genes that control developmental processes using growth equations, has proven powerful for addressing questions about the roles of genes in development. By treating phenotypic formation as a cohesive system using differential equations, a different approach-systems mapping-dissects the system into interconnected elements and then map genes that determine a web of interactions among these elements, facilitating our understanding of the genetic machineries for phenotypic development. Here, we argue that genetic mapping can play a more important role in studying the genotype-phenotype relationship by filling the gaps in the biochemical and regulatory process from DNA to end-point phenotype. We describe a new framework, named network mapping, to study the genetic architecture of complex traits by integrating the regulatory networks that cause a high-order phenotype. Network mapping makes use of a system of differential equations to quantify the rule by which transcriptional, proteomic and metabolomic components interact with each other to organize into a functional whole. The synthesis of functional mapping, systems mapping and network mapping provides a novel avenue to decipher a comprehensive picture of the genetic landscape of complex phenotypes that underlie economically and biomedically important traits.

A multivalent three-point linkage analysis model of autotetraploids

Because of its widespread occurrence and role in shaping evolutionary processes in the biological kingdom, especially in plants, polyploidy has been increasingly studied from cytological to molecular levels. By inferring gene order, gene distances and gene homology, linkage mapping with molecular markers has proven powerful for investigating genome structure and organization. Here we review and assess a general statistical model for three-point linkage analysis in autotetraploids by integrating double reduction, a phenomenon that commonly occurs in autopolyploids whose chromosomes are derived from a single ancestral species. This model does not require any assumption on the distribution of the occurrence of double reduction and can handle the complexity of multilocus linkage in terms of crossover interference. Implemented with the expectation-maximization (EM) algorithms, the model can estimate and test the recombination fractions between less informative dominant markers, thus facilitating its practical implications for any autopolyploids in most of which inexpensive dominant markers are still used for their genetic and evolutionary studies. The model was applied to reanalyze a published data in tetraploid switchgrass, validating its practical usefulness and utilization.

Delivering systems pharmacogenomics towards precision medicine through mathematics

The latest developments of pharmacology in the post-genomic era foster the emergence of new biomarkers that represent the future of drug targets. To identify these biomarkers, we need a major shift from traditional genomic analyses alone, moving the focus towards systems approaches to elucidating genetic variation in biochemical pathways of drug response. Is there any general model that can accelerate this shift via a merger of systems biology and pharmacogenomics? Here we describe a statistical framework for mapping dynamic genes that affect drug response by incorporating its pharmacokinetic and pharmacodynamic pathways. This framework is expanded to shed light on the mechanistic and therapeutic differences of drug response based on pharmacogenetic information, coupled with genomic, proteomic and metabolic data, allowing novel therapeutic targets and genetic biomarkers to be characterized and utilized for drug discovery.

A quantitative model of transcriptional differentiation driving host–pathogen interactions

Despite our expanding knowledge about the biochemistry of gene regulation involved in host-pathogen interactions, a quantitative understanding of this process at a transcriptional level is still limited. We devise and assess a computational framework that can address this question. This framework is founded on a mixture model-based likelihood, equipped with functionality to cluster genes per dynamic and functional changes of gene expression within an interconnected system composed of the host and pathogen. If genes from the host and pathogen are clustered in the same group due to a similar pattern of dynamic profiles, they are likely to be reciprocally co-evolving. If genes from the two organisms are clustered in different groups, this means that they experience strong host-pathogen interactions. The framework can test the rates of change for individual gene clusters during pathogenic infection and quantify their impacts on host-pathogen interactions. The framework was validated by a pathological study of poplar leaves infected by fungal Marssonina brunnea in which co-evolving and interactive genes that determine poplar-fungus interactions are identified. The new framework should find its wide application to studying host-pathogen interactions for any other interconnected systems.

Dynamic modeling of genes controlling cancer stem cell proliferation

The growing evidence that cancer originates from stem cells (SC) holds a great promise to eliminate this disease by designing specific drug therapies for removing cancer SC. Translation of this knowledge into predictive tests for the clinic is hampered due to the lack of methods to discriminate cancer SC from non-cancer SC. Here, we address this issue by describing a conceptual strategy for identifying the genetic origins of cancer SC. The strategy incorporates a high-dimensional group of differential equations that characterizes the proliferation, differentiation, and reprogramming of cancer SC in a dynamic cellular and molecular system. The deployment of robust mathematical models will help uncover and explain many still unknown aspects of cell behavior, tissue function, and network organization related to the formation and division of cancer SC. The statistical method developed allows biologically meaningful hypotheses about the genetic control mechanisms of carcinogenesis and metastasis to be tested in a quantitative manner.

Systems mapping of genes controlling chemotherapeutic drug efficiency for cancer stem cells

Cancer can be controlled effectively by using chemotherapeutic drugs to inhibit cancer stem cells, but there is considerable inter-patient variability regarding how these cells respond to drug intervention. Here, we describe a statistical framework for mapping genes that control tumor responses to chemotherapeutic drugs as well as the efficacy of treatments in arresting tumor growth. The framework integrates the mathematical aspects of the cancer stem cell hypothesis into genetic association studies, equipped with a capacity to quantify the magnitude and pattern of genetic effects on the kinetic decline of cancer stem cells in response to therapy. By quantifying how specific genes and their interactions govern drug response, the model provides essential information to tailor personalized drugs for individual patients.

A block mixture model to map eQTLs for gene clustering and networking

To study how genes function in a cellular and physiological process, a general procedure is to classify gene expression profiles into categories based on their similarity and reconstruct a regulatory network for functional elements. However, this procedure has not been implemented with the genetic mechanisms that underlie the organization of gene clusters and networks, despite much effort made to map expression quantitative trait loci (eQTLs) that affect the expression of individual genes. Here we address this issue by developing a computational approach that integrates gene clustering and network reconstruction with genetic mapping into a unifying framework. The approach can not only identify specific eQTLs that control how genes are clustered and organized toward biological functions, but also enable the investigation of the biological mechanisms that individual eQTLs perturb in a signaling pathway. We applied the new approach to characterize the effects of eQTLs on the structure and organization of gene clusters in Caenorhabditis elegans. This study provides the first characterization, to our knowledge, of the effects of genetic variants on the regulatory network of gene expression. The approach developed can also facilitate the genetic dissection of other dynamic processes, including development, physiology and disease progression in any organisms.

A Computational Algorithm for Functional Clustering of Proteome Dynamics During Development

Phenotypic traits, such as seed development, are a consequence of complex biochemical interactions among genes, proteins and metabolites, but the underlying mechanisms that operate in a coordinated and sequential manner remain elusive. Here, we address this issue by developing a computational algorithm to monitor proteome changes during the course of trait development. The algorithm is built within the mixture-model framework in which each mixture component is modeled by a specific group of proteins that display a similar temporal pattern of expression in trait development. A nonparametric approach based on Legendre orthogonal polynomials was used to fit dynamic changes of protein expression, increasing the power and flexibility of protein clustering. By analyzing a dataset of proteomic dynamics during early embryogenesis of the Chinese fir, the algorithm has successfully identified several distinct types of proteins that coordinate with each other to determine seed development in this forest tree commercially and environmentally important to China. The algorithm will find its immediate applications for the characterization of mechanistic underpinnings for any other biological processes in which protein abundance plays a key role.