Motivation: Clustering of protein–protein interaction networks is one of the most common approaches for predicting functional modules, protein complexes and protein functions. But, how well does clustering perform at these tasks?Results: We develop a general framework to assess how well computationally derived clusters in physical interactomes overlap functional modules derived via the Gene Ontology (GO). Using this framework, we evaluate six diverse network clustering algorithms using Saccharomyces cerevisiae and show that (i) the performances of these algorithms can differ substantially when run on the same network and (ii) their relative performances change depending upon the topological characteristics of the network under consideration. For the specific task of function prediction in S.cerevisiae, we demonstrate that, surprisingly, a simple non-clustering guilt-by-association approach outperforms widely used clustering-based approaches that annotate a protein with the overrepresented biological process and cellular component terms in its cluster; this is true over the range of clustering algorithms considered. Further analysis parameterizes performance based on the number of annotated proteins, and suggests when clustering approaches should be used for interactome functional analyses. Overall our results suggest a re-examination of when and how clustering approaches should be applied to physical interactomes, and establishes guidelines by which novel clustering approaches for biological networks should be justified and evaluated with respect to functional analysis.Contact: msingh@cs.princeton.eduSupplementary information: Supplementary data are available at Bioinformatics online.

BackgroundPiwi-interacting RNAs (piRNAs) are a recently discovered class of small non-coding RNAs whose best-understood function is to repress mobile element (ME) activity in animal germline. To date, nearly all piRNA studies have been conducted in model organisms and little is known about piRNA diversity, target specificity and biological function in human.ResultsHere we performed high-throughput sequencing of piRNAs from three human adult testis samples. We found that more than 81% of the ~17 million putative piRNAs mapped to ~6,000 piRNA-producing genomic clusters using a relaxed definition of clusters. A set of human protein-coding genes produces a relatively large amount of putative piRNAs from their 3’UTRs, and are significantly enriched for certain biological processes, suggestive of non-random sampling by the piRNA biogenesis machinery. Up to 16% of putative piRNAs mapped to a few hundred annotated long non-coding RNA (lncRNA) genes, suggesting that some lncRNA genes can act as piRNA precursors. Among major ME families, young families of LTR and endogenous retroviruses have a greater association with putative piRNAs than other MEs. In addition, piRNAs preferentially mapped to specific regions in the consensus sequences of several ME (sub)families and some piRNA mapping peaks showed patterns consistent with the “ping-pong” cycle of piRNA targeting and amplification.ConclusionsOverall our data provide a comprehensive analysis and improved annotation of human piRNAs in adult human testes and shed new light into the relationship of piRNAs with protein-coding genes, lncRNAs, and mobile genetic elements in human.Electronic supplementary materialThe online version of this article (doi:10.1186/1471-2164-15-545) contains supplementary material, which is available to authorized users.

Mutator phenotype tumors provide unique opportunities to unravel malignant progression because of various gene alterations acquired during clonal tumor evolution. Gastric carcinomas, which have been known to show frequent genetic instability, would be composed of initial gene alterations shared by most tumor areas and subsequent alterations restricted to particular tumor sites. To analyse the timing of genetic events, we examined separate sites of tumor tissue obtained from a given gastric carcinoma patient with microsatellite instability (MSI). Our study included 95 normal/tumor area pairs from 25 patients. Six of the 25 patients (24%) demonstrated various levels of MSI ranging from 7% (two of 30) to 97% (28 of 29) of markers tested in multiple tumor sites. Of the six patients, ®ve manifested frameshift mutations in a tract of ten deoxyadenosines within transforming growth factor b receptor type II and four demonstrated frameshift mutations in a tract of eight deoxyguanosines within BAX. These mutations were common to all tumor sites regardless of the various level of MSI phenotype, indicating initial events. Two of the six patients exhibited frameshift mutations in mononucleotide repeats of mismatch repair genes, hMSH3 and hMSH6, and the insulin-like growth factor II receptor in restricted tumor areas, indicating additional alterations. Insulin-like growth factor II receptor mutations appear to be caused by hMSH3 and hMSH6 mutations because the former mutations were con®ned to tumor portions with the latter two mismatch repair lesions. These results provide genetic progression evidence for gastric carcinomas of the mutator pathway. In this pathway, mismatch repair insuciency initially targets mononucleotide tracts of transforming growth factor b receptor type II and BAX. During tumorigenesis, primary mismatch repair failure may give rise to the secondary mismatch repair lesions, frameshift mutations of hMSH3 and hMSH6, which result in another tumorigenic mutation in the insulin-like growth factor II receptor.

