Punctuated evolution and transitional hybrid network in an ancestral cell cycle of fungi.

Although cell cycle control is an ancient, conserved, and essential process, some core animal and fungal cell cycle regulators share no more sequence identity than non-homologous proteins. Here, we show that evolution along the fungal lineage was punctuated by the early acquisition and entrainment of the SBF transcription factor through horizontal gene transfer. Cell cycle evolution in the fungal ancestor then proceeded through a hybrid network containing both SBF and its ancestral animal counterpart E2F, which is still maintained in many basal fungi. We hypothesize that a virally-derived SBF may have initially hijacked cell cycle control by activating transcription via the cis-regulatory elements targeted by the ancestral cell cycle regulator E2F, much like extant viral oncogenes. Consistent with this hypothesis, we show that SBF can regulate promoters with E2F binding sites in budding yeast.

Evaluation of predicted network modules in yeast metabolism using NMR-based metabolite profiling

Genome-scale metabolic models promise important insights into cell function. However, the definition of pathways and functional network modules within these models, and in the biochemical literature in general, is often based on intuitive reasoning. Although mathematical methods have been proposed to identify modules, which are defined as groups of reactions with correlated fluxes, there is a need for experimental verification. We show here that multivariate statistical analysis of the NMR-derived intra- and extracellular metabolite profiles of single-gene deletion mutants in specific metabolic pathways in the yeast Saccharomyces cerevisiae identified outliers whose profiles were markedly different from those of the other mutants in their respective pathways. Application of flux coupling analysis to a metabolic model of this yeast showed that the deleted gene in an outlying mutant encoded an enzyme that was not part of the same functional network module as the other enzymes in the pathway. We suggest that metabolomic methods such as this, which do not require any knowledge of how a gene deletion might perturb the metabolic network, provide an empirical method for validating and ultimately refining the predicted network structure.

Information processing in the transcriptional regulatory network of yeast: Functional robustness

Background: Gene networks are considered to represent various aspects of molecular biological systems meaningfully because they naturally provide a systems perspective of molecular interactions. In this respect, the functional understanding of the transcriptional regulatory network is considered as key to elucidate the functional organization of an organism.

Identification of gene networks associated with acute myeloid leukemia by comparative molecular methylation and expression profiling

Around 80% of acute myeloid leukemia (AML) patients achieve a complete remission, however many will relapse and ultimately die of their disease. The association between karyotype and prognosis has been studied extensively and identified patient cohorts as having favourable [e.g. t(8; 21), inv (16)/t(16; 16), t(15; 17)], intermediate [e.g. cytogenetically normal (NK-AML)] or adverse risk [e.g. complex karyotypes]. Previous studies have shown that gene expression profiling signatures can classify the sub-types of AML, although few reports have shown a similar feature by using methylation markers. The global methylation patterns in 19 diagnostic AML samples were investigated using the Methylated CpG Island Amplification Microarray (MCAM) method and CpG island microarrays containing 12,000 CpG sites. The first analysis, comparing favourable and intermediate cytogenetic risk groups, revealed significantly differentially methylated CpG sites (594 CpG islands) between the two subgroups. Mutations in the NPM1 gene occur at a high frequency (40%) within the NK-AML subgroup and are associated with a more favourable prognosis in these patients. A second analysis comparing the NPM1 mutant and wild-type research study subjects again identified distinct methylation profiles between these two subgroups. Network and pathway analysis revealed possible molecular mechanisms associated with the different risk and/or mutation sub-groups. This may result in a better classification of the risk groups, improved monitoring targets, or the identification of novel molecular therapies.

Inferring the conservative causal core of gene regulatory networks

Background
Inferring gene regulatory networks from large-scale expression data is an important problem that received much attention in recent years. These networks have the potential to gain insights into causal molecular interactions of biological processes. Hence, from a methodological point of view, reliable estimation methods based on observational data are needed to approach this problem practically.

