Zhong Sheng Sun

mirTools: microRNA profiling and discovery based on high-throughput sequencing

miRNAs are small, non-coding RNA that negatively regulate gene expression at post-transcriptional level, which play crucial roles in various physiological and pathological processes, such as development and tumorigenesis. Although deep sequencing technologies have been applied to investigate various small RNA transcriptomes, their computational methods are far away from maturation as compared to microarray-based approaches. In this study, a comprehensive web server mirTools was developed to allow researchers to comprehensively characterize small RNA transcriptome. With the aid of mirTools, users can: (i) filter low-quality reads and 3/5′ adapters from raw sequenced data; (ii) align large-scale short reads to the reference genome and explore their length distribution; (iii) classify small RNA candidates into known categories, such as known miRNAs, non-coding RNA, genomic repeats and coding sequences; (iv) provide detailed annotation information for known miRNAs, such as miRNA/miRNA*, absolute/relative reads count and the most abundant tag; (v) predict novel miRNAs that have not been characterized before; and (vi) identify differentially expressed miRNAs between samples based on two different counting strategies: total read tag counts and the most abundant tag counts. We believe that the integration of multiple computational approaches in mirTools will greatly facilitate current microRNA researches in multiple ways. mirTools can be accessed at http://centre.bioinformatics.zj.cn/mirtools/ and http://59.79.168.90/mirtools.

Genes with de novo mutations are shared by four neuropsychiatric disorders discovered from NPdenovo database

Currently, many studies on neuropsychiatric disorders have utilized massive trio-based whole-exome sequencing (WES) and whole-genome sequencing (WGS) to identify numerous de novo mutations (DNMs). Here, we retrieved 17 104 DNMs from 3555 trios across four neuropsychiatric disorders: autism spectrum disorder, epileptic encephalopathy, intellectual disability and schizophrenia, in addition to unaffected siblings (control), from 36 studies by WES/WGS. After eliminating non-exonic variants, we focused on 3334 exonic DNMs for evaluation of their association with these diseases. Our results revealed a higher prevalence of DNMs in the probands of all four disorders compared with the one in the controls (P<1.3 × 10−7). The elevated DNM frequency is dominated by loss-of-function/deleterious single-nucleotide variants and frameshift indels (that is, extreme mutations, P<4.5 × 10−5). With extensive annotation of these ‘extreme’ mutations, we prioritized 764 candidate genes in these four disorders. A combined analysis of Gene Ontology, microRNA targets and transcription factor targets revealed shared biological process and non-coding regulatory elements of candidate genes in the pathology of neuropsychiatric disorders. In addition, weighted gene co-expression network analysis of human laminar-specific neocortical expression data showed that candidate genes are convergent on eight shared modules with specific layer enrichment and biological process features. Furthermore, we identified that 53 candidate genes are associated with more than one disorder (P<0.000001), suggesting a possibly shared genetic etiology underlying these disorders. Particularly, DNMs of the SCN2A gene are frequently occurred across all four disorders. Finally, we constructed a freely available NPdenovo database, which provides a comprehensive catalog of the DNMs identified in neuropsychiatric disorders.

Characterization of the Small RNA Transcriptomes of Androgen Dependent and Independent Prostate Cancer Cell Line by Deep Sequencing

Given the important roles of miRNA in post-transcriptional regulation and its implications for cancer, characterization of miRNA facilitates us to uncover molecular mechanisms underlying the progression of androgen-independent prostate cancer (PCa). The emergence of next-generation sequencing technologies has dramatically changed the speed of all aspects of sequencing in a rapid and cost-effective fashion, which can permit an unbiased, quantitive and in-depth investigation of small RNA transcriptome. In this study, we used high-throughput Illumina sequencing to comprehensively represent the full complement of individual small RNA and to characterize miRNA expression profiles in both the androgen dependent and independent Pca cell line. At least 83 miRNAs are significantly differentially expressed, of which 41 are up-regulated and 42 are down-regulated, indicating these miRNAs may be involved in the transition of LNCaP to an androgen-independent phenotype. In addition, we have identified 43 novel miRNAs from the androgen dependent and independent PCa library and 3 of them are specific to the androgen-independent PCa. Function annotation of target genes indicated that most of these differentially expressed miRNAs tend to target genes involved in signal transduction and cell communication, epically the MAPK signaling pathway. The small RNA transcriptomes obtained in this study provide considerable insights into a better understanding of the expression and function of small RNAs in the development of androgen-independent prostate cancer.

