Genome-wide association study identifies a new susceptibility locus in PLA2G4C for Multiple System Atrophy

To elucidate the molecular basis of multiple system atrophy (MSA), a neurodegenerative disease, we conducted a genome-wide association study (GWAS) in a Japanese MSA case/control series followed by replication studies in Japanese, Korean, Chinese, European and North American samples. In the GWAS stage rs2303744 on chromosome 19 showed a suggestive association ( P = 6.5 × 10 ^-7 ) that was replicated in additional Japanese samples ( P = 2.9 × 10 ^-6 . OR = 1.58; 95% confidence interval, 1.30 to 1.91), and then confirmed as highly significant in a meta-analysis of East Asian population data ( P = 5.0 × 10 ^-15 . Odds ratio= 1.49; 95% CI 1.35 to 1.72). The association of rs2303744 with MSA remained significant in combined European/North American samples ( P =0.023. Odds ratio=1.14; 95% CI 1.02 to 1.28) despite allele frequencies being quite different between these populations. rs2303744 leads to an amino acid substitution in PLA2G4C that encodes the cPLA2γ lysophospholipase/transacylase. The cPLA2γ-Ile143 isoform encoded by the MSA risk allele has significantly decreased transacylase activity compared with the alternate cPLA2γ-Val143 isoform that may perturb membrane phospholipids and α-synuclein biology.

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

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

Muscle transcriptomics shows overexpression of cadherin 1 in inclusion body myositis

Objective

This study aimed to elucidate the molecular features of inclusion body myositis (IBM).

Methods

We performed RNA sequencing analysis of muscle biopsy samples from 67 participants, consisting of 58 myositis patients with the pathological finding of CD8‐positive T cells invading non‐necrotic muscle fibers expressing major histocompatibility complex class I (43 IBM, 6 polymyositis, and 9 unclassifiable myositis), and 9 controls.

Results

Cluster analysis, principal component analysis, and pathway analysis showed that differentially expressed genes and pathways identified in IBM and polymyositis were mostly comparable. However, pathways related to cell adhesion molecules were upregulated in IBM as compared with polymyositis and controls (p < 0.01). Notably, CDH1, which encodes the epidermal cell junction protein cadherin 1, was overexpressed in the muscles of IBM, which was validated by another RNA sequencing dataset from previous publications. Western blotting confirmed the presence of mature cadherin 1 protein in the muscles of IBM. Immunohistochemical staining confirmed the positivity for anti‐cadherin 1 antibody in the muscles of IBM, whereas there was no muscle fiber positive for anti‐cadherin 1 antibody in immune‐mediated necrotizing myopathy, antisynthetase syndrome, and controls. The fibers stained with anti‐cadherin 1 antibody did not have rimmed vacuoles or abnormal protein accumulation. Experimental skeletal muscle regeneration and differentiation systems showed that CDH1 is expressed during skeletal muscle regeneration and differentiation.

Interpretation

CDH1 was detected as a differentially expressed gene, and immunohistochemistry showed that cadherin 1 exists in the muscles of IBM, whereas it was rarely seen in those of other idiopathic inflammatory myopathies. Cadherin 1 upregulation in muscle could provide a valuable clue to the pathological mechanisms of IBM. ANN NEUROL 2022;91:317–328

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

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

Stentless Interventional Procedure Using Rotational Atherectomy and Drug-Coated Balloon for Noncalcified De Novo Lesions

Background

Several recent reports have shown that a stentless interventional procedure using rotational atherectomy followed by drug-coated balloon (DCB) treatment (RA/DCB) is a potent revascularization therapy for calcified de novo lesions even in the new-generation drug-eluting stent era; however, the role of the RA/DCB procedure for noncalcified de novo lesions remains unclear.

Methods

A total of 47 consecutive patients (53 lesions) who underwent RA/DCB for coronary de novo lesions were enrolled. According to the presence or absence of severe calcification at target lesions on fluoroscopy, the 47 patients were divided into the noncalcified cases (n = 12) and the calcified cases (n = 35), and the 53 lesions were divided into the noncalcified lesions (n = 14) and the calcified lesions (n = 39).

Results

The noncalcified cases tended to have a higher frequency of bleeding risk and had a significantly lower prevalence of dual antiplatelet therapy compared with the calcified cases. The main lesion-specific factors for the RA/DCB procedure among the noncalcified lesions were presence of left circumflex coronary artery ostial lesion. The final burr size, DCB diameter used, and angiographic success rate did not significantly differ between the 2 groups. The noncalcified lesions had a larger reference diameter and a shorter lesion length than the calcified lesions, whereas acute gain and late lumen loss did not differ between the 2 groups. Nine-month clinical outcomes were comparable between the 2 groups.