Results
In this paper, we introduce a novel gene regulatory network inference (GRNI) algorithm, called C3NET. We compare C3NET with four well known methods, ARACNE, CLR, MRNET and RN, conducting in-depth numerical ensemble simulations and demonstrate also for biological expression data from E. coli that C3NET performs consistently better than the best known GRNI methods in the literature. In addition, it has also a low computational complexity. Since C3NET is based on estimates of mutual information values in conjunction with a maximization step, our numerical investigations demonstrate that our inference algorithm exploits causal structural information in the data efficiently.

Conclusions
For systems biology to succeed in the long run, it is of crucial importance to establish methods that extract large-scale gene networks from high-throughput data that reflect the underlying causal interactions among genes or gene products. Our method can contribute to this endeavor by demonstrating that an inference algorithm with a neat design permits not only a more intuitive and possibly biological interpretation of its working mechanism but can also result in superior results.

Analysis of HSC activity and compensatory Hox gene expression profile in Hoxb cluster mutant fetal liver cells

Overexpression of Hoxb4 in bone marrow cells promotes expansion of hematopoietic stem cell (HSC) populations in vivo and in vitro, indicating that this homeoprotein can activate the genetic program that determines self-renewal. However, this function cannot be solely attributed to Hoxb4 because Hoxb4(-/-) mice are viable and have an apparently normal HSC number. Quantitative polymerase chain reaction analysis showed that Hoxb4(-/-) c-Kit(+) fetal liver cells expressed moderately higher levels of several Hoxb cluster genes than control cells, raising the possibility that normal HSC activity in Hoxb4(-/-) mice is due to a compensatory up-regulation of other Hoxb genes. In this study, we investigated the competitive repopulation potential of HSCs lacking Hoxb4 alone, or in conjunction with 8 other Hoxb genes. Our results show that Hoxb4(-/-) and Hoxb1-b9(-/-) fetal liver cells retain full competitive repopulation potential and the ability to regenerate all myeloid and lymphoid lineages. Quantitative Hox gene expression profiling in purified c-KIt(+) Hoxb1-bg(-/-) fetal liver cells revealed an interaction between the Hoxa, b, and c clusters with variation in expression levels of Hoxa4, -a11, and -c4. Together, these studies show a complex network of genetic interactions between several Hox genes in primitive hematopoietic cells and demonstrate that HSCs lacking up to 30% of the active Hox genes remain fully competent.

Global down-regulation of HOX gene expression in PML-RAR alpha(+) acute promyelocytic leukemia identified by small-array real-time PCR

Acute promyelocytic leukemia (APL) is associated with a reciprocal and balanced translocation involving the retinoic acid receptor-alpha (RARalpha). All-trans retinoic acid (ATRA) is used to treat APL and is a potent morphogen that regulates HOX gene expression in embryogenesis and organogenesis. HOX genes are also involved in hematopoiesis and leukemogenesis. Thirty-nine mammalian HOX genes have been identified and classified into 13 paralogous groups clustered on 4 chromosomes. They encode a complex net-Work of transcription regulatory proteins whose precise targets remain poorly understood. The overall function of the network appears to be dictated by gene dosage. To investigate the mechanisms involved in HOX gene regulation in hematopoiesis and leukemogenesis by precise measurement of individual HOX genes, a small-array real-time HOX (SMART-HOX) quantitative polymerase chain reaction (PCR) platform was designed and validated. Application of SMART-HOX to 16 APL bone marrow samples revealed a global down-regulation of 26 HOX genes compared with normal controls. HOX gene expression was also altered during differentiation induced by ATRA in the PML-RARalpha(+) NB4 cell line. PML-RARalpha, fusion proteins have been reported to act as part of a repressor complex during myelold cell differentiation, and a model linking HOX gene expression to this PML-RARalpha repressor complex is now proposed.

Activated oncogenic pathways and therapeutic targets in extranodal nasal-type NK/T cell lymphoma revealed by gene expression profiling

