Detecting separate time scales in genetic expression data.

BACKGROUND: Biological processes occur on a vast range of time scales, and many of them occur concurrently. As a result, system-wide measurements of gene expression have the potential to capture many of these processes simultaneously. The challenge however, is to separate these processes and time scales in the data. In many cases the number of processes and their time scales is unknown. This issue is particularly relevant to developmental biologists, who are interested in processes such as growth, segmentation and differentiation, which can all take place simultaneously, but on different time scales. RESULTS: We introduce a flexible and statistically rigorous method for detecting different time scales in time-series gene expression data, by identifying expression patterns that are temporally shifted between replicate datasets. We apply our approach to a Saccharomyces cerevisiae cell-cycle dataset and an Arabidopsis thaliana root developmental dataset. In both datasets our method successfully detects processes operating on several different time scales. Furthermore we show that many of these time scales can be associated with particular biological functions. CONCLUSIONS: The spatiotemporal modules identified by our method suggest the presence of multiple biological processes, acting at distinct time scales in both the Arabidopsis root and yeast. Using similar large-scale expression datasets, the identification of biological processes acting at multiple time scales in many organisms is now possible.

Yeast screens identify the RNA polymerase II CTD and SPT5 as relevant targets of BRCA1 interaction.

BRCA1 has been implicated in numerous DNA repair pathways that maintain genome integrity, however the function responsible for its tumor suppressor activity in breast cancer remains obscure. To identify the most highly conserved of the many BRCA1 functions, we screened the evolutionarily distant eukaryote Saccharomyces cerevisiae for mutants that suppressed the G1 checkpoint arrest and lethality induced following heterologous BRCA1 expression. A genome-wide screen in the diploid deletion collection combined with a screen of ionizing radiation sensitive gene deletions identified mutants that permit growth in the presence of BRCA1. These genes delineate a metabolic mRNA pathway that temporally links transcription elongation (SPT4, SPT5, CTK1, DEF1) to nucleopore-mediated mRNA export (ASM4, MLP1, MLP2, NUP2, NUP53, NUP120, NUP133, NUP170, NUP188, POM34) and cytoplasmic mRNA decay at P-bodies (CCR4, DHH1). Strikingly, BRCA1 interacted with the phosphorylated RNA polymerase II (RNAPII) carboxy terminal domain (P-CTD), phosphorylated in the pattern specified by the CTDK-I kinase, to induce DEF1-dependent cleavage and accumulation of a RNAPII fragment containing the P-CTD. Significantly, breast cancer associated BRCT domain defects in BRCA1 that suppressed P-CTD cleavage and lethality in yeast also suppressed the physical interaction of BRCA1 with human SPT5 in breast epithelial cells, thus confirming SPT5 as a relevant target of BRCA1 interaction. Furthermore, enhanced P-CTD cleavage was observed in both yeast and human breast cells following UV-irradiation indicating a conserved eukaryotic damage response. Moreover, P-CTD cleavage in breast epithelial cells was BRCA1-dependent since damage-induced P-CTD cleavage was only observed in the mutant BRCA1 cell line HCC1937 following ectopic expression of wild type BRCA1. Finally, BRCA1, SPT5 and hyperphosphorylated RPB1 form a complex that was rapidly degraded following MMS treatment in wild type but not BRCA1 mutant breast cells. These results extend the mechanistic links between BRCA1 and transcriptional consequences in response to DNA damage and suggest an important role for RNAPII P-CTD cleavage in BRCA1-mediated cancer suppression.

Comparison of pattern detection methods in microarray time series of the segmentation clock.

While genome-wide gene expression data are generated at an increasing rate, the repertoire of approaches for pattern discovery in these data is still limited. Identifying subtle patterns of interest in large amounts of data (tens of thousands of profiles) associated with a certain level of noise remains a challenge. A microarray time series was recently generated to study the transcriptional program of the mouse segmentation clock, a biological oscillator associated with the periodic formation of the segments of the body axis. A method related to Fourier analysis, the Lomb-Scargle periodogram, was used to detect periodic profiles in the dataset, leading to the identification of a novel set of cyclic genes associated with the segmentation clock. Here, we applied to the same microarray time series dataset four distinct mathematical methods to identify significant patterns in gene expression profiles. These methods are called: Phase consistency, Address reduction, Cyclohedron test and Stable persistence, and are based on different conceptual frameworks that are either hypothesis- or data-driven. Some of the methods, unlike Fourier transforms, are not dependent on the assumption of periodicity of the pattern of interest. Remarkably, these methods identified blindly the expression profiles of known cyclic genes as the most significant patterns in the dataset. Many candidate genes predicted by more than one approach appeared to be true positive cyclic genes and will be of particular interest for future research. In addition, these methods predicted novel candidate cyclic genes that were consistent with previous biological knowledge and experimental validation in mouse embryos. Our results demonstrate the utility of these novel pattern detection strategies, notably for detection of periodic profiles, and suggest that combining several distinct mathematical approaches to analyze microarray datasets is a valuable strategy for identifying genes that exhibit novel, interesting transcriptional patterns.

