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.

Genome Report: chromosome-scale genome assembly of the African spiny mouse (Acomys cahirinus).

There is increasing interest in the African spiny mouse (Acomys cahirinus) as a model organism because of its ability for regeneration of tissue after injury in skin, muscle, and internal organs such as the kidneys. A high-quality reference genome is needed to better understand these regenerative properties at the molecular level. Here, we present an improved reference genome for A. cahirinus generated from long Nanopore sequencing reads. We confirm the quality of our annotations using RNA sequencing data from 4 different tissues. Our genome is of higher contiguity and quality than previously reported genomes from this species and will facilitate ongoing efforts to better understand the regenerative properties of this organism.

GENOME REPORT: Chromosome-scale genome assembly of the African spiny mouse ( Acomys cahirinus ).

There is increasing interest in the African spiny mouse ( Acomys cahirinus ) as a model organism because of its ability for regeneration of tissue after injury in skin, muscle, and internal organs such as the kidneys. A high-quality reference genome is needed to better understand these regenerative properties at the molecular level. Here, we present an improved reference genome for A. cahirinus generated from long Nanopore sequencing reads. We confirm the quality of our annotations using RNA sequencing data from four different tissues. Our genome is of higher contiguity and quality than previously reported genomes from this species and will facilitate ongoing efforts to better understand the regenerative properties of this organism.

Mammalian organ regeneration in spiny mice.

Fibrosis-driven solid organ failure is a major world-wide health burden with few therapeutic options. Spiny mice (genus: Acomys) are terrestrial mammals that regenerate severe skin wounds without fibrotic scars to evade predators. Recent studies have shown that spiny mice also regenerate acute ischemic and traumatic injuries to kidney, heart, spinal cord, and skeletal muscle. A common feature of this evolved wound healing response is a lack of formation of fibrotic scar tissue that degrades organ function, inhibits regeneration, and leads to organ failure. Complex tissue regeneration is an extremely rare property among mammalian species. In this article, we discuss the evidence that Acomys represents an emerging model organism that offers a unique opportunity for the biomedical community to investigate and clinically translate molecular mechanisms of scarless wound healing and regeneration of organ function in a mammalian species.

The arginine methyltransferase Carm1 is necessary for heart development.

To discover genes implicated in human congenital disorders, we performed ENU mutagenesis in the mouse and screened for mutations affecting embryonic development. In this work, we report defects of heart development in mice homozygous for a mutation of coactivator-associated arginine methyltransferase 1 (Carm1). While Carm1 has been extensively studied, it has never been previously associated with a role in heart development. Phenotype analysis combining histology and microcomputed tomography imaging shows a range of cardiac defects. Most notably, many affected midgestation embryos appear to have cardiac rupture and hemorrhaging in the thorax. Mice that survive to late gestation show a variety of cardiac defects, including ventricular septal defects, double outlet right ventricle, and persistent truncus arteriosus. Transcriptome analyses of the mutant embryos by mRNA-seq reveal the perturbation of several genes involved in cardiac morphogenesis and muscle development and function. In addition, we observe the mislocalization of cardiac neural crest cells at E12.5 in the outflow tract. The cardiac phenotype of Carm1 mutant embryos is similar to that of Pax3 null mutants, and PAX3 is a putative target of CARM1. However, our analysis does not support the hypothesis that developmental defects in Carm1 mutant embryos are primarily due to a functional defect of PAX3.

Wound healing and regeneration in spiny mice (Acomys cahirinus).

The winds of Patagonia are referred to by locals as "The Broom of God" because they sweep away the less fit species that cannot survive there. Fitness as an evolutionary trait has been considered as fundamental for many aspects of morphogenesis and behavior in metazoans. Yet, it has not received much attention in the area of wound healing, despite the obvious relevance of this polygenic trait to an organism's survival in nature. In this chapter, we review the evidence that the rodent species Acomys cahirinus is an emerging mammalian model system that has evolved a non-typical (for mammals) wound healing response that offers unique opportunities for the study of organ regeneration without fibrosis in an adult mammalian species.

Spiny mice activate unique transcriptional programs after severe kidney injury regenerating organ function without fibrosis.

