The modulative effects of microcurrent electrical nerve stimulation on diabetic mice.

Diabetes (one of non-communicable diseases) is serious due to its complications, such like, cardiovascular ailments, neuropathy, nephropathy, retinopathy, wound gangrene and sexual impotence. Diabetes and associated chronic conditions are rapidly emerging as major health problems. In clinical, there were different drugs for diabetes treatment on different mechanisms. However, there were limited studies on the efficacy of electric stimulations on diabetes therapeutic application. In current study, we try to evaluate the effect of microcurrent electrical nerve stimulator (MENS) on diabetes modulation as an alternative medicine. A total of 36 male ICR mice of 6 weeks old were randomly divided into 4 groups [1] Control, [2] MENS only, [3] DM, [4] DM with MENS. During 8 weeks treatments, the diabetes-associated assessments included body weight, diet utilization, blood glucose measurement, other biochemistries and histopathological observations. The diabetes animal model induced by STZ had 180 mg/dl fasting blood glucose (GLU-AC) before MENS intervention. After 3 and 6 weeks administration, the GLU-AC of DM+MENS group significantly decreased 31.97% and 50.82% (P < 0.0001), respectively, as compared to DM group and the OGTT also demonstrated the similar significant results. The diabetic syndromes of polydipsia and polyphagia were also significantly ameliorated by MENS intervention. In other biochemical indexes, the glycated hemoglobin (HbA1c), hyperinsulinemia, liver functions (AST & ALT) and kidneys function (BUN & Creatinine) were also significantly mitigated by MENS under diabetes model. The histological observation also showed the MENS administration improved the diabetes-related pathological characteristics in liver, kidney and pancreas tissues. Our results suggest that administration of MENS could significantly improve diabetes animal model on blood sugar homeostasis, diabetic polydipsia, biochemistries, and tissue damage. In the health conditions, the MENS didn't exist obvious side effects on assessments. Therefore, the MENS could be potential on alternative medicine or supportive applications to future DM therapeutics.

Prevalence of intestinal parasitic infections among school children in capital areas of the Democratic Republic of São Tomé and Príncipe, West Africa.

BACKGROUND: Although the Democratic Republic of Sao Tome and Principe (DRSTP) has undertaken school children-based deworming programs against intestinal parasitic infections (IPIs) using a single dose of mebendazole annually since 2005, it remains unclear as to the outcome to date. The present study intends to investigate the recent IPIs status among school children living in capital areas of the DRSTP. METHODS: A total of 252 school children (121 boys and 131 girls) of grades 4 and 5 from 4 primary schools located in the capital areas participated in the present study and their fresh fecal specimens were examined for the presence of any parasites using the merthiolate-iodine-formaldehyde concentration method as conducted. RESULTS: The overall prevalence of IPIs was 64.7% (163/ 252). No significant gender difference in prevalence between boys (67.8%) and girls (61.8%) was found (p = 0.3). The majority of school children were infected with a single species of parasite (55.8%). Altogether, 12 different intestinal parasite species were identified in DRSTP school children, of which 9 species were pathogenic and the remaining 3 were non-pathogenic. CONCLUSION: Improving the detection method, sanitation facilities and personal hygiene as well as utilizing combined drugs are all important measures to greatly reduce IPIs in DRSTP school children.

Seroprevalence and Associated Risk Factors of Toxocariasis among College Students in Taipei City, Taiwan.

BACKGROUND: Infection by Toxocara spp. is known to be significantly associated with partial epilepsy. It has become popular for people to raise dogs/cats as pets and consume roasted meat/viscera, and the status of Toxocara spp. infection, epilepsy awareness, and associated risk factors among the general population are currently unknown in Taiwan. METHODS: A seroepidemiological investigation among 203 college students (CSs), consisting of 110 males and 93 females with an average age of 21.5 ± 1.2 years, was conducted in 2009 in Taipei City. A Western blot analysis based on excretory-secretory antigens derived from Toxocara canis larvae (TcESs) was applied to determine the positivity of serum immunoglobulin G antibodies. A self-administered questionnaire was also given to obtain information about demographic characteristics, epilepsy awareness, and risk factors. A logistic regression model was applied for the statistical analysis using SPSS software. RESULTS: The overall seropositive rate of Toxocara spp. infection was 8.4% (17/203). As to epilepsy awareness, a non-significantly higher seroprevalence was found in CSs who claimed to "know" about epilepsy compared to those who did not know (P > 0.05). CONCLUSIONS: It appears that appropriate educational programs are urgently needed to provide correct knowledge related to the prevention and control measures against Toxocara spp. infections to avoid potential threats by this parasite to the general population in Taiwan.

