Defining the Design Principles of Skin Epidermis Postnatal Growth.

During embryonic and postnatal development, organs and tissues grow steadily to achieve their final size at the end of puberty. However, little is known about the cellular dynamics that mediate postnatal growth. By combining in vivo clonal lineage tracing, proliferation kinetics, single-cell transcriptomics, and in vitro micro-pattern experiments, we resolved the cellular dynamics taking place during postnatal skin epidermis expansion. Our data revealed that harmonious growth is engineered by a single population of developmental progenitors presenting a fixed fate imbalance of self-renewing divisions with an ever-decreasing proliferation rate. Single-cell RNA sequencing revealed that epidermal developmental progenitors form a more uniform population compared with adult stem and progenitor cells. Finally, we found that the spatial pattern of cell division orientation is dictated locally by the underlying collagen fiber orientation. Our results uncover a simple design principle of organ growth where progenitors and differentiated cells expand in harmony with their surrounding tissues.

Methods and applications for single-cell and spatial multi-omics.

The joint analysis of the genome, epigenome, transcriptome, proteome and/or metabolome from single cells is transforming our understanding of cell biology in health and disease. In less than a decade, the field has seen tremendous technological revolutions that enable crucial new insights into the interplay between intracellular and intercellular molecular mechanisms that govern development, physiology and pathogenesis. In this Review, we highlight advances in the fast-developing field of single-cell and spatial multi-omics technologies (also known as multimodal omics approaches), and the computational strategies needed to integrate information across these molecular layers. We demonstrate their impact on fundamental cell biology and translational research, discuss current challenges and provide an outlook to the future.

Dedicated macrophages organize and maintain the enteric nervous system.

Correct development and maturation of the enteric nervous system (ENS) is critical for survival1. At birth, the ENS is immature and requires considerable refinement to exert its functions in adulthood2. Here we demonstrate that resident macrophages of the muscularis externa (MMϕ) refine the ENS early in life by pruning synapses and phagocytosing enteric neurons. Depletion of MMϕ before weaning disrupts this process and results in abnormal intestinal transit. After weaning, MMϕ continue to interact closely with the ENS and acquire a neurosupportive phenotype. The latter is instructed by transforming growth factor-ß produced by the ENS; depletion of the ENS and disruption of transforming growth factor-ß signalling result in a decrease in neuron-associated MMϕ associated with loss of enteric neurons and altered intestinal transit. These findings introduce a new reciprocal cell-cell communication responsible for maintenance of the ENS and indicate that the ENS, similarly to the brain, is shaped and maintained by a dedicated population of resident macrophages that adapts its phenotype and transcriptome to the timely needs of the ENS niche.

The extracellular matrix dictates regional competence for tumour initiation.

The skin epidermis is constantly renewed throughout life1,2. Disruption of the balance between renewal and differentiation can lead to uncontrolled growth and tumour initiation3. However, the ways in which oncogenic mutations affect the balance between renewal and differentiation and lead to clonal expansion, cell competition, tissue colonization and tumour development are unknown. Here, through multidisciplinary approaches that combine in vivo clonal analysis using intravital microscopy, single-cell analysis and functional analysis, we show how SmoM2-a constitutively active oncogenic mutant version of Smoothened (SMO) that induces the development of basal cell carcinoma-affects clonal competition and tumour initiation in real time. We found that expressing SmoM2 in the ear epidermis of mice induced clonal expansion together with tumour initiation and invasion. By contrast, expressing SmoM2 in the back-skin epidermis led to a clonal expansion that induced lateral cell competition without dermal invasion and tumour formation. Single-cell analysis showed that oncogene expression was associated with a cellular reprogramming of adult interfollicular cells into an embryonic hair follicle progenitor (EHFP) state in the ear but not in the back skin. Comparisons between the ear and the back skin revealed that the dermis has a very different composition in these two skin types, with increased stiffness and a denser collagen I network in the back skin. Decreasing the expression of collagen I in the back skin through treatment with collagenase, chronic UV exposure or natural ageing overcame the natural resistance of back-skin basal cells to undergoing EHFP reprogramming and tumour initiation after SmoM2 expression. Altogether, our study shows that the composition of the extracellular matrix regulates how susceptible different regions of the body are to tumour initiation and invasion.