Fibrosis-driven solid organ failure is an enormous burden on global health. Spiny mice (Acomys) are terrestrial mammals that can regenerate severe skin wounds without scars to avoid predation. Whether spiny mice also regenerate internal organ injuries is unknown. Here, we show that despite equivalent acute obstructive or ischemic kidney injury, spiny mice fully regenerate nephron structure and organ function without fibrosis, whereas C57Bl/6 or CD1 mice progress to complete organ failure with extensive renal fibrosis. Two mechanisms for vertebrate regeneration have been proposed that emphasize either extrinsic (pro-regenerative macrophages) or intrinsic (surviving cells of the organ itself) controls. Comparative transcriptome analysis revealed that the Acomys genome appears poised at the time of injury to initiate regeneration by surviving kidney cells, whereas macrophage accumulation was not detected until about day 7. Thus, we provide evidence for rapid activation of a gene expression signature for regenerative wound healing in the spiny mouse kidney.

Reelin Mediates Hippocampal Cajal-Retzius Cell Positioning and Infrapyramidal Blade Morphogenesis.

We have previously described hypomorphic reelin (Reln) mutant mice, RelnCTRdel, in which the morphology of the dentate gyrus is distinct from that seen in reeler mice. In the RelnCTRdel mutant, the infrapyramidal blade of the dentate gyrus fails to extend, while the suprapyramidal blade forms with a relatively compact granule neuron layer. Underlying this defect, we now report several developmental anomalies in the RelnCTRdel dentate gyrus. Most strikingly, the distribution of Cajal-Retzius cells was aberrant; Cajal-Retzius neurons were increased in the suprapyramidal blade, but were greatly reduced along the subpial surface of the prospective infrapyramidal blade. We also observed multiple abnormalities of the fimbriodentate junction. Firstly, progenitor cells were distributed abnormally; the "neurogenic cluster" at the fimbriodentate junction was absent, lacking the normal accumulation of Tbr2-positive intermediate progenitors. However, the number of dividing cells in the dentate gyrus was not generally decreased. Secondly, a defect of secondary glial scaffold formation, limited to the infrapyramidal blade, was observed. The densely radiating glial fibers characteristic of the normal fimbriodentate junction were absent in mutants. These fibers might be required for migration of progenitors, which may account for the failure of neurogenic cluster formation. These findings suggest the importance of the secondary scaffold and neurogenic cluster of the fimbriodentate junction in morphogenesis of the mammalian dentate gyrus. Our study provides direct genetic evidence showing that normal RELN function is required for Cajal-Retzius cell positioning in the dentate gyrus, and for formation of the fimbriodentate junction to promote infrapyramidal blade extension.

Differential requirements of tubulin genes in mammalian forebrain development.

Tubulin genes encode a series of homologous proteins used to construct microtubules which are essential for multiple cellular processes. Neural development is particularly reliant on functional microtubule structures. Tubulin genes comprise a large family of genes with very high sequence similarity between multiple family members. Human genetics has demonstrated that a large spectrum of cortical malformations are associated with de novo heterozygous mutations in tubulin genes. However, the absolute requirement for many of these genes in development and disease has not been previously tested in genetic loss of function models. Here we directly test the requirement for Tuba1a, Tubb2a and Tubb2b in the mouse by deleting each gene individually using CRISPR-Cas9 genome editing. We show that loss of Tubb2a or Tubb2b does not impair survival but does lead to relatively mild cortical malformation phenotypes. In contrast, loss of Tuba1a is perinatal lethal and leads to significant forebrain dysmorphology. We also present a novel mouse ENU allele of Tuba1a with phenotypes similar to the null allele. This demonstrates the requirements for each of the tubulin genes and levels of functional redundancy are quite different throughout the gene family. The ability of the mouse to survive in the absence of some tubulin genes known to cause disease in humans suggests future intervention strategies for these devastating tubulinopathy diseases.

Loss of function of Colgalt1 disrupts collagen post-translational modification and causes musculoskeletal defects.

