Specific mammalian genes functionally and dynamically associate together within the nucleus. Yet, how an array of many genes along the chromosome sequence can be spatially organized and folded together is unknown. We investigated the 3D structure of a well-annotated, highly conserved 4.3-Mb region on mouse chromosome 14 that contains four clusters of genes separated by gene "deserts." In nuclei, this region forms multiple, nonrandom "higher order" structures. These structures are based on the gene distribution pattern in primary sequence and are marked by preferential associations among multiple gene clusters. Associating gene clusters represent expressed chromatin, but their aggregation is not simply dependent on ongoing transcription. In chromosomes with aggregated gene clusters, gene deserts preferentially align with the nuclear periphery, providing evidence for chromosomal region architecture by specific associations with functional nuclear domains. Together, these data suggest dynamic, probabilistic 3D folding states for a contiguous megabase-scale chromosomal region, supporting the diverse activities of multiple genes and their conserved primary sequence organization.

The generation of a comprehensive catalog of null alleles covering all protein-coding genes is the goal of the International Mouse Phenotyping Consortium. Over the past 20 years, significant progress has been made towards achieving this goal through the combined efforts of many large-scale programs that built an embryonic stem cell resource to generate knockout mice and more recently employed CRISPR/Cas9-based mutagenesis to delete critical regions predicted to result in frameshift mutations, thus, ablating gene function. The IMPC initiative builds on prior and ongoing work by individual research groups creating gene knockouts in the mouse. Here, we analyze the collective efforts focusing on the combined null allele resource resulting from strains developed by the research community and large-scale production programs. Based upon this pooled analysis, we examine the remaining fraction of protein-coding genes focusing on clearly defined mouse-human orthologs as the highest priority for completing the mutant mouse null resource. In summary, we find that there are less than 3400 mouse-human orthologs remaining in the genome without a targeted null allele that can be further prioritized to achieve our overall goal of the complete functional annotation of the protein-coding portion of a mammalian genome.

The purpose of this study was to determine whether implementation of a multicomponent organizational intervention can produce significant change in diabetes care and outcomes in community primary care practices.

Primary cilia are nearly ubiquitous organelles that transduce molecular and mechanical signals. Although the basic structure of the cilium and the cadre of genes that contribute to ciliary formation and function (the ciliome) are believed to be evolutionarily conserved, the presentation of ciliopathies with narrow, tissue-specific phenotypes and distinct molecular readouts suggests that an unappreciated heterogeneity exists within this organelle. Here, we provide a searchable transcriptomic resource for a curated primary ciliome, detailing various subgroups of differentially expressed genes within the ciliome that display tissue and temporal specificity. Genes within the differentially expressed ciliome exhibited a lower level of functional constraint across species, suggesting organism and cell-specific function adaptation. The biological relevance of ciliary heterogeneity was functionally validated by using Cas9 gene-editing to disrupt ciliary genes that displayed dynamic gene expression profiles during osteogenic differentiation of multipotent neural crest cells. Collectively, this novel primary cilia-focused resource will allow researchers to explore longstanding questions related to how tissue and cell-type specific functions and ciliary heterogeneity may contribute to the range of phenotypes associated with ciliopathies.

Approximately one-third of all mammalian genes are essential for life. Phenotypes resulting from knockouts of these genes in mice have provided tremendous insight into gene function and congenital disorders. As part of the International Mouse Phenotyping Consortium effort to generate and phenotypically characterize 5,000 knockout mouse lines, here we identify 410 lethal genes during the production of the first 1,751 unique gene knockouts. Using a standardized phenotyping platform that incorporates high-resolution 3D imaging, we identify phenotypes at multiple time points for previously uncharacterized genes and additional phenotypes for genes with previously reported mutant phenotypes. Unexpectedly, our analysis reveals that incomplete penetrance and variable expressivity are common even on a defined genetic background. In addition, we show that human disease genes are enriched for essential genes, thus providing a dataset that facilitates the prioritization and validation of mutations identified in clinical sequencing efforts.