Gene Expression Profiling and Pathway Network Analysis Predicts a Novel Antitumor Function for a Botanical-Derived Drug, PG2.

PG2 is a botanical drug that is mostly composed of Astragalus polysaccharides (APS). Its role in hematopoiesis and relieving cancer-related fatigue has recently been clinically investigated in cancer patients. However, systematic analyses of its functions are still limited. The aim of this study was to use microarray-based expression profiling to evaluate the quality and consistency of PG2 from three different product batches and to study biological mechanisms of PG2. An integrative molecular analysis approach has been designed to examine significant PG2-induced signatures in HL-60 leukemia cells. A quantitative analysis of gene expression signatures was conducted for PG2 by hierarchical clustering of correlation coefficients. The results showed that PG2 product batches were consistent and of high quality. These batches were also functionally equivalent to each other with regard to how they modulated the immune and hematopoietic systems. Within the PG2 signature, there were five genes associated with doxorubicin: IL-8, MDM4, BCL2, PRODH2, and BIRC5. Moreover, the combination of PG2 and doxorubicin had a synergistic effect on induced cell death in HL-60 cells. Together with the bioinformatics-based approach, gene expression profiling provided a quantitative measurement for the quality and consistency of herbal medicines and revealed new roles (e.g., immune modulation) for PG2 in cancer treatment.

Differential network biology reveals a positive correlation between a novel protein-protein interaction and cancer cells migration.

This paper introduces a differential network biology for discovering tumor migration. We applied statistical methods to prioritize PPI candidates and an in situ proximity ligation assay to verify 67 endogenous PPIs among 21 interlinked pathways in two hepatocellular carcinoma (HCC) cells, Huh7 (minimally migratory cells) and Mahlavu (highly migratory cells). Differential network biology analysis was applied to determine the novel interaction, CRKL-FLT1, has a high centrality ranking, and the expression of this interaction is strongly correlated with the migratory ability of HCC and other cancer cell lines. Knockdown of CRKL and FLT1 in HCC cells leads to a decrease in cell migration. This study demonstrated that functional exploration of a disease network with differential network in interlinked pathways via PPIs can be used to discover tumor migration.

Genotyping of human neutrophil antigens (HNA) from whole genome sequencing data.

BACKGROUND: Neutrophil antigens are involved in a variety of clinical conditions including transfusion-related acute lung injury (TRALI) and other transfusion-related diseases. Recently, there are five characterized groups of human neutrophil antigen (HNA) systems, the HNA1 to 5. Characterization of all neutrophil antigens from whole genome sequencing (WGS) data may be accomplished for revealing complete genotyping formats of neutrophil antigens collectively at genome level with molecular variations which may respectively be revealed with available genotyping techniques for neutrophil antigens conventionally. RESULTS: We developed a computing method for the genotyping of human neutrophil antigens. Six samples from two families, available from the 1000 Genomes projects, were used for a HNA typing test. There are 500 ~ 3000 reads per sample filtered from the adopted human WGS datasets in order for identifying single nucleotide polymorphisms (SNPs) of neutrophil antigens. The visualization of read alignment shows that the yield reads from WGS dataset are enough to cover all of the SNP loci for the antigen system: HNA1, HNA3, HNA4 and HNA5. Consequently, our implemented Bioinformatics tool successfully revealed HNA types on all of the six samples including sequence-based typing (SBT) as well as PCR sequence-specific oligonucleotide probes (SSOP), PCR sequence-specific primers (SSP) and PCR restriction fragment length polymorphism (RFLP) along with parentage possibility. CONCLUSIONS: The next-generation sequencing technology strives to deliver affordable and non-biased sequencing results, hence the complete genotyping formats of HNA may be reported collectively from mining the output data of WGS. The study shows the feasibility of HNA genotyping through new WGS technologies. Our proposed algorithmic methodology is implemented in a HNATyping software package with user's guide available to the public at http://sourceforge.net/projects/hnatyping/.

