The Mass-Longevity Triangle: Pareto Optimality and the Geometry of Life-History Trait Space

When organisms need to perform multiple tasks they face a fundamental tradeoff: no phenotype can be optimal at all tasks. This situation was recently analyzed using Pareto optimality, showing that tradeoffs between tasks lead to phenotypes distributed on low dimensional polygons in trait space. The vertices of these polygons are archetypes—phenotypes optimal at a single task. This theory was applied to examples from animal morphology and gene expression. Here we ask whether Pareto optimality theory can apply to life history traits, which include longevity, fecundity and mass. To comprehensively explore the geometry of life history trait space, we analyze a dataset of life history traits of 2105 endothermic species. We find that, to a first approximation, life history traits fall on a triangle in log-mass log-longevity space. The vertices of the triangle suggest three archetypal strategies, exemplified by bats, shrews and whales, with specialists near the vertices and generalists in the middle of the triangle. To a second approximation, the data lies in a tetrahedron, whose extra vertex above the mass-longevity triangle suggests a fourth strategy related to carnivory. Each animal species can thus be placed in a coordinate system according to its distance from the archetypes, which may be useful for genome-scale comparative studies of mammalian aging and other biological aspects. We further demonstrate that Pareto optimality can explain a range of previous studies which found animal and plant phenotypes which lie in triangles in trait space. This study demonstrates the applicability of multi-objective optimization principles to understand life history traits and to infer archetypal strategies that suggest why some mammalian species live much longer than others of similar mass.     

A Mechanistic Collective Cell Model for Epithelial Colony Growth and Contact Inhibition

We present a mechanistic hybrid continuum-discrete model to simulate the dynamics of epithelial cell colonies. Collective cell dynamics are modeled using continuum equations that capture plastic, viscoelastic, and elastic deformations in the clusters while providing single-cell resolution. The continuum equations can be viewed as a coarse-grained version of previously developed discrete models that treat epithelial clusters as a two-dimensional network of vertices or stochastic interacting particles and follow the framework of dynamic density functional theory appropriately modified to account for cell size and shape variability. The discrete component of the model implements cell division and thus influences cell size and shape that couple to the continuum component. The model is validated against recent in vitro studies of epithelial cell colonies using Madin-Darby canine kidney cells. In good agreement with experiments, we find that mechanical interactions and constraints on the local expansion of cell size cause inhibition of cell motion and reductive cell division. This leads to successively smaller cells and a transition from exponential to quadratic growth of the colony that is associated with a constant-thickness rim of growing cells at the cluster edge, as well as the emergence of short-range ordering and solid-like behavior. A detailed analysis of the model reveals a scale invariance of the growth and provides insight into the generation of stresses and their influence on the dynamics of the colonies. Compared to previous models, our approach has several advantages: it is independent of dimension, it can be parameterized using classical elastic properties (Poisson’s ratio and Young’s modulus), and it can easily be extended to incorporate multiple cell types and general substrate geometries.     

Positional specificity within plant cells

Plant cells are clearly definable polyhedra composed of discrete facets, edges and vertices. It is proposed that these three elements of geometry along with the mitotic spindle and certain cytoplasmic particles constitute an ensemble of structural entities with morphogenetic capacities. Facets, edges and vertices are programmed to assume various states as dictated by the nucleus, and according to the combination of these states within a cell, the cell will express specific patterns of growth, thickening and separation from adjacent cells.A number of examples of cell types in plants and their origins from generalized types of cells are examined by this model and it is recognized that each cell type does not represent a unique behavior but rather expresses a particular combination of changes from a total of about 20 fundamental processes. These five elements of form are viewed as composed of polymers with stereospecific receptors except that the elements collectively are an architectural continuum having both structural and formative features.     

The geometry of the Pareto front in biological phenotype space

When organisms perform a single task, selection leads to phenotypes that maximize performance at that task. When organisms need to perform multiple tasks, a trade‐off arises because no phenotype can optimize all tasks. Recent work addressed this question, and assumed that the performance at each task decays with distance in trait space from the best phenotype at that task. Under this assumption, the best‐fitness solutions (termed the Pareto front) lie on simple low‐dimensional shapes in trait space: line segments, triangles and other polygons. The vertices of these polygons are specialists at a single task. Here, we generalize this finding, by considering performance functions of general form, not necessarily functions that decay monotonically with distance from their peak. We find that, except for performance functions with highly eccentric contours, simple shapes in phenotype space are still found, but with mildly curving edges instead of straight ones. In a wide range of systems, complex data on multiple quantitative traits, which might be expected to fill a high‐dimensional phenotype space, is predicted instead to collapse onto low‐dimensional shapes; phenotypes near the vertices of these shapes are predicted to be specialists, and can thus suggest which tasks may be at play.     

