Biomolecular network motif counting and discovery by color coding

Protein-protein interaction (PPI) networks of many organisms share global topological features such as degree distribution, k-hop reachability, betweenness and closeness. Yet, some of these networks can differ significantly from the others in terms of local structures: e.g. the number of specific network motifs can vary significantly among PPI networks. Counting the number of network motifs provides a major challenge to compare biomolecular networks. Recently developed algorithms have been able to count the number of induced occurrences of subgraphs with k < or = 7 vertices. Yet no practical algorithm exists for counting non-induced occurrences, or counting subgraphs with k > or = 8 vertices. Counting non-induced occurrences of network motifs is not only challenging but also quite desirable as available PPI networks include several false interactions and miss many others. In this article, we show how to apply the 'color coding' technique for counting non-induced occurrences of subgraph topologies in the form of trees and bounded treewidth subgraphs. Our algorithm can count all occurrences of motif G' with k vertices in a network G with n vertices in time polynomial with n, provided k = O(log n). We use our algorithm to obtain 'treelet' distributions for k < or = 10 of available PPI networks of unicellular organisms (Saccharomyces cerevisiae Escherichia coli and Helicobacter Pyloris), which are all quite similar, and a multicellular organism (Caenorhabditis elegans) which is significantly different. Furthermore, the treelet distribution of the unicellular organisms are similar to that obtained by the 'duplication model' but are quite different from that of the 'preferential attachment model'. The treelet distribution is robust w.r.t. sparsification with bait/edge coverage of 70% but differences can be observed when bait/edge coverage drops to 50%.     

Optimal pooling for genome re-sequencing with ultra-high-throughput short-read technologies

New generation sequencing technologies offer unique opportunities and challenges for re-sequencing studies. In this article, we focus on re-sequencing experiments using the Solexa technology, based on bacterial artificial chromosome (BAC) clones, and address an experimental design problem. In these specific experiments, approximate coordinates of the BACs on a reference genome are known, and fine-scale differences between the BAC sequences and the reference are of interest. The high-throughput characteristics of the sequencing technology makes it possible to multiplex BAC sequencing experiments by pooling BACs for a cost-effective operation. However, the way BACs are pooled in such re-sequencing experiments has an effect on the downstream analysis of the generated data, mostly due to subsequences common to multiple BACs. The experimental design strategy we develop in this article offers combinatorial solutions based on approximation algorithms for the well-known max n-cut problem and the related max n-section problem on hypergraphs. Our algorithms, when applied to a number of sample cases give more than a 2-fold performance improvement over random partitioning.     

Whole genome sequencing of Turkish genomes reveals functional private alleles and impact of genetic interactions with Europe,Asia and Africa

Background

Turkey is a crossroads of major population movements throughout history and has been a hotspot of cultural interactions. Several studies have investigated the complex population history of Turkey through a limited set of genetic markers. However, to date, there have been no studies to assess the genetic variation at the whole genome level using whole genome sequencing. Here, we present whole genome sequences of 16 Turkish individuals resequenced at high coverage (32 × -48×).

Results

We show that the genetic variation of the contemporary Turkish population clusters with South European populations, as expected, but also shows signatures of relatively recent contribution from ancestral East Asian populations. In addition, we document a significant enrichment of non-synonymous private alleles, consistent with recent observations in European populations. A number of variants associated with skin color and total cholesterol levels show frequency differentiation between the Turkish populations and European populations. Furthermore, we have analyzed the 17q21.31 inversion polymorphism region (MAPT locus) and found increased allele frequency of 31.25% for H1/H2 inversion polymorphism when compared to European populations that show about 25% of allele frequency.

Conclusion

This study provides the first map of common genetic variation from 16 western Asian individuals and thus helps fill an important geographical gap in analyzing natural human variation and human migration. Our data will help develop population-specific experimental designs for studies investigating disease associations and demographic history in Turkey.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-963) contains supplementary material, which is available to authorized users.     

Cortical spreading depression confounds concentration-dependent pial arteriolar dilation during N-methyl-D-aspartate superfusion

Pial arterioles do not express N-methyl-D-aspartate (NMDA) receptors but dilate in response to topical NMDA application. We explored the mechanism underlying NMDA-mediated responses in murine pial arterioles (11-31 microm), using a closed cranial window preparation, and found that arteriolar dilation was not concentration dependent. Pial arteriolar diameter abruptly increased within 3 min of superfusing 50 or 100 microM NMDA. Dilation reached a peak within 1 min (46 +/- 14%) and then declined to a plateau (28 +/- 13%) for the duration of superfusion. Whereas a higher concentration (200 microM) did not produce further dilation, lower concentrations (1-10 microM) did not dilate the arterioles at all. MK-801 (10 microM) abrogated the dilation response, whereas Nomega-nitro-L-arginine (1 mM) attenuated the peak and abolished the sustained dilation during NMDA superfusion. We determined that NMDA-induced pial arteriolar responses were evoked by cortical spreading depression, because abrupt vasodilation during 50 or 100 microM NMDA superfusion was associated with a large negative slow potential shift and electrocorticogram suppression that spread from the superfusion window to distant cortical areas. Our data suggest that the responses of pial arterioles to NMDA are caused in part by neurovascular coupling due to cortical spreading depression.     

MELK-T1, a small-molecule inhibitor of protein kinase MELK,decreases DNA-damage tolerance in proliferating cancer cells

Maternal embryonic leucine zipper kinase (MELK), a serine/threonine protein kinase, has oncogenic properties and is overexpressed in many cancer cells. The oncogenic function of MELK is attributed to its capacity to disable critical cell-cycle checkpoints and reduce replication stress. Most functional studies have relied on the use of siRNA/shRNA-mediated gene silencing. In the present study, we have explored the biological function of MELK using MELK-T1, a novel and selective small-molecule inhibitor. Strikingly, MELK-T1 triggered a rapid and proteasome-dependent degradation of the MELK protein. Treatment of MCF-7 (Michigan Cancer Foundation-7) breast adenocarcinoma cells with MELK-T1 induced the accumulation of stalled replication forks and double-strand breaks that culminated in a replicative senescence phenotype. This phenotype correlated with a rapid and long-lasting ataxia telangiectasia-mutated (ATM) activation and phosphorylation of checkpoint kinase 2 (CHK2). Furthermore, MELK-T1 induced a strong phosphorylation of p53 (cellular tumour antigen p53), a prolonged up-regulation of p21 (cyclin-dependent kinase inhibitor 1) and a down-regulation of FOXM1 (Forkhead Box M1) target genes. Our data indicate that MELK is a key stimulator of proliferation by its ability to increase the threshold for DNA-damage tolerance (DDT). Thus, targeting MELK by the inhibition of both its catalytic activity and its protein stability might sensitize tumours to DNA-damaging agents or radiation therapy by lowering the DNA-damage threshold.