We performed comprehensive genome-wide gene expression profiling (GEP) of extranodal nasal-type natural killer/T-cell lymphoma (NKTL) using formalin-fixed, paraffin-embedded tissue (n = 9) and NK cell lines (n = 5) in comparison with normal NK cells, with the objective of understanding the oncogenic pathways involved in the pathogenesis of NKTL and to identify potential therapeutic targets. Pathway and network analysis of genes differentially expressed between NKTL and normal NK cells revealed significant enrichment for cell cycle-related genes and pathways, such as PLK1, CDK1, and Aurora-A. Furthermore, our results demonstrated a pro-proliferative and anti-apoptotic phenotype in NKTL characterized by activation of Myc and nuclear factor kappa B (NF-kappa B), and deregulation of p53. In corroboration with GEP findings, a significant percentage of NKTLs (n = 33) overexpressed c-Myc (45.4%), p53 (87.9%), and NF-kappa B p50 (67.7%) on immunohistochemistry using a tissue microarray containing 33 NKTL samples. Notably, overexpression of survivin was observed in 97% of cases. Based on our findings, we propose a model of NKTL pathogenesis where deregulation of p53 together with activation of Myc and NF-kappa B, possibly driven by EBV LMP-1, results in the cumulative up-regulation of survivin. Down-regulation of survivin with Terameprocol (EM-1421, a survivin inhibitor) results in reduced cell viability and increased apoptosis in tumour cells, suggesting that targeting survivin may be a potential novel therapeutic strategy in NKTL. Copyright (C) 2011 Pathological Society of Great Britain and Ireland. Published by John Wiley & Sons, Ltd.

Network biology: a direct approach to study biological function

In this paper we discuss the dualism of gene networks and their role in systems biology. We argue that gene networks ( 1) can serve as a conceptual framework, forming a fundamental level of a phenomenological description, and ( 2) are a means to represent and analyze data. The latter point does not only allow a systems analysis but is even amenable for a direct approach to study biological function. Here we focus on the clarity of our main arguments and conceptual meaning of gene networks, rather than the causal inference of gene networks from data. (C) 2010 John Wiley & Sons, Inc. WIREs Syst Biol Med 2011 3 379-391 DOI: 10.1002/wsbm.134

Rod photoreceptor loss in Rho(-/-) mice reduces retinal hypoxia and hypoxia-regulated gene expression

PURPOSE. This study was conducted to evaluate whether regions of the retinal neuropile become hypoxic during periods of high oxygen consumption and whether depletion of the outer retina reduces hypoxia and related changes in gene expression.

METHODS. Retinas from rhodopsin knockout (Rho(-/-)) mice were evaluated along with those of wild-type (WT) control animals. Retinas were also examined at the end of 12-hour dark or light periods, and a separate group was treated with L-cis-diltiazem at the beginning of a 12-hour dark period. Hypoxia was assessed by deposition of hypoxyprobe (HP) and HP-protein adducts were localized by immunohistochemistry and quantified using ELISA. Also, hypoxia-regulated gene expression and transcriptional activity were assessed alongside vascular density.

RESULTS. Hypoxia was observed in the inner nuclear and ganglion cell layers in WT retina and was significantly reduced in Rho (-/-) mice (P < 0.05). Retinal hypoxia was significantly increased during dark adaptation in WT mice (P < 0.05), whereas no change was observed in Rho(-/-) or with L-cis-diltiazem-treated WT mice. Hypoxia-inducible factor (HIF)-1 alpha DNA-binding and VEGF mRNA expression in Rho(-/-) retina was significantly reduced in unison with outer retinal depletion (P < 0.05). Retina from the Rho(-/-) mice displayed an extensive intraretinal vascular network after 6 months, although there was evidence that capillary density was depleted in comparison with that in WT retinas.

CONCLUSIONS. Relative hypoxia occurs in the inner retina especially during dark adaptation. Photoreceptor loss reduces retinal oxygen usage and hypoxia which corresponds with attenuation of the retinal microvasculature. These studies suggest that in normal physiological conditions and diurnal cycles the adult retina exists in a state of borderline hypoxia, making this tissue particularly susceptible to even subtle reductions in perfusion.

Global down-regulation of HOX gene expression in PML-RARalpha + acute promyelocytic leukemia identified by small-array real-time PCR

