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.

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.

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.

Genome-wide screening reveals the genetic basis of mammalian embryonic eye development.

BACKGROUND: Microphthalmia, anophthalmia, and coloboma (MAC) spectrum disease encompasses a group of eye malformations which play a role in childhood visual impairment. Although the predominant cause of eye malformations is known to be heritable in nature, with 80% of cases displaying loss-of-function mutations in the ocular developmental genes OTX2 or SOX2, the genetic abnormalities underlying the remaining cases of MAC are incompletely understood. This study intended to identify the novel genes and pathways required for early eye development. Additionally, pathways involved in eye formation during embryogenesis are also incompletely understood. This study aims to identify the novel genes and pathways required for early eye development through systematic forward screening of the mammalian genome. RESULTS: Query of the International Mouse Phenotyping Consortium (IMPC) database (data release 17.0, August 01, 2022) identified 74 unique knockout lines (genes) with genetically associated eye defects in mouse embryos. The vast majority of eye abnormalities were small or absent eyes, findings most relevant to MAC spectrum disease in humans. A literature search showed that 27 of the 74 lines had previously published knockout mouse models, of which only 15 had ocular defects identified in the original publications. These 12 previously published gene knockouts with no reported ocular abnormalities and the 47 unpublished knockouts with ocular abnormalities identified by the IMPC represent 59 genes not previously associated with early eye development in mice. Of these 59, we identified 19 genes with a reported human eye phenotype. Overall, mining of the IMPC data yielded 40 previously unimplicated genes linked to mammalian eye development. Bioinformatic analysis showed that several of the IMPC genes colocalized to several protein anabolic and pluripotency pathways in early eye development. Of note, our analysis suggests that the serine-glycine pathway producing glycine, a mitochondrial one-carbon donator to folate one-carbon metabolism (FOCM), is essential for eye formation. CONCLUSIONS: Using genome-wide phenotype screening of single-gene knockout mouse lines, STRING analysis, and bioinformatic methods, this study identified genes heretofore unassociated with MAC phenotypes providing models to research novel molecular and cellular mechanisms involved in eye development. These findings have the potential to hasten the diagnosis and treatment of this congenital blinding disease.

Inflammation Impacts Androgen Receptor Signaling in Basal Prostate Stem Cells Through Interleukin 1 Receptor Antagonist.

The majority of patients with benign prostate hyperplasia (BPH) exhibit chronic prostate inflammation and the extent of inflammation correlates with the severity of symptoms. How inflammation contributes to prostate enlargement and/or BPH symptoms and the underlying mechanisms are not clearly understood. We established a unique mouse model Prostate Ovalbumin Expressing Transgenic 3 (POET3) that mimics chronic non-bacterial prostatitis in men to study the role of inflammation in prostate hyperplasia. After the injection of ovalbumin peptide-specific T cells, POET3 prostates exhibited an influx of inflammatory cells and an increase in pro-inflammatory cytokines that led to epithelial and stromal hyperplasia. We have previously demonstrated with the POET3 model that inflammation expands the basal prostate stem cell (bPSC) population and promotes bPSC differentiation in organoid cultures. In this study, we investigated the mechanisms underlying the impact of inflammation on bPSC. We found that AR activity was enhanced in inflamed bPSC and was essential for bPSC differentiation in organoid cultures. Most importantly, we identified, for the first time, interleukin 1 receptor antagonist (IL-1RA) as a key regulator of AR in basal stem cells. IL-1RA was one of the top genes upregulated by inflammation and inhibition of IL-1RA abrogated the enhanced AR nuclear accumulation and activity in organoids derived from inflamed bPSC. The mirroring effects of IL-1RA recombinant protein and IL-1α neutralizing antibody suggest that IL-1RA may function by antagonizing IL-1α inhibition of AR expression. Furthermore, we established a lineage tracing model to follow bPSC during inflammation and under castrate conditions. We found that inflammation induced bPSC proliferation and differentiation into luminal cells even under castrate conditions, indicating that AR activation driven by inflammation in bPSC is sufficient for their proliferation and differentiation under androgen-deprived conditions. However, proliferation of the differentiated bPSC in the luminal layer significantly diminished with castration, suggesting inflammation may not maintain AR activity in stromal cells, as stromal cells deprived of androgen after castration could no longer provide paracrine growth factors essential for luminal proliferation. Taken together, we have discovered novel mechanisms through which inflammation modulates AR signaling in bPSC and induces bPSC luminal differentiation that contributes to prostate hyperplasia.

