Early embryonic mutations reveal dynamics of somatic and germ cell lineages in mice.

De novo mutations accumulate with zygotic cell divisions. However, the occurrence of these mutations and the way they are inherited by somatic cells and germ cells remain unclear. Here, we present a novel method to reconstruct cell lineages. We identified mosaic mutations in mice using deep whole-genome sequencing and reconstructed embryonic cell lineages based on the variant allele frequencies of the mutations. The reconstructed trees were confirmed using nuclear transfer experiments and the genotyping of approximately 50 offspring of each tree. The most detailed tree had 32 terminal nodes and showed cell divisions from the fertilized egg to germ cell- and somatic cell-specific lineages, indicating at least five independent cell lineages that would be selected as founders of the primordial germ cells. The contributions of each lineage to germ cells and offspring varied widely. At the emergence of the germ cell-specific lineages, 10-15 embryonic mutations had accumulated, suggesting that the pregastrulation mutation rate is 1.0 mutation per mitosis. Subsequent mutation rates were 0.7 for germ cells and 13.2 for tail fibroblasts. Our results show a new framework to assess embryonic lineages; further, we suggest an evolutionary strategy for preserving heterogeneity owing to postzygotic mutations in offspring.

Mice with endogenous TDP-43 mutations exhibit gain of splicing function and characteristics of amyotrophic lateral sclerosis.

TDP-43 (encoded by the gene TARDBP) is an RNA binding protein central to the pathogenesis of amyotrophic lateral sclerosis (ALS). However, how TARDBP mutations trigger pathogenesis remains unknown. Here, we use novel mouse mutants carrying point mutations in endogenous Tardbp to dissect TDP-43 function at physiological levels both in vitro and in vivo Interestingly, we find that mutations within the C-terminal domain of TDP-43 lead to a gain of splicing function. Using two different strains, we are able to separate TDP-43 loss- and gain-of-function effects. TDP-43 gain-of-function effects in these mice reveal a novel category of splicing events controlled by TDP-43, referred to as "skiptic" exons, in which skipping of constitutive exons causes changes in gene expression. In vivo, this gain-of-function mutation in endogenous Tardbp causes an adult-onset neuromuscular phenotype accompanied by motor neuron loss and neurodegenerative changes. Furthermore, we have validated the splicing gain-of-function and skiptic exons in ALS patient-derived cells. Our findings provide a novel pathogenic mechanism and highlight how TDP-43 gain of function and loss of function affect RNA processing differently, suggesting they may act at different disease stages.

Establishment of mouse line showing inducible priapism-like phenotypes.

Purpose: Penile research is expected to reveal new targets for treatment and prevention of the complex mechanisms of its disorder including erectile dysfunction (ED). Thus, analyses of the molecular processes of penile ED and continuous erection as priapism are essential issues of reproductive medicine. Methods: By performing mouse N-ethyl-N-nitrosourea mutagenesis and exome sequencing, we established a novel mouse line displaying protruded genitalia phenotype (PGP; priapism-like phenotype) and identified a novel Pitpna gene mutation for PGP. Extensive histological analyses on the Pitpna mutant and intracavernous pressure measurement (ICP) and liquid chromatography-electrospray ionization tandem mass spectrometry (LC-ESI/MS)/MS analyses were performed. Results: We evaluated the role of phospholipids during erection for the first time and showed the mutants of inducible phenotypes of priapism. Moreover, quantitative analysis using LC-ESI/MS/MS revealed that the level of phosphatidylinositol (PI) was significantly lower in the mutant penile samples. These results imply that PI may contribute to penile erection by PITPα. Conclusions: Our findings suggest that the current mutant is a mouse model for priapism and abnormalities in PI signaling pathways through PITPα may lead to priapism providing an attractive novel therapeutic target in its treatment.

Crim1C140S mutant mice reveal the importance of cysteine 140 in the internal region 1 of CRIM1 for its physiological functions.

