A Comprehensive Atlas of AAV Tropism in the Mouse.

Gene therapy with Adeno-Associated Viral (AAV) vectors requires knowledge of their tropism within the body. Here we analyze the tropism of ten naturally occurring AAV serotypes (AAV3B, AAV4, AAV5, AAV6, AAV7, AAV8, AAV9, AAVrh8, AAVrh10 and AAVrh74) following systemic delivery into male and female mice. A transgene expressing ZsGreen and Cre recombinase was used to identify transduction in a cell-dependent manner based on fluorescence. Cre-driven activation of tdTomato fluorescence offered superior sensitivity for transduced cells. All serotypes except AAV3B and AAV4 had high liver tropism. Fluorescence activation revealed transduction of unexpected tissues, including adrenals, testes and ovaries. Rare transduced cells within tissues were also readily visualized. Biodistribution of AAV genomes correlated with fluorescence, except in immune tissues. AAV4 was found to have a pan-endothelial tropism while also targeting pancreatic beta cells. This public resource enables selection of the best AAV serotypes for basic science and preclinical applications in mice.

Amphiphilic shuttle peptide delivers base editor ribonucleoprotein to correct the CFTR R553X mutation in well-differentiated airway epithelial cells.

Base editing could correct nonsense mutations that cause cystic fibrosis (CF), but clinical development is limited by the lack of delivery methods that efficiently breach the barriers presented by airway epithelia. Here, we present a novel amphiphilic shuttle peptide based on the previously reported S10 peptide that substantially improved base editor ribonucleoprotein (RNP) delivery. Studies of the S10 secondary structure revealed that the alpha-helix formed by the endosomal leakage domain (ELD), but not the cell penetrating peptide (CPP), was functionally important for delivery. By isolating and extending the ELD, we created a novel shuttle peptide, termed S237. While S237 achieved lower delivery of green fluorescent protein, it outperformed S10 at Cas9 RNP delivery to cultured human airway epithelial cells and to pig airway epithelia in vivo, possibly due to its lower net charge. In well-differentiated primary human airway epithelial cell cultures, S237 achieved a 4.6-fold increase in base editor RNP delivery, correcting up to 9.4% of the cystic fibrosis transmembrane conductance regulator (CFTR) R553X allele and restoring CFTR channel function close to non-CF levels. These findings deepen the understanding of peptide-mediated delivery and offer a translational approach for base editor RNP delivery for CF airway disease.

Impact of essential genes on the success of genome editing experiments generating 3313 new genetically engineered mouse lines.

The International Mouse Phenotyping Consortium (IMPC) systematically produces and phenotypes mouse lines with presumptive null mutations to provide insight into gene function. The IMPC now uses the programmable RNA-guided nuclease Cas9 for its increased capacity and flexibility to efficiently generate null alleles in the C57BL/6N strain. In addition to being a valuable novel and accessible research resource, the production of 3313 knockout mouse lines using comparable protocols provides a rich dataset to analyze experimental and biological variables affecting in vivo gene engineering with Cas9. Mouse line production has two critical steps - generation of founders with the desired allele and germline transmission (GLT) of that allele from founders to offspring. A systematic evaluation of the variables impacting success rates identified gene essentiality as the primary factor influencing successful production of null alleles. Collectively, our findings provide best practice recommendations for using Cas9 to generate alleles in mouse essential genes, many of which are orthologs of genes linked to human disease.

A graph theoretical approach to experimental prioritization in genome-scale investigations.

The goal of systems biology is to gain a network level understanding of how gene interactions influence biological states, and ultimately inform upon human disease. Given the scale and scope of systems biology studies, resource constraints often limit researchers when validating genome-wide phenomena and potentially lead to an incomplete understanding of the underlying mechanisms. Further, prioritization strategies are often biased towards known entities (e.g. previously studied genes/proteins with commercially available reagents), and other technical issues that limit experimental breadth. Here, heterogeneous biological information is modeled as an association graph to which a high-performance minimum dominating set solver is applied to maximize coverage across the graph, and thus increase the breadth of experimentation. First, we tested our model on retrieval of existing gene functional annotations and demonstrated that minimum dominating set returns more diverse terms when compared to other computational methods. Next, we utilized our heterogenous network and minimum dominating set solver to assist in the process of identifying understudied genes to be interrogated by the International Mouse Phenotyping Consortium. Using an unbiased algorithmic strategy, poorly studied genes are prioritized from the remaining thousands of genes yet to be characterized. This method is tunable and extensible with the potential to incorporate additional user-defined prioritizing information. The minimum dominating set approach can be applied to any biological network in order to identify a tractable subset of features to test experimentally or to assist in prioritizing candidate genes associated with human disease.