SeqEntropy: genome-wide assessment of repeats for short read sequencing.

BACKGROUND: Recent studies on genome assembly from short-read sequencing data reported the limitation of this technology to reconstruct the entire genome even at very high depth coverage. We investigated the limitation from the perspective of information theory to evaluate the effect of repeats on short-read genome assembly using idealized (error-free) reads at different lengths. METHODOLOGY/PRINCIPAL FINDINGS: We define a metric H(k) to be the entropy of sequencing reads at a read length k and use the relative loss of entropy ΔH(k) to measure the impact of repeats for the reconstruction of whole-genome from sequences of length k. In our experiments, we found that entropy loss correlates well with de-novo assembly coverage of a genome, and a score of ΔH(k)>1% indicates a severe loss in genome reconstruction fidelity. The minimal read lengths to achieve ΔH(k)<1% are different for various organisms and are independent of the genome size. For example, in order to meet the threshold of ΔH(k)<1%, a read length of 60 bp is needed for the sequencing of human genome (3.2 10(9) bp) and 320 bp for the sequencing of fruit fly (1.8×10(8) bp). We also calculated the ΔH(k) scores for 2725 prokaryotic chromosomes and plasmids at several read lengths. Our results indicate that the levels of repeats in different genomes are diverse and the entropy of sequencing reads provides a measurement for the repeat structures. CONCLUSIONS/SIGNIFICANCE: The proposed entropy-based measurement, which can be calculated in seconds to minutes in most cases, provides a rapid quantitative evaluation on the limitation of idealized short-read genome sequencing. Moreover, the calculation can be parallelized to scale up to large euakryotic genomes. This approach may be useful to tune the sequencing parameters to achieve better genome assemblies when a closely related genome is already available.

EBARDenovo: highly accurate de novo assembly of RNA-Seq with efficient chimera-detection.

MOTIVATION: High-accuracy de novo assembly of the short sequencing reads from RNA-Seq technology is very challenging. We introduce a de novo assembly algorithm, EBARDenovo, which stands for Extension, Bridging And Repeat-sensing Denovo. This algorithm uses an efficient chimera-detection function to abrogate the effect of aberrant chimeric reads in RNA-Seq data. RESULTS: EBARDenovo resolves the complications of RNA-Seq assembly arising from sequencing errors, repetitive sequences and aberrant chimeric amplicons. In a series of assembly experiments, our algorithm is the most accurate among the examined programs, including de Bruijn graph assemblers, Trinity and Oases. AVAILABILITY AND IMPLEMENTATION: EBARDenovo is available at http://ebardenovo.sourceforge.net/. This software package (with patent pending) is free of charge for academic use only. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Analysis of protein-protein interactions in cross-talk pathways reveals CRKL protein as a novel prognostic marker in hepatocellular carcinoma.

Deciphering the network of signaling pathways in cancer via protein-protein interactions (PPIs) at the cellular level is a promising approach but remains incomplete. We used an in situ proximity ligation assay to identify and quantify 67 endogenous PPIs among 21 interlinked pathways in two hepatocellular carcinoma (HCC) cells, Huh7 (minimally migratory cells) and Mahlavu (highly migratory cells). We then applied a differential network biology analysis and determined that the novel interaction, CRKL-FLT1, has a high centrality ranking, and the expression of this interaction is strongly correlated with the migratory ability of HCC and other cancer cell lines. Knockdown of CRKL and FLT1 in HCC cells leads to a decrease in cell migration via ERK signaling and the epithelial-mesenchymal transition process. Our immunohistochemical analysis shows high expression levels of the CRKL and CRKL-FLT1 pair that strongly correlate with reduced disease-free and overall survival in HCC patient samples, and a multivariate analysis further established CRKL and the CRKL-FLT1 as novel prognosis markers. This study demonstrated that functional exploration of a disease network with interlinked pathways via PPIs can be used to discover novel biomarkers.

