The details in the distributions: why and how to study phenotypic variability

SIEVE: One-stop differential expression, variability, and skewness analyses using RNA-Seq data

Motivation

RNA-Seq data analysis is commonly biased towards detecting differentially expressed genes and insufficiently conveys the complexity of gene expression changes between biological conditions. This bias arises because discrete models of RNA-Seq count data cannot fully characterize the mean, variance, and skewness of gene expression distribution using independent model parameters. A unified framework that simultaneously tests for differential expression, variability, and skewness is needed to realize the full potential of RNA-Seq data analysis in a systems biology context.

Results

We present SIEVE, a statistical methodology that provides the desired unified framework. SIEVE embraces a compositional data analysis framework that transforms discrete RNA-Seq counts to a continuous form with a distribution that is well-fitted by a skew-normal distribution. Simulation results show that SIEVE controls the false discovery rate and probability of Type II error better than existing methods for differential expression analysis. Analysis of the Mayo RNA-Seq dataset for Alzheimer's disease using SIEVE reveals that a gene set with significant expression difference in mean, standard deviation and skewness between the control and the Alzheimer's disease group strongly predicts a subject's disease state. Furthermore, functional enrichment analysis shows that relying solely on differentially expressed genes detects only a segment of a much broader spectrum of biological aspects associated with Alzheimer's disease. The latter aspects can only be revealed using genes that show differential variability and skewness. Thus, SIEVE enables fresh perspectives for understanding the intricate changes in gene expression that occur in complex diseases.

Availability

The SIEVE R package and source codes are available at https://github.com/Divo-Lee/SIEVE .

Detecting genetic effects on phenotype variability to capture gene-by-environment interactions: a systematic method comparison

Genetically associated phenotypic variability has been widely observed across organisms and traits, including in humans. Both gene-gene and gene-environment interactions can lead to an increase in genetically associated phenotypic variability. Therefore, detecting the underlying genetic variants, or variance Quantitative Trait Loci (vQTLs), can provide novel insights into complex traits. Established approaches to detect vQTLs apply different methodologies from variance-only approaches to mean-variance joint tests, but a comprehensive comparison of these methods is lacking. Here, we review available methods to detect vQTLs in humans, carry out a simulation study to assess their performance under different biological scenarios of gene-environment interactions, and apply the optimal approaches for vQTL identification to gene expression data. Overall, with a minor allele frequency (MAF) of less than 0.2, the squared residual value linear model (SVLM) and the deviation regression model (DRM) are optimal when the data follow normal and non-normal distributions, respectively. In addition, the Brown-Forsythe (BF) test is one of the optimal methods when the MAF is 0.2 or larger, irrespective of phenotype distribution. Additionally, a larger sample size and more balanced sample distribution in different exposure categories increase the power of BF, SVLM, and DRM. Our results highlight vQTL detection methods that perform optimally under realistic simulation settings and show that their relative performance depends on the phenotype distribution, allele frequency, sample size, and the type of exposure in the interaction model underlying the vQTL.

Application of unimodal probability distribution models for morphological phenotyping of budding yeast

Morphological phenotyping of the budding yeast Saccharomyces cerevisiae has helped to greatly clarify the functions of genes and increase our understanding of cellular functional networks. It is necessary to understand cell morphology and perform quantitative morphological analysis (QMA) but assigning precise values to morphological phenotypes has been challenging. We recently developed the Unimodal Morphological Data image analysis pipeline for this purpose. All true values can be estimated theoretically by applying an appropriate probability distribution if the distribution of experimental values follows a unimodal pattern. This reliable pipeline allows several downstream analyses, including detection of subtle morphological differences, selection of mutant strains with similar morphology, clustering based on morphology, and study of morphological diversity. In addition to basic research, morphological analyses of yeast cells can also be used in applied research to monitor breeding and fermentation processes and control the fermentation activity of yeast cells.

High throughput microscopy and single cell phenotypic image-based analysis in toxicology and drug discovery

Measuring single cell responses to the universe of chemicals (drugs, natural products, environmental toxicants etc.) is of paramount importance to human health as phenotypic variability in sensing stimuli is a hallmark of biology that is considered during high throughput screening. One of the ways to approach this problem is via high throughput, microscopy-based assays coupled with multi-dimensional single cell analysis methods. Here, we will summarize some of the efforts in this vast and growing field, focusing on phenotypic screens (e.g., Cell Painting), single cell analytics and quality control, with particular attention to environmental toxicology and drug screening. We will discuss advantages and limitations of high throughput assays with various end points and levels of complexity.

Acute stress reduces population-level metabolic and proteomic variation

BACKGROUND

Variation in omics data due to intrinsic biological stochasticity is often viewed as a challenging and undesirable feature of complex systems analyses. In fact, numerous statistical methods are utilized to minimize the variation among biological replicates.

RESULTS

We demonstrate that the common statistics relative standard deviation (RSD) and coefficient of variation (CV), which are often used for quality control or part of a larger pipeline in omics analyses, can also be used as a metric of a physiological stress response. Using an approach we term Replicate Variation Analysis (RVA), we demonstrate that acute physiological stress leads to feature-wide canalization of CV profiles of metabolomes and proteomes across biological replicates. Canalization is the repression of variation between replicates, which increases phenotypic similarity. Multiple in-house mass spectrometry omics datasets in addition to publicly available data were analyzed to assess changes in CV profiles in plants, animals, and microorganisms. In addition, proteomics data sets were evaluated utilizing RVA to identify functionality of reduced CV proteins.

CONCLUSIONS

RVA provides a foundation for understanding omics level shifts that occur in response to cellular stress. This approach to data analysis helps characterize stress response and recovery, and could be deployed to detect populations under stress, monitor health status, and conduct environmental monitoring.

Stochastic Fluctuations Drive Non-genetic Evolution of Proliferation in Clonal Cancer Cell Populations

Evolutionary dynamics allows us to understand many changes happening in a broad variety of biological systems, ranging from individuals to complete ecosystems. It is also behind a number of remarkable organizational changes that happen during the natural history of cancers. These reflect tumour heterogeneity, which is present at all cellular levels, including the genome, proteome and phenome, shaping its development and interrelation with its environment. An intriguing observation in different cohorts of oncological patients is that tumours exhibit an increased proliferation as the disease progresses, while the timescales involved are apparently too short for the fixation of sufficient driver mutations to promote explosive growth. Here, we discuss how phenotypic plasticity, emerging from a single genotype, may play a key role and provide a ground for a continuous acceleration of the proliferation rate of clonal populations with time. We address this question by combining the analysis of real-time growth of non-small-cell lung carcinoma cells (N-H460) together with stochastic and deterministic mathematical models that capture proliferation trait heterogeneity in clonal populations to elucidate the contribution of phenotypic transitions on tumour growth dynamics.

Genetic architecture of natural variation of cardiac performance from flies to humans

Deciphering the genetic architecture of human cardiac disorders is of fundamental importance but their underlying complexity is a major hurdle. We investigated the natural variation of cardiac performance in the sequenced inbred lines of the Drosophila Genetic Reference Panel (DGRP). Genome-wide associations studies (GWAS) identified genetic networks associated with natural variation of cardiac traits which were used to gain insights as to the molecular and cellular processes affected. Non-coding variants that we identified were used to map potential regulatory non-coding regions, which in turn were employed to predict transcription factors (TFs) binding sites. Cognate TFs, many of which themselves bear polymorphisms associated with variations of cardiac performance, were also validated by heart-specific knockdown. Additionally, we showed that the natural variations associated with variability in cardiac performance affect a set of genes overlapping those associated with average traits but through different variants in the same genes. Furthermore, we showed that phenotypic variability was also associated with natural variation of gene regulatory networks. More importantly, we documented correlations between genes associated with cardiac phenotypes in both flies and humans, which supports a conserved genetic architecture regulating adult cardiac function from arthropods to mammals. Specifically, roles for PAX9 and EGR2 in the regulation of the cardiac rhythm were established in both models, illustrating that the characteristics of natural variations in cardiac function identified in Drosophila can accelerate discovery in humans.