Analysis of large-scale human genomic data has yielded unexplained mutations known to cause severe disease in healthy individuals. Here, we report the unexpected recovery of a rare dominant lethal mutation in TPM1, a sarcomeric actin-binding protein, in eight individuals with large atrial septal defect (ASD) in a five-generation pedigree. Mice with Tpm1 mutation exhibit early embryonic lethality with disrupted myofibril assembly and no heartbeat. However, patient-induced pluripotent-stem-cell-derived cardiomyocytes show normal beating with mild myofilament defect, indicating disease suppression. A variant in TLN2, another myofilament actin-binding protein, is identified as a candidate suppressor. Mouse CRISPR knock-in (KI) of both the TLN2 and TPM1 variants rescues heart beating, with near-term fetuses exhibiting large ASD. Thus, the role of TPM1 in ASD pathogenesis unfolds with suppression of its embryonic lethality by protective TLN2 variant. These findings provide evidence that genetic resiliency can arise with genetic suppression of a deleterious mutation.

Despite a common understanding that Gli TFs are utilized to convey a Hh morphogen gradient, genetic analyses suggest craniofacial development does not completely fit this paradigm. Using the mouse model (Mus musculus), we demonstrated that rather than being driven by a Hh threshold, robust Gli3 transcriptional activity during skeletal and glossal development required interaction with the basic helix-loop-helix TF Hand2. Not only did genetic and expression data support a co-factorial relationship, but genomic analysis revealed that Gli3 and Hand2 were enriched at regulatory elements for genes essential for mandibular patterning and development. Interestingly, motif analysis at sites co-occupied by Gli3 and Hand2 uncovered mandibular-specific, low-affinity, 'divergent' Gli-binding motifs (dGBMs). Functional validation revealed these dGBMs conveyed synergistic activation of Gli targets essential for mandibular patterning and development. In summary, this work elucidates a novel, sequence-dependent mechanism for Gli transcriptional activity within the craniofacial complex that is independent of a graded Hh signal.

Genome editing with CRISPR-associated (Cas) proteins holds exceptional promise for "correcting" variants causing genetic disease. To realize this promise, off-target genomic changes cannot occur during the editing process. Here, we use whole genome sequencing to compare the genomes of 50 Cas9-edited founder mice to 28 untreated control mice to assess the occurrence of S. pyogenes Cas9-induced off-target mutagenesis. Computational analysis of whole-genome sequencing data detects 26 unique sequence variants at 23 predicted off-target sites for 18/163 guides used. While computationally detected variants are identified in 30% (15/50) of Cas9 gene-edited founder animals, only 38% (10/26) of the variants in 8/15 founders validate by Sanger sequencing. In vitro assays for Cas9 off-target activity identify only two unpredicted off-target sites present in genome sequencing data. In total, only 4.9% (8/163) of guides tested have detectable off-target activity, a rate of 0.2 Cas9 off-target mutations per founder analyzed. In comparison, we observe ~1,100 unique variants in each mouse regardless of genome exposure to Cas9 indicating off-target variants comprise a small fraction of genetic heterogeneity in Cas9-edited mice. These findings will inform future design and use of Cas9-edited animal models as well as provide context for evaluating off-target potential in genetically diverse patient populations.

Genome editing using the CRISPR/Cas9 RNA-guided endonuclease system has rapidly become a driving force for discovery in modern biomedical research. This simple yet elegant system has been widely used to generate both loss-of-function alleles and precision knock-in mutations using single-stranded donor oligonucleotides. Our CRISPRtools platform supports both of these applications in order to facilitate the use of CRISPR/Cas9. While there are several tools that facilitate CRISPR/Cas9 design and screen for potential off-target sites, the process is typically performed sequentially on single genes, limiting scalability for large-scale programs. Here, the design principle underlying gene ablation is based upon using paired guides flanking a critical region/exon of interest to create deletions. Guide pairs are rank ordered based upon published efficiency scores and off-target analyses, and reported in a concise format for downstream implementation. The exon deletion strategy simplifies characterization of founder animals and is the strategy employed for the majority of knockouts in the mouse. In proof-of-principle experiments, the effectiveness of this approach is demonstrated using microinjection and electroporation to introduce CRISPR/Cas9 components into mouse zygotes to delete critical exons.

