A prebiotic diet modulates microglial states and motor deficits in α-synuclein overexpressing mice.

Parkinson's disease (PD) is a movement disorder characterized by neuroinflammation, α-synuclein pathology, and neurodegeneration. Most cases of PD are non-hereditary, suggesting a strong role for environmental factors, and it has been speculated that disease may originate in peripheral tissues such as the gastrointestinal (GI) tract before affecting the brain. The gut microbiome is altered in PD and may impact motor and GI symptoms as indicated by animal studies, although mechanisms of gut-brain interactions remain incompletely defined. Intestinal bacteria ferment dietary fibers into short-chain fatty acids, with fecal levels of these molecules differing between PD and healthy controls and in mouse models. Among other effects, dietary microbial metabolites can modulate activation of microglia, brain-resident immune cells implicated in PD. We therefore investigated whether a fiber-rich diet influences microglial function in α-synuclein overexpressing (ASO) mice, a preclinical model with PD-like symptoms and pathology. Feeding a prebiotic high-fiber diet attenuates motor deficits and reduces α-synuclein aggregation in the substantia nigra of mice. Concomitantly, the gut microbiome of ASO mice adopts a profile correlated with health upon prebiotic treatment, which also reduces microglial activation. Single-cell RNA-seq analysis of microglia from the substantia nigra and striatum uncovers increased pro-inflammatory signaling and reduced homeostatic responses in ASO mice compared to wild-type counterparts on standard diets. However, prebiotic feeding reverses pathogenic microglial states in ASO mice and promotes expansion of protective disease-associated macrophage (DAM) subsets of microglia. Notably, depletion of microglia using a CSF1R inhibitor eliminates the beneficial effects of prebiotics by restoring motor deficits to ASO mice despite feeding a prebiotic diet. These studies uncover a novel microglia-dependent interaction between diet and motor symptoms in mice, findings that may have implications for neuroinflammation and PD.

Mapping hormone-regulated cell-cell interaction networks in the human breast at single-cell resolution.

The rise and fall of estrogen and progesterone across menstrual cycles and during pregnancy regulates breast development and modifies cancer risk. How these hormones impact each cell type in the breast remains poorly understood because they act indirectly through paracrine networks. Using single-cell analysis of premenopausal breast tissue, we reveal a network of coordinated transcriptional programs representing the tissue-level response to changing hormone levels. Our computational approach, DECIPHER-seq, leverages person-to-person variability in breast composition and cell state to uncover programs that co-vary across individuals. We use differences in cell-type proportions to infer a subset of programs that arise from direct cell-cell interactions regulated by hormones. Further, we demonstrate that prior pregnancy and obesity modify hormone responsiveness through distinct mechanisms: obesity reduces the proportion of hormone-responsive cells, whereas pregnancy dampens the direct response of these cells to hormones. Together, these results provide a comprehensive map of the cycling human breast.

CloudPred: Predicting Patient Phenotypes From Single-cell RNA-seq.

Single-cell RNA sequencing (scRNA-seq) has the potential to provide powerful, high-resolution signatures to inform disease prognosis and precision medicine. This paper takes an important first step towards this goal by developing an interpretable machine learning algorithm, CloudPred, to predict individuals' disease phenotypes from their scRNA-seq data. Predicting phenotype from scRNA-seq is challenging for standard machine learning methods-the number of cells measured can vary by orders of magnitude across individuals and the cell populations are also highly heterogeneous. Typical analysis creates pseudo-bulk samples which are biased toward prior annotations and also lose the single cell resolution. CloudPred addresses these challenges via a novel end-to-end differentiable learning algorithm which is coupled with a biologically informed mixture of cell types model. CloudPred automatically infers the cell subpopulation that are salient for the phenotype without prior annotations. We developed a systematic simulation platform to evaluate the performance of CloudPred and several alternative methods we propose, and find that CloudPred outperforms the alternative methods across several settings. We further validated CloudPred on a real scRNA-seq dataset of 142 lupus patients and controls. CloudPred achieves AUROC of 0.98 while identifying a specific subpopulation of CD4 T cells whose presence is highly indicative of lupus. CloudPred is a powerful new framework to predict clinical phenotypes from scRNA-seq data and to identify relevant cells.

Bimodal function of chromatin remodeler Hmga1 in neural crest induction and Wnt-dependent emigration.

