Meander: visually exploring the structural variome using space-filling curves

The introduction of next generation sequencing methods in genome studies has made it possible to shift research from a gene-centric approach to a genome wide view. Although methods and tools to detect single nucleotide polymorphisms are becoming more mature, methods to identify and visualize structural variation (SV) are still in their infancy. Most genome browsers can only compare a given sequence to a reference genome; therefore, direct comparison of multiple individuals still remains a challenge. Therefore, the implementation of efficient approaches to explore and visualize SVs and directly compare two or more individuals is desirable. In this article, we present a visualization approach that uses space-filling Hilbert curves to explore SVs based on both read-depth and pair-end information. An interactive open-source Java application, called Meander, implements the proposed methodology, and its functionality is demonstrated using two cases. With Meander, users can explore variations at different levels of resolution and simultaneously compare up to four different individuals against a common reference. The application was developed using Java version 1.6 and Processing.org and can be run on any platform. It can be found at http://homes.esat.kuleuven.be/~bioiuser/meander.     

Phosphorylation of CK1delta: identification of Ser370 as the major phosphorylation site targeted by PKA in vitro and in vivo

The involvement of CK1 (casein kinase 1) delta in the regulation of multiple cellular processes implies a tight regulation of its activity on many different levels. At the protein level, reversible phosphorylation plays an important role in modulating the activity of CK1delta. In the present study, we show that PKA (cAMP-dependent protein kinase), Akt (protein kinase B), CLK2 (CDC-like kinase 2) and PKC (protein kinase C) alpha all phosphorylate CK1delta. PKA was identified as the major cellular CK1deltaCK (CK1delta C-terminal-targeted protein kinase) for the phosphorylation of CK1delta in vitro and in vivo. This was implied by the following evidence: PKA was detectable in the CK1deltaCK peak fraction of fractionated MiaPaCa-2 cell extracts, PKA shared nearly identical kinetic properties with those of CK1deltaCK, and both PKA and CK1deltaCK phosphorylated CK1delta at Ser370 in vitro. Furthermore, phosphorylation of CK1delta by PKA decreased substrate phosphorylation of CK1delta in vitro. Mutation of Ser370 to alanine increased the phosphorylation affinity of CK1delta for beta-casein and the GST (gluthatione S-transferase)-p53 1-64 fusion protein in vitro and enhanced the formation of an ectopic dorsal axis during Xenopus laevis development. Anchoring of PKA and CK1delta to centrosomes was mediated by AKAP (A-kinase-anchoring protein) 450. Interestingly, pre-incubation of MiaPaCa-2 cells with the synthetic peptide St-Ht31, which prevents binding between AKAP450 and the regulatory subunit RII of PKA, resulted in a 6-fold increase in the activity of CK1delta. In summary, we conclude that PKA phosphorylates CK1delta, predominantly at Ser370 in vitro and in vivo, and that site-specific phosphorylation of CK1delta by PKA plays an important role in modulating CK1delta-dependent processes.     

Relationship of CSF neurotransmitter metabolite levels to disease severity and disability in multiple sclerosis

Axonal degeneration and brain tissue loss occur during disease progression in multiple sclerosis (MS) and are expected to influence neurotransmitter activities, with consequences on neurologic and psychiatric symptomatology. We searched for relationships of disease duration, disability, and severity of MS patients to CSF levels of the major metabolites of noradrenaline, dopamine, and serotonin, MHPG, methoxyhydroxyphenylglycol (MHPG), homovanillic acid, and 5-hydroxyindoleacetic acid (5-HIAA), respectively, in 39 patients with relapsing–remitting (RR) MS in remission, and 26 patients with progressive (PR) MS. Disability and Disease Severity were assessed by the Expanded Disability Status Scale (EDSS) and the Multiple Sclerosis Severity Score (MSSS). Compared with the levels of 50 control subjects, MHPG levels were not different in either MS group, correlated negatively to duration of illness and number of relapses in the RRMS group, but not to EDSS score or to MSSS. Homovanillic acid levels were significantly lower only in the PRMS group, with a negative correlation to duration of illness, and a strong negative correlation to EDSS score, but not to MSSS. 5-HIAA was significantly lower in both RRMS and PRMS groups. In the RRMS group, 5-HIAA levels were negatively related to EDSS and to MSSS. Multiple regression analyses revealed a significant association of MHPG to duration of illness, and a strong negative association of 5-HIAA to MSSS rather than to EDSS. The strong negative correlation of MSSS to CSF 5-HIAA levels in RRMS group of patients indicates that deficits in central serotonergic activity are related to the rate of disability accumulation in RRMS, and could be linked to the reported reduction of disease activity by serotonergic drugs.     

