Fuzzy methods for analysis of microarrays and networks

Bioinformatics involves analyses of biological data such as DNA sequences, microarrays and protein-protein interaction (PPI) networks. Its two main objectives are the identification of genes or proteins and the prediction of their functions. Biological data often contain uncertain and imprecise information. Fuzzy theory provides useful tools to deal with this type of information, hence has played an important role in analyses of biological data. In this thesis, we aim to develop some new fuzzy techniques and apply them on DNA microarrays and PPI networks. We will focus on three problems: (1) clustering of microarrays; (2) identification of disease-associated genes in microarrays; and (3) identification of protein complexes in PPI networks. The first part of the thesis aims to detect, by the fuzzy C-means (FCM) method, clustering structures in DNA microarrays corrupted by noise. Because of the presence of noise, some clustering structures found in random data may not have any biological significance. In this part, we propose to combine the FCM with the empirical mode decomposition (EMD) for clustering microarray data. The purpose of EMD is to reduce, preferably to remove, the effect of noise, resulting in what is known as denoised data. We call this method the fuzzy C-means method with empirical mode decomposition (FCM-EMD). We applied this method on yeast and serum microarrays, and the silhouette values are used for assessment of the quality of clustering. The results indicate that the clustering structures of denoised data are more reasonable, implying that genes have tighter association with their clusters. Furthermore we found that the estimation of the fuzzy parameter m, which is a difficult step, can be avoided to some extent by analysing denoised microarray data. The second part aims to identify disease-associated genes from DNA microarray data which are generated under different conditions, e.g., patients and normal people. We developed a type-2 fuzzy membership (FM) function for identification of diseaseassociated genes. This approach is applied to diabetes and lung cancer data, and a comparison with the original FM test was carried out. Among the ten best-ranked genes of diabetes identified by the type-2 FM test, seven genes have been confirmed as diabetes-associated genes according to gene description information in Gene Bank and the published literature. An additional gene is further identified. Among the ten best-ranked genes identified in lung cancer data, seven are confirmed that they are associated with lung cancer or its treatment. The type-2 FM-d values are significantly different, which makes the identifications more convincing than the original FM test. The third part of the thesis aims to identify protein complexes in large interaction networks. Identification of protein complexes is crucial to understand the principles of cellular organisation and to predict protein functions. In this part, we proposed a novel method which combines the fuzzy clustering method and interaction probability to identify the overlapping and non-overlapping community structures in PPI networks, then to detect protein complexes in these sub-networks. Our method is based on both the fuzzy relation model and the graph model. We applied the method on several PPI networks and compared with a popular protein complex identification method, the clique percolation method. For the same data, we detected more protein complexes. We also applied our method on two social networks. The results showed our method works well for detecting sub-networks and give a reasonable understanding of these communities.

A highly divergent South African geminivirus species illuminates the ancient evolutionary history of this family

Background. We have characterised a new highly divergent geminivirus species, Eragrostis curvula streak virus (ECSV), found infecting a hardy perennial South African wild grass. ECSV represents a new genus-level geminivirus lineage, and has a mixture of features normally associated with other specific geminivirus genera. Results. Whereas the ECSV genome is predicted to express a replication associated protein (Rep) from an unspliced complementary strand transcript that is most similar to those of begomoviruses, curtoviruses and topocuviruses, its Rep also contains what is apparently a canonical retinoblastoma related protein interaction motif such as that found in mastreviruses. Similarly, while ECSV has the same unusual TAAGATTCC virion strand replication origin nonanucleotide found in another recently described divergent geminivirus, Beet curly top Iran virus (BCTIV), the rest of the transcription and replication origin is structurally more similar to those found in begomoviruses and curtoviruses than it is to those found in BCTIV and mastreviruses. ECSV also has what might be a homologue of the begomovirus transcription activator protein gene found in begomoviruses, a mastrevirus-like coat protein gene and two intergenic regions. Conclusion. Although it superficially resembles a chimaera of geminiviruses from different genera, the ECSV genome is not obviously recombinant, implying that the features it shares with other geminiviruses are those that were probably present within the last common ancestor of these viruses. In addition to inferring how the ancestral geminivirus genome may have looked, we use the discovery of ECSV to refine various hypotheses regarding the recombinant origins of the major geminivirus lineages. © 2009 Varsani et al; licensee BioMed Central Ltd.