Acute promyelocytic leukemia (APL) is associated with a reciprocal and balanced translocation involving the retinoic acid receptor-alpha (RARalpha). All-trans retinoic acid (ATRA) is used to treat APL and is a potent morphogen that regulates HOX gene expression in embryogenesis and organogenesis. HOX genes are also involved in hematopoiesis and leukemogenesis. Thirty-nine mammalian HOX genes have been identified and classified into 13 paralogous groups clustered on 4 chromosomes. They encode a complex network of transcription regulatory proteins whose precise targets remain poorly understood. The overall function of the network appears to be dictated by gene dosage. To investigate the mechanisms involved in HOX gene regulation in hematopoiesis and leukemogenesis by precise measurement of individual HOX genes, a small-array real-time HOX (SMART-HOX) quantitative polymerase chain reaction (PCR) platform was designed and validated. Application of SMART-HOX to 16 APL bone marrow samples revealed a global down-regulation of 26 HOX genes compared with normal controls. HOX gene expression was also altered during differentiation induced by ATRA in the PML-RARalpha(+) NB4 cell line. PML-RARalpha fusion proteins have been reported to act as part of a repressor complex during myeloid cell differentiation, and a model linking HOX gene expression to this PML-RARalpha repressor complex is now proposed.

Regulatory network of lipopolysaccharide O-antigen biosynthesis in Yersinia enterocolitica includes cell envelope-dependent signals

Lipopolysaccharide (LPS) is a glycolipid present in the outer membrane of all Gram-negative bacteria, and it is one of the signature molecules recognized by the receptors of the innate immune system. In addition to its lipid A portion (the endotoxin), its O-chain polysaccharide (the O-antigen) plays a critical role in the bacterium-host interplay and, in a number of bacterial pathogens, it is a virulence factor. We present evidence that, in Yersinia enterocolitica serotype O:8, a complex signalling network regulates O-antigen expression in response to temperature. Northern blotting and reporter fusion analyses indicated that temperature regulates the O-antigen expression at the transcriptional level. Promoter cloning showed that the O-antigen gene cluster contains two transcriptional units under the control of promoters P(wb1) and P(wb2). The activity of both promoters is under temperature regulation and is repressed in bacteria grown at 37 degrees C. We demonstrate that the RosA/RosB efflux pump/potassium antiporter system and Wzz, the O-antigen chain length determinant, are indirectly involved in the regulation mainly affecting the activity of promoter P(wb2). The rosAB transcription, under the control of P(ros), is activated at 37 degrees C, and P(wb2) is repressed through the signals generated by the RosAB system activation, i.e. decreased [K+] and increased [H+]. The wzz transcription is under the control of P(wb2), and we show that, at 37 degrees C, overexpression of Wzz downregulates slightly the P(wb1) and P(wb2) activities and more strongly the P(ros) activity, with the net result that more O-antigen is produced. Finally, we demonstrate that overexpression of Wzz causes membrane stress that activates the CpxAR two-component signal transduction system.

The Interlaboratory Robustness Of Next-Generation Sequencing (IRON) Study Phase II: Deep-Sequencing Analyses Of Hematological Malignancies Performed In 8,867 Cases By An International Network Involving 27 Laboratories

Introduction: Amplicon deep-sequencing using second-generation sequencing technology is an innovative molecular diagnostic technique and enables a highly-sensitive detection of mutations. As an international consortium we had investigated previously the robustness, precision, and reproducibility of 454 amplicon next-generation sequencing (NGS) across 10 laboratories from 8 countries (Leukemia, 2011;25:1840-8).

Aims: In Phase II of the study, we established distinct working groups for various hematological malignancies, i.e. acute myeloid leukemia (AML), acute lymphoblastic leukemia (ALL), chronic lymphocytic leukemia (CLL), chronic myelogenous leukemia (CML), myelodysplastic syndromes (MDS), myeloproliferative neoplasms (MPN), and multiple myeloma. Currently, 27 laboratories from 13 countries are part of this research consortium. In total, 74 gene targets were selected by the working groups and amplicons were developed for a NGS deep-sequencing assay (454 Life Sciences, Branford, CT). A data analysis pipeline was developed to standardize mutation interpretation both for accessing raw data (Roche Amplicon Variant Analyzer, 454 Life Sciences) and variant interpretation (Sequence Pilot, JSI Medical Systems, Kippenheim, Germany).

