Deciphering the Palimpsest: Studying the Relationship Between Morphological Integration and Phenotypic Covariation

Organisms represent a complex arrangement of anatomical structures and individuated parts that must maintain functional associations through development. This integration of variation between functionally related body parts and the modular organization of development are fundamental determinants of their evolvability. This is because integration results in the expression of coordinated variation that can create preferred directions for evolutionary change, while modularity enables variation in a group of traits or regions to accumulate without deleterious effects on other aspects of the organism. Using our own work on both model systems (e.g., lab mice, avians) and natural populations of rodents and primates, we explore in this paper the relationship between patterns of phenotypic covariation and the developmental determinants of integration that those patterns are assumed to reflect. We show that integration cannot be reliably studied through phenotypic covariance patterns alone and argue that the relationship between phenotypic covariation and integration is obscured in two ways. One is the superimposition of multiple determinants of covariance in complex systems and the other is the dependence of covariation structure on variances in covariance-generating processes. As a consequence, we argue that the direct study of the developmental determinants of integration in model systems is necessary to fully interpret patterns of covariation in natural populations, to link covariation patterns to the processes that generate them, and to understand their significance for evolutionary explanation.     

Allometric covariation: a hallmark behavior of plants and leaves

Size is one of the most important axes of variation among plants. As such, plant biologists have long searched for unifying principles that can explain how matter and energy flux and organ partitioning scale with plant size. Several recent models have proposed a universal biophysical basis for numerous scaling phenomena in plants based on vascular network geometry. Here, we review statistical analyses of several large-scale plant datasets that demonstrate that a true hallmark of plant form variability is systematic covariation among traits. This covariation is constrained by allometries that combine and trade off with one another, rather than any single universal allometric scaling exponent for a trait or suite of traits. Further, we show that covariation can be successfully modeled using network approaches that allow for species-specific designs in plants and geometric approaches that constrain relationships among economic traits in leaves. Finally, we report large-scale efforts utilizing semi-automated software tools that quantify physical networks and can inform our attempts to link vascular network structure to plant form and function. Collectively, this work highlights how the linking of morphology, biomass partitioning and the structure of physical distribution networks can improve our empirical and theoretical understanding of important drivers of plant functional diversity.     

Morphological variation and covariation in mandibular molars of platyrrhine primates

Molars are highly integrated biological structures that have been used for inferring evolutionary relationships among taxa. However, parallel and convergent morphological traits can be affected by developmental and functional constraints. Here, we analyze molar shapes of platyrrhines in order to explore if platyrrhine molar diversity reflects homogeneous patterns of molar variation and covariation. We digitized 30 landmarks on mandibular first and second molars of 418 extant and 11 fossil platyrrhine specimens to determine the degree of integration of both molars when treated as a single module. We combined morphological and phylogenetic data to investigate the phylogenetic signal and to visualize the history of molar shape changes. All platyrrhine taxa show a common shape pattern suggesting that a relatively low degree of phenotypic variation is caused by convergent evolution, although molar shape carries significant phylogenetic signal. Atelidae and Pitheciidae show high levels of integration with low variation between the two molars, whereas the Cebinae/Saimiriinae, and especially Callitrichinae, show greater variation between molars and trend toward a modular organization. We hypothesize that biomechanical constraints of the masticatory apparatus, and the dietary profile of each taxon are the main factors that determine high covariation in molars. In contrast, low molar shape covariation may result from the fact that each molar exhibits a distinct ecological signal, as molars can be exposed to distinct occlusal loadings during food processing, suggesting that different selective pressures on molars can reduce overall molar integration.     

Biological network design strategies: discovery through dynamic optimization