Receptor tyrosine kinases and their activation in melanoma

Receptor tyrosine kinases (RTKs) and their downstream signalling pathways have long been hypothesized to play key roles in melanoma development. A decade ago, evidence was derived largely from animal models, RTK expression studies and detection of activated RAS isoforms in a small fraction of melanomas. Predictions that overexpression of specific RTKs implied increased kinase activity and that some RTKs would show activating mutations in melanoma were largely untested. However, technological advances including rapid gene sequencing, siRNA methods and phospho-RTK arrays now give a more complete picture. Mutated forms of RTK genes including KIT, ERBB4, the EPH and FGFR families and others are known in melanoma. Additional over- or underexpressed RTKs and also protein tyrosine phosphatases (PTPs) have been reported, and activities measured. Complex interactions between RTKs and PTPs are implicated in the abnormal signalling driving aberrant growth and survival in malignant melanocytes, and indeed in normal melanocytic signalling including the response to ultraviolet radiation. Kinases are considered druggable targets, so characterization of global RTK activity in melanoma should assist the rational development of tyrosine kinase inhibitors for clinical use. © 2011 John Wiley & Sons A/S.

Epidermal growth factor receptors and cyclooxygenase-2 in the pathogenesis of non-small cell lung cancer : potential targets for chemoprevention and systemic therapy

The epidermal growth factor receptor (EGFR) is part of a family of plasma membrane receptor tyrosine kinases that control many important cellular functions, from growth and proliferation to cell death. Cyclooxygenase (COX)-2 is an enzyme which catalyses the conversion of arachidonic acid to prostagladins and thromboxane. It is induced by various inflammatory stimuli, including the pro-inflammatory cytokines, Interleukin (IL)-1β, Tumour Necrosis Factor (TNF)-α and IL-2. Both EGFR and COX-2 are over-expressed in non-small cell lung cancer (NSCLC) and have been implicated in the early stages of tumourigenesis. This paper considers their roles in the development and progression of lung cancer, their potential interactions, and reviews the recent progress in cancer therapies that are directed toward these targets. An increasing body of evidence suggests that selective inhibitors of both EGFR and COX-2 are potential therapeutic agents for the treatment of NSCLC, in the adjuvant, metastatic and chemopreventative settings. © 2002 Elsevier Science Ireland Ltd. All rights reserved.

miR-451 regulates zebrafish erythroid maturation in vivo via its target gata2

We demonstrate that in zebrafish, the microRNA miR-451 plays a crucial role in promoting erythroid maturation, in part via its target transcript gata2. Zebrafish miR-144 and miR-451 are processed from a single precursor transcript selectively expressed in erythrocytes. In contrast to other hematopoietic mutants, the ze-brafish mutant meunier (mnr) showed intact erythroid specification but diminished miR-144/451 expression. Although erythropoiesis initiated normally in mnr, erythrocyte maturation was morphologically retarded. Morpholino knockdown of miR-451 increased erythrocyte immaturity in wild-type embryos, and miR-451 RNA duplexes partially rescued erythroid maturation in mnr, demonstrating a requirement and role for miR-451 in erythro-cyte maturation. mnr provided a selectively miR-144/451-deficient background, facilitating studies to discern miRNA function and validate candidate targets. Among computer-predicted miR-451 targets potentially mediating these biologic effects, the pro-stem cell transcription factor gata2 was an attractive candidate. In vivo reporter assays validated the predicted miR-451/gata2-3'UTR interaction, gata2 down-regulation was delayed in miR-451-knockdown and mnr embryos, and gata2 knockdown partially restored erythroid maturation in mnr, collectively confirming gata2down-regulation as pivotal for miR-451-driven erythroid maturation. These studies define a new genetic pathway promoting erythroid maturation (mnr/miR-451/gata2) and provide a rare example of partial rescue of a mutant phenotype solely by miRNA overexpression. © 2009 by The American Society of Hematology.

The roles of plant dsRNA-binding proteins in RNAi-like pathways