Multiple endocytic pathways of G protein-coupled receptors delineated by GIT1 sensitivity.

Recently, we identified a GTPase-activating protein for the ADP ribosylation factor family of small GTP-binding proteins that we call GIT1. This protein initially was identified as an interacting partner for the G protein-coupled receptor kinases, and its overexpression was found to affect signaling and internalization of the prototypical beta(2)-adrenergic receptor. Here, we report that GIT1 overexpression regulates internalization of numerous, but not all, G protein-coupled receptors. The specificity of the GIT1 effect is not related to the type of G protein to which a receptor couples, but rather to the endocytic route it uses. GIT1 only affects the function of G protein-coupled receptors that are internalized through the clathrin-coated pit pathway in a beta-arrestin- and dynamin-sensitive manner. Furthermore, the GIT1 effect is not limited to G protein-coupled receptors because overexpression of this protein also affects internalization of the epidermal growth factor receptor. However, constitutive agonist-independent internalization is not regulated by GIT1, because transferrin uptake is not affected by GIT1 overexpression. Thus, GIT1 is a protein involved in regulating the function of signaling receptors internalized through the clathrin pathway and can be used as a diagnostic tool for defining the endocytic pathway of a receptor.

beta2-Adrenergic receptor regulation by GIT1, a G protein-coupled receptor kinase-associated ADP ribosylation factor GTPase-activating protein.

G protein-coupled receptor activation leads to the membrane recruitment and activation of G protein-coupled receptor kinases, which phosphorylate receptors and lead to their inactivation. We have identified a novel G protein-coupled receptor kinase-interacting protein, GIT1, that is a GTPase-activating protein (GAP) for the ADP ribosylation factor (ARF) family of small GTP-binding proteins. Overexpression of GIT1 leads to reduced beta2-adrenergic receptor signaling and increased receptor phosphorylation, which result from reduced receptor internalization and resensitization. These cellular effects of GIT1 require its intact ARF GAP activity and do not reflect regulation of GRK kinase activity. These results suggest an essential role for ARF proteins in regulating beta2-adrenergic receptor endocytosis. Moreover, they provide a mechanism for integration of receptor activation and endocytosis through regulation of ARF protein activation by GRK-mediated recruitment of the GIT1 ARF GAP to the plasma membrane.

Wee1-regulated apoptosis mediated by the crk adaptor protein in Xenopus egg extracts.

Many of the biochemical reactions of apoptotic cell death, including mitochondrial cytochrome c release and caspase activation, can be reconstituted in cell-free extracts derived from Xenopus eggs. In addition, because caspase activation does not occur until the egg extract has been incubated for several hours on the bench, upstream signaling processes occurring before full apoptosis are rendered accessible to biochemical manipulation. We reported previously that the adaptor protein Crk is required for apoptotic signaling in egg extracts (Evans, E.K., W. Lu, S.L. Strum, B.J. Mayer, and S. Kornbluth. 1997. EMBO (Eur. Mol. Biol. Organ.) J. 16:230-241). Moreover, we demonstrated that removal of Crk Src homology (SH)2 or SH3 interactors from the extracts prevented apoptosis. We now report the finding that the relevant Crk SH2-interacting protein, important for apoptotic signaling in the extract, is the well-known cell cycle regulator, Wee1. We have demonstrated a specific interaction between tyrosine-phosphorylated Wee1 and the Crk SH2 domain and have shown that recombinant Wee1 can restore apoptosis to an extract depleted of SH2 interactors. Moreover, exogenous Wee1 accelerated apoptosis in egg extracts, and this acceleration was largely dependent on the presence of endogenous Crk protein. As other Cdk inhibitors, such as roscovitine and Myt1, did not act like Wee1 to accelerate apoptosis, we propose that Wee1-Crk complexes signal in a novel apoptotic pathway, which may be unrelated to Wee1's role as a cell cycle regulator.