In a screen for organogenesis defects in N-ethyl-N-nitrosourea (ENU)-induced mutant mice, we discovered a line carrying a mutation in Colgalt1 [collagen beta(1-O)galactosyltransferase type 1], which is required for proper galactosylation of hydroxylysine residues in a number of collagens. Colgalt1 mutant embryos have not been previously characterized; here, we show that they exhibit skeletal and muscular defects. Analysis of mutant-derived embryonic fibroblasts reveals that COLGALT1 acts on collagen IV and VI, and, while collagen VI appears stable and its secretion is not affected, collagen IV accumulates inside of cells and within the extracellular matrix, possibly due to instability and increased degradation. We also generated mutant zebrafish that do not express the duplicated orthologs of mammalian Colgalt1 The double-homozygote mutants have muscle defects; they are viable through the larvae stage but do not survive to 10 days post-fertilization. We hypothesize that the Colgalt1 mutant could serve as a model of a human connective tissue disorder and/or congenital muscular dystrophy or myopathy.

Applicability of the Mutation-Selection Balance Model to Population Genetics of Heterozygous Protein-Truncating Variants in Humans.

The fate of alleles in the human population is believed to be highly affected by the stochastic force of genetic drift. Estimation of the strength of natural selection in humans generally necessitates a careful modeling of drift including complex effects of the population history and structure. Protein-truncating variants (PTVs) are expected to evolve under strong purifying selection and to have a relatively high per-gene mutation rate. Thus, it is appealing to model the population genetics of PTVs under a simple deterministic mutation-selection balance, as has been proposed earlier (Cassa et al. 2017). Here, we investigated the limits of this approximation using both computer simulations and data-driven approaches. Our simulations rely on a model of demographic history estimated from 33,370 individual exomes of the Non-Finnish European subset of the ExAC data set (Lek et al. 2016). Additionally, we compared the African and European subset of the ExAC study and analyzed de novo PTVs. We show that the mutation-selection balance model is applicable to the majority of human genes, but not to genes under the weakest selection.

Optimizing Genomic Methods for Mapping and Identification of Candidate Variants in ENU Mutagenesis Screens Using Inbred Mice.

Positional cloning of ENU-induced mutations has traditionally relied on analysis of polymorphic variation between two strains. In contrast, the application of whole-genome sequencing (WGS) has enabled gene discovery in mutant lines maintained on an inbred genetic background. This approach utilizes genetic variation derived from ENU-induced variants for mapping and reduces the likelihood of phenotypic variation, making it an ideal method for genetic modifier screening. Here, we describe the results of such a screen, wherein we determined the minimal number of mutant genomic DNA samples to include in our analyses and improved the sensitivity of our screen by individually barcoding each genomic DNA library. We present several unique cases to illustrate this approach's efficacy, including the discovery of two distinct mutations that generate essentially identical mutant phenotypes, the ascertainment of a non-ENU-induced candidate variant through homozygosity mapping, and an approach for the identification of putative dominant genetic modifiers.

Estimating the selective effects of heterozygous protein-truncating variants from human exome data.

The evolutionary cost of gene loss is a central question in genetics and has been investigated in model organisms and human cell lines. In humans, tolerance of the loss of one or both functional copies of a gene is related to the gene's causal role in disease. However, estimates of the selection and dominance coefficients in humans have been elusive. Here we analyze exome sequence data from 60,706 individuals to make genome-wide estimates of selection against heterozygous loss of gene function. Using this distribution of selection coefficients for heterozygous protein-truncating variants (PTVs), we provide corresponding Bayesian estimates for individual genes. We find that genes under the strongest selection are enriched in embryonic lethal mouse knockouts, Mendelian disease-associated genes, and regulators of transcription. Screening by essentiality, we find a large set of genes under strong selection that are likely to have crucial functions but have not yet been thoroughly characterized.

Generation of a multipurpose Prdm16 mouse allele by targeted gene trapping.

Gene trap mutagenesis is a powerful tool to create loss-of-function mutations in mice and other model organisms. Modifications of traditional gene trap cassettes, including addition of conditional features in the form of Flip-excision (FlEx) arrays to enable directional gene trap cassette inversions by Cre and Flpe site-specific recombinases, greatly enhanced their experimental potential. By taking advantage of these conditional gene trap cassettes, we developed a generic strategy for generating conditional mutations and validated this strategy in mice carrying a multipurpose allele of the Prdm16 transcription factor gene. We demonstrate that the gene trap insertion creates a null mutation replicating the Pierre Robin sequence-type cleft palate phenotype of other Prdm16 mutant mice. Consecutive breeding to Flpe and Emx1IREScre deleter mice spatially restricted Prdm16 loss to regions of the forebrain expressing the homeobox gene Emx1, demonstrating the utility of the technology for the analysis of tissue-specific gene functions.

