Studying stochastic systems biology of the cell with single-cell genomics data.
Ontology highlight
ABSTRACT:
SUBMITTER: Gorin G
PROVIDER: S-EPMC10245677 | biostudies-literature | 2023 May
REPOSITORIES: biostudies-literature
Publications
Studying stochastic systems biology of the cell with single-cell genomics data.
Gorin Gennady G Vastola John J JJ Pachter Lior L
bioRxiv : the preprint server for biology 20230529
Recent experimental developments in genome-wide RNA quantification hold considerable promise for systems biology. However, rigorously probing the biology of living cells requires a unified mathematical framework that accounts for single-molecule biological stochasticity in the context of technical variation associated with genomics assays. We review models for a variety of RNA transcription processes, as well as the encapsulation and library construction steps of microfluidics-based single-cell ...[more]
Similar Datasets
Project description:Recent experimental developments in genome-wide RNA quantification hold considerable promise for systems biology. However, rigorously probing the biology of living cells requires a unified mathematical framework that accounts for single-molecule biological stochasticity in the context of technical variation associated with genomics assays. We review models for a variety of RNA transcription processes, as well as the encapsulation and library construction steps of microfluidics-based single-cell RNA sequencing, and present a framework to integrate these phenomena by the manipulation of generating functions. Finally, we use simulated scenarios and biological data to illustrate the implications and applications of the approach.
| S-EPMC10725240 | biostudies-literature
Project description:Cellular circuits sense the environment, process signals, and compute decisions using networks of interacting proteins. To model such a system, the abundance of each activated protein species can be described as a stochastic function of the abundance of other proteins. High-dimensional single-cell technologies, such as mass cytometry, offer an opportunity to characterize signaling circuit-wide. However, the challenge of developing and applying computational approaches to interpret such complex data remains. Here, we developed computational methods, based on established statistical concepts, to characterize signaling network relationships by quantifying the strengths of network edges and deriving signaling response functions. In comparing signaling between naïve and antigen-exposed CD4(+) T lymphocytes, we find that although these two cell subtypes had similarly wired networks, naïve cells transmitted more information along a key signaling cascade than did antigen-exposed cells. We validated our characterization on mice lacking the extracellular-regulated mitogen-activated protein kinase (MAPK) ERK2, which showed stronger influence of pERK on pS6 (phosphorylated-ribosomal protein S6), in naïve cells as compared with antigen-exposed cells, as predicted. We demonstrate that by using cell-to-cell variation inherent in single-cell data, we can derive response functions underlying molecular circuits and drive the understanding of how cells process signals.
| S-EPMC4334155 | biostudies-literature
Project description:Cells communicate with each other via receptor-ligand interactions. Here, we describe lentiviral-mediated cell entry by engineered receptor-ligand interaction (ENTER) to display ligand proteins, deliver payloads, and record receptor specificity. We optimize ENTER to decode interactions between T cell receptor (TCR)-MHC peptides, antibody-antigen, and other receptor-ligand pairs. A viral presentation strategy allows ENTER to capture interactions between B cell receptor and any antigen. We engineer ENTER to deliver genetic payloads to antigen-specific T or B cells to selectively modulate cellular behavior in mixed populations. Single-cell readout of ENTER by RNA sequencing (ENTER-seq) enables multiplexed enumeration of antigen specificities, TCR clonality, cell type, and states of individual T cells. ENTER-seq of CMV-seropositive patient blood samples reveals the viral epitopes that drive effector memory T cell differentiation and inter-clonal vs. intra-clonal phenotypic diversity targeting the same epitope. ENTER technology enables systematic discovery of receptor specificity, linkage to cell fates, and antigen-specific cargo delivery.
| S-EPMC9789208 | biostudies-literature
Project description:The cellular abundance of proteins can vary even between isogenic single cells. This variability between single-cell protein levels can have regulatory roles, such as controlling cell fate during apoptosis induction or the proliferation/quiescence decision. Here, we review examples connecting protein levels and their dynamics in single cells to cellular functions. Such findings were made possible by the introduction of antibodies, and subsequently fluorescent proteins, for tracking protein levels in single cells. However, in heterogeneous cell populations, such as tumors or differentiating stem cells, cellular decisions are controlled by hundreds, even thousands of proteins acting in concert. Characterizing such complex systems demands measurements of thousands of proteins across thousands of single cells. This demand has inspired the development of new methods for single-cell protein analysis, and we discuss their trade-offs, with an emphasis on their specificity and coverage. We finish by highlighting the potential of emerging mass-spec methods to enable systems-level measurement of single-cell proteomes with unprecedented coverage and specificity. Combining such methods with methods for quantitating the transcriptomes and metabolomes of single cells will provide essential data for advancing quantitative systems biology.
| S-EPMC6204083 | biostudies-literature
Project description:Gene expression is controlled by interactions between trans-regulatory factors and cis-regulatory DNA sequences, and these interactions constitute the essential functional linkages of gene regulatory networks (GRNs). Validation of GRN models requires experimental cis-regulatory tests of predicted linkages to authenticate their identities and proposed functions. However, cis-regulatory analysis is, at present, at a severe bottleneck in genomic system biology because of the demanding experimental methodologies currently in use for discovering cis-regulatory modules (CRMs), in the genome, and for measuring their activities. Here we demonstrate a high-throughput approach to both discovery and quantitative characterization of CRMs. The unique aspect is use of DNA sequence tags to "barcode" CRM expression constructs, which can then be mixed, injected together into sea urchin eggs, and subsequently deconvolved. This method has increased the rate of cis-regulatory analysis by >100-fold compared with conventional one-by-one reporter assays. The utility of the DNA-tag reporters was demonstrated by the rapid discovery of 81 active CRMs from 37 previously unexplored sea urchin genes. We then obtained simultaneous high-resolution temporal characterization of the regulatory activities of more than 80 CRMs. On average 2-3 CRMs were discovered per gene. Comparison of endogenous gene expression profiles with those of the CRMs recovered from each gene showed that, for most cases, at least one CRM is active in each phase of endogenous expression, suggesting that CRM recovery was comprehensive. This approach will qualitatively alter the practice of GRN construction as well as validation, and will impact many additional areas of regulatory system biology.
| S-EPMC2840491 | biostudies-literature
Project description:Zika virus (ZIKV) is responsible for an ongoing and intensifying epidemic in the Western Hemisphere. We examined the complete predicted proteomes, glycomes, and selectomes of 33 ZIKV strains representing temporally diverse members of the African lineage, the Asian lineage, and the current outbreak in the Americas. Derivation of the complete selectome is an 'omics' approach to identify distinct evolutionary pressures acting on different features of an organism. Employment of the M8 model did not show evidence of global diversifying selection acting on the ZIKV polyprotein; however, a mixed effect model of evolution showed strong evidence (P<0.05) for episodic diversifying selection acting on specific sites. Single nucleotide polymorphisms (SNPs) were predictably frequent across strains relative to the derived consensus sequence. None of the 9 published detection procedures utilize targets that share 100% identity across the 33 strains examined, indicating that ZIKV escape from molecular detection is predictable. The predicted O-linked glycome showed marked diversity across strains; however, the N-linked glycome was highly stable. All Asian and American strains examined were predicted to include glycosylation of E protein ASN154, a modification proposed to mediate neurotropism, whereas the modification was not predicted for African strains. SNP diversity, episodic diversifying selection, and differential glycosylation, particularly of ASN154, may have major biological implications for ZIKV disease. Taken together, the systems biology perspective of ZIKV indicates: a.) The recently emergent Asian/American N-glycotype is mediating the new and emerging neuropathogenic potential of ZIKV; and b.) further divergence at specific sites is predictable as endemnicity is established in the Americas.
| S-EPMC5008700 | biostudies-literature
Project description:One motivation of systems biology research is to understand gene functions and interactions from functional genomics data such as that derived from microarrays. Up-to-date structural and functional annotations of genes are an essential foundation of systems biology modeling. We propose that the first essential step in any systems biology modeling of functional genomics data, especially for species with recently sequenced genomes, is gene structural and functional re-annotation. To demonstrate the impact of such re-annotation, we structurally and functionally re-annotated a microarray developed, and previously used, as a tool for disease research. We quantified the impact of this re-annotation on the array based on the total numbers of structural- and functional-annotations, the Gene Annotation Quality (GAQ) score, and canonical pathway coverage. We next quantified the impact of re-annotation on systems biology modeling using a previously published experiment that used this microarray. We show that re-annotation improves the quantity and quality of structural- and functional-annotations, allows a more comprehensive Gene Ontology based modeling, and improves pathway coverage for both the whole array and a differentially expressed mRNA subset. Our results also demonstrate that re-annotation can result in a different knowledge outcome derived from previous published research findings. We propose that, because of this, re-annotation should be considered to be an essential first step for deriving value from functional genomics data.
| S-EPMC2871057 | biostudies-literature
Project description:Mathematical models are routinely calibrated to experimental data, with goals ranging from building predictive models to quantifying parameters that cannot be measured. Whether or not reliable parameter estimates are obtainable from the available data can easily be overlooked. Such issues of parameter identifiability have important ramifications for both the predictive power of a model, and the mechanistic insight that can be obtained. Identifiability analysis is well-established for deterministic, ordinary differential equation (ODE) models, but there are no commonly adopted methods for analysing identifiability in stochastic models. We provide an accessible introduction to identifiability analysis and demonstrate how existing ideas for analysis of ODE models can be applied to stochastic differential equation (SDE) models through four practical case studies. To assess structural identifiability, we study ODEs that describe the statistical moments of the stochastic process using open-source software tools. Using practically motivated synthetic data and Markov chain Monte Carlo methods, we assess parameter identifiability in the context of available data. Our analysis shows that SDE models can often extract more information about parameters than deterministic descriptions. All code used to perform the analysis is available on Github.
| S-EPMC7811582 | biostudies-literature
Project description:Humans are not autonomous entities. We are all living in a complex environment, interacting not only with our peers, but as true holobionts; we are also very much in interaction with our coexisting microbial ecosystems living on and especially within us, in the intestine. Intestinal microorganisms, often collectively referred to as intestinal microbiota, contribute significantly to our daily energy uptake by breaking down complex carbohydrates into simple sugars, which are fermented to short-chain fatty acids and subsequently absorbed by human cells. They also have an impact on our immune system, by suppressing or enhancing the growth of malevolent and beneficial microbes. Our lifestyle can have a large influence on this ecosystem. What and how much we consume can tip the ecological balance in the intestine. A "western diet" containing mainly processed food will have a different effect on our health than a balanced diet fortified with pre- and probiotics. In recent years, new technologies have emerged, which made a more detailed understanding of microbial communities and ecosystems feasible. This includes progress in the sequencing of PCR-amplified phylogenetic marker genes as well as the collective microbial metagenome and metatranscriptome, allowing us to determine with an increasing level of detail, which microbial species are in the microbiota, understand what these microorganisms do and how they respond to changes in lifestyle and diet. These new technologies also include the use of synthetic and in vitro systems, which allow us to study the impact of substrates and addition of specific microbes to microbial communities at a high level of detail, and enable us to gather quantitative data for modelling purposes. Here, we will review the current state of microbiome research, summarizing the computational methodologies in this area and highlighting possible outcomes for personalized nutrition and medicine.
| S-EPMC5840735 | biostudies-literature
Project description:We describe a workflow for preprocessing a wide variety of single-cell genomics data types. The approach is based on parsing of machine-readable seqspec assay specifications to customize inputs for kb-python, which uses kallisto and bustools to catalog reads, error correct barcodes, and count reads. The universal preprocessing method is implemented in the Python package cellatlas that is available for download at: https://github.com/cellatlas/cellatlas/.
| S-EPMC10515959 | biostudies-literature