Charles Warden

Heterogeneous Genomic Molecular Clocks in Primates

Using data from primates, we show that molecular clocks in sites that have been part of a CpG dinucleotide in recent past (CpG sites) and non-CpG sites are of markedly different nature, reflecting differences in their molecular origins. Notably, single nucleotide substitutions at non-CpG sites show clear generation-time dependency, indicating that most of these substitutions occur by errors during DNA replication. On the other hand, substitutions at CpG sites occur relatively constantly over time, as expected from their primary origin due to methylation. Therefore, molecular clocks are heterogeneous even within a genome. Furthermore, we propose that varying frequencies of CpG dinucleotides in different genomic regions may have contributed significantly to conflicting earlier results on rate constancy of mammalian molecular clock. Our conclusion that different regions of genomes follow different molecular clocks should be considered when inferring divergence times using molecular data and in phylogenetic analysis.

Genome-Wide Mutagenesis Reveals That ORF7 Is a Novel VZV Skin-Tropic Factor

The Varicella Zoster Virus (VZV) is a ubiquitous human alpha-herpesvirus that is the causative agent of chicken pox and shingles. Although an attenuated VZV vaccine (v-Oka) has been widely used in children in the United States, chicken pox outbreaks are still seen, and the shingles vaccine only reduces the risk of shingles by 50%. Therefore, VZV still remains an important public health concern. Knowledge of VZV replication and pathogenesis remains limited due to its highly cell-associated nature in cultured cells, the difficulty of generating recombinant viruses, and VZVs almost exclusive tropism for human cells and tissues. In order to circumvent these hurdles, we cloned the entire VZV (p-Oka) genome into a bacterial artificial chromosome that included a dual-reporter system (GFP and luciferase reporter genes). We used PCR-based mutagenesis and the homologous recombination system in the E. coli to individually delete each of the genomes 70 unique ORFs. The collection of viral mutants obtained was systematically examined both in MeWo cells and in cultured human fetal skin organ samples. We use our genome-wide deletion library to provide novel functional annotations to 51% of the VZV proteome. We found 44 out of 70 VZV ORFs to be essential for viral replication. Among the 26 non-essential ORF deletion mutants, eight have discernable growth defects in MeWo. Interestingly, four ORFs were found to be required for viral replication in skin organ cultures, but not in MeWo cells, suggesting their potential roles as skin tropism factors. One of the genes (ORF7) has never been described as a skin tropic factor. The global profiling of the VZV genome gives further insights into the replication and pathogenesis of this virus, which can lead to improved prevention and therapy of chicken pox and shingles.

Deleterious effects of endocrine disruptors are corrected in the mammalian germline by epigenome reprogramming

BackgroundExposure to environmental endocrine-disrupting chemicals during pregnancy reportedly causes transgenerationally inherited reproductive defects. We hypothesized that to affect the grandchild, endocrine-disrupting chemicals must alter the epigenome of the germ cells of the in utero-exposed G1 male fetus. Additionally, to affect the great-grandchild, the aberration must persist in the germ cells of the unexposed G2 grandchild.ResultsHere, we treat gestating female mice with vinclozolin, bisphenol A, or di-(2-ethylhexyl)phthalate during the time when global de novo DNA methylation and imprint establishment occurs in the germ cells of the G1 male fetus. We map genome-wide features in purified G1 and G2 prospermatogonia, in order to detect immediate and persistent epigenetic aberrations, respectively. We detect changes in transcription and methylation in the G1 germline immediately after endocrine-disrupting chemicals exposure, but changes do not persist into the G2 germline. Additional analysis of genomic imprints shows no persistent aberrations in DNA methylation at the differentially methylated regions of imprinted genes between the G1 and G2 prospermatogonia, or in the allele-specific transcription of imprinted genes between the G2 and G3 soma.ConclusionsOur results suggest that endocrine-disrupting chemicals exert direct epigenetic effects in exposed fetal germ cells, which are corrected by reprogramming events in the next generation. Avoiding transgenerational inheritance of environmentally-caused epigenetic aberrations may have played an evolutionary role in the development of dual waves of global epigenome reprogramming in mammals.

COHCAP: an integrative genomic pipeline for single-nucleotide resolution DNA methylation analysis

