Two genes on A/J chromosome 18 are associated with susceptibility to Staphylococcus aureus infection by combined microarray and QTL analyses.

Although it has recently been shown that A/J mice are highly susceptible to Staphylococcus aureus sepsis as compared to C57BL/6J, the specific genes responsible for this differential phenotype are unknown. Using chromosome substitution strains (CSS), we found that loci on chromosomes 8, 11, and 18 influence susceptibility to S. aureus sepsis in A/J mice. We then used two candidate gene selection strategies to identify genes on these three chromosomes associated with S. aureus susceptibility, and targeted genes identified by both gene selection strategies. First, we used whole genome transcription profiling to identify 191 (56 on chr. 8, 100 on chr. 11, and 35 on chr. 18) genes on our three chromosomes of interest that are differentially expressed between S. aureus-infected A/J and C57BL/6J. Second, we identified two significant quantitative trait loci (QTL) for survival post-infection on chr. 18 using N(2) backcross mice (F(1) [C18A]xC57BL/6J). Ten genes on chr. 18 (March3, Cep120, Chmp1b, Dcp2, Dtwd2, Isoc1, Lman1, Spire1, Tnfaip8, and Seh1l) mapped to the two significant QTL regions and were also identified by the expression array selection strategy. Using real-time PCR, 6 of these 10 genes (Chmp1b, Dtwd2, Isoc1, Lman1, Tnfaip8, and Seh1l) showed significantly different expression levels between S. aureus-infected A/J and C57BL/6J. For two (Tnfaip8 and Seh1l) of these 6 genes, siRNA-mediated knockdown of gene expression in S. aureus-challenged RAW264.7 macrophages induced significant changes in the cytokine response (IL-1 beta and GM-CSF) compared to negative controls. These cytokine response changes were consistent with those seen in S. aureus-challenged peritoneal macrophages from CSS 18 mice (which contain A/J chromosome 18 but are otherwise C57BL/6J), but not C57BL/6J mice. These findings suggest that two genes, Tnfaip8 and Seh1l, may contribute to susceptibility to S. aureus in A/J mice, and represent promising candidates for human genetic susceptibility studies.

The Trim39 ubiquitin ligase inhibits APC/CCdh1-mediated degradation of the Bax activator MOAP-1.

Proapoptotic Bcl-2 family members, such as Bax, promote release of cytochrome c from mitochondria, leading to caspase activation and cell death. It was previously reported that modulator of apoptosis protein 1 (MOAP-1), an enhancer of Bax activation induced by DNA damage, is stabilized by Trim39, a protein of unknown function. In this paper, we show that MOAP-1 is a novel substrate of the anaphase-promoting complex (APC/C(Cdh1)) ubiquitin ligase. The influence of Trim39 on MOAP-1 levels stems from the ability of Trim39 (a RING domain E3 ligase) to directly inhibit APC/C(Cdh1)-mediated protein ubiquitylation. Accordingly, small interfering ribonucleic acid-mediated knockdown of Cdh1 stabilized MOAP-1, thereby enhancing etoposide-induced Bax activation and apoptosis. These data identify Trim39 as a novel APC/C regulator and provide an unexpected link between the APC/C and apoptotic regulation via MOAP-1.

Tyrosine phosphorylation of G protein alpha subunits by pp60c-src.

A number of lines of evidence suggest that cross-talk exists between the cellular signal transduction pathways involving tyrosine phosphorylation catalyzed by members of the pp60c-src kinase family and those mediated by guanine nucleotide regulatory proteins (G proteins). In this study, we explore the possibility that direct interactions between pp60c-src and G proteins may occur with functional consequences. Preparations of pp60c-src isolated by immunoprecipitation phosphorylate on tyrosine residues the purified G-protein alpha subunits (G alpha) of several heterotrimeric G proteins. Phosphorylation is highly dependent on G-protein conformation, and G alpha(GDP) uncomplexed by beta gamma subunits appears to be the preferred substrate. In functional studies, phosphorylation of stimulatory G alpha (G alpha s) modestly increases the rate of binding of guanosine 5'-[gamma-[35S]thio]triphosphate to Gs as well as the receptor-stimulated steady-state rate of GTP hydrolysis by Gs. Heterotrimeric G proteins may represent a previously unappreciated class of potential substrates for pp60c-src.

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.