C-Terminal Region Truncation of RELN Disrupts an Interaction with VLDLR, Causing Abnormal Development of the Cerebral Cortex and Hippocampus.

We discovered a hypomorphic reelin (Reln) mutant with abnormal cortical lamination and no cerebellar hypoplasia. This mutant, RelnCTRdel, carries a chemically induced splice-site mutation that truncates the C-terminal region (CTR) domain of RELN protein and displays remarkably distinct phenotypes from reeler The mutant does not have an inverted cortex, but cortical neurons overmigrate and invade the marginal zone, which are characteristics similar to a phenotype seen in the cerebral cortex of Vldlrnull mice. The dentate gyrus shows a novel phenotype: the infrapyramidal blade is absent, while the suprapyramidal blade is present and laminated. Genetic epistasis analysis showed that RelnCTRdel/Apoer2null double homozygotes have phenotypes akin to those of reeler mutants, while RelnCTRdel/Vldlrnull mice do not. Given that the receptor double knock-out mice resemble reeler mutants, we infer that RelnCTRdel/Apoer2null double homozygotes have both receptor pathways disrupted. This suggests that CTR-truncation disrupts an interaction with VLDLR (very low-density lipoprotein receptor), while the APOER2 signaling pathway remains active, which accounts for the hypomorphic phenotype in RelnCTRdel mice. A RELN-binding assay confirms that CTR truncation significantly decreases RELN binding to VLDLR, but not to APOER2. Together, the in vitro and in vivo results demonstrate that the CTR domain confers receptor-binding specificity of RELN. SIGNIFICANCE STATEMENT: Reelin signaling is important for brain development and is associated with human type II lissencephaly. Reln mutations in mice and humans are usually associated with cerebellar hypoplasia. A new Reln mutant with a truncation of the C-terminal region (CTR) domain shows that Reln mutation can cause abnormal phenotypes in the cortex and hippocampus without cerebellar hypoplasia. Genetic analysis suggested that CTR truncation disrupts an interaction with the RELN receptor VLDLR (very low-density lipoprotein receptor); this was confirmed by a RELN-binding assay. This result provides a mechanistic explanation for the hypomorphic phenotype of the CTR-deletion mutant, and further suggests that Reln mutations may cause more subtle forms of human brain malformation than classic lissencephalies.

Dysfunction of intraflagellar transport-A causes hyperphagia-induced obesity and metabolic syndrome.

Primary cilia extend from the plasma membrane of most vertebrate cells and mediate signaling pathways. Ciliary dysfunction underlies ciliopathies, which are genetic syndromes that manifest multiple clinical features, including renal cystic disease and obesity. THM1 (also termed TTC21B or IFT139) encodes a component of the intraflagellar transport-A complex and mutations in THM1 have been identified in 5% of individuals with ciliopathies. Consistent with this, deletion of murine Thm1 during late embryonic development results in cystic kidney disease. Here, we report that deletion of murine Thm1 during adulthood results in obesity, diabetes, hypertension and fatty liver disease, with gender differences in susceptibility to weight gain and metabolic dysfunction. Pair-feeding of Thm1 conditional knock-out mice relative to control littermates prevented the obesity and related disorders, indicating that hyperphagia caused the obese phenotype. Thm1 ablation resulted in increased localization of adenylyl cyclase III in primary cilia that were shortened, with bulbous distal tips on neurons of the hypothalamic arcuate nucleus, an integrative center for signals that regulate feeding and activity. In pre-obese Thm1 conditional knock-out mice, expression of anorexogenic pro-opiomelanocortin (Pomc) was decreased by 50% in the arcuate nucleus, which likely caused the hyperphagia. Fasting of Thm1 conditional knock-out mice did not alter Pomc nor orexogenic agouti-related neuropeptide (Agrp) expression, suggesting impaired sensing of changes in peripheral signals. Together, these data indicate that the Thm1-mutant ciliary defect diminishes sensitivity to feeding signals, which alters appetite regulation and leads to hyperphagia, obesity and metabolic disease.