During gastrulation, neural crest cells are specified at the neural plate border, as characterized by Pax7 expression. Using single-cell RNA sequencing coupled with high-resolution in situ hybridization to identify novel transcriptional regulators, we show that chromatin remodeler Hmga1 is highly expressed prior to specification and maintained in migrating chick neural crest cells. Temporally controlled CRISPR-Cas9-mediated knockouts uncovered two distinct functions of Hmga1 in neural crest development. At the neural plate border, Hmga1 regulates Pax7-dependent neural crest lineage specification. At premigratory stages, a second role manifests where Hmga1 loss reduces cranial crest emigration from the dorsal neural tube independent of Pax7. Interestingly, this is rescued by stabilized ß-catenin, thus implicating Hmga1 as a canonical Wnt activator. Together, our results show that Hmga1 functions in a bimodal manner during neural crest development to regulate specification at the neural plate border, and subsequent emigration from the neural tube via canonical Wnt signaling.

The neural plate is a structure that serves as the basis for the brain and central nervous system during the development of animals with a backbone. In particular, the tissues at the border of the neural plate become the neural crest, a group of highly mobile cells that can specialize to form nerves and parts of the face. The exact molecular mechanisms that allow the crest to emerge are still unknown. The protein Hmga1 alters how genes are packaged and organized inside cells, which in turn influences how genes are switched on and off. Here, Gandhi et al. studied how Hmga1 helps to shape the neural crest in developing chicken embryos. To do so, they harnessed a genetic tool called CRISPR-Cas9, and deleted the gene that encodes Hmga1 at specific developmental stages. This manipulation highlighted two periods where Hmga1 is active. First, Hmga1 helped to define neural crest cells at the neural plate border by activating a gene called pax7. Then, at a later stage, Hmga1 allowed these cells to move to other parts of the body by triggering the Wnt communication system. Failure for the neural crest to develop properly causes birth defects and cancers such as melanoma and childhood neuroblastoma, highlighting the need to better understand how this structure is formed. In addition, a better grasp of the roles of Hmga1 in healthy development could help to appreciate how it participates in a range of adult cancers.

Highly multiplexed single-cell RNA-seq by DNA oligonucleotide tagging of cellular proteins.

We describe a universal sample multiplexing method for single-cell RNA sequencing in which fixed cells are chemically labeled by attaching identifying DNA oligonucleotides to cellular proteins. Analysis of a 96-plex perturbation experiment revealed changes in cell population structure and transcriptional states that cannot be discerned from bulk measurements, establishing an efficient method for surveying cell populations from large experiments or clinical samples with the depth and resolution of single-cell RNA sequencing.

Integrated measurement of intracellular proteins and transcripts in single cells.

Biological function arises from the interplay of proteins, transcripts, and metabolites. An ongoing revolution in miniaturization technologies has created tools to analyze any one of these species in single cells, thus resolving the heterogeneity of tissues previously invisible to bulk measurements. An emerging frontier is single cell multi-omics, which is the measurement of multiple classes of analytes from single cells. Here, we combine bead-based transcriptomics with microchip-based proteomics to measure intracellular proteins and transcripts from single cells and defined small numbers of cells. The transcripts and proteins are independently measured by sequencing and fluorescent immunoassays respectively, to preserve their optimal measurement modes, and linked by encoding the physical address locations of the cells into digital sequencing space using spatially patterned DNA barcodes. We resolve cell-type-specific protein and transcript signatures and present a path forward to scaling the platform to high-throughput.

Adult Neurogenesis Is Sustained by Symmetric Self-Renewal and Differentiation.

Somatic stem cells have been identified in multiple adult tissues. Whether self-renewal occurs symmetrically or asymmetrically is key to understanding long-term stem cell maintenance and generation of progeny for cell replacement. In the adult mouse brain, neural stem cells (NSCs) (B1 cells) are retained in the walls of the lateral ventricles (ventricular-subventricular zone [V-SVZ]). The mechanism of B1 cell retention into adulthood for lifelong neurogenesis is unknown. Using multiple clonal labeling techniques, we show that the vast majority of B1 cells divide symmetrically. Whereas 20%-30% symmetrically self-renew and can remain in the niche for several months before generating neurons, 70%-80% undergo consuming divisions generating progeny, resulting in the depletion of B1 cells over time. This cellular mechanism decouples self-renewal from the generation of progeny. Limited rounds of symmetric self-renewal and consuming symmetric differentiation divisions can explain the levels of neurogenesis observed throughout life.

Diffusion as a Ruler: Modeling Kinesin Diffusion as a Length Sensor for Intraflagellar Transport.