Pharmacological targeting of netrin-1 inhibits EMT in cancer.

Epithelial-to-mesenchymal transition (EMT) regulates tumour initiation, progression, metastasis and resistance to anti-cancer therapy1-7. Although great progress has been made in understanding the role of EMT and its regulatory mechanisms in cancer, no therapeutic strategy to pharmacologically target EMT has been identified. Here we found that netrin-1 is upregulated in a primary mouse model of skin squamous cell carcinoma (SCC) exhibiting spontaneous EMT. Pharmacological inhibition of netrin-1 by administration of NP137, a netrin-1-blocking monoclonal antibody currently used in clinical trials in human cancer (ClinicalTrials.gov identifier NCT02977195 ), decreased the proportion of EMT tumour cells in skin SCC, decreased the number of metastases and increased the sensitivity of tumour cells to chemotherapy. Single-cell RNA sequencing revealed the presence of different EMT states, including epithelial, early and late hybrid EMT, and full EMT states, in control SCC. By contrast, administration of NP137 prevented the progression of cancer cells towards a late EMT state and sustained tumour epithelial states. Short hairpin RNA knockdown of netrin-1 and its receptor UNC5B in EPCAM+ tumour cells inhibited EMT in vitro in the absence of stromal cells and regulated a common gene signature that promotes tumour epithelial state and restricts EMT. To assess the relevance of these findings to human cancers, we treated mice transplanted with the A549 human cancer cell line-which undergoes EMT following TGFß1 administration8,9-with NP137. Netrin-1 inhibition decreased EMT in these transplanted A549 cells. Together, our results identify a pharmacological strategy for targeting EMT in cancer, opening up novel therapeutic interventions for anti-cancer therapy.

Mechanisms of stretch-mediated skin expansion at single-cell resolution.

The ability of the skin to grow in response to stretching has been exploited in reconstructive surgery1. Although the response of epidermal cells to stretching has been studied in vitro2,3, it remains unclear how mechanical forces affect their behaviour in vivo. Here we develop a mouse model in which the consequences of stretching on skin epidermis can be studied at single-cell resolution. Using a multidisciplinary approach that combines clonal analysis with quantitative modelling and single-cell RNA sequencing, we show that stretching induces skin expansion by creating a transient bias in the renewal activity of epidermal stem cells, while a second subpopulation of basal progenitors remains committed to differentiation. Transcriptional and chromatin profiling identifies how cell states and gene-regulatory networks are modulated by stretching. Using pharmacological inhibitors and mouse mutants, we define the step-by-step mechanisms that control stretch-mediated tissue expansion at single-cell resolution in vivo.

Heterotypic cell-cell communication regulates glandular stem cell multipotency.

Glandular epithelia, including the mammary and prostate glands, are composed of basal cells (BCs) and luminal cells (LCs)1,2. Many glandular epithelia develop from multipotent basal stem cells (BSCs) that are replaced in adult life by distinct pools of unipotent stem cells1,3-8. However, adult unipotent BSCs can reactivate multipotency under regenerative conditions and upon oncogene expression3,9-13. This suggests that an active mechanism restricts BSC multipotency under normal physiological conditions, although the nature of this mechanism is unknown. Here we show that the ablation of LCs reactivates the multipotency of BSCs from multiple epithelia both in vivo in mice and in vitro in organoids. Bulk and single-cell RNA sequencing revealed that, after LC ablation, BSCs activate a hybrid basal and luminal cell differentiation program before giving rise to LCs-reminiscent of the genetic program that regulates multipotency during embryonic development7. By predicting ligand-receptor pairs from single-cell data14, we find that TNF-which is secreted by LCs-restricts BC multipotency under normal physiological conditions. By contrast, the Notch, Wnt and EGFR pathways were activated in BSCs and their progeny after LC ablation; blocking these pathways, or stimulating the TNF pathway, inhibited regeneration-induced BC multipotency. Our study demonstrates that heterotypic communication between LCs and BCs is essential to maintain lineage fidelity in glandular epithelial stem cells.

Adipose tissue macrophage dysfunction is associated with a breach of vascular integrity in NASH.