Parental ecological history can differentially modulate parental age effects on offspring physiological traits in Drosophila

Parents adjust their reproductive investment over their lifespan based on their condition, age, and social environment, creating the potential for inter-generational effects to differentially affect offspring physiology. To date, however, little is known about how social environments experienced by parents throughout development and adulthood influence the effect of parental age on the expression of life-history traits in the offspring. Here, I collected data on Drosophila melanogaster offspring traits (i.e., body weight, water content, and lipid reserves) from populations where either mothers, fathers both, or neither parents experienced different social environments during development (larval crowding) and adulthood. Parental treatment modulated parental age effects on offspring lipid reserves but did not influence parental age effects on offspring water content. Importantly, parents in social environments where all individuals were raised in uncrowded larval densities produced daughters and sons lighter than parental treatments which produced the heaviest offspring. The peak in offspring body weight was delayed relative to the peak in parental reproductive success, but more strongly so for daughters from parental treatments where some or all males in the parental social environments were raised in crowded larval densities (irrespective of their social context), suggesting a potential father-to-daughter effect. Overall, the findings of this study reveal that parental ecological history (here, developmental and adult social environments) can modulate the effects of parental age at reproduction on the expression of offspring traits.

Variation and Variability in Drosophila Grooming Behavior

Behavioral differences can be observed between species or populations (variation) or between individuals in a genetically similar population (variability). Here, we investigate genetic differences as a possible source of variation and variability in Drosophila grooming. Grooming confers survival and social benefits. Grooming features of five Drosophila species exposed to a dust irritant were analyzed. Aspects of grooming behavior, such as anterior to posterior progression, were conserved between and within species. However, significant differences in activity levels, proportion of time spent in different cleaning movements, and grooming syntax were identified between species. All species tested showed individual variability in the order and duration of action sequences. Genetic diversity was not found to correlate with grooming variability within a species: melanogaster flies bred to increase or decrease genetic heterogeneity exhibited similar variability in grooming syntax. Individual flies observed on consecutive days also showed grooming sequence variability. Standardization of sensory input using optogenetics reduced but did not eliminate this variability. In aggregate, these data suggest that sequence variability may be a conserved feature of grooming behavior itself. These results also demonstrate that large genetic differences result in distinguishable grooming phenotypes (variation), but that genetic heterogeneity within a population does not necessarily correspond to an increase in the range of grooming behavior (variability).

Evaluation of the phenotypic and genomic background of variability based on litter size of Large White pigs

BACKGROUND

The genetic background of trait variability has captured the interest of ecologists and animal breeders because the genes that control it could be involved in buffering various environmental effects. Phenotypic variability of a given trait can be assessed by studying the heterogeneity of the residual variance, and the quantitative trait loci (QTL) that are involved in the control of this variability are described as variance QTL (vQTL). This study focuses on litter size (total number born, TNB) and its variability in a Large White pig population. The variability of TNB was evaluated either using a simple method, i.e. analysis of the log-transformed variance of residuals (LnVar), or the more complex double hierarchical generalized linear model (DHGLM). We also performed a single-SNP (single nucleotide polymorphism) genome-wide association study (GWAS). To our knowledge, this is only the second study that reports vQTL for litter size in pigs and the first one that shows GWAS results when using two methods to evaluate variability of TNB: LnVar and DHGLM.

RESULTS

Based on LnVar, three candidate vQTL regions were detected, on Sus scrofa chromosomes (SSC) 1, 7, and 18, which comprised 18 SNPs. Based on the DHGLM, three candidate vQTL regions were detected, i.e. two on SSC7 and one on SSC11, which comprised 32 SNPs. Only one candidate vQTL region overlapped between the two methods, on SSC7, which also contained the most significant SNP. Within this vQTL region, two candidate genes were identified, ADGRF1, which is involved in neurodevelopment of the brain, and ADGRF5, which is involved in the function of the respiratory system and in vascularization. The correlation between estimated breeding values based on the two methods was 0.86. Three-fold cross-validation indicated that DHGLM yielded EBV that were much more accurate and had better prediction of missing observations than LnVar.

CONCLUSIONS

The results indicated that the LnVar and DHGLM methods resulted in genetically different traits. Based on their validation, we recommend the use of DHGLM over the simpler method of log-transformed variance of residuals. These conclusions can be useful for future studies on the evaluation of the variability of any trait in any species.

Single-cell image analysis to explore cell-to-cell heterogeneity in isogenic populations

Single-cell image analysis provides a powerful approach for studying cell-to-cell heterogeneity, which is an important attribute of isogenic cell populations, from microbial cultures to individual cells in multicellular organisms. This phenotypic variability must be explained at a mechanistic level if biologists are to fully understand cellular function and address the genotype-to-phenotype relationship. Variability in single-cell phenotypes is obscured by bulk readouts or averaging of phenotypes from individual cells in a sample; thus, single-cell image analysis enables a higher resolution view of cellular function. Here, we consider examples of both small- and large-scale studies carried out with isogenic cell populations assessed by fluorescence microscopy, and we illustrate the advantages, challenges, and the promise of quantitative single-cell image analysis.

Multiomic Big Data Analysis Challenges: Increasing Confidence in the Interpretation of Artificial Intelligence Assessments

The need for holistic molecular measurements to better understand disease initiation, development, diagnosis, and therapy has led to an increasing number of multiomic analyses. The wealth of information available from multiomic assessments, however, requires both the evaluation and interpretation of extremely large data sets, limiting analysis throughput and ease of adoption. Computational methods utilizing artificial intelligence (AI) provide the most promising way to address these challenges, yet despite the conceptual benefits of AI and its successful application in singular omic studies, the widespread use of AI in multiomic studies remains limited. Here, we discuss present and future capabilities of AI techniques in multiomic studies while introducing analytical checks and balances to validate the computational conclusions.

Developmental stochasticity and variation in floral phyllotaxis

Floral phyllotaxis is a relatively robust phenotype; trimerous and pentamerous arrangements are widely observed in monocots and core eudicots. Conversely, it also shows variability in some angiosperm clades such as 'ANA' grade (Amborellales, Nymphaeales, and Austrobaileyales), magnoliids, and Ranunculales. Regardless of the phylogenetic relationship, however, phyllotactic pattern formation appears to be a common process. What are the causes of the variability in floral phyllotaxis and how has the variation of floral phyllotaxis contributed to floral diversity? In this review, I summarize recent progress in studies on two related fields to develop answers to these questions. First, it is known that molecular and cellular stochasticity are inevitably found in biological systems, including plant development. Organisms deal with molecular stochasticity in several ways, such as dampening noise through gene networks or maintaining function through cellular redundancy. Recent studies on molecular and cellular stochasticity suggest that stochasticity is not always detrimental to plants and that it is also essential in development. Second, studies on vegetative and inflorescence phyllotaxis have shown that plants often exhibit variability and flexibility in phenotypes. Three types of phyllotaxis variations are observed, namely, fluctuation around the mean, transition between regular patterns, and a transient irregular organ arrangement called permutation. Computer models have demonstrated that stochasticity in the phyllotactic pattern formation plays a role in pattern transitions and irregularities. Variations are also found in the number and positioning of floral organs, although it is not known whether such variations provide any functional advantages. Two ways of diversification may be involved in angiosperm floral evolution: precise regulation of organ position and identity that leads to further specialization of organs and organ redundancy that leads to flexibility in floral phyllotaxis.