Evolution of networks and sequences in eukaryotic cell cycle control.

The molecular networks regulating the G1-S transition in budding yeast and mammals are strikingly similar in network structure. However, many of the individual proteins performing similar network roles appear to have unrelated amino acid sequences, suggesting either extremely rapid sequence evolution, or true polyphyly of proteins carrying out identical network roles. A yeast/mammal comparison suggests that network topology, and its associated dynamic properties, rather than regulatory proteins themselves may be the most important elements conserved through evolution. However, recent deep phylogenetic studies show that fungal and animal lineages are relatively closely related in the opisthokont branch of eukaryotes. The presence in plants of cell cycle regulators such as Rb, E2F and cyclins A and D, that appear lost in yeast, suggests cell cycle control in the last common ancestor of the eukaryotes was implemented with this set of regulatory proteins. Forward genetics in non-opisthokonts, such as plants or their green algal relatives, will provide direct information on cell cycle control in these organisms, and may elucidate the potentially more complex cell cycle control network of the last common eukaryotic ancestor.

Computational Inference of Genome-Wide Protein-DNA Interactions Using High-Throughput Genomic Data

Transcriptional regulation has been studied intensively in recent decades. One important aspect of this regulation is the interaction between regulatory proteins, such as transcription factors (TF) and nucleosomes, and the genome. Different high-throughput techniques have been invented to map these interactions genome-wide, including ChIP-based methods (ChIP-chip, ChIP-seq, etc.), nuclease digestion methods (DNase-seq, MNase-seq, etc.), and others. However, a single experimental technique often only provides partial and noisy information about the whole picture of protein-DNA interactions. Therefore, the overarching goal of this dissertation is to provide computational developments for jointly modeling different experimental datasets to achieve a holistic inference on the protein-DNA interaction landscape.

We first present a computational framework that can incorporate the protein binding information in MNase-seq data into a thermodynamic model of protein-DNA interaction. We use a correlation-based objective function to model the MNase-seq data and a Markov chain Monte Carlo method to maximize the function. Our results show that the inferred protein-DNA interaction landscape is concordant with the MNase-seq data and provides a mechanistic explanation for the experimentally collected MNase-seq fragments. Our framework is flexible and can easily incorporate other data sources. To demonstrate this flexibility, we use prior distributions to integrate experimentally measured protein concentrations.

We also study the ability of DNase-seq data to position nucleosomes. Traditionally, DNase-seq has only been widely used to identify DNase hypersensitive sites, which tend to be open chromatin regulatory regions devoid of nucleosomes. We reveal for the first time that DNase-seq datasets also contain substantial information about nucleosome translational positioning, and that existing DNase-seq data can be used to infer nucleosome positions with high accuracy. We develop a Bayes-factor-based nucleosome scoring method to position nucleosomes using DNase-seq data. Our approach utilizes several effective strategies to extract nucleosome positioning signals from the noisy DNase-seq data, including jointly modeling data points across the nucleosome body and explicitly modeling the quadratic and oscillatory DNase I digestion pattern on nucleosomes. We show that our DNase-seq-based nucleosome map is highly consistent with previous high-resolution maps. We also show that the oscillatory DNase I digestion pattern is useful in revealing the nucleosome rotational context around TF binding sites.

Finally, we present a state-space model (SSM) for jointly modeling different kinds of genomic data to provide an accurate view of the protein-DNA interaction landscape. We also provide an efficient expectation-maximization algorithm to learn model parameters from data. We first show in simulation studies that the SSM can effectively recover underlying true protein binding configurations. We then apply the SSM to model real genomic data (both DNase-seq and MNase-seq data). Through incrementally increasing the types of genomic data in the SSM, we show that different data types can contribute complementary information for the inference of protein binding landscape and that the most accurate inference comes from modeling all available datasets.

This dissertation provides a foundation for future research by taking a step toward the genome-wide inference of protein-DNA interaction landscape through data integration.

A C-terminal motif found in the beta2-adrenergic receptor, P2Y1 receptor and cystic fibrosis transmembrane conductance regulator determines binding to the Na+/H+ exchanger regulatory factor family of PDZ proteins.

