Exploring the genetic diversity of the Japanese population: Insights from a large-scale whole genome sequencing analysis.

The Japanese archipelago is a terminal location for human migration, and the contemporary Japanese people represent a unique population whose genomic diversity has been shaped by multiple migrations from Eurasia. We analyzed the genomic characteristics that define the genetic makeup of the modern Japanese population from a population genetics perspective from the genomic data of 9,287 samples obtained by high-coverage whole-genome sequencing (WGS) by the National Center Biobank Network. The dataset comprised populations from the Ryukyu Islands and other parts of the Japanese archipelago (Hondo). The Hondo population underwent two episodes of population decline during the Jomon period, corresponding to the Late Neolithic, and the Edo period, corresponding to the Early Modern era, while the Ryukyu population experienced a population decline during the shell midden period of the Late Neolithic in this region. Haplotype analysis suggested increased allele frequencies for genes related to alcohol and fatty acid metabolism, which were reported as loci that had experienced positive natural selection. Two genes related to alcohol metabolism were found to be 12,500 years out of phase with the time when they began to increase in the allele frequency; this finding indicates that the genomic diversity of Japanese people has been shaped by events closely related to agriculture and food production.

Expansion of GGC Repeat in GIPC1 Is Associated with Oculopharyngodistal Myopathy.

Oculopharyngodistal myopathy (OPDM) is an adult-onset inherited neuromuscular disorder characterized by progressive ptosis, external ophthalmoplegia, and weakness of the masseter, facial, pharyngeal, and distal limb muscles. The myopathological features are presence of rimmed vacuoles (RVs) in the muscle fibers and myopathic changes of differing severity. Inheritance is variable, with either putative autosomal-dominant or autosomal-recessive pattern. Here, using a comprehensive strategy combining whole-genome sequencing (WGS), long-read whole-genome sequencing (LRS), linkage analysis, repeat-primed polymerase chain reaction (RP-PCR), and fluorescence amplicon length analysis polymerase chain reaction (AL-PCR), we identified an abnormal GGC repeat expansion in the 5' UTR of GIPC1 in one out of four families and three sporadic case subjects from a Chinese OPDM cohort. Expanded GGC repeats were further confirmed as the cause of OPDM in an additional 2 out of 4 families and 6 out of 13 sporadic Chinese individuals with OPDM, as well as 7 out of 194 unrelated Japanese individuals with OPDM. Methylation, qRT-PCR, and western blot analysis indicated that GIPC1 mRNA levels were increased while protein levels were unaltered in OPDM-affected individuals. RNA sequencing indicated p53 signaling, vascular smooth muscle contraction, ubiquitin-mediated proteolysis, and ribosome pathways were involved in the pathogenic mechanisms of OPDM-affected individuals with GGC repeat expansion in GIPC1. This study provides further evidence that OPDM is associated with GGC repeat expansions in distinct genes and highly suggests that expanded GGC repeat units are essential in the pathogenesis of OPDM, regardless of the genes in which the expanded repeats are located.

Pathogenic Variants in the Myosin Chaperone UNC-45B Cause Progressive Myopathy with Eccentric Cores.

The myosin-directed chaperone UNC-45B is essential for sarcomeric organization and muscle function from Caenorhabditis elegans to humans. The pathological impact of UNC-45B in muscle disease remained elusive. We report ten individuals with bi-allelic variants in UNC45B who exhibit childhood-onset progressive muscle weakness. We identified a common UNC45B variant that acts as a complex hypomorph splice variant. Purified UNC-45B mutants showed changes in folding and solubility. In situ localization studies further demonstrated reduced expression of mutant UNC-45B in muscle combined with abnormal localization away from the A-band towards the Z-disk of the sarcomere. The physiological relevance of these observations was investigated in C. elegans by transgenic expression of conserved UNC-45 missense variants, which showed impaired myosin binding for one and defective muscle function for three. Together, our results demonstrate that UNC-45B impairment manifests as a chaperonopathy with progressive muscle pathology, which discovers the previously unknown conserved role of UNC-45B in myofibrillar organization.

Distinctive chaperonopathy in skeletal muscle associated with the dominant variant in DNAJB4.

DnaJ homolog, subfamily B, member 4, a member of the heat shock protein 40 chaperones encoded by DNAJB4, is highly expressed in myofibers. We identified a heterozygous c.270 T > A (p.F90L) variant in DNAJB4 in a family with a dominantly inherited distal myopathy, in which affected members have specific features on muscle pathology represented by the presence of cytoplasmic inclusions and the accumulation of desmin, p62, HSP70, and DNAJB4 predominantly in type 1 fibers. Both Dnajb4F90L knockin and knockout mice developed muscle weakness and recapitulated the patient muscle pathology in the soleus muscle, where DNAJB4 has the highest expression. These data indicate that the identified variant is causative, resulting in defective chaperone function and selective muscle degeneration in specific muscle fibers. This study demonstrates the importance of DNAJB4 in skeletal muscle proteostasis by identifying the associated chaperonopathy.