Improvement of ENU Mutagenesis Efficiency Using Serial Injection and Mismatch Repair Deficiency Mice.

ENU mutagenesis is a powerful method for generating novel lines of mice that are informative with respect to both fundamental biological processes and human disease. Rapid developments in genomic technology have made the task of identifying causal mutations by positional cloning remarkably efficient. One limitation of this approach remains the mutation frequency achievable using standard treatment protocols, which currently generate approximately 1-2 sequence changes per megabase when optimized. In this study we used two strategies to attempt to increase the number of mutations induced by ENU treatment. One approach employed mice carrying a mutation in the DNA repair enzyme Msh6. The second strategy involved injection of ENU to successive generations of mice. To evaluate the number of ENU-induced mutations, single mice or pooled samples were analyzed using whole exome sequencing. The results showed that there is considerable variability in the induced mutation frequency using these approaches, but an overall increase in ENU-induced variants from one generation to another was observed. The analysis of the mice deficient for Msh6 also showed an increase in the ENU-induced variants compared to the wild-type ENU-treated mice. However, in both cases the increase in ENU-induced mutation frequency was modest.

Mutations in Dnaaf1 and Lrrc48 Cause Hydrocephalus, Laterality Defects, and Sinusitis in Mice.

We have previously described a forward genetic screen in mice for abnormalities of brain development. Characterization of two hydrocephalus mutants by whole-exome sequencing after whole-genome SNP mapping revealed novel recessive mutations in Dnaaf1 and Lrrc48 Mouse mutants of these two genes have not been previously reported. The Dnaaf1 mutant carries a mutation at the splice donor site of exon 4, which results in abnormal transcripts. The Lrrc48 mutation is a missense mutation at a highly conserved leucine residue, which is also associated with a decrease in Lrrc48 transcription. Both Dnaaf1 and Lrrc48 belong to a leucine-rich repeat-containing protein family and are components of the ciliary axoneme. Their Chlamydomonas orthologs are known to be required for normal ciliary beat frequency or flagellar waveform, respectively. Some Dnaaf1 or Lrrc48 homozygote mutants displayed laterality defects, suggesting a motile cilia defect in the embryonic node. Mucus accumulation and neutrophil infiltration in the maxillary sinuses suggested sinusitis. Dnaaf1 mutants showed postnatal lethality, and none survived to weaning age. Lrrc48 mutants survive to adulthood, but had male infertility. ARL13B immunostaining showed the presence of motile cilia in the mutants, and the distal distribution of DNAH9 in the axoneme of upper airway motile cilia appeared normal. The phenotypic abnormalities suggest that mutations in Dnaaf1 and Lrrc48 cause defects in motile cilia function.

Reduction of ciliary length through pharmacologic or genetic inhibition of CDK5 attenuates polycystic kidney disease in a model of nephronophthisis.

Polycystic kidney diseases (PKDs) comprise a subgroup of ciliopathies characterized by the formation of fluid-filled kidney cysts and progression to end-stage renal disease. A mechanistic understanding of cystogenesis is crucial for the development of viable therapeutic options. Here, we identify CDK5, a kinase active in post mitotic cells, as a new and important mediator of PKD progression. We show that long-lasting attenuation of PKD in the juvenile cystic kidneys (jck) mouse model of nephronophthisis by pharmacological inhibition of CDK5 using either R-roscovitine or S-CR8 is accompanied by sustained shortening of cilia and a more normal epithelial phenotype, suggesting this treatment results in a reprogramming of cellular differentiation. Also, a knock down of Cdk5 in jck cells using small interfering RNA results in significant shortening of ciliary length, similar to what we observed with R-roscovitine. Finally, conditional inactivation of Cdk5 in the jck mice significantly attenuates cystic disease progression and is associated with shortening of ciliary length as well as restoration of cellular differentiation. Our results suggest that CDK5 may regulate ciliary length by affecting tubulin dynamics via its substrate collapsin response mediator protein 2. Taken together, our data support therapeutic approaches aimed at restoration of ciliogenesis and cellular differentiation as a promising strategy for the treatment of renal cystic diseases.