The advent of next-generation sequencing technologies is accompanied with the development of many whole-genome sequence assembly methods and software, especially for de novo fragment assembly. Due to the poor knowledge about the applicability and performance of these software tools, choosing a befitting assembler becomes a tough task. Here, we provide the information of adaptivity for each program, then above all, compare the performance of eight distinct tools against eight groups of simulated datasets from Solexa sequencing platform. Considering the computational time, maximum random access memory (RAM) occupancy, assembly accuracy and integrity, our study indicate that string-based assemblers, overlap-layout-consensus (OLC) assemblers are well-suited for very short reads and longer reads of small genomes respectively. For large datasets of more than hundred millions of short reads, De Bruijn graph-based assemblers would be more appropriate. In terms of software implementation, string-based assemblers are superior to graph-based ones, of which SOAPdenovo is complex for the creation of configuration file. Our comparison study will assist researchers in selecting a well-suited assembler and offer essential information for the improvement of existing assemblers or the developing of novel assemblers.

The origin of new genes with novel functions creates genetic and phenotypic diversity in organisms. To acquire functional roles, new genes must integrate into ancestral gene-gene interaction (GGI) networks. The mechanisms by which new genes are integrated into ancestral networks, and their evolutionary significance, are yet to be characterized. Herein, we present a study investigating the rates and patterns of new gene-driven evolution of GGI networks in the human and mouse genomes. We examine the network topological and functional evolution of new genes that originated at various stages in the human and mouse lineages by constructing and analyzing three different GGI datasets. We find a large number of new genes integrated into GGI networks throughout vertebrate evolution. These genes experienced a gradual integration process into GGI networks, starting on the network periphery and gradually becoming highly connected hubs, and acquiring pleiotropic and essential functions. We identify a few human lineage-specific hub genes that have evolved brain development-related functions. Finally, we explore the possible underlying mechanisms driving the GGI network evolution and the observed patterns of new gene integration process. Our results unveil a remarkable network topological integration process of new genes: over 5000 new genes were integrated into the ancestral GGI networks of human and mouse; new genes gradually acquire increasing number of gene partners; some human-specific genes evolved into hub structure with critical phenotypic effects. Our data cast new conceptual insights into the evolution of genetic networks.

MicroRNAs (miRNAs) are a class of non-coding regulatory RNAs approximately 22 nucleotides in length that play a role in a wide range of biological processes. Abnormal miRNA function has been implicated in various human cancers including prostate cancer (PCa). Altered miRNA expression may serve as a biomarker for cancer diagnosis and treatment. However, limited data are available on the role of cancer-specific miRNAs. Integrative computational bioinformatics approaches are effective for the detection of potential outlier miRNAs in cancer.

Screening key microRNAs for castration-resistant prostate cancer based on miRNA/mRNA functional synergistic network

Background: Clear cell renal cell carcinoma (ccRCC) represents the most invasive and common adult kidney neoplasm. Mounting evidence suggests that microRNAs (miRNAs) are important regulators of gene expression. But their function in tumourigenesis in this tumour type remains elusive. With the development of high throughput technologies such as microarrays and NGS, aberrant miRNA expression has been widely observed in ccRCC. Systematic and integrative analysis of multiple microRNA expression datasets may reveal potential mechanisms by which microRNAs contribute to ccRCC pathogenesis. Methods: We collected 5 public microRNA expression datasets in ccRCC versus non-matching normal renal tissues from GEO database and published literatures. We analyzed these data sets with an integrated bioinformatics framework to identify expression signatures. The framework incorporates a novel statistic method for abnormal gene expression detection and an in-house developed predictor to assess the regulatory activity of microRNAs. We then mapped target genes of DE-miRNAs to different databases, such as GO, KEGG, GeneGo etc, for functional enrichment analysis. Results: Using this framework we identified a consistent panel of eleven deregulated miRNAs shared by five independent datasets that can distinguish normal kidney tissues from ccRCC. After comparison with 3 RNA-seq based microRNA profiling studies, we found that our data correlated well with the results of next generation sequencing. We also discovered 14 novel molecular pathways that are likely to play a role in the tumourigenesis of ccRCC. Conclusions: The integrative framework described in this paper greatly improves the inter-dataset consistency of microRNA expression signatures. Consensus expression profile should be identified at pathway or network level to address the heterogeneity of cancer. The DE-miRNA signature and novel pathways identified herein could provide potential biomarkers for ccRCC that await further validation.

Next-generation sequencing (NGS) technology has rapidly advanced and generated the massive data volumes. To align and map the NGS data, biologists often randomly select a number of aligners without concerning their suitable feature, high performance, and high accuracy as well as sequence variations and polymorphisms existing on reference genome. This study aims to systematically evaluate and compare the capability of multiple aligners for NGS data analysis. To explore this capability, we firstly performed alignment algorithms comparison and classification. We further used long-read and short-read datasets from both real-life andin silicoNGS data for comparative analysis and evaluation of these aligners focusing on three criteria, namely, application-specific alignment feature, computational performance, and alignment accuracy. Our study demonstrated the overall evaluation and comparison of multiple aligners for NGS data analysis. This serves as an important guiding resource for biologists to gain further insight into suitable selection of aligners for specific and broad applications.