Intranuclear inclusions in skin biopsies are not limited to neuronal intranuclear inclusion disease but can also be seen in oculopharyngodistal myopathy.

AIMS: Oculopharyngodistal myopathy (OPDM) is caused by the expansion of CGG repeats in NOTCH2NLC (OPDM_NOTCH2NLC) GIPC1 (OPDM_GIPC1), or LRP12 (OPDM_LRP12). Neuronal intranuclear inclusion disease (NIID) is clinically distinct from OPDM but is also caused by the expansion of CGG repeats in NOTCH2NLC, which may be an indicator of intranuclear inclusion in skin biopsy. We investigated the presence of intranuclear inclusions in skin biopsies from patients with OPDM and muscle diseases with a similar pathology to evaluate whether they will have similar diagnostic findings on skin biopsy. METHODS: We analysed the frequency of p62-positive intranuclear inclusions in sweat gland cells, adipocytes and fibroblasts in skin biopsy samples from patients with OPDM (OPDM_NOTCH2NLC [n = 2], OPDM_GIPC1 [n = 6] and OPDM_LRP12 [n = 3]), NIID (n = 1), OPMD (n = 1), IBM (n = 4) and GNE myopathy (n = 2). RESULTS: The p62-postive intranuclear inclusions were observed in all three cell types in both patients with OPDM_NOTCH2NLC and a patient with NIID, in at least one cell type in all six patients with OPDM_GIPC1, and all in three cell types in one of the three patients with OPDM_LRP12. These findings were not observed in patients with OPMD, IBM or GNE myopathy. CONCLUSION: Intranuclear inclusions in skin biopsy samples are not specific to NIID and are found in all three types of genetically confirmed OPDM, suggesting that the underlying mechanism of OPDM may be similar to NIID, regardless of causative genes.

Whole genome sequencing of 45 Japanese patients with intellectual disability.

Intellectual disability (ID) is characterized by significant limitations in both intellectual functioning and adaptive behaviors, originating before the age of 18 years. However, the genetic etiologies of ID are still incompletely elucidated due to the wide range of clinical and genetic heterogeneity. Whole genome sequencing (WGS) has been applied as a single-step clinical diagnostic tool for ID because it detects genetic variations with a wide range of resolution from single nucleotide variants (SNVs) to structural variants (SVs). To explore the causative genes for ID, we employed WGS in 45 patients from 44 unrelated Japanese families and performed a stepwise screening approach focusing on the coding variants in the genes. Here, we report 12 pathogenic and likely pathogenic variants: seven heterozygous variants of ADNP, SATB2, ANKRD11, PTEN, TCF4, SPAST, and KCNA2, three hemizygous variants of SMS, SLC6A8, and IQSEC2, and one homozygous variant in AGTPBP1. Of these, four were considered novel. Furthermore, a novel 76 kb deletion containing exons 1 and 2 in DYRK1A was identified. We confirmed the clinical and genetic heterogeneity and high frequency of de novo causative variants (8/12, 66.7%). This is the first report of WGS analysis in Japanese patients with ID. Our results would provide insight into the correlation between novel variants and expanded phenotypes of the disease.

COX6A2 variants cause a muscle-specific cytochrome c oxidase deficiency.

OBJECTIVE: Cytochrome c oxidase (COX) deficiency is a major mitochondrial respiratory chain defect that has vast genetic and phenotypic heterogeneity. This study aims to identify novel causative genes of COX deficiency with only striated muscle-specific symptoms. METHODS: Whole exome sequencing was performed in 2 unrelated individuals who were diagnosed with congenital myopathy and presented COX deficiency in muscle pathology. We assessed the COX6A2 variants using measurements of enzymatic activities and assembly of mitochondrial respiratory chain complexes in the samples from the patients and knockout mice. RESULTS: Both patients presented muscle weakness and hypotonia in 4 limbs along with facial muscle weakness. One patient had cardiomyopathy. Neither patient exhibited involvement from other organs. Whole exome sequencing identified biallelic missense variants in COX6A2, which is expressed only in the skeletal muscle and heart. The variants detected were homozygous c.117C > A (p.Ser39Arg) and compound heterozygous c.117C > A (p.Ser39Arg) and c.127T > C (p.Cys43Arg). We found specific reductions in complex IV activities in the skeletal muscle of both individuals. Assembly of complex IV and its supercomplex formation were impaired in the muscle. INTERPRETATION: This study indicates that biallelic variants in COX6A2 cause a striated muscle-specific form of COX deficiency. ANN NEUROL 2019;86:193-202.

