Language：English   Publishing type：Research paper (scientific journal)  

CEP55 encodes centrosomal protein 55 kDa, which plays a crucial role in mitosis, particularly cytokinesis. Biallelic CEP55 variants cause MARCH syndrome (multinucleated neurons, anhydramnios, renal dysplasia, cerebellar hypoplasia and hydranencephaly). Here, we describe a Japanese family with two affected siblings harboring novel compound heterozygous CEP55 variants, NM_001127182: c.[1357 C > T];[1358 G > A] p.[(Arg453Cys)];[(Arg453His)]. Both presented clinically with typical lethal MARCH syndrome. Although a combination of missense and nonsense variants has been reported previously, this is the first report of biallelic missense CEP55 variants. These variants biallelically affected the same amino acid, Arg453, in the last 40 amino acids of CEP55. These residues are functionally important for CEP55 localization to the midbody during cell division, and may be associated with severe clinical outcomes. More cases of pathogenic CEP55 variants are needed to establish the genotype-phenotype correlation.

DOI： 10.1038/s10038-024-01298-7

PubMed

researchmap

DOI： 10.20944/preprints202409.1911.v1

researchmap

Language：English   Publishing type：Research paper (scientific journal)  

Technologies for detecting structural variation (SV) have advanced with the advent of long-read sequencing, which enables the validation of SV at a nucleotide level. Optical genome mapping (OGM), a technology based on physical mapping, can also provide comprehensive SVs analysis. We applied long-read whole genome sequencing (LRWGS) to accurately reconstruct breakpoint (BP) segments in a patient with complex chromosome 6q rearrangements that remained elusive by conventional karyotyping. Although all BPs were precisely identified by LRWGS, there were two possible ways to construct the BP segments in terms of their orders and orientations. Thus, we also used OGM analysis. Notably, OGM recognized entire inversions exceeding 500 kb in size, which LRWGS could not characterize. Consequently, here we successfully unveil the full genomic structure of this complex chromosomal 6q rearrangement and cryptic SVs through combined long-molecule genomic analyses, showcasing how LRWGS and OGM can complement each other in SV analysis.

DOI： 10.1016/j.ygeno.2024.110894

PubMed

researchmap

Language：English   Publishing type：Research paper (scientific journal)  

SLC5A6 encodes the sodium-dependent multivitamin transporter, a transmembrane protein that uptakes biotin, pantothenic acid, and lipoic acid. Biallelic SLC5A6 variants cause sodium-dependent multivitamin transporter deficiency (SMVTD) and childhood-onset biotin-responsive peripheral motor neuropathy (COMNB), which both respond well to replacement therapy with the above three nutrients. SMVTD usually presents with various symptoms in multiple organs, such as gastrointestinal hemorrhage, brain atrophy, and global developmental delay, at birth or in infancy. Without nutrient replacement therapy, SMVTD can be lethal in early childhood. COMNB is clinically milder and has a later onset than SMVTD, at approximately 10 years of age. COMNB symptoms are mostly limited to peripheral motor neuropathy. Here we report three patients from one Japanese family harboring novel compound heterozygous missense variants in SLC5A6, namely NM_021095.4:c.[221C>T];[642G>C] p.[(Ser74Phe)];[(Gln214His)]. Both variants were predicted to be deleterious through multiple lines of evidence, including amino acid conservation, in silico predictions of pathogenicity, and protein structure considerations. Drosophila analysis also showed c.221C>T to be pathogenic. All three patients had congenital brain cysts on neonatal cranial imaging, but no other morphological abnormalities. They also had a mild motor developmental delay that almost completely resolved despite no treatment. In terms of severity, their phenotypes were intermediate between SMVTD and COMNB. From these findings we propose a new SLC5A6-related disorder, spontaneously remitting developmental delay with brain cysts (SRDDBC) whose phenotypic severity is between that of SMVTD and COMNB. Further clinical and genetic evidence is needed to support our suggestion.

DOI： 10.1038/s10038-023-01206-5

PubMed

researchmap

Language：English   Publishing type：Research paper (scientific journal)  

The gene for ATP binding cassette subfamily A member 2 (ABCA2) is located at chromosome 9q34.3. Biallelic ABCA2 variants lead to intellectual developmental disorder with poor growth and with or without seizures or ataxia (IDPOGSA). In this study, we identified novel compound heterozygous ABCA2 variants (NM_001606.5:c.[5300-17C>A];[6379C>T]) by whole exome sequencing in a 28-year-old Korean female patient with intellectual disability. These variants included intronic and nonsense variants of paternal and maternal origin, respectively, and are absent from gnomAD. SpliceAI predicted that the intron variant creates a cryptic acceptor site. Reverse transcription-PCR using RNA extracted from a lymphoblastoid cell line of the patient confirmed two aberrant transcripts. Her clinical features are compatible with those of IDPOGSA.

DOI： 10.1038/s10038-024-01219-8

PubMed

researchmap

Language：English   Publishing type：Research paper (scientific journal)   Publisher：Springer Science and Business Media LLC  

DOI： 10.1038/s10038-023-01217-2

researchmap

Other Link： https://www.nature.com/articles/s10038-023-01217-2

Authorship：Lead author,　Corresponding author   Publishing type：Research paper (scientific journal)   Publisher：Oxford University Press (OUP)  

Abstract

SARS-CoV-2 is the cause of the current worldwide pandemic of severe acute respiratory syndrome. The change of nucleotide composition of the SARS-CoV-2 genome is crucial for understanding the spread and transmission dynamics of the virus because viral nucleotide sequences are essential in identifying viral strains. Recent studies have shown that cytosine (C) to uracil (U) substitutions are overrepresented in SARS-CoV-2 genome sequences. These asymmetric substitutions between C and U indicate that traditional time-reversible substitution models cannot be applied to the evolution of SARS-CoV-2 sequences. Thus, we develop a new time-irreversible model of nucleotide substitutions to estimate the substitution rates in SARS-CoV-2 genomes. We investigated the number of nucleotide substitutions among the 7862 genomic sequences of SARS-CoV-2 registered in the Global Initiative on Sharing All Influenza Data (GISAID) that have been sampled from all over the world. Using the new method, the substitution rates in SARS-CoV-2 genomes were estimated. The C-to-U substitution rates of SARS-CoV-2 were estimated to be 1.95 × 10−3 ± 4.88 × 10−4 per site per year, compared with 1.48 × 10−4 ± 7.42 × 10−5 per site per year for all other types of substitutions.

DOI： 10.1093/nargab/lqae009

researchmap

Language：English   Publishing type：Research paper (scientific journal)  

DOI： 10.1038/s10038-023-01209-2

PubMed

researchmap

Language：English   Publishing type：Research paper (scientific journal)  

BACKGROUND: Charcot-Marie-Tooth disease type 1A (CMT1A) is one of the most common hereditary peripheral neuropathies caused by duplication of 1.5 Mb genome region including PMP22 gene. We aimed to correct the duplication in human CMT1A patient-derived iPS cells (CMT1A-iPSCs) by genome editing and intended to analyze the effect on Schwann cells differentiated from CMT1A-iPSCs. METHODS: We designed multiple gRNAs targeting a unique sequence present at two sites that sandwich only a single copy of duplicated peripheral myelin protein 22 (PMP22) genes, and selected one of them (gRNA3) from screening their efficiencies by T7E1 mismatch detection assay. AAV2-hSaCas9-gRNAedit was generated by subcloning gRNA3 into pX601-AAV-CMV plasmid, and the genome editing AAV vector was infected to CMT1A-iPSCs or CMT1A-iPSC-derived Schwann cell precursors. The effect of the genome editing AAV vector on myelination was evaluated by co-immunostaining of myelin basic protein (MBP), a marker of mature myelin, and microtubule-associated protein  2(MAP2), a marker of neurites or by electron microscopy. RESULTS: Here we show that infection of CMT1A-iPS cells (iPSCs) with AAV2-hSaCas9-gRNAedit expressing both hSaCas9 and gRNA targeting the tandem repeat sequence decreased PMP22 gene duplication by 20-40%. Infection of CMT1A-iPSC-derived Schwann cell precursors with AAV2-hSaCas9-gRNAedit normalized PMP22 mRNA and PMP22 protein expression levels, and also ameliorated increased apoptosis and impaired myelination in CMT1A-iPSC-derived Schwann cells. CONCLUSIONS: In vivo transfer of AAV2-hSaCas9-gRNAedit to peripheral nerves could be a potential therapeutic modality for CMT1A patient after careful examinations of toxicity including off-target mutations.