An important question in cell biology is whether cells are able to measure size, either whole cell size or organelle size. Perhaps cells have an internal chemical representation of size that can be used to precisely regulate growth, or perhaps size is just an accident that emerges due to constraint of nutrients. The eukaryotic flagellum is an ideal model for studying size sensing and control because its linear geometry makes it essentially one-dimensional, greatly simplifying mathematical modeling. The assembly of flagella is regulated by intraflagellar transport (IFT), in which kinesin motors carry cargo adaptors for flagellar proteins along the flagellum and then deposit them at the tip, lengthening the flagellum. The rate at which IFT motors are recruited to begin transport into the flagellum is anticorrelated with the flagellar length, implying some kind of communication between the base and the tip and possibly indicating that cells contain some mechanism for measuring flagellar length. Although it is possible to imagine many complex scenarios in which additional signaling molecules sense length and carry feedback signals to the cell body to control IFT, might the already-known components of the IFT system be sufficient to allow length dependence of IFT? Here we investigate a model in which the anterograde kinesin motors unbind after cargo delivery, diffuse back to the base, and are subsequently reused to power entry of new IFT trains into the flagellum. By mathematically modeling and simulating such a system, we are able to show that the diffusion time of the motors can in principle be sufficient to serve as a proxy for length measurement. We found that the diffusion model can not only achieve a stable steady-state length without the addition of any other signaling molecules or pathways, but also is able to produce the anticorrelation between length and IFT recruitment rate that has been observed in quantitative imaging studies.

Transient Thresholding: A Mechanism Enabling Noncooperative Transcriptional Circuitry to Form a Switch.

Threshold generation in fate-selection circuits is often achieved through deterministic bistability, which requires cooperativity (i.e., nonlinear activation) and associated hysteresis. However, the Tat positive-feedback loop that controls HIV's fate decision between replication and proviral latency lacks self-cooperativity and deterministic bistability. Absent cooperativity, it is unclear how HIV can temporarily remain in an off-state long enough for the kinetically slower epigenetic silencing mechanisms to act-expression fluctuations should rapidly trigger active positive feedback and replication, precluding establishment of latency. Here, using flow cytometry and single-cell imaging, we find that the Tat circuit exhibits a transient activation threshold. This threshold largely disappears after ∼40 h-accounting for the lack of deterministic bistability-and promoter activation shortens the lifetime of this transient threshold. Continuous differential equation models do not recapitulate this phenomenon. However, chemical reaction (master equation) models where the transcriptional transactivator and promoter toggle between inactive and active states can recapitulate the phenomenon because they intrinsically create a single-molecule threshold transiently requiring excess molecules in the inactive state to achieve at least one molecule (rather than a continuous fractional value) in the active state. Given the widespread nature of promoter toggling and transcription factor modifications, transient thresholds may be a general feature of inducible promoters.

In vitro patterning of pluripotent stem cell-derived intestine recapitulates in vivo human development.

The intestine plays a central role in digestion, nutrient absorption and metabolism, with individual regions of the intestine having distinct functional roles. Many examples of region-specific gene expression in the adult intestine are known, but how intestinal regional identity is established during development is a largely unresolved issue. Here, we have identified several genes that are expressed in a region-specific manner in the developing human intestine. Using human embryonic stem cell-derived intestinal organoids, we demonstrate that the duration of exposure to active FGF and WNT signaling controls regional identity. Short-term exposure to FGF4 and CHIR99021 (a GSK3ß inhibitor that stabilizes ß-catenin) resulted in organoids with gene expression patterns similar to developing human duodenum, whereas longer exposure resulted in organoids similar to ileum. When region-specific organoids were transplanted into immunocompromised mice, duodenum-like organoids and ileum-like organoids retained their regional identity, demonstrating that regional identity of organoids is stable after initial patterning occurs. This work provides insights into the mechanisms that control regional specification of the developing human intestine and provides new tools for basic and translational research.

Signaling Boundary Conditions Drive Self-Organization of Human "Gastruloids".

Morphogenesis requires tissues to sense and respond to their geometry. In this issue of Developmental Cell, Etoc et al. (2016) show that a confined colony of human embryonic stem cells can spontaneously sense its boundary, generating a self-organized TGF-ß signaling gradient that patterns it into a tissue resembling a gastrulating embryo.

SOX2 O-GlcNAcylation alters its protein-protein interactions and genomic occupancy to modulate gene expression in pluripotent cells.