COHCAP (City of Hope CpG Island Analysis Pipeline) is an algorithm to analyze single-nucleotide resolution DNA methylation data produced by either an Illumina methylation array or targeted bisulfite sequencing. The goal of the COHCAP algorithm is to identify CpG islands that show a consistent pattern of methylation among CpG sites. COHCAP is currently the only DNA methylation package that provides integration with gene expression data to identify a subset of CpG islands that are most likely to regulate downstream gene expression, and it can generate lists of differentially methylated CpG islands with ∼50% concordance with gene expression from both cell line data and heterogeneous patient data. For example, this article describes known breast cancer biomarkers (such as estrogen receptor) with a negative correlation between DNA methylation and gene expression. COHCAP also provides visualization for quality control metrics, regions of differential methylation and correlation between methylation and gene expression. This software is freely available at https://sourceforge.net/projects/cohcap/.

Herpesvirus BACs: Past, Present, and Future

The herpesviridae are a large family of DNA viruses with large and complicated genomes. Genetic manipulation and the generation of recombinant viruses have been extremely difficult. However, herpesvirus bacterial artificial chromosomes (BACs) that were developed approximately 10 years ago have become useful and powerful genetic tools for generating recombinant viruses to study the biology and pathogenesis of herpesviruses. For example, BAC-directed deletion mutants are commonly used to determine the function and essentiality of viral genes. In this paper, we discuss the creation of herpesvirus BACs, functional analyses of herpesvirus mutants, and future applications for studies of herpesviruses. We describe commonly used methods to create and mutate herpesvirus BACs (such as site-directed mutagenesis and transposon mutagenesis). We also evaluate the potential future uses of viral BACs, including vaccine development and gene therapy.

Altered Expression of Polycomb Group Genes in Glioblastoma Multiforme

The Polycomb group (PcG) proteins play a critical role in histone mediated epigenetics which has been implicated in the malignant evolution of glioblastoma multiforme (GBM). By systematically interrogating The Cancer Genome Atlas (TCGA), we discovered widespread aberrant expression of the PcG members in GBM samples compared to normal brain. The most striking differences were upregulation of EZH2, PHF19, CBX8 and PHC2 and downregulation of CBX7, CBX6, EZH1 and RYBP. Interestingly, changes in EZH2, PHF19, CBX7, CBX6 and EZH1 occurred progressively as astrocytoma grade increased. We validated the aberrant expression of CBX6, CBX7, CBX8 and EZH2 in GBM cell lines by Western blotting and qRT-PCR, and further the aberrant expression of CBX6 in GBM tissue samples by immunohistochemical staining. To determine if there was functional significance to the diminished CBX6 levels in GBM, CBX6 was overexpressed in GBM cells resulting in decreased proliferative capacity. In conclusion, aberrant expression of PcG proteins in GBMs may play a role in the development or maintenance of the malignancy.

Glioma IL13Rα2 is associated with mesenchymal signature gene expression and poor patient prognosis.

A major challenge for successful immunotherapy against glioma is the identification and characterization of validated targets. We have taken a bioinformatics approach towards understanding the biological context of IL-13 receptor α2 (IL13Rα2) expression in brain tumors, and its functional significance for patient survival. Querying multiple gene expression databases, we show that IL13Rα2 expression increases with glioma malignancy grade, and expression for high-grade tumors is bimodal, with approximately 58% of WHO grade IV gliomas over-expressing this receptor. By several measures, IL13Rα2 expression in patient samples and low-passage primary glioma lines most consistently correlates with the expression of signature genes defining mesenchymal subclass tumors and negatively correlates with proneural signature genes as defined by two studies. Positive associations were also noted with proliferative signature genes, whereas no consistent associations were found with either classical or neural signature genes. Probing the potential functional consequences of this mesenchymal association through IPA analysis suggests that IL13Rα2 expression is associated with activation of proinflammatory and immune pathways characteristic of mesenchymal subclass tumors. In addition, survival analyses indicate that IL13Rα2 over-expression is associated with poor patient prognosis, a single gene correlation ranking IL13Rα2 in the top ~1% of total gene expression probes with regard to survival association with WHO IV gliomas. This study better defines the functional consequences of IL13Rα2 expression by demonstrating association with mesenchymal signature gene expression and poor patient prognosis. It thus highlights the utility of IL13Rα2 as a therapeutic target, and helps define patient populations most likely to respond to immunotherapy in present and future clinical trials.

Inhibition of the proliferation of acquired aromatase inhibitor-resistant breast cancer cells by histone deacetylase inhibitor LBH589 (panobinostat)

