Jinyu Wu

Characterization of the stress associated microRNAs in Glycine max by deep sequencing

BackgroundPlants involved in highly complex and well-coordinated systems have evolved a considerable degree of developmental plasticity, thus minimizing the damage caused by stress. MicroRNAs (miRNAs) have recently emerged as key regulators in gene regulation, developmental processes and stress tolerance in plants.ResultsIn this study, soybean miRNAs associated with stress responses (drought, salinity, and alkalinity) have been identified and analyzed in combination with deep sequencing technology and in-depth bioinformatics analysis. One hundred and thirty three conserved miRNAs representing 95 miRNA families were expressed in soybeans under three treatments. In addition, 71, 50, and 45 miRNAs are either uniquely or differently expressed under drought, salinity, and alkalinity, respectively, suggesting that many miRNAs are inducible and are differentially expressed in response to certain stress.ConclusionOur study has important implications for further identification of gene regulation under abiotic stresses and significantly contributes a complete profile of miRNAs in Glycine max.

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.

Gene expression profiling of soybean leaves and roots under salt, saline-alkali and drought stress by high-throughput Illumina sequencing.

Salt, saline-alkali and drought stresses are major environmental constraints for the production and yield of soybean worldwide. To identify genes responsible for stress tolerance, the transcriptional profiles of genes in leaves and roots of seedlings (two-leaf stage) of the soybean inbred line HJ-1 were examined after 48 h under various stress conditions; salt (120 mM NaCl), saline-alkali (70 mM NaCl and 50mM NaHCO(3)) and drought (2% PEG 8000). Gene expression at the transcriptional level was investigated using high-throughput Illumina sequencing technology and bioinformatics tools. Under salt, saline-alkali and drought stress, 874, 1897, and 535 genes, respectively, were up-regulated in leaves, and 1822, 1731 and 1690 genes, respectively, were up-regulated in roots, compared with expression in the corresponding organ in control plants. Comparisons among salt, saline-alkali and drought stress yielded similar results in terms of the percentage of genes classified into each GO category. Moreover, 69 genes differentially expressed in both organs with similar expression patterns clustered together in the taxonomic tree across all conditions. Furthermore, comparison of gene expression among salt, saline-alkali and drought treated plants revealed that genes associated with calcium-signaling and nucleic acid pathways were up-regulated in the responses to all three stresses, indicating a degree of cross-talk among these pathways. These results could provide new insights into the stress tolerance mechanisms of soybean.

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.

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

Sequencing and Genetic Variation of Multidrug Resistance Plasmids in Klebsiella pneumoniae

Background The development of multidrug resistance is a major problem in the treatment of pathogenic microorganisms by distinct antimicrobial agents. Characterizing the genetic variation among plasmids from different bacterial species or strains is a key step towards understanding the mechanism of virulence and their evolution. Results We applied a deep sequencing approach to 206 clinical strains of Klebsiella pneumoniae collected from 2002 to 2008 to understand the genetic variation of multidrug resistance plasmids, and to reveal the dynamic change of drug resistance over time. First, we sequenced three plasmids (70 Kb, 94 Kb, and 147 Kb) from a clonal strain of K. pneumoniae using Sanger sequencing. Using the Illumina sequencing technology, we obtained more than 17 million of short reads from two pooled plasmid samples. We mapped these short reads to the three reference plasmid sequences, and identified a large number of single nucleotide polymorphisms (SNPs) in these pooled plasmids. Many of these SNPs are present in drug-resistance genes. We also found that a significant fraction of short reads could not be mapped to the reference sequences, indicating a high degree of genetic variation among the collection of K. pneumoniae isolates. Moreover, we identified that plasmid conjugative transfer genes and antibiotic resistance genes are more likely to suffer from positive selection, as indicated by the elevated rates of nonsynonymous substitution. Conclusion These data represent the first large-scale study of genetic variation in multidrug resistance plasmids and provide insight into the mechanisms of plasmid diversification and the genetic basis of antibiotic resistance.

De Novo Transcriptome of Safflower and the Identification of Putative Genes for Oleosin and the Biosynthesis of Flavonoids

Safflower (Carthamus tinctorius L.) is one of the most extensively used oil crops in the world. However, little is known about how its compounds are synthesized at the genetic level. In this study, Solexa-based deep sequencing on seed, leaf and petal of safflower produced a de novo transcriptome consisting of 153,769 unigenes. We annotated 82,916 of the unigenes with gene annotation and assigned functional terms and specific pathways to a subset of them. Metabolic pathway analysis revealed that 23 unigenes were predicted to be responsible for the biosynthesis of flavonoids and 8 were characterized as seed-specific oleosins. In addition, a large number of differentially expressed unigenes, for example, those annotated as participating in anthocyanin and chalcone synthesis, were predicted to be involved in flavonoid biosynthesis pathways. In conclusion, the de novo transcriptome investigation of the unique transcripts provided candidate gene resources for studying oleosin-coding genes and for investigating genes related to flavonoid biosynthesis and metabolism in safflower.