DOI： 10.1038/s43856-023-00400-y

PubMed

researchmap

Language：English   Publishing type：Research paper (scientific journal)  

Benign adult familial myoclonic epilepsy type 1 (BAFME1) is an autosomal dominant, adult-onset neurological disease caused by SAMD12 repeat expansion. In BAFME1, anticipation, such as the earlier onset of tremor and/or seizures in the next generation, was reported. This could be explained by intergenerational repeat instability, leading to larger expansions in successive generations. We report a four-generation BAFME1-affected family with anticipation. Using Nanopore long-read sequencing, detailed information regarding the sizes, configurations, and compositions of the expanded SAMD12 repeats across generations was obtained. Unexpectedly, a grandmother-mother-daughter triad showed similar repeat structures but with slight repeat expansions, despite quite variable age of onset of seizures (range: 52-14 years old), implying a complex relationship between the SAMD12 repeat expansion sequence and anticipation. This study suggests that different factor(s) from repeat expansion could modify the anticipation in BAFME1.

DOI： 10.1038/s10038-023-01187-5

PubMed

researchmap

Publishing type：Research paper (scientific journal)   Publisher：Life Science Alliance, LLC  

We discovered biallelic intragenic structural variations (SVs) inFGF12by applying long-read whole genome sequencing to an exome-negative patient with developmental and epileptic encephalopathy (DEE). We also found another DEE patient carrying a biallelic (homozygous) single-nucleotide variant (SNV) inFGF12that was detected by exome sequencing.FGF12heterozygous recurrent missense variants with gain-of-function or heterozygous entire duplication ofFGF12are known causes of epilepsy, but biallelic SNVs/SVs have never been described.FGF12encodes intracellular proteins interacting with the C-terminal domain of the alpha subunit of voltage-gated sodium channels 1.2, 1.5, and 1.6, promoting excitability by delaying fast inactivation of the channels. To validate the molecular pathomechanisms of these biallelicFGF12SVs/SNV, highly sensitive gene expression analyses using lymphoblastoid cells from the patient with biallelic SVs, structural considerations, andDrosophilain vivo functional analysis of the SNV were performed, confirming loss-of-function. Our study highlights the importance of small SVs in Mendelian disorders, which may be overlooked by exome sequencing but can be detected efficiently by long-read whole genome sequencing, providing new insights into the pathomechanisms of human diseases.

DOI： 10.26508/lsa.202302025

PubMed

researchmap

Language：English   Publishing type：Research paper (scientific journal)  

RAC1 at 7p22.1 encodes a RAC family small GTPase that regulates actin cytoskeleton organization and intracellular signaling pathways. Pathogenic RAC1 variants result in developmental delay and multiple anomalies. Here, exome sequencing identified a rare de novo RAC1 variant [NM_018890.4:c.118T > C p.(Tyr40His)] in a male patient. Fetal ultrasonography indicated the patient to have multiple anomalies, including persistent left superior vena cava, total anomalous pulmonary venous return, esophageal atresia, scoliosis, and right-hand polydactyly. After birth, craniofacial dysmorphism and esophagobronchial fistula were confirmed and VACTERL association was suspected. One day after birth, the patient died of respiratory failure caused by tracheal aplasia type III. The molecular mechanisms of pathogenic RAC1 variants remain largely unclear; therefore, we biochemically examined the pathophysiological significance of RAC1-p.Tyr40His by focusing on the best characterized downstream effector of RAC1, PAK1, which activates Hedgehog signaling. RAC1-p.Tyr40His interacted minimally with PAK1, and did not enable PAK1 activation. Variants in the RAC1 Switch II region consistently activate downstream signals, whereas the p.Tyr40His variant at the RAC1-PAK1 binding site and adjacent to the Switch I region may deactivate the signals. It is important to accumulate data from individuals with different RAC1 variants to gain a full understanding of their varied clinical presentations.

DOI： 10.1038/s41598-023-36381-0

PubMed

researchmap

Language：English   Publishing type：Research paper (scientific journal)  

Hereditary spastic paraplegias (HSPs) are a heterogeneous group of neurodegenerative disorders characterized by progressive spasticity and weakness in the lower extremities. To date, a total of 88 types of SPG are known. To diagnose HSP, multiple technologies, including microarray, direct sequencing, multiplex ligation-dependent probe amplification, and short-read next-generation sequencing, are often chosen based on the frequency of HSP subtypes. Exome sequencing (ES) is commonly used. We used ES to analyze ten cases of HSP from eight families. We identified pathogenic variants in three cases (from three different families); however, we were unable to determine the cause of the other seven cases using ES. We therefore applied long-read sequencing to the seven undetermined HSP cases (from five families). We detected intragenic deletions within the SPAST gene in four families, and a deletion within PSEN1 in the remaining family. The size of the deletion ranged from 4.7 to 12.5 kb and involved 1-7 exons. All deletions were entirely included in one long read. We retrospectively performed an ES-based copy number variation analysis focusing on pathogenic deletions, but were not able to accurately detect these deletions. This study demonstrated the efficiency of long-read sequencing in detecting intragenic pathogenic deletions in ES-negative HSP patients.

DOI： 10.1038/s10038-023-01170-0

PubMed

researchmap

Language：English   Publishing type：Research paper (scientific journal)  

Tandem repeats (TRs) are one of the largest sources of polymorphism, and their length is associated with gene regulation. Although previous studies reported several tandem repeats regulating gene splicing in cis (spl-TRs), no large-scale study has been conducted. In this study, we established a genome-wide catalog of 9537 spl-TRs with a total of 58,290 significant TR-splicing associations across 49 tissues (false discovery rate 5%) by using Genotype-Tissue expression (GTex) Project data. Regression models explaining splicing variation by using spl-TRs and other flanking variants suggest that at least some of the spl-TRs directly modulate splicing. In our catalog, two spl-TRs are known loci for repeat expansion diseases, spinocerebellar ataxia 6 (SCA6) and 12 (SCA12). Splicing alterations by these spl-TRs were compatible with those observed in SCA6 and SCA12. Thus, our comprehensive spl-TR catalog may help elucidate the pathomechanism of genetic diseases.

DOI： 10.1101/gr.277335.122

PubMed

researchmap

Language：English   Publishing type：Research paper (scientific journal)   Publisher：Springer Science and Business Media LLC  

Abstract

The Drosophila behavior/human splicing protein family is involved in numerous steps of gene regulation. In humans, this family consists of three proteins: SFPQ, PSPC1, and NONO. Hemizygous loss-of-function (LoF) variants in NONO cause a developmental delay with several complications (e.g., distinctive facial features, cardiac symptoms, and skeletal symptoms) in an X-linked recessive manner. Most of the reported variants have been LoF variants, and two missense variants have been reported as likely deleterious but with no functional validation. We report three individuals from two families harboring an identical missense variant that is located in the nuclear localization signal, NONO: NM_001145408.2:c.1375C > G p.(Pro459Ala). All of them were male and the variant was inherited from their asymptomatic mothers. Individual 1 was diagnosed with developmental delay and cardiac phenotypes (ventricular tachycardia and dilated cardiomyopathy), which overlapped with the features of reported individuals having NONO LoF variants. Individuals 2 and 3 were monozygotic twins. Unlike in Individual 1, developmental delay with autistic features was the only symptom found in them. A fly experiment and cell localization experiment showed that the NONO variant impaired its proper intranuclear localization, leading to mild LoF. Our findings suggest that deleterious NONO missense variants should be taken into consideration when whole-exome sequencing is performed on male individuals with developmental delay with or without cardiac symptoms.

