IQCELL: A platform for predicting the effect of gene perturbations on developmental trajectories using single-cell RNA-seq data.

The increasing availability of single-cell RNA-sequencing (scRNA-seq) data from various developmental systems provides the opportunity to infer gene regulatory networks (GRNs) directly from data. Herein we describe IQCELL, a platform to infer, simulate, and study executable logical GRNs directly from scRNA-seq data. Such executable GRNs allow simulation of fundamental hypotheses governing developmental programs and help accelerate the design of strategies to control stem cell fate. We first describe the architecture of IQCELL. Next, we apply IQCELL to scRNA-seq datasets from early mouse T-cell and red blood cell development, and show that the platform can infer overall over 74% of causal gene interactions previously reported from decades of research. We will also show that dynamic simulations of the generated GRN qualitatively recapitulate the effects of known gene perturbations. Finally, we implement an IQCELL gene selection pipeline that allows us to identify candidate genes, without prior knowledge. We demonstrate that GRN simulations based on the inferred set yield results similar to the original curated lists. In summary, the IQCELL platform offers a versatile tool to infer, simulate, and study executable GRNs in dynamic biological systems.

How do the object-file and physical-reasoning systems interact? Evidence from priming effects with object arrays or novel labels.

How do infants reason about simple physical events such as containment, tube, and support events? According to the two-system model, two cognitive systems, the object-file (OF) and physical-reasoning (PR) systems, work together to guide infants' responses to these events. When an event begins, the OF system sends categorical information about the objects and their arrangements to the PR system. This system then categorizes the event, assigns event roles to the objects, and taps the OF system for information about features previously identified as causally relevant for the event category selected. All of the categorical and featural information included in the event's representation is interpreted by the PR system's domain knowledge, which includes core principles such as persistence and gravity. The present research tested a novel prediction of the model: If the OF system could be primed to also send, at the beginning of an event, information about an as-yet-unidentified feature, the PR system would then interpret this information using its core principles, allowing infants to detect core violations involving the feature earlier than they normally would. We examined this prediction using two types of priming manipulations directed at the OF system, object arrays and novel labels. In six experiments, infants aged 7-13 months (N = 304) were tested using different event categories and as-yet-unidentified features (color in containment events, height in tube events, and proportional distribution in support events) as well as different tasks (violation-of-expectation and action tasks). In each case, infants who were effectively primed reasoned successfully about the as-yet-unidentified feature, sometimes as early as six months before they would typically do so. These converging results provide strong support for the two-system model and for the claim that uncovering how the OF and PR systems represent and exchange information is essential for understanding how infants respond to physical events.

What's New to You? Preschoolers' Partner-Specific Online Processing of Disfluency.

Speech disfluencies (e.g., "Point to thee um turtle") can signal that a speaker is about to refer to something difficult to name. In two experiments, we found evidence that 4-year-olds, like adults, flexibly interpret a particular partner's disfluency based on their estimate of that partner's knowledge, derived from the preceding conversation. In entrainment trials, children established partner-specific shared knowledge of names for tangram pictures with one or two adult interlocutors. In each test trial, an adult named one of two visible tangrams either fluently or disfluently while children's eye-movements were monitored. We manipulated speaker knowledge in the test trials. In Experiment 1, the test-trial speaker was the same speaker from entrainment or a naïve experimenter; in Experiment 2, the test-trial speaker had been one of the child's partners in entrainment and had seen half of the tangrams (either animal or vehicle tangrams). When hearing disfluent expressions, children looked more at a tangram that was unfamiliar from the speaker's perspective; this systematic disfluency effect disappeared in Experiment 1 when the speaker was entirely naïve, and depended on each speaker's entrainment experience in Experiment 2. These findings show that 4-year-olds can keep track of two different partners' knowledge states, and use this information to determine what should be difficult for a particular partner to name, doing so efficiently enough to guide online interpretation of disfluent speech.

An efficient method for generation of bi-allelic null mutant mouse embryonic stem cells and its application for investigating epigenetic modifiers.

