CRISPR-Based Editing Reveals Edge-Specific Effects in Biological Networks

Modeling Control of Supercoiling Dynamics and Transcription Using DNA-Binding Proteins

Nearly all natural and synthetic gene networks rely on the fundamental process of transcription to enact biological feedback, genetic programs, and living circuitry. In this work, we investigate the efficacy of controlling transcription using a new biophysical mechanism, control of localized supercoiling near a gene of interest. We postulate a basic reaction network model for describing the general phenomenon of transcription and introduce a separate set of equations to describe the dynamics of supercoiling. We show that supercoiling and transcription introduce a shared reaction flux term in the model dynamics and illustrate how the modulation of supercoiling can be used to control transcription rates. We show the supercoiling-transcription model can be written as a nonlinear state-space model, with a radial basis function nonlinearity to capture the empirical relationship between supercoiling and transcription rates. We show the system admits a single, globally exponentially stable equilibrium point. Notably, we show that mRNA steady-state levels can be controlled directly by increasing a length-scale parameter for genetic spacing. Finally, we build a mathematical model to explore the use of a DNA binding protein, to define programmable boundary conditions on supercoiling propagation, which we show can be used to control transcriptional bursting or pulsatile transcriptional response. We show there exists a stabilizing control law for mRNA tracking, using the method of control Lyapunov functions and illustrate these results with numerical simulations.

Reduction of Tumor Growth with RNA-Targeting Treatment of the NAB2-STAT6 Fusion Transcript in Solitary Fibrous Tumor Models

Solitary fibrous tumor (SFT) is a rare soft-tissue sarcoma. This nonhereditary cancer is the result of an environmental intrachromosomal gene fusion between NAB2 and STAT6 on chromosome 12, which fuses the activation domain of STAT6 with the repression domain of NAB2. Currently there is not an approved chemotherapy regimen for SFTs. The best response on available pharmaceuticals is a partial response or stable disease for several months. The purpose of this study is to investigate the potential of RNA-based therapies for the treatment of SFTs. Specifically, in vitro SFT cell models were engineered to harbor the characteristic NAB2-STAT6 fusion using the CRISPR/SpCas9 system. Cell migration as well as multiple cancer-related signaling pathways were increased in the engineered cells as compared to the fusion-absent parent cells. The SFT cell models were then used for evaluating the targeting efficacies of NAB2-STAT6 fusion-specific antisense oligonucleotides (ASOs) and CRISPR/CasRx systems. Our results showed that fusion specific ASO treatments caused a 58% reduction in expression of fusion transcripts and a 22% reduction in cell proliferation after 72 h in vitro. Similarly, the AAV2-mediated CRISPR/CasRx system led to a 59% reduction in fusion transcript expressions in vitro, and a 55% reduction in xenograft growth after 29 days ex vivo.

Genetic physical unclonable functions in human cells

A physical unclonable function (PUF) is a physical entity that provides a measurable output that can be used as a unique and irreproducible identifier for the artifact wherein it is embedded. Popularized by the electronics industry, silicon PUFs leverage the inherent physical variations of semiconductor manufacturing to establish intrinsic security primitives for attesting integrated circuits. Owing to the stochastic nature of these variations, photolithographically manufactured silicon PUFs are impossible to reproduce (thus unclonable). Inspired by the success of silicon PUFs, we sought to create the first generation of genetic PUFs in human cells. We demonstrate that these PUFs are robust (i.e., they repeatedly produce the same output), unique (i.e., they do not coincide with any other identically produced PUF), and unclonable (i.e., they are virtually impossible to replicate). Furthermore, we demonstrate that CRISPR-engineered PUFs (CRISPR-PUFs) can serve as a foundational principle for establishing provenance attestation protocols.

Machine learning-based approaches for identifying human blood cells harboring CRISPR-mediated fetal chromatin domain ablations

Two common hemoglobinopathies, sickle cell disease (SCD) and β-thalassemia, arise from genetic mutations within the β-globin gene. In this work, we identified a 500-bp motif (Fetal Chromatin Domain, FCD) upstream of human ϒ-globin locus and showed that the removal of this motif using CRISPR technology reactivates the expression of ϒ-globin. Next, we present two different cell morphology-based machine learning approaches that can be used identify human blood cells (KU-812) that harbor CRISPR-mediated FCD genetic modifications. Three candidate models from the first approach, which uses multilayer perceptron algorithm (MLP 20-26, MLP26-18, and MLP 30-26) and flow cytometry-derived cellular data, yielded 0.83 precision, 0.80 recall, 0.82 accuracy, and 0.90 area under the ROC (receiver operating characteristic) curve when predicting the edited cells. In comparison, the candidate model from the second approach, which uses deep learning (T2D5) and DIC microscopy-derived imaging data, performed with less accuracy (0.80) and ROC AUC (0.87). We envision that equivalent machine learning-based models can complement currently available genotyping protocols for specific genetic modifications which result in morphological changes in human cells.