Transposable elements (TEs) are one of the most important features of genome architecture, so their evolution and relationship with host defense mechanisms have been topics of intense study, especially in model systems such as Drosophila melanogaster. Recently, a novel small RNA-based defense mechanism in animals called the Piwi-interacting RNA (piRNA) pathway was discovered to form an adaptive defense mechanism against TEs. To investigate the relationship between piRNA and TE content between strains of a species, we sequenced piRNAs from 16 inbred lines of D. melanogaster from the Drosophila Genetic Reference Panel. Instead of a global correlation of piRNA expression and TE content, we found evidence for a host response through de novo piRNA production from novel TE insertions. Although approximately 20% of novel TE insertions induced de novo piRNA production, the abundance of de novo piRNAs was low and did not markedly affect the global pool of ovarian piRNAs. Our results provide new insights into the evolution of TEs and the piRNA system in an important model organism.

Numerous studies have suggested that hub proteins in the S. cerevisiae physical interaction network are more likely to be essential than other proteins. The proposed reasons underlying this observed relationship between topology and functioning have been subject to some controversy, with recent work suggesting that it arises due to the participation of hub proteins in essential complexes and processes. However, do these essential modules themselves have distinct network characteristics, and how do their essential proteins differ in their topological properties from their non-essential proteins? We aimed to advance our understanding of protein essentiality by analyzing proteins, complexes and processes within their broader functional context and by considering physical interactions both within and across complexes and biological processes. In agreement with the view that essentiality is a modular property, we found that the number of intracomplex or intraprocess interactions that a protein has is a better indicator of its essentiality than its overall number of interactions. Moreover, we found that within an essential complex, its essential proteins have on average more interactions, especially intracomplex interactions, than its non-essential proteins. Finally, we built a module-level interaction network and found that essential complexes and processes tend to have higher interaction degrees in this network than non-essential complexes and processes; that is, they exhibit a larger amount of functional cross-talk than their non-essential counterparts.

Epigenomic data from ENCODE can be used to associate specific combinations of chromatin marks with regulatory elements in the human genome. Hidden Markov models and the expectation-maximization (EM) algorithm are often used to analyze epigenomic data. However, the EM algorithm can have overfitting problems in data sets where the chromatin states show high class-imbalance and it is often slow to converge. Here we use spectral learning instead of EM and find that our software Spectacle overcame these problems. Furthermore, Spectacle is able to find enhancer subtypes not found by ChromHMM but strongly enriched in GWAS SNPs. Spectacle is available at https://github.com/jiminsong/Spectacle.Electronic supplementary materialThe online version of this article (doi:10.1186/s13059-015-0598-0) contains supplementary material, which is available to authorized users.

No abstract

PIWI-interacting RNAs (piRNAs) and their associated PIWI proteins play an important role in repressing transposable elements in animal germlines. However, little is known about the function of PIWI proteins and piRNAs in the developing brain. Here, we investigated the role of an important PIWI family member, Piwi-like protein 1 (Piwil1; also known as Miwi in mouse) in the developing mouse neocortex. Using a Piwil1 knock-out (Piwil1 KO) mouse strain, we found that Piwil1 is essential for several steps of neocorticogenesis, including neocortical cell cycle, neuron migration and dendritogenesis. Piwil1 deletion resulted in increased cell cycle re-entry at embryonic day 17 (E17) when predominantly intracortically projecting neurons are being produced. Prenatal Piwil1 deletion increased the number of Pax6+ radial glia at postnatal day 0 (P0). Furthermore, Piwil1 deletion disrupted migration of Satb2+ neurons within deep layers at E17, P0 and P7. Satb2+ neurons showed increased co-localization with Bcl11b (also known as Ctip2), marker of subcortically projecting neurons. Piwil1 knockouts had disrupted neocortical circuitry represented by thinning of the corpus callosum and altered dendritogenesis. We further investigated if Piwil1 deletion disrupted expression levels of neocortical piRNAs by small RNA-sequencing in neocortex. We did not find differential expression of piRNAs in the neocortices of Piwil1 KO, while differences were observed in other Piwil1 KO tissues. This result suggests that Piwil1 may act independently of piRNAs and have novel roles in higher cognitive centers, such as neocortex. In addition, we report a screen of piRNAs derived from tRNA fragments in developing neocortices. Our result is the first report of selective subsets of piRNAs and tRNA fragments in developing prenatal neocortices and helps clarify some outstanding questions about the role of the piRNA pathway in the brain.peer-reviewed)