BACKGROUND & AIMS: In non-alcoholic fatty liver disease (NAFLD), monocytes infiltrate visceral adipose tissue promoting local and hepatic inflammation. However, it remains unclear what drives inflammation and how the immune landscape in adipose tissue differs across the NAFLD severity spectrum. We aimed to assess adipose tissue macrophage (ATM) heterogeneity in a NAFLD cohort. METHODS: Visceral adipose tissue macrophages from lean and obese patients, stratified by NAFLD phenotypes, underwent single-cell RNA sequencing. Adipose tissue vascular integrity and breaching was assessed on a protein level via immunohistochemistry and immunofluorescence to determine targets of interest. RESULTS: We discovered multiple ATM populations, including resident vasculature-associated macrophages (ResVAMs) and distinct metabolically active macrophages (MMacs). Using trajectory analysis, we show that ResVAMs and MMacs are replenished by a common transitional macrophage (TransMac) subtype and that, during NASH, MMacs are not effectively replenished by TransMac precursors. We postulate an accessory role for MMacs and ResVAMs in protecting the adipose tissue vascular barrier, since they both interact with endothelial cells and localize around the vasculature. However, across the NAFLD severity spectrum, alterations occur in these subsets that parallel an adipose tissue vasculature breach characterized by albumin extravasation into the perivascular tissue. CONCLUSIONS: NAFLD-related macrophage dysfunction coincides with a loss of adipose tissue vascular integrity, providing a plausible mechanism by which tissue inflammation is perpetuated in adipose tissue and downstream in the liver. IMPACT AND IMPLICATIONS: Our study describes for the first time the myeloid cell landscape in human visceral adipose tissue at single-cell level within a cohort of well-characterized patients with non-alcoholic fatty liver disease. We report unique non-alcoholic steatohepatitis-specific transcriptional changes within metabolically active macrophages (MMacs) and resident vasculature-associated macrophages (ResVAMs) and we demonstrate their spatial location surrounding the vasculature. These dysfunctional transcriptional macrophage states coincided with the loss of adipose tissue vascular integrity, providing a plausible mechanism by which tissue inflammation is perpetuated in adipose tissue and downstream in the liver. Our study provides a theoretical basis for new therapeutic strategies to be directed towards reinstating the endogenous metabolic, homeostatic and cytoprotective functions of ResVAMs and MMacs, including their role in protecting vascular integrity.

Detecting cryptic clinically relevant structural variation in exome-sequencing data increases diagnostic yield for developmental disorders.

Structural variation (SV) describes a broad class of genetic variation greater than 50 bp in size. SVs can cause a wide range of genetic diseases and are prevalent in rare developmental disorders (DDs). Individuals presenting with DDs are often referred for diagnostic testing with chromosomal microarrays (CMAs) to identify large copy-number variants (CNVs) and/or with single-gene, gene-panel, or exome sequencing (ES) to identify single-nucleotide variants, small insertions/deletions, and CNVs. However, individuals with pathogenic SVs undetectable by conventional analysis often remain undiagnosed. Consequently, we have developed the tool InDelible, which interrogates short-read sequencing data for split-read clusters characteristic of SV breakpoints. We applied InDelible to 13,438 probands with severe DDs recruited as part of the Deciphering Developmental Disorders (DDD) study and discovered 63 rare, damaging variants in genes previously associated with DDs missed by standard SNV, indel, or CNV discovery approaches. Clinical review of these 63 variants determined that about half (30/63) were plausibly pathogenic. InDelible was particularly effective at ascertaining variants between 21 and 500 bp in size and increased the total number of potentially pathogenic variants identified by DDD in this size range by 42.9%. Of particular interest were seven confirmed de novo variants in MECP2, which represent 35.0% of all de novo protein-truncating variants in MECP2 among DDD study participants. InDelible provides a framework for the discovery of pathogenic SVs that are most likely missed by standard analytical workflows and has the potential to improve the diagnostic yield of ES across a broad range of genetic diseases.

Differential alternative splicing analysis links variation in ZRSR2 to a novel type of oral-facial-digital syndrome.