Dicers are associated with double-stranded RNA-binding proteins (dsRBPs) in animals. In the plant, Arabidopsis, there are four dicer-like (DCL) proteins and five potential dsRBPs. These DCLs act redundantly and hierarchically. However, we show there is little or no redundancy or hierarchy amongst the DRBs in their DCL interactions. DCL1 operates exclusively with DRB1 to produce micro (mi)RNAs, DCL4 operates exclusively with DRB4 to produce trans-acting (ta) siRNAs and 21nt siRNAs from viral RNA. DCL2 and DCL3 produce viral siRNAs without requiring assistance from any dsRBP. DRB2, DRB3 and DRB5 appear unnecessary for mi-, tasi-, viral si-, or heterochromatinising siRNA production but act redundantly in a developmental pathway. © 2008 Federation of European Biochemical Societies.

Different in vitro activation methods for latent transforming growth factors (TGF)–β : considerable exogenous factors to promote higher mesenchymal-origin cell proliferation in a bioprocessing platform

Regenerative medicine includes two efficient techniques, namely tissue-engineering and cell-based therapy in order to repair tissue damage efficiently. Most importantly, huge numbers of autologous cells are required to deal these practices. Nevertheless, primary cells, from autologous tissue, grow very slowly while culturing in vitro; moreover, they lose their natural characteristics over prolonged culturing period. Transforming growth factors-beta (TGF-β) is a ubiquitous protein found biologically in its latent form, which prevents it from eliciting a response until conversion to its active form. In active form, TGF-β acts as a proliferative agent in many cell lines of mesenchymal origin in vitro. This article reviews on some of the important activation methods-physiochemical, enzyme-mediated, non-specific protein interaction mediated, and drug-induced- of TGF-β, which may be established as exogenous factors to be used in culturing medium to obtain extensive proliferation of primary cells.

Direct and quantitative detection of bacteriophage by "hearing" surface detachment using a quartz crystal microbalance

We show that it is possible to detect specifically adsorbed bacteriophage directly by breaking the interactions between proteins displayed on the phage coat and ligands immobilized on the surface of a quartz crystal microbalance (QCM). This is achieved through increasing the amplitude of oscillation of the QCM surface and sensitively detecting the acoustic emission produced when the bacteriophage detaches from the surface. There is no interference from nonspecifically adsorbed phage. The detection is quantitative over at least 5 orders of magnitude and is sensitive enough to detect as few as 20 phage. The method has potential as a sensitive and low-cost method for virus detection.

The binding of nuclear factors to the as-1 element in the CaMV 35S promoter is affected by cytosine methylation in vitro

Transcriptional gene silencing (TCS) is often associated with an increased level of cytosine methylation in the affected promoters. The effect of methylation of the cauliflower mosaic virus (CaMV) 35S promoter sequence on its binding to factors present in the nuclei was analyzed by electrophoretic mobility shift assays using extracts of petunia flowers. Specific DNA-protein interactions were detected in the region of the CaMV 35S promoter that contains the as-1 element and the region between -345 and -208. The binding of protein factor(s) to the as-1 element was influenced by cytosine methylation, whereas the binding to the region between -345 and -208 was unaffected. The results suggest that cytosine methylation of the as-1 element potentially affects the activity of the CaMV 35S promoter. © Georg Thieme Verlag KG Stuttgart.

PRMT2 and RORγ expression are associated with breast cancer survival outcomes

Protein arginine methyltransferases (PRMTs) methylate arginine residues on histones and target transcription factors that play critical roles in many cellular processes, including gene transcription, mRNA splicing, proliferation, and differentiation. Recent studies have linked PRMT-dependent epigenetic marks and modifications to carcinogenesis and metastasis in cancer. However, the role of PRMT2-dependent signaling in breast cancer remains obscure. We demonstrate PRMT2 mRNA expression was significantly decreased in breast cancer relative to normal breast. Gene expression profiling, Ingenuity and protein-protein interaction network analysis after PRMT2-short interfering RNA transfection into MCF-7 cells, revealed that PRMT2-dependent gene expression is involved in cell-cycle regulation and checkpoint control, chromosomal instability, DNA repair, and carcinogenesis. For example, PRMT2 depletion achieved the following: 1) increased p21 and decreased cyclinD1 expression in (several) breast cancer cell lines, 2) decreased cell migration, 3) induced an increase in nucleotide excision repair and homologous recombination DNA repair, and 4) increased the probability of distance metastasis free survival (DMFS). The expression of PRMT2 and retinoid-related orphan receptor-γ (RORγ) is inversely correlated in estrogen receptor-positive breast cancer and increased RORγ expression increases DMFS. Furthermore, we found decreased expression of the PRMT2-dependent signature is significantly associated with increased probability of DMFS. Finally, weighted gene coexpression network analysis demonstrated a significant correlation between PRMT2-dependent genes and cell-cycle checkpoint, kinetochore, and DNA repair circuits. Strikingly, these PRMT2-dependent circuits are correlated with pan-cancer metagene signatures associated with epithelial-mesenchymal transition and chromosomal instability. This study demonstrates the role and significant correlation between a histone methyltransferase (PRMT2)-dependent signature, RORγ, the cell-cycle regulation, DNA repair circuits, and breast cancer survival outcomes.