High-Throughput Live Imaging of Microcolonies to Measure Heterogeneity in Growth and Gene Expression

Precise measurements of between- and within-strain heterogeneity in microbial growth rates are essential for understanding genetic and environmental inputs into stress tolerance, pathogenicity, and other key components of fitness. This manuscript describes a microscope-based assay that tracks approximately 10⁵ Saccharomyces cerevisiae microcolonies per experiment. After automated time-lapse imaging of yeast immobilized in a multiwell plate, microcolony growth rates are easily analyzed with custom image-analysis software. For each microcolony, expression and localization of fluorescent proteins and survival of acute stress can also be monitored. This assay allows precise estimation of strains' average growth rates, as well as comprehensive measurement of heterogeneity in growth, gene expression, and stress tolerance within clonal populations.

Urn models for stochastic gene expression yield intuitive insights into the probability distributions of single-cell mRNA and protein counts

Fitting the probability mass functions from analytical solutions of stochastic models of gene expression to the single-cell count distributions of mRNA and protein molecules can yield valuable insights into mechanisms underlying gene expression. Solutions of chemical master equations are available for various kinetic schemes but, even for the basic ON-OFF genetic switch, they take complex forms with generating functions given as hypergeometric functions. Interpretation of gene expression dynamics in terms of bursts is not consistent with the complete range of parameters for these functions. Physical insights into the probability mass functions are essential to ensure proper interpretations but are lacking for models considering genetic switches. To fill this gap, we develop urn models for stochastic gene expression. We sample RNA polymerases or ribosomes from a master urn, which represents the cytosol, and assign them to recipient urns of two or more colors, which represent time intervals in which no switching occurs. Colors of the recipient urns represent sub-systems of the promoter states, and the assignments to urns of a specific color represent gene expression. We use elementary principles of discrete probability theory to solve a range of kinetic models without feedback, including the Peccoud-Ycart model, the Shahrezaei-Swain model, and models with an arbitrary number of promoter states. In the last case, we obtain a novel result for the protein distribution. For activated genes, we show that transcriptional lapses, which are events of gene inactivation for short time intervals separated by long active intervals, quantify the transcriptional dynamics better than bursts. We show that the intuition gained from our urn models may also be useful in understanding existing solutions for models with feedback. We contrast our models with urn models for related distributions, discuss a generalization of the Delaporte distribution for single-cell data analysis, and highlight the limitations of our models.

Advancing Antimicrobial Resistance Research Through Quantitative Modeling and Synthetic Biology

Antimicrobial resistance (AMR) is an emerging global health crisis that is undermining advances in modern medicine and, if unmitigated, threatens to kill 10 million people per year worldwide by 2050. Research over the last decade has demonstrated that the differences between genetically identical cells in the same environment can lead to drug resistance. Fluctuations in gene expression, modulated by gene regulatory networks, can lead to non-genetic heterogeneity that results in the fractional killing of microbial populations causing drug therapies to fail; this non-genetic drug resistance can enhance the probability of acquiring genetic drug resistance mutations. Mathematical models of gene networks can elucidate general principles underlying drug resistance, predict the evolution of resistance, and guide drug resistance experiments in the laboratory. Cells genetically engineered to carry synthetic gene networks regulating drug resistance genes allow for controlled, quantitative experiments on the role of non-genetic heterogeneity in the development of drug resistance. In this perspective article, we emphasize the contributions that mathematical, computational, and synthetic gene network models play in advancing our understanding of AMR to discover effective therapies against drug-resistant infections.

Heritable gene expression variability and stochasticity govern clonal heterogeneity in circadian period

A ubiquitous feature of the circadian clock across life forms is its organization as a network of cellular oscillators, with individual cellular oscillators within the network often exhibiting considerable heterogeneity in their intrinsic periods. The interaction of coupling and heterogeneity in circadian clock networks is hypothesized to influence clock’s entrainability, but our knowledge of mechanisms governing period heterogeneity within circadian clock networks remains largely elusive. In this study, we aimed to explore the principles that underlie intercellular period variation in circadian clock networks (clonal period heterogeneity). To this end, we employed a laboratory selection approach and derived a panel of 25 clonal cell populations exhibiting circadian periods ranging from 22 to 28 h. We report that a single parent clone can produce progeny clones with a wide distribution of circadian periods, and this heterogeneity, in addition to being stochastically driven, has a heritable component. By quantifying the expression of 20 circadian clock and clock-associated genes across our clone panel, we found that inheritance of expression patterns in at least three clock genes might govern clonal period heterogeneity in circadian clock networks. Furthermore, we provide evidence suggesting that heritable epigenetic variation in gene expression regulation might underlie period heterogeneity.

How do genetically identical cells exhibit a different circadian phenotype? This study reveals that a single parent clone can produce progeny with a wide distribution of circadian periods and that this heterogeneity, in addition to being stochastically driven, has a heritable component, likely via heritable epigenetic variation in gene expression regulation.

FEATURE SELECTION FOR GENERALIZED VARYING COEFFICIENT MIXED-EFFECT MODELS WITH APPLICATION TO OBESITY GWAS

Motivated by an empirical analysis of data from a genome-wide association study on obesity, measured by the body mass index (BMI), we propose a two-step gene-detection procedure for generalized varying coefficient mixed-effects models with ultrahigh dimensional covariates. The proposed procedure selects significant single nucleotide polymorphisms (SNPs) impacting the mean BMI trend, some of which have already been biologically proven to be "fat genes." The method also discovers SNPs that significantly influence the age-dependent variability of BMI. The proposed procedure takes into account individual variations of genetic effects and can also be directly applied to longitudinal data with continuous, binary or count responses. We employ Monte Carlo simulation studies to assess the performance of the proposed method and further carry out causal inference for the selected SNPs.

Inter- and intra-animal variation in the integrative properties of stellate cells in the medial entorhinal cortex

Distinctions between cell types underpin organizational principles for nervous system function. Functional variation also exists between neurons of the same type. This is exemplified by correspondence between grid cell spatial scales and the synaptic integrative properties of stellate cells (SCs) in the medial entorhinal cortex. However, we know little about how functional variability is structured either within or between individuals. Using ex-vivo patch-clamp recordings from up to 55 SCs per mouse, we found that integrative properties vary between mice and, in contrast to the modularity of grid cell spatial scales, have a continuous dorsoventral organization. Our results constrain mechanisms for modular grid firing and provide evidence for inter-animal phenotypic variability among neurons of the same type. We suggest that neuron type properties are tuned to circuit-level set points that vary within and between animals.

Single-Cell Expression Variability Implies Cell Function