Suppression of DNA-damage checkpoint signaling by Rsk-mediated phosphorylation of Mre11.

Ataxia telangiectasia mutant (ATM) is an S/T-Q-directed kinase that is critical for the cellular response to double-stranded breaks (DSBs) in DNA. Following DNA damage, ATM is activated and recruited by the MRN protein complex [meiotic recombination 11 (Mre11)/DNA repair protein Rad50/Nijmegen breakage syndrome 1 proteins] to sites of DNA damage where ATM phosphorylates multiple substrates to trigger cell-cycle arrest. In cancer cells, this regulation may be faulty, and cell division may proceed even in the presence of damaged DNA. We show here that the ribosomal s6 kinase (Rsk), often elevated in cancers, can suppress DSB-induced ATM activation in both Xenopus egg extracts and human tumor cell lines. In analyzing each step in ATM activation, we have found that Rsk targets loading of MRN complex components onto DNA at DSB sites. Rsk can phosphorylate the Mre11 protein directly at S676 both in vitro and in intact cells and thereby can inhibit the binding of Mre11 to DNA with DSBs. Accordingly, mutation of S676 to Ala can reverse inhibition of the response to DSBs by Rsk. Collectively, these data point to Mre11 as an important locus of Rsk-mediated checkpoint inhibition acting upstream of ATM activation.

Ubiquitylation of p53 by the APC/C inhibitor Trim39.

Tripartite motif 39 (Trim39) is a RING domain-containing E3 ubiquitin ligase able to inhibit the anaphase-promoting complex (APC/C) directly. Through analysis of Trim39 function in p53-positive and p53-negative cells, we have found, surprisingly, that p53-positive cells lacking Trim39 could not traverse the G1/S transition. This effect did not result from disinhibition of the APC/C. Moreover, although Trim39 loss inhibited etoposide-induced apoptosis in p53-negative cells, apoptosis was enhanced by Trim39 knockdown in p53-positive cells. Furthermore, we show here that the Trim39 can directly bind and ubiquitylate p53 in vitro and in vivo, leading to p53 degradation. Depletion of Trim39 significantly increased p53 protein levels and cell growth retardation in multiple cell lines. We found that the relative importance of Trim39 and the well-characterized p53-directed E3 ligase, murine double minute 2 (MDM2), varied between cell types. In cells that were relatively insensitive to the MDM2 inhibitor, nutlin-3a, apoptosis could be markedly enhanced by siRNA directed against Trim39. As such, Trim39 may serve as a potential therapeutic target in tumors with WT p53 when MDM2 inhibition is insufficient to elevate p53 levels and apoptosis.

Regulation of DNA Double Strand Break Response

To ensure genomic integrity, dividing cells implement multiple checkpoint pathways during the course of the cell cycle. In response to DNA damage, cells may either halt the progression of the cycle (cell cycle arrest) or undergo apoptosis. This choice depends on the extent of damage and the cell's capacity for DNA repair. Cell cycle arrest induced by double-stranded DNA breaks relies on the activation of the ataxia-telangiectasia (ATM) protein kinase, which phosphorylates cell cycle effectors (e.g., Chk2 and p53) to inhibit cell cycle progression. ATM is an S/T-Q directed kinase that is critical for the cellular response to double-stranded DNA breaks. Following DNA damage, ATM is activated and recruited to sites of DNA damage by the MRN protein complex (Mre11-Rad50-Nbs1 proteins) where ATM phosphorylates multiple substrates to trigger a cell cycle arrest. In cancer cells, this regulation may be faulty and cell division may proceed even in the presence of damaged DNA. We show here that the RSK kinase, often elevated in cancers, can suppress DSB-induced ATM activation in both Xenopus egg extracts and human tumor cell lines. In analyzing each step in ATM activation, we have found that RSK disrupts the binding of the MRN complex to DSB DNA. RSK can directly phosphorylate the Mre11 protein at Ser 676 both in vitro and in intact cells and can thereby inhibit loading of Mre11 onto DSB DNA. Accordingly, mutation of Ser 676 to Ala can reverse inhibition of the DSB response by RSK. Collectively, these data point to Mre11 as an important locus of RSK-mediated checkpoint inhibition acting upstream of ATM activation.

