Single-cell, whole-embryo phenotyping of mammalian developmental disorders.

Mouse models are a critical tool for studying human diseases, particularly developmental disorders1. However, conventional approaches for phenotyping may fail to detect subtle defects throughout the developing mouse2. Here we set out to establish single-cell RNA sequencing of the whole embryo as a scalable platform for the systematic phenotyping of mouse genetic models. We applied combinatorial indexing-based single-cell RNA sequencing3 to profile 101 embryos of 22 mutant and 4 wild-type genotypes at embryonic day 13.5, altogether profiling more than 1.6 million nuclei. The 22 mutants represent a range of anticipated phenotypic severities, from established multisystem disorders to deletions of individual regulatory regions4,5. We developed and applied several analytical frameworks for detecting differences in composition and/or gene expression across 52 cell types or trajectories. Some mutants exhibit changes in dozens of trajectories whereas others exhibit changes in only a few cell types. We also identify differences between widely used wild-type strains, compare phenotyping of gain- versus loss-of-function mutants and characterize deletions of topological associating domain boundaries. Notably, some changes are shared among mutants, suggesting that developmental pleiotropy might be 'decomposable' through further scaling of this approach. Overall, our findings show how single-cell profiling of whole embryos can enable the systematic molecular and cellular phenotypic characterization of mouse mutants with unprecedented breadth and resolution.

STIGMA: Single-cell tissue-specific gene prioritization using machine learning.

Clinical exome and genome sequencing have revolutionized the understanding of human disease genetics. Yet many genes remain functionally uncharacterized, complicating the establishment of causal disease links for genetic variants. While several scoring methods have been devised to prioritize these candidate genes, these methods fall short of capturing the expression heterogeneity across cell subpopulations within tissues. Here, we introduce single-cell tissue-specific gene prioritization using machine learning (STIGMA), an approach that leverages single-cell RNA-seq (scRNA-seq) data to prioritize candidate genes associated with rare congenital diseases. STIGMA prioritizes genes by learning the temporal dynamics of gene expression across cell types during healthy organogenesis. To assess the efficacy of our framework, we applied STIGMA to mouse limb and human fetal heart scRNA-seq datasets. In a cohort of individuals with congenital limb malformation, STIGMA prioritized 469 variants in 345 genes, with UBA2 as a notable example. For congenital heart defects, we detected 34 genes harboring nonsynonymous de novo variants (nsDNVs) in two or more individuals from a set of 7,958 individuals, including the ortholog of Prdm1, which is associated with hypoplastic left ventricle and hypoplastic aortic arch. Overall, our findings demonstrate that STIGMA effectively prioritizes tissue-specific candidate genes by utilizing single-cell transcriptome data. The ability to capture the heterogeneity of gene expression across cell populations makes STIGMA a powerful tool for the discovery of disease-associated genes and facilitates the identification of causal variants underlying human genetic disorders.

Single-cell RNA-based phenotyping reveals a pivotal role of thyroid hormone receptor alpha for hypothalamic development.

Thyroid hormone and its receptor TRα1 play an important role in brain development. Several animal models have been used to investigate this function, including mice heterozygous for the TRα1R384C mutation, which confers receptor-mediated hypothyroidism. These mice display abnormalities in several autonomic functions, which was partially attributed to a developmental defect in hypothalamic parvalbumin neurons. However, whether other cell types in the hypothalamus are similarly affected remains unknown. Here, we used single-nucleus RNA sequencing to obtain an unbiased view on the importance of TRα1 for hypothalamic development and cellular diversity. Our data show that defective TRα1 signaling has surprisingly little effect on the development of hypothalamic neuronal populations, but it heavily affects hypothalamic oligodendrocytes. Using selective reactivation of the mutant TRα1 during specific developmental periods, we find that early postnatal thyroid hormone action seems to be crucial for proper hypothalamic oligodendrocyte maturation. Taken together, our findings underline the well-known importance of postnatal thyroid health for brain development and provide an unbiased roadmap for the identification of cellular targets of TRα1 action in mouse hypothalamic development.