Cysteine-rich transmembrane bone morphogenetic protein regulator 1 (CRIM1) is a type I transmembrane protein involved in the organogenesis of many tissues via its interactions with growth factors including BMP, TGF-ß, and VEGF. In this study, we used whole-exome sequencing and linkage analysis to identify a novel Crim1 mutant allele generated by ENU mutagenesis in mice. This allele is a missense mutation that causes a cysteine-to-serine substitution at position 140, and is referred to as Crim1C140S. In addition to the previously reported phenotypes in Crim1 mutants, Crim1C140S homozygous mice exhibited several novel phenotypes, including dwarfism, enlarged seminal vesicles, and rectal prolapse. In vitro analyses showed that Crim1C140S mutation affected the formation of CRIM1 complexes and decreased the amount of the overexpressed CRIM1 proteins in the cell culture supernatants. Cys140 is located in the internal region 1 (IR1) of the N-terminal extracellular region of CRIM1 and resides outside any identified functional domains. Inference of the domain architecture suggested that the Crim1C140S mutation disturbs an intramolecular disulfide bond in IR1, leading to the protein instability and the functional defects of CRIM1. Crim1C140S highlights the functional importance of the IR1, and Crim1C140S mice should serve as a valuable model for investigating the functions of CRIM1 that are unidentified as yet.

The linear ubiquitin-specific deubiquitinase gumby regulates angiogenesis.

A complex interaction of signalling events, including the Wnt pathway, regulates sprouting of blood vessels from pre-existing vasculature during angiogenesis. Here we show that two distinct mutations in the (uro)chordate-specific gumby (also called Fam105b) gene cause an embryonic angiogenic phenotype in gumby mice. Gumby interacts with disheveled 2 (DVL2), is expressed in canonical Wnt-responsive endothelial cells and encodes an ovarian tumour domain class of deubiquitinase that specifically cleaves linear ubiquitin linkages. A crystal structure of gumby in complex with linear diubiquitin reveals how the identified mutations adversely affect substrate binding and catalytic function in line with the severity of their angiogenic phenotypes. Gumby interacts with HOIP (also called RNF31), a key component of the linear ubiquitin assembly complex, and decreases linear ubiquitination and activation of NF-κB-dependent transcription. This work provides support for the biological importance of linear (de)ubiquitination in angiogenesis, craniofacial and neural development and in modulating Wnt signalling.

Mouse models for ROS1-fusion-positive lung cancers and their application to the analysis of multikinase inhibitor efficiency.

ROS1-fusion genes, resulting from chromosomal rearrangement, have been reported in 1-2% of human non-small cell lung cancer cases. More than 10 distinct ROS1-fusion genes, including break-point variants, have been identified to date. In this study, to investigate the in vivo oncogenic activities of one of the most frequently detected fusions, CD74-ROS1, as well as another SDC4-ROS1 fusion that has also been reported in several studies, we generated transgenic (TG) mouse strains that express either of the two ROS1-fusion genes specifically in lung alveolar type II cells. Mice in all TG lines developed tumorigenic nodules in the lung, and a few strains of both TG mouse lines demonstrated early-onset nodule development (multiple tumor lesions present in the lung at 2-4 weeks after birth); therefore, these two strains were selected for further investigation. Tumors developed progressively in the untreated TG mice of both lines, whereas those receiving oral administration of an ALK/MET/ROS1 inhibitor, crizotinib, and an ALK/ROS1 inhibitor, ASP3026, showed marked reduction in the tumor burden. Collectively, these data suggest that each of these two ROS1-fusion genes acts as a driver for the pathogenesis of lung adenocarcinoma in vivo The TG mice developed in this study are expected to serve as valuable tools for exploring novel therapeutic agents against ROS1-fusion-positive lung cancer.

Effects of background mutations and single nucleotide polymorphisms (SNPs) on the Disc1 L100P behavioral phenotype associated with schizophrenia in mice.

