Somatic GJA4 gain-of-function mutation in orbital cavernous venous malformations.

Orbital cavernous venous malformation (OCVM) is a sporadic vascular anomaly of uncertain etiology characterized by abnormally dilated vascular channels. Here, we identify a somatic missense mutation, c.121G > T (p.Gly41Cys) in GJA4, which encodes a transmembrane protein that is a component of gap junctions and hemichannels in the vascular system, in OCVM tissues from 25/26 (96.2%) individuals with OCVM. GJA4 expression was detected in OCVM tissue including endothelial cells and the stroma, through immunohistochemistry. Within OCVM tissue, the mutation allele frequency was higher in endothelial cell-enriched fractions obtained using magnetic-activated cell sorting. Whole-cell voltage clamp analysis in Xenopus oocytes revealed that GJA4 c.121G > T (p.Gly41Cys) is a gain-of-function mutation that leads to the formation of a hyperactive hemichannel. Overexpression of the mutant protein in human umbilical vein endothelial cells led to a loss of cellular integrity, which was rescued by carbenoxolone, a non-specific gap junction/hemichannel inhibitor. Our data suggest that GJA4 c.121G > T (p.Gly41Cys) is a potential driver gene mutation for OCVM. We propose that hyperactive hemichannel plays a role in the development of this vascular phenotype.

Muscle Transcriptomics Shows Overexpression of Cadherin 1 in Inclusion Body Myositis.

OBJECTIVE: This study aimed to elucidate the molecular features of inclusion body myositis (IBM). METHODS: We performed RNA sequencing analysis of muscle biopsy samples from 67 participants, consisting of 58 myositis patients with the pathological finding of CD8-positive T cells invading non-necrotic muscle fibers expressing major histocompatibility complex class I (43 IBM, 6 polymyositis, and 9 unclassifiable myositis), and 9 controls. RESULTS: Cluster analysis, principal component analysis, and pathway analysis showed that differentially expressed genes and pathways identified in IBM and polymyositis were mostly comparable. However, pathways related to cell adhesion molecules were upregulated in IBM as compared with polymyositis and controls (p < 0.01). Notably, CDH1, which encodes the epidermal cell junction protein cadherin 1, was overexpressed in the muscles of IBM, which was validated by another RNA sequencing dataset from previous publications. Western blotting confirmed the presence of mature cadherin 1 protein in the muscles of IBM. Immunohistochemical staining confirmed the positivity for anti-cadherin 1 antibody in the muscles of IBM, whereas there was no muscle fiber positive for anti-cadherin 1 antibody in immune-mediated necrotizing myopathy, antisynthetase syndrome, and controls. The fibers stained with anti-cadherin 1 antibody did not have rimmed vacuoles or abnormal protein accumulation. Experimental skeletal muscle regeneration and differentiation systems showed that CDH1 is expressed during skeletal muscle regeneration and differentiation. INTERPRETATION: CDH1 was detected as a differentially expressed gene, and immunohistochemistry showed that cadherin 1 exists in the muscles of IBM, whereas it was rarely seen in those of other idiopathic inflammatory myopathies. Cadherin 1 upregulation in muscle could provide a valuable clue to the pathological mechanisms of IBM. ANN NEUROL 2022;91:317-328.

Hamster PIWI proteins bind to piRNAs with stage-specific size variations during oocyte maturation.

In animal gonads, transposable elements are actively repressed to preserve genome integrity through the PIWI-interacting RNA (piRNA) pathway. In mice, piRNAs are abundantly expressed in male germ cells, and form effector complexes with three distinct PIWIs. The depletion of individual Piwi genes causes male-specific sterility with no discernible phenotype in female mice. Unlike mice, most other mammals have four PIWI genes, some of which are expressed in the ovary. Here, purification of PIWI complexes from oocytes of the golden hamster revealed that the size of the PIWIL1-associated piRNAs changed during oocyte maturation. In contrast, PIWIL3, an ovary-specific PIWI in most mammals, associates with short piRNAs only in metaphase II oocytes, which coincides with intense phosphorylation of the protein. An improved high-quality genome assembly and annotation revealed that PIWIL1- and PIWIL3-associated piRNAs appear to share the 5'-ends of common piRNA precursors and are mostly derived from unannotated sequences with a diminished contribution from TE-derived sequences, most of which correspond to endogenous retroviruses. Our findings show the complex and dynamic nature of biogenesis of piRNAs in hamster oocytes, and together with the new genome sequence generated, serve as the foundation for developing useful models to study the piRNA pathway in mammalian oocytes.