Conclusions

Under drug-eluting stent-unsuitable clinical or lesion conditions, acute and midterm outcomes of RA/DCB for noncalcified de novo lesions might be comparable with those for calcified de novo lesions.

Clinical usefulness of multigene screening with phenotype-driven bioinformatics analysis for the diagnosis of patients with monogenic diabetes or severe insulin resistance

AIMS

Monogenic diabetes is clinically heterogeneous and differs from common forms of diabetes (type 1 and 2). We aimed to investigate the clinical usefulness of a comprehensive genetic testing system, comprised of targeted next-generation sequencing (NGS) with phenotype-driven bioinformatics analysis in patients with monogenic diabetes, which uses patient genotypic and phenotypic data to prioritize potentially causal variants.

METHODS

We performed targeted NGS of 383 genes associated with monogenic diabetes or common forms of diabetes in 13 Japanese patients with suspected (n = 10) or previously diagnosed (n = 3) monogenic diabetes or severe insulin resistance. We performed in silico structural analysis and phenotype-driven bioinformatics analysis of candidate variants from NGS data.

RESULTS

Among the patients suspected having monogenic diabetes or insulin resistance, we diagnosed 3 patients as subtypes of monogenic diabetes due to disease-associated variants of INSR, LMNA, and HNF1B. Additionally, in 3 other patients, we detected rare variants with potential phenotypic effects. Notably, we identified a novel missense variant in TBC1D4 and an MC4R variant, which together may cause a mixed phenotype of severe insulin resistance.

CONCLUSIONS

This comprehensive approach could assist in the early diagnosis of patients with monogenic diabetes and facilitate the provision of tailored therapy.

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

BACKGROUND

Although 60% of patients with de novo neurofibromatosis type 2 (NF2) are presumed to have mosaic NF2, the actual diagnostic rate of this condition remains low at around 20% because of the existing difficulties in detecting NF2 variants with low variant allele frequency (VAF). Here, we examined the correlation between the genotype and phenotype of mosaic NF2 after improving the diagnostic rate of mosaic NF2.

METHODS

We performed targeted deep sequencing of 36 genes including NF2 using DNA samples from multiple tissues (blood, buccal mucosa, hair follicle and tumour) of 53 patients with de novo NF2 and elucidated their genotype-phenotype correlation.

RESULTS

Twenty-four patients (45.2%) had the NF2 germline variant, and 20 patients with NF2 (37.7%) had mosaic NF2. The mosaic NF2 phenotype was significantly different from that in patients with NF2 germline variant in terms of distribution of NF2-related disease, tumour growth rate and hearing outcome. The behaviour of schwannoma correlated to the extent of VAF with NF2 variant in normal tissues unlike meningioma.

CONCLUSION

We have improved the diagnostic rate of mosaic NF2 compared with that of previous studies by targeted deep sequencing of DNA from multiple tissues. Many atypical patients with NF2 diagnosed with 'unilateral vestibular schwannoma' or 'multiple meningiomas' presumably have mosaic NF2. Finally, we suggest that the highly diverse phenotype of NF2 could result not only from the type and location of NF2 variant but also the extent of VAF in the NF2 variant within normal tissue DNA.

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

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

Comprehensive investigation of RNF213 nonsynonymous variants associated with intracranial artery stenosis

Intracranial artery stenosis (ICAS) is the most common cause of ischemic stroke worldwide. RNF213 single nucleotide variant c.14429G > A (p.Arg4810Lys, rs112735431) was recently reported to be associated with ICAS in East Asians. However, the disease susceptibility of other RNF213 variants has not been clarified. This study comprehensively investigated ICAS-associated RNF213 variants in a pool of 168 Japanese ICAS patients and 1,194 control subjects. We found 138 nonsynonymous germline variants by target resequencing of all coding exons in RNF213. Association study between ICAS patients and control subjects revealed that only p.Arg4810Lys had significant association with ICAS (P = 1.5 × 10^-28, odds ratio = 29.3, 95% confidence interval 15.31-56.2 [dominant model]). Fourteen of 138 variants were rare variants detected in ICAS patients not harboring p.Arg4810Lys variant. Two of these rare variants (p.Cys118Arg and p.Leu2356Phe) consistent with variants previously reported in moyamoya disease patients characterized by stenosis of intracranial artery and association with RNF213, and three rare variants (p.Ser193Gly, p.Val1817Leu, and p.Asp3329Tyr) were found neither in control subjects and Single Nucleotide Polymorphism Database. The present findings may improve our understanding of the genetic background of intracranial artery stenosis.

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

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