PURPOSE: Oral-facial-digital (OFD) syndromes are genetically heterogeneous developmental disorders, caused by pathogenic variants in genes involved in primary cilia formation and function. We identified a previously undescribed type of OFD with brain anomalies, ranging from alobar holoprosencephaly to pituitary anomalies, in 6 unrelated families. METHODS: Exome sequencing of affected probands was supplemented with alternative splicing analysis in patient and control lymphoblastoid and fibroblast cell lines, and primary cilia structure analysis in patient fibroblasts. RESULTS: In 1 family with 2 affected males, we identified a germline variant in the last exon of ZRSR2, NM_005089.4:c.1211_1212del NP_005080.1:p.(Gly404GlufsTer23), whereas 7 affected males from 5 unrelated families were hemizygous for the ZRSR2 variant NM_005089.4:c.1207_1208del NP_005080.1:p.(Arg403GlyfsTer24), either occurring de novo or inherited in an X-linked recessive pattern. ZRSR2, located on chromosome Xp22.2, encodes a splicing factor of the minor spliceosome complex, which recognizes minor introns, representing 0.35% of human introns. Patient samples showed significant enrichment of minor intron retention. Among differentially spliced targets are ciliopathy-related genes, such as TMEM107 and CIBAR1. Primary fibroblasts containing the NM_005089.4:c.1207_1208del ZRSR2 variant had abnormally elongated cilia, confirming an association between defective U12-type intron splicing, OFD and abnormal primary cilia formation. CONCLUSION: We introduce a novel type of OFD associated with elongated cilia and differential splicing of minor intron-containing genes due to germline variation in ZRSR2.

De novo mutations in regulatory elements in neurodevelopmental disorders.

We previously estimated that 42% of patients with severe developmental disorders carry pathogenic de novo mutations in coding sequences. The role of de novo mutations in regulatory elements affecting genes associated with developmental disorders, or other genes, has been essentially unexplored. We identified de novo mutations in three classes of putative regulatory elements in almost 8,000 patients with developmental disorders. Here we show that de novo mutations in highly evolutionarily conserved fetal brain-active elements are significantly and specifically enriched in neurodevelopmental disorders. We identified a significant twofold enrichment of recurrently mutated elements. We estimate that, genome-wide, 1-3% of patients without a diagnostic coding variant carry pathogenic de novo mutations in fetal brain-active regulatory elements and that only 0.15% of all possible mutations within highly conserved fetal brain-active elements cause neurodevelopmental disorders with a dominant mechanism. Our findings represent a robust estimate of the contribution of de novo mutations in regulatory elements to this genetically heterogeneous set of disorders, and emphasize the importance of combining functional and evolutionary evidence to identify regulatory causes of genetic disorders.

Identification of the tumour transition states occurring during EMT.

In cancer, the epithelial-to-mesenchymal transition (EMT) is associated with tumour stemness, metastasis and resistance to therapy. It has recently been proposed that, rather than being a binary process, EMT occurs through distinct intermediate states. However, there is no direct in vivo evidence for this idea. Here we screen a large panel of cell surface markers in skin and mammary primary tumours, and identify the existence of multiple tumour subpopulations associated with different EMT stages: from epithelial to completely mesenchymal states, passing through intermediate hybrid states. Although all EMT subpopulations presented similar tumour-propagating cell capacity, they displayed differences in cellular plasticity, invasiveness and metastatic potential. Their transcriptional and epigenetic landscapes identify the underlying gene regulatory networks, transcription factors and signalling pathways that control these different EMT transition states. Finally, these tumour subpopulations are localized in different niches that differentially regulate EMT transition states.

Integrative analysis of genomic variants reveals new associations of candidate haploinsufficient genes with congenital heart disease.

Numerous genetic studies have established a role for rare genomic variants in Congenital Heart Disease (CHD) at the copy number variation (CNV) and de novo variant (DNV) level. To identify novel haploinsufficient CHD disease genes, we performed an integrative analysis of CNVs and DNVs identified in probands with CHD including cases with sporadic thoracic aortic aneurysm. We assembled CNV data from 7,958 cases and 14,082 controls and performed a gene-wise analysis of the burden of rare genomic deletions in cases versus controls. In addition, we performed variation rate testing for DNVs identified in 2,489 parent-offspring trios. Our analysis revealed 21 genes which were significantly affected by rare CNVs and/or DNVs in probands. Fourteen of these genes have previously been associated with CHD while the remaining genes (FEZ1, MYO16, ARID1B, NALCN, WAC, KDM5B and WHSC1) have only been associated in small cases series or show new associations with CHD. In addition, a systems level analysis revealed affected protein-protein interaction networks involved in Notch signaling pathway, heart morphogenesis, DNA repair and cilia/centrosome function. Taken together, this approach highlights the importance of re-analyzing existing datasets to strengthen disease association and identify novel disease genes and pathways.