Recompleting the Caenorhabditis elegans genome.

Caenorhabditis elegans was the first multicellular eukaryotic genome sequenced to apparent completion. Although this assembly employed a standard C. elegans strain (N2), it used sequence data from several laboratories, with DNA propagated in bacteria and yeast. Thus, the N2 assembly has many differences from any C. elegans available today. To provide a more accurate C. elegans genome, we performed long-read assembly of VC2010, a modern strain derived from N2. Our VC2010 assembly has 99.98% identity to N2 but with an additional 1.8 Mb including tandem repeat expansions and genome duplications. For 116 structural discrepancies between N2 and VC2010, 97 structures matching VC2010 (84%) were also found in two outgroup strains, implying deficiencies in N2. Over 98% of N2 genes encoded unchanged products in VC2010; moreover, we predicted ≥53 new genes in VC2010. The recompleted genome of C. elegans should be a valuable resource for genetics, genomics, and systems biology.

Targeted deep sequencing of DNA from multiple tissue types improves the diagnostic rate and reveals a highly diverse phenotype of mosaic neurofibromatosis type 2.

BACKGROUND: Although 60% of patients with de novo neurofibromatosis type 2 (NF2) are presumed to have mosaic NF2, the actual diagnostic rate of this condition remains low at around 20% because of the existing difficulties in detecting NF2 variants with low variant allele frequency (VAF). Here, we examined the correlation between the genotype and phenotype of mosaic NF2 after improving the diagnostic rate of mosaic NF2. METHODS: We performed targeted deep sequencing of 36 genes including NF2 using DNA samples from multiple tissues (blood, buccal mucosa, hair follicle and tumour) of 53 patients with de novo NF2 and elucidated their genotype-phenotype correlation. RESULTS: Twenty-four patients (45.2%) had the NF2 germline variant, and 20 patients with NF2 (37.7%) had mosaic NF2. The mosaic NF2 phenotype was significantly different from that in patients with NF2 germline variant in terms of distribution of NF2-related disease, tumour growth rate and hearing outcome. The behaviour of schwannoma correlated to the extent of VAF with NF2 variant in normal tissues unlike meningioma. CONCLUSION: We have improved the diagnostic rate of mosaic NF2 compared with that of previous studies by targeted deep sequencing of DNA from multiple tissues. Many atypical patients with NF2 diagnosed with 'unilateral vestibular schwannoma' or 'multiple meningiomas' presumably have mosaic NF2. Finally, we suggest that the highly diverse phenotype of NF2 could result not only from the type and location of NF2 variant but also the extent of VAF in the NF2 variant within normal tissue DNA.

Splice-site mutations in KIF5A in the Japanese case series of amyotrophic lateral sclerosis.

Our objective was to investigate the frequency of KIF5A variants in amyotrophic lateral sclerosis (ALS) and the clinical characteristics of familial ALS (FALS) associated with variants in KIF5A. Whole-exome sequence analysis was performed for a Japanese series of 43 families with FALS and 444 patients with sporadic ALS (SALS), in whom causative variants had not been identified. We compared the frequencies of rare variants (MAF < 0.01) in KIF5A, including missense and loss of function (LoF) variants, between ALS and control subjects (n = 1163). Clinical characteristics of patients with FALS carrying pathogenic variants in KIF5A were also described. LoF variants were identified only in the probands of two families with FALS, both of which were 3' splice-site variants leading to exon skipping and an altered C-terminal domain, located in the mutational hotspot causing FALS, and were considered to be pathogenic for FALS. Rare missense variants in KIF5A were identified in five patients with SALS (1.13%) and 11 control subjects (0.95%, carrier frequency), which were not significantly different. Consequently, the pathogenic LoF variants in KIF5A accounted for 2.1% of all FALS families in this study. These patients suffered from ALS characteristically associated with the predominant involvement of upper motor neuron. In conclusion, we identified two pathogenic splice-site variants in KIF5A in the probands in two Japanese families with FALS, which altered the C-terminal region of KIF5A. Our findings broaden the phenotype spectrum of ALS associated with variants in KIF5A in the Japanese series.