Improving laboratory animal genetic reporting: LAG-R guidelines.

The biomedical research community addresses reproducibility challenges in animal studies through standardized nomenclature, improved experimental design, transparent reporting, data sharing, and centralized repositories. The ARRIVE guidelines outline documentation standards for laboratory animals in experiments, but genetic information is often incomplete. To remedy this, we propose the Laboratory Animal Genetic Reporting (LAG-R) framework. LAG-R aims to document animals' genetic makeup in scientific publications, providing essential details for replication and appropriate model use. While verifying complete genetic compositions may be impractical, better reporting and validation efforts enhance reliability of research. LAG-R standardization will bolster reproducibility, peer review, and overall scientific rigor.

A frameshift mutation in the murine Prkra gene causes dystonia and exhibits abnormal cerebellar development and reduced eIF2α phosphorylation.

Mutations in Prkra gene, which encodes PACT/RAX cause early onset primary dystonia DYT-PRKRA, a movement disorder that disrupts coordinated muscle movements. PACT/RAX activates protein kinase R (PKR, aka EIF2AK2) by a direct interaction in response to cellular stressors to mediate phosphorylation of the α subunit of the eukaryotic translation initiation factor 2 (eIF2α). Mice homozygous for a naturally arisen, recessively inherited frameshift mutation, Prkra lear-5J exhibit progressive dystonia. In the present study, we investigate the biochemical and developmental consequences of the Prkra lear-5J mutation. Our results indicate that the truncated PACT/RAX protein retains its ability to interact with PKR, however, it inhibits PKR activation. Furthermore, mice homozygous for the mutation have abnormalities in the cerebellar development as well as a severe lack of dendritic arborization of Purkinje neurons. Additionally, reduced eIF2α phosphorylation is noted in the cerebellums and Purkinje neurons of the homozygous Prkra lear-5J mice. These results indicate that PACT/RAX mediated regulation of PKR activity and eIF2α phosphorylation plays a role in cerebellar development and contributes to the dystonia phenotype resulting from this mutation.

Caspase-11 drives macrophage hyperinflammation in models of Polg-related mitochondrial disease.

Mitochondrial diseases (MtD) represent a significant public health challenge due to their heterogenous clinical presentation, often severe and progressive symptoms, and the lack of effective therapies. Environmental exposures, such bacterial and viral infection, can further compromise mitochondrial function and exacerbate the progression of MtD. Infections in MtD patients more frequently progress to sepsis, pneumonia, and other detrimental inflammatory endpoints. However, the underlying immune alterations that enhance immunopathology in MtD remain unclear, constituting a key gap in knowledge that complicates treatment and increases mortality in this population. Here we employ in vitro and in vivo approaches to clarify the molecular and cellular basis for innate immune hyperactivity in models of polymerase gamma (Polg)-related MtD. We reveal that type I interferon (IFN-I)-mediated upregulation of caspase-11 and guanylate-binding proteins (GBPs) increase macrophage sensing of the opportunistic microbe Pseudomonas aeruginosa (PA) in Polg mutant mice. Furthermore, we show that excessive macrophage cytokine secretion and pyroptotic cell death contribute to lung inflammation and morbidity after infection with PA. Our work sheds new light on innate immune dysregulation in MtD and reveals potential targets for limiting infection- and inflammation-related complications in Polg-related MtD.

Lethal phenotypes in Mendelian disorders.

PURPOSE: Existing resources that characterize the essentiality status of genes are based on either proliferation assessment in human cell lines, viability evaluation in mouse knockouts, or constraint metrics derived from human population sequencing studies. Several repositories document phenotypic annotations for rare disorders; however, there is a lack of comprehensive reporting on lethal phenotypes. METHODS: We queried Online Mendelian Inheritance in Man for terms related to lethality and classified all Mendelian genes according to the earliest age of death recorded for the associated disorders, from prenatal death to no reports of premature death. We characterized the genes across these lethality categories, examined the evidence on viability from mouse models and explored how this information could be used for novel gene discovery. RESULTS: We developed the Lethal Phenotypes Portal to showcase this curated catalog of human essential genes. Differences in the mode of inheritance, physiological systems affected, and disease class were found for genes in different lethality categories, as well as discrepancies between the lethal phenotypes observed in mouse and human. CONCLUSION: We anticipate that this resource will aid clinicians in the diagnosis of early lethal conditions and assist researchers in investigating the properties that make these genes essential for human development.

