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.

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.

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.

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.

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.

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.

Pharmacological characterization of a recombinant, fluorescent somatostatin receptor agonist.

Somatostatin (SST) is a peptide neurotransmitter/hormone found in several mammalian tissue types. Apart from its natural importance, labeled SST/analogues are utilized in clinical applications such as targeting/diagnosis of neuroendocrine tumors. We report on the development and characterization of a novel, recombinant, fluorescent somatostatin analogue that has potential to elucidate somatostatin-activated cell signaling. SST was genetically fused with a monomeric-red fluorescent protein (mRFP) as the fluorescent label. The attachment of SST to mRFP had no detectable effect on its fluorescent properties. This analogue's potency to activate the endogenous and transfected somatostatin receptors was characterized using assays of membrane potential and Ca(2+) mobilization and immunocytochemistry. SST-mRFP was found to be an effective somatostatin receptor agonist, able to trigger the membrane hyperpolarization, mobilization of the intracellular Ca(2+) and receptor-ligand internalization in cells expressing somatostatin receptors. This complex represents a novel optical reporter due to its red emission spectral band suitable for in vivo imaging and tracking of the somatostatin receptor signaling pathways, affording higher resolution and sensitivity than those of the state-of-the-art radiolabeling bioassays.

Functionalization and Bioconjugation of Nanoruby for Long-Term, Ultrasensitive Imaging of Μu-Opioid Receptors.

Sensitive and long-term fluorescence imaging of G-protein-coupled receptors enables exploration of molecular level details of these therapeutically relevant proteins, including their expression, localization, signaling, and intracellular trafficking. In this context, labeling these receptors with bright and photostable fluorescent probes is necessary to overcome current imaging problems such as optical background and photobleaching. Here, we describe the procedures to functionalize nanoruby (and other similar nanoparticles) with NeutrAvidin (a streptavidin analog) and to apply this bioconjugate for ultrasensitive, long-term imaging of µ-opioid receptors heterologously expressed in AtT-20 cells. The receptor targeting is mediated via a biotinylated primary antibody, rendering this methodology extendable to other G-protein-coupled or, more generally, cell-surface receptors. Nanoruby-based time-gated imaging enables indefinitely long visualization of single particles even in high-autofluorescence media, such as serum, by completely suppressing autofluorescence and any laser backscatter.

Lifetime-Engineered Ruby Nanoparticles (Tau-Rubies) for Multiplexed Imaging of µ-Opioid Receptors.

To address the growing demand for simultaneous imaging of multiple biomarkers in highly scattering media such as organotypic cell cultures, we introduce a new type of photoluminescent nanomaterial termed "tau-ruby" composed of ruby nanocrystals (Al2O3:Cr3+) with tunable emission lifetime. The lifetime tuning range from 2.4 to 3.2 ms was achieved by varying the Cr3+ dopant concentration from 0.8% to 0.2%, affording facile implementation of background-free detection. We developed inexpensive scalable production of tau-ruby characterized by bright emission, narrow spectrum (693 ± 2 nm), and virtually unlimited photostability upon excitation with affordable excitation/detection sources, noncytotoxic and insensitive to microenvironmental fluctuations. By functionalizing the surface of tau-rubies with targeting antibodies, we obtained different biomarkers suitable for multiplexed lifetime imaging. As a proof of principle, three tau-ruby bioprobes, characterized by three mean lifetimes, were deployed to label three µ-opioid receptor species expressed on transfected cancer cells, each fused to a unique epitope, so that three types of cells were lifetime-encoded. Robust decoding of photoluminescent signals that report on each cell type was achieved by using a home-built lifetime imaging system and resulted in high-contrast multiplexed lifetime imaging of the cells.

Influence of FRET and fluorescent protein maturation on the quantification of binding affinity with dual-channel fluorescence cross-correlation spectroscopy.