DOI： 10.1038/s41598-023-27770-6

researchmap

Other Link： https://www.nature.com/articles/s41598-023-27770-6

Language：English   Publishing type：Research paper (scientific journal)  

TNNI2 at 11p15.5 encodes troponin I2, fast skeletal type, which is a member of the troponin I gene family and a component of the troponin complex. Distal arthrogryposis (DA) is characterized by congenital limb contractures without primary neurological or muscular effects. DA is inherited in an autosomal dominant fashion and is clinically and genetically heterogeneous. Exome sequencing identified a causative variant in TNNI2 [NM_003282.4:c.532T>C p.(Phe178Leu)] in a Japanese girl with typical DA2b. Interestingly, the familial study using Sanger sequencing suggested a mosaic variant in her healthy father. Subsequent targeted amplicon-based deep sequencing detected the TNNI2 variant with variant allele frequencies of 9.4-17.7% in genomic DNA derived from peripheral blood leukocytes, saliva, hair, and nails in the father. We confirmed a disease-causing variant in TNNI2 in the proband inherited from her asymptomatic father with its somatic variant. Our case demonstrates that careful clinical and genetic evaluation is required in DA.

DOI： 10.1038/s10038-022-01117-x

PubMed

researchmap

Language：English   Publishing type：Research paper (scientific journal)  

Pontocerebellar hypoplasia (PCH) is currently classified into 16 subgroups. Using mostly next-generation sequencing, pathogenic variants have been identified in as many as 24 PCH-associated genes. PCH type 8 (PCH8) is a rare heterogeneous disorder. Its clinical presentation includes severe development delay, increased muscle tone, microcephaly, and magnetic resonance imaging (MRI) abnormalities such as reduced cerebral white matter, a thin corpus callosum, and brainstem and cerebellar hypoplasia. To date, only two variants in the CHMP1A gene (MIM: 164010), NM_002768.5: c.88 C > T (p.Glu30*) and c.28-13 G > A, have been identified homozygously in seven patients with PCH8 from four families (MIM: 614961). CHMP1A is a subunit of the endosomal sorting complex required for transport III (ESCRT-III), which regulates the formation and release of extracellular vesicles. Biallelic CHMP1A loss of function impairs the ESCRT-III-mediated release of extracellular vesicles, which causes impaired progenitor proliferation in the developing brain. Herein, we report a patient with PCH8 who had a homozygous CHMP1A variant, c.122delA (p.Asn41Metfs*2), which arose from segmental uniparental disomy. Although our patient had similar MRI findings to those of previously reported patients, with no progression, we report some novel neurological and developmental findings that expand our knowledge of the clinical consequences associated with CHMP1A variants.

DOI： 10.1038/s10038-022-01098-x

PubMed

researchmap

Language：English   Publishing type：Research paper (scientific journal)  

PURPOSE: Cerebellar hypoplasia and atrophy (CBHA) in children is an extremely heterogeneous group of disorders, but few comprehensive genetic studies have been reported. Comprehensive genetic analysis of CBHA patients may help differentiating atrophy and hypoplasia and potentially improve their prognostic aspects. METHODS: Patients with CBHA in 176 families were genetically examined using exome sequencing. Patients with disease-causing variants were clinically evaluated. RESULTS: Disease-causing variants were identified in 96 of the 176 families (54.5%). After excluding 6 families, 48 patients from 42 families were categorized as having syndromic associations with CBHA, whereas the remaining 51 patients from 48 families had isolated CBHA. In 51 patients, 26 aberrant genes were identified, of which, 20 (76.9%) caused disease in 1 family each. The most prevalent genes were CACNA1A, ITPR1, and KIF1A. Of the 26 aberrant genes, 21 and 1 were functionally annotated to atrophy and hypoplasia, respectively. CBHA+S was more clinically severe than CBHA-S. Notably, ARG1 and FOLR1 variants were identified in 2 families, leading to medical treatments. CONCLUSION: A wide genetic and clinical diversity of CBHA was revealed through exome sequencing in this cohort, which highlights the importance of comprehensive genetic analyses. Furthermore, molecular-based treatment was available for 2 families.

DOI： 10.1016/j.gim.2022.08.007

PubMed

researchmap

Language：English   Publishing type：Research paper (scientific journal)  

We developed a diagnostic method for repeat expansion diseases using a long-read sequencer to improve currently available, low throughput diagnostic methods. We employed the real-time target enrichment system of the nanopore GridION sequencer using the adaptive sampling option, in which software-based target assignment is available without prior sample enrichment, and built an analysis pipeline that prioritized the disease-causing loci. Twenty-two patients with various neurological and neuromuscular diseases, including 12 with genetically diagnosed repeat expansion diseases and 10 manifesting cerebellar ataxia, but without genetic diagnosis, were analyzed. We first sequenced the 12 molecularly diagnosed patients and accurately confirmed expanded repeats in all with uniform depth of coverage across the loci. Next, we applied our method and a conventional method to 10 molecularly undiagnosed patients. Our method corrected inaccurate diagnoses of two patients by the conventional method. Our method is superior to conventional diagnostic methods in terms of speed, accuracy, and comprehensiveness.

DOI： 10.1038/s41525-022-00331-y

PubMed

researchmap

Language：English   Publishing type：Research paper (scientific journal)   Publisher：Wiley  

<jats:title>Abstract</jats:title><jats:p>Recent advances in sequencing technologies enable genome‐wide analyses for thousands of individuals. The sequential kernel association test (SKAT) is a widely used method to test for associations between a phenotype and a set of rare variants. As the sample size of human genetics studies increases, the computational time required to calculate a kernel is becoming more and more problematic. In this study, a new method to obtain kernel statistics without calculating a kernel matrix is proposed. A simple method for the computation of two kernel statistics, namely, a kernel statistic based on a genetic relationship matrix (GRM) and one based on an identity by state (IBS) matrix, are proposed. By using this method, calculation of the kernel statistics can be conducted using vector calculation without matrix calculation. The proposed method enables one to conduct SKAT for large samples of human genetics.</jats:p>

DOI： 10.1002/ggn2.202100066

PubMed

researchmap

Language：English   Publishing type：Research paper (scientific journal)  

We report two patients with autosomal dominant neuronal intranuclear inclusion disease (NIID) harboring the biallelic GGC repeat expansion in NOTCH2NLC to uncover the impact of repeat expansion zygosity on the clinical phenotype. The zygosity of the entire NOTCH2NLC GGC repeat expansion and DNA methylation were comprehensively evaluated using fluorescent amplicon length PCR (AL-PCR), Southern blotting and targeted long-read sequencing, and detailed genetic/epigenetic and clinical features were described. In AL-PCR, we could not recognize the wild-type allele in both patients. Targeted long-read sequencing revealed that one patient harbored a homozygous repeat expansion. The other patient harbored compound heterozygous repeat expansions. The GGC repeats and the nearest CpG island were hypomethylated in all expanded alleles in both patients. Both patients harboring the biallelic GGC repeat expansion showed a typical dementia-dominant NIID phenotype. In conclusion, the biallelic GGC repeat expansion in two typical NIID patients indicated that NOTCH2NLC-related diseases could be completely dominant.

DOI： 10.1016/j.ygeno.2022.110469

PubMed

researchmap

Language：English   Publishing type：Research paper (scientific journal)  

Recent studies suggest that transcript isoforms significantly overlap (approximately 60%) between brain tissue and Epstein-Barr virus-transformed lymphoblastoid cell lines (LCLs). Interestingly, 14 cohesion-related genes with variants that cause Cornelia de Lange Syndrome (CdLS) are highly expressed in the brain and LCLs. In this context, we first performed RNA sequencing of LCLs from 22 solved (with pathogenic variants) and 19 unsolved (with no confirmed variants) CdLS cases. Next, an RNA sequencing pipeline was developed using solved cases with two different methods: short variant analysis (for single-nucleotide and indel variants) and aberrant splicing detection analysis. Then, 19 unsolved cases were subsequently applied to our pipeline, and four pathogenic variants in NIPBL (one inframe deletion and three intronic variants) were newly identified. Two of three intronic variants were located at Alu elements in deep-intronic regions, creating cryptic exons. RNA sequencing with LCLs was useful for identifying hidden variants in exome-negative cases.