Mutations of ANK3 identified by exome sequencing are associated with Autism susceptibility

Autism spectrum disorders (ASDs) are common neurodevelopmental disorders with a strong genetic etiology. However, due to the extreme genetic heterogeneity of ASDs, traditional approaches for gene discovery are challenging. Next‐generation sequencing technologies offer an opportunity to accelerate the identification of the genetic causes of ASDs. Here, we report the results of whole‐exome sequence in a cohort of 20 ASD patients. By extensive bioinformatic analysis, we identified novel mutations in seven genes that are implicated in synaptic function and neurodevelopment. After sequencing an additional 47 ASD samples, we identified three different missense mutations in ANK3 in four unrelated ASD patients, one of which, c.4705T>G (p.S1569A), is a de novo mutation. Given the fact that ANK3 has been shown to strongly associate with schizophrenia and bipolar disorder, our findings support an association between ANK3 mutations and ASD susceptibility and imply a shared molecular pathophysiology between ASDs and other neuropsychiatric disorders. Hum Mutat 33:1635–1638, 2012.

Reversal of Cocaine-Conditioned Place Preference through Methyl Supplementation in Mice: Altering Global DNA Methylation in the Prefrontal Cortex

Analysis of global methylation in cells has revealed correlations between overall DNA methylation status and some biological states. Recent studies suggest that epigenetic regulation through DNA methylation could be responsible for neuroadaptations induced by addictive drugs. However, there is no investigation to determine global DNA methylation status following repeated exposure to addictive drugs. Using mice conditioned place preference (CPP) procedure, we measured global DNA methylation level in the nucleus accumbens (NAc) and the prefrontal cortex (PFC) associated with drug rewarding effects. We found that cocaine-, but not morphine- or food-CPP training decreased global DNA methylation in the PFC. Chronic treatment with methionine, a methyl donor, for 25 consecutive days prior to and during CPP training inhibited the establishment of cocaine, but not morphine or food CPP. We also found that both mRNA and protein level of DNMT (DNA methytransferase) 3b in the PFC were downregulated following the establishment of cocaine CPP, and the downregulation could be reversed by repeated administration of methionine. Our study indicates a crucial role of global PFC DNA hypomethylation in the rewarding effects of cocaine. Reversal of global DNA hypomethylation could significantly attenuate the rewarding effects induced by cocaine. Our results suggest that methionine may have become a potential therapeutic target to treat cocaine addiction.

mirTools 2.0 for non-coding RNA discovery, profiling, and functional annotation based on high-throughput sequencing

Next-generation sequencing has been widely applied to understand the complexity of non-coding RNAs (ncRNAs) in a cost-effective way. In this study, we developed mirTools 2.0, an updated version of mirTools 1.0, which includes the following new features. (1) From miRNA discovery in mirTools 1.0, mirTools 2.0 allows users to detect and profile various types of ncRNAs, such as miRNA, tRNA, snRNA, snoRNA, rRNA, and piRNA. (2) From miRNA profiling in mirTools 1.0, mirTools 2.0 allows users to identify miRNA-targeted genes and performs detailed functional annotation of miRNA targets, including Gene Ontology, KEGG pathway and protein-protein interaction. (3) From comparison of two samples for differentially expressed miRNAs in mirTools 1.0, mirTools 2.0 allows users to detect differentially expressed ncRNAs between two experimental groups or among multiple samples. (4) Other significant improvements include strategies used to detect novel miRNAs and piRNAs, more taxonomy categories to discover more known miRNAs and a stand-alone version of mirTools 2.0. In conclusion, we believe that mirTools 2.0 (122.228.158.106/mr2_dev and centre.bioinformatics.zj.cn/mr2_dev) will provide researchers with more detailed insight into small RNA transcriptomes.

MagicViewer: integrated solution for next-generation sequencing data visualization and genetic variation detection and annotation