Heterotaxy results from a failure to establish normal left-right asymmetry early in embryonic development. By whole-exome sequencing, whole-genome sequencing and high-throughput cohort resequencing, we identified recessive mutations in MMP21 (encoding matrix metallopeptidase 21) in nine index cases with heterotaxy. In addition, Mmp21-mutant mice and mmp21-morphant zebrafish displayed heterotaxy and abnormal cardiac looping, respectively, suggesting a new role for extracellular matrix remodeling in the establishment of laterality in vertebrates.

Congenital heart disease (CHD) affects up to 1% of live births. Although a genetic etiology is indicated by an increased recurrence risk, sporadic occurrence suggests that CHD genetics is complex. Here, we show that hypoplastic left heart syndrome (HLHS), a severe CHD, is multigenic and genetically heterogeneous. Using mouse forward genetics, we report what is, to our knowledge, the first isolation of HLHS mutant mice and identification of genes causing HLHS. Mutations from seven HLHS mouse lines showed multigenic enrichment in ten human chromosome regions linked to HLHS. Mutations in Sap130 and Pcdha9, genes not previously associated with CHD, were validated by CRISPR-Cas9 genome editing in mice as being digenic causes of HLHS. We also identified one subject with HLHS with SAP130 and PCDHA13 mutations. Mouse and zebrafish modeling showed that Sap130 mediates left ventricular hypoplasia, whereas Pcdha9 increases penetrance of aortic valve abnormalities, both signature HLHS defects. These findings show that HLHS can arise genetically in a combinatorial fashion, thus providing a new paradigm for the complex genetics of CHD.

The identification of causal variants in sequencing studies remains a considerable challenge that can be partially addressed by new gene-specific knowledge. Here, we integrate measures of how essential a gene is to supporting life, as inferred from viability and phenotyping screens performed on knockout mice by the International Mouse Phenotyping Consortium and essentiality screens carried out on human cell lines. We propose a cross-species gene classification across the Full Spectrum of Intolerance to Loss-of-function (FUSIL) and demonstrate that genes in five mutually exclusive FUSIL categories have differing biological properties. Most notably, Mendelian disease genes, particularly those associated with developmental disorders, are highly overrepresented among genes non-essential for cell survival but required for organism development. After screening developmental disorder cases from three independent disease sequencing consortia, we identify potentially pathogenic variants in genes not previously associated with rare diseases. We therefore propose FUSIL as an efficient approach for disease gene discovery.

Precise cis-regulatory control of gene expression is essential for normal embryogenesis and tissue development. The BMP antagonist Gremlin1 (Grem1) is a key node in the signalling system that coordinately controls limb bud development. Here, we use mouse reverse genetics to identify the enhancers in the Grem1 genomic landscape and the underlying cis-regulatory logics that orchestrate the spatio-temporal Grem1 expression dynamics during limb bud development. We establish that transcript levels are controlled in an additive manner while spatial regulation requires synergistic interactions among multiple enhancers. Disrupting these interactions shows that altered spatial regulation rather than reduced Grem1 transcript levels prefigures digit fusions and loss. Two of the enhancers are evolutionary ancient and highly conserved from basal fishes to mammals. Analysing these enhancers from different species reveal the substantial spatial plasticity in Grem1 regulation in tetrapods and basal fishes, which provides insights into the fin-to-limb transition and evolutionary diversification of pentadactyl limbs.

Advancing our understanding of embryonic development is heavily dependent on identification of novel pathways or regulators. Although genome-wide techniques such as RNA sequencing are ideally suited for discovering novel candidate genes, they are unable to yield spatially resolved information in embryos or tissues. Microscopy-based approaches, using in situ hybridization, for example, can provide spatial information about gene expression, but are limited to analyzing one or a few genes at a time. Here, we present a method where we combine traditional histological techniques with low-input RNA sequencing and mathematical image reconstruction to generate a high-resolution genome-wide 3D atlas of gene expression in the zebrafish embryo at three developmental stages. Importantly, our technique enables searching for genes that are expressed in specific spatial patterns without manual image annotation. We envision broad applicability of RNA tomography as an accurate and sensitive approach for spatially resolved transcriptomics in whole embryos and dissected organs.