Correction: Integrative analysis of genomic variants reveals new associations of candidate haploinsufficient genes with congenital heart disease.

[This corrects the article DOI: 10.1371/journal.pgen.1009679.].

Detection of structural mosaicism from targeted and whole-genome sequencing data.

Structural mosaic abnormalities are large post-zygotic mutations present in a subset of cells and have been implicated in developmental disorders and cancer. Such mutations have been conventionally assessed in clinical diagnostics using cytogenetic or microarray testing. Modern disease studies rely heavily on exome sequencing, yet an adequate method for the detection of structural mosaicism using targeted sequencing data is lacking. Here, we present a method, called MrMosaic, to detect structural mosaic abnormalities using deviations in allele fraction and read coverage from next-generation sequencing data. Whole-exome sequencing (WES) and whole-genome sequencing (WGS) simulations were used to calculate detection performance across a range of mosaic event sizes, types, clonalities, and sequencing depths. The tool was applied to 4911 patients with undiagnosed developmental disorders, and 11 events among nine patients were detected. For eight of these 11 events, mosaicism was observed in saliva but not blood, suggesting that assaying blood alone would miss a large fraction, possibly >50%, of mosaic diagnostic chromosomal rearrangements.

Mutations in Either TUBB or MAPRE2 Cause Circumferential Skin Creases Kunze Type.

Circumferential skin creases Kunze type (CSC-KT) is a specific congenital entity with an unknown genetic cause. The disease phenotype comprises characteristic circumferential skin creases accompanied by intellectual disability, a cleft palate, short stature, and dysmorphic features. Here, we report that mutations in either MAPRE2 or TUBB underlie the genetic origin of this syndrome. MAPRE2 encodes a member of the microtubule end-binding family of proteins that bind to the guanosine triphosphate cap at growing microtubule plus ends, and TUBB encodes a ß-tubulin isotype that is expressed abundantly in the developing brain. Functional analyses of the TUBB mutants show multiple defects in the chaperone-dependent tubulin heterodimer folding and assembly pathway that leads to a compromised yield of native heterodimers. The TUBB mutations also have an impact on microtubule dynamics. For MAPRE2, we show that the mutations result in enhanced MAPRE2 binding to microtubules, implying an increased dwell time at microtubule plus ends. Further, in vivo analysis of MAPRE2 mutations in a zebrafish model of craniofacial development shows that the variants most likely perturb the patterning of branchial arches, either through excessive activity (under a recessive paradigm) or through haploinsufficiency (dominant de novo paradigm). Taken together, our data add CSC-KT to the growing list of tubulinopathies and highlight how multiple inheritance paradigms can affect dosage-sensitive biological systems so as to result in the same clinical defect.

Making new genetic diagnoses with old data: iterative reanalysis and reporting from genome-wide data in 1,133 families with developmental disorders.

PURPOSE: Given the rapid pace of discovery in rare disease genomics, it is likely that improvements in diagnostic yield can be made by systematically reanalyzing previously generated genomic sequence data in light of new knowledge. METHODS: We tested this hypothesis in the United Kingdom-wide Deciphering Developmental Disorders study, where in 2014 we reported a diagnostic yield of 27% through whole-exome sequencing of 1,133 children with severe developmental disorders and their parents. We reanalyzed existing data using improved variant calling methodologies, novel variant detection algorithms, updated variant annotation, evidence-based filtering strategies, and newly discovered disease-associated genes. RESULTS: We are now able to diagnose an additional 182 individuals, taking our overall diagnostic yield to 454/1,133 (40%), and another 43 (4%) have a finding of uncertain clinical significance. The majority of these new diagnoses are due to novel developmental disorder-associated genes discovered since our original publication. CONCLUSION: This study highlights the importance of coupling large-scale research with clinical practice, and of discussing the possibility of iterative reanalysis and recontact with patients and health professionals at an early stage. We estimate that implementing parent-offspring whole-exome sequencing as a first-line diagnostic test for developmental disorders would diagnose >50% of patients.

Beegle: from literature mining to disease-gene discovery.