Input-specific control of interneuron numbers in nascent striatal networks.

The assembly of functional neuronal circuits requires appropriate numbers of distinct classes of neurons, but the mechanisms through which their relative proportions are established remain poorly defined. Investigating the mouse striatum, we found that the two most prominent subtypes of striatal interneurons, parvalbumin-expressing (PV+) GABAergic and cholinergic (ChAT+) interneurons, undergo extensive programmed cell death between the first and second postnatal weeks. Remarkably, the survival of PV+ and ChAT+ interneurons is regulated by distinct mechanisms mediated by their specific afferent connectivity. While long-range cortical inputs control PV+ interneuron survival, ChAT+ interneuron survival is regulated by local input from the medium spiny neurons. Our results identify input-specific circuit mechanisms that operate during the period of programmed cell death to establish the final number of interneurons in nascent striatal networks.

Spindle reorientation in response to mechanical stress is an emergent property of the spindle positioning mechanisms.

Proper orientation of the mitotic spindle plays a crucial role in embryos, during tissue development, and in adults, where it functions to dissipate mechanical stress to maintain tissue integrity and homeostasis. While mitotic spindles have been shown to reorient in response to external mechanical stresses, the subcellular cues that mediate spindle reorientation remain unclear. Here, we used a combination of optogenetics and computational modeling to investigate how mitotic spindles respond to inhomogeneous tension within the actomyosin cortex. Strikingly, we found that the optogenetic activation of RhoA only influences spindle orientation when it is induced at both poles of the cell. Under these conditions, the sudden local increase in cortical tension induced by RhoA activation reduces pulling forces exerted by cortical regulators on astral microtubules. This leads to a perturbation of the balance of torques exerted on the spindle, which causes it to rotate. Thus, spindle rotation in response to mechanical stress is an emergent phenomenon arising from the interaction between the spindle positioning machinery and the cell cortex.

Position effects at the FGF8 locus are associated with femoral hypoplasia.

Copy-number variations (CNVs) are a common cause of congenital limb malformations and are interpreted primarily on the basis of their effect on gene dosage. However, recent studies show that CNVs also influence the 3D genome chromatin organization. The functional interpretation of whether a phenotype is the result of gene dosage or a regulatory position effect remains challenging. Here, we report on two unrelated families with individuals affected by bilateral hypoplasia of the femoral bones, both harboring de novo duplications on chromosome 10q24.32. The ∼0.5 Mb duplications include FGF8, a key regulator of limb development and several limb enhancer elements. To functionally characterize these variants, we analyzed the local chromatin architecture in the affected individuals' cells and re-engineered the duplications in mice by using CRISPR-Cas9 genome editing. We found that the duplications were associated with ectopic chromatin contacts and increased FGF8 expression. Transgenic mice carrying the heterozygous tandem duplication including Fgf8 exhibited proximal shortening of the limbs, resembling the human phenotype. To evaluate whether the phenotype was a result of gene dosage, we generated another transgenic mice line, carrying the duplication on one allele and a concurrent Fgf8 deletion on the other allele, as a control. Surprisingly, the same malformations were observed. Capture Hi-C experiments revealed ectopic interaction with the duplicated region and Fgf8, indicating a position effect. In summary, we show that duplications at the FGF8 locus are associated with femoral hypoplasia and that the phenotype is most likely the result of position effects altering FGF8 expression rather than gene dosage effects.

Pyramidal cell regulation of interneuron survival sculpts cortical networks.