Hedgehog signaling controls pattern formation in many vertebrate tissues. The downstream effectors of the pathway are the bifunctional Gli transcription factors, which, depending on hedgehog concentration, act as either transcriptional activators or repressors. Quantitatively understanding the interplay between Gli activator and repressor forms for patterning complex tissues is an open challenge. Here, we describe a reductionist mathematical model for how Gli activators and repressors are integrated in space and time to regulate transcriptional outputs of hedgehog signaling, using the pathway readouts Gli1 and Ptch1 as a model system. Spatially resolved measurements of absolute transcript numbers for these genes allow us to infer spatiotemporal variations of Gli activator and repressor levels. We validate our model by successfully predicting expression changes of Gli1 and Ptch1 in mutants at different developmental stages and in different tissues. Our results provide a starting point for understanding gene regulation by bifunctional transcription factors during mammalian development.

The Hedgehog signaling pathway is part of the ancient developmental-evolutionary animal toolkit. Frequently co-opted to pattern new structures, the pathway is conserved among eumetazoans yet flexible and pleiotropic in its effects. The Hedgehog receptor, Patched, is transcriptionally activated by Hedgehog, providing essential negative feedback in all tissues. Our locus-wide dissections of the cis-regulatory landscapes of fly patched and mouse Ptch1 reveal abundant, diverse enhancers with stage- and tissue-specific expression patterns. The seemingly simple, constitutive Hedgehog response of patched/Ptch1 is driven by a complex regulatory architecture, with batteries of context-specific enhancers engaged in promoter-specific interactions to tune signaling individually in each tissue, without disturbing patterning elsewhere. This structure-one of the oldest cis-regulatory features discovered in animal genomes-explains how patched/Ptch1 can drive dramatic adaptations in animal morphology while maintaining its essential core function. It may also suggest a general model for the evolutionary flexibility of conserved regulators and pathways.

Nephron progenitor number determines nephron endowment; a reduced nephron count is linked to the onset of kidney disease. Several transcriptional regulators including Six2, Wt1, Osr1, Sall1, Eya1, Pax2, and Hox11 paralogues are required for specification and/or maintenance of nephron progenitors. However, little is known about the regulatory intersection of these players. Here, we have mapped nephron progenitor-specific transcriptional networks of Six2, Hoxd11, Osr1, and Wt1. We identified 373 multi-factor associated 'regulatory hotspots' around genes closely associated with progenitor programs. To examine their functional significance, we deleted 'hotspot' enhancer elements for Six2 and Wnt4. Removal of the distal enhancer for Six2 leads to a ~40% reduction in Six2 expression. When combined with a Six2 null allele, progeny display a premature depletion of nephron progenitors. Loss of the Wnt4 enhancer led to a significant reduction of Wnt4 expression in renal vesicles and a mildly hypoplastic kidney, a phenotype also enhanced in combination with a Wnt4 null mutation. To explore the regulatory landscape that supports proper target gene expression, we performed CTCF ChIP-seq to identify insulator-boundary regions. One such putative boundary lies between the Six2 and Six3 loci. Evidence for the functional significance of this boundary was obtained by deep sequencing of the radiation-induced Brachyrrhine (Br) mutant allele. We identified an inversion of the Six2/Six3 locus around the CTCF-bound boundary, removing Six2 from its distal enhancer regulation, but placed next to Six3 enhancer elements which support ectopic Six2 expression in the lens where Six3 is normally expressed. Six3 is now predicted to fall under control of the Six2 distal enhancer. Consistent with this view, we observed ectopic Six3 in nephron progenitors. 4C-seq supports the model for Six2 distal enhancer interactions in wild-type and Br/+ mouse kidneys. Together, these data expand our view of the regulatory genome and regulatory landscape underpinning mammalian nephrogenesis.