A single-cell time-lapse of mouse prenatal development from gastrula to birth.

The house mouse (Mus musculus) is an exceptional model system, combining genetic tractability with close evolutionary affinity to humans1,2. Mouse gestation lasts only 3 weeks, during which the genome orchestrates the astonishing transformation of a single-cell zygote into a free-living pup composed of more than 500 million cells. Here, to establish a global framework for exploring mammalian development, we applied optimized single-cell combinatorial indexing3 to profile the transcriptional states of 12.4 million nuclei from 83 embryos, precisely staged at 2- to 6-hour intervals spanning late gastrulation (embryonic day 8) to birth (postnatal day 0). From these data, we annotate hundreds of cell types and explore the ontogenesis of the posterior embryo during somitogenesis and of kidney, mesenchyme, retina and early neurons. We leverage the temporal resolution and sampling depth of these whole-embryo snapshots, together with published data4-8 from earlier timepoints, to construct a rooted tree of cell-type relationships that spans the entirety of prenatal development, from zygote to birth. Throughout this tree, we systematically nominate genes encoding transcription factors and other proteins as candidate drivers of the in vivo differentiation of hundreds of cell types. Remarkably, the most marked temporal shifts in cell states are observed within one hour of birth and presumably underlie the massive physiological adaptations that must accompany the successful transition of a mammalian fetus to life outside the womb.

Discovery and validation of genes driving drug-intake and related behavioral traits in mice.

Substance use disorders are heritable disorders characterized by compulsive drug use, the biological mechanisms for which remain largely unknown. Genetic correlations reveal that predisposing drug-naïve phenotypes, including anxiety, depression, novelty preference and sensation seeking, are predictive of drug-use phenotypes, thereby implicating shared genetic mechanisms. High-throughput behavioral screening in knockout (KO) mice allows efficient discovery of the function of genes. We used this strategy in two rounds of candidate prioritization in which we identified 33 drug-use candidate genes based upon predisposing drug-naïve phenotypes and ultimately validated the perturbation of 22 genes as causal drivers of substance intake. We selected 19/221 KO strains (8.5%) that had a difference from control on at least one drug-naïve predictive behavioral phenotype and determined that 15/19 (~80%) affected the consumption or preference for alcohol, methamphetamine or both. No mutant exhibited a difference in nicotine consumption or preference which was possibly confounded with saccharin. In the second round of prioritization, we employed a multivariate approach to identify outliers and performed validation using methamphetamine two-bottle choice and ethanol drinking-in-the-dark protocols. We identified 15/401 KO strains (3.7%, which included one gene from the first cohort) that differed most from controls for the predisposing phenotypes. 8 of 15 gene deletions (53%) affected intake or preference for alcohol, methamphetamine or both. Using multivariate and bioinformatic analyses, we observed multiple relations between predisposing behaviors and drug intake, revealing many distinct biobehavioral processes underlying these relationships. The set of mouse models identified in this study can be used to characterize these addiction-related processes further.

Lethal phenotypes in Mendelian disorders.

Essential genes are those whose function is required for cell proliferation and/or organism survival. A gene's intolerance to loss-of-function can be allocated within a spectrum, as opposed to being considered a binary feature, since this function might be essential at different stages of development, genetic backgrounds or other contexts. Existing resources that collect and characterise the essentiality status of genes are based on either proliferation assessment in human cell lines, embryonic and postnatal viability evaluation in different model organisms, and gene metrics such as intolerance to variation scores derived from human population sequencing studies. There are also several repositories available that document phenotypic annotations for rare disorders in humans such as the Online Mendelian Inheritance in Man (OMIM) and the Human Phenotype Ontology (HPO) knowledgebases. This raises the prospect of being able to use clinical data, including lethality as the most severe phenotypic manifestation, to further our characterisation of gene essentiality. Here we queried OMIM for terms related to lethality and classified all Mendelian genes into categories, according to the earliest age of death recorded for the associated disorders, from prenatal death to no reports of premature death. To showcase this curated catalogue of human essential genes, we developed the Lethal Phenotypes Portal (https://lethalphenotypes.research.its.qmul.ac.uk), where we also explore the relationships between these lethality categories, constraint metrics and viability in cell lines and mouse. Further analysis of the genes in these categories reveals differences in the mode of inheritance of the associated disorders, physiological systems affected and disease class. We highlight how the phenotypic similarity between genes in the same lethality category combined with gene family/group information can be used for novel disease gene discovery. Finally, we explore the overlaps and discrepancies between the lethal phenotypes observed in mouse and human and discuss potential explanations that include differences in transcriptional regulation, functional compensation and molecular disease mechanisms. We anticipate that this resource will aid clinicians in the diagnosis of early lethal conditions and assist researchers in investigating the properties that make these genes essential for human development.