The transcription factor SOX2 is central in establishing and maintaining pluripotency. The processes that modulate SOX2 activity to promote pluripotency are not well understood. Here, we show SOX2 is O-GlcNAc modified in its transactivation domain during reprogramming and in mouse embryonic stem cells (mESCs). Upon induction of differentiation SOX2 O-GlcNAcylation at serine 248 is decreased. Replacing wild type with an O-GlcNAc-deficient SOX2 (S248A) increases reprogramming efficiency. ESCs with O-GlcNAc-deficient SOX2 exhibit alterations in gene expression. This change correlates with altered protein-protein interactions and genomic occupancy of the O-GlcNAc-deficient SOX2 compared to wild type. In addition, SOX2 O-GlcNAcylation impairs the SOX2-PARP1 interaction, which has been shown to regulate ESC self-renewal. These findings show that SOX2 activity is modulated by O-GlcNAc, and provide a novel regulatory mechanism for this crucial pluripotency transcription factor.

Engineering Customized Cell Sensing and Response Behaviors Using Synthetic Notch Receptors.

The Notch protein is one of the most mechanistically direct transmembrane receptors-the intracellular domain contains a transcriptional regulator that is released from the membrane when engagement of the cognate extracellular ligand induces intramembrane proteolysis. We find that chimeric forms of Notch, in which both the extracellular sensor module and the intracellular transcriptional module are replaced with heterologous protein domains, can serve as a general platform for generating novel cell-cell contact signaling pathways. Synthetic Notch (synNotch) pathways can drive user-defined functional responses in diverse mammalian cell types. Because individual synNotch pathways do not share common signaling intermediates, the pathways are functionally orthogonal. Thus, multiple synNotch receptors can be used in the same cell to achieve combinatorial integration of environmental cues, including Boolean response programs, multi-cellular signaling cascades, and self-organized cellular patterns. SynNotch receptors provide extraordinary flexibility in engineering cells with customized sensing/response behaviors to user-specified extracellular cues.

Transcriptome-wide Analysis Reveals Hallmarks of Human Intestine Development and Maturation In Vitro and In Vivo.

Human intestinal organoids (HIOs) are a tissue culture model in which small intestine-like tissue is generated from pluripotent stem cells. By carrying out unsupervised hierarchical clustering of RNA-sequencing data, we demonstrate that HIOs most closely resemble human fetal intestine. We observed that genes involved in digestive tract development are enriched in both fetal intestine and HIOs compared to adult tissue, whereas genes related to digestive function and Paneth cell host defense are expressed at higher levels in adult intestine. Our study also revealed that the intestinal stem cell marker OLFM4 is expressed at very low levels in fetal intestine and in HIOs, but is robust in adult crypts. We validated our findings using in vivo transplantation to show that HIOs become more adult-like after transplantation. Our study emphasizes important maturation events that occur in the intestine during human development and demonstrates that HIOs can be used to model fetal-to-adult maturation.

A strategy for tissue self-organization that is robust to cellular heterogeneity and plasticity.

Developing tissues contain motile populations of cells that can self-organize into spatially ordered tissues based on differences in their interfacial surface energies. However, it is unclear how self-organization by this mechanism remains robust when interfacial energies become heterogeneous in either time or space. The ducts and acini of the human mammary gland are prototypical heterogeneous and dynamic tissues comprising two concentrically arranged cell types. To investigate the consequences of cellular heterogeneity and plasticity on cell positioning in the mammary gland, we reconstituted its self-organization from aggregates of primary cells in vitro. We find that self-organization is dominated by the interfacial energy of the tissue-ECM boundary, rather than by differential homo- and heterotypic energies of cell-cell interaction. Surprisingly, interactions with the tissue-ECM boundary are binary, in that only one cell type interacts appreciably with the boundary. Using mathematical modeling and cell-type-specific knockdown of key regulators of cell-cell cohesion, we show that this strategy of self-organization is robust to severe perturbations affecting cell-cell contact formation. We also find that this mechanism of self-organization is conserved in the human prostate. Therefore, a binary interfacial interaction with the tissue boundary provides a flexible and generalizable strategy for forming and maintaining the structure of two-component tissues that exhibit abundant heterogeneity and plasticity. Our model also predicts that mutations affecting binary cell-ECM interactions are catastrophic and could contribute to loss of tissue architecture in diseases such as breast cancer.

The rational parameterization theorem for multisite post-translational modification systems.