Protein-protein interactions at the plasma membrane mediate transmembrane signaling. Dual-channel fluorescence cross-correlation spectroscopy (dc-FCCS) is a method with which these interactions can be quantified in a cellular context. However, factors such as incomplete maturation of fluorescent proteins, spectral crosstalk, and fluorescence resonance energy transfer (FRET) affect quantification. Some of these can be corrected or accounted for during data analysis and/or interpretation. Here, we experimentally and analytically demonstrate that it is difficult to correct the error caused due to FRET when applying dc-FCCS to measure binding affinity or bound molecular concentrations. Additionally, the presence of dark fluorescent proteins due to incomplete maturation introduces further errors, which too cannot be corrected in the presence of FRET. Based on simulations, we find that modalities such as pulse-interleaved excitation FCCS do not eliminate FRET-induced errors. Finally, we demonstrate that the detrimental effect of FRET can be eliminated with careful experimental design when applying dc-FCCS to quantify protein-protein interactions at the plasma membrane of living cells.

Force Measurements for Cancer Cells.

During cytoskeleton remodeling, cancer cells generate force at the plasma membrane that originates from chemical motors (e.g., actin). This force (pN) and its time course reflect the on and off-rates of the motors. We describe the design and calibration of a force-measuring device (i.e., optical tweezers) that is used to monitor this force and its time course at the edge of a cell, with particular emphasis on the temporal resolution of the instrument.

Development of Bright and Biocompatible Nanoruby and Its Application to Background-Free Time-Gated Imaging of G-Protein-Coupled Receptors.

At the forefront of developing fluorescent probes for biological imaging applications are enhancements aimed at increasing their brightness, contrast, and photostability, especially toward demanding applications of single-molecule detection. In comparison with existing probes, nanorubies exhibit unlimited photostability and a long emission lifetime (∼4 ms), which enable continuous imaging at single-particle sensitivity in highly scattering and fluorescent biological specimens. However, their wide application as fluorescence probes has so far been hindered by the absence of facile methods for scaled-up high-volume production and molecularly specific targeting. The present work encompasses the large-scale production of colloidally stable nanoruby particles, the demonstration of their biofunctionality and negligible cytotoxicity, as well as the validation of its use for targeted biomolecular imaging. In addition, optical characteristics of nanorubies are found to be comparable or superior to those of state-of-the-art quantum dots. Protocols of reproducible and robust coupling of functional proteins to the nanoruby surface are also presented. As an example, NeutrAvidin-coupled nanoruby show excellent affinity and specificity to µ-opioid receptors in fixed and live cells, allowing wide-field imaging of G-protein coupled receptors with single-particle sensitivity.

In vivo imaging of endothelial cell adhesion molecule expression after radiosurgery in an animal model of arteriovenous malformation.

Focussed radiosurgery may provide a means of inducing molecular changes on the luminal surface of diseased endothelium to allow targeted delivery of novel therapeutic compounds. We investigated the potential of ionizing radiation to induce surface expression of intercellular adhesion molecule 1 (ICAM-1) and vascular cell adhesion molecule 1 (VCAM-1) on endothelial cells (EC) in vitro and in vivo, to assess their suitability as vascular targets in irradiated arteriovenous malformations (AVMs). Cultured brain microvascular EC were irradiated by linear accelerator at single doses of 0, 5, 15 or 25 Gy and expression of ICAM-1 and VCAM-1 measured by qRT-PCR, Western, ELISA and immunocytochemistry. In vivo, near-infrared (NIR) fluorescence optical imaging using Xenolight 750-conjugated ICAM-1 or VCAM-1 antibodies examined luminal biodistribution over 84 days in a rat AVM model after Gamma Knife surgery at a single 15 Gy dose. ICAM-1 and VCAM-1 were minimally expressed on untreated EC in vitro. Doses of 15 and 25 Gy stimulated expression equally; 5 Gy was not different from the unirradiated. In vivo, normal vessels did not bind or retain the fluorescent probes, however binding was significant in AVM vessels. No additive increases in probe binding were found in response to radiosurgery at a dose of 15 Gy. In summary, radiation induces adhesion molecule expression in vitro but elevated baseline levels in AVM vessels precludes further induction in vivo. These molecules may be suitable targets in irradiated vessels without hemodynamic derangement, but not AVMs. These findings demonstrate the importance of using flow-modulated, pre-clinical animal models for validating candidate proteins for vascular targeting in irradiated AVMs.