BACKGROUND: Disrupted-in-schizophrenia 1 (DISC1) is a promising candidate susceptibility gene for psychiatric disorders, including schizophrenia, bipolar disorder and major depression. Several previous studies reported that mice with N-ethyl-N-nitrosourea (ENU)-induced L100P mutation in Disc1 showed some schizophrenia-related behavioral phenotypes. This line originally carried several thousands of ENU-induced point mutations in the C57BL/6 J strain and single nucleotide polymorphisms (SNPs) from the DBA/2 J inbred strain. METHODS: To investigate the effect of Disc1 L100P, background mutations and SNPs on phenotypic characterization, we performed behavioral analyses to better understand phenotypes of Disc1 L100P mice and comprehensive genetic analyses using whole-exome resequencing and SNP panels to map ENU-induced mutations and strain-specific SNPs, respectively. RESULTS: We found no differences in spontaneous or methamphetamine-induced locomotor activity, sociability or social novelty preference among Disc1 L100P/L100P, L100P/+ mutants and wild-type littermates. Whole-exome resequencing of the original G1 mouse identified 117 ENU-induced variants, including Disc1 L100P per se. Two females and three males from the congenic L100P strain after backcrossing to C57BL/6 J were deposited to RIKEN BioResource Center in 2008. We genotyped them with DBA/2 J × C57BL/6 J SNPs and found a number of the checked SNPs still remained. CONCLUSION: These results suggest that causal attribution of the discrepancy in behavioral phenotypes to the Disc1 L100P mutant mouse line existing among different research groups needs to be cautiously investigated in further study by taking into account the effect(s) of other ENU-induced mutations and/or SNPs from DBA/2 J.

Tumor cell enrichment by tissue suspension improves sensitivity to copy number variation in diffuse gastric cancer with low tumor content.

The detection of copy number variations (CNVs) and somatic mutations in cancer is important for the selection of specific drugs for patients with cancer. In cancers with sporadic tumor cells, low tumor content prevents the accurate detection of somatic alterations using targeted sequencing. To efficiently identify CNVs, we performed tumor cell enrichment using tissue suspensions of formalin-fixed paraffin-embedded (FFPE) tissue sections with low tumor cell content. Tumor-enriched and residual fractions were separated from FFPE tissue suspensions of intestinal and diffuse-type gastric cancers containing sporadic tumor cells, and targeted sequencing was performed on 225 cancer-related genes. Sequencing of a targeted panel of cancer-related genes using tumor-enriched fractions increased the number of detectable CNVs and the copy number of amplified genes. Furthermore, CNV analysis using the normal cell-enriched residual fraction as a reference for CNV scoring allowed targeted sequencing to detect CNV characteristics of diffuse-type gastric cancer with low tumor content. Our approach improves the CNV detection rate in targeted sequencing with tumor enrichment and the accuracy of CNV detection in archival samples without paired blood.

Tumor cell enrichment by tissue suspension enables detection of mutations with low variant allele frequency and estimation of germline mutations.

Targeted sequencing offers an opportunity to select specific drugs for cancer patients based on alterations in their genome. However, accurate sequencing cannot be performed in cancers harboring diffuse tumor cells because of low tumor content. We performed tumor cell enrichment using tissue suspension of formalin-fixed, paraffin-embedded (FFPE) tissue sections with low tumor cell content. The enriched fractions were used to efficiently identify mutations by sequencing a target panel of cancer-related genes. Tumor-enriched and residual fractions were isolated from FFPE tissue sections of intestinal and diffuse gastric cancers harboring diffuse tumor cells and DNA of suitable quality was isolated for next-generation sequencing. Sequencing of a target panel of cancer-related genes using the tumor-enriched fraction increased the number of detectable mutations and variant allele frequency. Furthermore, mutation analysis of DNA isolated from tumor-enriched and residual fractions allowed us to estimate germline mutations without a blood reference. This approach of tumor cell enrichment will not only enhance the success rate of target panel sequencing, but can also improve the accuracy of detection of somatic mutations in archived specimens.

Serine racemase is associated with schizophrenia susceptibility in humans and in a mouse model.