Mouse embryonic stem (ES) cells are a popular model system to study biological processes, though uncovering recessive phenotypes requires inactivating both alleles. Building upon resources from the International Knockout Mouse Consortium (IKMC), we developed a targeting vector for second allele inactivation in conditional-ready IKMC 'knockout-first' ES cell lines. We applied our technology to several epigenetic regulators, recovering bi-allelic targeted clones with a high efficiency of 60% and used Flp recombinase to restore expression in two null cell lines to demonstrate how our system confirms causality through mutant phenotype reversion. We designed our strategy to select against re-targeting the 'knockout-first' allele and identify essential genes in ES cells, including the histone methyltransferase Setdb1. For confirmation, we exploited the flexibility of our system, enabling tamoxifen inducible conditional gene ablation while controlling for genetic background and tamoxifen effects. Setdb1 ablated ES cells exhibit severe growth inhibition, which is not rescued by exogenous Nanog expression or culturing in naive pluripotency '2i' media, suggesting that the self-renewal defect is mediated through pluripotency network independent pathways. Our strategy to generate null mutant mouse ES cells is applicable to thousands of genes and repurposes existing IKMC Intermediate Vectors.

A microfluidic device for characterizing nuclear deformations.

Cell nuclei experience and respond to a wide range of forces, both in vivo and in vitro. In order to characterize the nuclear response to physical stress, we developed a microfluidic chip and used it to apply mechanical stress to live cells and measure their nuclear deformability. The device design is optimized for the detection of both nucleus and cytoplasm, which can then be conveniently quantified using a custom-written Matlab program. We measured nuclear sizes and strains of embryonic stem cells, for which we observed negative Poisson ratios in the nuclei. In addition, we were able to detect changes in the nuclear response after treatment with actin depolymerizing and chromatin decondensing agents. Finally, we showed that the device can be used for biologically relevant high-resolution confocal imaging of cells under compression. Thus, the device presented here allows for accurate physical phenotyping at high throughput and has the potential to be applied to a range of cell types.

The interaction of PRC2 with RNA or chromatin is mutually antagonistic.

Polycomb repressive complex 2 (PRC2) modifies chromatin to maintain genes in a repressed state during development. PRC2 is primarily associated with CpG islands at repressed genes and also possesses RNA binding activity. However, the RNAs that bind PRC2 in cells, the subunits that mediate these interactions, and the role of RNA in PRC2 recruitment to chromatin all remain unclear. By performing iCLIP for PRC2 in comparison with other RNA binding proteins, we show here that PRC2 binds nascent RNA at essentially all active genes. Although interacting with RNA promiscuously, PRC2 binding is enriched at specific locations within RNAs, primarily exon-intron boundaries and the 3' UTR. Deletion of other PRC2 subunits reveals that SUZ12 is sufficient to establish this RNA binding profile. Contrary to prevailing models, we also demonstrate that the interaction of PRC2 with RNA or chromatin is mutually antagonistic in cells and in vitro. RNA degradation in cells triggers PRC2 recruitment to CpG islands at active genes. Correspondingly, the release of PRC2 from chromatin in cells increases RNA binding. Consistent with this, RNA and nucleosomes compete for PRC2 binding in vitro. We propose that RNA prevents PRC2 recruitment to chromatin at active genes and that mutual antagonism between RNA and chromatin underlies the pattern of PRC2 chromatin association across the genome.

Auxetic nuclei in embryonic stem cells exiting pluripotency.

Embryonic stem cells (ESCs) self-renew in a state of naïve pluripotency in which they are competent to generate all somatic cells. It has been hypothesized that, before irreversibly committing, ESCs pass through at least one metastable transition state. This transition would represent a gateway for differentiation and reprogramming of somatic cells. Here, we show that during the transition, the nuclei of ESCs are auxetic: they exhibit a cross-sectional expansion when stretched and a cross-sectional contraction when compressed, and their stiffness increases under compression. We also show that the auxetic phenotype of transition ESC nuclei is driven at least in part by global chromatin decondensation. Through the regulation of molecular turnover in the differentiating nucleus by external forces, auxeticity could be a key element in mechanotransduction. Our findings highlight the importance of nuclear structure in the regulation of differentiation and reprogramming.

Chromatin states in pluripotent, differentiated, and reprogrammed cells.

The pluripotent state of embryonic stem cells is maintained by a core network of transcription factors, and by chromatin remodelling factors that support an environment permissive for transcription. Polycomb and trithorax Group proteins enable 'bivalent' chromatin to be established at lineage-specific genes within pluripotent cells that is thought to poise genes for rapid activation upon induction of differentiation. As differentiation proceeds, chromatin condenses and there is a genome-wide increase in the abundance of repressive histone modifications, alterations in the subnuclear organisation of particular genomic regions, and changes in DNA methylation profiles within genes. Reprogramming of somatic cells provides a platform to investigate the role of chromatin-based factors in establishing and maintaining pluripotency.