Loss-of-function variants in NEK1 are associated with an increased risk of sporadic ALS in the Japanese population.

Loss-of-function (LoF) variants in NEK1 have recently been reported to be associated with amyotrophic lateral sclerosis (ALS). In this study, we investigated the association of NEK1 LoF variants with an increased risk of sporadic ALS (SALS) and the clinical characteristics of patients with SALS carrying LoF variants in a Japanese case series. Whole-exome sequencing analysis was performed for a series of 446 SALS patients in whom pathogenic variants in familial ALS-causative genes have not been identified and 1163 healthy control subjects in our Japanese series. We evaluated LoF variants, defined as nonsense, splice-site disrupting single-nucleotide variants (SNVs), or short insertion/deletion (indel) variants predicted to cause frameshifts in NEK1. We identified seven NEK1 LoF variants in patients with SALS (1.57%), whereas only one was identified in control subjects (0.086%) (P = 0.00073, Fisher's exact test). This finding is consistent with those in recent reports from other regions in the world. In conclusion, we demonstrated that NEK1 LoF variants are also associated with an increased risk of SALS in the Japanese population.

Clinical and molecular genetic characterization of two female patients harboring the Xq27.3q28 deletion with different ratios of X chromosome inactivation.

A heterozygous deletion at Xq27.3q28 including FMR1, AFF2, and IDS causing intellectual disability and characteristic facial features is very rare in females, with only 10 patients having been reported. Here, we examined two female patients with different clinical features harboring the Xq27.3q28 deletion and determined the chromosomal breakpoints. Moreover, we assessed the X chromosome inactivation (XCI) in peripheral blood from both patients. Both patients had an almost overlapping deletion at Xq27.3q28, however, the more severe patient (Patient 1) showed skewed XCI of the normal X chromosome (79:21) whereas the milder patient (Patient 2) showed random XCI. Therefore, deletion at Xq27.3q28 critically affected brain development, and the ratio of XCI of the normal X chromosome greatly affected the clinical characteristics of patients with deletion at Xq27.3q28. As the chromosomal breakpoints were determined, we analyzed a change in chromatin domains termed topologically associated domains (TADs) using published Hi-C data on the Xq27.3q28 region, and found that only patient 1 had a possibility of a drastic change in TADs. The altered chromatin topologies on the Xq27.3q28 region might affect the clinical features of patient 1 by changing the expression of genes just outside the deletion and/or the XCI establishment during embryogenesis resulting in skewed XCI.

Integrative analysis of genomic alterations in triple-negative breast cancer in association with homologous recombination deficiency.

Triple-negative breast cancer (TNBC) cells do not express estrogen receptors, progesterone receptors, or human epidermal growth factor receptor 2. Currently, apart from poly ADP-ribose polymerase inhibitors, there are few effective therapeutic options for this type of cancer. Here, we present comprehensive characterization of the genetic alterations in TNBC performed by high coverage whole genome sequencing together with transcriptome and whole exome sequencing. Silencing of the BRCA1 gene impaired the homologous recombination pathway in a subset of TNBCs, which exhibited similar phenotypes to tumors with BRCA1 mutations; they harbored many structural variations (SVs) with relative enrichment for tandem duplication. Clonal analysis suggested that TP53 mutations and methylation of CpG dinucleotides in the BRCA1 promoter were early events of carcinogenesis. SVs were associated with driver oncogenic events such as amplification of MYC, NOTCH2, or NOTCH3 and affected tumor suppressor genes including RB1, PTEN, and KMT2C. Furthermore, we identified putative TGFA enhancer regions. Recurrent SVs that affected the TGFA enhancer region led to enhanced expression of the TGFA oncogene that encodes one of the high affinity ligands for epidermal growth factor receptor. We also identified a variety of oncogenes that could transform 3T3 mouse fibroblasts, suggesting that individual TNBC tumors may undergo a unique driver event that can be targetable. Thus, we revealed several features of TNBC with clinically important implications.

Ataxic phenotype with altered CaV3.1 channel property in a mouse model for spinocerebellar ataxia 42.

Spinocerebellar ataxia 42 (SCA42) is a neurodegenerative disorder recently shown to be caused by c.5144G > A (p.Arg1715His) mutation in CACNA1G, which encodes the T-type voltage-gated calcium channel CaV3.1. Here, we describe a large Japanese family with SCA42. Postmortem pathological examination revealed severe cerebellar degeneration with prominent Purkinje cell loss without ubiquitin accumulation in an SCA42 patient. To determine whether this mutation causes ataxic symptoms and neurodegeneration, we generated knock-in mice harboring c.5168G > A (p.Arg1723His) mutation in Cacna1g, corresponding to the mutation identified in the SCA42 family. Both heterozygous and homozygous mutants developed an ataxic phenotype from the age of 11-20 weeks and showed Purkinje cell loss at 50 weeks old. Degenerative change of Purkinje cells and atrophic thinning of the molecular layer were conspicuous in homozygous knock-in mice. Electrophysiological analysis of Purkinje cells using acute cerebellar slices from young mice showed that the point mutation altered the voltage dependence of CaV3.1 channel activation and reduced the rebound action potentials after hyperpolarization, although it did not significantly affect the basic properties of synaptic transmission onto Purkinje cells. Finally, we revealed that the resonance of membrane potential of neurons in the inferior olivary nucleus was decreased in knock-in mice, which indicates that p.Arg1723His CaV3.1 mutation affects climbing fiber signaling to Purkinje cells. Altogether, our study shows not only that a point mutation in CACNA1G causes an ataxic phenotype and Purkinje cell degeneration in a mouse model, but also that the electrophysiological abnormalities at an early stage of SCA42 precede Purkinje cell loss.

Burden of rare variants in causative genes for amyotrophic lateral sclerosis (ALS) accelerates age at onset of ALS.

OBJECTIVES: To evaluate the burden of rare variants in the causative genes for amyotrophic lateral sclerosis (ALS) on the age at onset of ALS in a Japanese case series. METHODS: We conducted whole-exome sequencing analysis of 89 families with familial ALS (FALS) and 410 patients with sporadic ALS (SALS) to identify known pathogenic mutations or rare functionally predicted deleterious variants in the causative genes for ALS. Rare variants (minor allele frequency <1%) with scaled Combined Annotation-Dependent Depletion score >20 were defined as rare functionally predicted deleterious variants. The patients with ALS were classified on the basis of the number of pathogenic and/or rare functionally predicted deleterious variants, and the age at onset was compared among the classified groups. RESULTS: Whole-exome sequencing analysis revealed known pathogenic mutations or rare functionally predicted deleterious variants in causative genes for ALS in 56 families with FALS (62.9%) and 87 patients with SALS (21.2%). Such variants in multiple genes were identified in seven probands with FALS and eight patients with SALS. The ages at onset in the patients with ALS with multiple variants were significantly earlier than those in other patients with ALS. Even when the patients with known pathogenic mutations were excluded, a significantly earlier onset of the disease was still observed in patients with multiple rare functionally predicted deleterious variants. CONCLUSIONS: A substantial number of patients carried rare variants in multiple genes, and the burden of rare variants in the known causative genes for ALS affects the age at onset in the Japanese ALS series.

Mutations in COA7 cause spinocerebellar ataxia with axonal neuropathy.

Several genes related to mitochondrial functions have been identified as causative genes of neuropathy or ataxia. Cytochrome c oxidase assembly factor 7 (COA7) may have a role in assembling mitochondrial respiratory chain complexes that function in oxidative phosphorylation. Here we identified four unrelated patients with recessive mutations in COA7 among a Japanese case series of 1396 patients with Charcot-Marie-Tooth disease (CMT) or other inherited peripheral neuropathies, including complex forms of CMT. We also found that all four patients had characteristic neurological features of peripheral neuropathy and ataxia with cerebellar atrophy, and some patients showed leukoencephalopathy or spinal cord atrophy on MRI scans. Validated mutations were located at highly conserved residues among different species and segregated with the disease in each family. Nerve conduction studies showed axonal sensorimotor neuropathy. Sural nerve biopsies showed chronic axonal degeneration with a marked loss of large and medium myelinated fibres. An immunohistochemical assay with an anti-COA7 antibody in the sural nerve from the control patient showed the positive expression of COA7 in the cytoplasm of Schwann cells. We also observed mildly elevated serum creatine kinase levels in all patients and the presence of a few ragged-red fibres and some cytochrome c oxidase-negative fibres in a muscle biopsy obtained from one patient, which was suggestive of subclinical mitochondrial myopathy. Mitochondrial respiratory chain enzyme assay in skin fibroblasts from the three patients showed a definitive decrease in complex I or complex IV. Immunocytochemical analysis of subcellular localization in HeLa cells indicated that mutant COA7 proteins as well as wild-type COA7 were localized in mitochondria, which suggests that mutant COA7 does not affect the mitochondrial recruitment and may affect the stability or localization of COA7 interaction partners in the mitochondria. In addition, Drosophila COA7 (dCOA7) knockdown models showed rough eye phenotype, reduced lifespan, impaired locomotive ability and shortened synaptic branches of motor neurons. Our results suggest that loss-of-function COA7 mutation is responsible for the phenotype of the presented patients, and this new entity of disease would be referred to as spinocerebellar ataxia with axonal neuropathy type 3.

Novel De Novo KCND3 Mutation in a Japanese Patient with Intellectual Disability, Cerebellar Ataxia, Myoclonus, and Dystonia.

Spinocerebellar ataxia 19/22 (SCA19/22) is a rare type of autosomal dominant SCA that was previously described in 11 families. We report the case of a 30-year-old Japanese man presenting with intellectual disability, early onset cerebellar ataxia, myoclonus, and dystonia without a family history. MRI showed cerebellar atrophy, and electroencephalograms showed paroxysmal sharp waves during hyperventilation and photic stimulation. Trio whole-exome sequencing analysis of DNA samples from the patient and his parents revealed a de novo novel missense mutation (c.1150G>A, p.G384S) in KCND3, the causative gene of SCA19/22, substituting for evolutionally conserved glycine. The mutation was predicted to be functionally deleterious by bioinformatic analysis. Although pure cerebellar ataxia is the most common clinical feature in SCA19/22 families, extracerebellar symptoms including intellectual disability and myoclonus are reported in a limited number of families, suggesting a genotype-phenotype correlation for particular mutations. Although autosomal recessive diseases are more common in patients with early onset sporadic cerebellar ataxia, the present study emphasizes that such a possibility of de novo mutation should be considered.

Analysis of a novel gene, Sdgc, reveals sex chromosome-dependent differences of medaka germ cells prior to gonad formation.

In vertebrates that have been examined to date, the sexual identity of germ cells is determined by the sex of gonadal somatic cells. In the teleost fish medaka, a sex-determination gene on the Y chromosome, DMY/dmrt1bY, is expressed in gonadal somatic cells and regulates the sexual identity of germ cells. Here, we report a novel mechanism by which sex chromosomes cell-autonomously confer sexually different characters upon germ cells prior to gonad formation in a genetically sex-determined species. We have identified a novel gene, Sdgc (sex chromosome-dependent differential expression in germ cells), whose transcripts are highly enriched in early XY germ cells. Chimeric analysis revealed that sexually different expression of Sdgc is controlled in a germ cell-autonomous manner by the number of Y chromosomes. Unexpectedly, DMY/dmrt1bY was expressed in germ cells prior to gonad formation, but knockdown and overexpression of DMY/dmrt1bY did not affect Sdgc expression. We also found that XX and XY germ cells isolated before the onset of DMY/dmrt1bY expression in gonadal somatic cells behaved differently in vitro and were affected by Sdgc. Sdgc maps close to the sex-determination locus, and recombination around the two loci appears to be repressed. Our results provide important insights into the acquisition and plasticity of sexual differences at the cellular level even prior to the developmental stage of sex determination.

AgIn: measuring the landscape of CpG methylation of individual repetitive elements.

MOTIVATION: Determining the methylation state of regions with high copy numbers is challenging for second-generation sequencing, because the read length is insufficient to map reads uniquely, especially when repetitive regions are long and nearly identical to each other. Single-molecule real-time (SMRT) sequencing is a promising method for observing such regions, because it is not vulnerable to GC bias, it produces long read lengths, and its kinetic information is sensitive to DNA modifications. RESULTS: We propose a novel linear-time algorithm that combines the kinetic information for neighboring CpG sites and increases the confidence in identifying the methylation states of those sites. Using a practical read coverage of ∼30-fold from an inbred strain medaka (Oryzias latipes), we observed that both the sensitivity and precision of our method on individual CpG sites were ∼93.7%. We also observed a high correlation coefficient (R = 0.884) between our method and bisulfite sequencing, and for 92.0% of CpG sites, methylation levels ranging over [0,1] were in concordance within an acceptable difference 0.25. Using this method, we characterized the landscape of the methylation status of repetitive elements, such as LINEs, in the human genome, thereby revealing the strong correlation between CpG density and hypomethylation and detecting hypomethylation hot spots of LTRs and LINEs. We uncovered the methylation states for nearly identical active transposons, two novel LINE insertions of identity ∼99% and length 6050 base pairs (bp) in the human genome, and 16 Tol2 elements of identity >99.8% and length 4682 bp in the medaka genome. AVAILABILITY AND IMPLEMENTATION: AgIn (Aggregate on Intervals) is available at: https://github.com/hacone/AgIn CONTACT: ysuzuki@cb.k.u-tokyo.ac.jp or moris@cb.k.u-tokyo.ac.jp SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

TBCD may be a causal gene in progressive neurodegenerative encephalopathy with atypical infantile spinal muscular atrophy.

Spinal muscular atrophy (SMA) is an autosomal recessive neurodegenerative disorder caused by survival motor neuron gene mutations. Variant forms of SMA accompanied by additional clinical presentations have been classified as atypical SMA and are thought to be caused by variants in as yet unidentified causative genes. Here, we presented the clinical findings of two siblings with an SMA variant followed by progressive cerebral atrophy, and the results of whole-exome sequencing analyses of the family quartet that was performed to identify potential causative variants. We identified two candidate homozygous missense variants, R942Q in the tubulin-folding cofactor D (TBCD) gene and H250Q in the bromo-adjacent homology domain and coiled-coil containing 1 (BAHCC1) gene, located on chromosome 17q25.3 with an interval of 1.4 Mbp. The in silico analysis of both variants suggested that TBCD rather than BAHCC1 was likely the pathogenic gene (TBCD sensitivity, 0.68; specificity, 0.97; BAHCC1 sensitivity, 1.00; specificity, 0.00). Thus, our results show that TBCD is a likely novel candidate gene for atypical SMA with progressive cerebral atrophy. TBCD is predicted to have important functions on tubulin integrity in motor neurons as well as in the central nervous system.

Mutations in MME cause an autosomal-recessive Charcot-Marie-Tooth disease type 2.

OBJECTIVE: The objective of this study was to identify new causes of Charcot-Marie-Tooth (CMT) disease in patients with autosomal-recessive (AR) CMT. METHODS: To efficiently identify novel causative genes for AR-CMT, we analyzed 303 unrelated Japanese patients with CMT using whole-exome sequencing and extracted recessive variants/genes shared among multiple patients. We performed mutation screening of the newly identified membrane metalloendopeptidase (MME) gene in 354 additional patients with CMT. We clinically, genetically, pathologically, and radiologically examined 10 patients with the MME mutation. RESULTS: We identified recessive mutations in MME in 10 patients. The MME gene encodes neprilysin (NEP), which is well known to be one of the most prominent beta-amyloid (Aß)-degrading enzymes. All patients had a similar phenotype consistent with late-onset axonal neuropathy. They showed muscle weakness, atrophy, and sensory disturbance in the lower extremities. All the MME mutations could be loss-of-function mutations, and we confirmed a lack/decrease of NEP protein expression in a peripheral nerve. No patients showed symptoms of dementia, and 1 patient showed no excess Aß in Pittsburgh compound-B positron emission tomography imaging. INTERPRETATION: Our results indicate that loss-of-function MME mutations are the most frequent cause of adult-onset AR-CMT2 in Japan, and we propose that this new disease should be termed AR-CMT2T. A loss-of-function MME mutation did not cause early-onset Alzheimer's disease. Identifying the MME mutation responsible for AR-CMT could improve the rate of molecular diagnosis and the understanding of the molecular mechanisms of CMT.

Identification of candidate genes involved in the etiology of sporadic Tourette syndrome by exome sequencing.

Tourette Syndrome (TS) is a neurodevelopmental disorder characterized by chronic motor and vocal tics. Although there is a large genetic contribution, the genetic architecture of TS remains unclear. Exome sequencing has successfully revealed the contribution of de novo mutations in sporadic cases with neuropsychiatric disorders such as autism and schizophrenia. Here, using exome sequencing, we investigated de novo mutations in individuals with sporadic TS to identify novel risk loci and elucidate the genetic background of TS. Exome analysis was conducted for sporadic TS cases: nine trio families and one quartet family with concordant twins were investigated. Missense mutations were evaluated using functional prediction algorithms, and their population frequencies were calculated based on three public databases. Gene expression patterns in the brain were analyzed using the BrainSpan Developmental Transcriptome. Thirty de novo mutations, including four synonymous and four missense mutations, were identified. Among the missense mutations, one in the rapamycin-insensitive companion of mammalian target of rapamycin (RICTOR)-coding gene (rs140964083: G > A, found in one proband) was predicted to be hazardous. In the three public databases analyzed, variants in the same SNP locus were absent, and variants in the same gene were either absent or present at an extremely low frequency (3/5,008), indicating the rarity of hazardous RICTOR mutations in the general population. The de novo variant of RICTOR may be implicated in the development of sporadic TS, and RICTOR is a novel candidate factor for TS etiology.

ERBB4 mutations that disrupt the neuregulin-ErbB4 pathway cause amyotrophic lateral sclerosis type 19.

Amyotrophic lateral sclerosis (ALS) is a devastating neurological disorder characterized by the degeneration of motor neurons and typically results in death within 3-5 years from onset. Familial ALS (FALS) comprises 5%-10% of ALS cases, and the identification of genes associated with FALS is indispensable to elucidating the molecular pathogenesis. We identified a Japanese family affected by late-onset, autosomal-dominant ALS in which mutations in genes known to be associated with FALS were excluded. A whole- genome sequencing and parametric linkage analysis under the assumption of an autosomal-dominant mode of inheritance with incomplete penetrance revealed the mutation c.2780G>A (p. Arg927Gln) in ERBB4. An extensive mutational analysis revealed the same mutation in a Canadian individual with familial ALS and a de novo mutation, c.3823C>T (p. Arg1275Trp), in a Japanese simplex case. These amino acid substitutions involve amino acids highly conserved among species, are predicted as probably damaging, and are located within a tyrosine kinase domain (p. Arg927Gln) or a C-terminal domain (p. Arg1275Trp), both of which mediate essential functions of ErbB4 as a receptor tyrosine kinase. Functional analysis revealed that these mutations led to a reduced autophosphorylation of ErbB4 upon neuregulin-1 (NRG-1) stimulation. Clinical presentations of the individuals with mutations were characterized by the involvement of both upper and lower motor neurons, a lack of obvious cognitive dysfunction, and relatively slow progression. This study indicates that disruption of the neuregulin-ErbB4 pathway is involved in the pathogenesis of ALS and potentially paves the way for the development of innovative therapeutic strategies such using NRGs or their agonists to upregulate ErbB4 functions.

Human genetic variation database, a reference database of genetic variations in the Japanese population.

Whole-genome and -exome resequencing using next-generation sequencers is a powerful approach for identifying genomic variations that are associated with diseases. However, systematic strategies for prioritizing causative variants from many candidates to explain the disease phenotype are still far from being established, because the population-specific frequency spectrum of genetic variation has not been characterized. Here, we have collected exomic genetic variation from 1208 Japanese individuals through a collaborative effort, and aggregated the data into a prevailing catalog. In total, we identified 156 622 previously unreported variants. The allele frequencies for the majority (88.8%) were lower than 0.5% in allele frequency and predicted to be functionally deleterious. In addition, we have constructed a Japanese-specific major allele reference genome by which the number of unique mapping of the short reads in our data has increased 0.045% on average. Our results illustrate the importance of constructing an ethnicity-specific reference genome for identifying rare variants. All the collected data were centralized to a newly developed database to serve as useful resources for exploring pathogenic variations. Public access to the database is available at http://www.genome.med.kyoto-u.ac.jp/SnpDB/.