Mathematical modeling of cancer cell invasion of tissue: biological insight from mathematical analysis and computational simulation

The ability of cancer cells to break out of tissue compartments and invade locally gives solid tumours a defining deadly characteristic. One of the first steps of invasion is the remodelling of the surrounding tissue or extracellular matrix (ECM) and a major part of this process is the over-expression of proteolytic enzymes, such as the urokinase-type plasminogen activator (uPA) and matrix metalloproteinases (MMPs), by the cancer cells to break down ECM proteins. Degradation of the matrix enables the cancer cells to migrate through the tissue and subsequently to spread to secondary sites in the body, a process known as metastasis. In this paper we undertake an analysis of a mathematical model of cancer cell invasion of tissue, or ECM, which focuses on the role of the urokinase plasminogen activation system. The model consists of a system of five reaction-diffusion-taxis partial differential equations describing the interactions between cancer cells, uPA, uPA inhibitors, plasmin and the host tissue. Cancer cells react chemotactically and haptotactically to the spatio-temporal effects of the uPA system. The results obtained from computational simulations carried out on the model equations produce dynamic heterogeneous spatio-temporal solutions and using linear stability analysis we show that this is caused by a taxis-driven instability of a spatially homogeneous steady-state. Finally we consider the biological implications of the model results, draw parallels with clinical samples and laboratory based models of cancer cell invasion using three-dimensional invasion assay, and go on to discuss future development of the model.     

Extensive polymorphism of the major histocompatibility complex DRA gene in Balkan donkeys: perspectives on selection and genealogy

The major histocompatibility complex (MHC) contains genes important for immune response in mammals, and these genes exhibit high polymorphism and diversity. The DRA gene, a member of the MHC class II family, is highly conserved across a large number of mammalian species, but it displays exceptionally rich sequence variations in Equidae members. We analyzed allelic polymorphism of the DRA locus in 248 donkeys sampled across the Balkan Peninsula (Albania, Bulgaria, Croatia, Macedonia, Greece and Montenegro). Five known alleles and two new alleles were identified. The new allele Eqas‐DRA*0601 was found to carry a synonymous mutation, and new allele Eqas‐DRA*0701, a non‐synonymous mutation. We further analyzed the historical selection and allele genealogy at the DRA locus in equids. Signals of positive selection obtained by various tests were ambiguous. A conservative conclusion is that DRA polymorphism occurred relatively recently and that positive selection has been acting on the DRA locus for a relatively brief period.     

Revealing Pathway Dynamics in Heart Diseases by Analyzing Multiple Differential Networks