Analysis of schizophrenia and hepatocellular carcinoma genetic network with corresponding modularity and pathways: novel insights to the immune system.

BACKGROUND: Schizophrenic patients show lower incidences of cancer, implicating schizophrenia may be a protective factor against cancer. To study the genetic correlation between the two diseases, a specific PPI network was constructed with candidate genes of both schizophrenia and hepatocellular carcinoma. The network, designated schizophrenia-hepatocellular carcinoma network (SHCN), was analysed and cliques were identified as potential functional modules or complexes. The findings were compared with information from pathway databases such as KEGG, Reactome, PID and ConsensusPathDB. RESULTS: The functions of mediator genes from SHCN show immune system and cell cycle regulation have important roles in the eitology mechanism of schizophrenia. For example, the over-expressing schizophrenia candidate genes, SIRPB1, SYK and LCK, are responsible for signal transduction in cytokine production; immune responses involving IL-2 and TREM-1/DAP12 pathways are relevant for the etiology mechanism of schizophrenia. Novel treatments were proposed by searching the target genes of FDA approved drugs with genes in potential protein complexes and pathways. It was found that Vitamin A, retinoid acid and a few other immune response agents modulated by RARA and LCK genes may be potential treatments for both schizophrenia and hepatocellular carcinoma. CONCLUSIONS: This is the first study showing specific mediator genes in the SHCN which may suppress tumors. We also show that the schizophrenic protein interactions and modulation with cancer implicates the importance of immune system for etiology of schizophrenia.

The impact of normalization and phylogenetic information on estimating the distance for metagenomes.

Metagenomics enables the study of unculturable microorganisms in different environments directly. Discriminating between the compositional differences of metagenomes is an important and challenging problem. Several distance functions have been proposed to estimate the differences based on functional profiles or taxonomic distributions; however, the strengths and limitations of such functions are still unclear. Initially, we analyzed three well-known distance functions and found very little difference between them in the clustering of samples. This motivated us to incorporate suitable normalizations and phylogenetic information into the functions so that we could cluster samples from both real and synthetic data sets. The results indicate significant improvement in sample clustering over that derived by rank-based normalization with phylogenetic information, regardless of whether the samples are from real or synthetic microbiomes. Furthermore, our findings suggest that considering suitable normalizations and phylogenetic information is essential when designing distance functions for estimating the differences between metagenomes. We conclude that incorporating rank-based normalization with phylogenetic information into the distance functions helps achieve reliable clustering results.

Quantitative assessment of mitochondrial DNA copies from whole genome sequencing.

BACKGROUND: Mitochondrial dysfunction is associated with various aging diseases. The copy number of mtDNA in human cells may therefore be a potential biomarker for diagnostics of aging. Here we propose a new computational method for the accurate assessment of mtDNA copies from whole genome sequencing data. RESULTS: Two families of the human whole genome sequencing datasets from the HapMap and the 1000 Genomes projects were used for the accurate counting of mitochondrial DNA copy numbers. The results revealed the parental mitochondrial DNA copy numbers are significantly lower than that of their children in these samples. There are 8%~21% more copies of mtDNA in samples from the children than from their parents. The experiment demonstrated the possible correlations between the quantity of mitochondrial DNA and aging-related diseases. CONCLUSIONS: Since the next-generation sequencing technology strives to deliver affordable and non-biased sequencing results, accurate assessment of mtDNA copy numbers can be achieved effectively from the output of whole genome sequencing. We implemented the method as a software package MitoCounter with the source code and user's guide available to the public at http://sourceforge.net/projects/mitocounter/.

Gene expression-based chemical genomics identifies potential therapeutic drugs in hepatocellular carcinoma.