The phosphorylation of Mre11 on Ser 676 is antagonized by phosphatases. Here, we screened for phosphatases that target this site and identified PP5 as a candidate. This finding is consistent with the fact that PP5 is required for the ATM-mediated DNA damage response, indicating that PP5 may promote DSB-induced, ATM-dependent DNA damage response by targeting Mre11 upstream of ATM.

Measuring fast gene dynamics in single cells with time-lapse luminescence microscopy.

Time-lapse fluorescence microscopy is an important tool for measuring in vivo gene dynamics in single cells. However, fluorescent proteins are limited by slow chromophore maturation times and the cellular autofluorescence or phototoxicity that arises from light excitation. An alternative is luciferase, an enzyme that emits photons and is active upon folding. The photon flux per luciferase is significantly lower than that for fluorescent proteins. Thus time-lapse luminescence microscopy has been successfully used to track gene dynamics only in larger organisms and for slower processes, for which more total photons can be collected in one exposure. Here we tested green, yellow, and red beetle luciferases and optimized substrate conditions for in vivo luminescence. By combining time-lapse luminescence microscopy with a microfluidic device, we tracked the dynamics of cell cycle genes in single yeast with subminute exposure times over many generations. Our method was faster and in cells with much smaller volumes than previous work. Fluorescence of an optimized reporter (Venus) lagged luminescence by 15-20 min, which is consistent with its known rate of chromophore maturation in yeast. Our work demonstrates that luciferases are better than fluorescent proteins at faithfully tracking the underlying gene expression.

Nucleolar organization, ribosomal DNA array stability, and acrocentric chromosome integrity are linked to telomere function.

The short arms of the ten acrocentric human chromosomes share several repetitive DNAs, including ribosomal RNA genes (rDNA). The rDNA arrays correspond to nucleolar organizing regions that coalesce each cell cycle to form the nucleolus. Telomere disruption by expressing a mutant version of telomere binding protein TRF2 (dnTRF2) causes non-random acrocentric fusions, as well as large-scale nucleolar defects. The mechanisms responsible for acrocentric chromosome sensitivity to dysfunctional telomeres are unclear. In this study, we show that TRF2 normally associates with the nucleolus and rDNA. However, when telomeres are crippled by dnTRF2 or RNAi knockdown of TRF2, gross nucleolar and chromosomal changes occur. We used the controllable dnTRF2 system to precisely dissect the timing and progression of nucleolar and chromosomal instability induced by telomere dysfunction, demonstrating that nucleolar changes precede the DNA damage and morphological changes that occur at acrocentric short arms. The rDNA repeat arrays on the short arms decondense, and are coated by RNA polymerase I transcription binding factor UBF, physically linking acrocentrics to one another as they become fusogenic. These results highlight the importance of telomere function in nucleolar stability and structural integrity of acrocentric chromosomes, particularly the rDNA arrays. Telomeric stress is widely accepted to cause DNA damage at chromosome ends, but our findings suggest that it also disrupts chromosome structure beyond the telomere region, specifically within the rDNA arrays located on acrocentric chromosomes. These results have relevance for Robertsonian translocation formation in humans and mechanisms by which acrocentric-acrocentric fusions are promoted by DNA damage and repair.

Integrated pathway and epistasis analysis reveals interactive effect of genetic variants at TERF1 and AFAP1L2 loci on melanoma risk.

Genome-wide association studies (GWASs) have characterized 13 loci associated with melanoma, which only account for a small part of melanoma risk. To identify new genes with too small an effect to be detected individually but which collectively influence melanoma risk and/or show interactive effects, we used a two-step analysis strategy including pathway analysis of genome-wide SNP data, in a first step, and epistasis analysis within significant pathways, in a second step. Pathway analysis, using the gene-set enrichment analysis (GSEA) approach and the gene ontology (GO) database, was applied to the outcomes of MELARISK (3,976 subjects) and MDACC (2,827 subjects) GWASs. Cross-gene SNP-SNP interaction analysis within melanoma-associated GOs was performed using the INTERSNP software. Five GO categories were significantly enriched in genes associated with melanoma (false discovery rate ≤ 5% in both studies): response to light stimulus, regulation of mitotic cell cycle, induction of programmed cell death, cytokine activity and oxidative phosphorylation. Epistasis analysis, within each of the five significant GOs, showed significant evidence for interaction for one SNP pair at TERF1 and AFAP1L2 loci (pmeta-int  = 2.0 × 10(-7) , which met both the pathway and overall multiple-testing corrected thresholds that are equal to 9.8 × 10(-7) and 2.0 × 10(-7) , respectively) and suggestive evidence for another pair involving correlated SNPs at the same loci (pmeta-int  = 3.6 × 10(-6) ). This interaction has important biological relevance given the key role of TERF1 in telomere biology and the reported physical interaction between TERF1 and AFAP1L2 proteins. This finding brings a novel piece of evidence for the emerging role of telomere dysfunction into melanoma development.