Focused ultrasound-mediated brain genome editing.

Gene editing in the brain has been challenging because of the restricted transport imposed by the blood-brain barrier (BBB). Current approaches mainly rely on local injection to bypass the BBB. However, such administration is highly invasive and not amenable to treating certain delicate regions of the brain. We demonstrate a safe and effective gene editing technique by using focused ultrasound (FUS) to transiently open the BBB for the transport of intravenously delivered CRISPR/Cas9 machinery to the brain.

DISCOVERY AND VALIDATION OF GENES DRIVING DRUG-INTAKE AND RELATED BEHAVIORAL TRAITS IN MICE.

Substance use disorders (SUDs) are heritable disorders characterized by compulsive drug use, but the biological mechanisms driving addiction remain largely unknown. Genetic correlations reveal that predisposing drug-naïve phenotypes, including anxiety, depression, novelty preference, and sensation seeking, are predictive of drug-use phenotypes, implicating shared genetic mechanisms. Because of this relationship, high-throughput behavioral screening of predictive phenotypes in knockout (KO) mice allows efficient discovery of genes likely to be involved in drug use. We used this strategy in two rounds of screening in which we identified 33 drug-use candidate genes and ultimately validated the perturbation of 22 of these genes as causal drivers of substance intake. In our initial round of screening, we employed the two-bottle-choice paradigms to assess alcohol, methamphetamine, and nicotine intake. We identified 19 KO strains that were extreme responders on at least one predictive phenotype. Thirteen of the 19 gene deletions (68%) significantly affected alcohol use three methamphetamine use, and two both. In the second round of screening, we employed a multivariate approach to identify outliers and performed validation using methamphetamine two-bottle choice and ethanol drinking-in-the-dark protocols. We identified 15 KO strains that were extreme responders across the predisposing drug-naïve phenotypes. Eight of the 15 gene deletions (53%) significantly affected intake or preference for three alcohol, eight methamphetamine or three both (3). We observed multiple relations between predisposing behaviors and drug intake, revealing many distinct biobehavioral processes underlying these relationships. The set of mouse models identified in this study can be used to characterize these addiction-related processes further.

Comprehensive ECG reference intervals in C57BL/6N substrains provide a generalizable guide for cardiac electrophysiology studies in mice.

Reference ranges provide a powerful tool for diagnostic decision-making in clinical medicine and are enormously valuable for understanding normality in pre-clinical scientific research that uses in vivo models. As yet, there are no published reference ranges for electrocardiography (ECG) in the laboratory mouse. The first mouse-specific reference ranges for the assessment of electrical conduction are reported herein generated from an ECG dataset of unprecedented scale. International Mouse Phenotyping Consortium data from over 26,000 conscious or anesthetized C57BL/6N wildtype control mice were stratified by sex and age to develop robust ECG reference ranges. Interesting findings include that heart rate and key elements from the ECG waveform (RR-, PR-, ST-, QT-interval, QT corrected, and QRS complex) demonstrate minimal sexual dimorphism. As expected, anesthesia induces a decrease in heart rate and was shown for both inhalation (isoflurane) and injectable (tribromoethanol) anesthesia. In the absence of pharmacological, environmental, or genetic challenges, we did not observe major age-related ECG changes in C57BL/6N-inbred mice as the differences in the reference ranges of 12-week-old compared to 62-week-old mice were negligible. The generalizability of the C57BL/6N substrain reference ranges was demonstrated by comparison with ECG data from a wide range of non-IMPC studies. The close overlap in data from a wide range of mouse strains suggests that the C57BL/6N-based reference ranges can be used as a robust and comprehensive indicator of normality. We report a unique ECG reference resource of fundamental importance for any experimental study of cardiac function in mice.