Abnormal N-methyl-d-aspartate receptor (NMDAR) function has been implicated in the pathophysiology of schizophrenia. d-serine is an important NMDAR modulator, and to elucidate the role of the d-serine synthesis enzyme serine racemase (Srr) in schizophrenia, we identified and characterized mice with an ENU-induced mutation that results in a complete loss of Srr activity and dramatically reduced d-serine levels. Mutant mice displayed behaviors relevant to schizophrenia, including impairments in prepulse inhibition, sociability and spatial discrimination. Behavioral deficits were exacerbated by an NMDAR antagonist and ameliorated by d-serine or the atypical antipsychotic clozapine. Expression profiling revealed that the Srr mutation influenced several genes that have been linked to schizophrenia and cognitive ability. Transcript levels altered by the Srr mutation were also normalized by d-serine or clozapine treatment. Furthermore, analysis of SRR genetic variants in humans identified a robust association with schizophrenia. This study demonstrates that aberrant Srr function and diminished d-serine may contribute to schizophrenia pathogenesis.

A Comprehensive Mouse Transcriptomic BodyMap across 17 Tissues by RNA-seq.

The mouse has been widely used as a model organism for studying human diseases and for evaluating drug safety and efficacy. Many diseases and drug effects exhibit tissue specificity that may be reflected by tissue-specific gene-expression profiles. Here we construct a comprehensive mouse transcriptomic BodyMap across 17 tissues of six-weeks old C57BL/6JJcl mice using RNA-seq. We find different expression patterns between protein-coding and non-coding genes. Liver expressed the least complex transcriptomes, that is, the smallest number of genes detected in liver across all 17 tissues, whereas testis and ovary harbor more complex transcriptomes than other tissues. We report a comprehensive list of tissue-specific genes across 17 tissues, along with a list of 4,781 housekeeping genes in mouse. In addition, we propose a list of 27 consistently and highly expressed genes that can be used as reference controls in expression-profiling analysis. Our study provides a unique resource of mouse gene-expression profiles, which is helpful for further biomedical research.

An ENU-induced p.C225S missense mutation in the mouse Tgfb1 gene does not cause Camurati-Engelmann disease-like skeletal phenotypes.

Camurati-Engelmann disease (CED) is a rare sclerosing bone disorder in humans with autosomal dominant inheritance. Mutations in the gene (TGFB1) that encodes transforming growth factor-ß1 (TGF-ß1) are causative for CED. TGF-ß1 signaling is enhanced by the CED-causing mutations. In this study, we performed Tgfb1 mutation screening in an ENU-mutagenized mouse genomic DNA library. We identified a missense mutation in which cysteine was substituted by serine at position 225 (p.C225S), that corresponded to the CED-causing mutation (p.C225R). TGF-ß1 mutant protein carrying p.C225S was secreted normally into the extracellular space. Reporter gene assays showed that the p.C225S mutants enhanced TGF-ß signaling at the same level as p.C225R mutants. We generated p.C225S homozygous mice and confirmed that the mature TGF-ß1 levels in the culture supernatants of the calvarial cells from the homozygotes were significantly higher than those from wild-type mice. Although the skull and femur are sclerotic in CED, these phenotypes were not observed in p.C225S homozygous mice. These results suggest that human and mouse bone tissue react differently to TGF-ß1. These findings are useful to pharmacological studies using mouse models in developing drugs that will target TGF-ß signaling.

Identification of genes that express in response to light exposure and express rhythmically in a circadian manner in the mouse suprachiasmatic nucleus.

Most biological phenomena, including behavior and metabolic pathways, are governed by an internal clock system that is circadian (i.e., with a period of approximately 24 h) and is reset by light exposure from outside. In order to understand the molecular basis of the resetting mechanism of the clock, we attempted to isolate light-inducible transcripts in the suprachiasmatic nucleus, where the master clock resides, using a new gene expression profiling procedure. We identified 87 such transcripts, successfully cloned 60 of them and confirmed their light inducibility. Six of the 60 were already known to be light inducible and 17 are protein-coding transcripts registered in the public database that were not known to be light inducible. Induction is subjective night specific in most of the transcripts. Interestingly, 6 of the transcripts exhibit rhythmic expression in a circadian manner in the suprachiasmatic nucleus.