A Novel de novo KIF1A Mutation in a Patient with Autism, Hyperactivity, Epilepsy, Sensory Disturbance, and Spastic Paraplegia

Heterozygous mutations in KIF1A have been reported to cause syndromic intellectual disability or pure spastic paraplegia. However, their genotype-phenotype correlations have not been fully elucidated. We herein report a man with autism and hyperactivity along with sensory disturbance and spastic paraplegia, carrying a novel de novo mutation in KIF1A [c.37C>T (p.R13C)]. Autism and hyperactivity have only previously been reported in a patient with c.38 G>A (R13H) mutation. This case suggests that alterations in this arginine at codon 13 might lead to a common clinical spectrum and further expand the genetic and clinical spectra associated with KIF1A mutations.

Long-read metagenomic exploration of extrachromosomal mobile genetic elements in the human gut

BACKGROUND

Elucidating the ecological and biological identity of extrachromosomal mobile genetic elements (eMGEs), such as plasmids and bacteriophages, in the human gut remains challenging due to their high complexity and diversity.

RESULTS

Here, we show efficient identification of eMGEs as complete circular or linear contigs from PacBio long-read metagenomic data. De novo assembly of PacBio long reads from 12 faecal samples generated 82 eMGE contigs (2.5~666.7-kb), which were classified as 71 plasmids and 11 bacteriophages, including 58 novel plasmids and six bacteriophages, and complete genomes of five diverse crAssphages with terminal direct repeats. In a dataset of 413 gut metagenomes from five countries, many of the identified plasmids were highly abundant and prevalent. The ratio of gut plasmids by our plasmid data is more than twice that in the public database. Plasmids outnumbered bacterial chromosomes three to one on average in this metagenomic dataset. Host prediction suggested that Bacteroidetes-associated plasmids predominated, regardless of microbial abundance. The analysis found several plasmid-enriched functions, such as inorganic ion transport, while antibiotic resistance genes were harboured mostly in low-abundance Proteobacteria-associated plasmids.

CONCLUSIONS

Overall, long-read metagenomics provided an efficient approach for unravelling the complete structure of human gut eMGEs, particularly plasmids.

Mutations in COA7 cause spinocerebellar ataxia with axonal neuropathy

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

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

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

Recompleting the Caenorhabditis elegans genome

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

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

OBJECTIVES

To evaluate the burden of rare variants in the causative genes for amyotrophic lateral sclerosis (ALS) on the age at onset of ALS in a Japanese case series.

METHODS

We conducted whole-exome sequencing analysis of 89 families with familial ALS (FALS) and 410 patients with sporadic ALS (SALS) to identify known pathogenic mutations or rare functionally predicted deleterious variants in the causative genes for ALS. Rare variants (minor allele frequency <1%) with scaled Combined Annotation-Dependent Depletion score >20 were defined as rare functionally predicted deleterious variants. The patients with ALS were classified on the basis of the number of pathogenic and/or rare functionally predicted deleterious variants, and the age at onset was compared among the classified groups.

RESULTS

Whole-exome sequencing analysis revealed known pathogenic mutations or rare functionally predicted deleterious variants in causative genes for ALS in 56 families with FALS (62.9%) and 87 patients with SALS (21.2%). Such variants in multiple genes were identified in seven probands with FALS and eight patients with SALS. The ages at onset in the patients with ALS with multiple variants were significantly earlier than those in other patients with ALS. Even when the patients with known pathogenic mutations were excluded, a significantly earlier onset of the disease was still observed in patients with multiple rare functionally predicted deleterious variants.

CONCLUSIONS

A substantial number of patients carried rare variants in multiple genes, and the burden of rare variants in the known causative genes for ALS affects the age at onset in the Japanese ALS series.

Neuroimaging, genetic, and enzymatic study in a Japanese family with a GBA gross deletion

INTRODUCTION

Glucocerebrosidase gene (GBA) variants are associated with Parkinson's disease (PD) and dementia with Lewy bodies (DLB). The molecular mechanisms underlying these diseases with GBA variants, however, are not well understood. In order to determine the effect of a deletion mutation in GBA, we performed a neuroimaging, genetic, and enzymatic study in a Japanese family with a gross deletion of exons 3 to 11 in GBA.

METHODS