Hepatocellular carcinoma (HCC) is an aggressive tumor with a poor prognosis. Currently, only sorafenib is approved by the FDA for advanced HCC treatment; therefore, there is an urgent need to discover candidate therapeutic drugs for HCC. We hypothesized that if a drug signature could reverse, at least in part, the gene expression signature of HCC, it might have the potential to inhibit HCC-related pathways and thereby treat HCC. To test this hypothesis, we first built an integrative platform, the "Encyclopedia of Hepatocellular Carcinoma genes Online 2", dubbed EHCO2, to systematically collect, organize and compare the publicly available data from HCC studies. The resulting collection includes a total of 4,020 genes. To systematically query the Connectivity Map (CMap), which includes 6,100 drug-mediated expression profiles, we further designed various gene signature selection and enrichment methods, including a randomization technique, majority vote, and clique analysis. Subsequently, 28 out of 50 prioritized drugs, including tanespimycin, trichostatin A, thioguanosine, and several anti-psychotic drugs with anti-tumor activities, were validated via MTT cell viability assays and clonogenic assays in HCC cell lines. To accelerate their future clinical use, possibly through drug-repurposing, we selected two well-established drugs to test in mice, chlorpromazine and trifluoperazine. Both drugs inhibited orthotopic liver tumor growth. In conclusion, we successfully discovered and validated existing drugs for potential HCC therapeutic use with the pipeline of Connectivity Map analysis and lab verification, thereby suggesting the usefulness of this procedure to accelerate drug repurposing for HCC treatment.

MetaABC--an integrated metagenomics platform for data adjustment, binning and clustering.

SUMMARY: MetaABC is a metagenomic platform that integrates several binning tools coupled with methods for removing artifacts, analyzing unassigned reads and controlling sampling biases. It allows users to arrive at a better interpretation via series of distinct combinations of analysis tools. After execution, MetaABC provides outputs in various visual formats such as tables, pie and bar charts as well as clustering result diagrams. AVAILABILITY: MetaABC source code and documentation are available at http://bits2.iis.sinica.edu.tw/MetaABC/ CONTACT: dywang@gate.sinica.edu.tw; hktsai@iis.sinica.edu.tw SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

A quantitative analysis of monochromaticity in genetic interaction networks.

BACKGROUND: A genetic interaction refers to the deviation of phenotypes from the expected when perturbing two genes simultaneously. Studying genetic interactions help clarify relationships between genes, such as compensation and masking, and identify gene groups of functional modules. Recently, several genome-scale experiments for measuring quantitative (positive and negative) genetic interactions have been conducted. The results revealed that genes in the same module usually interact with each other in a consistent way (pure positive or negative); this phenomenon was designated as monochromaticity. Monochromaticity might be the underlying principle that can be utilized to unveil the modularity of cellular networks. However, no appropriate quantitative measurement for this phenomenon has been proposed. RESULTS: In this study, we propose the monochromatic index (MCI), which is able to quantitatively evaluate the monochromaticity of potential functional modules of genes, and the MCI was used to study genetic landscapes in different cellular subsystems. We demonstrated that MCI not only amend the deficiencies of MP-score but also properly incorporate the background effect. The results showed that not only within-complex but also between-complex connections present significant monochromatic tendency. Furthermore, we also found that significantly higher proportion of protein complexes are connected by negative genetic interactions in metabolic network, while transcription and translation system adopts relatively even number of positive and negative genetic interactions to link protein complexes. CONCLUSION: In summary, we demonstrate that MCI improves deficiencies suffered by MP-score, and can be used to evaluate monochromaticity in a quantitative manner. In addition, it also helps to unveil features of genetic landscapes in different cellular subsystems. Moreover, MCI can be easily applied to data produced by different types of genetic interaction methodologies such as Synthetic Genetic Array (SGA), and epistatic miniarray profile (E-MAP).

Construction and analysis of the protein-protein interaction networks for schizophrenia, bipolar disorder, and major depression.