Aromatase inhibitors (AIs) are important drugs for treating postmenopausal patients with hormone receptor-positive breast cancer. However, acquired resistance to AI therapies is a significant problem. Our study has revealed that the histone deacetylase inhibitor LBH589 treatment abrogated growth of AI-resistant cells in vitro and in vivo, causing cell cycle G2/M arrest and induced apoptosis. LBH589 treatment also reduced the level of NF-κB1 which is overexpressed when AI resistance develops. Analyzing paired tumor specimens from 12 patients, we found that NF-κB1 expression was increased in recurrent AI-resistant tumors as compared to the paired primary tumors before AI treatment. This finding was consistent with up-regulated NF-κB1 expression seen in a collection of well-established AI-resistant cell lines. Furthermore, knockdown of NF-κB1 expression significantly suppressed the proliferation of AI-resistant cells. Treatment of AI-resistant cell lines with LBH589 suppressed NF-κB1 mRNA and protein expression. In addition, LBH589 treatment abrogated growth of AI-resistant tumors in mice, and was associated with significantly decreased levels of NF-κB1 in tumors. In all, our findings strongly support further investigation of LBH589 as a novel therapeutic strategy for patients with AI-resistant breast cancer, in part by suppressing the NF-κB1 pathway.

S100B Promotes Glioma Growth through Chemoattraction of Myeloid-Derived Macrophages

Purpose: S100B is member of a multigenic family of Ca2+-binding proteins, which is overexpressed by gliomas. Recently, we showed that low concentrations of S100B attenuated microglia activation through the induction of Stat3. We hypothesized that overexpression of S100B in gliomas could promote tumor growth by modulating the activity of tumor-associated macrophages (TAM). Experimental Design: We stably transfected GL261 glioma cell lines with constructs that overexpressed (S100Bhigh) or underexpressed (S100Blow) S100B and compared their growth characteristics to intracranial wild-type (S100Bwt) tumors. Results: Downregulation of S100B in gliomas had no impact on cell division in vitro but abrogated tumor growth in vivo. Interestingly, compared to S100Blow tumors, S100Bwt and S100Bhigh intracranial gliomas exhibited higher infiltration of TAMs, stronger inflammatory cytokine expression, and increased vascularity. To identify the potential mechanisms involved, the expression of the S100B receptor, receptor for advanced glycation end products (RAGE), was evaluated in gliomas. Although S100B expression induced RAGE in vivo, RAGE ablation in mice did not significantly inhibit TAM infiltration into gliomas, suggesting that other pathways were involved in this process. To evaluate other mechanisms responsible for TAM chemoattraction, we then examined chemokine pathways and found that C-C motif ligand 2 (CCL2) was upregulated in S100Bhigh tumors. Furthermore, analysis of The Cancer Genome Atlass glioma data bank showed a positive correlation between S100B and CCL2 expression in human proneural and neural glioma subtypes, supporting our finding. Conclusions: These observations suggest that S100B promotes glioma growth by TAM chemoattraction through upregulation of CCL2 and introduces the potential utility of S100B inhibitors for glioma therapy. Clin Cancer Res; 19(14); 3764–75. ©2013 AACR.

RNA-seq Reveals Aurora Kinase Driven-mTOR Pathway Activation in Patients with Sarcomatoid Metastatic Renal Cell Carcinoma

Sarcomatoid metastatic renal cell carcinoma (mRCC) is associated with a poor prognosis, and the biology of the disease has been inadequately characterized. RNA sequencing (RNA-seq) was performed on adjacent benign, clear cell, and sarcomatoid components from clinical specimens with sarcomatoid mRCC. M phase and cell-cycle pathways were enriched in sarcomatoid versus adjacent clear cell components, suggesting greater cell proliferation. The expression of aurora kinase A (AURKA) was increased as part of these pathways, and its increased expression was validated by quantitative PCR (qPCR). Immunohistochemical (IHC) analysis revealed that AURKA levels were increased in sarcomatoid tissue compared with their benign or clear cell parts. The increase in AURKA correlated with increased mTOR pathway activity, as evidenced by increased expression of phosphorylated mTOR (S2448) and ribosomal protein S6K (T389). When AURKA was stably expressed in a RCC cell line (Renca), it resulted in increased expression and activity of mTOR, suggesting that overexpression of AURKA can activate the mTOR pathway. These results warrant the analysis of a larger clinical cohort and suggest that targeting AURKA and/or mTOR in patients with sarcomatoid mRCC should be explored. Implications: Comparative RNA-seq of adjacent sarcomatoid and clear cell histology of RCC indicates a proliferative phenotype and increased AURKA-dependent activation of mTOR signaling in sarcomatoid RCC, which could be targeted by available agents. Mol Cancer Res; 13(1); 130–7. ©2014 AACR.