As single-cell RNA sequencing (scRNA-seq) data becomes widely available, cell-to-cell variability in gene expression, or single-cell expression variability (scEV), has been increasingly appreciated. However, it remains unclear whether this variability is functionally important and, if so, what are its implications for multi-cellular organisms. Here, we analyzed multiple scRNA-seq data sets from lymphoblastoid cell lines (LCLs), lung airway epithelial cells (LAECs), and dermal fibroblasts (DFs) and, for each cell type, selected a group of homogenous cells with highly similar expression profiles. We estimated the scEV levels for genes after correcting the mean-variance dependency in that data and identified 465, 466, and 364 highly variable genes (HVGs) in LCLs, LAECs, and DFs, respectively. Functions of these HVGs were found to be enriched with those biological processes precisely relevant to the corresponding cell type’s function, from which the scRNA-seq data used to identify HVGs were generated—e.g., cytokine signaling pathways were enriched in HVGs identified in LCLs, collagen formation in LAECs, and keratinization in DFs. We repeated the same analysis with scRNA-seq data from induced pluripotent stem cells (iPSCs) and identified only 79 HVGs with no statistically significant enriched functions; the overall scEV in iPSCs was of negligible magnitude. Our results support the “variation is function” hypothesis, arguing that scEV is required for cell type-specific, higher-level system function. Thus, quantifying and characterizing scEV are of importance for our understating of normal and pathological cellular processes.

Complex Patterns of Cannabinoid Alkyl Side-Chain Inheritance in Cannabis

The cannabinoid alkyl side-chain represents an important pharmacophore, where genetic targeting of alkyl homologs has the potential to provide enhanced forms of Cannabis for biopharmaceutical manufacture. Delta(9)-tetrahydrocannabinolic acid (THCA) and cannabidiolic acid (CBDA) synthase genes govern dicyclic (CBDA) and tricyclic (THCA) cannabinoid composition. However, the inheritance of alkyl side-chain length has not been resolved, and few studies have investigated the contributions and interactions between cannabinoid synthesis pathway loci. To examine the inheritance of chemical phenotype (chemotype), THCAS and CBDAS genotypes were scored and alkyl cannabinoid segregation analysed in 210 F₂ progeny derived from a cross between two Cannabis chemotypes divergent for alkyl and cyclic cannabinoids. Inheritance patterns of F₂ progeny were non-Gaussian and deviated from Mendelian expectations. However, discrete alkyl cannabinoid segregation patterns consistent with digenic as well as epistatic modes of inheritance were observed among F₂ THCAS and CBDAS genotypes. These results suggest linkage between cannabinoid pathway loci and highlight the need for further detailed characterisation of cannabinoid inheritance to facilitate metabolic engineering of chemically elite germplasm.

Split and Merge Watershed: a two-step method for cell segmentation in fluorescence microscopy images

The development of advanced techniques in medical imaging has allowed scanning of the human body to microscopic levels, making research on cell behavior more complex and more in-depth. Recent studies have focused on cellular heterogeneity since cell-to-cell differences are always present in the cell population and this variability contains valuable information. However, identifying each cell is not an easy task because, in the images acquired from the microscope, there are clusters of cells that are touching one another. Therefore, the segmentation stage is a problem of considerable difficulty in cell image processing. Although several methods for cell segmentation are described in the literature, they have drawbacks in terms of over-segmentation, under-segmentation or misidentification. Consequently, our main motivation in studying cell segmentation was to develop a new method to achieve a good tradeoff between accurately identifying all relevant elements and not inserting segmentation artifacts. This article presents a new method for cell segmentation in fluorescence microscopy images. The proposed approach combines the well-known Marker-Controlled Watershed algorithm (MC-Watershed) with a new, two-step method based on Watershed, Split and Merge Watershed (SM-Watershed): in the first step, or split phase, the algorithm identifies the clusters using inherent characteristics of the cell, such as size and convexity, and separates them using watershed. In the second step, or the merge stage, it identifies the over-segmented regions using proper features of the cells and eliminates the divisions. Before applying our two-step method, the input image is first preprocessed, and the MC-Watershed algorithm is used to generate an initial segmented image. However, this initial result may not be suitable for subsequent tasks, such as cell count or feature extraction, because not all cells are separated, and some cells may be mistakenly confused with the background. Thus, our proposal corrects this issue with its two-step process, reaching a high performance, a suitable tradeoff between over-segmentation and under-segmentation and preserving the shape of the cell, without the need of any labeled data or relying on machine learning processes. The latter is advantageous over state-of-the-art techniques that in order to achieve similar results require labeled data, which may not be available for all of the domains. Two cell datasets were used to validate this approach, and the results were compared with other methods in the literature, using traditional metrics and quality visual assessment. We obtained 90% of average visual accuracy and an F-index higher than 80%. This proposal outperforms other techniques for cell separation, achieving an acceptable balance between over-segmentation and under-segmentation, which makes it suitable for several applications in cell identification, such as virus infection analysis, high-content cell screening, drug discovery, and morphometry.

Isolating and quantifying the role of developmental noise in generating phenotypic variation

Genotypic variation, environmental variation, and their interaction may produce variation in the developmental process and cause phenotypic differences among individuals. Developmental noise, which arises during development from stochasticity in cellular and molecular processes when genotype and environment are fixed, also contributes to phenotypic variation. While evolutionary biology has long focused on teasing apart the relative contribution of genes and environment to phenotypic variation, our understanding of the role of developmental noise has lagged due to technical difficulties in directly measuring the contribution of developmental noise. The influence of developmental noise is likely underestimated in studies of phenotypic variation due to intrinsic mechanisms within organisms that stabilize phenotypes and decrease variation. Since we are just beginning to appreciate the extent to which phenotypic variation due to stochasticity is potentially adaptive, the contribution of developmental noise to phenotypic variation must be separated and measured to fully understand its role in evolution. Here, we show that variation in the component of the developmental process corresponding to environmental and genetic factors (here treated together as a unit called the LALI-type) versus the contribution of developmental noise, can be distinguished for leopard gecko (Eublepharis macularius) head color patterns using mathematical simulations that model the role of random variation (corresponding to developmental noise) in patterning. Specifically, we modified the parameters of simulations corresponding to variation in the LALI-type to generate the full range of phenotypic variation in color pattern seen on the heads of eight leopard geckos. We observed that over the range of these parameters, variation in color pattern due to LALI-type variation exceeds that due to developmental noise in the studied gecko cohort. However, the effect of developmental noise on patterning is also substantial. Our approach addresses one of the major goals of evolutionary biology: to quantify the role of stochasticity in shaping phenotypic variation.

Magnesium Flux Modulates Ribosomes to Increase Bacterial Survival

Bacteria exhibit cell-to-cell variability in their resilience to stress, for example, following antibiotic exposure. Higher resilience is typically ascribed to "dormant" non-growing cellular states. Here, by measuring membrane potential dynamics of Bacillus subtilis cells, we show that actively growing bacteria can cope with ribosome-targeting antibiotics through an alternative mechanism based on ion flux modulation. Specifically, we observed two types of cellular behavior: growth-defective cells exhibited a mathematically predicted transient increase in membrane potential (hyperpolarization), followed by cell death, whereas growing cells lacked hyperpolarization events and showed elevated survival. Using structural perturbations of the ribosome and proteomic analysis, we uncovered that stress resilience arises from magnesium influx, which prevents hyperpolarization. Thus, ion flux modulation provides a distinct mechanism to cope with ribosomal stress. These results suggest new approaches to increase the effectiveness of ribosome-targeting antibiotics and reveal an intriguing connection between ribosomes and the membrane potential, two fundamental properties of cells.

The rise of the distributions: why non-normality is important for understanding the transcriptome and beyond

The application of statistics has been instrumental in clarifying our understanding of the genome. While insights have been derived for almost all levels of genome function, most importantly, statistics has had the greatest impact on improving our knowledge of transcriptional regulation. But the drive to extract the most meaningful inferences from big data can often force us to overlook the fundamental role that statistics plays, and specifically, the basic assumptions that we make about big data. Normality is a statistical property that is often swept up into an assumption that we may or may not be consciously aware of making. This review highlights the inherent value of non-normal distributions to big data analysis by discussing use cases of non-normality that focus on gene expression data. Collectively, these examples help to motivate the premise of why at this stage, now more than ever, non-normality is important for learning about gene regulation, transcriptomics, and more.