Genotoxicity of inorganic mercury salts based on disturbed microtubule function

This study investigated the hypothesis that the chromosomal genotoxicity of inorganic mercury results from interaction(s) with cytoskeletal proteins. Effects of Hg2+ salts on functional activities of tubulin and kinesin were investigated by determining tubulin assembly and kinesin-driven motility in cell-free systems. Hg2+ inhibits microtubule assembly at concentrations above 1 μM, and inhibition is complete at about 10 μM. In this range, the tubulin assembly is fully (up to 6 μM) or partially (∼6-10 μM) reversible. The inhibition of tubulin assembly by mercury is independent of the anion, chloride or nitrate. The no-observed-effect- concentration for inhibition of microtubule assembly in vitro was 1 μM Hg2+, the IC50 5.8 μM. Mercury(II) salts at the IC 50 concentrations partly inhibiting tubulin assembly did not cause the formation of aberrant microtubule structures. Effects of mercury salts on the functionality of the microtubule motility apparatus were studied with the motor protein kinesin. By using a "gliding assay" mimicking intracellular movement and transport processes in vitro, HgCl2 affected the gliding velocity of paclitaxel-stabilised microtubules in a clear dose-dependent manner. An apparent effect is detected at a concentration of 0.1 μM and a complete inhibition is reached at 1 μM. Cytotoxicity of mercury chloride was studied in V79 cells using neutral red uptake, showing an influence above 17 μM HgCl2. Between 15 and 20 μM HgCl2 there was a steep increase in cell toxicity. Both mercury chloride and mercury nitrate induced micronuclei concentration-dependently, starting at concentrations above 0.01 μM. CREST analyses on micronuclei formation in V79 cells demonstrated both clastogenic (CREST-negative) and aneugenic effects of Hg2+, with some preponderance of aneugenicity. A morphological effect of high Hg2+ concentrations (100 μM HgCl2) on the microtubule cytoskeleton was verified in V79 cells by immuno-fluorescence staining. The overall data are consistent with the concept that the chromosomal genotoxicity could be due to interaction of Hg2+ with the motor protein kinesin mediating cellular transport processes. Interactions of Hg 2+ with the tubulin shown by in vitro investigations could also partly influence intracellular microtubule functions leading, together with the effects on the kinesin, to an impaired chromosome distribution as shown by the micronucleus test.

Biological monitoring in workers in a nitrobenzene reduction plant : Haemoglobin versus serum albumin adducts

The high priority of monitoring workers exposed to nitrobenzene is a consequence of clear findings of experimental carcinogenicity of nitrobenzene and the associated evaluations by the International Agency for Research on Cancer. Eighty male employees of a nitrobenzene reduction plant, with potential skin contact with nitrobenzene and aniline, participated in a current medical surveillance programme. Blood samples were routinely taken and analysed for aniline, 4-aminodiphenyl (4-ADP) and benzidine adducts of haemoglobin (Hb) and human serum albumin (HSA). Also, levels of methaemoglobin (Met-Hb) and of carbon monoxide haemoglobin (CO-Hb) were monitored. Effects of smoking were straightforward. Using the rank sum test of Wilcoxon, we found that very clear-cut and statistically significant smoking effects (about 3-fold increases) were apparent on CO-Hb (P = 0.00085) and on the Hb adduct of 4-ADP (P = 0.0006). The mean aniline-Hb adduct level in smokers was 1.5 times higher than in non-smokers; the significance (P = 0.05375) was close to the 5% level. The strongest correlation was evident between the Hb and HSA adducts of aniline (rs = 0.846). Less pronounced correlations (but with P values < 0.02) appeared between aniline-Hb and 4-ADP-Hb adducts (rs = 0.388), between 4-ADP and 4-ADP-HSA adducts (rs = 0.373), and between 4-ADP-Hb and aniline-HSA adducts (rs = 0.275). In view of the proposal for additional use of the aniline-HSA adduct for biological monitoring, particularly in cases of acute overexposures or poisonings, the strong correlation of the Hb and HSA conjugates is noteworthy; the ratio aniline-HSA:aniline-Hb was 1:42 for the entire cohort.