A sensitive transcriptome analysis method that can detect unknown transcripts.

We have developed an AFLP-based gene expression profiling method called 'high coverage expression profiling' (HiCEP) analysis. By making improvements to the selective PCR technique we have reduced the rate of false positive peaks to approximately 4% and consequently the number of peaks, including overlapping peaks, has been markedly decreased. As a result we can determine the relationship between peaks and original transcripts unequivocally. This will make it practical to prepare a database of all peaks, allowing gene assignment without having to isolate individual peaks. This precise selection also enables us to easily clone peaks of interest and predict the corresponding gene for each peak in some species. The procedure is highly reproducible and sensitive enough to detect even a 1.2-fold difference in gene expression. Most importantly, the low false positive rate enables us to analyze gene expression with wide coverage by means of four instead of six nucleotide recognition site restriction enzymes for fingerprinting mRNAs. Therefore, the method detects 70-80% of all transcripts, including non-coding transcripts, unknown and known genes. Moreover, the method requires no sequence information and so is applicable even to eukaryotes for which there is no genome information available.

Illegitimate translation causes unexpected gene expression from on-target out-of-frame alleles created by CRISPR-Cas9.

CRISPR-Cas9 is efficient enough to knock out both alleles directly by introducing out-of-frame mutations. We succeeded in making biallelic on-target frameshift mutations of the endogenous Gli3 gene; however, the GLI3 protein was expressed in all six of the established cell lines carrying homozygous out-of-frame mutations. We developed a dual-tagged expression vector and proved that illegitimate translation (ITL) was the cause of the unexpected Gli3 expression. Thus, gene expression must be examined even if designed on-target out-of-frame sequences are introduced by genome editing. In addition, it is highly recommended to pre-examine the occurrence of ITL in vitro prior to the design and construction of any genome-editing vectors. In vitro assay systems such as the dual-tagged ITL assay system developed in this study should aid the identification and elucidation of ITL-based human diseases and gene expression.

A normalization strategy applied to HiCEP (an AFLP-based expression profiling) analysis: toward the strict alignment of valid fragments across electrophoretic patterns.

BACKGROUND: Gene expression analysis based on comparison of electrophoretic patterns is strongly dependent on the accuracy of DNA fragment sizing. The current normalization strategy based on molecular weight markers has limited accuracy because marker peaks are often masked by intense peaks nearby. Cumulative errors in fragment lengths cause problems in the alignment of same-length fragments across different electropherograms, especially for small fragments (< 100 bp). For accurate comparison of electrophoretic patterns, further inspection and normalization of electrophoretic data after fragment sizing by conventional strategies is needed. RESULTS: Here we describe a method for the normalization of a set of time-course electrophoretic data to be compared. The method uses Gaussian curves fitted to the complex peak mixtures in each electropherogram. It searches for target ranges for which patterns are dissimilar to the other patterns (called "dissimilar ranges") and for references (a kind of mean or typical pattern) in the set of resultant approximate patterns. It then constructs the optimal normalized pattern whose correlation coefficient against the reference in the range achieves the highest value among various combinations of candidates. We applied the procedure to time-course electrophoretic data produced by HiCEP, an AFLP-based expression profiling method which can detect a slight expression change in DNA fragments. We obtained dissimilar ranges whose electrophoretic patterns were obviously different from the reference and as expected, most of the fragments in the detected ranges were short (< 100 bp). The normalized electrophoretic patterns also agreed well with reference patterns. CONCLUSION: The normalization strategy presented here demonstrates the importance of pre-processing before electrophoretic signal comparison, and we anticipate its usefulness especially for temporal expression analysis by the electrophoretic method.

T396I mutation of mouse Sufu reduces the stability and activity of Gli3 repressor.