RNA sequencing solved the most common but unrecognized NEB pathogenic variant in Japanese nemaline myopathy.

PURPOSE: The diagnostic rate for Mendelian diseases by exome sequencing (ES) is typically 20-40%. The low rate is partly because ES misses deep-intronic or synonymous variants leading to aberrant splicing. In this study, we aimed to apply RNA sequencing (RNA-seq) to efficiently detect the aberrant splicings and their related variants. METHODS: Aberrant splicing in biopsied muscles from six nemaline myopathy (NM) cases unresolved by ES were analyzed with RNA-seq. Variants related to detected aberrant splicing events were analyzed with Sanger sequencing. Detected variants were screened in NM patients unresolved by ES. RESULTS: We identified a novel deep-intronic NEB pathogenic variant, c.1569+339A>G in one case, and another novel synonymous NEB pathogenic variant, c.24684G>C (p.Ser8228Ser) in three cases. The c.24684G>C variant was observed to be the most frequent among all NEB pathogenic variants in normal Japanese populations with a frequency of 1 in 178 (20 alleles in 3552 individuals), but was previously unrecognized. Expanded screening of the variant identified it in a further four previously unsolved nemaline myopathy cases. CONCLUSION: These results indicated that RNA-seq may be able to solve a large proportion of previously undiagnosed muscle diseases.

A Functional SNP in BNC2 Is Associated with Adolescent Idiopathic Scoliosis.

Adolescent idiopathic scoliosis (AIS) is the most common spinal deformity. We previously conducted a genome-wide association study (GWAS) and detected two loci associated with AIS. To identify additional loci, we extended our GWAS by increasing the number of cohorts (2,109 affected subjects and 11,140 control subjects in total) and conducting a whole-genome imputation. Through the extended GWAS and replication studies using independent Japanese and Chinese populations, we identified a susceptibility locus on chromosome 9p22.2 (p = 2.46 × 10(-13); odds ratio = 1.21). The most significantly associated SNPs were in intron 3 of BNC2, which encodes a zinc finger transcription factor, basonuclin-2. Expression quantitative trait loci data suggested that the associated SNPs have the potential to regulate the BNC2 transcriptional activity and that the susceptibility alleles increase BNC2 expression. We identified a functional SNP, rs10738445 in BNC2, whose susceptibility allele showed both higher binding to a transcription factor, YY1 (yin and yang 1), and higher BNC2 enhancer activity than the non-susceptibility allele. BNC2 overexpression produced body curvature in developing zebrafish in a gene-dosage-dependent manner. Our results suggest that increased BNC2 expression is implicated in the etiology of AIS.

Compound Heterozygosity for Null Mutations and a Common Hypomorphic Risk Haplotype in TBX6 Causes Congenital Scoliosis.

Congenital scoliosis (CS) occurs as a result of vertebral malformations and has an incidence of 0.5-1/1,000 births. Recently, TBX6 on chromosome 16p11.2 was reported as a disease gene for CS; about 10% of Chinese CS patients were compound heterozygotes for rare null mutations and a common haplotype defined by three SNPs in TBX6. All patients had hemivertebrae. We recruited 94 Japanese CS patients, investigated the TBX6 locus for both mutations and the risk haplotype, examined transcriptional activities of mutant TBX6 in vitro, and evaluated clinical and radiographic features. We identified TBX6 null mutations in nine patients, including a missense mutation that had a loss of function in vitro. All had the risk haplotype in the opposite allele. One of the mutations showed dominant negative effect. Although all Chinese patients had one or more hemivertebrae, two Japanese patients did not have hemivertebra. The compound heterozygosity of null mutations and the common risk haplotype in TBX6 also causes CS in Japanese patients with similar incidence. Hemivertebra was not a specific type of spinal malformation in TBX6-associated CS (TACS). A heterozygous TBX6 loss-of-function mutation has been reported in a family with autosomal-dominant spondylocostal dysostosis, but it may represent a spectrum of the same disease with TACS.

Axial spondylometaphyseal dysplasia is also caused by NEK1 mutations.