Loss-of-function Additional sex combs like 1 mutations disrupt hematopoiesis but do not cause severe myelodysplasia or leukemia.

The Additional sex combs like 1 (Asxl1) gene is 1 of 3 mammalian homologs of the Additional sex combs (Asx) gene of Drosophila. Asx is unusual because it is required to maintain both activation and silencing of Hox genes in flies and mice. Asxl proteins are characterized by an amino terminal homology domain, by interaction domains for nuclear receptors, and by a C-terminal plant homeodomain protein-protein interaction domain. A recent study of patients with myelodysplastic syndrome (MDS) and chronic myelomonocytic leukemia (CMML) revealed a high incidence of truncation mutations that would delete the PHD domain of ASXL1. Here, we show that Asxl1 is expressed in all hematopoietic cell fractions analyzed. Asxl1 knockout mice exhibit defects in frequency of differentiation of lymphoid and myeloid progenitors, but not in multipotent progenitors. We do not detect effects on hematopoietic stem cells, or in peripheral blood. Notably, we do not detect severe myelodysplastic phenotypes or leukemia in this loss-of-function model. We conclude that Asxl1 is needed for normal hematopoiesis. The mild phenotypes observed may be because other Asxl genes have redundant function with Asxl1, or alternatively, MDS or oncogenic phenotypes may result from gain-of-function Asxl mutations caused by genomic amplification, gene fusion, or truncation of Asxl1.

Characterization of Asxl1, a murine homolog of Additional sex combs, and analysis of the Asx-like gene family.

The Additional sex combs (Asx) gene of Drosophila is required to maintain homeotic gene activation and silencing. Here we characterize the three murine homologs of Asx: Additional sex combs-like (Asxl1, Asxl2, and Asxl3) and identify conserved sequence features. The predicted amino acid sequence of Asxl1 has 16% identity and 40% similarity to Drosophila Asx, and 74% identity and 81% similarity to human ASXL1. Murine Asxl1 contains two regions highly conserved with Drosophila Asx: 1) a conserved amino terminal region of unknown function, termed the ASX homology domain (ASXH) which contains two nuclear receptor (NR) co-regulator binding motifs; and 2) a conserved C-terminal PHD domain. The mammalian Asx-like predicted proteins possess additional conserved sequence features of unknown function not present in Drosophila Asx, including three more NR co-regulator binding motifs. Asxl1and Asxl2 are expressed as multiple transcripts, at varying levels, in adult tissues and in embryonic stem cells analyzed by Northern blot, and exhibit similar expression patterns suggesting they may be co-regulated. Whole mount RNA in situ hybridization revealed that Asxl1 is also expressed in 10.5-11.0 dpc mouse embryos.

Maintenance of gene expression patterns.

In development, cells pass on established gene expression patterns to daughter cells over multiple rounds of cell division. The cellular memory of the gene expression state is termed maintenance, and the proteins required for this process are termed maintenance proteins. The best characterized are proteins of the Polycomb and trithorax Groups that are required for silencing and maintenance of activation of target loci, respectively. These proteins act through DNA elements termed maintenance elements. Here, we re-examine the genetics and molecular biology of maintenance proteins. We discuss molecular models for the maintenance of activation and silencing, and the establishment of epigenetic marks, and suggest that maintenance proteins may play a role in propagating the mark through DNA synthesis.

A human homolog of Additional sex combs, ADDITIONAL SEX COMBS-LIKE 1, maps to chromosome 20q11.

The Additional sex combs (Asx) gene of Drosophila is an Enhancer of trithorax and Polycomb (ETP), and is required to maintain activation and silencing of homeotic loci. The molecular basis of this dual function is not understood. Here, we identify a human homolog of Asx, termed ADDITIONAL SEX COMBS-LIKE 1 (ASXL1). Overall, the amino acid sequence of ASXL1 has 21% identity and 41% similarity to Drosophila ASX. The ASXL1 protein contains a 118 amino acid conserved amino terminal region of unknown function that we term the ASX homology domain (ASXH) that contains two LXXLL consensus sequences for nuclear receptor binding. ASXL1 also contains a conserved C-terminal cysteine cluster that is a variant of the PHD domain. Three ASXL1 transcripts of differing size are detected, which are widely expressed in adult tissues. ASXL1 maps to chromosome 20q11, a region frequently amplified in human tumors. Interestingly, ASXL1 is overexpressed in cell lines derived from carcinomas.