Post-translational modification of proteins plays a central role in cellular regulation but its study has been hampered by the exponential increase in substrate modification forms ("modforms") with increasing numbers of sites. We consider here biochemical networks arising from post-translational modification under mass-action kinetics, allowing for multiple substrates, having different types of modification (phosphorylation, methylation, acetylation, etc.) on multiple sites, acted upon by multiple forward and reverse enzymes (in total number L), using general enzymatic mechanisms. These assumptions are substantially more general than in previous studies. We show that the steady-state modform concentrations constitute an algebraic variety that can be parameterized by rational functions of the L free enzyme concentrations, with coefficients which are rational functions of the rate constants. The parameterization allows steady states to be calculated by solving L algebraic equations, a dramatic reduction compared to simulating an exponentially large number of differential equations. This complexity collapse enables analysis in contexts that were previously intractable and leads to biological predictions that we review. Our results lay a foundation for the systems biology of post-translational modification and suggest deeper connections between biochemical networks and algebraic geometry.

Unlimited multistability in multisite phosphorylation systems.

Reversible phosphorylation on serine, threonine and tyrosine is the most widely studied posttranslational modification of proteins. The number of phosphorylated sites on a protein (n) shows a significant increase from prokaryotes, with n /= 150 sites. Multisite phosphorylation has many roles and site conservation indicates that increasing numbers of sites cannot be due merely to promiscuous phosphorylation. A substrate with n sites has an exponential number (2(n)) of phospho-forms and individual phospho-forms may have distinct biological effects. The distribution of these phospho-forms and how this distribution is regulated have remained unknown. Here we show that, when kinase and phosphatase act in opposition on a multisite substrate, the system can exhibit distinct stable phospho-form distributions at steady state and that the maximum number of such distributions increases with n. Whereas some stable distributions are focused on a single phospho-form, others are more diffuse, giving the phospho-proteome the potential to behave as a fluid regulatory network able to encode information and flexibly respond to varying demands. Such plasticity may underlie complex information processing in eukaryotic cells and suggests a functional advantage in having many sites. Our results follow from the unusual geometry of the steady-state phospho-form concentrations, which we show to constitute a rational algebraic curve, irrespective of n. We thereby reduce the complexity of calculating steady states from simulating 3 x 2(n) differential equations to solving two algebraic equations, while treating parameters symbolically. We anticipate that these methods can be extended to systems with multiple substrates and multiple enzymes catalysing different modifications, as found in posttranslational modification 'codes' such as the histone code. Whereas simulations struggle with exponentially increasing molecular complexity, mathematical methods of the kind developed here can provide a new language in which to articulate the principles of cellular information processing.

Programming with models: modularity and abstraction provide powerful capabilities for systems biology.

Mathematical models are increasingly used to understand how phenotypes emerge from systems of molecular interactions. However, their current construction as monolithic sets of equations presents a fundamental barrier to progress. Overcoming this requires modularity, enabling sub-systems to be specified independently and combined incrementally, and abstraction, enabling generic properties of biological processes to be specified independently of specific instances. These, in turn, require models to be represented as programs rather than as datatypes. Programmable modularity and abstraction enables libraries of modules to be created, which can be instantiated and reused repeatedly in different contexts with different components. We have developed a computational infrastructure that accomplishes this. We show here why such capabilities are needed, what is required to implement them and what can be accomplished with them that could not be done previously.

Influence of common orthodontic appliances on the diagnostic quality of cranial magnetic resonance images.

INTRODUCTION: The aim of this study was to evaluate cranial magnetic resonance (MR) image distortion caused by various orthodontic brackets. METHODS: Ten subjects received 5 consecutive cranial MR scans. A control scan was conducted with Essix trays (GAC International, Bohemia, NY) fitted over the maxillary and mandibular teeth. Four experimental MR scans of the head were conducted with plastic, ceramic, titanium, and stainless steel brackets incorporated into the Essix tray material. Each MR scan consisted of 4 sequences: sagittal T1-weighted spin echo (T1 sagittal), axial T2-weighted spin echo (T2 axial), gradient echo, and diffusion-weighted imaging. Three board-certified neuroradiologists examined the MR images for distortion in predetermined regions of the head. RESULTS: The paired Wilcoxon signed rank test showed a statistically significant difference between the mean distortion scores of stainless steel brackets and the mean distortion scores of the other experimental MR scans (P <0.0001). Interrater and intrarater agreement was high (kappa statistic and associated 95% confidence intervals). CONCLUSIONS: The study showed that plastic, ceramic, and titanium brackets cause minimal distortion of cranial MR images (similar to the control). On the other hand, stainless steel brackets cause significant distortion, rendering several cranial regions nondiagnostic. Areas with the most distortion were the body of the mandible, the hard palate, the base of the tongue, the globes, the nasopharynx, and the frontal lobes. In general, the closer the stainless steel appliance was to a specific anatomic region, the greater the distortion of the MR image.