DOI： 10.1016/j.ygeno.2022.110468

PubMed

researchmap

Language：English   Publishing type：Research paper (scientific journal)  

AIMS: The molecular genetic mechanisms underlying postoperative nausea and vomiting (PONV) in the brain have not been fully elucidated. This study aimed to determine the changes in whole transcriptome in the nucleus of the solitary tract (NTS) in an animal model of PONV, to screen a drug candidate and to elucidate the molecular genetic mechanisms of PONV development. METHODS: Twenty-one female musk shrews were assigned into three groups: the Surgery group (shrew PONV model, n = 9), the Sham group (n = 6), and the Naïve group (n = 6). In behavioral studies, the main outcome was the number of emetic episodes. In genetic experiments, changes in the transcriptome in the NTS were measured. In a separate study, 12 shrews were used to verify the candidate mechanism underlying PONV. RESULTS: A median of six emetic episodes occurred in both the Sham and Surgery groups. Whole-transcriptome analysis indicated the inhibition of the GABAB receptor-mediated signaling pathway in the PONV model. Baclofen (GABAB receptor agonist) administration eliminated emetic behaviors in the shrew PONV model. CONCLUSIONS: Our findings suggest that the GABAB receptor-mediated signaling pathway is involved in emesis and that baclofen may be a novel therapeutic or prophylactic agent for PONV.

DOI： 10.1111/cns.13823

PubMed

researchmap

Language：English   Publishing type：Research paper (scientific journal)   Publisher：{MDPI} {AG}  

Hypouricemia is recognized as a rare disorder, defined as a serum uric acid level of 2.0 mg/dL or less. Hypouricemia is divided into an overexcretion type and an underproduction type. The former typical disease is xanthinuria, and the latter is renal hypouricemia (RHUC). The frequency of nephrogenic hypouricemia due to a deficiency of URAT1 is high in Japan, accounting for most asymptomatic and persistent cases of hypouricemia. RHUC results in a high risk of exercise-induced acute kidney injury and urolithiasis. It is vital to promote research on RHUC, as this will lead not only to the elucidation of its pathophysiology but also to the development of new treatments for gout and hyperuricemia.

DOI： 10.3390/biomedicines10030652

PubMed

researchmap

Publishing type：Research paper (scientific journal)   Publisher：Cold Spring Harbor Laboratory  

<jats:title>Abstract</jats:title><jats:p>Principal component analyses (PCAs) are often used to visualize patterns of genetic variation in human populations. Previous studies showed a close correspondence between genetic and geographic distances. In such PCAs, the principal components are eigenvectors of the data’s variance-covariance matrix, which is obtained by a genetic relationship matrix (GRM). However, it is difficult to apply GRM to multiallelic sites. In this paper, I showed that a PCA from GRM is equivalent to multidimensional scaling (MDS) from nucleotide differences. Therefore, a PCA can be conducted using nucleotide differences. The new method provided in this study provides a straightforward method to predict the effects of different demographic processes on genetic diversity.</jats:p>

DOI： 10.1101/2022.03.02.482744

researchmap

Language：English   Publishing type：Research paper (scientific journal)  

To reveal gene-environment interactions underlying common diseases and estimate the risk for common diseases, the Tohoku Medical Megabank (TMM) project has conducted prospective cohort studies and genomic and multiomics analyses. To establish an integrated biobank, we developed an integrated database called “dbTMM” that incorporates both the individual cohort/clinical data and the genome/multiomics data of 157,191 participants in the Tohoku Medical Megabank project. To our knowledge, dbTMM is the first database to store individual whole-genome data on a variant-by-variant basis as well as cohort/clinical data for over one hundred thousand participants in a prospective cohort study. dbTMM enables us to stratify our cohort by both genome-wide genetic factors and environmental factors, and it provides a research and development platform that enables prospective analysis of large-scale data from genome cohorts.

DOI： 10.1038/s41439-021-00175-5

Scopus

researchmap

Authorship：Lead author,　Corresponding author   Language：Japanese   Publishing type：Research paper (scientific journal)  

DOI： 10.14867/gnamtsunyo.45.1_23

researchmap

Language：English   Publishing type：Research paper (scientific journal)   Publisher：Nature Publishing Group  

Amyotrophic lateral sclerosis (ALS) is a devastating progressive motor neuron disease that affects people of all ethnicities. Approximately 90% of ALS cases are sporadic and thought to have multifactorial pathogenesis. To understand the genetics of sporadic ALS, we conducted a genome-wide association study using 1,173 sporadic ALS cases and 8,925 controls in a Japanese population. A combined meta-analysis of our Japanese cohort with individuals of European ancestry revealed a significant association at the ACSL5 locus (top SNP p = 2.97 × 10-8). We validated the association with ACSL5 in a replication study with a Chinese population and an independent Japanese population (1941 ALS cases, 3821 controls; top SNP p = 1.82 × 10-4). In the combined meta-analysis, the intronic ACSL5 SNP rs3736947 showed the strongest association (p = 7.81 × 10-11). Using a gene-based analysis of the full multi-ethnic dataset, we uncovered additional genes significantly associated with ALS: ERGIC1, RAPGEF5, FNBP1, and ATXN3. These results advance our understanding of the genetic basis of sporadic ALS.

DOI： 10.1038/s42003-020-01251-2

researchmap

Language：English   Publishing type：Research paper (scientific journal)  

Over the past few decades, obesity has become a public health issue of global concern. Even though disparities exist between human populations, e.g., the higher liver fat content of the Japanese despite a lower body mass index (BMI), studies on the genetics of obesity still largely focus on populations of European descent, leading to a dearth of genetic data on non-European populations. In this context, this study aimed to establish a broad picture of the genetic attributes of the Japanese population, by examining a representative sample of 18,889 individuals participating in the Tohoku Medical Megabank Project cohort. We applied linear mixed model methods to 17 traits related to obesity and associated diseases to estimate the heritabilities explained by common genetic variants and the genetic correlations between each pair of traits. These analyses allowed us to quantify the SNP heritability of health indicators such as BMI (0.248 ± 0.032) and HDL cholesterol (0.324 ± 0.031), and to provide one of the few estimates of the SNP heritability of cystatin C in unrelated individuals (0.260 ± 0.025). We discuss potential differences between the Japanese and people of European ancestry with respect to the genetic correlations between urinary biomarkers and adiposity traits, for which large estimates were obtained. For instance, the genetic correlations between urine potassium level and the values for weight, BMI, waist circumference, and waist-to-height ratio ranged from 0.290 to 0.559, much higher than the corresponding estimates in the UK Biobank.

DOI： 10.1534/g3.120.401242

PubMed

researchmap

Authorship：Lead author,　Corresponding author   Language：English   Publishing type：Research paper (scientific journal)  

Gout is a common arthritis caused by monosodium urate crystals. The heritability of serum urate levels is estimated to be 30‒70%; however, common genetic variants account for only 7.9% of the variance in serum urate levels. This discrepancy is an example of "missing heritability." The "missing heritability" suggests that variants associated with uric acid levels are yet to be found. By using genomic sequences of the ToMMo cohort, we identified rare variants of the SLC22A12 gene that affect the urate transport activity of URAT1. URAT1 is a transporter protein encoded by the SLC22A12 gene. We grouped the participants with variants affecting urate uptake by URAT1 and analyzed the variance of serum urate levels. The results showed that the heritability explained by the SLC22A12 variants of men and women exceeds 10%, suggesting that rare variants underlie a substantial portion of the "missing heritability" of serum urate levels.

DOI： 10.1534/genetics.119.303006

PubMed

researchmap

Authorship：Lead author,　Corresponding author   Language：English   Publishing type：Research paper (scientific journal)   Publisher：Springer Science and Business Media {LLC}  