Disease-gene identification is a challenging process that has multiple applications within functional genomics and personalized medicine. Typically, this process involves both finding genes known to be associated with the disease (through literature search) and carrying out preliminary experiments or screens (e.g. linkage or association studies, copy number analyses, expression profiling) to determine a set of promising candidates for experimental validation. This requires extensive time and monetary resources. We describe Beegle, an online search and discovery engine that attempts to simplify this process by automating the typical approaches. It starts by mining the literature to quickly extract a set of genes known to be linked with a given query, then it integrates the learning methodology of Endeavour (a gene prioritization tool) to train a genomic model and rank a set of candidate genes to generate novel hypotheses. In a realistic evaluation setup, Beegle has an average recall of 84% in the top 100 returned genes as a search engine, which improves the discovery engine by 12.6% in the top 5% prioritized genes. Beegle is publicly available at http://beegle.esat.kuleuven.be/.

Clinical and molecular consequences of disease-associated de novo mutations in SATB2.

PURPOSE: To characterize features associated with de novo mutations affecting SATB2 function in individuals ascertained on the basis of intellectual disability. METHODS: Twenty previously unreported individuals with 19 different SATB2 mutations (11 loss-of-function and 8 missense variants) were studied. Fibroblasts were used to measure mutant protein production. Subcellular localization and mobility of wild-type and mutant SATB2 were assessed using fluorescently tagged protein. RESULTS: Recurrent clinical features included neurodevelopmental impairment (19/19), absent/near absent speech (16/19), normal somatic growth (17/19), cleft palate (9/19), drooling (12/19), and dental anomalies (8/19). Six of eight missense variants clustered in the first CUT domain. Sibling recurrence due to gonadal mosaicism was seen in one family. A nonsense mutation in the last exon resulted in production of a truncated protein retaining all three DNA-binding domains. SATB2 nuclear mobility was mutation-dependent; p.Arg389Cys in CUT1 increased mobility and both p.Gly515Ser in CUT2 and p.Gln566Lys between CUT2 and HOX reduced mobility. The clinical features in individuals with missense variants were indistinguishable from those with loss of function. CONCLUSION: SATB2 haploinsufficiency is a common cause of syndromic intellectual disability. When mutant SATB2 protein is produced, the protein appears functionally inactive with a disrupted pattern of chromatin or matrix association.Genet Med advance online publication 02 February 2017.

Genetic diagnosis of developmental disorders in the DDD study: a scalable analysis of genome-wide research data.

BACKGROUND: Human genome sequencing has transformed our understanding of genomic variation and its relevance to health and disease, and is now starting to enter clinical practice for the diagnosis of rare diseases. The question of whether and how some categories of genomic findings should be shared with individual research participants is currently a topic of international debate, and development of robust analytical workflows to identify and communicate clinically relevant variants is paramount. METHODS: The Deciphering Developmental Disorders (DDD) study has developed a UK-wide patient recruitment network involving over 180 clinicians across all 24 regional genetics services, and has performed genome-wide microarray and whole exome sequencing on children with undiagnosed developmental disorders and their parents. After data analysis, pertinent genomic variants were returned to individual research participants via their local clinical genetics team. FINDINGS: Around 80,000 genomic variants were identified from exome sequencing and microarray analysis in each individual, of which on average 400 were rare and predicted to be protein altering. By focusing only on de novo and segregating variants in known developmental disorder genes, we achieved a diagnostic yield of 27% among 1133 previously investigated yet undiagnosed children with developmental disorders, whilst minimising incidental findings. In families with developmentally normal parents, whole exome sequencing of the child and both parents resulted in a 10-fold reduction in the number of potential causal variants that needed clinical evaluation compared to sequencing only the child. Most diagnostic variants identified in known genes were novel and not present in current databases of known disease variation. INTERPRETATION: Implementation of a robust translational genomics workflow is achievable within a large-scale rare disease research study to allow feedback of potentially diagnostic findings to clinicians and research participants. Systematic recording of relevant clinical data, curation of a gene-phenotype knowledge base, and development of clinical decision support software are needed in addition to automated exclusion of almost all variants, which is crucial for scalable prioritisation and review of possible diagnostic variants. However, the resource requirements of development and maintenance of a clinical reporting system within a research setting are substantial. FUNDING: Health Innovation Challenge Fund, a parallel funding partnership between the Wellcome Trust and the UK Department of Health.