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/.