To explore the effect of aging, a cohort study is being performed on Sado Island, which is located in the Sea of Japan. Sado Island is close to the eastern coast of Japan, yet its population speaks the western Japanese dialect. Consequently, the genetic background of the population of Sado Island is of interest. Based on Nei's genetic distance, we compared the allele frequencies of people from Sado Island to those of people from Nagahama and Miyagi, which are located in the western and northeastern parts of Honshu, respectively. The results showed that the populations of Miyagi and Nagahama are genetically closer to each other than to the population of Sado Island. Because the Sado and Honshu Islands are isolated by a channel, it is possible that genetic drift occurred within Sado Island, which would explain the uniqueness of the people of this region.

DOI： 10.1038/s41439-019-0058-6

PubMed

researchmap

Language：English   Publishing type：Research paper (scientific journal)  

In recent genome analyses, population-specific reference panels have indicated important. However, reference panels based on short-read sequencing data do not sufficiently cover long insertions. Therefore, the nature of long insertions has not been well documented. Here, we assembled a Japanese genome using single-molecule real-time sequencing data and characterized insertions found in the assembled genome. We identified 3691 insertions ranging from 100 bps to ~10,000 bps in the assembled genome relative to the international reference sequence (GRCh38). To validate and characterize these insertions, we mapped short-reads from 1070 Japanese individuals and 728 individuals from eight other populations to insertions integrated into GRCh38. With this result, we constructed JRGv1 (Japanese Reference Genome version 1) by integrating the 903 verified insertions, totaling 1,086,173 bases, shared by at least two Japanese individuals into GRCh38. We also constructed decoyJRGv1 by concatenating 3559 verified insertions, totaling 2,536,870 bases, shared by at least two Japanese individuals or by six other assemblies. This assembly improved the alignment ratio by 0.4% on average. These results demonstrate the importance of refining the reference assembly and creating a population-specific reference genome. JRGv1 and decoyJRGv1 are available at the JRG website.

DOI： 10.1038/s41439-019-0057-7

PubMed

researchmap

Language：English   Publishing type：Research paper (scientific journal)   Publisher：Springer Nature  

DOI： 10.1038/s41397-017-0010-4

Scopus

researchmap

Language：English   Publishing type：Research paper (scientific journal)   Publisher：Elsevier {BV}  

Forensic DNA typing is widely used to identify missing persons and plays a central role in forensic profiling. DNA typing usually uses capillary electrophoresis fragment analysis of PCR amplification products to detect the length of short tandem repeat (STR) markers. Here, we analyzed whole genome data from 1,070 Japanese individuals generated using massively parallel short-read sequencing of 162 paired-end bases. We have analyzed 843,473 STR loci with two to six basepair repeat units and cataloged highly polymorphic STR loci in the Japanese population. To evaluate the performance of the cataloged STR loci, we compared 23 STR loci, widely used in forensic DNA typing, with capillary electrophoresis based STR genotyping results in the Japanese population. Seventeen loci had high correlations and high call rates. The other six loci had low call rates or low correlations due to either the limitations of short-read sequencing technology, the bioinformatics tool used, or the complexity of repeat patterns. With these analyses, we have also purified the suitable 218 STR loci with four basepair repeat units and 53 loci with five basepair repeat units both for short read sequencing and PCR based technologies, which would be candidates to the actual forensic DNA typing in Japanese population.

DOI： 10.1016/j.heliyon.2018.e00625

Scopus

PubMed

researchmap

Language：English   Publishing type：Research paper (scientific journal)  

Background<br />
In the estimation of repeat numbers in a short tandem repeat (STR) region from high-throughput sequencing data, two types of strategies are mainly taken: a strategy based on counting repeat patterns included in sequence reads spanning the region and a strategy based on estimating the difference between the actual insert size and the insert size inferred from paired-end reads. The quality of sequence alignment is crucial, especially in the former approaches although usual alignment methods have difficulty in STR regions due to insertions and deletions caused by the variations of repeat numbers.<br />
<br />
Results<br />
We proposed a new dynamic programming based realignment method named STR-realigner that considers repeat patterns in STR regions as prior knowledge. By allowing the size change of repeat patterns with low penalty in STR regions, accurate realignment is expected. For the performance evaluation, publicly available STR variant calling tools were applied to three types of aligned reads: synthetically generated sequencing reads aligned with BWA-MEM, those realigned with STR-realigner, those realigned with ReviSTER, and those realigned with GATK IndelRealigner. From the comparison of root mean squared errors between estimated and true STR region size, the results for the dataset realigned with STR-realigner are better than those for other cases. For real data analysis, we used a real sequencing dataset from Illumina HiSeq 2000 for a parent-offspring trio. RepeatSeq and lobSTR were applied to the sequence reads for these individuals aligned with BWA-MEM, those realigned with STR-realigner, ReviSTER, and GATK IndelRealigner. STR-realigner shows the best performance in terms of consistency of the size of estimated STR regions in Mendelian inheritance. Root mean squared error values were also calculated from the comparison of these estimated results with STR region sizes obtained from high coverage PacBio sequencing data, and the results from the realigned sequencing data with STR-realigner showed the least (the best) root mean squared error value.<br />
<br />
Conclusions<br />
The effectiveness of the proposed realignment method for STR regions was verified from the comparison with an existing method on both simulation datasets and real whole genome sequencing dataset.

DOI： 10.1186/s12864-016-3294-x

Scopus

PubMed

researchmap

Language：English   Publishing type：Research paper (scientific journal)  

Background<br />
Genome-wide association studies have revealed associations between single-nucleotide polymorphisms (SNPs) and phenotypes such as disease symptoms and drug tolerance. To address the small sample size for rare variants, association studies tend to group gene or pathway level variants and evaluate the effect on the set of variants. One of such strategies, known as the sequential kernel association test (SKAT), is a widely used collapsing method. However, the reported p-values from SKAT tend to be biased because the asymptotic property of the statistic is used to calculate the p-value. Although this bias can be corrected by applying permutation procedures for the test statistics, the computational cost of obtaining p-values with high resolution is prohibitive.<br />
<br />
Results<br />
To address this problem, we devise an adaptive SKAT procedure termed AP-SKAT that efficiently classifies significant SNP sets and ranks them according to the permuted p-values. Our procedure adaptively stops the permutation test when the significance level is outside some confidence interval of the estimated p-value for a binomial distribution. To evaluate the performance, we first compare the power and sample size calculation and the type I error rates estimate of SKAT, SKAT-O, and the proposed procedure using genotype data in the SKAT R package and from 1000 Genome Project. Through computational experiments using whole genome sequencing and SNP array data, we show that our proposed procedure is highly efficient and has comparable accuracy to the standard procedure.<br />
<br />
Conclusions<br />
For several types of genetic data, the developed procedure could achieve competitive power and sample size under small and large sample size conditions with controlling considerable type I error rates, and estimate p-values of significant SNP sets that are consistent with those estimated by the standard permutation test within a realistic time. This demonstrates that the procedure is sufficiently powerful for recent whole genome sequencing and SNP array data with increasing numbers of phenotypes. Additionally, this procedure can be used in other association tests by employing alternative methods to calculate the statistics.

DOI： 10.1186/s12864-016-3094-3

Scopus

PubMed

researchmap

Authorship：Lead author,　Last author,　Corresponding author   Language：English   Publishing type：Research paper (scientific journal)   Publisher：Elsevier {BV}  

Repetitive sequences in a genome cause mapping errors, especially in case of short reads, because of the presence of similar or identical sequences. Distribution of repetitive sequences in a genome must be studied to distinguish between mappable and unmappable regions. Previous studies showed that short tandem repeats (STRs) are clustered in the terminal regions of chromosomes in the human genome. It is an open question whether formation of STRs in the terminal regions of chromosomes occurs only in humans. The present study investigated the distribution of STRs in the human, cow, mouse, chicken, and lizard genomes. In this study, it was shown that STRs were concentrated in the terminal regions of chromosomes not only in the human genome, but also in mouse, cow, chicken, and lizard genomes. The results suggested that the mechanism through which STRs are shared by amniotes in which mammals, birds and lizards are included. Thus, we must be careful when the genomic sequences at chromosome ends of amniotes by using the next generation sequencers.