BACKGROUND: Schizophrenia, bipolar disorder, and major depression are devastating mental diseases, each with distinctive yet overlapping epidemiologic characteristics. Microarray and proteomics data have revealed genes which expressed abnormally in patients. Several single nucleotide polymorphisms (SNPs) and mutations are associated with one or more of the three diseases. Nevertheless, there are few studies on the interactions among the disease-associated genes and proteins. RESULTS: This study, for the first time, incorporated microarray and protein-protein interaction (PPI) databases to construct the PPI network of abnormally expressed genes in postmortem brain samples of schizophrenia, bipolar disorder, and major depression patients. The samples were collected from Brodmann area (BA) 10 of the prefrontal cortex. Abnormally expressed disease genes were selected by t-tests comparing the disease and control samples. These genes were involved in housekeeping functions (e.g. translation, transcription, energy conversion, and metabolism), in brain specific functions (e.g. signal transduction, neuron cell differentiation, and cytoskeleton), or in stress responses (e.g. heat shocks and biotic stress).The diseases were interconnected through several "switchboard"-like nodes in the PPI network or shared abnormally expressed genes. A "core" functional module which consisted of a tightly knitted sub-network of clique-5 and -4s was also observed. These cliques were formed by 12 genes highly expressed in both disease and control samples. CONCLUSIONS: Several previously unidentified disease marker genes and drug targets, such as SBNO2 (schizophrenia), SEC24C (bipolar disorder), and SRRT (major depression), were identified based on statistical and topological analyses of the PPI network. The shared or interconnecting marker genes may explain the shared symptoms of the studied diseases. Furthermore, the "switchboard" genes, such as APP, UBC, and YWHAZ, are proposed as potential targets for developing new treatments due to their functional and topological significance.

Inferring a transcriptional regulatory network of the cytokinesis-related genes by network component analysis.

BACKGROUND: Network Component Analysis (NCA) is a network structure-driven framework for deducing regulatory signal dynamics. In contrast to principal component analysis, which can be employed to select the high-variance genes, NCA makes use of the connectivity structure from transcriptional regulatory networks to infer dynamics of transcription factor activities. Using the budding yeast Saccharomyces cerevisiae as a model system, we aim to deduce regulatory actions of cytokinesis-related genes, using precise spatial proximity (midbody) and/or temporal synchronicity (cytokinesis) to avoid full-scale computation from genome-wide databases. RESULTS: NCA was applied to infer regulatory actions of transcription factor activity from microarray data and partial transcription factor-gene connectivity information for cytokinesis-related genes, which were a subset of genome-wide datasets. No literature has so far discussed the inferred results through NCA are independent of the scale of the gene expression dataset. To avoid full-scale computation from genome-wide databases, four cytokinesis-related gene cases were selected for NCA by running computational analysis over the transcription factor database to confirm the approach being scale-free. The inferred dynamics of transcription factor activity through NCA were independent of the scale of the data matrix selected from the four cytokinesis-related gene sets. Moreover, the inferred regulatory actions were nearly identical to published observations for the selected cytokinesis-related genes in the budding yeast; namely, Mcm1, Ndd1, and Fkh2, which form a transcription factor complex to control expression of the CLB2 cluster (i.e. BUD4, CHS2, IQG1, and CDC5). CONCLUSION: In this study, using S. cerevisiae as a model system, NCA was successfully applied to infer similar regulatory actions of transcription factor activities from two various microarray databases and several partial transcription factor-gene connectivity datasets for selected cytokinesis-related genes independent of data sizes. The regulated action for four selected cytokinesis-related genes (BUD4, CHS2, IQG1, and CDC5) belongs to the M-phase or M/G1 phase, consistent with the empirical observations that in S. cerevisiae, the Mcm1-Ndd1-Fkh2 transcription factor complex can regulate expression of the cytokinesis-related genes BUD4, CHS2, IQG1, and CDC5. Since Bud4, Iqg1, and Cdc5 are highly conserved between human and yeast, results obtained from NCA for cytokinesis in the budding yeast can lead to a suggestion that human cells should have the transcription regulator(s) as the budding yeast Mcm1-Ndd1-Fkh2 transcription factor complex in controlling occurrence of cytokinesis.