An important challenge in systems biology is the inherent complexity of biological network models, which complicates the task of relating network structure to function and of understanding the conceptual design principles by which a given network operates. Here we investigate an approach to analyze the relationship between a network structure and its function using the framework of optimization. A common feature found in a variety of biochemical networks involves the opposition of a pair of enzymatic chemical modification reactions such as phosphorylation-dephosphorylation or methylation-demethylation. The modification pair frequently adjusts biochemical properties of its target, such as activating and deactivating function. We applied optimization methodology to study a reversible modification network unit commonly found in signal transduction systems, and we explored the use of this methodology to discover design principles. The results demonstrate that different sets of rate constants used to parameterize the same network topology represent different compromises made in the resulting network operating characteristics. Moreover, the same topology can be used to encode different strategies for achieving performance goals. The ability to adopt multiple strategies may lead to significantly improved performance across a range of conditions through rate modulation or evolutionary processes. The optimization framework explored here is a practical approach to support the discovery of design principles in biological networks.     

Module-Based Association Analysis for Omics Data with Network Structure

Module-based analysis (MBA) aims to evaluate the effect of a group of biological elements sharing common features, such as SNPs in the same gene or metabolites in the same pathways, and has become an attractive alternative to traditional single bio-element approaches. Because bio-elements regulate and interact with each other as part of network, incorporating network structure information can more precisely model the biological effects, enhance the ability to detect true associations, and facilitate our understanding of the underlying biological mechanisms. How-ever, most MBA methods ignore the network structure information, which depicts the interaction and regulation relationship among basic functional units in biology system. We construct the con-nectivity kernel and the topology kernel to capture the relationship among bio-elements in a mod-ule, and use a kernel machine framework to evaluate the joint effect of bio-elements. Our proposed kernel machine approach directly incorporates network structure so to enhance the study effi-ciency; it can assess interactions among modules, account covariates, and is computational effi-cient. Through simulation studies and real data application, we demonstrate that the proposed network-based methods can have markedly better power than the approaches ignoring network information under a range of scenarios.     

DISSECTING CORRELATED CHARACTERS: ADAPTIVE ASPECTS OF PHENOTYPIC COVARIATION IN MELANIZATION PATTERN OF PIERIS BUTTERFLIES

In this study we address the question of how much of the covariation among phenotypic characters observed in natural populations is adaptive. We examine covariation among a set of phenotypic characters that describe the wing-melanization pattern of Pieris butterflies. Previous functional analyses of thermoregulatory performance allow us to predict a priori whether and how different wing melanic characters should be correlated. We quantify and analyze the variation in the wing-melanization pattern within species for a series of Pieris populations from relatively cool environments in North America and compare these results with the predictions based on our adaptive hypothesis. We consider adaptive covariation both for biogeographic variation among populations and for seasonal polyphenism (phenotypic plasticity) within populations. Our hypothesis correctly predicts many of the qualitative features of covariation in melanization among major regions of the wings, at the level of biogeographic variation among populations, for both males and females of Pieris occidentalis. When within-population variation is considered, agreement with the adaptive predictions varies considerably in different populations for both P. occidentalis and P. napi males and females. Agreement for P. napi, particularly the females, is generally poorer than for P. occidentalis. In both species, there is a consistent difference in melanization pattern between alpine and arctic sites; this difference is discussed in relation to the differences in the radiative environment between these two types of “cold” habitats. Our results suggest that some important aspects of phenotypic correlation among wing melanic characters in Pieris are adaptive. We emphasize the important distinction between covariation and co-occurrence of characters, and we discuss these results in relation to the extensive biogeographic variation and phenotypic plasticity (seasonal polyphenism) in Pieris wing-melanization patterns.     

Mistranslation-induced protein misfolding as a dominant constraint on coding-sequence evolution