DOI： 10.1016/j.genrep.2016.01.004

Scopus

researchmap

Language：Japanese   Publishing type：Research paper (conference, symposium, etc.)  

researchmap

researchmap

Language：Japanese   Publishing type：Research paper (conference, symposium, etc.)  

researchmap

Authorship：Lead author,　Corresponding author   Language：English   Publishing type：Research paper (scientific journal)  

DNA sequencers that can conduct real-time sequencing from a single polymerase molecule are known as third-generation sequencers. Third-generation sequencers enable sequencing of reads that are several kilobases long. However, the raw data generated from third-generation sequencers are known to be error-prone. Because of sequencing errors, it is difficult to identify which genes are homologous to the reads obtained using third-generation sequencers. In this study, a new method for homology search algorithm, PAFFT, is developed. This method is the extension of the MAFFT algorithm which was used for multiple alignments. PAFFT detects global homology rather than local homology so that homologous regions can be detected even when the error rate of sequencing is high. PAFFT will boost application of third-generation sequencers.

DOI： 10.1016/j.ygeno.2015.09.005

PubMed

researchmap

Authorship：Lead author,　Last author,　Corresponding author   Language：English   Publishing type：Research paper (scientific journal)  

Programmable sequence-specific endonucleases are powerful tools for<br />
genome alteration with high precision. For example, the CRISPR system is an<br />
efficient tool for genome engineering in eukaryotic cells by simply specifying a<br />
20-bp targeting sequence within its guide RNA. When studying large genomes,<br />
however, the design of target sequences is complicated by the redundancy of<br />
sequences. The distribution of unique sequences in the genome is of interest.<br />
In this paper, I describe the development of a novel method, UF, for detecting<br />
unique 20-bp sequences in entire genomes. UF stands for “Unique Finder”.<br />
By using UF, the distribution of unique sequences in the human genome<br />
was investigated. It was found that 60% of the human genome is unique on<br />
average. However, non-unique regions of human genome are concentrated<br />
on centromeres and terminal regions of the chromosomes. The proportions of<br />
unique sequences are about 80% in the rest part of the genome. The program<br />
for obtaining unique sequences is available at https://sourceforge.jp/projects/<br />
parallelgwas/releases/

researchmap

Authorship：Lead author,　Last author,　Corresponding author   Language：English   Publishing type：Research paper (scientific journal)  

Next-generation sequencing platforms generate short (50-150 bp) reads that can be mapped onto the reference genome. Repetitive sequences in the genome, because of the presence of similar or identical sequences, cause mapping errors in the case of the short reads. By filtering short reads with repeats, mapping will be improved. I developed RF. RF is a new method that filters short reads with tandem repeats. A scoring scheme was developed that assigned higher scores to regions with tandem repeats and lower scores to regions without tandem repeats. In this study, RF was applied to filter out short reads with repeats, before short reads were mapped onto the same genomic contig by using a short read-mapping program. The result suggests RF improved the proportion of correctly mapped short reads on filtering the repeats. RF is a useful tool for reducing mapping errors of short reads onto reference genomes.

DOI： 10.1016/j.ygeno.2013.03.002

Scopus

researchmap

Authorship：Lead author,　Corresponding author   Language：English   Publishing type：Research paper (scientific journal)  

Efficient determination of evolutionary distances is important for the correct reconstruction of phylogenetic trees. The performance of the pooled distance required for reconstructing a phylogenetic tree can be improved by applying large weights to appropriate distances for reconstructing phylogenetic trees and small weights to inappropriate distances. We developed two weighting methods, the modified Tajima–Takezaki method and the modified least-squares method, for reconstructing phylogenetic trees from multiple loci. By computer simulations, we found that both of the new methods were more efficient in reconstructing correct topologies than the no-weight method. Hence, we reconstructed hominoid phylogenetic trees from mitochondrial DNA using our new methods, and found that the levels of bootstrap support were significantly increased by the modified Tajima–Takezaki and by the modified least-squares method.

DOI： 10.1007/s00239-012-9513-4

Scopus

PubMed

researchmap

Authorship：Lead author,　Last author,　Corresponding author  

DOI： 10.5772/34923

researchmap

Language：English   Publishing type：Research paper (scientific journal)  

Many of the genes that control photosynthesis are carried in the chloroplast. These genes differ among species. However, evidence has yet to be reported revealing the involvement of organelle genes in the initial stages of plant speciation. To elucidate the molecular basis of aquatic plant speciation, we focused on the unique plant species Chara braunii C. C. Gmel. that inhabits both shallow and deep freshwater habitats and exhibits habitat‐based dimorphism of chloroplast DNA (cpDNA). Here, we examined the “shallow” and “deep” subpopulations of C. braunii using two nuclear DNA (nDNA) markers and cpDNA. Genetic differentiation between the two subpopulations was measured in both nDNA and cpDNA regions, although phylogenetic analyses suggested nuclear gene flow between subpopulations. Neutrality tests based on Tajima’s D demonstrated diversifying selection acting on organelle DNA regions. Furthermore, both “shallow” and “deep” haplotypes of cpDNA detected in cultures originating from bottom soils of three deep environments suggested that migration of oospores (dormant zygotes) between the two habitats occurs irrespective of the complete habitat‐based dimorphism of cpDNA from field‐collected vegetative thalli. Therefore, the two subpopulations are highly selected by their different aquatic habitats and show prezygotic isolation, which represents an initial process of speciation affected by ecologically based divergent selection of organelle genes.

DOI： 10.1111/j.1529-8817.2011.01037.x

researchmap

Authorship：Lead author,　Corresponding author   Language：English   Publishing type：Research paper (scientific journal)  

Background<br />
Codon bias is a phenomenon that refers to the differences in the frequencies of synonymous codons among different genes. In many organisms, natural selection is considered to be a cause of codon bias because codon usage in highly expressed genes is biased toward optimal codons. Methods have previously been developed to predict the expression level of genes from their nucleotide sequences, which is based on the observation that synonymous codon usage shows an overall bias toward a few codons called major codons. However, the relationship between codon bias and gene expression level, as proposed by the translation-selection model, is less evident in mammals.<br />
<br />
Findings<br />
We investigated the correlations between the expression levels of 1,182 mouse genes and amino acid composition, as well as between gene expression and codon preference. We found that a weak but significant correlation exists between gene expression levels and amino acid composition in mouse. In total, less than 10% of variation of expression levels is explained by amino acid components. We found the effect of codon preference on gene expression was weaker than the effect of amino acid composition, because no significant correlations were observed with respect to codon preference.<br />
<br />
Conclusion<br />
These results suggest that it is difficult to predict expression level from amino acid components or from codon bias in mouse.

DOI： 10.1186/1756-0500-4-20

researchmap

Authorship：Lead author,　Corresponding author   Language：English   Publishing type：Research paper (scientific journal)  

BACKGROUND: Use of missing genotype imputations and haplotype reconstructions are valuable in genome-wide association studies (GWASs). By modeling the patterns of linkage disequilibrium in a reference panel, genotypes not directly measured in the study samples can be imputed and used for GWASs. Since millions of single nucleotide polymorphisms need to be imputed in a GWAS, faster methods for genotype imputation and haplotype reconstruction are required. RESULTS: We developed a program package for parallel computation of genotype imputation and haplotype reconstruction. Our program package, ParaHaplo 3.0, is intended for use in workstation clusters using the Intel Message Passing Interface. We compared the performance of ParaHaplo 3.0 on the Japanese in Tokyo, Japan and Han Chinese in Beijing, and Chinese in the HapMap dataset. A parallel version of ParaHaplo 3.0 can conduct genotype imputation 20 times faster than a non-parallel version of ParaHaplo. CONCLUSIONS: ParaHaplo 3.0 is an invaluable tool for conducting haplotype-based GWASs. The need for faster genotype imputation and haplotype reconstruction using parallel computing will become increasingly important as the data sizes of such projects continue to increase. ParaHaplo executable binaries and program sources are available at http://en.sourceforge.jp/projects/parallelgwas/releases/.