From midbody protein-protein interaction network construction to novel regulators in cytokinesis.

Midbody, a transient organelle-like structure, is known as central for abscission and is indispensable for termination of cytokinesis. Here, we used the midbody proteome inventories to construct the potential midbody protein-protein interaction (PPI) network. To delineate novel regulators participating in cytokinesis, the z-score, a standard statistic score, rather than hub degree was implemented to prioritize the novel hubs. Of these hubs, KIAA0133, SEPT1, KIAA1377, and CRMP-1 were localized to the midbody, whereas HTR3A and ICAM2 were associated with the cleavage furrow as examined by immunofluorescence. Knockdown of SEPT1 and KIAA1377 resulted in increasing numbers of cytokinesis defect cells, suggesting these newly identified hubs play critical roles in cytokinesis progression. Moreover, ectopic expression of CRMP-1 mutant in which Aurora-A phosphorylation sites have been replaced with Ala results in a cytokinesis defect. This subproteome network construction not only sheds light on the intimate interactions of the midbody proteomes, but also prioritizes novel hubs or protein complexes that may govern the process of cytokinesis.

Cliques in mitotic spindle network bring kinetochore-associated complexes to form dependence pathway.

The mitotic spindle is an essential molecular machine for chromosome segregation during mitosis. Achieving a better understanding of its organization at the topological level remains a daunting task. To determine the functional connections among 137 mitotic spindle proteins, a protein-protein interaction network among queries was constructed. Many hub proteins, which connect more than one query and serve as highly plausible candidates for expanding the mitotic spindle proteome, are ranked by conventional degree centrality and a new subnetwork specificity score. Evaluation of the ranking results by literature reviews and empirical verification of SEPT6, a novel top-ranked hub, suggests that the subnetwork specificity score could enrich for putative spindle-related proteins. Topological analysis of this expanded network shows the presence of 30 3-cliques and six 4-cliques (fully connected subgraphs) that, respectively, reside in eight kinetochore-associated complexes, of which seven are evolution conserved. Notably, these complexes strikingly form dependence pathways for the assembly of the kinetochore complex. These analyses indicate the feasibility of using network topology, i.e. cliques, to uncover novel pathways to accelerate our understanding of potential biological processes.

POINeT: protein interactome with sub-network analysis and hub prioritization.

BACKGROUND: Protein-protein interactions (PPIs) are critical to every aspect of biological processes. Expansion of all PPIs from a set of given queries often results in a complex PPI network lacking spatiotemporal consideration. Moreover, the reliability of available PPI resources, which consist of low- and high-throughput data, for network construction remains a significant challenge. Even though a number of software tools are available to facilitate PPI network analysis, an integrated tool is crucial to alleviate the burden on querying across multiple web servers and software tools. RESULTS: We have constructed an integrated web service, POINeT, to simplify the process of PPI searching, analysis, and visualization. POINeT merges PPI and tissue-specific expression data from multiple resources. The tissue-specific PPIs and the numbers of research papers supporting the PPIs can be filtered with user-adjustable threshold values and are dynamically updated in the viewer. The network constructed in POINeT can be readily analyzed with, for example, the built-in centrality calculation module and an integrated network viewer. Nodes in global networks can also be ranked and filtered using various network analysis formulas, i.e., centralities. To prioritize the sub-network, we developed a ranking filtered method (S3) to uncover potential novel mediators in the midbody network. Several examples are provided to illustrate the functionality of POINeT. The network constructed from four schizophrenia risk markers suggests that EXOC4 might be a novel marker for this disease. Finally, a liver-specific PPI network has been filtered with adult and fetal liver expression profiles. CONCLUSION: The functionalities provided by POINeT are highly improved compared to previous version of POINT. POINeT enables the identification and ranking of potential novel genes involved in a sub-network. Combining with tissue-specific gene expression profiles, PPIs specific to selected tissues can be revealed. The straightforward interface of POINeT makes PPI search and analysis just a few clicks away. The modular design permits further functional enhancement without hampering the simplicity. POINeT is available at (http://poinet.bioinformatics.tw/).