Mean-Variance QTL Mapping Identifies Novel QTL for Circadian Activity and Exploratory Behavior in Mice

We illustrate, through two case studies, that "mean-variance QTL mapping"-QTL mapping that models effects on the mean and the variance simultaneously-can discover QTL that traditional interval mapping cannot. Mean-variance QTL mapping is based on the double generalized linear model, which extends the standard linear model used in interval mapping by incorporating not only a set of genetic and covariate effects for mean but also set of such effects for the residual variance. Its potential for use in QTL mapping has been described previously, but it remains underutilized, with certain key advantages undemonstrated until now. In the first case study, a reduced complexity intercross of C57BL/6J and C57BL/6N mice examining circadian behavior, our reanalysis detected a mean-controlling QTL for circadian wheel running activity that interval mapping did not; mean-variance QTL mapping was more powerful than interval mapping at the QTL because it accounted for the fact that mice homozygous for the C57BL/6N allele had less residual variance than other mice. In the second case study, an intercross between C57BL/6J and C58/J mice examining anxiety-like behaviors, our reanalysis detected a variance-controlling QTL for rearing behavior; interval mapping did not identify this QTL because it does not target variance QTL. We believe that the results of these reanalyses, which in other respects largely replicated the original findings, support the use of mean-variance QTL mapping as standard practice.

Redundancy, Feedback, and Robustness in the Arabidopsis thaliana BZR/BEH Gene Family

Organismal development is remarkably robust, tolerating stochastic errors to produce consistent, so-called canalized adult phenotypes. The mechanistic underpinnings of developmental robustness are poorly understood, but recent studies implicate certain features of genetic networks such as functional redundancy, connectivity, and feedback. Here, we examine the BZR/BEH gene family, whose function contributes to embryonic stem development in the plant Arabidopsis thaliana, to test current assumptions on functional redundancy and trait robustness. Our analyses of BZR/BEH gene mutants and mutant combinations revealed that functional redundancy among these gene family members is not necessary for trait robustness. Connectivity is another commonly cited determinant of robustness; however, we found no correlation between connectivity among gene family members or their connectivity with other transcription factors and effects on developmental robustness. Instead, our data suggest that BEH4, the earliest diverged family member, modulates developmental robustness. We present evidence indicating that regulatory cross-talk among gene family members is integrated by BEH4 to promote wild-type levels of developmental robustness. Further, the chaperone HSP90, a known determinant of developmental robustness, appears to act via BEH4 in maintaining robustness of embryonic stem length. In summary, we demonstrate that even among closely related transcription factors, trait robustness can arise through the activity of a single gene family member, challenging common assumptions about the molecular underpinnings of robustness.

Control of nongenetic heterogeneity in growth rate and stress tolerance of Saccharomyces cerevisiae by cyclic AMP-regulated transcription factors

Genetically identical cells exhibit extensive phenotypic variation even under constant and benign conditions. This so-called nongenetic heterogeneity has important clinical implications: within tumors and microbial infections, cells show nongenetic heterogeneity in growth rate and in susceptibility to drugs or stress. The budding yeast, Saccharomyces cerevisiae, shows a similar form of nongenetic heterogeneity in which growth rate correlates positively with susceptibility to acute heat stress at the single-cell level. Using genetic and chemical perturbations, combined with high-throughput single-cell assays of yeast growth and gene expression, we show here that heterogeneity in intracellular cyclic AMP (cAMP) levels acting through the conserved Ras/cAMP/protein kinase A (PKA) pathway and its target transcription factors, Msn2 and Msn4, underlies this nongenetic heterogeneity. Lower levels of cAMP correspond to slower growth, as shown by direct comparison of cAMP concentration in subpopulations enriched for slower vs. faster growing cells. Concordantly, an endogenous reporter of this pathway’s activity correlates with growth in individual cells. The paralogs Msn2 and Msn4 differ in their roles in nongenetic heterogeneity in a way that demonstrates slow growth and stress tolerance are not inevitably linked. Heterogeneity in growth rate requires each, whereas only Msn2 is required for heterogeneity in expression of Tsl1, a subunit of trehalose synthase that contributes to acute-stress tolerance. Perturbing nongenetic heterogeneity by mutating genes in this pathway, or by culturing wild-type cells with the cell-permeable cAMP analog 8-bromo-cAMP or the PKA inhibitor H89, significantly impacts survival of acute heat stress. Perturbations that increase intracellular cAMP levels reduce the slower-growing subpopulation and increase susceptibility to acute heat stress, whereas PKA inhibition slows growth and decreases susceptibility to acute heat stress. Loss of Msn2 reduces, but does not completely eliminate, the correlation in individual cells between growth rate and acute-stress survival, suggesting a major role for the Msn2 pathway in nongenetic heterogeneity but also a residual benefit of slow growth. Our results shed light on the genetic control of nongenetic heterogeneity and suggest a possible means of defeating bet-hedging pathogens or tumor cells by making them more uniformly susceptible to treatment.

Nongenetic heterogeneity exists when a trait differs among individuals that have identical genotypes and environments. A clonal population can maximize its long-term success in an uncertain environment by diversifying its phenotypes via nongenetic heterogeneity: the currently unfavored ones may become the favored ones when conditions change. Nongenetic heterogeneity has clinical relevance. For example, populations of tumor cells or infectious microbes show cell-to-cell differences in growth and in drug or stress tolerance. This heterogeneity hampers efficient treatment and can potentiate harmful evolution of a tumor or pathogen. We show that in budding yeast, heterogeneity in intracellular cyclic AMP levels acting through the conserved Ras/cAMP/protein kinase A (PKA) pathway and its target transcription factors, Msn2 and Msn4, underlies the nongenetic heterogeneity of both single-cell growth rate and acute heat-stress tolerance. Perturbations of this pathway significantly affect population survival upon acute heat stress. These results illuminate a mechanism of nongenetic heterogeneity and suggest the potential value of antitumor or antifungal treatment strategies that target nongenetic heterogeneity to render the tumor or pathogen population more uniformly susceptible to a second drug that aims to kill.

A sibling method for identifying vQTLs

The propensity of a trait to vary within a population may have evolutionary, ecological, or clinical significance. In the present study we deploy sibling models to offer a novel and unbiased way to ascertain loci associated with the extent to which phenotypes vary (variance-controlling quantitative trait loci, or vQTLs). Previous methods for vQTL-mapping either exclude genetically related individuals or treat genetic relatedness among individuals as a complicating factor addressed by adjusting estimates for non-independence in phenotypes. The present method uses genetic relatedness as a tool to obtain unbiased estimates of variance effects rather than as a nuisance. The family-based approach, which utilizes random variation between siblings in minor allele counts at a locus, also allows controls for parental genotype, mean effects, and non-linear (dominance) effects that may spuriously appear to generate variation. Simulations show that the approach performs equally well as two existing methods (squared Z-score and DGLM) in controlling type I error rates when there is no unobserved confounding, and performs significantly better than these methods in the presence of small degrees of confounding. Using height and BMI as empirical applications, we investigate SNPs that alter within-family variation in height and BMI, as well as pathways that appear to be enriched. One significant SNP for BMI variability, in the MAST4 gene, replicated. Pathway analysis revealed one gene set, encoding members of several signaling pathways related to gap junction function, which appears significantly enriched for associations with within-family height variation in both datasets (while not enriched in analysis of mean levels). We recommend approximating laboratory random assignment of genotype using family data and more careful attention to the possible conflation of mean and variance effects.