Complex neuronal circuitries such as those found in the mammalian cerebral cortex have evolved as balanced networks of excitatory and inhibitory neurons. Although the establishment of appropriate numbers of these cells is essential for brain function and behaviour, our understanding of this fundamental process is limited. Here we show that the survival of interneurons in mice depends on the activity of pyramidal cells in a critical window of postnatal development, during which excitatory synaptic input to individual interneurons predicts their survival or death. Pyramidal cells regulate interneuron survival through the negative modulation of PTEN signalling, which effectively drives interneuron cell death during this period. Our findings indicate that activity-dependent mechanisms dynamically adjust the number of inhibitory cells in nascent local cortical circuits, ultimately establishing the appropriate proportions of excitatory and inhibitory neurons in the cerebral cortex.

Comparative single-cell analysis of the adult heart and coronary vasculature.

The structure and function of the circulatory system, including the heart, have undergone substantial changes with the vertebrate evolution. Although the basic function of the heart is to pump blood through the body, its size, shape, speed, regeneration capacity, etc. vary considerably across species. Here, we address the differences among vertebrate hearts using a single-cell transcriptomics approach. Published datasets of macaque (Macaca fascicularis), mouse, and zebrafish hearts were integrated and compared to the human heart as a reference. While the three mammalian hearts integrated well, the zebrafish heart showed very little overlap with the other species. Our analysis revealed a mouse-specific cell subpopulation of ventricular cardiomyocytes (CM), represented by strikingly different expression patterns of specific genes related to high-energy metabolism. Interestingly, the observed differences between mouse and human CM coincided with actual biological differences between the two species. Smooth muscle and endothelial cells (EC) exhibited species-specific differences in clustering and gene expression, respectively, which we attribute to the tissues selected for sequencing, given different focuses of the original studies. Finally, we compared human and zebrafish heart-specific fibroblasts (FB) and identified a distinctively high expression of genes associated with heart regeneration following injury in zebrafish. Together, our results show that integration of numerous datasets of different species and different sequencing technologies is feasible and that this approach can identify species-specific differences and similarities in the heart.

The role of single-cell genomics in human genetics.

Single-cell sequencing is a powerful approach that can detect genetic alterations and their phenotypic consequences in the context of human development, with cellular resolution. Humans start out as single-cell zygotes and undergo fission and differentiation to develop into multicellular organisms. Before fertilisation and during development, the cellular genome acquires hundreds of mutations that propagate down the cell lineage. Whether germline or somatic in nature, some of these mutations may have significant genotypic impact and lead to diseased cellular phenotypes, either systemically or confined to a tissue. Single-cell sequencing enables the detection and monitoring of the genotype and the consequent molecular phenotypes at a cellular resolution. It offers powerful tools to compare the cellular lineage between 'normal' and 'diseased' conditions and to establish genotype-phenotype relationships. By preserving cellular heterogeneity, single-cell sequencing, unlike bulk-sequencing, allows the detection of even small, diseased subpopulations of cells within an otherwise normal tissue. Indeed, the characterisation of biopsies with cellular resolution can provide a mechanistic view of the disease. While single-cell approaches are currently used mainly in basic research, it can be expected that applications of these technologies in the clinic may aid the detection, diagnosis and eventually the treatment of rare genetic diseases as well as cancer. This review article provides an overview of the single-cell sequencing technologies in the context of human genetics, with an aim to empower clinicians to understand and interpret the single-cell sequencing data and analyses. We discuss the state-of-the-art experimental and analytical workflows and highlight current challenges/limitations. Notably, we focus on two prospective applications of the technology in human genetics, namely the annotation of the non-coding genome using single-cell functional genomics and the use of single-cell sequencing data for in silico variant prioritisation.

Transcriptional Alterations in X-Linked Dystonia-Parkinsonism Caused by the SVA Retrotransposon.