We performed [¹²³I] FP-CIT SPECT and [¹²³I] N-isopropyl-p-iodoamphetamine SPECT (IMP-SPECT), and determined GBA expression and glucocerebrosidase (GCase) activity in leukocytes in two GBA-associated PD patients and nine unaffected individuals (including four mutation carriers) in a Japanese family with a heterozygous gross deletion mutation in the GBA gene.

RESULTS

The two PD patients and two of the four clinically unaffected carriers showed decreased [¹²³I] FP-CIT uptake. IMP-SPECT showed a pattern like that in DLB in one patient. When we compared PD patients with GBA mutations with clinically unaffected carriers, there was a poor correlation between the development of PD and the expression level of GBA or GCase activity.

CONCLUSION

We confirmed the gross deletion mutation in the GBA gene, which appeared to be associated with the PD or reduced [¹²³I] FP-CIT in this family. However, since we cannot conclude whether a reduction of GCase activity is directly correlated with the pathogenesis of PD or not, longitudinal follow-up of this family is needed.

Centromere evolution and CpG methylation during vertebrate speciation

Centromeres and large-scale structural variants evolve and contribute to genome diversity during vertebrate speciation. Here, we perform de novo long-read genome assembly of three inbred medaka strains that are derived from geographically isolated subpopulations and undergo speciation. Using single-molecule real-time (SMRT) sequencing, we obtain three chromosome-mapped genomes of length ~734, ~678, and ~744Mbp with a resource of twenty-two centromeric regions of length 20-345kbp. Centromeres are positionally conserved among the three strains and even between four pairs of chromosomes that were duplicated by the teleost-specific whole-genome duplication 320-350 million years ago. The centromeres do not all evolve at a similar pace; rather, centromeric monomers in non-acrocentric chromosomes evolve significantly faster than those in acrocentric chromosomes. Using methylation sensitive SMRT reads, we uncover centromeres are mostly hypermethylated but have hypomethylated sub-regions that acquire unique sequence compositions independently. These findings reveal the potential of non-acrocentric centromere evolution to contribute to speciation.

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

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

Molecular epidemiological study of familial amyotrophic lateral sclerosis in Japanese population by whole-exome sequencing and identification of novel HNRNPA1 mutation

To elucidate the genetic epidemiology of familial amyotrophic lateral sclerosis (FALS) in the Japanese population, we conducted whole-exome sequencing analysis of 30 FALS families in whom causative mutations have not been identified in previous studies. Consequently, whole-exome sequencing analysis revealed novel mutations in HNRNPA1, TBK1, and VCP. Taken together with our previous results of mutational analyses by direct nucleotide sequencing analysis, a microarray-based resequencing method, or repeat-primed PCR analysis, causative mutations were identified in 41 of the 68 families (60.3%) with SOD1 being the most frequent cause of FALS (39.7%). Of the mutations identified in this study, a novel c.862/1018C>G (p.P288A/340A) mutation in HNRNPA1 located in the nuclear localization signal domain of hnRNPA1, enhances the recruitment of mutant hnRNPA1 into stress granules, indicating that an altered nuclear localization signal activity plays an essential role in amyotrophic lateral sclerosis pathogenesis.

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

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

Cancer can result from genetic alterations, some of which can be good drug targets. To reveal genetic alterations that provide important information for the development of ideal therapeutic strategies for triple-negative breast cancer (TNBC), TNBC tumor samples were subjected to comprehensive genomic analyses. We identified novel recurrent structural variations associated with enhanced expression of the TGFA gene that encodes one of the high affinity ligands for epidermal growth factor receptor (EGFR). Although TGFA expression is known to be elevated in a subset of TNBC tumors, this is the first report of the mechanistic basis of this phenomenon. It is of particular importance considering that anti-EGFR agents are possible therapeutic options for TNBC patients. Our study also revealed several features associated with “BRCAness”, which is critical for identification of patients who may be responsive to platinum agents and/or poly ADP-ribose polymerase inhibitors. Thus, the data presented in this report may advance our understanding of the pathogenesis of TNBC.

Ataxic form of autosomal recessive PEX10-related peroxisome biogenesis disorders with a novel compound heterozygous gene mutation and characteristic clinical phenotype