Epistasis and destabilizing mutations shape gene expression variability in humans via distinct modes of action

Increasing evidence shows that phenotypic variance is genetically determined, but the underlying mechanisms of genetic control over the variance remain obscure. Here, we conducted variance-association mapping analyses to identify expression variability QTLs (evQTLs), i.e. genomic loci associated with gene expression variance, in humans. We discovered that common genetic variants may contribute to increasing gene expression variance via two distinct modes of action—epistasis and destabilization. Specifically, epistasis explains a quarter of the identified evQTLs, of which the formation is attributed to the presence of ‘third-party’ eQTLs that influence the gene expression mean in a fraction, rather than the entire set, of sampled individuals. On the other hand, the destabilization model explains the other three-quarters of evQTLs, caused by mutations that disrupt the stability of the transcription process of genes. To show the destabilizing effect, we measured discordant gene expression between monozygotic twins, and estimated the stability of gene expression in single samples using repetitive qRT-PCR assays. The mutations that cause destabilizing evQTLs were found to be associated with more pronounced expression discordance between twin pairs and less stable gene expression in single samples. Together, our results suggest that common genetic variants work either interactively or independently to shape the variability of gene expression in humans. Our findings contribute to the understanding of the mechanisms of genetic control over phenotypic variance and may have implications for the development of variance-centred analytic methods for quantitative trait mapping.

Modeling Persistent Trends in Distributions

We present a nonparametric framework to model a short sequence of probability distributions that vary both due to underlying effects of sequential progression and confounding noise. To distinguish between these two types of variation and estimate the sequential-progression effects, our approach leverages an assumption that these effects follow a persistent trend. This work is motivated by the recent rise of single-cell RNA-sequencing experiments over a brief time course, which aim to identify genes relevant to the progression of a particular biological process across diverse cell populations. While classical statistical tools focus on scalar-response regression or order-agnostic differences between distributions, it is desirable in this setting to consider both the full distributions as well as the structure imposed by their ordering. We introduce a new regression model for ordinal covariates where responses are univariate distributions and the underlying relationship reflects consistent changes in the distributions over increasing levels of the covariate. This concept is formalized as a trend in distributions, which we define as an evolution that is linear under the Wasserstein metric. Implemented via a fast alternating projections algorithm, our method exhibits numerous strengths in simulations and analyses of single-cell gene expression data.

To Produce or to Survive: How Plastic Is Your Crop Stress Physiology?

Abiotic stress causes major crop losses and is considered a greater challenge than biotic stress. Comparisons of the number of published articles and patents regarding these different types of stresses, and the number of commercially released crops designed to tolerate different types of stresses, revealed a huge gap in the bench-to-field transfer rate of abiotic stress-tolerant crops, as compared to crops designed to tolerate biotic stress. These differences underscore the complexity of abiotic stress-response mechanisms. Here, we suggest that breeding programs favoring yield-related quantitative physiological traits (QPTs; e.g., photosynthesis rate or stomatal conductance) have canalized those QPTs at their highest levels. This has affected the sensitivity of those QPTs to changing environmental conditions and those traits have become less plastic. We also suggest that breeding pressure has had an asymmetric impact on different QPTs, depending on their sensitivity to environmental conditions and their interactions with other QPTs. We demonstrate this asymmetric impact on the regulation of whole-plant water balance, showing how plastic membrane water content, stomatal conductance and leaf hydraulic conductance interact to canalize whole-organ water content. We suggest that a QPT's plasticity is itself an important trait and that understanding this plasticity may help us to develop yield-optimized crops.

Stochastic loss and gain of symmetric divisions in the C. elegans epidermis perturbs robustness of stem cell number

Biological systems are subject to inherent stochasticity. Nevertheless, development is remarkably robust, ensuring the consistency of key phenotypic traits such as correct cell numbers in a certain tissue. It is currently unclear which genes modulate phenotypic variability, what their relationship is to core components of developmental gene networks, and what is the developmental basis of variable phenotypes. Here, we start addressing these questions using the robust number of Caenorhabditis elegans epidermal stem cells, known as seam cells, as a readout. We employ genetics, cell lineage tracing, and single molecule imaging to show that mutations in lin-22, a Hes-related basic helix-loop-helix (bHLH) transcription factor, increase seam cell number variability. We show that the increase in phenotypic variability is due to stochastic conversion of normally symmetric cell divisions to asymmetric and vice versa during development, which affect the terminal seam cell number in opposing directions. We demonstrate that LIN-22 acts within the epidermal gene network to antagonise the Wnt signalling pathway. However, lin-22 mutants exhibit cell-to-cell variability in Wnt pathway activation, which correlates with and may drive phenotypic variability. Our study demonstrates the feasibility to study phenotypic trait variance in tractable model organisms using unbiased mutagenesis screens.

Engineered Control of Genetic Variability Reveals Interplay among Quorum Sensing, Feedback Regulation, and Biochemical Noise

Stochastic fluctuations in gene expression trigger both beneficial and harmful consequences for cell behavior. Therefore, achieving a desired mean protein expression level while minimizing noise is of interest in many applications, including robust protein production systems in industrial biotechnology. Here, we consider a synthetic gene circuit combining intracellular negative feedback and cell-to-cell communication based on quorum sensing. Accounting for both intrinsic and extrinsic noise, stochastic simulations allow us to analyze the capability of the circuit to reduce noise strength as a function of its parameters. We obtain mean expression levels and noise strengths for all species under different scenarios, showing good agreement with system-wide available experimental data of protein abundance and noise in Escherichia coli. Our in silico experiments, validated by preliminary in vivo results, reveal significant noise attenuation in gene expression through the interplay between quorum sensing and negative feedback and highlight the differential role that they play in regard to intrinsic and extrinsic noise.

Instructed subsets or agile swarms: how T-helper cells may adaptively counter uncertainty with variability and plasticity

Over recent years, extensive phenotypic variability and plasticity have been revealed among the T-helper cells of the mammalian adaptive immune system, even within clonal lineages of identical antigen specificity. This challenges the conventional view that T-helper cells assort into functionally distinct subsets following differential instruction by the innate immune system. We argue that the adaptive value of coping with uncertainty can reconcile the 'instructed subset' framework with T-helper variability and plasticity. However, we also suggest that T-helper cells might better be understood as agile swarms engaged in collective decision-making to promote host fitness. With rigorous testing, the 'agile swarms' framework may illuminate how variable and plastic individual T-helper cells interact to create coherent immunity.

Heterogeneity in Pure Microbial Systems: Experimental Measurements and Modeling

Cellular heterogeneity influences bioprocess performance in ways that until date are not completely elucidated. In order to account for this phenomenon in the design and operation of bioprocesses, reliable analytical and mathematical descriptions are required. We present an overview of the single cell analysis, and the mathematical modeling frameworks that have potential to be used in bioprocess control and optimization, in particular for microbial processes. In order to be suitable for bioprocess monitoring, experimental methods need to be high throughput and to require relatively short processing time. One such method used successfully under dynamic conditions is flow cytometry. Population balance and individual based models are suitable modeling options, the latter one having in particular a good potential to integrate the various data collected through experimentation. This will be highly beneficial for appropriate process design and scale up as a more rigorous approach may prevent a priori unwanted performance losses. It will also help progressing synthetic biology applications to industrial scale.

Genetics of Genome-Wide Recombination Rate Evolution in Mice from an Isolated Island