X-linked dystonia-parkinsonism (XDP) is a severe neurodegenerative disorder that manifests as adult-onset dystonia combined with parkinsonism. A SINE-VNTR-Alu (SVA) retrotransposon inserted in an intron of the TAF1 gene reduces its expression and alters splicing in XDP patient-derived cells. As a consequence, increased levels of the TAF1 intron retention transcript TAF1-32i can be found in XDP cells as compared to healthy controls. Here, we investigate the sequence of the deep intronic region included in this transcript and show that it is also present in cells from healthy individuals, albeit in lower amounts than in XDP cells, and that it undergoes degradation by nonsense-mediated mRNA decay. Furthermore, we investigate epigenetic marks (e.g., DNA methylation and histone modifications) present in this intronic region and the spanning sequence. Finally, we show that the SVA evinces regulatory potential, as demonstrated by its ability to repress the TAF1 promoter in vitro. Our results enable a better understanding of the disease mechanisms underlying XDP and transcriptional alterations caused by SVA retrotransposons.

Intrinsic Efficacy of Opioid Ligands and Its Importance for Apparent Bias, Operational Analysis, and Therapeutic Window.

Evidence from several novel opioid agonists and knockout animals suggests that improved opioid therapeutic window, notably for analgesia versus respiratory depression, is a result of ligand bias downstream of activation of the µ-opioid receptor (MOR) toward G protein signaling and away from other pathways, such as arrestin recruitment. Here, we argue that published claims of opioid bias based on application of the operational model of agonism are frequently confounded by failure to consider the assumptions of the model. These include failure to account for intrinsic efficacy and ceiling effects in different pathways, distortions introduced by analysis of amplified (G protein) versus linear (arrestin) signaling mechanisms, and nonequilibrium effects in a dynamic signaling cascade. We show on both theoretical and experimental grounds that reduced intrinsic efficacy that is unbiased across different downstream pathways, when analyzed without due considerations, does produce apparent but erroneous MOR ligand bias toward G protein signaling, and the weaker the G protein partial agonism is the greater the apparent bias. Experimentally, such apparently G protein-biased opioids have been shown to exhibit low intrinsic efficacy for G protein signaling when ceiling effects are properly accounted for. Nevertheless, such agonists do display an improved therapeutic window for analgesia versus respiratory depression. Reduced intrinsic efficacy for G proteins rather than any supposed G protein bias provides a more plausible, sufficient explanation for the improved safety. Moreover, genetic models of G protein-biased opioid receptors and replication of previous knockout experiments suggest that reduced or abolished arrestin recruitment does not improve therapeutic window for MOR-induced analgesia versus respiratory depression. SIGNIFICANCE STATEMENT: Efforts to improve safety of µ-opioid analgesics have focused on agonists that show signaling bias for the G protein pathway versus other signaling pathways. This review provides theoretical and experimental evidence showing that failure to consider the assumptions of the operational model can lead to large distortions and overestimation of actual bias. We show that low intrinsic efficacy is a major determinant of these distortions, and pursuit of appropriately reduced intrinsic efficacy should guide development of safer opioids.

Electromechanics and Volume Dynamics in Nonexcitable Tissue Cells.

Cell volume regulation is fundamentally important in phenomena such as cell growth, proliferation, tissue homeostasis, and embryogenesis. How the cell size is set, maintained, and changed over a cell's lifetime is not well understood. In this work we focus on how the volume of nonexcitable tissue cells is coupled to the cell membrane electrical potential and the concentrations of membrane-permeable ions in the cell environment. Specifically, we demonstrate that a sudden cell depolarization using the whole-cell patch clamp results in a 50% increase in cell volume, whereas hyperpolarization results in a slight volume decrease. We find that cell volume can be partially controlled by changing the chloride or the sodium/potassium concentrations in the extracellular environment while maintaining a constant external osmotic pressure. Depletion of external chloride leads to a volume decrease in suspended HN31 cells. Introducing cells to a high-potassium solution causes volume increase up to 50%. Cell volume is also influenced by cortical tension: actin depolymerization leads to cell volume increase. We present an electrophysiology model of water dynamics driven by changes in membrane potential and the concentrations of permeable ions in the cells surrounding. The model quantitatively predicts that the cell volume is directly proportional to the intracellular protein content.