Peroxisome biogenesis factor 10 (PEX10) is involved in the import of peroxisomal matrix proteins, and the mutation of this gene causes 3 subtypes of peroxisome biogenesis disorders, namely Zellweger syndrome (severe), neonatal adrenoleukodystrophy (moderate) and an ataxic form (mild). Here, we report 3 siblings of the ataxic form with cerebellar ataxia, mild mental retardation, and 3 additional characteristic features: mydriasis, hyperreflexia and involuntary head movement. All 3 siblings are compound heterozygous for a previously reported mutation, c.2T>C (p.M1T), and a novel mutation, c.920G>A, causing a missense change (p.C307Y) located in the RING finger domain of PEX10. The present cases suggest that these PEX10 mutations involve not only cerebellar but also more multiple nervous systems including pupillary autonomic, pyramidal and extrapyramidal systems.

Modeling neurological diseases with induced pluripotent cells reprogrammed from immortalized lymphoblastoid cell lines

Patient-specific induced pluripotent stem cells (iPSCs) facilitate understanding of the etiology of diseases, discovery of new drugs and development of novel therapeutic interventions. A frequently used starting source of cells for generating iPSCs has been dermal fibroblasts (DFs) isolated from skin biopsies. However, there are also numerous repositories containing lymphoblastoid B-cell lines (LCLs) generated from a variety of patients. To date, this rich bioresource of LCLs has been underused for generating iPSCs, and its use would greatly expand the range of targeted diseases that could be studied by using patient-specific iPSCs. However, it remains unclear whether patient’s LCL-derived iPSCs (LiPSCs) can function as a disease model. Therefore, we generated Parkinson’s disease patient-specific LiPSCs and evaluated their utility as tools for modeling neurological diseases. We established iPSCs from two LCL clones, which were derived from a healthy donor and a patient carrying PARK2 mutations, by using existing non-integrating episomal protocols. Whole genome sequencing (WGS) and comparative genomic hybridization (CGH) analyses showed that the appearance of somatic variations in the genomes of the iPSCs did not vary substantially according to the original cell types (LCLs, T-cells and fibroblasts). Furthermore, LiPSCs could be differentiated into functional neurons by using the direct neurosphere conversion method (dNS method), and they showed several Parkinson’s disease phenotypes that were similar to those of DF-iPSCs. These data indicate that the global LCL repositories can be used as a resource for generating iPSCs and disease models. Thus, LCLs are the powerful tools for generating iPSCs and modeling neurological diseases.

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

Motivation: Determining the methylation state of regions with high copy numbers is challenging for second-generation sequencing, because the read length is insufficient to map reads uniquely, especially when repetitive regions are long and nearly identical to each other. Single-molecule real-time (SMRT) sequencing is a promising method for observing such regions, because it is not vulnerable to GC bias, it produces long read lengths, and its kinetic information is sensitive to DNA modifications.

Results: We propose a novel linear-time algorithm that combines the kinetic information for neighboring CpG sites and increases the confidence in identifying the methylation states of those sites. Using a practical read coverage of ∼30-fold from an inbred strain medaka (Oryzias latipes), we observed that both the sensitivity and precision of our method on individual CpG sites were ∼93.7%. We also observed a high correlation coefficient (R = 0.884) between our method and bisulfite sequencing, and for 92.0% of CpG sites, methylation levels ranging over [0,1] were in concordance within an acceptable difference 0.25. Using this method, we characterized the landscape of the methylation status of repetitive elements, such as LINEs, in the human genome, thereby revealing the strong correlation between CpG density and hypomethylation and detecting hypomethylation hot spots of LTRs and LINEs. We uncovered the methylation states for nearly identical active transposons, two novel LINE insertions of identity ∼99% and length 6050 base pairs (bp) in the human genome, and 16 Tol2 elements of identity >99.8% and length 4682 bp in the medaka genome.

Availability and Implementation: AgIn (Aggregate on Intervals) is available at: https://github.com/hacone/AgIn

Contact:ysuzuki@cb.k.u-tokyo.ac.jp or moris@cb.k.u-tokyo.ac.jp

Supplementary information:Supplementary data are available at Bioinformatics online.

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

Objective

The objective of this study was to identify new causes of Charcot–Marie–Tooth (CMT) disease in patients with autosomal‐recessive (AR) CMT.

Methods

To efficiently identify novel causative genes for AR‐CMT, we analyzed 303 unrelated Japanese patients with CMT using whole‐exome sequencing and extracted recessive variants/genes shared among multiple patients. We performed mutation screening of the newly identified membrane metalloendopeptidase (MME) gene in 354 additional patients with CMT. We clinically, genetically, pathologically, and radiologically examined 10 patients with the MME mutation.

Results

