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.     

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.     

Analysis of xyloglucan endotransglycosylase/hydrolase (XTH) genes from allotetraploid (Gossypium hirsutum) cotton and its diploid progenitors expressed during fiber elongation

Multiple cellular pathways have been shown to be involved during fiber initiation and elongation stages in the cultivated allotetraploid cotton (Gossypium hirsutum). The cell wall enzymes xyloglucan endotransglycosylase/hydrolases (XTH) have been reported to be associated with the biosynthesis of the cell wall and the growth of cotton fibers, probably regulating the plasticity of the primary cell wall. Among various cotton fiber cDNAs found to be preferentially expressed in cotton fibers, a xyloglucan endotransglycosylase (XTH) cDNA was significantly up-regulated during the elongation stage of cotton fiber development. In the present study, we isolated and characterized genomic clones encoding cotton XTH from cultivated cotton (Gossypium hirsutum) and its diploid progenitors (Gossypium arboreum and Gossypium raimondii), designated GhXTH1-1, GhXTH1-2, GaXTH1 and GrXTH, respectively. In addition, we isolated and characterized, by in silico methods, the putative promoter of XTH1 from Gossypium hirsutum. Sequence analysis revealed more than 50% homology to XTH's at the protein level. DNA gel blot hybridization indicated that at least two copies of GhXTH1 are present in Gossypium hirsutum whereas the diploid progenitor species Gossypium arboreum and Gossypium raimondii has only a single copy. Quantitative real-time PCR and high-resolution melting experiments indicated that in Gossypium hirsutum cultivars, in cotton fibers during early stages of fiber elongation specifically expressing only the GhXTH1-1 gene and expression levels of GhXTH1-1 in fibers varies among cultivars differing in fiber percentage and fiber length.     

Immune Response Gene Expression in Colorectal Cancer Carries Distinct Prognostic Implications According to Tissue,Stage and Site: A Prospective Retrospective Translational Study in the Context of a Hellenic Cooperative Oncology Group Randomised Trial

Background

Although host immune response is an emerging prognostic factor for colorectal cancer, there is no consensus on the optimal methodology, surrogate markers or tissue for study.

Patients and Methods

Tumour blocks were prospectively collected from 344 patients with stage II/III colorectal cancer (CRC) treated with adjuvant chemotherapy. Whole section lymphocytic infiltration was studied along with mRNA expression of CD3Z, CD8, CD4, CXCL9, CXCL13, IGHM, FOXP3, SNAI2 and ESR1 by qRT-qPCR in tissue microarray (TMA) cores from the centre of tumour, invasive margin and adjacent normal mucosa.

Results

Lymphocytic infiltration, deficient MMR (10.9%), KRAS (40.7%) and BRAF (4.9%) mutations or single mRNA gene expression were not prognostic. Tumour ESR1 gene expression (Hazard Ratio [HR] for relapse 2.33, 95% CI 1.35-4.02; HR for death 1.74, 95% CI 1.02-2.97) and absence of necrosis (HR for relapse 1.71, 95% CI 1.05-2.71; HR for death 1.98, 95% CI 1.14-3.43) were adverse prognostic features. We used CD3Z and CD8 expression in order to devise the mRNA-based Immune Score (mIS) and proceeded to partitioning analysis in 267 patients, with age, stage, tumour site (Right vs Left CRC), KRAS mutation and tumour mIS as input factors. Only in patients with stage III right-sided colon cancer, a low immune response was associated with inferior disease-free survival (mIS-low, HR for relapse 2.28, 95% CI 1.05-8.02). No prognostic significance was seen for tumour mIS in any other stage or site of CRC, or for a similar mIS score derived from adjacent normal mucosa. Independent adverse prognostic significance was retained in multivariable analysis for absence of necrosis, tumour ESR1 expression in all patients and low tumour mIS in stage III right-sided CRC.

Conclusions

In localised CRC, mRNA-based CD3Z/CD8 profiling of tumour immune response may have stage, site and tissue-specific prognostic significance, along with ESR1 expression.

Trial Registration

ANZCTR.org.au ACTRN12610000509066     

Gene expression and biochemical characterization of Azotobacter vinelandii cyclophilins and Protein Interaction Studies of the cytoplasmic isoform with dnaK and lpxH

The soil nitrogen-fixing bacterium Azotobacter vinelandii possesses two cyclophilins, comprising putative cytoplasmic and periplasmic isoforms, designated as AvPPIB and AvPPIA, respectively. Both recombinant cyclophilins have been purified and their peptidyl-prolyl cis/trans isomerase activity against Suc-Ala-Xaa-Pro-Phe-pNA synthetic peptides has been characterized. The substrate specificity of both cyclophilins is typical for bacterial cyclophilins, with Suc-Ala-Ala-Pro-Phe-pNA being the most rapidly catalyzed substrate. The cytoplasmic cyclophilin also displays a chaperone function in the citrate synthase thermal aggregation assay. Using real-time quantitative RT-PCR, we demonstrate that AvppiB is expressed under various physiological and growth conditions, mainly upregulated by acetate and downregulated by the stationary growth state, while AvppiA shows a tendency for downregulation under the tested conditions. Further, we identified chaperone protein dnaK and UDP-2, 3-diacylglucosamine hydrolase lpxH as probable interacting partners of AvPPIB and we demonstrate their physical interaction by coexpression studies. An increase in AvPPIB PPIase activity in the presence of AvdnaK and a decrease in the presence of AvlpxH further confirms each interaction. However, the PPIase activity does not seem to be essential for those interactions since AvPPIB active site mutants still interact with dnaK and lpxH, while their minor PPIase activity cannot be modulated by the interaction.     

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.     

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.