Estimating large-scale signaling networks through nested effect models with intervention effects from microarray data

Motivation: Targeted interventions using RNA interference in combination with the measurement of secondary effects with DNA microarrays can be used to computationally reverse engineer features of upstream non-transcriptional signaling cascades based on the nested structure of effects

Nuclear LEF1/TCF4 correlate with poor prognosis but not with nuclear β-catenin in cerebral metastasis of lung adenocarcinomas

An essential function of the transcription factors LEF1/TCF4 in cerebral metastases of lung adenocarcinomas has been described in mouse models, suggesting a WNT/β-catenin effect as potential mechanism. Their role in humans is still unclear, thus we analyzed LEF1, TCF4, β-catenin, and early stage prognostic markers in 25 adenocarcinoma brain metastases using immunohistochemistry (IHC). IHC revealed nuclear TCF4 in all adenocarcinoma samples, whereas only 36 % depicted nuclear LEF1 and nuclear β-catenin signals. Samples with nuclear LEF1 as well as high TCF4 (++++) expression were associated with a shorter survival (p = 0.01, HR = 6.68), while nuclear β-catenin had no significant impact on prognosis and did not significantly correlate with nuclear LEF1. High proliferation index Ki67 was associated with shorter survival in late-stage disease (p = 0.03, HR 3.27). Additionally, we generated a LEF1/TCF4 as well as an AXIN2 signature, the latter as representative of WNT/β-catenin activity, following a bioinformatics approach with a gene expression dataset of cerebral metastases in lung adenocarcinoma. To analyze the prognostic relevance in primary lung adenocarcinomas, we applied both signatures to a microarray dataset of 58 primary lung adenocarcinomas. Only the LEF1/TCF4 signature was able to separate clusters with impact on survival (p = 0.01, HR = 0.32). These clusters displayed diverging enrichment patterns of the cell cycle pathway. In conclusion, our data show that LEF1/TCF4, but not β-catenin, have prognostic relevance in primary and cerebrally metastasized human lung adenocarcinomas. In contrast to the previous in vivo findings, these results indicate that LEF1/TCF4 act independently of β-catenin in this setting

A new analysis approach of epidermal growth factor receptor pathway activation patterns provides insights into cetuximab resistance mechanisms in head and neck cancer

The pathways downstream of the epidermal growth factor receptor (EGFR) have often been implicated to play crucial roles in the development and progression of various cancer types. Different authors have proposed models in cell lines in which they study the modes of pathway activities after perturbation experiments. It is prudent to believe that a better understanding of these pathway activation patterns might lead to novel treatment concepts for cancer patients or at least allow a better stratification of patient collectives into different risk groups or into groups that might respond to different treatments. Traditionally, such analyses focused on the individual players of the pathways. More recently in the field of systems biology, a plethora of approaches that take a more holistic view on the signaling pathways and their downstream transcriptional targets has been developed. Fertig et al. have recently developed a new method to identify patterns and biological process activity from transcriptomics data, and they demonstrate the utility of this methodology to analyze gene expression activity downstream of the EGFR in head and neck squamous cell carcinoma to study cetuximab resistance. Please see related article: http://www.biomedcentral.com/1471-2164/13/16

The histone H2B monoubiquitination regulatory pathway is required for differentiation of multipotent stem cells.

Extensive changes in posttranslational histone modifications accompany the rewiring of the transcriptional program during stem cell differentiation. However, the mechanisms controlling the changes in specific chromatin modifications and their function during differentiation remain only poorly understood. We show that histone H2B monoubiquitination (H2Bub1) significantly increases during differentiation of human mesenchymal stem cells (hMSCs) and various lineage-committed precursor cells and in diverse organisms. Furthermore, the H2B ubiquitin ligase RNF40 is required for the induction of differentiation markers and transcriptional reprogramming of hMSCs. This function is dependent upon CDK9 and the WAC adaptor protein, which are required for H2B monoubiquitination. Finally, we show that RNF40 is required for the resolution of the H3K4me3/H3K27me3 bivalent poised state on lineage-specific genes during the transition from an inactive to an active chromatin conformation. Thus, these data indicate that H2Bub1 is required for maintaining multipotency of hMSCs and plays a central role in controlling stem cell differentiation