We identified recessive mutations in MME in 10 patients. The MME gene encodes neprilysin (NEP), which is well known to be one of the most prominent beta‐amyloid (Aβ)‐degrading enzymes. All patients had a similar phenotype consistent with late‐onset axonal neuropathy. They showed muscle weakness, atrophy, and sensory disturbance in the lower extremities. All the MME mutations could be loss‐of‐function mutations, and we confirmed a lack/decrease of NEP protein expression in a peripheral nerve. No patients showed symptoms of dementia, and 1 patient showed no excess Aβ in Pittsburgh compound‐B positron emission tomography imaging.

Interpretation

Our results indicate that loss‐of‐function MME mutations are the most frequent cause of adult‐onset AR‐CMT2 in Japan, and we propose that this new disease should be termed AR‐CMT2T. A loss‐of‐function MME mutation did not cause early‐onset Alzheimer's disease. Identifying the MME mutation responsible for AR‐CMT could improve the rate of molecular diagnosis and the understanding of the molecular mechanisms of CMT. Ann Neurol 2016;79:659–672

Atypical parkinsonism caused by Pro105Leu mutation of prion protein: A broad clinical spectrum

Objective:

To delineate molecular and clinical characteristics of 3 families with PRNP P105L mutation, a variant of Gerstmann-Sträussler-Scheinker syndrome whose main motor symptoms were parkinsonism and/or involuntary movements.

Methods:

The causative mutation was first determined in the affected patients of family 1 using whole-exome sequencing, and then mutational analysis was extended to families 2 and 3. The clinical features of the patients of these 3 families were summarized. Haplotype analysis was performed using high-density single nucleotide polymorphism array.

Results:

The whole-exome sequencing revealed that the heterozygous mutation c.314C>T (p.P105L) in PRNP was the only known pathogenic mutation shared by the 3 patients of the family with autosomal dominant parkinsonism. We further identified the same mutation in patients of the other 2 families with autosomal dominant parkinsonism and/or involuntary movements. The clinical features of our patients with PRNP P105L mutation included various motor symptoms such as parkinsonism and involuntary movements in addition to progressive dementia. The clinical features in part overlapped with those of other forms of inherited prion diseases, such as fatal familial insomnia and Huntington disease-like type 1. The patients with PRNP P105L mutation shared a haplotype spanning 7.1 Mb around PRNP, raising the possibility that the mutations in the patients originated from a common founder.

Conclusion:

Most of the patients presented with parkinsonism in addition to progressive dementia. Although spastic paraparesis has been emphasized as the main clinical feature, the clinical spectrum of patients with PRNP P105L is broader than expected.

Associations between nucleosome phasing, sequence asymmetry, and tissue-specific expression in a set of inbred Medaka species

BACKGROUND

Transcription start sites (TSSs) with pronounced and phased nucleosome arrays downstream and nucleosome-depleted regions upstream of TSSs are observed in various species.

RESULTS

We have characterized sequence variation and expression properties of this set of TSSs (which we call "Nucleocyclic TSSs") using germline and somatic cells of three medaka (Oryzias latipes) inbred isolates from different locations. We found nucleocyclic TSSs in medaka to be associated with higher gene expression and characterized by a clear boundary in sequence composition with potentially-nucleosome-destabilizing A/T-enrichment upstream (p < 10(-60)) and nucleosome- accommodating C/G-enrichment downstream (p < 10(-40)) that was highly conserved from an ancestor. A substantial genetic distance between the strains facilitated the in-depth analysis of patterns of fixed mutations, revealing a localization-specific equilibrium between the rates of distinct mutation categories that would serve to maintain the conserved sequence anisotropy around TSSs. Downstream of nucleocyclic TSSs, C to T, T to C, and other mutation rates on the sense strand increased around first nucleosome dyads and decreased around first linkers, which contrasted with genomewide mutational patterns around nucleosomes (p < 5 %). C to T rates are higher than G to A rates around nucleosome associated with germline nucleocyclic TSS sites (p < 5 %), potentially due to the asymmetric effect of transcription-coupled repair.

CONCLUSIONS

Our results demonstrate an atypical evolutionary process surrounding nucleocyclic TSSs.

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

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

A recurrent de novo FAM111A mutation causes Kenny-Caffey syndrome type 2