DOI： 10.1186/1751-0473-6-10

researchmap

Language：English   Publishing type：Research paper (scientific journal)  

DOI： 10.1186/1471-2164-12-397

researchmap

Authorship：Lead author,　Corresponding author   Language：English   Publishing type：Research paper (scientific journal)  

BACKGROUND: The use of haplotype-based association tests can improve the power of genome-wide association studies. Since the observed genotypes are unordered pairs of alleles, haplotype phase must be inferred. However, estimating haplotype phase is time consuming. When millions of single-nucleotide polymorphisms (SNPs) are analyzed in genome-wide association study, faster methods for haplotype estimation are required. METHODS: We developed a program package for parallel computation of haplotype estimation. Our program package, ParaHaplo 2.0, is intended for use in workstation clusters using the Intel Message Passing Interface (MPI). We compared the performance of our algorithm to that of the regular permutation test on both Japanese in Tokyo, Japan and Han Chinese in Beijing, China of the HapMap dataset. RESULTS: Parallel version of ParaHaplo 2.0 can estimate haplotypes 100 times faster than a non-parallel version of the ParaHaplo. CONCLUSION: ParaHaplo 2.0 is an invaluable tool for conducting haplotype-based genome-wide association studies (GWAS). The need for fast haplotype estimation using parallel computing will become increasingly important as the data sizes of such projects continue to increase. The executable binaries and program sources of ParaHaplo are available at the following address: http://en.sourceforge.jp/projects/parallelgwas/releases/

DOI： 10.1186/1751-0473-5-5

researchmap

Language：English   Publishing type：Research paper (scientific journal)  

DOI： 10.1186/1756-0381-3-6

researchmap

Authorship：Lead author,　Corresponding author   Language：English   Publishing type：Research paper (scientific journal)  

Background: Since more than a million single-nucleotide polymorphisms (SNPs) are analyzed in any given genome-wide association study (GWAS), performing multiple comparisons can be problematic. To cope with multiple-comparison problems in GWAS, haplotype-based algorithms were developed to correct for multiple comparisons at multiple SNP loci in linkage disequilibrium. A permutation test can also control problems inherent in multiple testing
however, both the calculation of exact probability and the execution of permutation tests are time-consuming. Faster methods for calculating exact probabilities and executing permutation tests are required. Methods: We developed a set of computer programs for the parallel computation of accurate P-values in haplotype-based GWAS. Our program, ParaHaplo, is intended for workstation clusters using the Intel Message Passing Interface (MPI). We compared the performance of our algorithm to that of the regular permutation test on JPT and CHB of HapMap. Results: ParaHaplo can detect smaller differences between 2 populations than SNP-based GWAS. We also found that parallel-computing techniques made ParaHaplo 100-fold faster than a non-parallel version of the program. Conclusion: ParaHaplo is a useful tool in conducting haplotype-based GWAS. Since the data sizes of such projects continue to increase, the use of fast computations with parallel computing--such as that used in ParaHaplo--will become increasingly important. The executable binaries and program sources of ParaHaplo are available at the following address: http://sourceforge.jp/projects/parallelgwas/?_sl=1. © 2009 Misawa and Kamatani
licensee BioMed Central Ltd.

DOI： 10.1186/1751-0473-4-7

Scopus

researchmap

Language：English   Publishing type：Research paper (scientific journal)  

DOI： 10.1016/j.ympev.2009.08.015

PubMed

researchmap

Language：English   Publishing type：Research paper (scientific journal)  

DOI： 10.1016/j.gene.2008.11.006

researchmap

Language：English   Publishing type：Research paper (scientific journal)   Publisher：Springer Science and Business Media LLC  

DOI： 10.1186/1471-2148-9-197

Web of Science

researchmap

Language：English   Publishing type：Research paper (scientific journal)   Publisher：Springer Science and Business Media LLC  

DOI： 10.1186/1471-2148-9-197

Web of Science

PubMed

researchmap

Language：English   Publishing type：Research paper (scientific journal)   Publisher：ACADEMIC PRESS INC ELSEVIER SCIENCE  

DOI： 10.1016/j.ympev.2008.08.018

Web of Science

Scopus

researchmap

Language：English   Publishing type：Research paper (scientific journal)  

DOI： 10.1016/j.ympev.2008.03.016

PubMed

researchmap

Language：English   Publishing type：Research paper (scientific journal)  

DOI： 10.1016/j.gene.2008.02.021

PubMed

researchmap

Language：English   Publishing type：Research paper (scientific journal)  

DOI： 10.1186/1471-2148-8-151

researchmap

Language：English   Publishing type：Research paper (scientific journal)  

DOI： 10.1007/s10038-008-0312-0

PubMed

researchmap

Language：English   Publishing type：Research paper (scientific journal)  

DOI： 10.1007/s00239-008-9141-1

researchmap

Language：English   Publisher：BLACKWELL PUBLISHING  

Web of Science

researchmap

Language：English   Publishing type：Research paper (scientific journal)   Publisher：OXFORD UNIV PRESS  

DOI： 10.1093/molbev/msm091

Web of Science

PubMed

researchmap

Language：Japanese   Publishing type：Research paper (scientific journal)  

Multiple sequence alignment is an important tool for computational analysis of nucleotide or amino acid sequences. It is also a challenging combinatorial optimization problem in computer science. As a large amount of sequence data is becoming available from genome and other large-scale sequencing projects, efficiency, as well as accuracy, is currently required for a multiple sequence alignment program. Several new programs are being developed aiming at improving both efficiency and accuracy. We overview the algorithms and performances of new programs including that by ourselves

DOI： 10.2142/biophys.46.312

researchmap

Language：English   Publishing type：Research paper (scientific journal)  

DOI： 10.1093/dnares/dsi001

PubMed

researchmap

Language：English   Publishing type：Research paper (scientific journal)  

DOI： 10.1006/jtbi.2003.3144

researchmap

Language：English   Publishing type：Research paper (scientific journal)  

DOI： 10.1016/S1055-7903(03)00079-4

researchmap

Language：English   Publishing type：Research paper (scientific journal)  

DOI： 10.1007/s00239-003-0039-7

researchmap

Language：English   Publishing type：Research paper (scientific journal)  

DOI： 10.1093/nar/gkf436

Scopus

PubMed

researchmap

Authorship：Lead author   Language：English   Publishing type：Research paper (scientific journal)   Publisher：SOC MOLECULAR BIOLOGY EVOLUTION  

A new simple method for classifying genes is proposed based on Klastorin's method. This method classifies genes into monophyletic groups which are made distinct from each other by evolutionary changes. The method is applicable as long as the phylogenetic tree of genes is obtained. There is a fast algorithm for obtaining the classification. A bootstrap test of a classification is also presented. As an example, we classified opsin genes. The classification obtained by this method is the same as the previous classification based on the function of opsins.

DOI： 10.1093/oxfordjournals.molbev.a026289

Web of Science

Scopus

PubMed

researchmap

Language：English   Publishing type：Research paper (scientific journal)   Publisher：ACADEMIC PRESS INC  

A combined data set of DNA sequences (6021 bp) from five protein-coding genes of the chloroplast genome (rbcL, atpB, psaA, psaB, and psbC genes) were analyzed for 42 strains representing 30 species of the colonial Volvocales (Volvox and its relatives) and 5 related species of green algae to deduce robust phylogenetic relationships within the colonial green flagellates. The 4-celled family Tetrabaenaceae was robustly resolved as the most basal group within the colonial Volvocales. The sequence data also suggested that all five volvocacean genera with 32 or more cells in a vegetative colony (all four of the anisogamous/oogamous genera, Eudorina, Platydorina, Pleodorina, and Volvox, plus the isogamous genus Yamagishiella) constituted a large monophyletic group, in which 2 Pleodorina species were positioned distally to 3 species of Volvox. Therefore, most of the evolution of the colonial Volvocales appears to constitute a gradual progression in colonial complexity and in types of sexual reproduction, as in the traditional volvocine lineage hypothesis, although reverse evolution must be considered for the origin of certain species of Pleodorina. Data presented here also provide robust support for a monophyletic family Goniaceae consisting of two genera: Gonium and Astrephomene. (C) 2000 Academic Press.

