Aposematic coloration offers an opportunity to explore the molecular mechanisms underlying canalization. In this study, the role of epigenetic regulation underlying robustness was explored in the aposematic coloration of the milkweed bug, Oncopeltus fasciatus. Polycomb (Pc) and Enhancer of zeste (E(z)), which encode components of the Polycomb repressive complex 1 (PRC1) and PRC2, respectively, and jing, which encodes a component of the PRC2.2 subcomplex, were knocked down in the fourth instar of O. fasciatus. Knockdown of these genes led to alterations in scutellar morphology and melanization. In particular, when Pc was knocked down, the adults developed a highly melanized abdomen, head and forewings at all temperatures examined. In contrast, the E(z) and jing knockdown led to increased plasticity of the dorsal forewing melanization across different temperatures. Moreover, jing knockdown adults exhibited increased plasticity in the dorsal melanization of the head and the thorax. These observations demonstrate that histone modifiers may play a key role during the process of canalization to confer robustness in the aposematic coloration.

Variability in quantitative traits has clinical, ecological, and evolutionary significance. Most genetic variants identified for complex quantitative traits have only a detectable effect on the mean of trait. We have developed the mean-variance test (MVtest) to simultaneously model the mean and log-variance of a quantitative trait as functions of genotypes and covariates by using estimating equations. The advantages of MVtest include the facts that it can detect effect modification, that multiple testing can follow conventional thresholds, that it is robust to non-normal outcomes, and that association statistics can be meta-analyzed. In simulations, we show control of type I error of MVtest over several alternatives. We identified 51 and 37 previously unreported associations for effects on blood-pressure variance and mean, respectively, in the UK Biobank. Transcriptome-wide association studies revealed 633 significant unique gene associations with blood-pressure mean variance. MVtest is broadly applicable to studies of complex quantitative traits and provides an important opportunity to detect novel loci.

Receptor-mediated signaling plays a central role in tissue regeneration, and it is dysregulated in disease. Here, we build a signaling-response map for a model regenerative human tissue: the airway epithelium. We analyzed the effect of 17 receptor-mediated signaling pathways on organotypic cultures to determine changes in abundance and phenotype of epithelial cell types. This map recapitulates the gamut of known airway epithelial signaling responses to these pathways. It defines convergent states induced by multiple ligands and diverse, ligand-specific responses in basal cell and secretory cell metaplasia. We show that loss of canonical differentiation induced by multiple pathways is associated with cell-cycle arrest, but that arrest is not sufficient to block differentiation. Using the signaling-response map, we show that a TGFB1-mediated response underlies specific aberrant cells found in multiple lung diseases and identify interferon responses in COVID-19 patient samples. Thus, we offer a framework enabling systematic evaluation of tissue signaling responses. A record of this paper\'s transparent peer review process is included in the supplemental information.

The concept of control is central to understanding and applications of biological network models. Some of their key structural features relate to control functions, through gene regulation, signaling, or metabolic mechanisms, and computational models need to encode these. Applications of models often focus on model-based control, such as in biomedicine or metabolic engineering. This paper presents an approach to model-based control that exploits two common features of biological networks, namely their modular structure and canalizing features of their regulatory mechanisms. The paper focuses on intracellular regulatory networks, represented by Boolean network models. A main result of this paper is that control strategies can be identified by focusing on one module at a time. This paper also presents a criterion based on canalizing features of the regulatory rules to identify modules that do not contribute to network control and can be excluded. For even moderately sized networks, finding global control inputs is computationally very challenging. The modular approach presented here leads to a highly efficient approach to solving this problem. This approach is applied to a published Boolean network model of blood cancer large granular lymphocyte (T-LGL) leukemia to identify a minimal control set that achieves a desired control objective.

Dysregulation of fatty acid metabolites can play a crucial role in the progression of complex diseases, such as cardiovascular disease, digestive diseases, and metabolic diseases. Metabolites can have either protective or risk effects on a disease; however, the details of such associations remain contentious. In this study, we demonstrate an integrative PheWAS approach to establish high confidence, causally suggestive of metabolite-disease associations for three fatty acid metabolites, namely, omega-3 fatty acids, omega-6 fatty acids, and docosahexaenoic acid, for 1,254 disease endpoints. Metabolite-disease associations were established if there was a concordant direction of effect and significance for metabolite level and genetic risk score for the metabolite. There was enrichment for metabolite associations with diseases of the respiratory system for omega-3 fatty acids, diseases of the circulatory system and endocrine system for omega-6 fatty acids, and diseases of the digestive system for docosahexaenoic acid. Upon performing Mendelian randomization on a subset of the outcomes, we identified 3, 6, and 15 significant diseases associated with omega-3 fatty acids, omega-6 fatty acids, and docosahexaenoic acid, respectively. We then demonstrate a class of prevalence-risk relationships indicative of (de)canalization of disease under high and low fatty acid metabolite levels. Finally, we show that the interaction between the metabolites and obesity demonstrates that the degree of protection afforded by fatty acid metabolites is strongly modulated by underlying metabolic health. This study evaluated the disease architectures of three polyunsaturated fatty acids (PUFAs), which were validated by several PheWAS modes of support. Our results not only highlight specific diseases associated with each metabolite but also disease group enrichments. In addition, we demonstrate an integrative PheWAS methodology that can be applied to other components of the human metabolome or other traits of interest. The results of this study can be used as an atlas to cross-compare genetic with non-genetic disease associations for the three PUFAs investigated. The findings can be explored through our R shiny app at https://pufa.biosci.gatech.edu.