Consensus nomenclature for the human ArfGAP domain-containing proteins.

At the FASEB summer research conference on "Arf Family GTPases", held in Il Ciocco, Italy in June, 2007, it became evident to researchers that our understanding of the family of Arf GTPase activating proteins (ArfGAPs) has grown exponentially in recent years. A common nomenclature for these genes and proteins will facilitate discovery of biological functions and possible connections to pathogenesis. Nearly 100 researchers were contacted to generate a consensus nomenclature for human ArfGAPs. This article describes the resulting consensus nomenclature and provides a brief description of each of the 10 subfamilies of 31 human genes encoding proteins containing the ArfGAP domain.

The Physarum polycephalum Genome Reveals Extensive Use of Prokaryotic Two-Component and Metazoan-Type Tyrosine Kinase Signaling.

Physarum polycephalum is a well-studied microbial eukaryote with unique experimental attributes relative to other experimental model organisms. It has a sophisticated life cycle with several distinct stages including amoebal, flagellated, and plasmodial cells. It is unusual in switching between open and closed mitosis according to specific life-cycle stages. Here we present the analysis of the genome of this enigmatic and important model organism and compare it with closely related species. The genome is littered with simple and complex repeats and the coding regions are frequently interrupted by introns with a mean size of 100 bases. Complemented with extensive transcriptome data, we define approximately 31,000 gene loci, providing unexpected insights into early eukaryote evolution. We describe extensive use of histidine kinase-based two-component systems and tyrosine kinase signaling, the presence of bacterial and plant type photoreceptors (phytochromes, cryptochrome, and phototropin) and of plant-type pentatricopeptide repeat proteins, as well as metabolic pathways, and a cell cycle control system typically found in more complex eukaryotes. Our analysis characterizes P. polycephalum as a prototypical eukaryote with features attributed to the last common ancestor of Amorphea, that is, the Amoebozoa and Opisthokonts. Specifically, the presence of tyrosine kinases in Acanthamoeba and Physarum as representatives of two distantly related subdivisions of Amoebozoa argues against the later emergence of tyrosine kinase signaling in the opisthokont lineage and also against the acquisition by horizontal gene transfer.

HDAC inhibitors cause site-specific chromatin remodeling at PU.1-bound enhancers in K562 cells

BACKGROUND: Small molecule inhibitors of histone deacetylases (HDACi) hold promise as anticancer agents for particular malignancies. However, clinical use is often confounded by toxicity, perhaps due to indiscriminate hyperacetylation of cellular proteins. Therefore, elucidating the mechanisms by which HDACi trigger differentiation, cell cycle arrest, or apoptosis of cancer cells could inform development of more targeted therapies. We used the myelogenous leukemia line K562 as a model of HDACi-induced differentiation to investigate chromatin accessibility (DNase-seq) and expression (RNA-seq) changes associated with this process. RESULTS: We identified several thousand specific regulatory elements [~10 % of total DNase I-hypersensitive (DHS) sites] that become significantly more or less accessible with sodium butyrate or suberanilohydroxamic acid treatment. Most of the differential DHS sites display hallmarks of enhancers, including being enriched for non-promoter regions, associating with nearby gene expression changes, and increasing luciferase reporter expression in K562 cells. Differential DHS sites were enriched for key hematopoietic lineage transcription factor motifs, including SPI1 (PU.1), a known pioneer factor. We found PU.1 increases binding at opened DHS sites with HDACi treatment by ChIP-seq, but PU.1 knockdown by shRNA fails to block the chromatin accessibility and expression changes. A machine-learning approach indicates H3K27me3 initially marks PU.1-bound sites that open with HDACi treatment, suggesting these sites are epigenetically poised. CONCLUSIONS: We find HDACi treatment of K562 cells results in site-specific chromatin remodeling at epigenetically poised regulatory elements. PU.1 shows evidence of a pioneer role in this process by marking poised enhancers but is not required for transcriptional activation.