Analysis of genome-wide knockout mouse database identifies candidate ciliopathy genes.

We searched a database of single-gene knockout (KO) mice produced by the International Mouse Phenotyping Consortium (IMPC) to identify candidate ciliopathy genes. We first screened for phenotypes in mouse lines with both ocular and renal or reproductive trait abnormalities. The STRING protein interaction tool was used to identify interactions between known cilia gene products and those encoded by the genes in individual knockout mouse strains in order to generate a list of "candidate ciliopathy genes." From this list, 32 genes encoded proteins predicted to interact with known ciliopathy proteins. Of these, 25 had no previously described roles in ciliary pathobiology. Histological and morphological evidence of phenotypes found in ciliopathies in knockout mouse lines are presented as examples (genes Abi2, Wdr62, Ap4e1, Dync1li1, and Prkab1). Phenotyping data and descriptions generated on IMPC mouse line are useful for mechanistic studies, target discovery, rare disease diagnosis, and preclinical therapeutic development trials. Here we demonstrate the effective use of the IMPC phenotype data to uncover genes with no previous role in ciliary biology, which may be clinically relevant for identification of novel disease genes implicated in ciliopathies.

Mendelian gene identification through mouse embryo viability screening.

BACKGROUND: The diagnostic rate of Mendelian disorders in sequencing studies continues to increase, along with the pace of novel disease gene discovery. However, variant interpretation in novel genes not currently associated with disease is particularly challenging and strategies combining gene functional evidence with approaches that evaluate the phenotypic similarities between patients and model organisms have proven successful. A full spectrum of intolerance to loss-of-function variation has been previously described, providing evidence that gene essentiality should not be considered as a simple and fixed binary property. METHODS: Here we further dissected this spectrum by assessing the embryonic stage at which homozygous loss-of-function results in lethality in mice from the International Mouse Phenotyping Consortium, classifying the set of lethal genes into one of three windows of lethality: early, mid, or late gestation lethal. We studied the correlation between these windows of lethality and various gene features including expression across development, paralogy and constraint metrics together with human disease phenotypes. We explored a gene similarity approach for novel gene discovery and investigated unsolved cases from the 100,000 Genomes Project. RESULTS: We found that genes in the early gestation lethal category have distinct characteristics and are enriched for genes linked with recessive forms of inherited metabolic disease. We identified several genes sharing multiple features with known biallelic forms of inborn errors of the metabolism and found signs of enrichment of biallelic predicted pathogenic variants among early gestation lethal genes in patients recruited under this disease category. We highlight two novel gene candidates with phenotypic overlap between the patients and the mouse knockouts. CONCLUSIONS: Information on the developmental period at which embryonic lethality occurs in the knockout mouse may be used for novel disease gene discovery that helps to prioritise variants in unsolved rare disease cases.

Identifying genetic determinants of inflammatory pain in mice using a large-scale gene-targeted screen.