New sequencing technologies, such as Roche 454, ABI SOLiD and Illumina, have been increasingly developed at an astounding pace with the advantages of high throughput, reduced time and cost. To satisfy the impending need for deciphering the large-scale data generated from next-generation sequencing, an integrated software MagicViewer is developed to easily visualize short read mapping, identify and annotate genetic variation based on the reference genome. MagicViewer provides a user-friendly environment in which large-scale short reads can be displayed in a zoomable interface under user-defined color scheme through an operating system-independent manner. Meanwhile, it also holds a versatile computational pipeline for genetic variation detection, filtration, annotation and visualization, providing details of search option, functional classification, subset selection, sequence association and primer design. In conclusion, MagicViewer is a sophisticated assembly visualization and genetic variation annotation tool for next-generation sequencing data, which can be widely used in a variety of sequencing-based researches, including genome re-sequencing and transcriptome studies. MagicViewer is freely available at http://bioinformatics.zj.cn/magicviewer/.

Recent Progress in CRISPR/Cas9 Technology.

The clustered regularly interspaced short palindromic repeats (CRISPR)/Cas9 system, a simple and efficient tool for genome editing, has experienced rapid progress in its technology and applicability in the past two years. Here, we review the recent advances in CRISPR/Cas9 technology and the ways that have been adopted to expand our capacity for precise genome manipulation. First, we introduce the mechanism of CRISPR/Cas9, including its biochemical and structural implications. Second, we highlight the latest improvements in the CRISPR/Cas9 system, especially Cas9 protein modifications for customization. Third, we review its current applications, in which the versatile CRISPR/Cas9 system was employed to edit the genome, epigenome, or RNA of various organisms. Although CRISPR/Cas9 allows convenient genome editing accompanied by many benefits, we should not ignore the significant ethical and biosafety concerns that it raises. Finally, we discuss the prospective applications and challenges of several promising techniques adapted from CRISPR/Cas9.

EpilepsyGene: a genetic resource for genes and mutations related to epilepsy

Epilepsy is one of the most prevalent chronic neurological disorders, afflicting about 3.5–6.5 per 1000 children and 10.8 per 1000 elderly people. With intensive effort made during the last two decades, numerous genes and mutations have been published to be associated with the disease. An organized resource integrating and annotating the ever-increasing genetic data will be imperative to acquire a global view of the cutting-edge in epilepsy research. Herein, we developed EpilepsyGene (http://61.152.91.49/EpilepsyGene). It contains cumulative to date 499 genes and 3931 variants associated with 331 clinical phenotypes collected from 818 publications. Furthermore, in-depth data mining was performed to gain insights into the understanding of the data, including functional annotation, gene prioritization, functional analysis of prioritized genes and overlap analysis focusing on the comorbidity. An intuitive web interface to search and browse the diversified genetic data was also developed to facilitate access to the data of interest. In general, EpilepsyGene is designed to be a central genetic database to provide the research community substantial convenience to uncover the genetic basis of epilepsy.

Antinociceptive synergistic effect of spinal mGluR2/3 antagonist and glial cells inhibitor on peripheral inflammation-induced mechanical hypersensitivity

Metabotropic glutamate receptor (mGluR) 2/3 is distributed in neurons and glial cells in many regions of the nervous system, but its role in nociceptive processing is unclear. In this study, we examined the mRNA expressions of mGluR2 and mGluR3, by real-time RT-PCR, in the spinal cord. We further investigated the possible involvement of mGluR2/3 and mechanisms underlying peripheral inflammatory pain induced by subcutaneous complete Freunds adjuvant (CFA) injection. We demonstrate that compared to the controls, the mRNA expression levels of mGluR2 and mGluR3 were significantly higher 4h after CFA injection. Functionally, blocking mGluR2/3 by their antagonist (2S)-2-amino-2-[(1S, 2S)-2-carboxycycloprop-1-yl]-3-(xanth-9-yl) propanoic acid (LY341495) alleviated the CFA-induced mechanical allodynia and the inhibitory effects were reversed by mGluR2/3 agonist (2R, 4R)-4-aminopyrrolidine-2,4-dicarboxylate ((2R, 4R)-APDC). In addition, a glial metabolism inhibitor dl-fluorocitric acid barium salt (fluorocitric acid) also inhibited the CFA-induced mechanical allodynia in a dose-dependent manner. Remarkably, simultaneous inhibition of mGluR2/3 and glial metabolism had synergistic effects. The co-administration of LY341495 and fluorocitric acid with minimal dosages produced significant more inhibition than the additive effects by the individual inhibitor alone. In summary, our data suggest that spinal mGluR2/3 contributes to the generation of mechanical allodynia induced by peripheral inflammation. We also suggest that involvement of mGluR2/3 in the communication between glial cells and neurons takes part in the processing of nociceptive information.