Diversity of Nontypeable Haemophilus influenzae Strains Colonizing Australian Aboriginal and Non-Aboriginal Children

Nontypeable Haemophilus influenzae (NTHI) strains are responsible for respiratory-related infections which cause a significant burden of disease in Australian children. We previously identified a disparity in NTHI culture-defined carriage rates between Aboriginal and non-Aboriginal children (42% versus 11%). The aim of this study was to use molecular techniques to accurately determine the true NTHI carriage rates (excluding other culture-identical Haemophilus spp.) and assess whether the NTHI strain diversity correlates with the disparity in NTHI carriage rates. NTHI isolates were cultured from 595 nasopharyngeal aspirates collected longitudinally from asymptomatic Aboriginal (n = 81) and non-Aboriginal (n = 76) children aged 0 to 2 years living in the Kalgoorlie-Boulder region, Western Australia. NTHI-specific 16S rRNA gene PCR and PCR ribotyping were conducted on these isolates. Confirmation of NTHI by 16S rRNA gene PCR corrected the NTHI carriage rates from 42% to 36% in Aboriginal children and from 11% to 9% in non-Aboriginal children. A total of 75 different NTHI ribotypes were identified, with 51% unique to Aboriginal children and 13% unique to non-Aboriginal children (P < 0.0001). The strain richness (proportion of different NTHI ribotypes) was similar for Aboriginal (19%, 65/346) and non-Aboriginal children (19%, 37/192) (P = 0.909). Persistent carriage of the same ribotype was rare in the two groups, but colonization with multiple NTHI strains was more common in Aboriginal children than in non-Aboriginal children. True NTHI carriage was less than that estimated by culture. The Aboriginal children were more likely to carry unique and multiple NTHI strains, which may contribute to the chronicity of NTHI colonization and subsequent diseas

Patterns of subnet usage reveal distinct scales of regulation in the transcriptional regulatory network of Escherichia coli

The set of regulatory interactions between genes, mediated by transcription factors, forms a species' transcriptional regulatory network (TRN). By comparing this network with measured gene expression data one can identify functional properties of the TRN and gain general insight into transcriptional control. We define the subnet of a node as the subgraph consisting of all nodes topologically downstream of the node, including itself. Using a large set of microarray expression data of the bacterium Escherichia coli, we find that the gene expression in different subnets exhibits a structured pattern in response to environmental changes and genotypic mutation. Subnets with less changes in their expression pattern have a higher fraction of feed-forward loop motifs and a lower fraction of small RNA targets within them. Our study implies that the TRN consists of several scales of regulatory organization: 1) subnets with more varying gene expression controlled by both transcription factors and post-transcriptional RNA regulation, and 2) subnets with less varying gene expression having more feed-forward loops and less post-transcriptional RNA regulation.Comment: 14 pages, 8 figures, to be published in PLoS Computational Biolog

Induction of microRNAs, mir-155, mir-222, mir-424 and mir-503, promotes monocytic differentiation through combinatorial regulation

Acute myeloid leukemia (AML) involves a block in terminal differentiation of the myeloid lineage and uncontrolled proliferation of a progenitor state. Using phorbol myristate acetate (PMA), it is possible to overcome this block in THP-1 cells (an M5-AML containing the MLL-MLLT3 fusion), resulting in differentiation to an adherent monocytic phenotype. As part of FANTOM4, we used microarrays to identify 23 microRNAs that are regulated by PMA. We identify four PMA-induced micro- RNAs (mir-155, mir-222, mir-424 and mir-503) that when overexpressed cause cell-cycle arrest and partial differentiation and when used in combination induce additional changes not seen by any individual microRNA. We further characterize these prodifferentiative microRNAs and show that mir-155 and mir-222 induce G2 arrest and apoptosis, respectively. We find mir-424 and mir-503 are derived from a polycistronic precursor mir-424-503 that is under repression by the MLL-MLLT3 leukemogenic fusion. Both of these microRNAs directly target cell-cycle regulators and induce G1 cell-cycle arrest when overexpressed in THP-1. We also find that the pro-differentiative mir-424 and mir-503 downregulate the anti-differentiative mir-9 by targeting a site in its primary transcript. Our study highlights the combinatorial effects of multiple microRNAs within cellular systems.Comment: 45 pages 5 figure