DOI： 10.1006/mpev.2000.0831

Web of Science

PubMed

researchmap

Language：English   Publishing type：Research paper (scientific journal)   Publisher：KLUWER ACADEMIC PUBL  

The amount and pattern of DNA polymorphism can give useful information on the maintenance mechanism of genetic variation at the DNA level. In this note we have shown the amount and pattern of DNA polymorphism expected under the neutral theory. The amount of DNA polymorphism can be estimated from the average number of nucleotide differences per site, the proportion of segregating sites, and so on. We have shown how to estimate theta from these quantities, where theta = 4N nu,N is the effective population size and nu is the mutation rate per site per generation. We have also shown the expectations of the nucleotide variation within and between allelic classes.

DOI： 10.1023/A:1017011631572

Web of Science

Scopus

researchmap

Publishing type：Research paper (scientific journal)   Publisher：Oxford University Press (OUP)  

Abstract

Knowing the amount of DNA polymorphism is essential to understand the mechanism of maintaining DNA polymorphism in a natural population. The amount of DNA polymorphism can be measured by the average number of nucleotide differences per site (π), the proportion of segregating (polymorphic) site (s) and the minimum number of mutations per site (s*). Since the latter two quantities depend on the sample size, θ is often used as a measure of the amount of DNA polymorphism, where θ = 4Nμ, N is the effective population size and μ is the neutral mutation rate per site per generation. It is known that θ estimated from π, s and s* under the infinite site model can be biased when the mutation rate varies among sites. We have therefore developed new methods for estimating θ under the finite site model. Using computer simulations, it has been shown that the new methods give almost unbiased estimates even when the mutation rate varies among sites substantially. Furthermore, we have also developed new statistics for testing neutrality by modifying Tajima's D statistic. Computer simulations suggest that the new test statistics can be used even when the mutation rate varies among sites.

DOI： 10.1093/genetics/147.4.1959

Web of Science

Scopus

researchmap

Language：Japanese  

researchmap

Language：English   Book type：Scholarly book

http://www.intechopen.com/books/point-mutation

researchmap

J-GLOBAL

researchmap

Authorship：Lead author,　Last author,　Corresponding author   Language：English   Publishing type：Research paper, summary (national, other academic conference)   Publisher：GENETICS SOC JAPAN  

Web of Science

researchmap

Language：English   Publishing type：Article, review, commentary, editorial, etc. (scientific journal)   Publisher：INT PHYCOLOGICAL SOC  

Web of Science

researchmap

Language：English   Publishing type：Article, review, commentary, editorial, etc. (scientific journal)   Publisher：INT PHYCOLOGICAL SOC  

Web of Science

researchmap

Event date： 2021.8

Language：Japanese   Presentation type：Poster presentation  

researchmap

Event date： 2021.8

Language：English   Presentation type：Poster presentation  

researchmap

Event date： 2020.9

Language：Japanese   Presentation type：Symposium, workshop panel (nominated)  

researchmap

Event date： 2020.2

Language：Japanese   Presentation type：Oral presentation (general)  

researchmap

Event date： 2017.9

Language：English  

Venue：Philadelphia, PA, USA  

researchmap

Event date： 2017.5

Language：Japanese   Presentation type：Poster presentation  

researchmap

Event date： 2016.10

Language：English  

Venue：Vancouver, Canada  

researchmap

Event date： 2015.9

Language：Japanese   Presentation type：Oral presentation (general)  

Venue：仙台市  

researchmap

Event date： 2014.9

Language：Japanese   Presentation type：Oral presentation (general)  

researchmap

Event date： 2014.8

Language：Japanese   Presentation type：Oral presentation (general)  

researchmap

Event date： 2013.11

Language：Japanese   Presentation type：Oral presentation (general)  

researchmap

Event date： 2008.8

Language：Japanese   Presentation type：Poster presentation  

researchmap

Event date： 2006.8

Language：Japanese   Presentation type：Poster presentation  

researchmap

Event date： 2005.12

Language：Japanese   Presentation type：Poster presentation  

researchmap

Event date： 2004.9

Language：Japanese   Presentation type：Oral presentation (general)  

researchmap

Event date： 2000.3

Language：English   Presentation type：Poster presentation  

researchmap

Event date： 1996.10

Language：Japanese   Presentation type：Oral presentation (general)  

researchmap

Language：English  

Venue：Yokohama, Japan  

researchmap

Language：Japanese  

Venue：京都市  

researchmap

Language：English   Presentation type：Poster presentation  

researchmap

Language：Japanese   Presentation type：Poster presentation  

Venue：Orlando, Florida  

researchmap

Language：Japanese   Presentation type：Poster presentation  

researchmap

Language：English  

Venue：仙台市  

researchmap

Language：English  

Venue：Baltimore, America  

researchmap

Language：English  

Venue：Sakyo-ku, Kyoto, Japan  

researchmap

Language：English   Presentation type：Poster presentation  

Venue：Chiyoda, Tokyo, Japan  

researchmap

Language：Japanese   Presentation type：Oral presentation (general)  

researchmap

Language：English   Presentation type：Poster presentation  

researchmap

Language：Japanese   Presentation type：Poster presentation  

researchmap

Language：English   Presentation type：Poster presentation  

Venue：Houston, TX  

researchmap

Language：Japanese   Presentation type：Poster presentation  

researchmap

Language：English   Presentation type：Poster presentation  

researchmap

Language：English   Presentation type：Oral presentation (general)  

researchmap

Language：Japanese   Presentation type：Oral presentation (general)  

researchmap

Language：English   Presentation type：Oral presentation (general)  

researchmap

Language：Japanese   Presentation type：Oral presentation (general)  

researchmap

Language：Japanese   Presentation type：Oral presentation (general)  

researchmap

Language：English   Presentation type：Oral presentation (general)  

researchmap

Language：English  

Venue：Atami, Japan  

researchmap

Language：Japanese   Presentation type：Oral presentation (general)  

researchmap

Language：Japanese   Presentation type：Oral presentation (general)  

researchmap

Language：Japanese   Presentation type：Oral presentation (general)  

researchmap

Language：Japanese   Presentation type：Poster presentation  

researchmap

Language：Japanese   Presentation type：Oral presentation (general)  

researchmap

Language：Japanese   Presentation type：Poster presentation  

researchmap

Language：English   Presentation type：Poster presentation  

Venue：San Francisco  

researchmap

Language：Japanese   Presentation type：Poster presentation  

researchmap

Language：English   Presentation type：Poster presentation  

researchmap

Language：English   Presentation type：Poster presentation  

researchmap

Language：English   Presentation type：Poster presentation  

researchmap

Language：Japanese   Presentation type：Symposium, workshop panel (nominated)  

researchmap

Language：Japanese   Presentation type：Poster presentation  

researchmap

Language：Japanese   Presentation type：Oral presentation (general)  

researchmap

Language：Japanese   Presentation type：Oral presentation (general)  

Venue：Yokohama  

researchmap

Language：English   Presentation type：Poster presentation  

Venue：Philadelphia, PA  

researchmap

Language：Japanese   Presentation type：Oral presentation (general)  

researchmap

researchmap

Language：Japanese   Presentation type：Oral presentation (general)  

Venue：Niigata  

researchmap

Language：English   Presentation type：Oral presentation (general)  

researchmap

Language：English   Presentation type：Oral presentation (general)  

researchmap

Language：Japanese   Presentation type：Oral presentation (general)  

researchmap

Language：Japanese   Presentation type：Poster presentation  

Venue：Sapporo  

researchmap

Language：Japanese   Presentation type：Oral presentation (general)  

researchmap

Language：Japanese   Presentation type：Poster presentation  

Venue：盛岡市  

researchmap