Recombination rate is a heritable quantitative trait that evolves despite the fundamentally conserved role that recombination plays in meiosis. Differences in recombination rate can alter the landscape of the genome and the genetic diversity of populations. Yet our understanding of the genetic basis of recombination rate evolution in nature remains limited. We used wild house mice (Mus musculus domesticus) from Gough Island (GI), which diverged recently from their mainland counterparts, to characterize the genetics of recombination rate evolution. We quantified genome-wide autosomal recombination rates by immunofluorescence cytology in spermatocytes from 240 F2 males generated from intercrosses between GI-derived mice and the wild-derived inbred strain WSB/EiJ. We identified four quantitative trait loci (QTL) responsible for inter-F2 variation in this trait, the strongest of which had effects that opposed the direction of the parental trait differences. Candidate genes and mutations for these QTL were identified by overlapping the detected intervals with whole-genome sequencing data and publicly available transcriptomic profiles from spermatocytes. Combined with existing studies, our findings suggest that genome-wide recombination rate divergence is not directional and its evolution within and between subspecies proceeds from distinct genetic loci.

The evolutionary advantage of heritable phenotypic heterogeneity

Phenotypic plasticity is an evolutionary driving force in diverse biological processes, including the adaptive immune system, the development of neoplasms, and the persistence of pathogens despite drug pressure. It is essential, therefore, to understand the evolutionary advantage of an allele that confers on cells the ability to express a range of phenotypes. Here, we study the fate of a new mutation that allows the expression of multiple phenotypic states, introduced into a finite population of individuals that can express only a single phenotype. We show that the advantage of such a mutation depends on the degree of phenotypic heritability between generations, called phenotypic memory. We analyze the fixation probability of the phenotypically plastic allele as a function of phenotypic memory, the variance of expressible phenotypes, the rate of environmental changes, and the population size. We find that the fate of a phenotypically plastic allele depends fundamentally on the environmental regime. In constant environments, plastic alleles are advantageous and their fixation probability increases with the degree of phenotypic memory. In periodically fluctuating environments, by contrast, there is an optimum phenotypic memory that maximizes the probability of the plastic allele's fixation. This same optimum memory also maximizes geometric mean fitness, in steady state. We interpret these results in the context of previous studies in an infinite-population framework. We also discuss the implications of our results for the design of therapies that can overcome persistence and, indirectly, drug resistance.

Resolving the Complex Genetic Basis of Phenotypic Variation and Variability of Cellular Growth

In all organisms, the majority of traits vary continuously between individuals. Explaining the genetic basis of quantitative trait variation requires comprehensively accounting for genetic and nongenetic factors as well as their interactions. The growth of microbial cells can be characterized by a lag duration, an exponential growth phase, and a stationary phase. Parameters that characterize these growth phases can vary among genotypes (phenotypic variation), environmental conditions (phenotypic plasticity), and among isogenic cells in a given environment (phenotypic variability). We used a high-throughput microscopy assay to map genetic loci determining variation in lag duration and exponential growth rate in growth rate-limiting and nonlimiting glucose concentrations, using segregants from a cross of two natural isolates of the budding yeast, Saccharomyces cerevisiae We find that some quantitative trait loci (QTL) are common between traits and environments whereas some are unique, exhibiting gene-by-environment interactions. Furthermore, whereas variation in the central tendency of growth rate or lag duration is explained by many additive loci, differences in phenotypic variability are primarily the result of genetic interactions. We used bulk segregant mapping to increase QTL resolution by performing whole-genome sequencing of complex mixtures of an advanced intercross mapping population grown in selective conditions using glucose-limited chemostats. We find that sequence variation in the high-affinity glucose transporter HXT7 contributes to variation in growth rate and lag duration. Allele replacements of the entire locus, as well as of a single polymorphic amino acid, reveal that the effect of variation in HXT7 depends on genetic, and allelic, background. Amplifications of HXT7 are frequently selected in experimental evolution in glucose-limited environments, but we find that HXT7 amplifications result in antagonistic pleiotropy that is absent in naturally occurring variants of HXT7 Our study highlights the complex nature of the genotype-to-phenotype map within and between environments.

Integrating mean and variance heterogeneities to identify differentially expressed genes

BACKGROUND

In functional genomics studies, tests on mean heterogeneity have been widely employed to identify differentially expressed genes with distinct mean expression levels under different experimental conditions. Variance heterogeneity (aka, the difference between condition-specific variances) of gene expression levels is simply neglected or calibrated for as an impediment. The mean heterogeneity in the expression level of a gene reflects one aspect of its distribution alteration; and variance heterogeneity induced by condition change may reflect another aspect. Change in condition may alter both mean and some higher-order characteristics of the distributions of expression levels of susceptible genes.

RESULTS

In this report, we put forth a conception of mean-variance differentially expressed (MVDE) genes, whose expression means and variances are sensitive to the change in experimental condition. We mathematically proved the null independence of existent mean heterogeneity tests and variance heterogeneity tests. Based on the independence, we proposed an integrative mean-variance test (IMVT) to combine gene-wise mean heterogeneity and variance heterogeneity induced by condition change. The IMVT outperformed its competitors under comprehensive simulations of normality and Laplace settings. For moderate samples, the IMVT well controlled type I error rates, and so did existent mean heterogeneity test (i.e., the Welch t test (WT), the moderated Welch t test (MWT)) and the procedure of separate tests on mean and variance heterogeneities (SMVT), but the likelihood ratio test (LRT) severely inflated type I error rates. In presence of variance heterogeneity, the IMVT appeared noticeably more powerful than all the valid mean heterogeneity tests. Application to the gene profiles of peripheral circulating B raised solid evidence of informative variance heterogeneity. After adjusting for background data structure, the IMVT replicated previous discoveries and identified novel experiment-wide significant MVDE genes.

CONCLUSIONS

Our results indicate tremendous potential gain of integrating informative variance heterogeneity after adjusting for global confounders and background data structure. The proposed informative integration test better summarizes the impacts of condition change on expression distributions of susceptible genes than do the existent competitors. Therefore, particular attention should be paid to explicitly exploit the variance heterogeneity induced by condition change in functional genomics analysis.

Selection Transforms the Landscape of Genetic Variation Interacting with Hsp90

The protein-folding chaperone Hsp90 has been proposed to buffer the phenotypic effects of mutations. The potential for Hsp90 and other putative buffers to increase robustness to mutation has had major impact on disease models, quantitative genetics, and evolutionary theory. But Hsp90 sometimes contradicts expectations for a buffer by potentiating rapid phenotypic changes that would otherwise not occur. Here, we quantify Hsp90’s ability to buffer or potentiate (i.e., diminish or enhance) the effects of genetic variation on single-cell morphological features in budding yeast. We corroborate reports that Hsp90 tends to buffer the effects of standing genetic variation in natural populations. However, we demonstrate that Hsp90 tends to have the opposite effect on genetic variation that has experienced reduced selection pressure. Specifically, Hsp90 tends to enhance, rather than diminish, the effects of spontaneous mutations and recombinations. This result implies that Hsp90 does not make phenotypes more robust to the effects of genetic perturbation. Instead, natural selection preferentially allows buffered alleles to persist and thereby creates the false impression that Hsp90 confers greater robustness.