GDSCalc: A Web-Based Application for Evaluating Discrete Graph Dynamical Systems

Discrete dynamical systems are used to model various realistic systems in network science, from social unrest in human populations to regulation in biological networks. A common approach is to model the agents of a system as vertices of a graph, and the pairwise interactions between agents as edges. Agents are in one of a finite set of states at each discrete time step and are assigned functions that describe how their states change based on neighborhood relations. Full characterization of state transitions of one system can give insights into fundamental behaviors of other dynamical systems. In this paper, we describe a discrete graph dynamical systems (GDSs) application called GDSCalc for computing and characterizing system dynamics. It is an open access system that is used through a web interface. We provide an overview of GDS theory. This theory is the basis of the web application; i.e., an understanding of GDS provides an understanding of the software features, while abstracting away implementation details. We present a set of illustrative examples to demonstrate its use in education and research. Finally, we compare GDSCalc with other discrete dynamical system software tools. Our perspective is that no single software tool will perform all computations that may be required by all users; tools typically have particular features that are more suitable for some tasks. We situate GDSCalc within this space of software tools.     

Single-Cell Gene Expression Analysis: Implications for Neurodegenerative and Neuropsychiatric Disorders

Technical and experimental advances in microaspiration techniques, RNA amplification, quantitative real-time polymerase chain reaction (qPCR), and cDNA microarray analysis have led to an increase in the number of studies of single-cell gene expression. In particular, the central nervous system (CNS) is an ideal structure to apply single-cell gene expression paradigms. Unlike an organ that is composed of one principal cell type, the brain contains a constellation of neuronal and noneuronal populations of cells. A goal is to sample gene expression from similar cell types within a defined region without potential contamination by expression profiles of adjacent neuronal subpopulations and noneuronal cells. The unprecedented resolution afforded by single-cell RNA analysis in combination with cDNA microarrays and qPCR-based analyses allows for relative gene expression level comparisons across cell types under different experimental conditions and disease states. The ability to analyze single cells is an important distinction from global and regional assessments of mRNA expression and can be applied to optimally prepared tissues from animal models as well as postmortem human brain tissues. This focused review illustrates the potential power of single-cell gene expression studies within the CNS in relation to neurodegenerative and neuropsychiatric disorders such as Alzheimer's disease (AD) and schizophrenia, respectively.     

Tools for Single-Cell Kinetic Analysis of Virus-Host Interactions

Measures of cellular gene expression or behavior, when performed on individual cells, inevitably reveal a diversity of behaviors and outcomes that can correlate with normal or diseased states. For virus infections, the potential diversity of outcomes are pushed to an extreme, where measures of infection reflect features of the specific infecting virus particle, the individual host cell, as well as interactions between viral and cellular components. Single-cell measures, while revealing, still often rely on specialized fluid handling capabilities, employ end-point measures, and remain labor-intensive to perform. To address these limitations, we consider a new microwell-based device that uses simple pipette-based fluid handling to isolate individual cells. Our design allows different experimental conditions to be implemented in a single device, permitting easier and more standardized protocols. Further, we utilize a recently reported dual-color fluorescent reporter system that provides dynamic readouts of viral and cellular gene expression during single-cell infections by vesicular stomatitis virus. In addition, we develop and show how free, open-source software can enable streamlined data management and batch image analysis. Here we validate the integration of the device and software using the reporter system to demonstrate unique single-cell dynamic measures of cellular responses to viral infection.     

scLM: Automatic Detection of Consensus Gene Clusters Across Multiple Single-cell Datasets