Animal personality has been shown to be influenced by both genetic and environmental factors and shaped by natural selection. Currently, little is known about mechanisms influencing the development of personality traits. This study examines the extent to which personality development is genetically influenced and/or environmentally responsive (plastic). We also investigated the role of evolutionary history, assessing whether personality traits could be canalized along a genetic and ecological divergence gradient. We tested the plastic potential of boldness in juveniles of five Icelandic Arctic charr morphs (Salvelinus alpinus), including two pairs of sympatric morphs, displaying various degrees of genetic and ecological divergence from the ancestral anadromous charr, split between treatments mimicking benthic versus pelagic feeding modalities. We show that differences in mean boldness are mostly affected by genetics. While the benthic treatment led to bolder individuals overall, the environmental effect was rather weak, suggesting that boldness lies under strong genetic influence with reduced plastic potential. Finally, we found hints of differences by morphs in boldness canalization through reduced variance and plasticity, and higher consistency in boldness within morphs. These findings provide new insights on how behavioural development may impact adaptive diversification.

Morphological analyses are critical to quantify phenotypic variation, identify taxa, inform phylogenetic relationships, and shed light on evolutionary patterns. This work is particularly important in groups that display great morphological disparity. Such is the case in geomyoid rodents, a group that includes 2 of the most species-rich families of rodents in North America: the Geomyidae (pocket gophers) and the Heteromyidae (kangaroo rats, pocket mice, and their relatives). We assessed variation in skull morphology (including both shape and size) among geomyoids to test the hypothesis that there are statistically significant differences in skull measurements at the family, genus, and species levels. Our sample includes 886 specimens representing all geomyoid genera and 39 species. We used the geometric mean to compare size across taxa. We used 14 measurements of the cranium and lower jaw normalized for size to compare shape among and within taxa. Our results show that skull measurements enable the distinction of geomyoids at the family, genus, and species levels. There is a larger amount of size variation within Geomyidae than within Heteromyidae. Our phylomorphospace analysis shows that the skull shape of the common ancestor of all geomyoids was more similar to the common ancestor of heteromyids than that of geomyids. Geomyid skulls display negative allometry whereas heteromyid skulls display positive allometry. Within heteromyids, dipodomyines, and non-dipodomyines show significantly different allometric patterns. Future analyses including fossils will be necessary to test our evolutionary hypotheses.

Canalization involves mutational robustness, the lack of phenotypic change as a result of genetic mutations. Given the large divergence in phenotype across species, understanding the relationship between high robustness and evolvability has been of interest to both theorists and experimentalists. Although canalization was originally proposed in the context of multicellular organisms, the effect of multicellularity and other classes of hierarchical organization on evolvability has not been considered by theoreticians. We address this issue using a Boolean population model with explicit representation of an environment in which individuals with explicit genotype and a hierarchical phenotype representing multicellularity evolve. Robustness is described by a single real number between zero and one which emerges from the genotype-phenotype map. We find that high robustness is favoured in constant environments, and lower robustness is favoured after environmental change. Multicellularity and hierarchical organization severely constrain robustness: peak evolvability occurs at an absolute level of robustness of about 0.99 compared with values of about 0.5 in a classical neutral network model. These constraints result in a sharp peak of evolvability in which the maximum is set by the fact that the fixation of adaptive mutations becomes more improbable as robustness decreases. When robustness is put under genetic control, robustness levels leading to maximum evolvability are selected for, but maximal relative fitness appears to require recombination.

Our understanding of plant biology has been revolutionized by modern genetics and biochemistry. However, biochemical genetics can be traced back to the foundation of Mendelian genetics; indeed, one of Mendel\'s milestone discoveries of seven characteristics of pea plants later came to be ascribed to a mutation in a starch branching enzyme. Here, we review both current and historical strategies for the elucidation of plant metabolic pathways and the genes that encode their component enzymes and regulators. We use this historical review to discuss a range of classical genetic phenomena including epistasis, canalization, and heterosis as viewed through the lens of contemporary high-throughput data obtained via the array of approaches currently adopted in multiomics studies.

Developmental robustness represents the ability of an organism to resist phenotypic variations despite environmental insults and inherent genetic variations. Derailment of developmental robustness leads to phenotypic variations that can get fixed in a population for many generations. Environmental pollution is a significant worldwide problem with detrimental consequences of human development. Understanding the genetic basis for how pollutants affect human development is critical for developing interventional therapies. Here, we report that environmental stress induced by hexavalent chromium, Cr(VI), a potent industrial pollutant, compromises developmental robustness, leading to phenotypic variations in the progeny. These phenotypic variations arise due to epigenetic instability and transposon activation in the somatic tissues of the progeny rather than novel genetic mutations and can be reduced by increasing the dosage of Piwi - a Piwi-interacting RNA-binding protein, in the ovary of the exposed mother. Significantly, the derailment of developmental robustness by Cr(VI) exposure leads to tumors in the progeny, and the predisposition to develop tumors is fixed in the population for at least three generations. Thus, we show for the first time that environmental pollution can derail developmental robustness and predispose the progeny of the exposed population to develop phenotypic variations and tumors.