Most biologists appreciate that natural selection filters new mutations (e.g., by eliminating deleterious ones), such that genetic variation in nature is biased. The idea that selection also skews the types of genetic interactions that exist in nature is less appreciated. For example, studies spanning diverse species have shown that the protein Hsp90, which helps other proteins to fold properly, tends to diminish the observable effects of genetic variation. This observation has led to the assumption that Hsp90 also buffers the effects of new mutations. This untested assumption has served as a rationale for cancer-treatment strategies and shaped our understanding of variation in complex traits. We measured the effects of new mutations on the shapes and sizes of individual yeast cells and found that Hsp90 does not tend to buffer these effects. Instead, Hsp90 interacts with new mutations in diverse ways, sometimes buffering, but more often enhancing mutational effects on cell shape and size. We conclude that selection preferentially allows buffered mutations to persist in natural populations. This result alters common perceptions about why cryptic (i.e., buffered) genetic variation exists and casts doubt on cancer-treatment strategies aiming to target presumed buffers of mutational effects.

Physiology of Saccharomyces cerevisiae strains isolated from Brazilian biomes: new insights into biodiversity and industrial applications

Fourteen indigenous Saccharomyces cerevisiae strains isolated from the barks of three tree species located in the Atlantic Rain Forest and Cerrado biomes in Brazil were genetically and physiologically compared to laboratory strains and to strains from the Brazilian fuel ethanol industry. Although no clear correlation could be found either between phenotype and isolation spot or between phenotype and genomic lineage, a set of indigenous strains with superior industrially relevant traits over commonly known industrial and laboratory strains was identified: strain UFMG-CM-Y257 has a very high specific growth rate on sucrose (0.57 ± 0.02 h^-1), high ethanol yield (1.65 ± 0.02 mol ethanol mol hexose equivalent^-1), high ethanol productivity (0.19 ± 0.00 mol L^-1 h^-1), high tolerance to acetic acid (10 g L^-1) and to high temperature (40°C). Strain UFMG-CM-Y260 displayed high ethanol yield (1.67 ± 0.13 mol ethanol mol hexose equivalent^-1), high tolerance to ethanol and to low pH, a trait which is important for non-aseptic industrial processes. Strain UFMG-CM-Y267 showed high tolerance to acetic acid and to high temperature (40°C), which is of particular interest to second generation industrial processes.

Neuromodulatory Regulation of Behavioral Individuality in Zebrafish

Inter-individual behavioral variation is thought to increase fitness and aid adaptation to environmental change, but the underlying mechanisms are poorly understood. We find that variation between individuals in neuromodulatory input contributes to individuality in short-term habituation of the zebrafish (Danio Rerio) acoustic startle response (ASR). ASR habituation varies greatly between individuals, but differences are stable over days and are heritable. Acoustic stimuli that activate ASR-command Mauthner cells also activate dorsal raphe nucleus (DRN) serotonergic neurons, which project to the vicinity of the Mauthner cells and their inputs. DRN neuron activity decreases during habituation in proportion to habituation and a genetic manipulation that reduces serotonin content in DRN neurons increases habituation, whereas serotonergic agonism or DRN activation with ChR2 reduces habituation. Finally, level of rundown of DRN activity co-segregates with extent of behavioral habituation across generations. Thus, variation between individuals in neuromodulatory input contributes to individuality in a core adaptive behavior. VIDEO ABSTRACT.

The developmental genetics of biological robustness

BACKGROUND

Living organisms are continuously confronted with perturbations, such as environmental changes that include fluctuations in temperature and nutrient availability, or genetic changes such as mutations. While some developmental systems are affected by such challenges and display variation in phenotypic traits, others continue consistently to produce invariable phenotypes despite perturbation. This ability of a living system to maintain an invariable phenotype in the face of perturbations is termed developmental robustness. Biological robustness is a phenomenon observed across phyla, and studying its mechanisms is central to deciphering the genotype-phenotype relationship. Recent work in yeast, animals and plants has shown that robustness is genetically controlled and has started to reveal the underlying mechinisms behind it.

SCOPE AND CONCLUSIONS

Studying biological robustness involves focusing on an important property of developmental traits, which is the phenotypic distribution within a population. This is often neglected because the vast majority of developmental biology studies instead focus on population aggregates, such as trait averages. By drawing on findings in animals and yeast, this Viewpoint considers how studies on plant developmental robustness may benefit from strict definitions of what is the developmental system of choice and what is the relevant perturbation, and also from clear distinctions between gene effects on the trait mean and the trait variance. Recent advances in quantitative developmental biology and high-throughput phenotyping now allow the design of targeted genetic screens to identify genes that amplify or restrict developmental trait variance and to study how variation propagates across different phenotypic levels in biological systems. The molecular characterization of more quantitative trait loci affecting trait variance will provide further insights into the evolution of genes modulating developmental robustness. The study of robustness mechanisms in closely related species will address whether mechanisms of robustness are evolutionarily conserved.

Applications and implications of the exponentially modified gamma distribution as a model for time variabilities related to cell proliferation and gene expression

A panel of published distributions of cell interdivision times (IDT) comprising 77 datasets related to 16 cell types, some studied under different conditions, was used to evaluate their conformance to the exponentially modified gamma distribution (EMGD) in comparison with distributions suggested for IDT data earlier. Lognormal, gamma, inverse Gaussian, and shifted Weibull and gamma distributions were found to be generally inferior to EMGD. Exponentially modified Gaussian (EMG) performed equally well. Although EMGD or EMG may be worse than some other distributions in specific cases, the reason that IDT distributions must be generated by a common mechanism of the cell cycle makes it unlikely that they differ essentially in different cell types. Therefore, exponentially modified peak functions, such as EMGD or EMG, are most appropriate if the use of a single distribution for IDT data is reasonable. EMGD is also shown to be the best descriptive tool for published data on the distribution of times between the bursts of mRNA synthesis at defined genes in single cells. EMG is inadequate to such data because its Gaussian component markedly extends to the negative time domain. The applicability of EMGD to comparable features of cells and genes behaviors are discussed to support the validity of the transition probability model and to relate the exponential component of EMGD to the times of cell dwelling in the restriction point of the cell cycle.

Genome-wide association study reveals novel loci for litter size and its variability in a Large White pig population

Background

In many traits, not only individual trait levels are under genetic control, but also the variation around that level. In other words, genotypes do not only differ in mean, but also in (residual) variation around the genotypic mean. New statistical methods facilitate gaining knowledge on the genetic architecture of complex traits such as phenotypic variability. Here we study litter size (total number born) and its variation in a Large White pig population using a Double Hierarchical Generalized Linear model, and perform a genome-wide association study using a Bayesian method.

Results

In total, 10 significant single nucleotide polymorphisms (SNPs) were detected for total number born (TNB) and 9 SNPs for variability of TNB (varTNB). Those SNPs explained 0.83 % of genetic variance in TNB and 1.44 % in varTNB. The most significant SNP for TNB was detected on Sus scrofa chromosome (SSC) 11. A possible candidate gene for TNB is ENOX1, which is involved in cell growth and survival. On SSC7, two possible candidate genes for varTNB are located. The first gene is coding a swine heat shock protein 90 (HSPCB = Hsp90), which is a well-studied gene stabilizing morphological traits in Drosophila and Arabidopsis. The second gene is VEGFA, which is activated in angiogenesis and vasculogenesis in the fetus. Furthermore, the genetic correlation between additive genetic effects on TNB and on its variation was 0.49. This indicates that the current selection to increase TNB will also increase the varTNB.

Conclusions

To the best of our knowledge, this is the first study reporting SNPs associated with variation of a trait in pigs. Detected genomic regions associated with varTNB can be used in genomic selection to decrease varTNB, which is highly desirable to avoid very small or very large litters in pigs. However, the percentage of variance explained by those regions was small. The SNPs detected in this study can be used as indication for regions in the Sus scrofa genome involved in maintaining low variability of litter size, but further studies are needed to identify the causative loci.