Results: We will report on the design, standardization, quality control aspects, landscape of mutations, as well as the prognostic and predictive utility of this assay in a cohort of 8,867 cases. Overall, 1,146 primer sequences were designed and tested. In detail, for example in AML, 924 cases had been screened for CEBPA mutations. RUNX1 mutations were analyzed in 1,888 cases applying the deep-sequencing read counts to study the stability of such mutations at relapse and their utility as a biomarker to detect residual disease. Analyses of DNMT3A (n=1,041) were focused to perform landscape investigations and to address the prognostic relevance. Additionally, this working group is focusing on TET2, ASXL1, and TP53 analyses. A novel prognostic model is being developed allowing stratification of AML into prognostic subgroups based on molecular markers only. In ALL, 1,124 pediatric and adult cases have been screened, including 763 assays for TP53 mutations both at diagnosis and relapse of ALL. Pediatric and adult leukemia expert labs developed additional content to study the mutation incidence of other B and T lineage markers such as IKZF1, JAK2, IL7R, PAX5, EP300, LEF1, CRLF2, PHF6, WT1, JAK1, PTEN, AKT1, IL7R, NOTCH1, CREBBP, or FBXW7. Further, the molecular landscape of CLL is changing rapidly. As such, a separate working group focused on analyses including NOTCH1, SF3B1, MYD88, XPO1, FBXW7 and BIRC3. Currently, 922 cases were screened to investigate the range of mutational burden of NOTCH1 mutations for their prognostic relevance. In MDS, RUNX1 mutation analyses were performed in 977 cases. The prognostic relevance of TP53 mutations in MDS was assessed in additional 327 cases, including isolated deletions of chromosome 5q. Next, content was developed targeting genes of the cellular splicing component, e.g. SF3B1, SRSF2, U2AF1, and ZRSR2. In BCR-ABL1-negative MPN, nine genes of interest (JAK2, MPL, TET2, CBL, KRAS, EZH2, IDH1, IDH2, ASXL1) have been analyzed in a cohort of 155 primary myelofibrosis cases searching for novel somatic mutations and addressing their relevance for disease progression and leukemia transformation. Moreover, an assay was developed and applied to CMML cases allowing the simultaneous analysis of 25 leukemia-associated target genes in a single sequencing run using just 20 ng of starting DNA. Finally, nine laboratories are studying CML, applying ultra-deep sequencing of the BCR-ABL1 tyrosine kinase domain. Analyses were performed on 615 cases investigating the dynamics of expansion of mutated clones under various tyrosine kinase inhibitor therapies.

Conclusion: Molecular characterization of hematological malignancies today requires high diagnostic sensitivity and specificity. As part of the IRON-II study, a network of laboratories analyzed a variety of disease entities applying amplicon-based NGS assays. Importantly, the consortium not only standardized assay design for disease-specific panels, but also achieved consensus on a common data analysis pipeline for mutation interpretation. Distinct working groups have been forged to address scientific tasks and in total 8,867 cases had been analyzed thus far.

Influence of the experimental design of gene expression studies on the inference of gene regulatory networks: environmental factors

The inference of gene regulatory networks gained within recent years a considerable interest in the biology and biomedical community. The purpose of this paper is to investigate the influence that environmental conditions can exhibit on the inference performance of network inference algorithms. Specifically, we study five network inference methods, Aracne, BC3NET, CLR, C3NET and MRNET, and compare the results for three different conditions: (I) observational gene expression data: normal environmental condition, (II) interventional gene expression data: growth in rich media, (III) interventional gene expression data: normal environmental condition interrupted by a positive spike-in stimulation. Overall, we find that different statistical inference methods lead to comparable, but condition-specific results. Further, our results suggest that non-steady-state data enhance the inferability of regulatory networks.