Human Prestin: A Candidate PE1 Protein Lacking Stringent Mass Spectrometric Evidence?

The evidence that any protein exists in the Human Proteome Project (HPP; protein evidence 1 or PE1) has revolved primarily (although not exclusively) around mass spectrometry (MS) (93% of PE1 proteins have MS evidence in the latest neXtProt release), with robust and stringent, well-curated metrics that have served the community well. This has led to a significant number of proteins still considered "missing" (i.e., PE2-4). Many PE2-4 proteins have MS evidence of unacceptable quality (small or not enough unitypic peptides and unacceptably high protein/peptide FDRs), transcriptomic, or antibody evidence. Here we use a Chromosome 7 PE2 example called Prestin to demonstrate that clear and robust criteria/metrics need to be developed for proteins that may not or cannot produce clear-cut MS evidence while possessing significant non-MS evidence, including disease-association data. Many of the PE2-4 proteins are inaccessible, spatiotemporally expressed in a limited way, or expressed at such a very low copy number as to be unable to be detected by current MS methodologies. We propose that the HPP community consider and lead a communal initiative to accelerate the discovery and characterization of these types of "missing" proteins.

Parallel pathways from motor and somatosensory cortex for controlling whisker movements in mice.

Mice can gather tactile sensory information by actively moving their whiskers to palpate objects in their immediate surroundings. Whisker sensory perception therefore requires integration of sensory and motor information, which occurs prominently in the neocortex. The signalling pathways from the neocortex for controlling whisker movements are currently poorly understood in mice. Here, we delineate two pathways, one originating from primary whisker somatosensory cortex (wS1) and the other from whisker motor cortex (wM1), that control qualitatively distinct movements of contralateral whiskers. Optogenetic stimulation of wS1 drove retraction of contralateral whiskers while stimulation of wM1 drove rhythmic whisker protraction. To map brainstem pathways connecting these cortical areas to whisker motor neurons, we used a combination of anterograde tracing using adenoassociated virus injected into neocortex and retrograde tracing using monosynaptic rabies virus injected into whisker muscles. Our data are consistent with wS1 driving whisker retraction by exciting glutamatergic premotor neurons in the rostral spinal trigeminal interpolaris nucleus, which in turn activate the motor neurons innervating the extrinsic retractor muscle nasolabialis. The rhythmic whisker protraction evoked by wM1 stimulation might be driven by excitation of excitatory and inhibitory premotor neurons in the brainstem reticular formation innervating both intrinsic and extrinsic muscles. Our data therefore begin to unravel the neuronal circuits linking the neocortex to whisker motor neurons.

A modular design of low-background bioassays based on a high-affinity molecular pair barstar:barnase.

High-affinity molecular pairs provide a convenient and flexible modular base for the design of molecular probes and protein/antigen assays. Specificity and sensitivity performance indicators of a bioassay critically depend on the dissociation constant (K(D)) of the molecular pair, with avidin:biotin being the state-of-the-art molecular pair (K(D) ∼ 1 fM) used almost universally for applications in the fields of nanotechnology and proteomics. In this paper, we present an alternative high-affinity protein pair, barstar:barnase (K(D) ∼ 10 fM), which addresses several shortfalls of the avidin:biotin system, including non-negligible background due to the non-specific binding. A quantitative assessment of the non-specific binding carried out using a model assay revealed inherent irreproducibility of the [strept]avidin:biotin-based assays, attributed to the avidin binding to solid phases, endogenous biotin molecules and serum proteins. On the other hand, the model assays assembled via a barstar:barnase protein linker proved to be immune to such non-specific binding, showing good prospects for high-sensitivity rare biomolecular event nanoproteomic assays.

Prediction by promoter logic in bacterial quorum sensing.