Strikingly consistent correlations between rates of coding-sequence evolution and gene expression levels are apparent across taxa, but the biological causes behind the selective pressures on coding-sequence evolution remain controversial. Here, we demonstrate conserved patterns of simple covariation between sequence evolution, codon usage, and mRNA level in E. coli, yeast, worm, fly, mouse, and human that suggest that all observed trends stem largely from a unified underlying selective pressure. In metazoans, these trends are strongest in tissues composed of neurons, whose structure and lifetime confer extreme sensitivity to protein misfolding. We propose, and demonstrate using a molecular-level evolutionary simulation, that selection against toxicity of misfolded proteins generated by ribosome errors suffices to create all of the observed covariation. The mechanistic model of molecular evolution that emerges yields testable biochemical predictions, calls into question the use of nonsynonymous-to-synonymous substitution ratios (Ka/Ks) to detect functional selection, and suggests how mistranslation may contribute to neurodegenerative disease.     

Genetics of body shape and armour variation in threespine sticklebacks

Patterns of genetic variation and covariation can influence the rate and direction of phenotypic evolution. We explored the possibility that the parallel morphological evolution seen in threespine stickleback (Gasterosteus aculeatus) populations colonizing freshwater environments is facilitated by patterns of genetic variation and covariation in the ancestral (marine) population. We estimated the genetic (G) and phenotypic (P) covariance matrices and directions of maximum additive genetic (g(max) ) and phenotypic (p(max) ) covariances of body shape and armour traits. Our results suggest a role for the ancestral G in explaining parallel morphological evolution in freshwater populations. We also found evidence of genetic constraints owing to the lack of variance in the ancestral G. Furthermore, strong genetic covariances and correlations among traits revealed that selective factors responsible for threespine stickleback body shape and armour divergence may be difficult to disentangle. The directions of g(max) and p(max) were correlated, but the correlations were not high enough to imply that phenotypic patterns of trait variation and covariation within populations are very informative of underlying genetic patterns.     

A METHOD FOR ASSESSING PHYLOGENETIC LEAST SQUARES MODELS FOR SHAPE AND OTHER HIGH‐DIMENSIONAL MULTIVARIATE DATA

Studies of evolutionary correlations commonly use phylogenetic regression (i.e., independent contrasts and phylogenetic generalized least squares) to assess trait covariation in a phylogenetic context. However, while this approach is appropriate for evaluating trends in one or a few traits, it is incapable of assessing patterns in highly multivariate data, as the large number of variables relative to sample size prohibits parametric test statistics from being computed. This poses serious limitations for comparative biologists, who must either simplify how they quantify phenotypic traits, or alter the biological hypotheses they wish to examine. In this article, I propose a new statistical procedure for performing ANOVA and regression models in a phylogenetic context that can accommodate high‐dimensional datasets. The approach is derived from the statistical equivalency between parametric methods using covariance matrices and methods based on distance matrices. Using simulations under Brownian motion, I show that the method displays appropriate Type I error rates and statistical power, whereas standard parametric procedures have decreasing power as data dimensionality increases. As such, the new procedure provides a useful means of assessing trait covariation across a set of taxa related by a phylogeny, enabling macroevolutionary biologists to test hypotheses of adaptation, and phenotypic change in high‐dimensional datasets.     

Understanding biological functions through molecular networks

The completion of genome sequences and subsequent high-throughput mapping of molecular networks have allowed us to study biology from the network perspective. Experimental, statistical and mathematical modeling approaches have been employed to study the structure, function and dynamics of molecular networks, and begin to reveal important links of various network properties to the functions of the biological systems. In agreement with these functional links, evolutionary selection of a network is apparently based on the function, rather than directly on the structure of the network. Dynamic modularity is one of the prominent features of molecular networks. Taking advantage of such a feature may simplify network-based biological studies through construction of process-specific modular networks and provide functional and mechanistic insights linking genotypic variations to complex traits or diseases, which is likely to be a key approach in the next wave of understanding complex human diseases. With the development of ready-to-use network analysis and modeling tools the networks approaches will be infused into everyday biological research in the near future.     

Competing to coordinate cell fate decisions: the MST2-Raf-1 signaling device