GOSim – an R-package for computation of information theoretic GO similarities between terms and gene products

<p>Abstract</p> <p>Background</p> <p>With the increased availability of high throughput data, such as DNA microarray data, researchers are capable of producing large amounts of biological data. During the analysis of such data often there is the need to further explore the similarity of genes not only with respect to their expression, but also with respect to their functional annotation which can be obtained from Gene Ontology (GO).</p> <p>Results</p> <p>We present the freely available software package <it>GOSim</it>, which allows to calculate the functional similarity of genes based on various information theoretic similarity concepts for GO terms. <it>GOSim </it>extends existing tools by providing additional lately developed functional similarity measures for genes. These can e.g. be used to cluster genes according to their biological function. Vice versa, they can also be used to evaluate the homogeneity of a given grouping of genes with respect to their GO annotation. <it>GOSim </it>hence provides the researcher with a flexible and powerful tool to combine knowledge stored in GO with experimental data. It can be seen as complementary to other tools that, for instance, search for significantly overrepresented GO terms within a given group of genes.</p> <p>Conclusion</p> <p><it>GOSim </it>is implemented as a package for the statistical computing environment <it>R </it>and is distributed under GPL within the CRAN project.</p

Accurate reconstruction of insertion-deletion histories by statistical phylogenetics

The Multiple Sequence Alignment (MSA) is a computational abstraction that represents a partial summary either of indel history, or of structural similarity. Taking the former view (indel history), it is possible to use formal automata theory to generalize the phylogenetic likelihood framework for finite substitution models (Dayhoff's probability matrices and Felsenstein's pruning algorithm) to arbitrary-length sequences. In this paper, we report results of a simulation-based benchmark of several methods for reconstruction of indel history. The methods tested include a relatively new algorithm for statistical marginalization of MSAs that sums over a stochastically-sampled ensemble of the most probable evolutionary histories. For mammalian evolutionary parameters on several different trees, the single most likely history sampled by our algorithm appears less biased than histories reconstructed by other MSA methods. The algorithm can also be used for alignment-free inference, where the MSA is explicitly summed out of the analysis. As an illustration of our method, we discuss reconstruction of the evolutionary histories of human protein-coding genes.Comment: 28 pages, 15 figures. arXiv admin note: text overlap with arXiv:1103.434

Cross-platform expression profiling demonstrates that SV40 small tumor antigen activates Notch, Hedgehog, and Wnt signaling in human cells

BACKGROUND: We previously analyzed human embryonic kidney (HEK) cell lines for the effects that simian virus 40 (SV40) small tumor antigen (ST) has on gene expression using Affymetrix U133 GeneChips. To cross-validate and extend our initial findings, we sought to compare the expression profiles of these cell lines using an alternative microarray platform. METHODS: We have analyzed matched cell lines with and without expression of SV40 ST using an Applied Biosystems (AB) microarray platform that uses single 60-mer oligonucleotides and single-color quantitative chemiluminescence for detection. RESULTS: While we were able to previously identify only 456 genes affected by ST with the Affymetrix platform, we identified 1927 individual genes with the AB platform. Additional technical replicates increased the number of identified genes to 3478 genes and confirmed the changes in 278 (61%) of our original set of 456 genes. Among the 3200 genes newly identified as affected by SV40 ST, we confirmed 20 by QRTPCR including several components of the Wnt, Notch, and Hedgehog signaling pathways, consistent with SV40 ST activation of these developmental pathways. While inhibitors of Notch activation had no effect on cell survival, cyclopamine had a potent killing effect on cells expressing SV40 ST. CONCLUSIONS: These data show that SV40 ST expression alters cell survival pathways to sensitize cells to the killing effect of Hedgehog pathway inhibitors