Hedgehog signaling is primarily transduced by two transcription factors: Gli2, which mainly acts as a full-length activator, and Gli3, which tends to be proteolytically processed from a full-length form (Gli3FL) to an N-terminal repressor (Gli3REP). Recent studies using a Sufu knockout mouse have indicated that Sufu is involved in regulating Gli2 and Gli3 activator and repressor activity at multiple steps of the signaling cascade; however, the mechanism of specific Gli2 and Gli3 regulation remains to be elucidated. In this study, we established an allelic series of ENU-induced mouse strains. Analysis of one of the missense alleles, SufuT396I, showed that Thr396 residue of Sufu played a key role in regulation of Gli3 activity. SufuT396I/T396I embryos exhibited severe polydactyly, which is indicative of compromised Gli3 activity. Concomitantly, significant quantitative reductions of unprocessed Gli3 (Gli3FL) and processed Gli3 (Gli3REP) were observed in vivo as well as in vitro. Genetic experiments showed that patterning defects in the limb buds of SufuT396I/T396I were rescued by a constitutive Gli3REP allele (Gli3∆699), strongly suggesting that SufuT396I reduced the truncated Gli3 repressor. In contrast, SufuT396I qualitatively exhibited no mutational effects on Gli2 regulation. Taken together, the results of this study show that the Thr396 residue of Sufu is specifically required for regulation of Gli3 but not Gli2. This implies a novel Sufu-mediated mechanism in which Gli2 activator and Gli3 repressor are differentially regulated.

Mouse dexamethasone-induced RAS protein 1 gene is expressed in a circadian rhythmic manner in the suprachiasmatic nucleus.

We identified the Dexamethasone-induced RAS protein 1 (Dexras1) gene as a cycling gene in the suprachiasmatic nucleus (SCN). Investigation of the whole brain using in situ hybridization demonstrated the localization of the expression of the gene in the SCN, thalamus, piriform cortex and hippocampus. However, rhythmic expression of the gene was observed only in the SCN. The rhythmic change in gene expression during 1 day was approximately five-fold, and the maximum expression was observed during subjective night. Real-time PCR using the SCN, paraventricular nucleus and cortex confirmed these results. Next, we analyzed the expression of the Dexras1 gene in the SCN of cryptochrome (Cry) 1 and 2 double knockout mice. We found that the rhythmic expression disappeared. The results indicate that Dexras1 rhythmicity and levels are dependent upon CRYs. This is the first time that the G protein, which may be involved in the input pathway, has been isolated as a cycling gene in the SCN.

Sequence analysis of 193.4 and 83.9 kbp of mouse and chicken genomic DNAs containing the entire Prkdc (DNA-PKcs) gene.

The catalytic subunit of DNA-dependent protein kinase plays critical roles in nonhomologous end joining in repair of DNA double-strand breaks and V(D)J recombination. In addition to the SCID phenotype, it has been suggested that the molecule contributes to the polymorphic variations in radiosensitivity and susceptibility to cancer in mouse strains. Here we show the nucleotide sequence of approximately 193-kbp and 84-kbp genomic regions encoding the entire Prkdc gene (also known as DNA-PKcs) in the mouse and chicken, respectively. A large retroposon was found in intron 51 in the mouse but not in the human or chicken. Comparative analyses of the genome strongly suggested that the region contains only two genes for Prkdc and Mcm4; however, several conserved sequences and cis elements were also predicted.

ß-CateninC429S mice exhibit sterility consequent to spatiotemporally sustained Wnt signalling in the internal genitalia.

Wnt/ß-catenin signalling regulates numerous developmental and homeostatic processes. Ctnnb1 (also known as ß-catenin) is the only protein that transmits signals from various Wnt ligands to downstream genes. In this study, we report that our newly established mouse strain, which harbours a Cys429 to Ser missense mutation in the ß-catenin gene, exhibited specific organ defects in contrast to mice with broadly functioning Wnt/ß-catenin signalling. Both homozygous mutant males and females produced normal gametes but were infertile because of abnormal seminal vesicle and vaginal morphogenesis. An ins-TOPGAL transgenic reporter spatiotemporally sustained Wnt/ß-catenin signalling during the corresponding organogenesis. Therefore, ß-catenin(C429S) should provide new insights into ß-catenin as a universal component of Wnt/ß-catenin signal transduction.