How do biochemical signaling pathways generate biological specificity? This question is fundamental to modern biology, and its enigma has been accentuated by the discovery that most proteins in signaling networks serve multifunctional roles. An answer to this question may lie in analyzing network properties rather than individual traits of proteins in order to elucidate design principles of biochemical networks that enable biological decision-making. We discuss how this is achieved in the MST2/Hippo-Raf-1 signaling network with the help of mathematical modeling and model-based analysis, which showed that competing protein interactions with affinities controlled by dynamic protein modifications can function as Boolean computing devices that determine cell fate decisions. In addition, we discuss areas of interest for future research and highlight how systems approaches would be of benefit.     

Fine-Scale Dissection of Functional Protein Network Organization by Statistical Network Analysis

Revealing organizational principles of biological networks is an important goal of systems biology. In this study, we sought to analyze the dynamic organizational principles within the protein interaction network by studying the characteristics of individual neighborhoods of proteins within the network based on their gene expression as well as protein-protein interaction patterns. By clustering proteins into distinct groups based on their neighborhood gene expression characteristics, we identify several significant trends in the dynamic organization of the protein interaction network. We show that proteins with distinct neighborhood gene expression characteristics are positioned in specific localities in the protein interaction network thereby playing specific roles in the dynamic network connectivity. Remarkably, our analysis reveals a neighborhood characteristic that corresponds to the most centrally located group of proteins within the network. Further, we show that the connectivity pattern displayed by this group is consistent with the notion of “rich club connectivity” in complex networks. Importantly, our findings are largely reproducible in networks constructed using independent and different datasets.     

Organization of Enzyme Concentration across the Metabolic Network in Cancer Cells

Rapid advances in mass spectrometry have allowed for estimates of absolute concentrations across entire proteomes, permitting the interrogation of many important biological questions. Here, we focus on a quantitative aspect of human cancer cell metabolism that has been limited by a paucity of available data on the abundance of metabolic enzymes. We integrate data from recent measurements of absolute protein concentration to analyze the statistics of protein abundance across the human metabolic network. At a global level, we find that the enzymes in glycolysis comprise approximately half of the total amount of metabolic proteins and can constitute up to 10% of the entire proteome. We then use this analysis to investigate several outstanding problems in cancer metabolism, including the diversion of glycolytic flux for biosynthesis, the relative contribution of nitrogen assimilating pathways, and the origin of cellular redox potential. We find many consistencies with current models, identify several inconsistencies, and find generalities that extend beyond current understanding. Together our results demonstrate that a relatively simple analysis of the abundance of metabolic enzymes was able to reveal many insights into the organization of the human cancer cell metabolic network.     

Toward scalable parts families for predictable design of biological circuits

Our current ability to engineer biological circuits is hindered by design cycles that are costly in terms of time and money, with constructs failing to operate as desired, or evolving away from the desired function once deployed. Synthetic biologists seek to understand biological design principles and use them to create technologies that increase the efficiency of the genetic engineering design cycle. Central to the approach is the creation of biological parts--encapsulated functions that can be composited together to create new pathways with predictable behaviors. We define five desirable characteristics of biological parts--independence, reliability, tunability, orthogonality and composability, and review studies of small natural and synthetic biological circuits that provide insights into each of these characteristics. We propose that the creation of appropriate sets of families of parts with these properties is a prerequisite for efficient, predictable engineering of new function in cells and will enable a large increase in the sophistication of genetic engineering applications.     

Discordant patterns of genetic and phenotypic differentiation in five grasshopper species codistributed across a microreserve network