In gene expression profiling studies, including single-cell RNA sequencing(sc RNA-seq)analyses, the identification and characterization of co-expressed genes provides critical information on cell identity and function. Gene co-expression clustering in sc RNA-seq data presents certain challenges. We show that commonly used methods for single-cell data are not capable of identifying co-expressed genes accurately, and produce results that substantially limit biological expectations of co-expressed genes. Herein, we present single-cell Latent-variable Model(sc LM), a gene coclustering algorithm tailored to single-cell data that performs well at detecting gene clusters with significant biologic context. Importantly, sc LM can simultaneously cluster multiple single-cell datasets, i.e., consensus clustering, enabling users to leverage single-cell data from multiple sources for novel comparative analysis. sc LM takes raw count data as input and preserves biological variation without being influenced by batch effects from multiple datasets. Results from both simulation data and experimental data demonstrate that sc LM outperforms the existing methods with considerably improved accuracy. To illustrate the biological insights of sc LM, we apply it to our in-house and public experimental sc RNA-seq datasets. sc LM identifies novel functional gene modules and refines cell states, which facilitates mechanism discovery and understanding of complex biosystems such as cancers. A user-friendly R package with all the key features of the sc LM method is available at https://github.com/QSong-github/sc LM.     

Cell type diversity in scallop adductor muscles revealed by single-cell RNA-Seq

Studies on cell atlas in marine invertebrates provide a better understanding of cell types, stem cell maintenance, and lineages of cell differentiation. To investigate the molecular features of various cell types in molluscan muscles, we performed single-cell RNA sequencing (scRNA-seq) to map cell types in scallop adductor muscles. We uncovered the cell type-specific features of 20 cell clusters defined by the expression of multiple specific molecular markers. These cell clusters are mainly classified into four broad classes, including mesenchymal stem cells, muscle cells, neurons, and haemolymph cells. In particular, we identified a diverse repertoire of neurons in the striated adductor muscle, but not in the smooth muscle. We further reconstructed the cell differentiation events using all the cell clusters by single-cell pseudotemporal trajectories. By integrating dual BrdU-PCNA immunodetection, neuron-specific staining and electron microscopy observation, we showed the spatial distribution of mesenchymal stem cells and neurons in striated adductor muscle of scallops. The present findings will not only be useful to address the cell type-specific gene expression profiles in scallop muscles, but also provide valuable resources for cross-species comparison of marine organisms.     

Maturation of cellular Fe-S proteins: an essential function of mitochondria

Iron-sulfur (Fe-S) cluster-containing proteins perform important tasks in catalysis, electron transfer and regulation of gene expression. In eukaryotes, mitochondria are the primary site of cluster formation of most Fe-S proteins. Assembly of the Fe-S clusters is mediated by the iron-sulphate cluster assembly (ISC) machinery consisting of some ten proteins.     

Single-cell exome sequencing and monoclonal evolution of a JAK2-negative myeloproliferative neoplasm

Tumor heterogeneity presents a challenge for inferring clonal evolution and driver gene identification. Here, we describe a method for analyzing the cancer genome at a single-cell nucleotide level. To perform our analyses, we first devised and validated a high-throughput whole-genome single-cell sequencing method using two lymphoblastoid cell line single cells. We then carried out whole-exome single-cell sequencing of 90 cells from a JAK2-negative myeloproliferative neoplasm patient. The sequencing data from 58 cells passed our quality control criteria, and these data indicated that this neoplasm represented a monoclonal evolution. We further identified essential thrombocythemia (ET)-related candidate mutations such as SESN2 and NTRK1, which may be involved in neoplasm progression. This pilot study allowed the initial characterization of the disease-related genetic architecture at the single-cell nucleotide level. Further, we established a single-cell sequencing method that opens the way for detailed analyses of a variety of tumor types, including those with high genetic complex between patients.     

INSISTC: Incorporating network structure information for single-cell type classification

Uncovering gene regulatory mechanisms in individual cells can provide insight into cell heterogeneity and function. Recent accumulated Single-Cell RNA-Seq data have made it possible to analyze gene regulation at single-cell resolution. Understanding cell-type-specific gene regulation can assist in more accurate cell type and state identification. Computational approaches utilizing such relationships are under development. Methods pioneering in integrating gene regulatory mechanism discovery with cell-type classification encounter challenges such as determine gene regulatory relationships and incorporate gene regulatory network structure. To fill this gap, we developed INSISTC, a computational method to incorporate gene regulatory network structure information for single-cell type classification. INSISTC is capable of identifying cell-type-specific gene regulatory mechanisms while performing single-cell type classification. INSISTC demonstrated its accuracy in cell type classification and its potential for providing insight into molecular mechanisms specific to individual cells. In comparison with the alternative methods, INSISTC demonstrated its complementary performance for gene regulation interpretation.