Quorum-sensing systems mediate chemical communication between bacterial cells, coordinating cell-density-dependent processes like biofilm formation and virulence-factor expression. In the proteobacterial LuxI/LuxR quorum sensing paradigm, a signaling molecule generated by an enzyme (LuxI) diffuses between cells and allosterically stimulates a transcriptional regulator (LuxR) to activate its cognate promoter (pR). By expressing either LuxI or LuxR in positive feedback from pR, these versatile systems can generate smooth (monostable) or abrupt (bistable) density-dependent responses to suit the ecological context. Here we combine theory and experiment to demonstrate that the promoter logic of pR - its measured activity as a function of LuxI and LuxR levels - contains all the biochemical information required to quantitatively predict the responses of such feedback loops. The interplay of promoter logic with feedback topology underlies the versatility of the LuxI/LuxR paradigm: LuxR and LuxI positive-feedback systems show dramatically different responses, while a dual positive/negative-feedback system displays synchronized oscillations. These results highlight the dual utility of promoter logic: to probe microscopic parameters and predict macroscopic phenotype.

Acute hypertensive response management in patients with acute stroke.

Acute elevation in blood pressure (acute hypertensive response) is commonly observed in the early period of both ischemic and hemorrhagic stroke. The management of acute hypertensive response depends upon the underlying intracranial pathology. Management of acute hypertensive response has been the focus of many trials and studies such as the SCAST trial, CHHIPS trial, COSSACS trial, INTERACT, and ATACH trial, which are discussed here. However, there were many limitations to these trials including time of presentation, different pathophysiology of ischemic strokes versus hemorrhagic strokes, and patient related factors. Ongoing clinical trials which take these limitations into account include ENCHANTED trial, ATACH II trial, INTERACT 2 trial, and ENOS trial. The results of these trials are expected to modify current guidelines of acute stroke, both ischemic and hemorrhagic, and potentially improve clinical outcomes and quality of life.

Nature of the amyloid-beta monomer and the monomer-oligomer equilibrium.

The monomer to oligomer transition initiates the aggregation and pathogenic transformation of Alzheimer amyloid-ß (Aß) peptide. However, the monomeric state of this aggregation-prone peptide has remained beyond the reach of most experimental techniques, and a quantitative understanding of this transition is yet to emerge. Here, we employ single-molecule level fluorescence tools to characterize the monomeric state and the monomer-oligomer transition at physiological concentrations in buffers mimicking the cerebrospinal fluid (CSF). Our measurements show that the monomer has a hydrodynamic radius of 0.9 ± 0.1 nm, which confirms the prediction made by some of the in silico studies. Surprisingly, at equilibrium, both Aß(40) and Aß(42) remain predominantly monomeric up to 3 µm, above which it forms large aggregates. This concentration is much higher than the estimated concentrations in the CSF of either normal or diseased brains. If Aß oligomers are present in the CSF and are the key agents in Alzheimer pathology, as is generally believed, then these must be released in the CSF as preformed entities. Although the oligomers are thermodynamically unstable, we find that a large kinetic barrier, which is mostly entropic in origin, strongly impedes their dissociation. Thermodynamic principles therefore allow the development of a pharmacological agent that can catalytically convert metastable oligomers into nontoxic monomers.

Single-cell sequencing: promises and challenges for human genetics.

Over the last decade, single-cell sequencing has transformed many fields. It has enabled the unbiased molecular phenotyping of even whole organisms with unprecedented cellular resolution. In the field of human genetics, where the phenotypic consequences of genetic and epigenetic alterations are of central concern, this transformative technology promises to functionally annotate every region in the human genome and all possible variants within them at a massive scale. In this review aimed at the clinicians in human genetics, we describe the current status of the field of single-cell sequencing and its role for human genetics, including how the technology works as well as how it is being applied to characterize and monitor diseases, to develop human cell atlases, and to annotate the genome.