The Na+/H+ exchanger regulatory factor (NHERF) binds to the tail of the beta2-adrenergic receptor and plays a role in adrenergic regulation of Na+/H+ exchange. NHERF contains two PDZ domains, the first of which is required for its interaction with the beta2 receptor. Mutagenesis studies of the beta2 receptor tail revealed that the optimal C-terminal motif for binding to the first PDZ domain of NHERF is D-S/T-x-L, a motif distinct from those recognized by other PDZ domains. The first PDZ domain of NHERF-2, a protein that is 52% identical to NHERF and also known as E3KARP, SIP-1, and TKA-1, exhibits binding preferences very similar to those of the first PDZ domain of NHERF. The delineation of the preferred binding motif for the first PDZ domain of the NHERF family of proteins allows for predictions for other proteins that may interact with NHERF or NHERF-2. For example, as would be predicted from the beta2 receptor tail mutagenesis studies, NHERF binds to the tail of the purinergic P2Y1 receptor, a seven-transmembrane receptor with an intracellular C-terminal tail ending in D-T-S-L. NHERF also binds to the tail of the cystic fibrosis transmembrane conductance regulator, which ends in D-T-R-L. Because the preferred binding motif of the first PDZ domain of the NHERF family of proteins is found at the C termini of a variety of intracellular proteins, NHERF and NHERF-2 may be multifunctional adaptor proteins involved in many previously unsuspected aspects of intracellular signaling.

MYC activity mitigates response to rapamycin in prostate cancer through eukaryotic initiation factor 4E-binding protein 1-mediated inhibition of autophagy.

Loss of PTEN and activation of phosphoinositide 3-kinase are commonly observed in advanced prostate cancer. Inhibition of mammalian target of rapamycin (mTOR), a downstream target of phosphoinositide 3-kinase signaling, results in cell cycle arrest and apoptosis in multiple in vitro and in vivo models of prostate cancer. However, single-agent use of mTOR inhibition has limited clinical success, and the identification of molecular events mitigating tumor response to mTOR inhibition remains a critical question. Here, using genetically engineered human prostate epithelial cells (PrEC), we show that MYC, a frequent target of genetic gain in prostate cancers, abrogates sensitivity to rapamycin by decreasing rapamycin-induced cytostasis and autophagy. Analysis of MYC and the mTOR pathway in human prostate tumors and PrEC showed selective increased expression of eukaryotic initiation factor 4E-binding protein 1 (4EBP1) with gain in MYC copy number or forced MYC expression, respectively. We have also found that MYC binds to regulatory regions of the 4EBP1 gene. Suppression of 4EBP1 expression resulted in resensitization of MYC-expressing PrEC to rapamycin and increased autophagy. Taken together, our findings suggest that MYC expression abrogates sensitivity to rapamycin through increased expression of 4EBP1 and reduced autophagy.

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.

A network of substrates of the E3 ubiquitin ligases MDM2 and HUWE1 control apoptosis independently of p53.

In the intrinsic pathway of apoptosis, cell-damaging signals promote the release of cytochrome c from mitochondria, triggering activation of the Apaf-1 and caspase-9 apoptosome. The ubiquitin E3 ligase MDM2 decreases the stability of the proapoptotic factor p53. We show that it also coordinated apoptotic events in a p53-independent manner by ubiquitylating the apoptosome activator CAS and the ubiquitin E3 ligase HUWE1. HUWE1 ubiquitylates the antiapoptotic factor Mcl-1, and we found that HUWE1 also ubiquitylated PP5 (protein phosphatase 5), which indirectly inhibited apoptosome activation. Breast cancers that are positive for the tyrosine receptor kinase HER2 (human epidermal growth factor receptor 2) tend to be highly aggressive. In HER2-positive breast cancer cells treated with the HER2 tyrosine kinase inhibitor lapatinib, MDM2 was degraded and HUWE1 was stabilized. In contrast, in breast cancer cells that acquired resistance to lapatinib, the abundance of MDM2 was not decreased and HUWE1 was degraded, which inhibited apoptosis, regardless of p53 status. MDM2 inhibition overcame lapatinib resistance in cells with either wild-type or mutant p53 and in xenograft models. These findings demonstrate broader, p53-independent roles for MDM2 and HUWE1 in apoptosis and specifically suggest the potential for therapy directed against MDM2 to overcome lapatinib resistance.

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.