Axial spondylometaphyseal dysplasia (axial SMD) is a unique form of SMD characterized by dysplasia of axial skeleton and retinal dystrophy. Recently, C21orf2 has been identified as the first disease gene for axial SMD; however, the presence of genetic heterogeneity is known. In this study, we identified NEK1 as the second disease gene for axial SMD. By whole-exome sequencing in a patient with axial SMD, we identified compound heterozygous mutations of NEK1, c.3107C>G (p.S1036*) and c.3830A>C (p.D1277A), which co-segregated in the family. NEK1 mutations have previously been found in three types of short rib thoracic dystrophy, which have no retinal dystrophy. The skeletal phenotype of our patient was milder than those of previously reported cases with NEK1 mutations and those with axial SMD harboring C21orf2 mutations. Phenotypes associated with NEK1 mutations are variable and the phenotype-genotype corelation in skeletal ciliopathies is challenging.

Identification of a novel LRRK1 mutation in a family with osteosclerotic metaphyseal dysplasia.

Osteosclerotic metaphyseal dysplasia (OSMD) is a rare skeletal dysplasia characterized by osteosclerotic metaphyses with osteopenic diaphyses of the long tubular bones. Our previous study identified a homozygous elongation mutation in leucine-rich repeat kinase 1 gene (LRRK1) in a patient with OSMD and showed that Lrrk1 knockout mice exhibited phenotypic similarity with OSMD. Here we report a second LRRK1 mutation in Indian sibs with OSMD. They had homozygous mutation (c.5971_5972insG) that produces an elongated mutant protein (p.A1991Gfs*31) similar to the first case. The sibs had normal stature, normal intelligence and recurrent fractures. The common radiographic feature was asymmetric and variable sclerosis of vertebral end plates, pelvic margin and metaphyses of tubular bones. One of the sibs had facial dysmorphisms, dentine abnormalities and acro-osteolysis. A comparison between the three OSMD cases with LRRK1 mutations with different ages suggested that the sclerotic lesions resolved with age. Our findings further support that LRRK1 would cause a subset of OSMD cases.

Novel and recurrent XYLT1 mutations in two Turkish families with Desbuquois dysplasia, type 2.

Desbuquois dysplasia (DBQD) is an autosomal recessive skeletal disorder characterized by growth retardation, joint laxity, short extremities, and progressive scoliosis. DBQD is classified into two types based on the presence (DBQD1) or absence (DBQD2) of characteristic hand abnormalities. CANT1 mutations have been reported in both DBQD1 and DBQD2. Recently, mutations in the gene encoding xylosyltransferase 1 (XYLT1) were identified in several families with DBQD2. In this study, we performed whole-exome sequencing in two Turkish families with DBQD2. We found a novel and a recurrent XYLT1 mutation in each family. The patients were homozygous for the mutations. Our results further support that XYLT1 is responsible for a major subset of DBQD2.

Identification of biallelic LRRK1 mutations in osteosclerotic metaphyseal dysplasia and evidence for locus heterogeneity.

BACKGROUND: Osteosclerotic metaphyseal dysplasia (OSMD) is a unique form of osteopetrosis characterised by severe osteosclerosis localised to the bone ends. The mode of inheritance is autosomal recessive. Its genetic basis is not known. OBJECTIVE: To identify the disease gene for OSMD. METHODS AND RESULTS: By whole exome sequencing in a boy with OSMD, we identified a homozygous 7 bp deletion (c.5938_5944delGAGTGGT) in the LRRK1 gene. His skeletal phenotype recapitulated that seen in the Lrrk1-deficient mouse. The shared skeletal hallmarks included severe sclerosis in the undermodelled metaphyses and epiphyseal margins of the tubular bones, costal ends, vertebral endplates and margins of the flat bones. The deletion is predicted to result in an elongated LRRK1 protein (p.E1980Afs*66) that lacks a part of its WD40 domains. In vitro functional studies using osteoclasts from Lrrk1-deficient mice showed that the deletion was a loss of function mutation. Genetic analysis of LRRK1 in two unrelated patients with OSMD suggested that OSMD is a genetically heterogeneous condition. CONCLUSIONS: This is the first study to identify the causative gene of OSMD. Our study provides evidence that LRRK1 plays a critical role in the regulation of bone mass in humans.

Mutations in B3GALT6, which encodes a glycosaminoglycan linker region enzyme, cause a spectrum of skeletal and connective tissue disorders.