Mutations in the gene encoding IFT dynein complex component WDR34 cause jeune asphyxiating thoracic dystrophy

Bidirectional (anterograde and retrograde) motor-based intraflagellar transport (IFT) governs cargo transport and delivery processes that are essential for primary cilia growth and maintenance and for hedgehog signaling functions. The IFT dynein-2 motor complex that regulates ciliary retrograde protein transport contains a heavy chain dynein ATPase/motor subunit, DYNC2H1, along with other less well functionally defined subunits. Deficiency of IFT proteins, including DYNC2H1, underlies a spectrum of skeletal ciliopathies. Here, by using exome sequencing and a targeted next-generation sequencing panel, we identified a total of 11 mutations in WDR34 in 9 families with the clinical diagnosis of Jeune syndrome (asphyxiating thoracic dystrophy). WDR34 encodes a WD40 repeat-containing protein orthologous to Chlamydomonas FAP133, a dynein intermediate chain associated with the retrograde intraflagellar transport motor. Three-dimensional protein modeling suggests that the identified mutations all affect residues critical for WDR34 protein-protein interactions. We find that WDR34 concentrates around the centrioles and basal bodies in mammalian cells, also showing axonemal staining. WDR34 coimmunoprecipitates with the dynein-1 light chain DYNLL1 in vitro, and mining of proteomics data suggests that WDR34 could represent a previously unrecognized link between the cytoplasmic dynein-1 and IFT dynein-2 motors. Together, these data show that WDR34 is critical for ciliary functions essential to normal development and survival, most probably as a previously unrecognized component of the mammalian dynein-IFT machinery.

Compound heterozygous mutations in RIPPLY2 associated with vertebral segmentation defects

Segmentation defects of the vertebrae (SDV) are caused by aberrant somite formation during embryogenesis and result in irregular formation of the vertebrae and ribs. The Notch signal transduction pathway plays a critical role in somite formation and patterning in model vertebrates. In humans, mutations in several genes involved in the Notch pathway are associated with SDV, with both autosomal recessive (MESP2, DLL3, LFNG, HES7) and autosomal dominant (TBX6) inheritance. However, many individuals with SDV do not carry mutations in these genes. Using whole-exome capture and massive parallel sequencing, we identified compound heterozygous mutations in RIPPLY2 in two brothers with multiple regional SDV, with appropriate familial segregation. One novel mutation (c.A238T:p.Arg80*) introduces a premature stop codon. In transiently transfected C2C12 mouse myoblasts, the RIPPLY2 mutant protein demonstrated impaired transcriptional repression activity compared with wild-type RIPPLY2 despite similar levels of expression. The other mutation (c.240-4T>G), with minor allele frequency <0.002, lies in the highly conserved splice site consensus sequence 5' to the terminal exon. Ripply2 has a well-established role in somitogenesis and vertebral column formation, interacting at both gene and protein levels with SDV-associated Mesp2 and Tbx6. We conclude that compound heterozygous mutations in RIPPLY2 are associated with SDV, a new gene for this condition. © The Author 2014.

Prediction of Candidate Primary Immunodeficiency Disease Genes Using a Support Vector Machine Learning Approach

Screening and early identification of primary immunodeficiency disease (PID) genes is a major challenge for physicians. Many resources have catalogued molecular alterations in known PID genes along with their associated clinical and immunological phenotypes. However, these resources do not assist in identifying candidate PID genes. We have recently developed a platform designated Resource of Asian PDIs, which hosts information pertaining to molecular alterations, protein-protein interaction networks, mouse studies and microarray gene expression profiling of all known PID genes. Using this resource as a discovery tool, we describe the development of an algorithm for prediction of candidate PID genes. Using a support vector machine learning approach, we have predicted 1442 candidate PID genes using 69 binary features of 148 known PID genes and 3162 non-PID genes as a training data set. The power of this approach is illustrated by the fact that six of the predicted genes have recently been experimentally confirmed to be PID genes. The remaining genes in this predicted data set represent attractive candidates for testing in patients where the etiology cannot be ascribed to any of the known PID genes.