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/ .

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 VZV's 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 genome's 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. The Varicella Zoster Virus (VZV) is the causative agent of chicken pox and shingles. The long-term efficacy of the current chickenpox vaccine is yet to be determined, and the current shingles vaccine fails to provide protective immunity for a substantial number of individuals. Shingles can also lead to post-herpetic neuralgia (PHN), a debilitating condition associated with an intractable pain that can linger for life. Therefore, VZV remains an important public health concern. We use growth-rate analysis of our genome-wide deletion library to determine the essentiality of all known VZV genes, including novel annotations for 51% of the VZV proteome. We also discovered a novel skin-tropic factor encoded by ORF7. Overall, our identification of genes essential for VZV replication and pathogenesis will serve as the basis for multiple in-depth genetic studies of VZV, which can lead to improved prevention and therapy of chicken pox and shingles. For example, essential genes may be appealing drug targets and genes whose deletion causes a substantial growth defect may be prospective candidates for novel live attenuated vaccines.

Triple negative breast cancer (TNBC) is aggressive with limited treatment options upon recurrence. Molecular discordance between primary and metastatic TNBC has been observed, but the degree of biological heterogeneity has not been fully explored. Furthermore, genomic evolution through treatment is poorly understood. In this study, we aim to characterize the genomic changes between paired primary and metastatic TNBCs through transcriptomic and genomic profiling, and to identify genomic alterations which may contribute to chemotherapy resistance. Genomic alterations and mRNA expression of 10 paired primary and metastatic TNBCs were determined through targeted sequencing, microarray analysis, and RNA sequencing. Commonly mutated genes, as well as differentially expressed and co-expressed genes were identified. We further explored the clinical relevance of differentially expressed genes between primary and metastatic tumors to patient survival using large public datasets. Through gene expression profiling, we observed a shift in TNBC subtype classifications between primary and metastatic TNBCs. A panel of eight cancer driver genes (CCNE1, TPX2, ELF3, FANCL, JAK2, GSK3B, CEP76, and SYK) were differentially expressed in recurrent TNBCs, and were also overexpressed in TCGA and METABRIC. CCNE1 and TPX2 were co-overexpressed in TNBCs. DNA mutation profiling showed that multiple mutations occurred in genes comprising a number of potentially targetable pathways including PI3K/AKT/mTOR, RAS/MAPK, cell cycle, and growth factor receptor signaling, reaffirming the wide heterogeneity of mechanisms driving TNBC. CCNE1 amplification was associated with poor overall survival in patients with metastatic TNBC. CCNE1 amplification may confer resistance to chemotherapy and is associated with poor overall survival in TNBC.

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. The rate at which mutations accumulate in a genome, referred as a “molecular clock,” is an instrumental tool in molecular evolution and phylogenetics. Different types of mutations occur via distinctive molecular pathways. In particular, while most mutations occur from errors in DNA replication, spontaneous deamination of methylated CpG dinucleotides is another important source of mutation in mammalian genomes. Molecular clock studies typically combined all types of mutations together. In this paper, the authors analyze molecular clocks of replication-origin and methylation-origin mutations separately. By utilizing high-quality sequence data from several primate species and fossil calibration, the authors demonstrate that the two types of mutations follow statistically different molecular clocks. Methylation-origin mutations accumulate relatively constantly over time, while replication-origin mutations scale with generation-times. Therefore, the genomic molecular clock, as a whole, is shaped by the molecular origins of mutations that have accumulated over time. The authors' results have direct implications on phylogenetic analyses, estimation of species divergence dates, and studies of the mechanisms and processes of evolution, where molecular clocks are imperative.

Menopause is a critical window of susceptibility for its sensitivity to endocrine disrupting chemicals due to the decline of endogenous estrogen. Using a surgical menopausal (ovariectomized) mouse model, we assessed how mammary tissue was affected by both 17β-estradiol (E2) and polybrominated diphenyl ethers (PBDEs). As flame retardants in household products, PBDEs are widely detected in human serum. During physiologically-relevant exposure to E2, PBDEs enhanced E2-mediated regrowth of mammary glands with terminal end bud-like structures. Analysis of mammary gland RNA revealed that PBDEs both augmented E2-facilitated gene expression and modulated immune regulation. Through single-cell RNA sequencing (scRNAseq) analysis, E2 was found to induce Pgr expression in both Esr1⁺ and Esr1⁻ luminal epithelial cells and Ccl2 expression in Esr1⁺ fibroblasts. PBDEs promote the E2-AREG-EGFR-M2 macrophage pathway. Our findings support that E2 + PBDE increases the risk of developing breast cancer through the expansion of estrogen-responsive luminal epithelial cells and immune modulation.

Ribonucleotide reductase (RR) is an essential enzyme involved in DNA synthesis. We hypothesized that RR subunit M2 (RRM2) might be a novel prognostic and predictive biomarker for estrogen receptor (ER)-negative breast cancers. Individual and pooled survival analyses were conducted on six independent large-scale breast cancer microarray data sets; and findings were validated on a human breast tissue set (ZJU set). Gene set enrichment analysis revealed that RRM2-high breast cancers were significantly enriched for expression of gene sets that increased in proliferation, invasiveness, undifferentiation, embryonic stem/progenitor-like phenotypes, and poor patient survival (p < 0.01). Independent and pooled analyses verified that increased RRM2 mRNA levels were associated with poor patient outcome in a dose-dependent manner. The prognostic power of RRM2 mRNA was comparable to multiple gene signatures, and it was superior to TNM stage. In ER-negative breast cancers, RRM2 showed more prognostic power than that in ER-positive breast cancers. Further analysis indicated that RRM2 was a more accurate prognostic biomarker for ER-negative breast cancers than the pathoclinical indicators and uPA. A new RR inhibitor, COH29, could significantly enhance the chemosensitivity to doxorubicin in ER-negative MDA-MB-231 cells, but not in ER-positive MCF-7 cells. RRM2 appears to be a promising prognostic biomarker and therapeutic target for ER-negative breast cancer patients.