Kenny-Caffey syndrome (KCS) is a rare dysmorphologic syndrome characterized by proportionate short stature, cortical thickening and medullary stenosis of tubular bones, delayed closure of anterior fontanelle, eye abnormalities, and hypoparathyroidism. The autosomal dominant form of KCS (KCS type 2 [KCS2]) is distinguished from the autosomal recessive form of KCS (KCS type 1 [KCS1]), which is caused by mutations of the tubulin-folding cofactor E (TBCE) gene, by the absence of mental retardation. In this study, we recruited four unrelated Japanese patients with typical sporadic KCS2, and performed exome sequencing in three patients and their parents to elucidate the molecular basis of KCS2. The possible candidate genes were explored by a de novo mutation detection method. A single gene, FAM111A (NM_001142519.1), was shared among three families. An identical missense mutation, R569H, was heterozygously detected in all three patients but not in the unaffected family members. This mutation was also found in an additional unrelated patient. These findings are in accordance with those of a recent independent report by a Swiss group that KCS2 is caused by a de novo mutation of FAM111A, and R569H is a hot spot mutation for KCS2. Although the function of FAM111A is not known, this study would provide evidence that FAM111A is a key molecule for normal bone development, height gain, and parathyroid hormone development and/or regulation.

Rapid detection of expanded short tandem repeats in personal genomics using hybrid sequencing

Motivation: Long expansions of short tandem repeats (STRs), i.e. DNA repeats of 2–6 nt, are associated with some genetic diseases. Cost-efficient high-throughput sequencing can quickly produce billions of short reads that would be useful for uncovering disease-associated STRs. However, enumerating STRs in short reads remains largely unexplored because of the difficulty in elucidating STRs much longer than 100 bp, the typical length of short reads.

Results: We propose ab initio procedures for sensing and locating long STRs promptly by using the frequency distribution of all STRs and paired-end read information. We validated the reproducibility of this method using biological replicates and used it to locate an STR associated with a brain disease (SCA31). Subsequently, we sequenced this STR site in 11 SCA31 samples using SMRT^TM sequencing (Pacific Biosciences), determined 2.3–3.1 kb sequences at nucleotide resolution and revealed that (TGGAA)- and (TAAAATAGAA)-repeat expansions determined the instability of the repeat expansions associated with SCA31. Our method could also identify common STRs, (AAAG)- and (AAAAG)-repeat expansions, which are remarkably expanded at four positions in an SCA31 sample. This is the first proposed method for rapidly finding disease-associated long STRs in personal genomes using hybrid sequencing of short and long reads.

Availability and implementation: Our TRhist software is available at http://trhist.gi.k.u-tokyo.ac.jp/.

Contact:moris@cb.k.u-tokyo.ac.jp

Supplementary information: Supplementary data are available at Bioinformatics online.

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

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

Exome analysis reveals a Japanese family with spinocerebellar ataxia, autosomal recessive 1

Spinocerebellar ataxia autosomal recessive 1 (SCAR1/AOA2) is clinically characterized by an early-onset progressive cerebellar ataxia with axonal neuropathy, ocular motor apraxia, and elevation of serum alpha-fetoprotein level. The disorder is caused by mutations in senataxin (SETX) gene. Here, we report a Japanese SCAR1/AOA2 family with a homozygous nonsense mutation (p.Q1441X) of SETX that was identified by exome sequencing. The family was previously reported as early-onset ataxia of undetermined cause. The present study emphasized the role of whole exome-sequence analysis to establish the molecular diagnosis of neurodegenerative disease presenting with diverse clinical presentations.

ATF6α/β-mediated adjustment of ER chaperone levels is essential for development of the notochord in medaka fish

ATF6α and ATF6β are membrane-bound transcription factors activated by regulated intramembrane proteolysis in response to endoplasmic reticulum (ER) stress to induce various ER quality control proteins. ATF6α- and ATF6β single-knockout mice develop normally, but ATF6α/β double knockout causes embryonic lethality, the reason for which is unknown. Here we show in medaka fish that ATF6α is primarily responsible for transcriptional induction of the major ER chaperone BiP and that ATF6α/β double knockout, but not ATF6α- or ATF6β single knockout, causes embryonic lethality, as in mice. Analyses of ER stress reporters reveal that ER stress occurs physiologically during medaka early embryonic development, particularly in the brain, otic vesicle, and notochord, resulting in ATF6α- and ATF6β-mediated induction of BiP, and that knockdown of the α1 chain of type VIII collagen reduces such ER stress. The absence of transcriptional induction of several ER chaperones in ATF6α/β double knockout causes more profound ER stress and impaired notochord development, which is partially rescued by overexpression of BiP. Thus ATF6α/β-mediated adjustment of chaperone levels to increased demands in the ER is essential for development of the notochord, which synthesizes and secretes large amounts of extracellular matrix proteins to serve as the body axis before formation of the vertebra.