Cell morphology-based machine learning models for human cell state classification

Herein, we implement and access machine learning architectures to ascertain models that differentiate healthy from apoptotic cells using exclusively forward (FSC) and side (SSC) scatter flow cytometry information. To generate training data, colorectal cancer HCT116 cells were subjected to miR-34a treatment and then classified using a conventional Annexin V/propidium iodide (PI)-staining assay. The apoptotic cells were defined as Annexin V-positive cells, which include early and late apoptotic cells, necrotic cells, as well as other dying or dead cells. In addition to fluorescent signal, we collected cell size and granularity information from the FSC and SSC parameters. Both parameters are subdivided into area, height, and width, thus providing a total of six numerical features that informed and trained our models. A collection of logistical regression, random forest, k-nearest neighbor, multilayer perceptron, and support vector machine was trained and tested for classification performance in predicting cell states using only the six aforementioned numerical features. Out of 1046 candidate models, a multilayer perceptron was chosen with 0.91 live precision, 0.93 live recall, 0.92 live f value and 0.97 live area under the ROC curve when applied on standardized data. We discuss and highlight differences in classifier performance and compare the results to the standard practice of forward and side scatter gating, typically performed to select cells based on size and/or complexity. We demonstrate that our model, a ready-to-use module for any flow cytometry-based analysis, can provide automated, reliable, and stain-free classification of healthy and apoptotic cells using exclusively size and granularity information.

Uncoupling gene expression noise along the central dogma using genome engineered human cell lines

Eukaryotic protein synthesis is an inherently stochastic process. This stochasticity stems not only from variations in cell content between cells but also from thermodynamic fluctuations in a single cell. Ultimately, these inherently stochastic processes manifest as noise in gene expression, where even genetically identical cells in the same environment exhibit variation in their protein abundances. In order to elucidate the underlying sources that contribute to gene expression noise, we quantify the contribution of each step within the process of protein synthesis along the central dogma. We uncouple gene expression at the transcriptional, translational, and post-translational level using custom engineered circuits stably integrated in human cells using CRISPR. We provide a generalized framework to approximate intrinsic and extrinsic noise in a population of cells expressing an unbalanced two-reporter system. Our decomposition shows that the majority of intrinsic fluctuations stem from transcription and that coupling the two genes along the central dogma forces the fluctuations to propagate and accumulate along the same path, resulting in increased observed global correlation between the products.

Molecular networks in Network Medicine: Development and applications

Network Medicine applies network science approaches to investigate disease pathogenesis. Many different analytical methods have been used to infer relevant molecular networks, including protein-protein interaction networks, correlation-based networks, gene regulatory networks, and Bayesian networks. Network Medicine applies these integrated approaches to Omics Big Data (including genetics, epigenetics, transcriptomics, metabolomics, and proteomics) using computational biology tools and, thereby, has the potential to provide improvements in the diagnosis, prognosis, and treatment of complex diseases. We discuss briefly the types of molecular data that are used in molecular network analyses, survey the analytical methods for inferring molecular networks, and review efforts to validate and visualize molecular networks. Successful applications of molecular network analysis have been reported in pulmonary arterial hypertension, coronary heart disease, diabetes mellitus, chronic lung diseases, and drug development. Important knowledge gaps in Network Medicine include incompleteness of the molecular interactome, challenges in identifying key genes within genetic association regions, and limited applications to human diseases. This article is categorized under: Models of Systems Properties and Processes > Mechanistic Models Translational, Genomic, and Systems Medicine > Translational Medicine Analytical and Computational Methods > Analytical Methods Analytical and Computational Methods > Computational Methods.

Coevolutionary Couplings Unravel PAM-Proximal Constraints of CRISPR-SpCas9

The clustered regularly interspaced short palindromic repeats (CRISPR) system, an immune system analog found in prokaryotes, allows a single-guide RNA to direct a CRISPR-associated protein (Cas) with combined helicase and nuclease activity to DNA. The presence of a specific protospacer adjacent motif (PAM) next to the DNA target site plays a crucial role in determining both efficacy and specificity of gene editing. Herein, we introduce a coevolutionary framework to computationally unveil nonobvious molecular interactions in CRISPR systems and experimentally probe their functional role. Specifically, we use direct coupling analysis, a statistical inference framework used to infer direct coevolutionary couplings, in the context of protein/nucleic acid interactions. Applied to Streptococcus pyogenes Cas9, a Hamiltonian metric obtained from coevolutionary relationships reveals, to our knowledge, novel PAM-proximal nucleotide preferences at the seventh position of S. pyogenes Cas9 PAM (5'-NGRNNNT-3'), which was experimentally confirmed by in vitro and functional assays in human cells. We show that coevolved and conserved interactions point to specific clues toward rationally engineering new generations of Cas9 systems and may eventually help decipher the diversity of this family of proteins.