ABSTRACT: Identifying the genetic determinants of pain is a scientific imperative given the magnitude of the global health burden that pain causes. Here, we report a genetic screen for nociception, performed under the auspices of the International Mouse Phenotyping Consortium. A biased set of 110 single-gene knockout mouse strains was screened for 1 or more nociception and hypersensitivity assays, including chemical nociception (formalin) and mechanical and thermal nociception (von Frey filaments and Hargreaves tests, respectively), with or without an inflammatory agent (complete Freund's adjuvant). We identified 13 single-gene knockout strains with altered nocifensive behavior in 1 or more assays. All these novel mouse models are openly available to the scientific community to study gene function. Two of the 13 genes (Gria1 and Htr3a) have been previously reported with nociception-related phenotypes in genetically engineered mouse strains and represent useful benchmarking standards. One of the 13 genes (Cnrip1) is known from human studies to play a role in pain modulation and the knockout mouse reported herein can be used to explore this function further. The remaining 10 genes (Abhd13, Alg6, BC048562, Cgnl1, Cp, Mmp16, Oxa1l, Tecpr2, Trim14, and Trim2) reveal novel pathways involved in nociception and may provide new knowledge to better understand genetic mechanisms of inflammatory pain and to serve as models for therapeutic target validation and drug development.

Mouse mutant phenotyping at scale reveals novel genes controlling bone mineral density.

The genetic landscape of diseases associated with changes in bone mineral density (BMD), such as osteoporosis, is only partially understood. Here, we explored data from 3,823 mutant mouse strains for BMD, a measure that is frequently altered in a range of bone pathologies, including osteoporosis. A total of 200 genes were found to significantly affect BMD. This pool of BMD genes comprised 141 genes with previously unknown functions in bone biology and was complementary to pools derived from recent human studies. Nineteen of the 141 genes also caused skeletal abnormalities. Examination of the BMD genes in osteoclasts and osteoblasts underscored BMD pathways, including vesicle transport, in these cells and together with in silico bone turnover studies resulted in the prioritization of candidate genes for further investigation. Overall, the results add novel pathophysiological and molecular insight into bone health and disease.

Machine learning-based automated phenotyping of inflammatory nocifensive behavior in mice.

The discovery and development of new and potentially nonaddictive pain therapeutics requires rapid, yet clinically relevant assays of nociception in preclinical models. A reliable and scalable automated scoring system for nocifensive behavior of mice in the formalin assay would dramatically lower the time and labor costs associated with experiments and reduce experimental variability. Here, we present a method that exploits machine learning techniques for video recordings that consists of three components: key point detection, per frame feature extraction using these key points, and classification of behavior using the GentleBoost algorithm. This approach to automation is flexible as different model classifiers or key points can be used with only small losses in accuracy. The adopted system identified the behavior of licking/biting of the hind paw with an accuracy that was comparable to a human observer (98% agreement) over 111 different short videos (total 284 min) at a resolution of 1 s. To test the system over longer experimental conditions, the responses of two inbred strains, C57BL/6NJ and C57BL/6J, were recorded over 90 min post formalin challenge. The automated system easily scored over 80 h of video and revealed strain differences in both response timing and amplitude. This machine learning scoring system provides the required accuracy, consistency, and ease of use that could make the formalin assay a feasible choice for large-scale genetic studies.

ADAD1 and ADAD2, testis-specific adenosine deaminase domain-containing proteins, are required for male fertility.

Adenosine-to-inosine RNA editing, a fundamental RNA modification, is regulated by adenosine deaminase (AD) domain containing proteins. Within the testis, RNA editing is catalyzed by ADARB1 and is regulated in a cell-type dependent manner. This study examined the role of two testis-specific AD domain proteins, ADAD1 and ADAD2, on testis RNA editing and male germ cell differentiation. ADAD1, previously shown to localize to round spermatids, and ADAD2 had distinct localization patterns with ADAD2 expressed predominantly in mid- to late-pachytene spermatocytes suggesting a role for both in meiotic and post-meiotic germ cell RNA editing. AD domain analysis showed the AD domain of both ADADs was likely catalytically inactive, similar to known negative regulators of RNA editing. To assess the impact of Adad mutation on male germ cell RNA editing, CRISPR-induced alleles of each were generated in mouse. Mutation of either Adad resulted in complete male sterility with Adad1 mutants displaying severe teratospermia and Adad2 mutant germ cells unable to progress beyond round spermatid. However, mutation of neither Adad1 nor Adad2 impacted RNA editing efficiency or site selection. Taken together, these results demonstrate ADAD1 and ADAD2 are essential regulators of male germ cell differentiation with molecular functions unrelated to A-to-I RNA editing.

GDNF promotes hair formation and cutaneous wound healing by targeting bulge stem cells.

Glial-cell-derived neurotrophic factor (GDNF) is a well-studied neuroregenerative factor; however, the degree to which it supports hair formation and skin wound repair is not known. By using a Gfra1 (GDNF family receptor alpha 1) knock-in reporter mouse line, GDNF signaling was found to occur within hair bulge stem cells (BSCs) during the initiation of the hair cycle and early stages of hair formation after depilation. Both recombinant and transgene overexpression of GDNF promoted BSC colony growth, hair formation, and skin repair after wounding through enhanced self-renewal of BSCs and commitment of BSC-derived progenitors into becoming epidermal cells at the injury site. Conditional ablation of Gfra1 among BSCs impaired the onset of the hair cycle, while conditional ablation of the GDNF family member signal transducer, Ret, within BSCs prevented the onset of the hair cycle and depilation-induced anagen development of hair follicles. Our findings reveal that GDNF promotes hair formation and wound repair and that bulge stem cells are critical mediators of both.

Human and mouse essentiality screens as a resource for disease gene discovery.

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.

An expanded mouse testis transcriptome and mass spectrometry defines novel proteins.

The testis transcriptome is exceptionally complex. Despite its complexity, previous testis transcriptome analyses relied on a reductive method for transcript identification, thus underestimating transcriptome complexity. We describe here a more complete testis transcriptome generated by combining Tuxedo, a reductive method, and spliced-RUM, a combinatorial transcript-building approach. Forty-two percent of the expanded testis transcriptome is composed of unannotated RNAs with novel isoforms of known genes and novel genes constituting 78 and 9.8% of the newly discovered transcripts, respectively. Across tissues, novel transcripts were predominantly expressed in the testis with the exception of novel isoforms which were also highly expressed in the adult ovary. Within the testis, novel isoform expression was distributed equally across all cell types while novel genes were predominantly expressed in meiotic and post-meiotic germ cells. The majority of novel isoforms retained their protein-coding potential while most novel genes had low protein-coding potential. However, a subset of novel genes had protein-coding potentials equivalent to known protein-coding genes. Shotgun mass spectrometry of round spermatid total protein identified unique peptides from four novel genes along with seven annotated non-coding RNAs. These analyses demonstrate the testis expresses a wide range of novel transcripts that give rise to novel proteins.

Identification of EOMES-expressing spermatogonial stem cells and their regulation by PLZF.

Long-term maintenance of spermatogenesis in mammals is supported by GDNF, an essential growth factor required for spermatogonial stem cell (SSC) self-renewal. Exploiting a transgenic GDNF overexpression model, which expands and normalizes the pool of undifferentiated spermatogonia between Plzf +/+ and Plzf lu/lu mice, we used RNAseq to identify a rare subpopulation of cells that express EOMES, a T-box transcription factor. Lineage tracing and busulfan challenge show that these are SSCs that contribute to steady state spermatogenesis as well as regeneration following chemical injury. EOMES+ SSCs have a lower proliferation index in wild-type than in Plzf lu/lu mice, suggesting that PLZF regulates their proliferative activity and that EOMES+ SSCs are lost through proliferative exhaustion in Plzf lu/lu mice. Single cell RNA sequencing of EOMES+ cells from Plzf +/+ and Plzf lu/lu mice support the conclusion that SSCs are hierarchical yet heterogeneous.

Robust mouse tracking in complex environments using neural networks.

The ability to track animals accurately is critical for behavioral experiments. For video-based assays, this is often accomplished by manipulating environmental conditions to increase contrast between the animal and the background in order to achieve proper foreground/background detection (segmentation). Modifying environmental conditions for experimental scalability opposes ethological relevance. The biobehavioral research community needs methods to monitor behaviors over long periods of time, under dynamic environmental conditions, and in animals that are genetically and behaviorally heterogeneous. To address this need, we applied a state-of-the-art neural network-based tracker for single mice. We compare three different neural network architectures across visually diverse mice and different environmental conditions. We find that an encoder-decoder segmentation neural network achieves high accuracy and speed with minimal training data. Furthermore, we provide a labeling interface, labeled training data, tuned hyperparameters, and a pretrained network for the behavior and neuroscience communities.

Erratum: Author Correction: Identification of genes required for eye development by high-throughput screening of mouse knockouts.

[This corrects the article DOI: 10.1038/s42003-018-0226-0.].

Large-scale discovery of mouse transgenic integration sites reveals frequent structural variation and insertional mutagenesis.

Transgenesis has been a mainstay of mouse genetics for over 30 yr, providing numerous models of human disease and critical genetic tools in widespread use today. Generated through the random integration of DNA fragments into the host genome, transgenesis can lead to insertional mutagenesis if a coding gene or an essential element is disrupted, and there is evidence that larger scale structural variation can accompany the integration. The insertion sites of only a tiny fraction of the thousands of transgenic lines in existence have been discovered and reported, due in part to limitations in the discovery tools. Targeted locus amplification (TLA) provides a robust and efficient means to identify both the insertion site and content of transgenes through deep sequencing of genomic loci linked to specific known transgene cassettes. Here, we report the first large-scale analysis of transgene insertion sites from 40 highly used transgenic mouse lines. We show that the transgenes disrupt the coding sequence of endogenous genes in half of the lines, frequently involving large deletions and/or structural variations at the insertion site. Furthermore, we identify a number of unexpected sequences in some of the transgenes, including undocumented cassettes and contaminating DNA fragments. We demonstrate that these transgene insertions can have phenotypic consequences, which could confound certain experiments, emphasizing the need for careful attention to control strategies. Together, these data show that transgenic alleles display a high rate of potentially confounding genetic events and highlight the need for careful characterization of each line to assure interpretable and reproducible experiments.