Development of heart diseases is driven by dynamic changes in both the activity and connectivity of gene pathways. Understanding these dynamic events is critical for understanding pathogenic mechanisms and development of effective treatment. Currently, there is a lack of computational methods that enable analysis of multiple gene networks, each of which exhibits differential activity compared to the network of the baseline/healthy condition. We describe the iMDM algorithm to identify both unique and shared gene modules across multiple differential co-expression networks, termed M-DMs (multiple differential modules). We applied iMDM to a time-course RNA-Seq dataset generated using a murine heart failure model generated on two genotypes. We showed that iMDM achieves higher accuracy in inferring gene modules compared to using single or multiple co-expression networks. We found that condition-specific M-DMs exhibit differential activities, mediate different biological processes, and are enriched for genes with known cardiovascular phenotypes. By analyzing M-DMs that are present in multiple conditions, we revealed dynamic changes in pathway activity and connectivity across heart failure conditions. We further showed that module dynamics were correlated with the dynamics of disease phenotypes during the development of heart failure. Thus, pathway dynamics is a powerful measure for understanding pathogenesis. iMDM provides a principled way to dissect the dynamics of gene pathways and its relationship to the dynamics of disease phenotype. With the exponential growth of omics data, our method can aid in generating systems-level insights into disease progression.     

Brugia malayi Antigen (BmA) Inhibits HIV-1 Trans-Infection but Neither BmA nor ES-62 Alter HIV-1 Infectivity of DC Induced CD4+ Th-Cells

One of the hallmarks of HIV-1 disease is the association of heightened CD4⁺ T-cell activation with HIV-1 replication. Parasitic helminths including filarial nematodes have evolved numerous and complex mechanisms to skew, dampen and evade human immune responses suggesting that HIV-1 infection may be modulated in co-infected individuals. Here we studied the effects of two filarial nematode products, adult worm antigen from Brugia malayi (BmA) and excretory-secretory product 62 (ES-62) from Acanthocheilonema viteae on HIV-1 infection in vitro. Neither BmA nor ES-62 influenced HIV-1 replication in CD4⁺ enriched T-cells, with either a CCR5- or CXCR4-using virus. BmA, but not ES-62, had the capacity to bind the C-type lectin dendritic cell-specific intercellular adhesion molecule-3-grabbing non-integrin (DC-SIGN) thereby inhibiting HIV-1 trans-infection of CD4⁺ enriched T-cells. As for their effect on DCs, neither BmA nor ES-62 could enhance or inhibit DC maturation as determined by CD83, CD86 and HLA-DR expression, or the production of IL-6, IL-10, IL-12 and TNF-α. As expected, due to the unaltered DC phenotype, no differences were found in CD4⁺ T helper (Th) cell phenotypes induced by DCs treated with either BmA or ES-62. Moreover, the HIV-1 susceptibility of the Th-cell populations induced by BmA or ES-62 exposed DCs was unaffected for both CCR5- and CXCR4-using HIV-1 viruses. In conclusion, although BmA has the potential capacity to interfere with HIV-1 transmission or initial viral dissemination through preventing the virus from interacting with DCs, no differences in the Th-cell polarizing capacity of DCs exposed to BmA or ES-62 were observed. Neither antigenic source demonstrated beneficial or detrimental effects on the HIV-1 susceptibility of CD4⁺ Th-cells induced by exposed DCs.     

One round of SELEX for the generation of DNA aptamers directed against KLK6

Kallikrein-related peptidase 6 (KLK6) is an active serine protease that has been implicated in common pathologies, including neurodegenerative disorders such as Parkinson and Alzheimer disease and certain types of cancer. Antibodies, either polyclonal or monoclonal, that exhibit specificity for distinct members of the extended kallikrein family, including KLK6, were developed. With the exception of KLK3/PSA, the identification and generation of aptamers, as potential new tools with improved characteristics demanded for therapeutic and diagnostic applications, has not been explored for KLKs. Here, we report for the first time the identification of novel DNA aptamers against KLK6 that were isolated using a modified systemic evolution of ligands by exponential enrichment technique. The identified aptamers were characterized using fluorescence spectroscopy, competition ELISA, and quartz crystal microbalance, and two aptamers (008 and 022) were found to exhibit high affinity (K(d) in the low nanomolar range) for KLK6. Aptamers were tested for their ability to bind to serum albumin, to demonstrate their specificity for their target, and the possible involvement of such proteins in the transport of aptamers into the bloodstream. The developed aptamers are expected to assist the development of novel diagnostic, biosensing, and therapeutic strategies.