Proteoglycans (PGs) are a major component of the extracellular matrix in many tissues and function as structural and regulatory molecules. PGs are composed of core proteins and glycosaminoglycan (GAG) side chains. The biosynthesis of GAGs starts with the linker region that consists of four sugar residues and is followed by repeating disaccharide units. By exome sequencing, we found that B3GALT6 encoding an enzyme involved in the biosynthesis of the GAG linker region is responsible for a severe skeletal dysplasia, spondyloepimetaphyseal dysplasia with joint laxity type 1 (SEMD-JL1). B3GALT6 loss-of-function mutations were found in individuals with SEMD-JL1 from seven families. In a subsequent candidate gene study based on the phenotypic similarity, we found that B3GALT6 is also responsible for a connective tissue disease, Ehlers-Danlos syndrome (progeroid form). Recessive loss-of-function mutations in B3GALT6 result in a spectrum of disorders affecting a broad range of skeletal and connective tissues characterized by lax skin, muscle hypotonia, joint dislocation, and spinal deformity. The pleiotropic phenotypes of the disorders indicate that B3GALT6 plays a critical role in a wide range of biological processes in various tissues, including skin, bone, cartilage, tendon, and ligament.

A rapid functional decline type of amyotrophic lateral sclerosis is linked to low expression of TTN.

OBJECTIVE: To classify the patterns of functional decline in patients with sporadic amyotrophic lateral sclerosis (ALS) and explore the genetic backgrounds that modified these patterns. METHODS: We included 465 patients with sporadic ALS in the analysis and clustered the longitudinal functional scores in the registered patients, using a mixture approach of a non-linear mixed-effects model. We conducted a genome-wide analysis of 572 983 single nucleotide polymorphisms (SNPs). We then assessed the association between the clusters of longitudinal functional scores and SNPs. RESULTS: We identified the following four clusters of longitudinal functional decline in the cases: a rapid decline cluster, an intermediate decline cluster, a sigmoidal decline cluster and a moderate decline cluster. We identified seven SNPs associated with the rapid decline cluster, using a recessive model (p=3.47-8.34×10(-8)). The OR for the probabilities of the rapid decline cluster ranged from 5.5 to 5.84. Homozygosity for the minor alleles in the seven SNPs, which constituted a linkage disequilibrium (LD) block, was associated with decreased expression of TTN (encoding Titin, a large sarcomere protein) in the expression quantitative trait loci database of a large-scale Japanese genetic variation database (p=8.6×10(-10)-1.1×10(-7)). TTN expression in immortalised lymphocyte lines was decreased in patients who were homozygous for the minor alleles compared with those who were homozygous for the major alleles (n=19 in each group, p=0.002). CONCLUSIONS: We detected an LD block associated with a rapid functional decline in patients with sporadic ALS, which is linked to decreased expression of TTN.

Novel DDR2 mutation identified by whole exome sequencing in a Moroccan patient with spondylo-meta-epiphyseal dysplasia, short limb-abnormal calcification type.

Spondylo-meta-epiphyseal dysplasia (SMED), short limb-abnormal calcification type (SMED, SL-AC), is a very rare autosomal recessive disorder with various skeletal changes characterized by premature calcification leading to severe disproportionate short stature. Twenty-two patients have been reported until now, but only five mutations (four missense and one splice-site) in the conserved sequence encoding the tyrosine kinase domain of the DDR2 gene has been identified. We report here a novel DDR2 missense mutation, c.370C > T (p.Arg124Trp) in a Moroccan girl with SMED, SL-AC, identified by whole exome sequencing. Our study has expanded the mutational spectrum of this rare disease and it has shown that exome sequencing is a powerful and cost-effective tool for the diagnosis of clinically heterogeneous disorders such as SMED.

Ectopic expression of Ptf1a induces spinal defects, urogenital defects, and anorectal malformations in Danforth's short tail mice.

Danforth's short tail (Sd) is a semidominant mutation on mouse chromosome 2, characterized by spinal defects, urogenital defects, and anorectal malformations. However, the gene responsible for the Sd phenotype was unknown. In this study, we identified the molecular basis of the Sd mutation. By positional cloning, we identified the insertion of an early transposon in the Sd candidate locus approximately 12-kb upstream of Ptf1a. We found that insertion of the transposon caused overexpression of three neighboring genes, Gm13344, Gm13336, and Ptf1a, in Sd mutant embryos and that the Sd phenotype was not caused by disruption of an as-yet-unknown gene in the candidate locus. Using multiple knockout and knock-in mouse models, we demonstrated that misexpression of Ptf1a, but not of Gm13344 or Gm13336, in the notochord, hindgut, cloaca, and mesonephros was sufficient to replicate the Sd phenotype. The ectopic expression of Ptf1a in the caudal embryo resulted in attenuated expression of Cdx2 and its downstream target genes T, Wnt3a, and Cyp26a1; we conclude that this is the molecular basis of the Sd phenotype. Analysis of Sd mutant mice will provide insight into the development of the spinal column, anus, and kidney.