CSF1R mutations identified in three families with autosomal dominantly inherited leukoencephalopathy

Genetic and phenotypic heterogeneities are considerably high in adult-onset leukoencephalopathy, in which comprehensive mutational analyses of the candidate genes by conventional methods are too laborious. We applied exome sequencing to conduct a comprehensive mutational analysis of genes for autosomal dominant leukoencephalopathies. Genomic DNA samples from four patients of three families with autosomal dominantly inherited adult-onset leukodystrophy were subjected to exome sequencing. On the basis of the results, 21 patients with adult-onset sporadic leukodystrophy and one patient with pathologically proven HDLS were additionally screened for CSF1R mutations. Exome sequencing identified heterozygous CSF1R mutations (p.I794T and p.R777W) in two families. I794T has recently been reported as a causative mutation for hereditary diffuse leukoencephalopathy with spheroids (HDLS), and R777W is a novel mutation. Although mutational analysis of CSF1R in 21 sporadic cases revealed no mutations, another novel CSF1R mutation, p.C653Y, was identified in one patient with autopsy-proven HDSL. These variants were located in the PTK domain where the causative mutations cluster. Functional prediction of the mutant CSF1R as well as cross-species conservation of the affected amino acids supports the notion that these variants are pathogenic for HDLS. Exome sequencing is useful for a comprehensive mutational analysis of causative genes for hereditary leukoencephalopathies, and CSF1R should be considered a candidate gene for patients with autosomal dominant leukoencephalopathies.

siDirect 2.0: updated software for designing functional siRNA with reduced seed-dependent off-target effect

BACKGROUND

RNA interference (RNAi), mediated by 21-nucleotide (nt)-length small interfering RNAs (siRNAs), is a powerful tool not only for studying gene function but also for therapeutic applications. RNAi, requiring perfect complementarity between the siRNA guide strand and the target mRNA, was believed to be extremely specific. However, a recent growing body of evidence has suggested that siRNA could down-regulate unintended genes whose transcripts possess complementarity to the 7-nt siRNA seed region. This off-target gene silencing may often provide incongruous results obtained from knockdown experiments, leading to misinterpretation. Thus, an efficient algorithm for designing functional siRNAs with minimal off-target effect based on the mechanistic features is considered of value.

RESULTS

We present siDirect 2.0, an update of our web-based software siDirect, which provides functional and off-target minimized siRNA design for mammalian RNAi. The previous version of our software designed functional siRNAs by considering the relationship between siRNA sequence and RNAi activity, and provided them along with the enumeration of potential off-target gene candidates by using a fast and sensitive homology search algorithm. In the new version, the siRNA design algorithm is extensively updated to eliminate off-target effects by reflecting our recent finding that the capability of siRNA to induce off-target effect is highly correlated to the thermodynamic stability, or the melting temperature (Tm), of the seed-target duplex, which is formed between the nucleotides positioned at 2-8 from the 5' end of the siRNA guide strand and its target mRNA. Selection of siRNAs with lower seed-target duplex stabilities (benchmark Tm < 21.5 degrees C) followed by the elimination of unrelated transcripts with nearly perfect match should minimize the off-target effects.

CONCLUSION

siDirect 2.0 provides functional, target-specific siRNA design with the updated algorithm which significantly reduces off-target silencing. When the candidate functional siRNAs could form seed-target duplexes with Tm values below 21.5 degrees C, and their 19-nt regions spanning positions 2-20 of both strands have at least two mismatches to any other non-targeted transcripts, siDirect 2.0 can design at least one qualified siRNA for >94% of human mRNA sequences in RefSeq. siDirect 2.0 is available at http://siDirect2.RNAi.jp/.

UTGB toolkit for personalized genome browsers

The advent of high-throughput DNA sequencers has increased the pace of collecting enormous amounts of genomic information, yielding billions of nucleotides on a weekly basis. This advance represents an improvement of two orders of magnitude over traditional Sanger sequencers in terms of the number of nucleotides per unit time, allowing even small groups of researchers to obtain huge volumes of genomic data over fairly short period. Consequently, a pressing need exists for the development of personalized genome browsers for analyzing these immense amounts of locally stored data. The UTGB (University of Tokyo Genome Browser) Toolkit is designed to meet three major requirements for personalization of genome browsers: easy installation of the system with minimum efforts, browsing locally stored data and rapid interactive design of web interfaces tailored to individual needs. The UTGB Toolkit is licensed under an open source license.

AVAILABILITY

The software is freely available at http://utgenome.org/.