The spontaneous mouse mutant low set ears (Lse) is caused by tandem duplication of Fgf3 and Fgf4.

The external ear develops from an organized convergence of ventrally migrating neural crest cells into the first and second branchial arches. Defects in external ear position are often symptomatic of complex syndromes such as Apert, Treacher-Collins, and Crouzon Syndrome. The low set ears (Lse) spontaneous mouse mutant is characterized by the dominant inheritance of a ventrally shifted external ear position and an abnormal external auditory meatus (EAM). We identified the causative mutation as a 148 Kb tandem duplication on Chromosome 7, which includes the entire coding sequences of Fgf3 and Fgf4. Duplications of FGF3 and FGF4 occur in 11q duplication syndrome in humans and are associated with craniofacial anomalies, among other features. Intercrosses of Lse-affected mice revealed perinatal lethality in homozygotes, and Lse/Lse embryos display additional phenotypes including polydactyly, abnormal eye morphology, and cleft secondary palate. The duplication results in increased Fgf3 and Fgf4 expression in the branchial arches and additional discrete domains in the developing embryo. This ectopic overexpression resulted in functional FGF signaling, demonstrated by increased Spry2 and Etv5 expression in overlapping domains of the developing arches. Finally, a genetic interaction between Fgf3/4 overexpression and Twist1, a regulator of skull suture development, resulted in perinatal lethality, cleft palate, and polydactyly in compound heterozygotes. These data indicate a role for Fgf3 and Fgf4 in external ear and palate development and provide a novel mouse model for further interrogation of the biological consequences of human FGF3/4 duplication.

Whole genome analysis for 163 gRNAs in Cas9-edited mice reveals minimal off-target activity.

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.

A single-cell transcriptional timelapse of mouse embryonic development, from gastrula to pup.

The house mouse, Mus musculus, is an exceptional model system, combining genetic tractability with close homology to human biology. Gestation in mouse development lasts just under three weeks, a period during which its genome orchestrates the astonishing transformation of a single cell zygote into a free-living pup composed of >500 million cells. Towards a global framework for exploring mammalian development, we applied single cell combinatorial indexing (sci-*) to profile the transcriptional states of 12.4 million nuclei from 83 precisely staged embryos spanning late gastrulation (embryonic day 8 or E8) to birth (postnatal day 0 or P0), with 2-hr temporal resolution during somitogenesis, 6-hr resolution through to birth, and 20-min resolution during the immediate postpartum period. From these data (E8 to P0), we annotate dozens of trajectories and hundreds of cell types and perform deeper analyses of the unfolding of the posterior embryo during somitogenesis as well as the ontogenesis of the kidney, mesenchyme, retina, and early neurons. Finally, we leverage the depth and temporal resolution of these whole embryo snapshots, together with other published data, to construct and curate a rooted tree of cell type relationships that spans mouse development from zygote to pup. Throughout this tree, we systematically nominate sets of transcription factors (TFs) and other genes as candidate drivers of the in vivo differentiation of hundreds of mammalian cell types. Remarkably, the most dramatic shifts in transcriptional state are observed in a restricted set of cell types in the hours immediately following birth, and presumably underlie the massive changes in physiology that must accompany the successful transition of a placental mammal to extrauterine life.

Identification of a heterogeneous and dynamic ciliome during embryonic development and cell differentiation.

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.

TMEM161B modulates radial glial scaffolding in neocortical development.

TMEM161B encodes an evolutionarily conserved widely expressed novel 8-pass transmembrane protein of unknown function in human. Here we identify TMEM161B homozygous hypomorphic missense variants in our recessive polymicrogyria (PMG) cohort. Patients carrying TMEM161B mutations exhibit striking neocortical PMG and intellectual disability. Tmem161b knockout mice fail to develop midline hemispheric cleavage, whereas knock-in of patient mutations and patient-derived brain organoids show defects in apical cell polarity and radial glial scaffolding. We found that TMEM161B modulates actin filopodia, functioning upstream of the Rho-GTPase CDC42. Our data link TMEM161B with human PMG, likely regulating radial glia apical polarity during neocortical development.