Conservation plans can be greatly improved when information on the evolutionary and demographic consequences of habitat fragmentation is available for several codistributed species. Here, we study spatial patterns of phenotypic and genetic variation among five grasshopper species that are codistributed across a network of microreserves but show remarkable differences in dispersal‐related morphology (body size and wing length), degree of habitat specialization and extent of fragmentation of their respective habitats in the study region. In particular, we tested the hypothesis that species with preferences for highly fragmented microhabitats show stronger genetic and phenotypic structure than codistributed generalist taxa inhabiting a continuous matrix of suitable habitat. We also hypothesized a higher resemblance of spatial patterns of genetic and phenotypic variability among species that have experienced a higher degree of habitat fragmentation due to their more similar responses to the parallel large‐scale destruction of their natural habitats. In partial agreement with our first hypothesis, we found that genetic structure, but not phenotypic differentiation, was higher in species linked to highly fragmented habitats. We did not find support for congruent patterns of phenotypic and genetic variability among any studied species, indicating that they show idiosyncratic evolutionary trajectories and distinctive demographic responses to habitat fragmentation across a common landscape. This suggests that conservation practices in networks of protected areas require detailed ecological and evolutionary information on target species to focus management efforts on those taxa that are more sensitive to the effects of habitat fragmentation.     

New insights into the physiology of natural foraging

The purpose of this symposium was to examine how foraging physiology is studied in the field across a diversity of species and habitats. While field studies are constrained by the relatively poor ability to control the experiment, the natural variability in both the environment and animal behavior provides insights into adaptation to change that are usually not tested in the laboratory. Talks in this session examined how foraging energy (both costs and gains) is partitioned over time. "Time," in this case, ranged from evolutionary time (how different animals are designed to most efficiently forage), to long, lifetime periods (development of foraging ability and growth), to short-duration feeding bouts, and ultimately to the minutes to hours following ingestion (metabolic and biochemical changes). From this diversity, two core themes emerged: that foraging strategies and behaviors are limited by physiology and biochemical processes and that time plays a central role in the organization of foraging behaviors and the physiological processes that underlie those behaviors.     

Codependency and mutual exclusivity for gene community detection from sparse single-cell transcriptome data

Single-cell RNA-seq (scRNA-seq) can be used to characterize cellular heterogeneity in thousands of cells. The reconstruction of a gene network based on coexpression patterns is a fundamental task in scRNA-seq analyses, and the mutual exclusivity of gene expression can be critical for understanding such heterogeneity. Here, we propose an approach for detecting communities from a genetic network constructed on the basis of coexpression properties. The community-based comparison of multiple coexpression networks enables the identification of functionally related gene clusters that cannot be fully captured through differential gene expression-based analysis. We also developed a novel metric referred to as the exclusively expressed index (EEI) that identifies mutually exclusive gene pairs from sparse scRNA-seq data. EEI quantifies and ranks the exclusive expression levels of all gene pairs from binary expression patterns while maintaining robustness against a low sequencing depth. We applied our methods to glioblastoma scRNA-seq data and found that gene communities were partially conserved after serum stimulation despite a considerable number of differentially expressed genes. We also demonstrate that the identification of mutually exclusive gene sets with EEI can improve the sensitivity of capturing cellular heterogeneity. Our methods complement existing approaches and provide new biological insights, even for a large, sparse dataset, in the single-cell analysis field.     

Subtle alterations in PCNA-partner interactions severely impair DNA replication and repair

The robustness of complex biological processes in the face of environmental and genetic perturbations is a key biological trait. However, while robustness has been extensively studied, little is known regarding the fragility of biological processes. Here, we have examined the susceptibility of DNA replication and repair processes mediated by the proliferating cell nuclear antigen (PCNA). Using protein directed evolution, biochemical, and genetic approaches, we have generated and characterized PCNA mutants with increased affinity for several key partners of the PCNA-partner network. We found that increases in PCNA-partner interaction affinities led to severe in vivo phenotypic defects. Surprisingly, such defects are much more severe than those induced by complete abolishment of the respective interactions. Thus, the subtle and tunable nature of these affinity perturbations produced different phenotypic effects than realized with traditional "on-off" analysis using gene knockouts. Our findings indicate that biological systems can be robust to one set of perturbations yet fragile to others.